Semiclassical Langevin dynamics of long-range dissipative time crystals

Resumen Técnico: Dinámica de Langevin Semiclásica de Cristales de Tiempo Disipativos de Largo Alcance

Planteamiento del Problema
Los cristales de tiempo de frontera disipativa (BTC, por sus siglas en inglés) representan una fase de la materia donde la simetría de traslación temporal continua o discreta se rompe espontáneamente en el límite termodinámico, resultando en oscilaciones colectivas persistentes. Mientras que los modelos idealizados con interacciones de rango infinito (totalmente colectivas) soportan un comportamiento de ciclo límite estable descrito por ecuaciones de campo medio, los sistemas reales suelen involucrar interacciones de rango finito y disipación local. Una cuestión crítica abierta es la robustez de la fase del cristal de tiempo cuando la estructura colectiva ideal se debilita. Específicamente, es necesario comprender cómo los efectos de tamaño finito influyen en la vida útil de estas oscilaciones e identificar los umbrales de rango de interacción (parametrizado por un exponente de ley de potencia $\alpha$) más allá de los cuales la fase del cristal de tiempo se disuelve en un estado estacionario. Las simulaciones numéricas exactas utilizando la ecuación maestra de Lindblad están limitadas a tamaños de sistema pequeños, lo que hace necesario otros métodos para investigar el escalamiento de tamaño finito en el límite termodinámico.

Metodología
Los autores desarrollan un enfoque de Langevin semiclasico basado en los marcos de la aproximación de Wigner truncada (TWA) y la TWA discreta (DTWA). En lugar de resolver la ecuación maestra de Lindblad completa, la dinámica se formula mediante ecuaciones de Langevin cuánticas para operadores de espín locales.

• Formalismo: El método reemplaza los productos de operadores por productos de sus valores esperados (aproximación de campo medio) mientras retiene términos de ruido estocástico para dar cuenta de las fluctuaciones cuánticas de primer orden generadas por los canales disipativos.
• Modelo de Espín-1/2: Para una cadena de espín-1/2 con operadores de Lindblad de ley de potencia (la disipación es de largo alcance, el Hamiltoniano es local), la dinámica se mapea a ecuaciones diferenciales estocásticas para los componentes de espín local $s_i^\mu$. El ruido se representa mediante ruidos blancos gaussianos reales.
• Modelo de Espín-1: Para un sistema de espín-1 con disipación local e interacciones de Hamiltoniano de largo alcance, la dinámica se formula utilizando una base de matrices de Gell-Mann (los generadores de $SU(3)$). Las variables locales son los valores esperados de estas matrices, $\lambda_i^\mu$.
• Implementación Numérica: Las ecuaciones estocásticas se resuelven utilizando una interpretación de Itô con un esquema de Runge-Kutta para el arrastre determinista y aumentos de Wiener gaussianos para el ruido. Crucialmente, los autores encuentran que para los modelos estudiados, el muestreo sobre diferentes condiciones iniciales (estándar en DTWA) es innecesario; una sola condición inicial es suficiente porque la memoria del estado inicial se pierde rápidamente en comparación con el periodo de oscilación.
• Diagnóstico: El orden del cristal de tiempo se diagnostica analizando el escalamiento de tamaño finito de la estabilidad de la oscilación. Para el espín-1/2, el enfoque está en la tasa de decaimiento de la envolvente de la oscilación. Para el espín-1, se utilizan dos sondas: el "tiempo de desviación" (cuánto tiempo la trayectoria de tamaño finito permanece cerca del ciclo límite de campo medio) y el escalamiento de la potencia del pico de Fourier dominante.

Resultados Clave

