PT-symmetry breaking phase transitions in an LMG dimer

Resumen Técnico: Transiciones de Fase de Ruptura de Simetría PT en un Dímero LMG

Planteamiento del Problema
Los sistemas cuánticos abiertos impulsados tanto por campos externos como por la disipación exhiben una fenomenología más rica que sus contrapartes en equilibrio, incluyendo la estabilización de fases imposibles en los estados fundamentales. Mientras que investigaciones previas han explorado las transiciones de fase disipativas en modelos como el modelo de Dicke abierto y el modelo Lipkin-Meshkov-Glick (LMG) disipativo, e investigado por separado las transiciones de ruptura de simetría Paridad-Tiempo (PT), la interacción entre estas distintas clases de transiciones no equilibrium permanece poco explorada. Específicamente, existe la necesidad de comprender cómo se comporta un sistema que combina un mecanismo de acoplamiento de tipo ruptura de simetría PT con una interacción LMG no lineal. Los autores pretenden construir un modelo que integre estos dos mecanismos para investigar la competencia resultante, la estructura de fases y los regímenes dinámicos emergentes.

Metodología
Los autores introducen un modelo consistente en un par de modelos LMG acoplados: uno disipativo y otro impulsado, vinculados mediante un Hamiltoniano con simetría PT. El sistema está descrito por una ecuación maestra que gobierna la matriz de densidad $\rho$ de dos espines colectivos (etiquetados como A y B) de longitud $S$.

• Construcción del Modelo: El Hamiltoniano total $H$ combina los Hamiltonianos LMG para cada espín ($H_A^{LMG} \otimes I_B + I_A \otimes H_B^{LMG}$) con un término de acoplamiento PT-simétrico ($H_{PT}$). La dinámica incluye superoperadores de Lindblad que representan ganancia en un espín y pérdida en el otro.
• Análisis de Campo Medio: Para analizar el sistema en el límite termodinámico ($S \to \infty$), los autores derivan un conjunto de seis ecuaciones diferenciales no lineales acopladas para los componentes de espín normalizados ($X, Y, Z$). Estas ecuaciones son exactas en el límite termodinámico para esta clase de modelos.
• Investigación Numérica: Los autores resuelven numéricamente las ecuaciones de campo medio para mapear el diagrama de fases en estado estacionario, identificando puntos fijos, ciclos límite y regiones caóticas. La estabilidad se determina mediante los autovalores de la matriz Jacobiana.
• Caracterización Dinámica: Para distinguir entre la dinámica regular y la caótica, los autores computan el exponente de Lyapunov máximo promediado sobre condiciones iniciales simétricas.
• Simulación Cuántica Exacta: Para validar las predicciones de campo medio, los autores realizan la diagonalización numérica exacta de la ecuación maestra para tamaños de espín finitos ($S$), trabajando en un espacio de Liouville de tamaño $(2S+1)^4$. También computan funciones Wigner de espín para comparar las representaciones de fase del estado cuántico con las trayectorias de campo medio.

Contribuciones Clave y Resultados

1. Diagrama de Fases Complejo: La competencia entre la transición de simetría PT (impulsada por el acoplamiento inter-espín $\lambda$) y la transición ferromagnética LMG (impulsada por la no linealidad $J$) genera un diagrama de fases con múltiples fases de estado estacionario que no están presentes en los modelos constituyentes.

   • Fase Normal: A valores bajos de $\lambda$ y $|J|$, el sistema exhibe un único punto fijo estable correspondiente a una magnetización escalonada (dominada por la disipación).
   • Fases LMG: A medida que $|J|$ aumenta, el sistema entra en regiones con dos puntos fijos estables, correspondientes a estados ferromagnéticos con ruptura de simetría (ordenamiento de espín-A o espín-B).
   • Régimen PT/Caótico: Para un $\lambda$ grande, el sistema pierde puntos fijos estables. La teoría de campo medio predice la emergencia de ciclos límite y, en regiones de parámetros específicos, dinámica caótica caracterizada por exponentes de Lyapunov máximos positivos.

2. Regímenes Dinámicos Emergentes: El artículo identifica que la interacción de las dos transiciones crea comportamientos dinámicos ausentes en los modelos individuales. Notablemente, existen regiones donde el análisis de campo medio predice puntos fijos estables, pero el exponente de Lyapunov es positivo, lo que indica una dinámica oscilatoria compleja dependiente de las condiciones iniciales (cuencas de atracción).

3. Validación mediante Escalamiento de Tamaño Finito: Los autores demuestran que las características predichas por la teoría de campo medio sobreviven en el modelo cuántico completo.

   • Los resultados de la diagonalización exacta para $S$ finito coinciden cualitativamente con las fronteras de fase de campo medio.
   • Los efectos de tamaño finito suavizan las transiciones agudas (por ejemplo, la transición PT en $\lambda = g/2$) y aseguran un operador de densidad de estado estacionario único incluso cuando la teoría de campo medio predice ciclos límite.
   • Las funciones Wigner de espín confirman que las distribuciones de probabilidad cuántica se alinean con los puntos fijos de campo medio (en la fase LMG) o con distribuciones difusas correspondientes a ciclos límite (en la fase PT).

Significancia y Reivindicaciones
El artículo afirma demostrar que el acoplamiento de sistemas con distintos mecanismos de ruptura de simetría disipativa puede generar fases dinámicas emergentes, incluyendo el caos y ciclos límite complejos, que están ausentes en cualquiera de los modelos constituyentes. Los autores proponen este trabajo como un marco general para diseñar fases en sistemas impulsados-disipativos mediante el acoplamiento de modelos con fases conocidas a través de interacciones que soportan distintas transiciones no equilibrium.

Los autores aseguran que su modelo es accesible utilizando las tecnologías experimentales actuales, citando realizaciones existentes de modelos tipo LMG en qudits superconductores, Dysprosium ultrafrío, QED de cavidad y QED de circuitos, combinadas con desarrollos en disipación diseñada. Sugieren que tales técnicas podrían ser relevantes para tecnologías cuánticas que requieren estados estacionarios diseñados o puntos críticos, como la computación cuántica disipativa, la detección cuántica y la termodinámica cuántica. El trabajo destaca el potencial para descubrir nuevas clases de universalidad y comportamientos de escalamiento de tamaño finito, aunque los autores señalan que determinar estas propiedades requeriría enfoques numéricos más especializados.