Fault-tolerant quantum algorithms for simulating atomic nuclei

Resumen Técnico: Algoritmos Cuánticos Tolerantes a Fallos para la Simulación de Núcleos Atómicos

Planteamiento del Problema
Si bien la computación cuántica tolerante a fallos (FTQC) ha experimentado avances significativos en la simulación de la estructura electrónica (química) y sistemas de materia condensada, la simulación de núcleos atómicos ha recibido comparativamente poca atención por parte de la comunidad de computación cuántica. La estructura nuclear presenta un problema de muchos cuerpos fermiónicos distinto que involucra protones y neutrones. A diferencia de los sistemas electrónicos, los núcleos carecen de un fondo iónico, son sistemas auto-ligados gobernados por la fuerza nuclear fuerte, y requieren la inclusión de interacciones de tres cuerpos (y potencialmente de orden superior) derivadas de la teoría de campo efectivo (EFT) quiral. Además, la interacción nuclear está menos comprendida que la interacción de Coulomb, y el escalamiento exponencial del espacio de Hilbert con el número de partículas plantea desafíos computacionales clásicos severos. A pesar de las similitudes con los problemas electrónicos y la importancia crítica de los datos nucleares para la astrofísica y la industria, el trabajo algorítmico cuántico existente se ha centrado mayoritariamente en heurísticas ruidosas de corto plazo en lugar de enfoques escalables y tolerantes a fallos.

Metodología
Los autores desarrollan y compilan algoritmos cuánticos específicos diseñados para hardware tolerante a fallos para simular la estructura nuclear. La metodología central involucra:

1. Hamiltonianos Objetivo: El estudio utiliza dos clases de Hamiltonianos nucleares:

   • Hamiltonianos del Modelo de Capas Fenomenológicos: Estos incluyen interacciones empíricas y fenomenológicas para núcleos específicos ($^{24}\text{Mg}$, $^{32}\text{Mg}$, y $^{219}\text{At}$). Estos modelos típicamente asumen un núcleo inerte (por ejemplo, $^{16}\text{O}$ o $^{208}\text{Pb}$) y generalmente contienen términos de uno y dos cuerpos.
   • Hamiltonianos del Modelo de Capas Sin Núcleo (No-Core Shell-Model): Derivados de la EFT Quiral, estos modelos no asumen un núcleo inerte y están construidos desde principios fundamentales. Incluyen explícitamente términos de interacción de uno, dos y tres cuerpos. Los autores utilizan el software NuHamil para generar estos Hamiltonianos, adecuados para núcleos ligeros hasta aproximadamente $^{40}\text{Ca}$.

2. Marco Algorítmico: Los autores emplean la estimación de fase cuántica (QPE) basada en qubitización. Este enfoque requiere codificar el Hamiltoniano $H$ como una combinación lineal de unitarias (LCU) utilizando dos oráculos primarios:

   • Preparación (Prep): Prepara un estado de superposición ponderado por los coeficientes de los términos del Hamiltoniano.
   • Selección (Sel): Un operador de caminata (walk operator) que selecciona los términos unitarios específicos basados en el estado preparado por Prep.

3. Optimización de Circuitos: Una parte significativa del trabajo se centra en optimizar los circuitos Sel y Prep para Hamiltonianos nucleares, abordando específicamente la transformación de Jordan-Wigner de los operadores fermiónicos.

   • Los autores introducen circuitos optimizados para términos de uno, dos y tres cuerpos que codifican directamente operadores $\sigma^{\pm}$ no unitarios sin expandirlos en cadenas de Pauli, un alejamiento de la técnica previa [71, 85].
   • Utilizan operaciones SwapUp y escaleras CNOT para manejar eficientemente las cadenas de Z de Jordan-Wigner.
   • Para los términos de tres cuerpos, construyen un circuito Sel unificado que codifica por bloques los términos de uno, dos y tres cuerpos simultáneamente, en lugar de tratarlos como componentes LCU separados.
   • Los autores también exploran una estrategia para Hamiltonianos de esquema M (donde los términos de dos cuerpos son bloque-diagonales por la proyección del momento angular total $M$), permitiendo la simulación de bloques $M$ individuales en lugar del Hamiltoniano completo para reducir los requisitos de recursos.

