Quantum Transport in Interacting Spin Chains: Exact Derivation of the GUE Tracy-Widom Distribution

Resumo Técnico: Transporte Quântico em Cadeias de Spin Interagentes: Derivação Exata da Distribuição GUE Tracy-Widom

Enunciado do Problema
O artigo aborda a questão de longa data sobre se a distribuição Tracy-Widom do GUE (Gaussian Unitary Ensemble), uma distribuição universal que caracteriza o maior autovalor de matrizes aleatórias e a corrente integrada na classe de universalidade de Kardar-Parisi-Zhang (KPZ), emerge na dinâmica de sistemas quânticos de muitos corpos interagentes. Embora esse comportamento tenha sido estabelecido anteriormente para o modelo XX não interagente (equivalente a férmions livres via transformação de Jordan-Wigner), sua presença em modelos interagentes como o modelo XXZ tem sido elusiva. Estudos numéricos anteriores no modelo XXZ produziram resultados conflitantes, com alguns sugerindo a ausência do comportamento GUE Tracy-Widom e outros indicando sua presença. Além disso, derivações analíticas exatas para modelos interagentes que não podem ser mapeados para férmions não interagentes eram escassas.

Metodologia
Os autores empregam uma abordagem de duas frentes combinando soluções analíticas exatas e simulações numéricas:

1. Derivação Analítica Exata (Modelo XXZ Dobrado):

   • Modelo: O estudo foca no modelo XXZ "dobrado", um Hamiltoniano efetivo derivado do modelo XXZ padrão no limite de grande anisotropia ($\Delta \gg 1$). Crucialmente, este modelo não pode ser mapeado para um Hamiltoniano de férmions não interagentes via transformação de Jordan-Wigner.
   • Estado Inicial: Um estado de parede de domínio alternante é considerado, onde spins para cima ocupam cada outro sítio em metade do sistema ($x \ge 2$).
   • Técnica: Os autores utilizam a técnica de Bethe ansatz, especificamente adaptando métodos desenvolvidos por Schütz, Tracy e Widom para o processo de exclusão simples assimétrico e o modelo XXZ.
   • Etapas da Derivação:
     • Eles derivam uma fórmula integral para a função de onda de muitos corpos $\Phi(x_1, \dots, x_N, t)$.
     • Usando o Bethe ansatz e propriedades do modelo de seis vértices (determinante de Izergin-Korepin), eles derivam uma fórmula determinante para a função de distribuição cumulativa complementar $F(x,t)$ de encontrar o spin para cima mais à esquerda no sítio $x$.
     • Ao tomar o limite $\Delta \to \infty$, eles simplificam as amplitudes de espalhamento para obter uma forma determinante específica compatível com a teoria de matrizes aleatórias.
     • A análise assintótica é realizada usando uma variável de escala $s$ definida por $x = 2 + \lfloor -t - s(t/2)^{1/3} \rfloor$. Isso envolve analisar o núcleo do determinante no limite de longo tempo ($t \to \infty$) usando técnicas relacionadas a operadores de Toeplitz e funções de Bessel, convergindo finalmente para o núcleo de Airy.

2. Investigação Numérica (Modelo XXZ Padrão):

   • Modelo: O modelo XXZ padrão com condições de contorno abertas e anisotropia finita $\Delta$.
   • Técnica: A Decimação de Bloco Evolutiva no Tempo (TEBD) é usada para simular a dinâmica em tempo real da equação de Schrödinger.
   • Análise: A probabilidade $P(x,t)$ de encontrar o spin para cima mais à esquerda é computada para vários valores de $\Delta$ (5, 10, 20). Os dados são reescalonados de acordo com os expoentes de escala KPZ para comparar contra a distribuição GUE Tracy-Widom.

Principais Contribuições e Resultados

1. Derivação Exata para Modelos Interagentes: A principal contribuição é a primeira derivação exata demonstrando que a distribuição de probabilidade do spin para cima mais à esquerda em um modelo de spin quântico interagente (o modelo XXZ dobrado) converge para a distribuição GUE Tracy-Widom ($F_2(s)$) no limite de longo tempo. Isso confirma que o comportamento GUE Tracy-Widom sobrevive em sistemas que não são mapeáveis para férmions não interagentes.
2. Conexão Matemática: Os autores estabelecem um elo matemático entre o modelo XXZ dobrado e o modelo XX. Eles mostram que a função de onda de muitos corpos do modelo XXZ dobrado possui uma estrutura determinante semelhante à do modelo XX, fornecendo uma origem matemática para a emergência da mesma distribuição universal, apesar da presença de interações.
3. Assinaturas Numéricas no XXZ: Para o modelo XXZ padrão com $\Delta$ finito, os resultados numéricos mostram que, embora a distribuição completa não combine perfeitamente com a curva GUE Tracy-Widom para $\Delta$ pequenos em tempos acessíveis, uma "assinatura universal" específica emerge. Especificamente, na cauda direita da distribuição ($s > 0$), a probabilidade $P(x,t)$ exibe uma curva independente de $\Delta$. Esta curva é caracterizada pelo núcleo de Airy diagonal $K_{Ai}(s,s)$.
4. Escalonamento de Convergência Rápida: O artigo introduz uma variável de escala refinada com um parâmetro $a=1/2$ (em vez do padrão $a=0$) para o reescalonamento de coordenadas. Testes numéricos demonstram que esta escolha leva a uma convergência significativamente mais rápida da distribuição de tempo finito para o limite assintótico de Tracy-Widom.

Significância e Alegações
O artigo afirma fornecer a primeira prova exata de que a distribuição GUE Tracy-Widom surge na dinâmica de um modelo quântico interagente que não pode ser reduzido a férmions livres. Isso resolve a ambiguidade sobre a universalidade da física KPZ em sistemas quânticos interagentes.

Os autores sugerem modestamente que suas descobertas oferecem uma base teórica para observar esses comportamentos universais em configurações experimentais atuais. Eles observam que experimentos recentes utilizando átomos ultra-frios e qubits supercondutores realizaram modelos do tipo XXZ e acessaram escalas de tempo onde os spins invertem múltiplas vezes (aproximadamente 5 vezes), o que se alinha com as escalas de tempo onde a assinatura GUE Tracy-Widom aparece em suas simulações numéricas para $\Delta$ grande.

O trabalho abre direções futuras para explorar distribuições Tracy-Widom em modelos não integráveis e outros sistemas integrantes, como o modelo de fase, e sugere que a hidrodinâmica generalizada e a teoria de flutuação macroscópica balística poderiam ser usadas para entender analiticamente o comportamento universal do modelo XXZ.