Nuclear scattering phase shifts with factorized geometric-time RODEO on a quantum processor

Resumen Técnico: Desfases de Dispersión Nuclear con RODEO de Tiempo Geométrico Factorizado en un Procesador Cuántico

Planteamiento del Problema
La reconstrucción de los desfases de dispersión nuclear a partir de espectros computados cuánticamente representa un desafío significativo. Si bien los algoritmos de resolución de problemas cuánticos (quantum eigensolvers) pueden estimar niveles de energía discretos en un sistema confinado (una "trampa"), conectar estos estados discretos con observables de espacio libre (dispersión en espacio libre) requiere un post-procesamiento clásico, específicamente la Expansión de Rango Efectivo Modificada (MERE, por sus siglas en la literatura). El rendimiento de un circuito cuántico en este contexto no puede juzgarse únicamente por el error de raíz cuadrática media (RMS) de los niveles de energía individuales. Los errores espectrales pueden ser amplificados, suprimidos o reorganizados durante la subsiguiente interpolación de energía y la extrapolación de confinamiento cero. Por lo tanto, un circuito cuántico con reducción de recursos debe evaluarse mediante su capacidad para preservar la información necesaria para la reconstrucción final de la dispersión, en lugar de basarse solo en la fidelidad espectral.

Metodología
Los autores reconstruyen los desfases de dispersión elástica s-wave neutrón-protón (np) utilizando datos de la procesadora cuántica IBM Aachen. El estudio compara dos implementaciones del algoritmo RODEO (Repetitive Operation for Discrete Energy Optimization):

1. Directo Dinámico R10: Una implementación base que utiliza diez ciclos con pasos de tiempo de cuantiles gaussianos fijos. Para cada energía de consulta, la evolución controlada se resintetiza y se compila de forma independiente.
2. FG–R6 de Tiempo Geométrico Factorizado: Una implementación comprimida que utiliza seis ciclos con pasos de tiempo muestreados geométricamente. Este método emplea una separación exacta de la fase de consulta, donde la evolución del Hamiltoniano se fija a una energía de referencia ($E_{ref}$) y se compila una sola vez. La fase dependiente de la consulta ($\Delta E = E_q - E_{ref}$) se vincula al qubit ancilla tras la compilación.

Ambos métodos utilizan un registro de tres qubits (dos qubits de sistema y un ancilla reciclado) en el backend IBM Aachen. El modelo físico es una interacción de pozo cuadrado s-wave np esquemática. Se utilizan Hamiltonianos efectivos de Okubo–Lee–Suzuki (OLS) de cuatro dimensiones, construidos clásicamente, para representar el sistema. El experimento mide cuatro niveles de energía positiva en cinco intensidades de trampa diferentes ($\omega = 3.6$ a $4.0$ MeV), lo que genera veinte entradas espectrales. Estas entradas se introducen en un análisis MERE clásico para extrapolar al espacio libre y calcular los desfases sobre una rejilla de energía de 0.1 a 30.0 MeV.

Contribuciones Clave

• Compresión de Circuitos mediante Factorización: El artículo introduce el circuito FG–R6, que combina el reúso de ancillas con una estrategia de evolución temporal factorizada. Al separar la fase de consulta de la evolución del Hamiltoniano, el circuito reduce el número de ciclos de evolución controlada de diez a seis y elimina la necesidad de recompilar el operador de evolución para cada energía de consulta.
• Evaluación de Hardware (Benchmarking): El estudio proporciona una evaluación completa derivada de hardware para la dispersión nuclear, yendo más allá de las pruebas de un solo nivel hacia un conjunto completo de veinte entradas espectrales requeridas para la reconstrucción en espacio libre.
• Análisis de la Relación Recurso vs. Exactitud: El trabajo demuestra que reducir la profundidad del circuito y el conteo de puertas de dos qubits no necesariamente degrada la exactitud del observable de dispersión final, incluso si el error RMS de los niveles de energía intermedios aumenta.

Resultados

• Reducción de Recursos: En comparación con la implementación Directa Dinámica R10, la implementación FG–R6 reduce la profundidad mediana del circuito compilado en un 33.1% (de 459 a 307) y el conteo de puertas de dos qubits en un 37.2% (de 113 a 71).
• Incertidumbre de Energía: El FG–R6 logra una incertidumbre de energía de disparo finito mediana menor (0.534 keV) en comparación con el R10 (0.610 keV), a pesar de tener una amplitud ajustada mediana mayor y un desplazamiento de centro RMS ($\Delta E_c$) mayor respecto a la energía de referencia (1.648 keV frente a 0.793 keV).
• Exactitud del Desfase de Dispersión: En la rejilla completa de 0.1–30.0 MeV, la desviación máxima del desfase de dispersión reconstruido respecto a la referencia de energía exacta de MERE disminuye significativamente de 0.842° (R10) a 0.163° (FG–R6). La desviación máxima respecto a la solución analítica del pozo cuadrado también cae de 0.981° a 0.156°.
• Rendimiento Dependiente de la Energía: La mejora es impulsada principalmente por la región de extrapolación de baja energía (por debajo de 10 MeV). En el subconjunto de interpolación común de 10.0–30.0 MeV, la implementación Directa R10 produce en realidad residuos ligeramente menores que el FG–R6.
• Sensibilidad de la Forma del Ajuste: El estudio observa que la sensibilidad a la forma del ajuste sigue siendo apreciable. La dispersión de los resultados entre los diferentes órdenes polinómicos para los ajustes de energía y trampa es sustancial (por ejemplo, un rango de 20.12° en la rejilla completa para R10 frente a 6.75° para FG–R6), lo que indica que la elección del modelo de extrapolación impacta significativamente el resultado final.

Significado y Reivindicaciones
El artículo afirma que esta evaluación de espacio reducido logra conectar con éxito la compresión de circuitos RODEO con los observables nucleares de continuo. La importancia principal radica en demostrar que el rendimiento del circuito debe evaluarse a través de la reconstrucción del observable de dispersión y su rango de energía específico, en lugar de basarse únicamente en los errores RMS espectrales. La implementación FG–R6 se presenta como un compromiso de exactitud-recurso validado para el espectro y la ventana específicos probados, en lugar de ser una solución de mínimo ciclo universal.

Los autores declaran explícitamente que estos resultados no demuestran ventaja cuántica ni escalabilidad de muchos cuerpos. Los datos de hardware no determinan por separado los efectos de la compilación de fase frente al muestreo temporal, ni establecen la reproducibilidad de ejecuciones independientes. Este trabajo sirve como punto de partida para extender estos métodos a modelos de mayor tamaño, interacciones realistas y la dinámica de reacciones de pocos cuerpos, tales como la dispersión resonante de ondas parciales superiores (por ejemplo, la dispersión neutrón–alfa) y sistemas que involucran interacciones de Coulomb de largo alcance.