Contribuciones Clave

• Primeras Estimaciones de Recursos: El artículo proporciona las primeras estimaciones conocidas de recursos cuánticos (en términos de conteo de puertas Toffoli y conteo de qubits lógicos) para la simulación de núcleos atómicos tolerante a fallos.
• Construcción de Algoritmos: La construcción de algoritmos de estimación de energía basados en qubitización adaptados específicamente para Hamiltonianos nucleares, incluyendo aquellos con interacciones de tres cuerpos.
• Mejoras de Circuitos: Los autores demuestran que sus construcciones de circuitos logran reducciones significativas en los costos de recursos en comparación con el arte previo. Específicamente, sus circuitos Sel logran una reducción de aproximadamente el 50% en los conteos de Toffoli para términos de un cuerpo, ~40% para términos de dos cuerpos, y una reducción total de ~70% al combinar términos de uno, dos y tres cuerpos en comparación con el enfoque de [71].
• Análisis Estratégico: El trabajo compara dos estrategias de simulación: la simulación directa del Hamiltoniano completo frente a la simulación de bloques $M$ individuales. El análisis muestra que el enfoque de bloques $M$ es más eficiente en recursos con respecto al conteo de qubits y los requisitos de tamaño físico de la computadora cuántica con corrección de errores.

Resultados
Los autores compilaron los algoritmos para casos nucleares específicos y reportaron las siguientes estimaciones de recursos (asumiendo un error de simulación objetivo $\epsilon = 0.01$ MeV):

• Núcleos del Modelo de Capas ($^{32}\text{Mg}$ y $^{219}\text{At}$): Las estimaciones de recursos para los Hamiltonianos de modelo de capas de $^{32}\text{Mg}$ y $^{219}\text{At}$ son comparables a las estimaciones recientes para la simulación del sitio activo de FeMoco (un estándar de referencia en química).
  • Para $^{32}\text{Mg}$, el circuito QPE coherente más grande requiere aproximadamente la mitad del número de qubits, pero un orden de magnitud más de puertas Toffoli que los mejores algoritmos para FeMoco.
  • Para $^{219}\text{At}$, el conteo de qubits es similar al de FeMoco, pero el conteo de Toffoli es, nuevamente, un orden de magnitud superior.
• Modelo de Capas Sin Núcleo ($^{40}\text{Ca}$): Para núcleos ligeros usando Hamiltonianos de EFT Quiral (hasta $\sim^{40}\text{Ca}$), los requisitos de recursos son significativamente más altos que para los casos de modelo de capas. Los autores señalan que la densidad del Hamiltoniano y la inclusión de términos de tres cuerpos hacen que estas simulaciones sean actualmente impracticables sin estrategias más personalizadas.

Significancia y Reivindicaciones
El artículo afirma que estos resultados representan un paso fundacional hacia la realización de simulaciones nucleares útiles en computadoras cuánticas. Al trazar paralelismos entre los problemas nucleares y electrónicos, resaltando al mismo tiempo los desafíos nucleares específicos (como las fuerzas de tres cuerpos y la auto-ligación), los autores pretenden fomentar la colaboración a largo plazo entre las comunidades de física nuclear y computación cuántica.

Los autores enmarcan modestamente su trabajo como una herramienta para generar estimaciones de recursos y demostrar la viabilidad algorítmica, en lugar de resolver preguntas específicas de física nuclear de inmediato. Enfatizan que, si bien los resultados del modelo de capas son prometedores y comparables con los referentes químicos, los resultados del modelo de capas sin núcleo indican que se requiere más innovación algorítmica antes de que estas simulaciones sean practicables en el hardware tolerante a fallos de la próxima generación. El objetivo final declarado es permitir el estudio de los límites de estabilidad nuclear, las formas nucleares y otras preguntas abiertas en la estructura nuclear que son actualmente inaccesibles para la computación clásica.