Full configuration interaction quantum Monte Carlo for accurate $\textit{ab initio}$ nuclear structure calculations: algorithms and calculation details

Resumen Técnico: Interacción de Configuración Completa mediante Monte Carlo Cuántico para Cálculos de Estructura Nuclear Ab Initio Precisos

Planteamiento del Problema
La teoría nuclear ab initio busca describir núcleos finitos a partir de Hamiltonianos microscópicos derivados de la teoría de campo efectivo quiral. Mientras que los métodos de expansión de muchos cuerpos (p. ej., Coupled CC, IMSRG, MBPT) extienden los cálculos a núcleos más pesados, estos dependen de truncamientos que introducen incertidumbres de muchos cuerpos residuales difíciles de cuantificar. Por el contrario, el Modelo de Capas Sin Núcleo (NCSM) proporciona acceso directo a la solución de la interacción de configuración completa (FCI), pero está limitado por el crecimiento exponencial del espacio de Hilbert. Existe una necesidad crítica de resolvedores estocásticos capaces de aproximarse al límite FCI en espacios de modelos grandes sin almacenar explícitamente la matriz del Hamiltoniano de muchos cuerpos, particularmente para sistemas con fuertes correlaciones de orden superior derivadas de la deformación, el agrupamiento (clustering) o estructuras de capa abierta.

Metodología
Este trabajo adapta el algoritmo de Interacción de Configuración Completa mediante Monte Carlo Cuántico (FCIQMC), desarrollado originalmente para la química cuántica, a los cálculos de estructura nuclear. El método emplea un enfoque de proyector estocástico en tiempo imaginario para resolver la ecuación de Schrödinger dentro de un sector de simetría elegido (número de partículas, paridad y proyección de momento angular).

Los componentes algorítmicos clave descritos incluyen:

• Propagación en Tiempo Imaginario: La función de onda se representa mediante una población de caminantes (walkers) distribuidos en determinantes de Slater. La propagación implica pasos de generación (spawning), muerte/clonación y aniquilación para establecer estocásticamente la estructura de signo fermiónica.
• Generación de Excitaciones: Para mejorar la eficiencia sobre la generación uniforme, los autores utilizan un generador de baño térmico (heat-bath) precalculado que muestrea preferencialmente excitaciones con elementos de matriz del Hamiltoniano grandes.
• Aproximación de Iniciador (i-FCIQMC): Para superar el problema de signo en grandes espacios de Hilbert, se restringe la generación desde determinantes débilmente ocupados. Esto introduce un sesgo de caminante finito que se reduce sistemáticamente a medida que aumenta la población de caminantes.
• Corrección de Desplazamiento Adaptativo (AS): Se asigna un desplazamiento local a los determinantes no iniciadores para acelerar la convergencia hacia el límite FCI, preservando al mismo tiempo el límite exacto a medida que el número de caminantes tiende a infinito.
• Estimadores de Observables:
  • Estimadores Proyectados/de Prueba: Utilizados para energías, empleando determinantes de referencia o funciones de onda de prueba compactas.
  • Estimadores Puros: Esenciales para observables no conmutantes como los radios. Estos se calculan mediante el muestreo de la Matriz de Densidad Reducida (RDM) utilizando dos réplicas de caminantes independientes para eliminar el sesgo inherente a los estimadores cuadráticos.
• Estados Excitados: El método se extiende a estados excitados mediante la propagación multi-estado ortogonalizada, donde múltiples conjuntos de caminantes se propagan en paralelo y se ortogonalizan en cada paso.

Contribuciones Clave

1. Adaptación del Formalismo: El artículo detalla la implementación específica de FCIQMC para Hamiltonianos nucleares, incluyendo el manejo de las fuerzas de tres nucleones y los esquemas específicos de generación de excitaciones requeridos para interacciones nucleares.
2. Validación contra FCI Determinista: En espacios de modelos pequeños ($e_{\text{max}}=2$) para $^4$He, $^8$Be, $^{12}$C y $^{16}$O, los resultados estocásticos de FCIQMC para energías y radios de protones al cuadrado ($R_p^2$) muestran que reproducen los valores exactos de FCI determinista dentro de las incertidumbres estadísticas.
3. Control de Error Sistemático: Los autores demuestran que el sesgo residual de caminante finito en las aproximaciones de iniciador puede controlarse y cuantificarse. Introducen una estrategia utilizando extrapolaciones de ley de potencia ($O(N_w) = O_\infty + \lambda N_w^{-\alpha}$) en grandes espacios de modelos ($e_{\text{max}}=10$). El exponente $\alpha$ está constreñido por requerir la consistencia entre la energía proyectada de convergencia rápida y la energía pura de convergencia más lenta.
4. Extensión a Estados Excitados: El trabajo demuestra la capacidad de calcular espectros de baja energía (específicamente para $^6$Li) resolviendo múltiples estados dentro del mismo sector de simetría sin diagonalización explícita.

Resultados

• Benchmarking: En el espacio $e_{\text{max}}=2$, las energías y radios de FCIQMC coinciden con los referentes de FCI. La corrección de desplazamiento adaptativo reduce significamente el sesgo de caminante finito y debilita la dependencia del umbral del iniciador ($n_{\text{init}}$).
• Grandes Espacios de Modelos ($e_{\text{max}}=10$):
  • Para $^4$He, el cálculo converge directamente con el aumento del número de caminantes ($N_w$), coincidiendo con los resultados de NCSM.
  • Para $^8$Be, $^{12}$C y $^{16}$O, la energía proyectada converge rápidamente, mientras que los estimadores puros basados en RDM (energía y radios) exhiben una convergencia más lenta. Al aplicar la extrapolación de ley de potencia restringida por la energía, los autores obtienen resultados que capturan correlaciones más allá de los métodos de expansión de muchos cuerpos truncados (MBPT, CC, IMSRG).
  • La estrategia de extrapolación se valida mostrando que los resultados son independientes del parámetro de desplazamiento de la corrección adaptativa ($\tilde{\Delta}$) y consistentes a través de diferentes umbrales de iniciador.
• Estados Excitados: Para $^6$Li, el método resuelve con éxito los cinco estados más bajos ($1^+, 3^+, 0^+, 2^+_1, 2^+_2$). El espectro de excitación está correctamente ordenado incluso con números de caminantes bajos, convergiendo más rápido que los resultados de NCSM con respecto al truncamiento de muchos cuerpos ($N_{\text{max}}$).

Significancia
El artículo posiciona a FCIQMC como un resolvedor estocástico sistemáticamente mejorable capaz de aproximarse al límite FCI para sistemas nucleares donde los métodos deterministas fallan debido al tamaño del espacio de Hilbert. Cumple dos funciones principales:

1. Resolvedor Preciso: Proporciona referencias de alta precisión para energías de estado fundamental y radios en grandes espacios de modelos, capturando altas correlaciones de orden que los métodos de expansión truncada omiten.
2. Herramienta de Referencia (Benchmark): Ofrece un medio para evaluar las incertidumbres residuales de muchos cuerpos de otros métodos de expansión de muchos cuerpos truncados al proporcionar una referencia cercana a la solución exacta.

El trabajo establece la base algorítmica y numérica para los resultados de física reportados en una Carta compañera [80], confirmando que FCIQMC puede tratar núcleos de capa cerrada y abierta dentro del mismo marco y extenderse a la espectroscopía de estados excitados. Los autores señalan que el trabajo futuro involucrará interacciones de tres nucleones explícitas, bases de partícula única más grandes y aplicaciones a observables colectivos y electrodébiles.