gerar uma matriz
N times N cujos elementos sejam números aleatórios gaussianos independentes com média zero e desvio padrão
sigma=1 ;
tomar a média de cada matriz com sua transposta para obter uma matriz simétrica.
Uma das propriedades mais notáveis que grandes matrizes aleatórias compartilham é a distribuição dos espaçamentos entre níveis.
(a) Gere um ensemble com
M=10000 matrizes do GOE de tamanhos
N=2 , 4 e 10. Para cada matriz, determine os autovalores
lambdan , ordenados do menor para o maior, e calcule a diferença
lambdan+1-lambdan entre autovalores consecutivos, mas apenas com
n=N / 2 . (Utilizar as diferenças entre todos os autovalores consecutivos distorceria um pouco a distribuição e dificultaria a visualização nas próximas etapas.) Para cada valor de
N , trace histogramas do espaçamento entre autovalores dividido pelo valor médio desses espaçamentos, utilizando caixas de largura suficientemente pequena para conseguir observar algumas flutuações. Você deve observar uma queda na altura dos histogramas quando o espaçamento tende a zero, que caracteriza o que se chama de repulsão entre níveis. (Se for utilizar o Mathematica, as funções RandomVariate, NormalDistribution, Transpose, Eigenvalues, Sort e Histogram são convenientes; veja a documentação do software.)
Para
N=2 , a distribuição de probabilidades dos espaçamentos pode ser calculada analiticamente. Vamos denotar a matriz por
H=(lla & b b & c) .
(b) Mostre que a diferença entre os autovalores de
H é
De .
© As distribuições de probabilidade de
a e
c são gaussianas com média zero e desvio padrão
sigmaa=sigmac=1 , ou seja
Pa(u)=Pc(u)=1raiz2 pi e^-u^2 / 2
Lembrando que a matriz
H é simetrizada por construção, e que portanto
b corresponde à média de dois números aleatórios gaussianos, cada um com média zero e desvio padrão 1, mostre que a distribuição de probabilidade de
b é
(d) Determine a distribuição de probabilidade de
De , ou seja, calcule
PDe .
Dica: efetue primeiro a transformação de variáveis
x=(c+a) / 2 e
y=(c-a) / 2 , integre sobre
x e em seguida efetue uma segunda transformação para variáveis polares
r e
theta tais que
y=r sen theta e
b=r cos theta .
(e) Reescale a distribuição do item anterior definindo
s=De para fazer o espaçamento médio igual a 1 e mostre que o resultado é
rho Wigner (s)=pi s2 e^-pi s^2 / 4,
que é conhecida como distribuição de espaçamentos de Wigner. Note que tanto
De quanto
s estão definidos apenas no intervalo
[0,+infty) .
(f) Apesar de ter sido deduzida para
N=2 , a distribuição de espaçamentos de Wigner é uma boa aproximação para os resultados obtidos com
N maior. Para verificar essa afirmação, trace os histogramas do item (a) junto com
rho Wigner (s) em gráficos para
N=4 e
N=10 .
(g) Por fim, repita o item (a) para um outro ensemble, agora com matrizes simétricas cujos elementos sejam variáveis aleatórias
pm 1 com igual probabilidade. As distribuições obtidas são quase universais, ou seja, são praticamente as mesmas para qualquer
N , como ocorreu no item (a)? Para
N=10 , trace o histograma obtido e superponha
rho Wigner (s) no mesmo gráfico. A concordância nesse caso parece razoável? O que isso sugere? Teste sua hipótese com outros ensembles de matrizes aleatórias simétricas que lhe ocorram.