Full Configuration Interaction Quantum Monte Carlo for Accurate $\textit{Ab Initio}$ Nuclear Structure Calculations

Resumen Técnico: Interacción de Configuración Completa mediante Monte Carlo de Funciones de Onda para Cálculos Estructurales Nucleares Ab Initio Precisos

Planteamiento del Problema
El desarrollo de métodos de muchos cuerpos ab initio para núcleos fuertemente interactuantes es un desafío central en la teoría nuclear moderna. Si bien la teoría de campo efectivo quiral (EFT) proporciona un marco sistemático para derivar interacciones nucleares a partir de la cromodinámica cuántica, resolver la ecuación de Schrödinger de muchos cuerpos para núcleos de masa media y pesada sigue siendo computacionalmente prohibitivo.

• Modelo de Capas Sin Núcleo (NCSM): Aunque es exitoso para núcleos ligeros, el NCSM está limitado por el crecimiento exponencial de la dimensión del espacio de modelo, que rápidamente supera la capacidad computacional de las supercomputadoras actuales (alcanzando $\sim 10^{12}$ dimensiones). En consecuencia, el NCSM depende en gran medida de extrapolaciones para sistemas más grandes.
• Expansiones de Muchos Cuerpos Truncadas: Métodos como la Teoría de Perturbaciones de Muchos Cuerpos (MBPT), el Grupo de Renormalización de Similitud In-Medio (IMSRG) y la teoría de Clústeres Acoplados (CC) ofrecen un escalamiento polinomial, permitiendo cálculos para núcleos más pesados. Sin embargo, sufren de errores de truncamiento. Lograr una alta precisión requiere llevar estas expansiones a órdenes superiores, lo que aumenta la complejidad matemática y el costo computacional. Además, proporcionar estimaciones fiables de las incertidumbres de muchos cuerpos en estos esquemas truncados es difícil, particularmente para sistemas con fuerte agrupamiento (clustering) o correlaciones colectivas donde la convergencia puede ser lenta.

Existe una necesidad crítica de un método que retenga la precisión de la Interacción de Configuración Completa (FCI) sin el almacenamiento explícito del Hamiltoniano, sirviendo como un referente para validar y guiar el desarrollo de métodos de expansión truncados.

Metodología
Los autores introducen el método de Interacción de Configuración Completa mediante Monte Carlo de Funciones de Onda (FCIQMC) al dominio de la estructura nuclear ab initio. El FCIQMC es un método de proyector estocástico que resuelve la ecuación de Schrödinger en tiempo imaginario:
$$-\frac{d}{d\tau} |\Psi(\tau)\rangle = (\hat{H} - E_0) |\Psi(\tau)\rangle$$
Las características clave de la implementación incluyen:

• Muestreo Estocástico: La función de onda se representa mediante un conjunto de "caminantes" (walkers, partículas con signo) distribuidos a través de determinantes de Slater en el espacio de configuración completo. El coeficiente de un determinante está determinado por la población neta de signos de los caminantes en dicho determinante.
• Pasos de Propagación: El conjunto de caminantes evoluciona a través de tres pasos estocásticos por cada paso de tiempo imaginario ($\Delta\tau$):
  1. Generación (Spawning): Muestrea los acoplamientos de la Hamiltoniana fuera de la diagonal ($H_{ij}$) para crear nuevos caminantes en determinantes conectados.
  2. Muerte/Clonación (Death/Cloning): Muestrea el término diagonal ($H_{ii}$) desplazado por un parámetro $S$, ajustando la población de caminantes en el determinante actual.
  3. Aniquilación: Cancela los caminantes de signos opuestos que residen en el mismo determinante, lo cual es crucial para superar el problema del signo de Fermión sin requerir una aproximación de nodo fijo o una función de onda guía.
• Aproximación de Iniciador: Para controlar el problema del signo en espacios de modelo grandes, los autores emplean la aproximación de iniciador combinada con una corrección de desplazamiento adaptativo. Esto permite estimaciones controlables y mejorables del límite FCI a medida que aumenta la población de caminantes ($N_w$).
• Configuración del Sistema: Los cálculos se realizaron para $^4$He, $^8$Be, $^{12}$C y $^{16}$O utilizando una interacción de dos nucleones quiral optimizada (N2LOopt) dentro de una base de oscilador armónico esférico ($\hbar\omega = 20$ MeV). No se aplicaron truncamientos de configuración en los cálculos de FCIQMC.

