NuLattice: Ab initio computations of atomic nuclei on lattices

Resumen Técnico de "NuLattice: Computaciones ab initio de núcleos atómicos en redes"

Planteamiento del Problema
Las computaciones ab initio de núcleos atómicos, que utilizan Hamiltonianos derivados de teorías de campo efectivo (EFT) de la cromodinámica cuántica, dependen tradicionalmente de bases de oscilador armónico. Si bien este enfoque es eficaz para muchas aplicaciones, sufre cuellos de botella computacionales significativos cuando se trata de grandes espacios de modelos o fuerzas de tres nucleones. Específicamente, la base del oscilador armónico "desordena" las interacciones nucleares locales de corto alcance, transformándolas en operadores no locales con representaciones matriciales densas. Para las fuerzas de tres nucleones, esto resulta en la necesidad de generar y almacenar decenas de terabytes de datos, haciendo que la construcción del Hamiltoniano sea más costosa computacionalmente que la resolución posterior del problema de muchos cuerpos. En consecuencia, estos métodos suelen requerir recursos de supercomputación masivos y son inaccesibles para muchos investigadores y educadores.

Metodología
El artículo presenta NuLattice, un paquete de software en Python diseñado para realizar computaciones ab initio utilizando una red espacial como la base de una sola partícula. Este enfoque aprovecha la naturaleza de corto alcance de las fuerzas nucleares (específicamente dentro de la teoría de campo efectivo sin piones en su orden principal) para mantener estructuras de datos dispersas.

• Formulación de la Red: La base de una sola partícula se define en una red cúbica con $L^3$ sitios. El Hamiltoniano consiste en energía cinética (acoplamiento de sitios vecinos), interacciones de contacto de dos cuerpos y de tres cuerpos. Debido a que las interacciones son locales (en el sitio o de vecino más cercano), el número de elementos de matriz no nulos escala linealmente con el número de sitios de la red ($O(D)$), en marcado contraste con la escala $O(D^3)$ o $O(D^5)$ encontrada en las bases de oscilador armónico para los términos de dos y tres cuerpos.
• Algoritmos: NuLattice implementa cuatro métodos computacionales primarios:
  1. Interacción de Configuración Completa (FCI): Utilizada para núcleos ligeros ($A=2,3,4$), empleando scipy.sparse para un almacenamiento y computación de autovalores eficientes.
  2. Hartree-Fock (HF): Implementada mediante operaciones de datos dispersos para computar la matriz de densidad y la energía.
  3. Coupled-Cluster (CCSD): El paquete emplea un estado de referencia basado en la base de la red en lugar de una base de Hartree-Fock para preservar la dispersión. Utiliza una aproximación de dos cuerpos de orden normal. El resolvedor utiliza el método de Newton inexacto con la Inversión Directa del Subespacio Iterativo (DIIS) para la convergencia. La implementación mezcla contracciones de tensores dispersos-densos y densos-densos, utilizando opt_einsum para la optimización.
  4. Grupo de Renormalización de Similitud In-Medium (IMSRG): Implementado al nivel IMSRG(2). Debido a la estructura evolutiva del Hamiltoniano durante el flujo, lo que complica los enfoques dispersos, la implementación actual utiliza tensores densos para tamaños de red pequeños ($L=2,3$).
• Interacciones: El código emplea actualmente la EFT sin piones de orden principal con términos de contacto de dos y tres cuerpos. Los parámetros se calibran con las energías del estado fundamental de $^2$H y $^4$He.

Contribuciones Clave

• Lanzamiento de Software: Los autores lanzan NuLattice como un paquete de Python de código abierto (BSD-3-clause), haciendo que los cálculos de estructura nuclear ab initio sean accesibles en computadoras portátiles y de escritorio estándar.
• Implementación Dispersa: Al explotar la localidad de las interacciones de la red, el paquete evita los cuellos de botella de almacenamiento y generación asociados con las fuerzas de tres nucleones en bases de oscilador armónico.
• Utilidad Educativa: El paquete está diseñado con funciones puras y documentación extensa (vía Sphinx) para facilitar su uso en investigación y educación, incluyendo su uso en escuelas de verano recientes para enseñar patrones de convergencia y renormalización.

Resultados
El artículo presenta resultados para núcleos ligeros ($^2$H, $^3$He, $^4$He, $^8$Be, $^{12}$C y $^{16}$O) en redes con un espaciado $a \approx 1.7\text{--}2.5$ fm.

• Núcleos Ligeros ($A=2,3,4$): Para $^2$H, CCSD produce la solución exacta, mientras que Hartree-Fock falla en converger. Para $^3$He y $^4$He, los resultados muestran que Hartree-Fock y Coupled-Cluster con Singulares (CCS) capturan la mayor parte de la energía del estado fundamental. CCSD proporciona solo mejoras marginales sobre CCS, lo que sugiere que para estos núcleos compactos con interacciones de rango cero, las excitaciones de una sola partícula dominan.
• Convergencia: La energía converge rápidamente con el tamaño de la red $L$. Para $^{16}$O, la energía cambia por menos de 0.2 MeV entre $L=4$ y $L=5$.
• IMSRG vs. CCSD: Los resultados de IMSRG(2) y CCSD son consistentes entre sí y cercanos a Hartree-Fock, lo que indica que las correlaciones de orden superior (más allá de singulares y dobles) contribuyen poco a la energía de enlace en este marco específico de rango cero.
• Limitaciones: El estudio confirma que la EFT sin piones de orden principal con interacciones de rango cero no logra ligar partículas $\alpha$ en núcleos más grandes como $^8$Be, $^{12}$C y $^{16}$O contra la desintegración $\alpha$. Esto se atribuye a la naturaleza de rango cero del potencial, que no puede ligar sistemas más allá de la partícula $\alpha$. Además, para sistemas débilmente ligados como $^3$He, los efectos de tamaño finito y las correlaciones de tres partículas-tres huecos son significativos, lo que conduce a mayores discrepancias entre los métodos aproximados (CCSD/IMSRG) y los referentes de FCI en redes pequeñas.

Significancia y Reivindicaciones
Los autores afirman que NuLattice representa un paso significativo hacia la democratización de la física nuclear ab initio. Al reducir el costo computacional mediante técnicas de red dispersas, el paquete permite que computaciones no triviales se realicen en computadoras portátiles, haciendo que estos métodos sean accesibles para un grupo más amplio de investigadores y educadores. El artículo enfatiza que, si bien la implementación actual se limita a la EFT sin piones de orden principal e interacciones de rango cero, el marco demuestra con éxito que los complejos resolvedores de muchos cuerpos (CCSD, IMSRG, FCI) pueden implementarse eficientemente en Python. Los autores invitan a la comunidad a contribuir al paquete, señalando que el trabajo futuro se centrará en interacciones de rango finito, EFT quiral, resolvedores de mayor precisión y el cálculo de estados excitados y observables más allá de las energías del estado fundamental.