Kardar-Parisi-Zhang dynamics in an open integrable system: beyond the spontaneous-symmetry-breaking ansatz

Resumen Técnico: Dinámica de Kardar-Parisi-Zhang en un sistema integrable abierto

Planteamiento del Problema
La universalidad de las leyes de escala dinámica es una piedra angular para comprender los sistemas de muchos cuerpos cuánticos, particularmente en entornos fuera del equilibrio. El progreso reciente en sistemas cuánticos abiertos ha dependido fuertemente del Ansatz de Ruptura Espontánea de Simetría (SSBA, por sus siglas en inglés) [20], un marco fenomenológico que utiliza funciones de onda variacionales de Feynman-Bijl para aproximar las excitaciones de baja energía cerca de los estados estacionarios. El SSBA predice un transporte difusivo universal (exponente dinámico $z=2$) para sistemas abiertos con interacciones locales y sin restricciones, capturando con éxito comportamientos en modelos de interacción de largo alcance y modelos con restricciones. Sin embargo, sigue sin estar claro bajo qué condiciones este ansatz falla al capturar la hidrodinámica emergente y, específicamente, si emergen dinámicas no triviales cuando la aproximación se rompe. El SSBA descuida términos en la ecuación de Lindblad que impulsan las funciones de onda variacionales fuera del subespacio variacional, dejando abierta la cuestión de si pueden surgir dinámicas anómalas en sistemas integrables abiertos.

Metodología
Los autores investigan el modelo B3 [28–30], un sistema cuántico abierto integrable de tipo Yang-Baxter gobernado por la ecuación de Lindblad. El estudio emplea una combinación de transformaciones analíticas y simulaciones numéricas:

1. Isomorfismo de Choi: La ecuación de Lindblad se mapea a un espacio de Hilbert duplicado, descomponiendo el superoperador de Liouville $\mathcal{L}$ en tres componentes: evoluciones no hermitianas independientes en cadenas "izquierda" y "derecha" ($A \otimes \mathbb{1}$ y $\mathbb{1} \otimes B$) y un término de acoplamiento entre cadenas ($l \otimes l^*$) que representa saltos cuánticos.
2. Análisis Variacional: Los autores analizan el SSBA dentro de este espacio duplicado, demostrando que los términos de cadena única ($A$ y $B$) impulsan el sistema fuera del subespacio variacional utilizado por el ansatz, mientras que el ansatz solo captura los efectos del acoplamiento entre cadenas.
3. Selección del Estado Inicial: El sistema se prepara en un estado inicial específico $\rho_0 = 2n_{j'} |\Psi\rangle\langle\Psi| n_{j'}$, donde $|\Psi\rangle$ es un estado producto que satisface $l^* |\Psi^*\rangle = 0$. Esta condición desacopla efectivamente las dos cadenas, haciendo que el acoplamiento entre cadenas sea dinámicamente irrelevante para la dinámica de transporte de carga.
4. Simulación Numérica: Se utiliza la Decimación de Bloques Evolutiva en el Tiempo (TEBD) para simular el transporte de carga para una cadena de longitud $N=256$ con una dimensión de enlace $\chi=64$ y un parámetro $\gamma=1.2$.
5. Mapeo Analítico: La dinámica de cadena única se mapea analíticamente a un modelo de cadena de espines XXZ asimétrica. El gap espectral de este Hamiltoniano efectivo se calcula utilizando el ansatz de Bethe para determinar el exponente dinámico.

Contribuciones Clave y Resultados

• Fallo del SSBA: El artículo demuestra que el SSBA falla al describir el transporte de carga en el modelo B3. Mientras que el SSBA predice un comportamiento difusivo ($z=2$) basado en el gap espectral del Lioviliano completo, la dinámica de carga real está gobernada por los términos de cadena única, los cuales el ansatz excluye.
• Emergencia de la Escala KPZ: Para estados iniciales apropiados, el modelo B3 exhibe un comportamiento de escala de Kardar-Parisi-Zhang (KPZ) ($z=3/2$) en lugar de difusión. Los resultados numéricos confirman que la función de correlación de densidad de carga $G(x,t)$ decae como $t^{-2/3}$ y colapsa sobre la función de escala universal de KPZ $f_{KPZ}(x/t^{2/3})$.
• Mecanismo de Desacoplamiento: El estudio revela que, bajo condiciones iniciales específicas, el acoplamiento entre cadenas (que genera el espectro difusivo predicho por el SSBA) se vuelve irrelevante. La dinámica está dominada por la evolución de la cadena única, la cual es equivalente a un modelo XXZ asimétrico.
• Origen de la XXZ Asimétrica: Se muestra que el operador de cadena única $A_{j,j+1}$ es equivalente al Hamiltoniano de una cadena XXZ asimétrica. Incluso cuando la tasa de salto se vuelve negativa (para $\gamma < 1$), el sistema se mapea a un proceso estadístico descrito por una ecuación de Burgers (el límite de campo medio de KPZ), produciendo el exponente $z=3/2$.
• Discrepancia del Gap Espectral: Los autores calculan el gap espectral del Hamiltoniano de cadena única efectiva ($H_{XXZ}$) mediante el ansatz de Bethe, encontrando un escalamiento de $\Delta \propto N^{-3/2}$. Esto contrasta con el escalamiento $N^{-2}$ del gap del Lioviliano completo, explicando por qué el SSBA (que sondea el Lioviliano completo) pierde la clase de universalidad KPZ.
• Tasas de Salto Negativas: El trabajo destaca que la tasa de salto efectiva en el modelo XXZ mapeado puede ser negativa, pero la constante de difusión permanece positiva, permitiendo que la ecuación de Burgers sirva como una teoría efectiva.

Significancia
El artículo sostiene que una comprensión simplista de la hidrodinámica de sistemas abiertos basada únicamente en la ruptura espontánea de simetría (SSB) es insuficiente. Establece que los sistemas abiertos integrables pueden exhibir dinámicas de KPZ, una clase de universalidad previamente asociada con modelos integrables aislados (como el imán de Heisenberg) o procesos estocásticos clásicos (ASEP), pero no típicamente con sistemas abiertos gobernados por el SSBA.

El trabajo proporciona una vía para conectar la mecánica estadística clásica con los sistemas cuánticos abiertos: si el acoplamiento entre cadenas desaparece y el operador de cadena única corresponde a un generador de un modelo estadístico, las funciones de correlación cuánticas coinciden con las funciones de Green de ese proceso clásico. Esto sugiere que el modelo B3 se reduce al Proceso de Exclusión Simple Asimétrica (ASEP) en el límite apropiado, ofreciendo un ejemplo concreto de cómo las herramientas estadísticas clásicas pueden aplicarse a la dinámica cuántica abierta. Los autores señalan que la emergencia de la escala KPZ incluso con tasas de salto negativas apunta a nuevos orígenes físicos para tales procesos, aunque no proponen implementaciones experimentales específicas más allá de notar que los avances en las plataformas de simulación cuántica (por ejemplo, arreglos de átomos neutros) podrían investigar estos fenómenos.