Circular Rosenzweig-Porter random matrix ensemble

Resumo Técnico: Ensemble Circular de Matrizes Aleatórias de Rosenzweig-Porter

Enunciado do Problema
O ensemble de matrizes aleatórias de Rosenzweig-Porter (RP) é um modelo fenomenológico bem estabelecido utilizado para descrever as estatísticas de níveis e a fractalidade dos autoestados através da transição de localização de muitos corpos (MBL) em sistemas quânticos interagentes, desordenados e estáticos. Em sistemas estáticos, a dinâmica é governada por Hamiltonianos hermitianos. No entanto, sistemas periodicamente acionados (Floquet) são descritos por operadores de Floquet unitários, onde os autovalores residem no círculo unitário no plano complexo. Embora a MBL tenha sido estendida para sistemas de Floquet, um análogo unitário (circular) direto do ensemble RP, capaz de capturar a mesma fenomenologia para sistemas periodicamente acionados, não havia sido construído. Este artigo aborda a questão de se tal análogo unitário pode ser definido e se ele compartilha as propriedades estatísticas fundamentais do ensemble RP original.

Metodologia
Os autores propõem uma construção do ensemble Circular de Rosenzweig-Porter (CRP) baseada na observação de que o ensemble RP padrão pode ser visto como o resultado de um processo de movimento browniano de Dyson (DBM) de tempo finito.

1. Construção: Os autores definem o ensemble CRP como o resultado de uma evolução estocástica de uma matriz unitária simétrica dependente do tempo $S(t)$ de dimensão $N$. A evolução segue um processo de movimento browniano de Dyson adaptado para matrizes unitárias:
   $$dS(t) = i\sqrt{dt} U^T(t) X U(t)$$
   onde $X$ é uma matriz do Ensemble Gaussiano Ortogonal (GOE) amostrada a cada passo de tempo infinitesimal, e $U(t)$ é derivada da decomposição $S(t) = U^T(t)U(t)$.
2. Condições Iniciais: O processo inicia com uma matriz diagonal $S(0) = \text{diag}(e^{i\theta_1(0)}, \dots, e^{i\theta_N(0)})$, onde as fases iniciais $\theta_n(0)$ são amostradas independentemente de uma distribuição uniforme sobre $[-\pi, \pi)$. Isso garante uma densidade de estados uniforme, consistente com operadores de Floquet genéricos.
3. Definição do Ensemble: O ensemble CRP é definido como o estado da matriz $S(t)$ em um tempo específico $t = \epsilon^2 / N^\gamma$, onde $\epsilon$ é a força de perturbação e $\gamma$ controla a força relativa entre os termos diagonais e fora da diagonal.
4. Implementação Numérica: Os autores avaliam numericamente o processo estocástico utilizando um algoritmo que atualiza a matriz via $S(t+dt) = S(t)e^{i\sqrt{dt}X}$. Esta abordagem preserva a unitariedade para passos de tempo finitos e evita decomposições explícitas de matrizes. Simulações foram realizadas para dimensões de matriz $N \leq 1000$, com média sobre pelo menos $10^6$ autovalores ou autoestados.

Principais Contribuições

• Definição do CRP: O artigo introduz o ensemble Circular de Rosenzweig-Porter como um análogo unitário do ensemble RP padrão, definido explicitamente através de um processo de movimento browniano de Dyson circular.
• Mapeamento Teórico: Os autores fornecem um argumento heurístico de que as estatísticas espectrais do ensemble CRP no meio do espectro devem coincidir com as do ensemble RP, desde que a densidade de estados seja escalada adequadamente (especificamente, mapeando o parâmetro RP $\epsilon$ para $\epsilon/\sqrt{2\pi}$ no contexto do CRP).
• Verificação Numérica: O estudo fornece evidências numéricas abrangentes caracterizando as propriedades estatísticas dos autovalores e autoestados do ensemble CRP.

Resultados
A investigação numérica confirma que o ensemble CRP exibe comportamentos estatísticos análogos ao ensemble RP em diferentes regimes do parâmetro $\gamma$:

1. Estatísticas de Espaçamento de Níveis: A razão média de espaçamentos de níveis consecutivos, $r$, serve como um diagnóstico para o caos quântico.
   • Para $\gamma < 2$, o sistema exibe estatísticas de Wigner-Dyson ($r \approx 0.530$), características de sistemas caóticos.
   • Para $\gamma \geq 2$, o sistema transita para estatísticas de Poisson ($r \approx 0.386$), características de sistemas integráveis.
   • Uma transição de fase ocorre no ponto crítico $\gamma_c = 2$ com um expoente crítico $\nu = 1$. Os dados mostram um colapso de escala de $r$ em função de $(\gamma - \gamma_c) \ln(N)^{1/\nu}$, consistente com o ensemble RP.
2. Fractalidade dos Autoestados: O inverso do índice de participação (IPR) dos autoestados foi analisado na base onde $S(0)$ é diagonal.
   • Para $\gamma \leq 1$, os autoestados são estendidos ($IPR_2 \sim N^{-1}$).
   • Para $\gamma \geq 2$, os autoestados são localizados ($IPR_2 \sim O(1)$).
   • No regime intermediário $1 < \gamma < 2$, os autoestados são fractais, escalando como $IPR_2 \sim N^{-(2-\gamma)}$.
3. Forma dos Autoestados: A distribuição das amplitudes dos autoestados $|\langle m | \psi_n \rangle|^2$ foi encontrada seguindo estatísticas de Breit-Wigner. Os dados alinham-se com a previsão teórica:
   $$|\langle m | \psi_n \rangle|^2 \sim \frac{1}{(\theta_n - \theta_m^{(0)})^2 + \Gamma^2}$$
   onde $\Gamma$ é a largura de espalhamento. Embora efeitos de tamanho finito e aproximações de contínuo levem a discrepâncias quantitativas leves, o acordo qualitativo e as diferenças estruturais distintas entre os regimes $\gamma = 0.5$ e $\gamma = 1.5$ são claramente observados.

Significado e Alegações
Os autores afirmam que o ensemble Circular de Rosenzweig-Porter serve como um modelo fenomenológico válido para as estatísticas de níveis e a fractalidade dos autoestados observadas através da transição de localização de muitos corpos em sistemas periodicamente acionados (Floquet). Ao estabelecer que o ensemble CRP compartilha propriedades estatísticas fundamentais com o ensemble RP, o trabalho fornece uma ferramenta teórica tratável para estudar fases não ergódicas em sistemas quânticos acionados.

O artigo enquadra modestamente este trabalho como um "primeiro passo" rumo à construção de generalizações mais amplas, como modelos incorporando multifractalidade ou outros modelos de Floquet com autoestados multifractais. Os autores sugerem que futuras generalizações poderiam ser alcançadas considerando processos estocásticos com incrementos correlacionados, semelhantes a propostas recentes para o ensemble RP padrão. Adicionalmente, eles notam a utilidade potencial do ensemble CRP como um bloco construtivo não maximamente aleatório em estudos de circuitos quânticos aleatórios.