Steady-state phase transition in one-dimensional quantum contact process

Planteamiento del Problema
El artículo investiga las transiciones de fase en estado estacionario del proceso de contacto cuántico (QCP) unidimensional, un modelo que describe sistemas de muchos cuerpos fuera del equilibrio cuántico. Mientras que el proceso de contacto clásico es un modelo bien establecido para fenómenos como la propagación de enfermedades y la difusión de información, su extensión cuántica introduce una compleja interacción entre la evolución temporal coherente y la disipación. Una controversia central en el campo concierne a la naturaleza de la transición de fase en el QCP unidimensional: específicamente, si la transición del fase absorbente (sitios vacíos) al fase activo (sitios ocupados) es continua o discontinua. Estudios previos han arrojado resultados contradictorios, con algunos sugiriendo transiciones continuas y otros insinuando discontinuidades, dejando sin resolver el orden de la transición en sistemas unidimensionales.

Metodología
Los autores emplean un enfoque teórico multifacético para resolver la naturaleza de la transición:

1. Análisis del Espectro de Liouvilliano: La evolución temporal del sistema se describe mediante una ecuación maestra de Lindblad. Los autores analizan la brecha de Liouvillian (la magnitud de la parte real negativa máxima de los autovalores de Liouvilliano) en el límite termodinámico para identificar transiciones de fase y metaestabilidad.
2. Aproximaciones de Campo Medio (Mean-Field):
   • Campo Medio de Sitio Único (MF): Utilizado como base para identificar soluciones de estado estacionario y comportamientos de bifurcación.
   • Campo Medio de Clúster (CMF): Para incorporar correlaciones de corto alcance, los autores utilizan el método CMF. Una innovación clave es la introducción de una nueva condición de autoconsistencia basada en un campo efectivo. Esta condición parametriza la interacción entre clústeres y permite a los autores capturar todas las posibles soluciones de estado estacionario, incluyendo las inestables, mientras evita que el sistema quede atrapado en estados metaestables durante la evolución temporal numérica.
3. Expansión de Clústeres Vinculados Numérica (LCE): Para investigar las correlaciones de largo alcance y el comportamiento de la susceptibilidad magnética ($\chi$) en el límite termodinámico, los autores implementan una LCE numérica. Este método suma las contribuciones de clústeres de tamaño creciente para extrapolar los resultados a un sistema infinito.

Resultados Clave

• Metaestabilidad y Brecha de Liouvilliano: El análisis del espectro de Liouvilliano revela el cierre de la brecha en $\Omega/\Gamma \approx 5.83$ en el límite termodinámico. La estructura del espectro indica la existencia de un estado metaestable cerca del punto de transición, caracterizado por una meseta de larga duración en la evolución temporal de la población antes de decaer al verdadero estado absorbente. Esto explica por qué las simulaciones previas podrían haber fallado al observar el verdadero estado estacionario si el tiempo de simulación era insuficiente.
• Transición de Fase Discontinua: Tanto el análisis de MF de sitio único como el de CMF revelan una bifurcación de silla-nodo del parámetro de orden (población en estado estacionario $\langle \bar{n} \rangle_{ss}$). A medida que varía el parámetro de control $\Omega$, el sistema exhibe una coexistencia de estados absorbentes y activos, donde la rama activa emerge vía un punto de silla-nodo ($S_1$) en lugar de una bifurcación de bifurcación (pitchfork) supercrítica. Este comportamiento es característico de una transición de fase discontinua (de primer orden).
• Extracción del Punto de Transición: Utilizando el método CMF con la nueva condición de autoconsistencia, los autores extraen el punto de transición para sistemas de tamaño infinito. Al extrapolar la dependencia del tamaño del clúster del punto de transición $\Omega_c(L)$, determinan que $\Omega_c/\Gamma \approx 6.04$, lo cual concuerda bien con la predicción de la brecha de Liouvilliano de $\approx 5.83$.
• Longitud de Correlación y Susceptibilidad: Los resultados de la LCE muestran que en la fase absorbente, los pesos de los clústeres decaen exponencialmente, indicando correlaciones de corto alcance. En la fase activa, el decaimiento sigue una ley de potencia. Crucialmente, la susceptibilidad magnética $\chi$ extrapolada exhibe un decrecimiento monotónico a través del punto de transición sin divergir. Esta falta de divergencia en la susceptibilidad y la longitud de correlación respalda adicionalmente la conclusión de que la transición es discontinua, distinguiéndola de las transiciones continuas donde tales cantidades típicamente divergen.

Significancia y Reivindicaciones
El artículo afirma resolver la controversia sobre la naturaleza de la transición de fase en el QCP unidimensional al demostrar que es una transición discontinua del fase absorbente al fase activo. Los autores enfatizan que la observación de esta discontinuidad depende de manejar correctamente la metaestabilidad del sistema, que puede atrapar las simulaciones numéricas en falsos estados estacionarios. La condición de autoconsistencia propuesta basada en campos efectivos se presenta como una herramienta robusta para sortear estas interferencias metaestables y capturar el panorama completo de las soluciones de estado estacionario.

Los resultados sugieren que las propiedades del estado estacionario están dominadas por correlaciones de corto alcance, validando el uso de aproximaciones de campo medio de clúster para este sistema. Los autores señalan que estos hallazgos teóricos podrían ser probados en simuladores cuánticos utilizando átomos de Rydberg. Concluyen identificando el estudio de procesos de contacto híbridos (que combinan ramificación/coagulación cuántica y clásica) y la extracción de exponentes críticos en tales sistemas como direcciones futuras prometedoras.