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

Resumen Técnico: Transporte Cuántico en Cadenas de Espín con Interacciones: Derivación Exacta de la Distribución de Tracy-Widom del GUE

Planteamiento del Problema
El artículo aborda la cuestión largamente debatida de si la distribución de Tracy-Widom del GUE (Conjunto Unitario Gaussiano), una distribución universal que caracteriza al mayor autovalor de matrices aleatorias y a la corriente integrada en la clase de universalidad de Kardar-Parisi-Zhang (KPZ), emerge en la dinámica de sistemas cuánticos de muchos cuerpos con interacciones. Si bien este comportamiento ya se había establecido para el modelo XX no interactuante (equivalente a fermiones libres mediante la transformación de Jordan-Wigner), su presencia en modelos con interacciones como el modelo XXZ ha permanecido esquiva. Estudios numéricos previos sobre el modelo XXZ arrojaron resultados contradictorios: algunos sugerían la ausencia del comportamiento de la distribución de Tracy-Widom del GUE y otros indicaban su presencia. Además, faltaban derivaciones analíticas exactas para modelos con interacciones que no pueden mapearse a fermiones no interactuantes.

Metodología
Los autores emplean un enfoque de dos vertientes que combina soluciones analíticas exactas y simulaciones numéricas:

1. Derivación Analítica Exacta (Modelo XXZ Plegado):

   • Modelo: El estudio se centra en el modelo XXZ "plegado" (folded), un Hamiltoniano efectivo derivado del modelo XXZ estándar en el límite de interacción anisotrópica grande ($\Delta \gg 1$). Crucialmente, este modelo no puede mapearse a un Hamiltoniano de fermiones no interactuantes mediante la transformación de Jordan-Wigner.
   • Estado Inicial: Se considera un estado de pared de dominio alternante, donde los espines "arriba" ocupan cada sitio alterno en la mitad del sistema ($x \ge 2$).
   • Técnica: Los autores utilizan la técnica de la aproximación de Bethe (Bethe ansatz), adaptando específicamente los métodos desarrollados por Schütz, Tracy y Widom para el proceso de exclusión simple asimétrica y el modelo XXZ.
   • Pasos de la Derivación:
     • Derivan una fórmula integral para la función de onda de muchos cuerpos $\Phi(x_1, \dots, x_N, t)$.
     • Utilizando la aproximación de Bethe y las propiedades del modelo de seis vértices (determinante de Izergin-Korepin), derivan una fórmula determinante de la función de distribución acumulada complementaria $F(x,t)$ para encontrar el espín "arriba" situado más a la izquierda en el sitio $x$.
     • Al tomar el límite $\Delta \to \infty$, simplifican las amplitudes de dispersión para obtener una forma determinante específica compatible con la teoría de matrices aleatorias.
     • Se realiza un análisis asintótico utilizando una variable de escala $s$ definida por $x = 2 + \lfloor -t - s(t/2)^{1/3} \rfloor$. Esto implica el análisis del núcleo del determinante en el límite de tiempo largo ($t \to \infty$) utilizando técnicas relacionadas con operadores de Toeplitz y funciones de Bessel, convergiendo finalmente al núcleo de Airy.

2. Investigación Numérica (Modelo XXZ Estándar):

   • Modelo: El modelo XXZ estándar con condiciones de contorno abiertas y anisotropía finita $\Delta$.
   • Técnica: Se utiliza la Decimación de Bloques Evolutiva en el Tiempo (TEBD) para simular la dinámica en tiempo real de la ecuación de Schrödinger.
   • Análisis: Se computa la probabilidad $P(x,t)$ de encontrar el espín "arriba" situado más a la izquierda para varios valores de $\Delta$ (5, 10, 20). Los datos se reescalan según los exponentes de escala KPZ para compararlos con la distribución de Tracy-Widom del GUE.

Contribuciones Clave y Resultados

1. Primera Derivación Exacta para Modelos con Interacción: La principal contribución es la primera derivación exacta que demuestra que la distribución de probabilidad del espín "arriba" situado más a la izquierda en un modelo de espín cuántico con interacción converge a la distribución de Tracy-Widom del GUE ($F_2(s)$) en el límite de tiempo largo. Esto confirma que el comportamiento de la distribución de Tracy-Widom del GUE sobrevive en sistemas que no son mapeables a fermiones no interactuantes.
2. Conexión Matemática: Los autores establecen un vínculo matemático entre el modelo XXZ plegado y el modelo XX. Demuestran que la función de onda de muchos cuerpos del modelo XXZ plegado posee una estructura determinante similar a la del modelo XX, proporcionando un origen matemático para la emergencia de la misma distribución universal a pesar de la presencia de interacciones.
3. Firmas Numéricas en el XXZ: Para el modelo XXZ estándar con $\Delta$ finito, los resultados numéricos muestran que, si bien la distribución completa no coincide perfectamente con la curva de Tracy-Widom del GUE para $\Delta$ pequeños en tiempos accesibles, emerge una "firma universal" específica. Concretamente, en la cola derecha de la distribución ($s > 0$), la probabilidad $P(x,t)$ exhibe una curva independiente de $\Delta$. Esta curva está caracterizada por el núcleo de Airy diagonal $K_{Ai}(s,s)$.
4. Escalamiento de Convergencia Rápida: El artículo introduce una variable de escala refinada con un parámetro $a=1/2$ (en lugar del estándar $a=0$) para el reescalamiento de la coordenada. Las pruebas numéricas demuestran que esta elección conduce a una convergencia significativamente más rápida de la distribución de tiempo finito hacia el límite asintótico de Tracy-Widom.

Significado y Reivindicaciones
El artículo afirma proporcionar la primera prueba exacta de que la distribución de Tracy-Widom del GUE surge en la dinámica de un modelo cuántico con interacción que no puede reducirse a fermiones libres. Esto resuelve la ambigüedad sobre la universalidad de la física KPZ en sistemas cuánticos con interacción.

Los autores sugieren modestamente que sus hallazgos ofrecen una base teórica para observar estos comportamientos universales en configuraciones experimentales actuales. Señalan que experimentos recientes utilizando átomos ultrafríos y cúbits superconductores han implementado modelos de tipo XXZ y han accedido a escalas de tiempo donde los espines cambian de orientación múltiples veces (aproximadamente 5 veces), lo cual se alinea con las escalas de tiempo donde aparece la firma de la distribución de Tracy-Widom del GUE en sus simulaciones numéricas para $\Delta$ grandes.

El trabajo abre futuras direcciones para explorar las distribuciones de Tracy-Widom en modelos no integrables y otros sistemas integrables, como el modelo de fase, y sugiere que la hidrodinámica generalizada y la teoría de fluctuación macroscópica balística podrían utilizarse para comprender analíticamente el comportamiento universal del modelo XXZ.