Beyond the imbalance: site-resolved dynamics probing resonances in many-body localization

Resumo Técnico: Além do Desequilíbrio: Sondagem de Dinâmicas Resolvidas por Sítio em Ressonâncias na Localização de Muitos Corpos

Problema
A Localização de Muitos Corpos (MBL) é um mecanismo que impede a termalização quântica em sistemas desordenados interagentes. Embora o diagnóstico padrão para MBL seja o "desequilíbrio" ($I_m(t)$) — uma medida espacialmente média da memória de um estado inicial de onda de densidade — este artigo argumenta que tal média espacial obscurece características microscópicas cruciais. Especificamente, os autores investigam se o observável de desequilíbrio padrão esconde fenômenos dinâmicos locais, como ressonâncias, que são críticos para entender a transição entre fases ergódicas e localizadas. O trabalho aborda as limitações de interpretar o decaimento do desequilíbrio de longo tempo apenas através de médias espaciais, observando que decai de forma ultra-lenta (lei de potência) em simulações numéricas, o que tem sido interpretado de forma controversa como sinais de uma eventual termalização.

Metodologia
Os autores focam no regime de desordem forte da cadeia de spins XXZ de campo aleatório ($L$ sítios, força de desordem $h$, interação $\Delta$). Eles empregam duas técnicas numéricas primárias:

1. Diagonalização Exata (ED): Utilizada para sistemas de até $L=18$ (no setor de magnetização zero) para computar autovetores exatos e limites de tempo infinito.
2. Evolução Temporal de Krylov: Utilizada para sistemas maiores ($L$ até 24) para simular a evolução unitária até $t \approx 10^4$.

Em vez de depender apenas do desequilíbrio padrão $I_m(t) = \frac{1}{L} \sum_j Z^m_j(t)$, o estudo centra-se no autocorreletor local resolvido por sítio $Z^m_j(t) = \langle m | \sigma^z_j(0) \sigma^z_j(t) | m \rangle$. Os autores analisam sistematicamente as distribuições dessas quantidades locais sob diferentes protocolos de média:

• Média espacial: Soma sobre os sítios $j$.
• Média do estado inicial: Soma sobre $N_s$ estados produto distintos (traço parcial) ou o espaço de Hilbert completo (traço total/temperatura infinita).
• Média de desordem: Soma sobre $N_r$ realizações de campos aleatórios.

Para interpretar os achados numéricos, os autores constroem um modelo de brinquedo analítico de poucos sítios. Este modelo trata as ressonâncias locais como clusters efetivos de poucos corpos (2 e 3 sítios) inseridos em um fundo de spins congelados (presos), permitindo a derivação de previsões analíticas para os limites de longo tempo de observáveis locais.

Contribuições Principais e Resultados

1. Mascaramento de Ressonâncias Locais pela Média Espacial:
   O estudo demonstra que o desequilíbrio padrão $I_m(t)$ exibe um patamar estável consistente com MBL, mas essa média esconde completamente o comportamento de sítios individuais. Enquanto a maioria dos sítios permanece polarizada (retendo memória do estado inicial), uma pequena fração de sítios exibe flutuações fortes e decai para zero. Esses sítios "flutuantes" correspondem a regiões locais onde os campos aleatórios são quase idênticos (ex: $h_j \approx h_{j+1}$), criando ressonâncias locais.

2. Estrutura de Múltiplos Picos nas Distribuições Locais:
   Ao analisar os histogramas dos autocorreletores locais de longo tempo ($P(Z^m_j)$ e $P(\mathcal{A}_j)$), os autores revelam uma estrutura rica obscurecida pelas médias.

   • A distribuição de observáveis resolvidos por sítio mostra um pico dominante próximo a 1 (spins congelados) e picos secundários distintos em valores menores (ex: $\approx 0,6$ para ressonâncias de 2 corpos e $\approx 0,46$ para ressonâncias de 3 corpos no caso de Heisenberg).
   • Esses picos secundários são assinaturas de ressonâncias de poucos corpos onde os spins invertem de forma coerente devido a níveis de energia quase degenerados.

