Variational non-gaussian approach to interacting spin-boson models

Resumen Técnico: Enfoque Variacional No Gaussiano para Modelos de Interacción Espín-Bosón

Planteamiento del Problema
Los modelos de interacción espín-bosón, donde grados de libertad cuánticos discretos (espines) se acoplan a uno o más modos bosónicos (fotones o fonones), son centrales en la ciencia cuántica contemporánea, incluyendo la QED de cavidad, la QED de circuitos y los sistemas espín-fonón. Un desafío primordial en la simulación de estos sistemas surge de la naturaleza de dimensión infinita del espacio de Hilbert bosónico combinada con las fuertes correlaciones de muchos cuerpos en el sector de espín. Los métodos numéricos convencionales enfrentan un compromiso: la truncación de la base de ocupación bosónica introduce sesgo, mientras que relajar el sector de espín hacia descripciones de campo medio falla en capturar el entrelazamiento esencial espín-bosón. Además, los enfoques directos de redes de tensores (como DMRG) aplicados al sistema compuesto completo suelen sufrir un crecimiento rápido en la dimensión de enlace debido a la necesidad de representar explícitamente el campo bosónico junto con los espines.

Metodología
Los autores proponen un marco variacional híbrido que desacopla la optimización del sector bosónico del sector de espín. El estado fundamental se representa dentro de una variedad variacional no gaussiana compacta definida por el ansatz:
$$|\psi_{\text{NGS}}\rangle = U_\lambda (U_{\text{GS}} |0_b\rangle \otimes |\phi_s\rangle)$$
donde:

1. Capa Bosónica Gaussiana ($U_{\text{GS}}$): Una transformación gaussiana que comprende operadores de desplazamiento ($\Delta_R$) y de compresión o squeezing ($\xi$) actuando sobre el vacío bosónico $|0_b\rangle$. Esto captura las fluctuaciones bosónicas y los desplazamientos de campo medio.
2. Transformación de Vestido o Dressing ($U_\lambda$): Un operador unitario que introduce correlaciones espín-bosónicas más allá del nivel de campo medio. Actúa como un desplazamiento dependiente del espín (transformación tipo polaron) controlado por un parámetro variacional $\lambda$. Esta capa mapea el acoplamiento lineal espín-bosón en interacciones efectivas espín-espín.
3. Sector de Espín Efectivo ($|\phi_s\rangle$): El estado de espín no es un simple estado producto, sino el estado fundamental exacto de un Hamiltoniano de espín efectivo ($H_{\text{eff}}$). Este Hamiltoniano se obtiene promediando el Hamiltoniano completo sobre el estado bosónico gaussiano y la transformación de vestido.

El método emplea un bucle de optimización autoconsistente:

• Para parámetros bosónicos fijos ($\Delta_R, \xi, \lambda$), el Hamiltoniano de espín efectivo se resuelve utilizando DMRG (Grupo de Renormalización de la Matriz de Densidad) para obtener $|\phi_s\rangle$ y los valores de expectación de espín.
• Estos observables de espín se utilizan para actualizar los parámetros variacionales bosónicos mediante la minimización del funcional de energía variacional total.
• El bucle itera hasta la convergencia de la energía y los observables seleccionados.

Este enfoque evita la truncación explícita del campo fotónico. En su lugar, los grados de libertad bosónicos se integran variacionalmente, mapeando el problema en un modelo de espín autoconsistente.

Contribuciones Clave y Resultados
El artículo evalúa este marco NGS-DMRG frente a los modelos Dicke y Dicke-Ising, comparando los resultados con simulaciones DMRG de referencia convergidas (que utilizan cortes de fotones explícitos).

1. Precisión y Mejora Variacional:

   • El enfoque no gaussiano produce consistentemente energías de estado fundamental más bajas que el límite gaussiano (campo medio), confirmando la mejora variacional.
   • Para el modelo de Dicke, el método reproduce con precisión la transición de fase cuántica (QPT) superradiante y el escalamiento de tamaño finito, coincidiendo con los resultados de DMRG de referencia con significativamente menos recursos computacionales.
   • En el modelo Dicke-Ising, el marco logra mapear un diagrama de fases complejo que presenta:
     • Transiciones continuas (de segundo orden) entre las fases normal y superradiante.
     • Transiciones discontinuas (de primer orden) y puntos trícriticos en el régimen ferromagnético.
     • Una estrecha fase intermedia de coexistencia de orden antiferromagnético y superradiancia (AFM-SP) en el régimen antiferromagnético, una característica predicha por la teoría de campo medio pero difícil de resolver numéricamente.
   • El método captura las fluctuaciones de tamaño finito y converge al límite termodinámico, donde el entrelazamiento espín-fotón desaparece, siendo consistente con los límites analíticos establecidos.

2. Eficiencia Computacional:

   • Reducción de la Dimensión de Enlace: Al tratar los bosones de forma variacional, el resolvedor DMRG opera sobre un Hamiltoniano efectivo de solo espines. Esto resulta en una dimensión de enlace MPS ($\chi$) significativamente menor en comparación con el DMRG de referencia, que debe manejar explícitamente la cadena compuesta espín-bosón.
   • Tiempo de Ejecución: A pesar de la sobrecarga del bucle autoconsistente, la reducción en la complejidad de la red de tensores conduce a una disminución sustancial en el tiempo total de ejecución. El bucle de optimización externo demuestra una convergencia exponencial.
   • La ventaja es más pronunciada en regímenes con fuertes correlaciones y fronteras de fase complejas (por ejemplo, la fase de coexistencia), donde el DMRG de referencia requiere dimensiones de enlace y cortes de fotones muy grandes.

Significado
El artículo afirma que este marco híbrido ofrece una vía algorítmica distinta respecto a los métodos de redes de tensores de optimización de base local (LBO). Al capturar los modos fotónicos a través de parámetros variacionales (desplazamiento, compresión, vestido) en lugar de una reconfiguración iterativa de la base, el método permite que el resolvedor DMRG se concentre exclusivamente en las correlaciones de espín.

La importancia radica en proporcionar un método numérico robusto y no perturbativo capaz de acceder a la física de espín-bosón fuertemente correlacionada sin eliminar o truncar arbitrariamente el campo bosónico. Los autores demuestran que este enfoque puede reproducir con precisión diagramas de fases complejos y comportamientos críticos en modelos como la cadena Dicke-Ising, logrando una alta fidelidad con un costo computacional reducido. Los autores proporcionan una implementación de código abierto que soporta acoplamientos espín-espín de rango arbitrario y desorden (disponible en: https://github.com/Jpedromend/NGS), con potenciales extensiones a sistemas de múltiples modos y geometrías de mayor dimensión mediante otros resolvedores (por ejemplo, PEPS, QMC).