Quantum chaos and late-time equipartition of symmetry-resolved Krylov complexity

Resumo Técnico: Caos Quântico e Equipartição de Longo Prazo da Complexidade de Krylov Resolvida por Simetria

Enunciado do Problema
A complexidade de Krylov ($C_K$) emergiu como um diagnóstico para o crescimento de operadores em sistemas de muitos corpos quânticos, tipicamente caracterizada por um crescimento exponencial inicial em sistemas caóticos, seguido por um regime linear e eventual saturação. No entanto, a interpretação do valor de saturação de longo prazo é complicada pela presença de simetrias exatas. Em sistemas onde o Hamiltoniano $H$ e o operador semente inicial $O$ comutam com uma carga conservada $P$, o espaço de Hilbert se decompõe em setores de simetria independentes. Embora esteja estabelecido que, em tempos iniciais, a complexidade de Krylov total relaciona-se com uma média ponderada das complexidades resolvidas por simetria (com pesos determinados pela projeção do operador inicial), a relação entre as complexidades total e resolvida em tempos tardios (apressando a saturação) permanecia incerta. Uma expectativa simples poderia sugerir que o valor de saturação escala com o quadrado da dimensão total do espaço de Hilbert ($D^2$), mas a presença de simetrias poderia alterar essa escala, potencialmente obscurecendo o poder de diagnóstico da $C_K$ para o caos quântico se não for devidamente resolvida.

Metodologia
Os autores empregam uma combinação de derivação analítica e simulação numérica para investigar o comportamento de longo prazo da complexidade de Krylov resolvida por simetria em sistemas caóticos de dimensão finita.

• Estrutura Analítica: O estudo utiliza o formalismo do espaço de Krylov gerado pelo superoperador de Liouville $\mathcal{L} = [H, \cdot]$. Ao decompor o espaço de Hilbert $\mathcal{H}$ em setores de simetria $\mathcal{H}_q$ com dimensões $d_q$, os autores analisam a dimensão do espaço de Krylov acessível dentro de cada setor ($K_q$) e o espaço total não resolvido ($K$). Eles assumem que, para Hamiltonianos caóticos genéricos, as diferenças de energia são não-degeneradas (exceto pelo modo zero universal) e que o operador semente possui sobreposição não nula com todas as frequências de Liouville permitidas pela simetria. Sob essas suposições, eles derivam as dimensões dos espaços de Krylov e aplicam a hipótese de deslocalização de longo prazo, onde a função de onda de Krylov torna-se uniformemente distribuída sobre a cadeia de Krylov acessível.
• Implementação Numérica: Os autores implementam o algoritmo de Lanczos (usando a iteração de Arnoldi para estabilidade numérica) para computar a complexidade de Krylov. Eles testam suas previsões analíticas em quatro modelos distintos:
  1. Modelo SYK Real: Um modelo de férmions de Majorana com simetria de paridade de férmion $Z_2$.
  2. Modelo SYK Complexo: Um modelo de férmions complexos com conservação de carga global $U(1)$.
  3. Modelo de Spin Bosônico: Um modelo de spin quadrático caótico com simetrias $Z_2 \times Z_2$.
  4. Cadeia de Ising de Campo Misto: Uma cadeia de spin não-integrável com simetria de reflexão.
     Eles também analisam um limite integrável da cadeia de Ising para contrastar comportamentos caóticos e integráveis.

Principais Contribuições e Resultados
O artigo estabelece que no regime de saturação ($t_{sat} \lesssim t \ll t_{rec}$), a complexidade de Krylov total torna-se aditiva sobre os setores de simetria:
$$C_K(t) \simeq \sum_q C_K^{(q)}(t)$$
Crucialmente, a contribuição fracionária de cada setor não é determinada pelos pesos de projeção iniciais (como no regime de tempo inicial) nem por partilha igual, mas sim pelas dimensões dos espaços de operadores dinamicamente acessíveis. Especificamente, a complexidade fracionária de longo prazo é governada por:
$$\frac{C_K^{(q)}}{C_K} \simeq \frac{d_q(d_q - 1)}{\sum_{q'} d_{q'}(d_{q'} - 1)} \approx \frac{d_q^2}{\sum_{q'} d_{q'}^2}$$
onde $d_q$ é a dimensão do espaço de Hilbert do setor $q$. Este fenômeno é denominado equipartição ponderada pela dimensão.

• Verificação Numérica:
  • No modelo SYK Real ($N=16$), os dois setores de paridade possuem dimensões iguais ($d_+ = d_-$), resultando em uma divisão igual ($50\%$) da complexidade de saturação, consistente com a escala $d_q^2$.
  • No modelo SYK Complexo ($N=7, 8$), os setores de carga possuem dimensões desiguais ($d_q = \binom{N}{q}$). Os resultados numéricos mostram que os setores com maiores $d_q$ contribuem com uma fração maior da complexidade total, seguindo precisamente o peso $d_q^2$.
  • O modelo de spin bosônico ($N=6$) com quatro setores iguais confirma a divisão igual ($25\%$).
  • A cadeia de Ising de campo misto ($N=8$) apresenta setores de reflexão com dimensões desiguais ($d_+ = 136, d_- = 120$). Os valores de saturação desviam-se de $50/50$ e, em vez disso, correspondem à razão prevista $d_+^2 / (d_+^2 + d_-^2) \approx 0.56$.
• Contraste Integrável: No limite integrável da cadeia de Ising, os autores demonstram que as degenerescências de Liouville e as regras de seleção reduzem as dimensões reais de Krylov ($K_q$) muito abaixo dos limites genéricos caóticos ($d_q(d_q-1)+1$). Consequentemente, a complexidade fracionária de longo prazo satura de acordo com a razão das verdadeiras dimensões de Krylov $(K_q-1)/(K_{full}-1)$ em vez da escala $d_q^2$, destacando que a regra de ponderação pela dimensão é específica para a dinâmica caótica, onde o espaço de Krylov é maximamente explorado.

Significância
O artigo argumenta que resolver as simetrias exatas é essencial para interpretar corretamente o valor de saturação da complexidade de Krylov como um diagnóstico do crescimento de operadores caóticos. Sem a resolução de simetria, poder-se-ia assumir incorretamente que a saturação escala com o quadrado da dimensão total do espaço de Hilbert ($D^2$). Em vez disso, os autores mostram que, para sistemas caóticos, a saturação escala com a soma dos quadrados das dimensões dos setores ($\sum d_q^2$), que é parametricamente menor que $D^2$ se os setores de simetria forem desiguais.

O trabalho fornece uma "relação genérica simples" entre as complexidades de Krylov total e resolvida por simetria em sistemas caóticos, distinguindo o regime de longo prazo (governado pelas dimensões dos espaços de operadores acessíveis) do regime de curto prazo (governado pelas projeções do estado inicial). Isso esclarece a interpretação da saturação da complexidade de Krylov, garantindo que o diagnóstico reflita a real extensão do crescimento de operadores dentro dos subespaços dinamicamente acessíveis, em vez de ser obscurecido por estruturas de blocos induzidas pela simetria. As descobertas sugerem que a resolução de simetria atua como um refinamento do diagnóstico de saturação, oferecendo uma caracterização mais nítida da dinâmica caótica.