3. Validação Analítica via Modelo de Brinquedo:
   Os autores derivam expressões de forma fechada para os limites de longo tempo dessas ressonâncias usando um Hamiltoniano efetivo de 2 e 3 sítios.

   • Para ressonâncias de 2 corpos com $\Delta=1$, a previsão analítica para a posição do pico é $A \approx \frac{5}{8} - \frac{3}{8(L-1)}$.
   • Para ressonâncias de 3 corpos, são previstos picos distintos para sítios centrais e de borda.
   • Essas previsões analíticas alinham-se quantitativamente com os histogramas numéricos, confirmando que os picos secundários surgem de clusters de ressonância de poucos corpos específicos.

4. Escalonamento de Tamanho Finito e Dependência do Estado Inicial:
   Um achado importante é a forte dependência do desequilíbrio de longo tempo em relação ao estado inicial e ao tamanho do sistema, impulsionada por essas ressonâncias:

   • Estados Aleatórios Genéricos: O desequilíbrio médio aproxima-se do limite termodinâmico por baixo (correções negativas de $1/L$). Isso é impulsionado pela ocorrência rara de sítios ressonantes que reduzem a média.
   • Estado Néel ($|\uparrow\downarrow\uparrow\downarrow\dots\rangle$): O desequilíbrio aproxima-se do limite termodinâmico por cima (correções positivas de $1/L$). Como o estado Néel consiste inteiramente de ligações invertíveis ($\uparrow\downarrow$), ele maximiza o acoplamento aos modos ressonantes, levando a um comportamento de escalonamento de tamanho finito diferente comparado a estados aleatórios.
   • Estado de Parede de Domínio: Mostra um comportamento intermediário com correções específicas de tamanho finito.

5. Detectabilidade Experimental:
   O artigo estabelece que esses efeitos de ressonância são detectáveis sob condições experimentais realistas. Os picos secundários na distribuição de autocorreletores locais podem ser resolvidos com um número relativamente pequeno de estados iniciais ($N_s \sim 20-50$) e janelas de tempo moderadas ($t \approx 500$), desde que a imagem resolvida por sítio esteja disponível.

Significância e Alegações
O artigo afirma que o "desequilíbrio" não é um probe universal de MBL, mas é altamente sensível às condições iniciais específicas e ao método de média. Os autores argumentam que:

• Probes Refinados são Necessários: Para entender plenamente a MBL e sua conexão com efeitos de regiões raras, deve-se ir além de quantidades espacialmente médias para dinâmicas resolvidas por sítio.
• Reinterpretação da Dinâmica Lenta: O decaimento lento do desequilíbrio observado em estudos numéricos de tempo finito (frequentemente partindo de estados Néel) pode não indicar necessariamente uma termalização eventual. Em vez disso, pode representar um crossover transiente de um valor de tamanho finito (aproximado por cima) para um valor assintótico termodinâmico não nulo menor.
• Papel das Ressonâncias Raras: Ressonâncias de poucos corpos de curto alcance, ocorrendo com probabilidade escalonando como $1/h$ (para 2 corpos) e $1/h^2$ (para 3 corpos), governam o escalonamento de tamanho finito do desequilíbrio no regime de forte desordem.
• Relevância Experimental: Com o advento de imagens de alta precisão resolvidas por sítio em átomos ultrafrios e circuitos supercondutores, as estruturas de múltiplos picos e os escalonamentos de tamanho finito dependentes do estado inicial previstos são diretamente testáveis.

O trabalho conclui que uma "abordagem refinada, resolvida por sítio" é essencial para distinguir as complexidades da dinâmica de MBL, particularmente a interação entre o comportamento localizado e as instabilidades de ressonância raras.