Real-time Krylov Diagonalisation for Open Quantum Systems

Resumen Técnico: Diagonalización de Krylov en Tiempo Real para Sistemas Cuánticos Abiertos

Enunciado del Problema
El artículo aborda el desafío de simular sistemas cuánticos abiertos descritos por ecuaciones maestras de Lindblad, centrándose específicamente en la estimación del hueco de Liouvilliano. Si bien la Diagonalización del Subespacio Cuántico (QSD) ha surgido como un algoritmo líder no tolerante a fallos para sistemas cerrados (hamiltonianos), su aplicación a sistemas disipativos no es trivial. Los sistemas abiertos están gobernados por superoperadores de Liouvilliano no hermíticos y no normales, lo que introduce diferencias estructurales como vectores propios bi-ortogonales, valores propios complejos y dinámicas no unitarias. Estas características complican los métodos estándar de subespacio de Krylov, particularmente en lo que respecta a la estabilidad numérica (mal condicionamiento) y la interpretación de filtros espectrales. Los autores buscan adaptar la QSD en tiempo real para estimar el hueco de Liouvilliano, un parámetro crítico que define el tiempo de relajación y la presencia de modos lentos en transiciones de fase disipativas.

Metodología
Los autores proponen una modificación de la Diagonalización del Subespacio Cuántico (QSD) en tiempo real para manejar dinámicas de Liouvilliano. La metodología central implica:

1. Construcción del Subespacio de Krylov: En lugar de evolucionar un estado bajo un hamiltoniano unitario $H$, el sistema evoluciona bajo el superoperador de Liouvilliano $\mathcal{L}$. Un subespacio de Krylov se genera aplicando potencias del operador de evolución temporal $e^{\mathcal{L}\tau}$ a un estado inicial $|\rho_0)$:
   $$\text{span}\{|\rho_0), e^{\mathcal{L}\tau}|\rho_0), \dots, e^{\mathcal{L}(D-1)\tau}|\rho_0)\}$$
2. Problema Generalizado de Valores Propios (GEVP): Dado que los circuitos cuánticos no pueden realizar reortogonalización sobre la marcha, el método resuelve un GEVP $\bar{H}v = \lambda Sv$. Para los Liouvillianos, la matriz de solapamiento $S$ y la matriz hamiltoniana $\bar{H}$ se construyen utilizando vectores propios izquierdos y derechos bi-ortogonales.
3. Interpretación del Filtro Espectral: El método se enmarca como la optimización de un filtro espectral. A diferencia de los sistemas hamiltonianos donde los estados propios son reales y la evolución temporal es oscilatoria (requiriendo un muestreo estricto de Nyquist), las dinámicas de Liouvilliano son disipativas. Los modos de decaimiento rápido giran hacia el interior en el plano complejo, permitiendo tiempos de evolución más largos $\tau$ sin artefactos de aliasing. Esto permite el aislamiento de modos lentos (valores propios con partes reales pequeñas) incluso con una dimensión de Krylov reducida $D$.
4. Manejo de la No Normalidad: Los autores reconocen que los operadores no normales conducen a una alta sensibilidad a las perturbaciones (teorema de Bauer-Ficke). La matriz de solapamiento $S$ se vuelve exponencialmente mal condicionada a medida que aumenta $D$. Para mitigar esto, el método se basa en la naturaleza disipativa del sistema para suprimir los modos rápidos, permitiendo un $\tau$ mayor y un $D$ menor, reduciendo así la inestabilidad numérica asociada con grandes dimensiones de Krylov.
5. Selección del Estado Inicial: Crucialmente, el estado inicial debe tener una proyección no nula sobre el sector de traza cero (coherencias o desequilibrios de población) para excitar los modos lentos. El sector de traza 1 (estado estacionario) no contribuye a la señal para la estimación del hueco en este marco.

Contribuciones Clave

• Extensión de QSD a Liouvillianos: El artículo proporciona un marco teórico para aplicar la QSD en tiempo real a superoperadores de Liouvilliano no hermíticos y no normales en forma de Lindblad.
• Estrategia de Estimación del Hueco: Demuestra cómo estimar el hueco de Liouvilliano ($\Delta = -\max \Re(\lambda_k)$) explotando la naturaleza disipativa del sistema para filtrar los modos de decaimiento rápido, eludiendo efectivamente algunas restricciones de Nyquist que se aplican a la evolución unitaria.
• Análisis de Estabilidad Numérica: El trabajo destaca la compensación entre el tiempo de evolución $\tau$ y la dimensión de Krylov $D$. Argumenta que aumentar $\tau$ (en lugar de $D$) es una estrategia viable para sistemas abiertos para mejorar la resolución del hueco manteniendo la estabilidad numérica, ya que los modos rápidos decaen naturalmente.
• Aplicación a Qubits de Gato de Kerr: El método se aplica a un resonador de Kerr superconductor impulsado por dos fotones (Qubit de Gato de Kerr), un sistema relevante para arquitecturas de qubits protegidos.

Resultados
Los autores aplicaron su método a un modelo simplificado de un Qubit de Gato de Kerr con pérdida de fotones individuales.

• Parámetros de Simulación: Las simulaciones se realizaron con un máximo de 30 fotones y una dimensión de Krylov $D=20$, resultando en una matriz de Krylov de $900 \times 20$.
• Reconstrucción del Hueco: El método reconstruyó con éxito el hueco de Liouvilliano a través de diferentes intensidades de impulsión. Los resultados mostraron que aumentar el tiempo de evolución $\tau$ mejoró la precisión de la estimación del hueco, particularmente en regímenes donde el conmutador de Liouvilliano (una medida de la no normalidad) era alto.
• Identificación de Modos Lentos: Las simulaciones confirmaron la aparición de un modo lento correspondiente al decaimiento de coherencias y desequilibrios de población. En el régimen de "ruido sesgado", el hueco se contrae exponencialmente con la amplitud del estado coherente ($\propto e^{-\alpha^2}$), indicando la formación de un subspace computacional protegido.
• Representación de Wigner: Los estados propios reconstruidos de la matriz de Krylov se utilizaron para visualizar la representación de Wigner del modo lento, mostrando la transición de un régimen amortiguado a un régimen con dos lóbulos distintos (ruptura de simetría) y de vuelta a un régimen comprimido a alta impulsión.

Significado y Afirmaciones
El artículo afirma demostrar la viabilidad de utilizar la QSD en tiempo real para estudiar sistemas cuánticos abiertos, específicamente para estimar el hueco de Liouvilliano en el contexto de transiciones de fase disipativas. Los autores posicionan este trabajo como una adaptación necesaria de algoritmos de la era NISQ para el diseño y la simulación de hardware cuántico, que intrínsecamente implica disipación.

El significado se enmarca en el potencial de utilizar computadoras cuánticas para simular modos lentos de sistemas cuánticos que son computacionalmente costosos de simular clásicamente, particularmente en el contexto de qubits protegidos como el Gato de Kerr. Los autores notan modestamente que, aunque los Qubits de Gato de Kerr individuales son simulables clásicamente, las redes de tales qubits podrían no serlo, sugiriendo la utilidad del método para sistemas más grandes y complejos. El trabajo se presenta como un paso hacia el "descubrimiento automatizado de hardware cuántico" al proporcionar herramientas para caracterizar las propiedades de ruido y la estabilidad de las arquitecturas de qubits propuestas. El artículo no afirma haber resuelto el problema general de la simulación de operadores no normales, sino que ofrece un enfoque específico y modificado para la estimación del hueco en sistemas disipativos.