Perturbative Analysis of Dark State Dynamics in Weakly Anharmonic Photon-Emitter Pairs

Resumen Técnico: Análisis Perturbativo de la Dinámica de Estados Oscuros en Pares de Emisores de Fotones Débilmente Anarmónicos

Enunciado del Problema
Los estados oscuros son estados cuánticos excitados en sistemas abiertos que se desacoplan de su entorno, impidiendo la emisión o absorción de fotones. Aunque estos estados ofrecen potencial para el almacenamiento de información cuántica libre de decoherencia y el entrelazamiento a larga distancia, su estabilidad se ve comprometida cuando el sistema incluye interacciones en el sitio (anarmonicidad). En el contexto de qubits transmon acoplados a una guía de ondas, modelados por el modelo de Bose-Hubbard, la inclusión de anarmonicidad ($U$) rompe el desacoplamiento perfecto de los estados oscuros, induciendo disipación.

El análisis teórico de estos sistemas se complica por su naturaleza no hermítica. El Hamiltoniano efectivo que gobierna la dinámica posee valores propios complejos y requiere el uso de vectores propios izquierdos y derechos biortogonales. A medida que aumenta la anarmonicidad, el sistema se acerca a puntos excepcionales donde los vectores propios coalescen, lo que conduce a inestabilidades numéricas y dificultades para rastrear las poblaciones de los estados, particularmente en los manifiestos de excitación superiores más allá del modelo de qubit. Los métodos existentes a menudo luchan por proporcionar una descripción analítica estable de estas dinámicas sin perder información a través de correcciones de renormalización.

Metodología
Los autores proponen un enfoque perturbativo para analizar el origen de la disipación en los estados oscuros, tratando la anarmonicidad débil como una perturbación a un sistema de dos osciladores armónicos acoplados disipativamente.

1. Formulación del Modelo: El sistema se describe utilizando un Hamiltoniano de Bose-Hubbard no hermítico efectivo. El Hamiltoniano incluye términos armónicos, un término de decaimiento colectivo (acoplamiento disipativo a través de una guía de ondas) y un término de interacción anarmónica ($U$).
2. Expansión Perturbativa: La matriz densidad $\hat{\rho}$ y los vectores de estado $|\psi\rangle$ se expanden en una serie de potencias con respecto a la fuerza de la perturbación $\epsilon$ (que representa la anarmonicidad).
   • La solución de orden cero corresponde al estado oscuro no perturbado en el límite armónico ($U=0$).
   • Se calculan correcciones de primer orden tanto para la energía como para la función de onda utilizando teoría de perturbaciones no degenerada adaptada para sistemas no hermíticos.
   • Se derivan correcciones de segundo orden, incorporando los efectos de los operadores de salto y la mezcla de estados.
3. Seguimiento de la Dinámica: Las matrices densidad corregidas (hasta segundo orden) se aplican a la ecuación maestra para rastrear la evolución dependiente del tiempo del sistema.
4. Validación: Los resultados perturbativos se comparan con la diagonalización numérica exacta del Hamiltoniano no hermítico y la evolución numérica estándar de la ecuación maestra. Para manejar la naturaleza biortogonal de los vectores propios cerca de los puntos excepcionales, los autores emplean un método de ortogonalización de Gram-Schmidt simultánea para vectores propios izquierdos y derechos durante la inicialización.

Contribuciones y Resultados Clave

• Mecanismo de Disipación: El análisis revela que la anarmonicidad induce disipación en los estados oscuros principalmente a través de un acoplamiento inducido entre el estado oscuro y un estado disipativo vecino (brillante) que comparte paridad colectiva similar.
• Correcciones Analíticas:
  • Primer Orden: La corrección de primer orden a la función de onda introduce una correlación de valor complejo entre el estado oscuro y el estado brillante. Esto crea una vía de transición que permite que la población se filtre desde el estado oscuro hacia canales disipativos.
  • Segundo Orden: La corrección de segundo orden a la energía revela un desplazamiento en la componente imaginaria, cuantificando directamente la tasa de disipación inducida. La corrección escala con el número de excitación $N$ como $\lambda_c(N) \propto -i N(N-1)U^2/16\gamma$.
• Comportamiento Dinámico:
  • El método perturbativo reproduce con éxito la dinámica de relajación numérica del estado oscuro hacia el estado fundamental con alta precisión para anarmonicidad débil ($U/\gamma < 0.5$).
  • La dinámica exhibe un estallido inicial de emisión de fotones (firma de tipo superradiante) antes de asentarse en un decaimiento casi exponencial. Este estallido surge porque la población se mezcla inicialmente en el estado brillante (que es máximamente disipativo) más rápido de lo que el estado decae.
  • Dependencia de la Paridad: El estado final de relajación depende de la paridad del número de excitación inicial. Los números de excitación pares se relajan al estado fundamental, mientras que los números de excitación impares se relajan y permanecen atrapados en el estado oscuro $N=1$, ya que el primer manifiesto de excitación no se ve afectado por la anarmonicidad en este modelo.
• Estabilidad: El enfoque perturbativo proporciona un método numéricamente estable y que preserva la traza para rastrear la dinámica de poblaciones, evitando las inestabilidades asociadas con la coalescencia de funciones de onda biortogonales cerca de los puntos excepcionales.

Significado y Afirmaciones
El artículo afirma que este marco perturbativo ofrece una alternativa robusta a la diagonalización numérica directa para estudiar sistemas cuánticos abiertos no hermíticos con anarmonicidad débil. Al proporcionar expresiones analíticas para las correcciones de la función de onda y la energía, el método permite la construcción de un "mapa" de vías de transición y correlaciones entre estados subradiantes y superradiantes.

Los autores declaran que este enfoque es particularmente útil para:

1. Comprender la ruptura de espacios libres de decoherencia en manifiestos de excitación superiores (más allá del límite de qubit).
2. Predecir la dinámica de poblaciones y los tiempos de coherencia para estados entrelazados en arreglos de transmon.
3. Informar esquemas de control cuántico en circuitos superconductores identificando cómo la anarmonicidad afecta la fidelidad del estado.

El trabajo es modesto en su alcance, reconociendo que las correcciones analíticas se desvían de los resultados numéricos a medida que la anarmonicidad aumenta hacia y más allá del punto excepcional, lo cual es consistente con las limitaciones de la teoría de perturbaciones. Los autores sugieren que el trabajo futuro podría involucrar escalar el sistema a arreglos más grandes e investigar otros factores de desfase, pero el artículo actual se centra estrictamente en validar el método perturbativo frente a técnicas numéricas establecidas para un sistema de dos transmon.