Resultos Clave
El estudio presenta energías del estado fundamental ($E$) y radios de carga ($R_{ch}$) para los núcleos seleccionados con incertidumbres de muchos cuerpos a nivel de sub-percentaje.

• Validación en Espacios Pequeños: En un espacio de modelo pequeño ($e_{max}=2$), los resultados de FCIQMC coincidieron con los resultados exactos de FCI obtenidos mediante diagonalización directa, validando la precisión del muestreo estocástico.
• Desempeño en Grandes Espacios de Modelo: En un espacio de modelo grande ($e_{max}=10$), donde la dimensión del espacio de Hilbert para $^{16}$O alcanza $\sim 10^{33}$ (haciendo imposible la diagonalización directa), el FCIQMC proporcionó resultados precisos.
  • Para $^4$He, los resultados de FCIQMC concordaron perfectamente con las extrapolaciones casi exactas de NCSM.
  • Para núcleos más pesados ($^8$Be, $^{12}$C, $^{16}$O), donde el NCSM falla en converger, el FCIQMC entregó resultados precisos.
• Evaluación de Métodos Truncados:
  • MBPT: Mostró las mayores desviaciones respecto a los resultados de FCIQMC.
  • IMSRG: Generalmente razonable, pero exhibió desviaciones significativas para $^8$Be, donde las correcciones de orden superior (IMSRG(3f2)+T) se desviaron aún más de los valores exactos, indicando la ausencia de correlaciones colectivas.
  • Clústeres Acoplados (CC): Mostró una mejora sistemática desde CCSD hasta CCSDT-3. Sin embargo, para el radio de carga de $^{16}$O, la inclusión de los triples iterativos principales condujo a una mayor desviación en comparación con FCIQMC, resaltando la dificultad de la convergencia en ciertos observables.
• Efectos de Agrupamiento (Clustering): Los resultados para $^8$Be, que posee una estructura molecular/de agrupamiento pronunciada, demostraron que las expansiones basadas en referencia (como IMSRG de una sola referencia) luchan por converger debido a la omisión de correlaciones colectivas de alto orden. El FCIQMC, al incluir correlaciones a todos los rangos de excitación dentro de la base, mantuvo la precisión para estos estados agrupados.

Significancia y Reivindicaciones
El artículo afirma que el FCIQMC establece un nuevo estándar para los cálculos de estructura nuclear ab initio:

1. Resolutor de Espacio Completo Estocástico: Es capaz de tratar sistemas más allá del alcance del NCSM convencional mediante el muestreo directo de la función de onda exacta en el espacio de configuración completo.
2. Incertidumbres Cuantificadas: A diferencia de los métodos de expansión truncada donde la estimación de la incertidumbre es difícil, el FCIQMC proporciona resultados con incertidumbres de muchos cuerpos cuantificadas y de sub-percentaje.
3. Referencia para la Teoría: Sus resultados sirven como referencias exactas para evaluar la precisión de las interacciones quirales y la convergencia de los métodos de escalamiento polinomial (MBPT, IMSRG, CC). El estudio demuestra explícitamente el papel crucial de las correlaciones de alto orden que a menudo faltan en los esquemas truncados.
4. Impacto Futuro: Los autores postulan que este método facilitará el desarrollo de la próxima generación de métodos de muchos cuerpos con mayor precisión y ayudará a resolver enigmas de larga data en la teoría nuclear, como la evolución de los radios de carga en la cadena de calcio y las correcciones de ruptura de isospín en decaimientos Fermi superpermitidos, donde los efectos de muchos cuerpos de alto orden son centrales.

El trabajo representa la primera aplicación de FCIQMC a cálculos ab initio precisos de núcleos finitos en interacciones fuertes, marcando un paso significativo hacia la comprensión de las fuerzas nucleares y las simetrías fundamentales con una precisión sin precedentes.