QwaveMPS: An efficient open-source Python package for simulating non-Markovian waveguide-QED using matrix product states

A continuación presento un resumen técnico detallado del artículo "QwaveMPS: An efficient open-source Python package for simulating non-Markovian waveguide-QED using matrix product states", traducido y estructurado en español.

1. El Problema

La Electrodinámica Cuántica en Guías de Ondas (Waveguide-QED) es una plataforma fundamental para estudiar interacciones luz-materia en sistemas cuasi-unidimensionales, donde emisores cuánticos (átomos, transmons, puntos cuánticos) se acoplan a un continuo de modos de campo cuantizado.

Sin embargo, modelar estos sistemas presenta desafíos significativos:

• Crecimiento exponencial del espacio de Hilbert: La combinación de no linealidades fuertes y retardos temporales (feedback) hace que la simulación exacta sea computacionalmente prohibitiva para métodos tradicionales.
• Limitaciones de las aproximaciones actuales: La mayoría de los enfoques existentes utilizan la aproximación de Markov, que ignora los efectos de retroalimentación retardada (time-delayed feedback) y asume acoplamientos instantáneos. Esto falla drásticamente en regímenes no markovianos.
• Restricciones de otros métodos: Técnicas como la teoría de scattering o métodos estocásticos cuánticos suelen estar limitados a excitaciones débiles o bajos números de fotones, lo que impide estudiar sistemas fuertemente no lineales o con múltiples cuantos.
• Falta de herramientas accesibles: Aunque existen métodos basados en Redes de Tensores (Tensor Networks), carecían de paquetes de software de código abierto y fácil uso específicos para la comunidad de óptica cuántica, a diferencia de herramientas como QuTiP (basado en ecuaciones maestras).

2. Metodología

El artículo introduce QwaveMPS, una biblioteca de Python de código abierto que utiliza el formalismo de Estados de Producto Matricial (MPS) para simular sistemas de Waveguide-QED.

• Discretización Temporal: En lugar de tratar el campo electromagnético como un continuo, el método discretiza la guía de ondas en "bines" de tiempo. Esto permite representar el campo como una cadena de espines o modos discretos.
• Redes de Tensores (TN): Se utiliza la descomposición en valores singulares (SVD) para descomponer el estado cuántico en un producto de tensores locales. Esto evita el crecimiento exponencial del espacio de Hilbert al capturar solo la entropía de entrelazamiento relevante.
• Operadores de Producto Matricial (MPO): La evolución temporal se realiza aplicando operadores de evolución (MPO) sobre los tensores del MPS.
  • En el régimen Markoviano, el operador actúa sobre el bin del sistema y el bin de tiempo actual.
  • En el régimen No Markoviano (con feedback), el operador se convierte en un MPO de 3 sitios que incluye el bin del sistema, el bin de tiempo actual y el bin de retroalimentación (retardo), tras aplicar operadores de intercambio (swap) para reordenar los tensores.
• Tratamiento Unificado: El paquete trata a los átomos cuantizados y a los fotones cuantizados en pie de igualdad, permitiendo simulaciones exactas sin aproximaciones de campo medio.

3. Contribuciones Clave del Paquete QwaveMPS

El paquete ofrece una interfaz amigable para construir, evolucionar y analizar estados cuánticos en sistemas 1D abiertos. Sus características principales incluyen:

• Versatilidad de Regímenes: Capacidad para simular tanto regímenes Markovianos como No Markovianos (con retardos de tiempo y espejos laterales).
• No Linealidad y Múltiples Cuantos: A diferencia de métodos lineales, QwaveMPS puede manejar sistemas fuertemente no lineales y estados con múltiples fotones (ej. estados de Fock con $n > 1$).
• Tipos de Excitación: Soporta:
  • Dinámica de vacío.
  • Bombeo clásico (ondas continuas CW y pulsos temporales).
  • Excitación cuántica directa mediante pulsos de estados de Fock (1, 2 o más fotones).
• Cálculo de Observables: Incluye funciones para calcular:
  • Dinámicas de población (átomos y flujo de fotones).
  • Funciones de correlación cuántica de dos tiempos ($G^{(1)}$, $G^{(2)}$).
  • Espectros de emisión (triplete de Mollow).
  • Entropía de entrelazamiento (sistema-campo y circuitos completos).
• Eficiencia Computacional: Diseñado para ejecutarse en estaciones de trabajo estándar e incluso portátiles, utilizando recursos mínimos de RAM (< 0.2 GB en ejemplos complejos).

4. Resultados y Ejemplos Demostrativos

Los autores validaron el paquete mediante una serie de ejemplos comparados con soluciones analíticas o solvers de ecuaciones diferenciales (ODE):

1. Decaimiento en Régimen Lineal:
   • Simulación de la desintegración de un TLS (sistema de dos niveles) en guías infinitas y semi-infinitas.
   • Demostración de efectos de feedback destructivo y constructivo, mostrando la conservación de cuantos y la dinámica de poblaciones con tiempos de ejecución de menos de 1 segundo.
2. Dinámica No Lineal (Dos Cuantos):
   • Estudio de dos TLSs inicialmente excitados. Se observó el atrapamiento de fotones en el bucle de retroalimentación y se calculó la entropía de entrelazamiento, mostrando desviaciones significativas respecto al límite Markoviano.
3. Bombeo Clásico (CW y Pulsos):
   • Análisis de un TLS bajo bombeo continuo, recuperando el triplete de Mollow y mostrando cómo el feedback no markoviano modifica las correlaciones y el espectro.
   • Simulación de pulsos gaussianos, revelando el atrapamiento de población y patrones de transmisión en correlaciones de segundo orden.
4. Pulsos de Estados de Fock:
   • Inyección directa de pulsos de 1 y 2 fotones. Se observó la forma característica de "pájaro" en la función de correlación de segundo orden para el caso de 2 fotones, un resultado difícil de obtener con otros métodos.

Rendimiento: Todos los ejemplos se ejecutaron en tiempos del orden de segundos (o menos de un minuto para correlaciones de dos tiempos complejas), demostrando una eficiencia superior a los métodos tradicionales de espacio de Hilbert completo.

5. Significado e Impacto

El trabajo de QwaveMPS es significativo por varias razones:

• Democratización de Métodos Avanzados: Lleva la potencia de las Redes de Tensores (MPS) a la comunidad de óptica cuántica, proporcionando una herramienta accesible que antes requería conocimientos profundos de teoría de tensores.
• Precisión en Regímenes No Markovianos: Permite estudiar fenómenos físicos reales (como la retroalimentación retardada en circuitos superconductores o guías de ondas nanofotónicas) que las aproximaciones de Markov no pueden capturar.
• Escalabilidad: Al evitar el crecimiento exponencial del espacio de Hilbert, permite simular sistemas con múltiples emisores y múltiples fotones, abriendo la puerta al estudio de correlaciones cuánticas complejas y dinámicas no lineales fuertes.
• Reproducibilidad: Al ser de código abierto, fomenta la reproducibilidad de estudios en Waveguide-QED y facilita la comparación de resultados entre diferentes grupos de investigación.

En conclusión, QwaveMPS cierra la brecha entre técnicas teóricas potentes y estudios prácticos, ofreciendo una solución eficiente y exacta para la próxima generación de dispositivos cuánticos basados en guías de ondas.