Esta tese tem como objetivo o estudo de aspectos dinâmicos e estatísticos da matéria escura em distribuições esféricas de massa. O fato de suas partículas constituintes interagirem gravitacionalmente mas não eletromagneticamente, e portanto sua evolução ser regida por interações de longo alcance, traz algumas complicações teóricas na descrição de suas propriedades nos termos da mecânica estatística, dificuldades compartilhadas com sistemas auto-gravitantes em geral. Para melhor compreender essas propriedades, estudamos as distribuições de matéria escura em três abordagens diferentes. Na primeira, utilizamos dados observacionais, utilizando lentes gravitacionais, em aglomerados de galáxias para comparar a performance de alguns modelos propostos para o perfil de densidade da matéria escura. Dividimos estes modelos em fenomenológicos ou teóricos. Dos primeiros, todos são capazes de descrever os dados observacionais com performance comparável. Entre os modelos teóricos estudados, o modelo chamado DARKexp descreve os dados tão bem quanto os primeiros. Numa segunda abordagem, utilizamos dados de simulações numéricas para testar uma função proposta para a distribuição de velocidades das partículas. Esta função inclui a anisotropia no campo de velocidades na chamada distribuição q-gaussiana. Comparamos a performance desta função com a da função gaussiana e concluímos que a primeira representa uma melhor descrição dos dados, mesmo levando em conta a introdução de um parâmetro extra, apesar de ainda apresentar algumas discrepâncias, especialmente nas regiões internas dos halos. Por fim, discutimos a possível relevância do conceito de indistinguibilidade na determinação dos estados de equilíbrio de sistemas auto-gravitantes em geral, propondo uma associação deste conceito com o nível de mistura do sistema. Implementamos esta associação numa análise combinatória e estudamos as conseqüências para a determinação da função distribuição e do perfil de densidades. Esta associação também levanta algumas dúvidas sobre a validade da equação de Vlasov durante o processo de relaxação violenta.