Low-depth simulation of non-Markovianity under quantum hardware noise

Resumen Técnico: Simulación de baja profundidad de no-markovianidad bajo ruido de hardware cuántico

Planteamiento del Problema
Simular sistemas cuánticos abiertos en computadoras cuánticas digitales es intrínsecamente desafiante porque los circuitos cuánticos soportan nativamente solo evoluciones unitarias, mientras que los sistemas abiertos experimentan dinámicas no unitarias debido a las interacciones con el entorno. El enfoque estándar, la dilatación de Stinespring, simula estas dinámicas acoplando el sistema a qubits auxiliares y eliminándolos mediante la traza parcial. Sin embargo, este método requiere topologías de interacción densas y un alto número de puertas de dos qubits (por ejemplo, CNOTs), que son significativamente más ruidosas que las puertas de un solo qubit. En el hardware de la era actual (near-term), la acumulación de errores de estas puertas auxiliares limita severamente la fidelidad y el tiempo de ejecución de las simulaciones, particularmente para evoluciones de largo tiempo o dinámicas no-markovianas complejas.

Metodología
El artículo propone aprovechar el método de mezcla de trayectorias (mezcla estadística de trayectorias de estados puros) para simular la dinámica de sistemas abiertos sin necesidad de qubits auxiliares. Este enfoque se aplica específicamente a canales unitarios mixtos, donde un canal cuántico $\Phi(\rho)$ puede expresarse como una mezcla estadística de evoluciones unitarias: $\Phi(\rho) = \sum_i p_i U_i \rho U_i^\dagger$.

Características metodológicas clave:

1. Probabilidades independientes del estado: Para los canales unitarios mixtos, la probabilidad de un salto unitario específico ($p_i$) es independiente del estado del sistema. Esto permite que la simulación proceda mediante la ejecución de múltiples circuitos independientes (trayectorias) y la mezcla estadística de sus resultados, evitando la necesidad de una costosa tomografía de estado para determinar dinámicamente las probabilidades de salto.
2. Muestreo de trayectorias: Aunque el número total de trayectorias posibles escala exponencialmente ($m^N$, donde $m$ es el número de operadores de Kraus y $N$ es el número de pasos de tiempo), los autores emplean un muestreo estocástico (enfoque Monte Carlo) para simular evoluciones de largo tiempo con un subconjunto de trayectorias estadísticamente suficiente.
3. Modelado de ruido: Las simulaciones utilizan un modelo de ruido realista calibrado para procesadores cuánticos modernos, incluyendo relajación/desfasaje $T_1$ y $T_2$ después de puertas de un solo qubit, y un canal de depolarización después de las puertas CNOT.
4. Línea base de comparación: El método de mezcla de trayectorias se compara con el enfoque tradicional de qubits auxiliares (dilatación de Stinespring) a través de varios escenarios.

Contribuciones Clave y Simulaciones
El artículo demuestra la viabilidad de este enfoque a través de cuatro escenarios de simulación específicos:

• Desfasaje Puro (Pure Dephasing): Un solo qubit que experimenta evolución libre y desfasaje puro. El método de mezcla de trayectorias (usando solo puertas de un solo qubit) mantiene una alta fidelidad de estado y sigue de cerca la dinámica sin ruido. En contraste, el método de qubits auxiliares (que requiere CNOTs) muestra una deriva significativa hacia el centro de la esfera de Bloch debido al ruido de las puertas.
• Dinámica No-Markoviana (Amortiguamiento de Amplitud y Desfasaje): Los autores simulan evoluciones no-markovianas mediante incrustaciones markovianas (Markovian embeddings), donde un sistema interactúa con un qubit de memoria.
  • Resultado: El método de mezcla de trayectorias ofrece consistentemente una mayor fidelidad que el método de qubits auxiliares, particularmente cuando se considera el estado conjunto sistema-memoria. Sin embargo, el artículo señala que para el sistema por sí solo, la ventaja de fidelidad es a veces marginal. Los autores sugieren que, debido a que el sistema y la memoria interactúan constantemente, el ruido de las puertas de dos qubits en el modelo de incrustación limita la ventaja del método de trayectoria en escenarios de muchos cuerpos complejos.
• Canales de Memoria: La simulación de un canal de desfasaje de memoria que actúa sobre dos qubits enviados en rápida sucesión. Este canal preserva el entrelazamiento en estados Bell.
  • Resultado: El método de mezcla de trayectorias preserva la concurrencia (entrelazamiento) casi perfectamente, mientras que el método de qubits auxiliares provoca un decaimiento significativo. Esto resalta una ventaja importante: el método de trayectoria evita las dos puertas CNOT requeridas para correlacionar el qubit auxiliar con ambos qubits del sistema.
• Hamiltonianos Impulsados Estocásticamente: Los autores simulan dinámicas no-markovianas impulsadas por un Hamiltoniano con Ruido de Telegrafista Aleatorio (RTN). Al descomponer la evolución en intervalos de ruido constante y mezclar las trayectorias, logran fidelidades extremadamente altas sin qubits auxiliares ni qubits de memoria.

Resultos
Utilizando un modelo de ruido calibrado para procesadores actuales ($T_1=232\,\mu s$, $T_2=127\,\mu s$, error de CNOT $\approx 2.4\%$), el estudio encuentra que:

• Fidelidad: El método de trayectoria mixta retiene consistentemente una mayor fidelidad de estado en comparación con el método de qubits auxiliares en los escenarios de desfasaje puro y canales de memoria.
• Preservación del Entrelazamiento: En las simulaciones de canales de memoria, el método de trayectoria logra preservar el entrelazamiento, mientras que el método auxiliar falla en hacerlo debido a la acumulación de errores de las puertas de dos qubits.
• Escalabilidad: Incluso con un muestreo limitado (1024 trayectorias para 25 pasos de tiempo), el método de trayectoria supera al método auxiliar, el cual sufre un decaimiento constante de la fidelidad a medida que aumenta el número de pasos de tiempo.

Significancia y Reivindicaciones
El artículo afirma que el método de mezcla de trayectorias ofrece un plano de implementación (blueprint) desplegable para simular efectos ambientales no locales en hardware cuántico de corto plazo. Su principal significancia radica en:

1. Reducción de Ruido: Al intercambiar puertas de entrelazamiento por una mezcla estadística de circuitos independientes de estado puro, el método reduce drásticamente la sobrecarga de ruido asociada con los qubits auxiliares.
2. Eficiencia de Recursos: Permite la simulación de evoluciones no-markovianas de largo tiempo con recursos limitados, evitando la "devastadora sobrecarga de ruido" de las puertas auxiliares.
3. Practicidad: Aunque el método requiere ejecutar múltiples circuitos, los autores argumentan que este es un intercambio favorable dado los actuales niveles de error de las puertas de dos qubits en los procesadores cuánticos.

Los autores mantienen la modestia respecto a la universalidad del método, señalando que para dinámicas de muchos cuerpos complejas donde las interacciones sistema-memoria son constantes, la ventaja puede ser menos pronunciada debido al ruido inherente de las interacciones de dos qubits requeridas en los modelos de incrustación. Sin embargo, para canales y escenarios de memoria específicos, el enfoque proporciona un camino claro hacia simulaciones de mayor fidelidad en el hardware existente.