Emergence of caustics in dynamics of the Kitaev model

Resumen Técnico: Emergencia de Caústicas en la Dinámica del Modelo de Kitaev

Planteamiento del Problema
El artículo investiga la dinámica de no equilibrio de sistemas cuánticos de muchos cuerpos cerrados, centrándose específicamente en la propagación de cuasipartículas en sistemas integrables bidimensionales (2D). Si bien las estructuras de cono de luz, acotadas por el límite de Lieb-Robinson, están bien establecidas en sistemas unidimensionales (1D) y sirven como firmas de caústicas cuánticas (análogas a singularidades en óptica geométrica), su comportamiento en modelos integrables 2D permanece menos explorado. Los autores tienen como objetivo identificar estas estructuras de cono de luz y caústicas cuánticas en el modelo de panal de Kitaev 2D, analizar su anisotropía espacial y determinar cómo los impulsos periódicos externos (ingeniería de Floquet) alteran estas dinámicas.

Metodología
El estudio emplea una combinación de marcos analíticos y simulaciones numéricas:

1. Mapeo del Modelo: El modelo de panal de Kitaev de espín-1/2 se mapea a un modelo fermiónico no interactuante utilizando la transformación de Jordan-Wigner. El sistema se analiza dentro del sector donde la constante de movimiento $\hat{D}_n = 1$.
2. Dinámica de Cuasipartículas: La dinámica de una sola cuasipartícula iniciada en un sitio específico se describe mediante una función de onda $\psi(r, t)$ expresada como una suma sobre modos de momento en la zona de Brillouin. La fase de esta función de onda está definida por un funcional de acción $S(k; r, t)$.
3. Identificación de Caústicas: Las caústicas cuánticas se identifican como los lugares de densidad infinita de rayos donde la amplitud de la función de onda se maximiza. Matemáticamente, estos corresponden a la coalescencia de puntos críticos (puntos de silla) del funcional de acción. Los autores resuelven las condiciones $\nabla_k S = 0$ y $\det(\nabla_k \nabla_k S) = 0$ para derivar las trayectorias exactas de estas caústicas.
4. Conducción Periódica: Para estudiar efectos de no equilibrio, los autores introducen protocolos dependientes del tiempo donde el parámetro de acoplamiento $J_3$ se modula periódicamente. Se analizan dos protocolos de impulso específicos:
   • Patadas de función $\delta$ aplicadas periódicamente.
   • Pulsos de onda cuadrada que alternan entre dos fuerzas de acoplamiento.
     La evolución se trata utilizando la teoría de Floquet, construyendo un Hamiltoniano efectivo de Floquet ($H_{eff}$) a partir del operador de evolución temporal sobre un periodo.

Contribuciones y Resultados Clave

1. Anisotropía Espacial en 2D: A diferencia de los sistemas 1D, la dinámica de cuasipartículas en el modelo de Kitaev 2D exhibe una fuerte anisotropía espacial. La amplitud de la función de onda se propaga a lo largo de direcciones específicas en el plano de la red, determinadas por la relación de los parámetros de acoplamiento ($J_2/J_1$). A medida que varía esta relación, cambia la dirección de máxima propagación, y la amplitud decae a medida que uno se aleja de este eje preferente.
2. Trayectorias Exactas de Caústicas y Límite de Lieb-Robinson: Al resolver las condiciones de coalescencia para puntos críticos, los autores derivan expresiones analíticas exactas para los envolventes de las caústicas. Estos envolventes definen la velocidad máxima de propagación de la información, recuperando efectivamente el límite de Lieb-Robinson para el sistema 2D.
   • En la fase con brecha ($J_3 \gg J_1, J_2$), los límites del cono de luz son lineales: $x_1 = 2J_1t$ y $x_2 = 2J_2t$.
   • En la fase sin brecha ($J_3 = J_1 = J_2$), los límites también forman una estructura de cono de luz donde el momento $k$ permanece real solo dentro de límites específicos.
3. Efecto de la Conducción Periódica: La introducción de un impulso dependiente del tiempo altera fundamentalmente la estructura de las caústicas.
   • Bajo ambos protocolos de patada $\delta$ y onda cuadrada, la propagación característica en forma de cono de luz observada en el sistema estático desaparece.
   • En su lugar, la amplitud de la función de onda exhibe una propagación local que decae con cada periodo de conducción.
   • La cuasipartícula queda confinada (localizada) a lo largo de líneas específicas (por ejemplo, $x_1 = 0$) en lugar de propagarse libremente. Esta localización persiste por más tiempo a medida que la frecuencia del impulso se acerca a cero ($\omega \to 0$).

Significado y Afirmaciones
El artículo afirma proporcionar el primer marco analítico para estudiar la dinámica de una sola cuasipartícula en el modelo de Kitaev 2D utilizando funciones de catástrofe. El significado principal radica en:

• Demostrar que los sistemas integrables 2D soportan estructuras similares a conos de luz que sirven como caústicas cuánticas, pero con una dependencia direccional (anisotropía) ausente en 1D.
• Proporcionar soluciones analíticas exactas para los envolventes de las caústicas, que corresponden al límite de Lieb-Robinson en 2D.
• Mostrar que los impulsos periódicos externos pueden suprimir completamente la estructura del cono de luz, conduciendo a la localización de cuasipartículas. Esto sugiere que la frecuencia del impulso es un parámetro de control crítico para manipular la propagación de información cuántica y las escalas de tiempo de medición en tales sistemas.

Los autores señalan que, si bien la realización del modelo de Kitaev en redes ópticas es un tema de investigación en curso, sus resultados teóricos ofrecen una base para rastrear la propagación de cuasipartículas y potencialmente controlar el movimiento de cuasipartículas individuales mediante campos externos en futuros montajes experimentales.