Feynman's clock and hierarchy-informed sampling for quantum error mitigation

Resumen Técnico: El Reloj de Feynman y el Muestreo Informado por la Jerarquía para la Mitigación de Errores Cuánticos

Planteamiento del Problema
Los algoritmos cuánticos de la era actual (NISQ) enfrentan limitaciones significativas debido a la naturaleza ruidosa del hardware, lo que perjudica la ejecución física de los circuitos cuánticos. Si bien la corrección de errores requiere dispositivos de gran escala que aún no están disponibles, es necesario contar con esquemas eficientes de mitigación de errores cuánticos (QEM) para aplicaciones de corto plazo. Los métodos existentes están fundamentalmente limitados, pero siguen siendo el medio principal para reducir el ruido cuántico. Este artículo aborda la necesidad de una técnica de QEM general capaz de mitigar el ruido derivado de la ejecución de circuitos cuánticos arbitrarios, extendiéndose más allá de las aplicaciones previas limitadas a simulaciones específicas de cadenas de espines.

Metodología
La técnica propuesta, una extensión del esquema de Mitigación por Muestreo Informado por BBGKY (BBGKY-ISM), opera a través de un proceso de tres etapas:

1. Mapeo de Circuito-a-Hamiltoniano mediante el Reloj de Feynman:
   Los autores mapean un circuito arbitrario de $N_Q$ qubits, compuesto por $N_G$ puertas, sobre la evolución temporal de un sistema cuántico gobernado por un Hamiltoniano independiente del tiempo. Esto utiliza la construcción de la computadora cuántica de Feynman, introduciendo un registro de reloj $H_C$ con $N_C = \lceil \log_2(N_G + 1) \rceil$ qubits auxiliares. El sistema evoluciona bajo el Hamiltoniano del reloj de Feynman $H$, donde los estados del reloj rastrean el progreso de la computación del circuito. La dinámica completa está descrita por $|\psi(t)\rangle = e^{-itH}|\psi(0)\rangle$.

2. Reconstrucción mediante Pruebas de Hadamard:
   En lugar de simular la evolución temporal mediante Trotterización (lo que amplificaría el ruido), el método reconstruye los valores de la esperanza temporal de interés utilizando pruebas de Hadamard sobre ejecuciones parciales del circuito original.

• Las mediciones Z en el registro de datos se reinterpretan como valores de la esperanza dependientes del tiempo de operadores Pauli-Z individuales, proyectados sobre el estado final del reloj $|N_G\rangle_C$.
• Estas cantidades se descomponen en un conjunto de "cantidades de interés" que involucran cadenas de Pauli.
• Las partes real e imaginaria de los elementos de matriz necesarios se miden mediante pruebas de Hadamard, evitando la profundidad de circuito profunda asociada con las rebanadas de Trotter.
• Los valores finales de la medición Z se extraen sumando estas cantidades de evolución temporal reconstruidas, ponderadas por la amplitud de probabilidad del estado del reloj, eliminando efectivamente la dependencia temporal.

3. Mitigación por Muestreo Informado por BBGKY:
   Las estimaciones ruidosas obtenidas de las pruebas de Hadamard se mitigan utilizando el esquema BBGKY-ISM.

• Distribución Informada por la Física: El método muestrea candidatos de mitigación de una distribución de probabilidad codificada con la dinámica sin ruido. Esta distribución está definida por una acción $S(\vec{x})$ que se minimiza cuando la configuración $\vec{x}$ aproxima la evolución ideal y sin ruido.
• Restricciones Jerárquicas: Las restricciones se derivan de una jerarquía de ecuaciones de movimiento de tipo BBGKY para los correladores de cadenas de Pauli. La jerarquía conecta las cantidades de interés con correladores de orden superior.
• Recocido Simulado (Simulated Annealing): Se utiliza un procedimiento de Monte Carlo de Cadenas de Markov (Recocido Simulado) para muestrear configuraciones que minimicen la acción. Un parámetro $r$ (radio) determina cuántas capas de la jerarquía se incluyen en las restricciones. Un $r$ mayor codifica una mayor parte de la dinámica ideal, estrechando las restricciones sobre el muestreo.

Contribuciones Clave

• Generalización de BBGLY-ISM: El artículo extiende el esquema BBGKY-ISM desde simulaciones específicas de cadenas de espines hacia circuitos cuánticos arbitrarios mediante el uso de la construcción del reloj de Feynman.
• Escalabilidad de Sobrecarga Eficiente: Los autores demuestran que tanto la sobrecarga clásica como la cuántica de este método escalan polinomialmente con el tamaño del circuito ($N_G$), el número de qubits ($N_Q$) y el número de pasos de tiempo ($N_T$). Específicamente, el número de pruebas de Hadamard requeridas escala como $N_Q \cdot \text{poly}(N_G)$.
• Naturaleza de Post-procesamiento: La técnica es enteramente un método de post-procesamiento, lo que permite acoplarla con esquemas de mitigación existentes o aplicarla después de la adquisición de datos sin modificar la ejecución del hardware cuántico.
• Reducción de Error Controlable: El método introduce un parámetro ajustable $r$ (el radio de la jerarquía) que permite a los usuarios controlar la intensidad de la mitigación.

Resultos
El método fue probado numéricamente en circuitos de preparación de estados Bell ajustables ($N_Q=2, N_G=2$) bajo un modelo de ruido que imita el procesador IBM Fez.

• Reducción Sistemática: Los resultados muestran una reducción sistemática del error cuántico ($\Delta' < \Delta$) en todos los parámetros de circuito y radios de jerarquía probados.
• Capacidad de Ajuste: La capacidad de reducción de error aumenta a medida que el radio de la jerarquía $r$ aumenta. Esto confirma que codificar una mayor parte de la jerarquía BBGKY restringe el muestreo de manera más efectiva hacia la dinámica sin ruido.
• Autoconsistencia: En un radio específico ($r=4$ en el caso probado), el sistema de ecuaciones se vuelve completamente autoconsistente ($z=1$). En este punto, el método se reduce a un cálculo numérico puramente clásico que recupera la dinámica exacta sin ruido, ignorando efectivamente la contribución del ruido cuántico.

Significancia y Reivindicaciones
El artículo afirma proporcionar una nueva técnica de mitigación de errores cuánticos de propósito general que es tanto eficiente (escalado polinomial) como controlable. Al mapear los circuitos a la dinámica de Hamiltoniano y aprovechar el muestreo informado por la física, el método ofrece una forma sistemática de reducir el ruido sin requerir recursos cuánticos adicionales más allá de la ejecución inicial del circuito y el post-procesamiento estándar. Los autores enfatizan que la efectividad de la técnica es ajustable mediante el radio de la jerarquía, lo que la hace adaptable a diferentes niveles de ruido y restricciones computacionales. El trabajo sugiere posibles aplicaciones futuras en algoritmos más complejos (por ejemplo, Grover, Shor) y ansatze variacionales más profundos, señalando su compatibilidad como una capa de post-procesamiento con otras estrategias de mitigación.