1. Modelo de Espín-1/2 (Disipación de Ley de Potencia):

   • En el límite totalmente colectivo ($\alpha=0$), el sistema exhibe oscilaciones robustas donde el amortiguamiento de tamaño finito se suprime a medida que el tamaño del sistema $L$ aumenta.
   • Para exponentes de ley de potencia finitos $\alpha$, la envolvente de la oscilación decae exponencialmente en el tiempo, pero la tasa de decaimiento $D$ escala algebraicamente con el tamaño del sistema: $D(L) \sim L^{-\beta}$.
   • Un exponente de escalamiento positivo $\beta$ indica que la tasa de decaimiento desaparece en el límite termodinámico, lo que significa un cristal de tiempo estable.
   • El exponente $\beta$ disminuye monótonamente a medida que $\alpha$ aumenta. El régimen de cristal de tiempo persiste para $\alpha \lesssim 1.2$. Notablemente, este umbral está ligeramente por encima de $\alpha=1$. Mientras que la teoría de campo medio predice la relajación hacia un punto fijo para $\alpha > 1$ (donde el campo medio ya no es exacto en el límite termodinámico), el enfoque de Langevin semiclasico predice oscilaciones persistentes en el rango $1 < \alpha \lesssim 1.2$. Esto sugiere que las fluctuaciones cuánticas, capturadas por el método, estabilizan el cristal de tiempo más allá del régimen de campo medio.

2. Modelo de Espín-1 (Disipación Local, Hamiltoniano de Largo Alcance):

   • El sistema muestra comportamiento de cristal de tiempo cuando la interacción coherente es suficientemente de largo alcance.
   • Tiempo de Desviación: Para $\alpha \leq 1$, el tiempo $t^*$ en el que la trayectoria de tamaño finito se desvía de la solución de campo medio escala como $t^*(L) \sim L^b$. El exponente $b$ es positivo para $\alpha$ pequeño (aprox. 0.37 para $\alpha=0$) y disminuye a medida que $\alpha$ se acerca a 1, desapareciendo cerca de $\alpha \approx 1.2$. Este análisis es estrictamente válido solo para $\alpha \leq 1$, donde la teoría de campo medio describe el límite termodinámico.
   • Escalamiento del Pico de Fourier: La potencia del pico de Fourier dominante escala como $P(L) \sim L^s$. El exponente $s$ es positivo para $\alpha \lesssim 1.2$, indicando comportamiento de cristal de tiempo. A diferencia del análisis del tiempo de desviación, esta sonda no depende de la suposición de que el límite termodinico sea de campo medio, lo que la hace robusta para $\alpha > 1$.
   • Ambas sondas identifican un umbral crítico para la desaparición del cristal de tiempo en el intervalo $\alpha \in [1, 1.5]$, consistente con los hallazgos previos de expansión de cumulantes.

Significancia y Reivindicaciones
El artículo afirma que la dinámica de Langevin semiclasica proporciona una sonda de tamaño finito útil para cristales de tiempo disipativos en regímenes donde las simulaciones exactas de Lindblad son computacionalmente intratables.

• Más allá del Campo Medio: El enfoque captura con éxito la robustez del orden de cristal de tiempo en regímenes donde la teoría de campo medio falla (específicamente para $\alpha > 1$ en el modelo de espín-1/2). Demuestra que la fase de cristal de tiempo se extiende más allá del límite estricto de campo medio, siempre que la disipación o la interacción sigan siendo suficientemente de largo alcance.
• Herramientas de Diagnóstico: El trabajo establece los exponentes de escalamiento de tamaño finito (tasa de decaimiento $\beta$, exponente de tiempo de desviación $b$ y exponente de pico de Fourier $s$) como diagnósticos directos de la robustez de la fase de cristal de tiempo.
• Validación Metodológica: Los resultados para el modelo de espín-1 se alinean con hallazgos previos de expansiones de cumulantes, validando el enfoque de Langevin semiclasico como una herramienta confiable para estudiar sistemas cuánticos abiertos con interacciones de largo alcance.
• Perspectiva Física: El estudio destaca que, si bien el límite termodinámico para $\alpha \leq 1$ es descrito por la dinámica de campo medio, el escalamiento de tamaño finito hacia este límite es no trivial y depende de $\alpha$. Además, el escalamiento algebraico del tiempo de desviación (en lugar de logarítmico) sugiere una dinámica regular y no caótica en estas fases de cristal de tiempo.