Nucleon-pair truncation of the shell model for medium-heavy nuclei

Resumen Técnico: Truncamiento de pares de nucleones del modelo de capas para núcleos de masa media-pesada

Planteamiento del Problema
El modelo de capas de interacción de configuración (SM) es un marco versátil para la estructura nuclear, capaz de generar estados de baja energía y correlaciones complejas de múltiples partículas. Sin embargo, para núcleos rotacionales de masa media y pesada, las dimensiones de la interacción de configuración completa (FCI) suelen exceder los límites computacionales actuales (típicamente $\sim 2-3 \times 10^{10}$). Aunque los esquemas de truncamiento existentes, como el esquema de seniority generalizado o la aproximación de pares de nucleones (NPA), reducen el tamaño de la base, enfrentan desafíos: los esquemas de seniority funcionan mejor para núcleos esféricos, y los cálculos de NPA que involucran pares de alto espín (pares G e I) siguen siendo computacionalmente onerosos para núcleos bien deformados. Por el contrario, los enfoques de campo medio como Hartree-Fock-Bogoliubov (HFB) o Bardeen-Cooper-Schrieffer (BCS) manejan la deformación y el apareamiento de manera eficiente, pero rompen la simetría rotacional y la conservación del número de partículas. La restauración de estas simetrías mediante proyección es a menudo computacionalmente prohibitiva, particularmente al realizar una "variación después de la proyección".

Metodología
Los autores proponen un esquema de truncamiento computacionalmente eficiente denominado BCS con número proyectado y conservado (PNBCS). La metodología procede en tres etapas principales:

1. Generación de la Base de Hartree-Fock (HF): Los autores generan estados de partícula única utilizando un cálculo HF dentro de una base de modelo de capas que impone la degeneración de Kramers. Esto asegura parejas de inversión temporal degeneradas sin imponer restricciones adicionales sobre la forma o la paridad.
2. BCS con Conservación del Número de Partículas (NBCS): Construyendo sobre los estados HF, los autores construyen un ansatz de condensado de pares variacional para núcleos de capa abierta (que contienen tanto protones como neutrones de valencia). A diferencia del BCS estándar, que viola el número de partículas, este enfoque utiliza un formalismo de conservación del número de partículas (NBCS) donde el estado fundamental es un condensado de $N$ pares. Los coeficientes de la estructura de los pares se optimizan minimizando la energía mediante un método de gradiente conjugado y fórmulas iterativas derivadas del principio variacional. Este paso incorpora explícitamente las correlaciones de apareamiento.
3. Proyección de Momento Angular (LAP): Dado que los estados HF y NBCS rompen la simetría rotacional, los autores aplican una proyección de álgebra lineal (LAP) para restaurar un buen momento angular. A diferencia de la integración numérica tradicional sobre ángulos de Euler, LAP resuelve un conjunto de ecuaciones lineales para extraer los elementos de la matriz del Hamiltoniano y de la norma. Se reporta que este enfoque es más de 10 veces más rápido que los métodos de cuadratura. El espectro final se obtiene diagonalizando el Hamiltoniano en el espacio definido por los estados proyectados.

Resultados Clave
El esquema PNBCS se aplicó a núcleos a través de las capas $pf$, $pf_{5/2}g_{9/2}$ y $sdg_{7/2}h_{11/2}$, cubriendo regiones transicionales y deformadas:

• Evolución de la Forma: En un modelo esquemático que varía la relación entre el apareamiento y las interacciones cuadrupolo-cuadrupolo, el PNBCS reprodujo con éxito la evolución desde los límites esféricos a los rotacionales. Proporcionó descripciones superiores de las energías de excitación y de las tasas de transición $B(E2)$ en comparación con el PHF (Hartree-Fock proyectado de momento angular), particularmente en regiones transicionales ($2.2 < R_{4/2} < 2.8$).
• Núcleos de Masa Media (capas $pf$y$pf_{5/2}g_{9/2}$): Para núcleos tales como $^{44,46,48}$Ti, $^{48,50}$Cr, $^{52}$Fe, $^{60,62,64}$Zn, $^{66,68}$Ge y $^{68}$Se, los resultados de PNBCS mostraron una buena concordancia con los cálculos de SM completo y con los datos experimentales.
  • Se encontró que la inclusión de las correlaciones de apareamiento reduce el momento de inercia (mejorando el espaciamiento de los niveles), mientras que la deformación (codificada en HF) dicta primordialmente las fuerzas de transición cuadrupolo eléctrico.
  • Coexistencia de Formas: En $^{66}$Ge y $^{68}$Ge, donde se observa coexistencia de formas (mínimos oblatos/triaxiales o triaxiales), los autores demostraron que tanto las correlaciones de apareamiento como la mezcla de configuraciones entre diferentes estados intrínsecos de HF son cruciales para reproducir las energías de las bandas laterales y los valores de $B(E2)$. En contraste, para $^{68}$Se, se encontró que la mezcla de configuraciones era menos crítica.
• Predicciones de Masa Pesada: Los autores aplicaron PNBCS a núcleos donde los cálculos de SM completo son actualmente infactibles debido a la dimensionalidad ($2 \times 10^{10}$ a $5 \times 10^{13}$). Proporcionaron predicciones para los estados de baja energía de $^{108,110}$Xe, $^{112,114}$Ba y $^{116,118,120}$Ce. Los resultados para $^{108,110}$Xe concordaron bien con los datos disponibles y los benchmarks de SM.

Significancia y Reivindicaciones
El artículo afirma que el PNBCS ofrece un esquema de truncamiento práctico y eficiente para el modelo de capas, dirigido específicamente a núcleos transicionales y deformados en las regiones de masa media y pesada. La significancia principal radica en su capacidad para:

1. Equilibrar Precisión y Eficiencia: Captura las correlaciones de apareamiento y la mezcla de configuraciones esenciales, manteniéndose computacionalmente rápido, evitando las pesadas integraciones multidimensionales requeridas por los métodos HFB proyectados o los métodos generales de condensado de pares.
2. Extender el Alcance: Permite el estudio de núcleos (por ejemplo, $^{112-120}$Ba/Ce) que están actualmente fuera del alcance de los cálculos de SM de interacción de configuración a gran escala.
3. Reproducir la Colectividad: El estudio demuestra que tanto las correlaciones de apareamiento como la mezcla de diferentes estados intrínsecos son clave para reproducir características colectivas y la coexistencia de formas.

Los autores señalan que, si bien la implementación actual utiliza estados de partícula única de HF como una base fija para la variación de NBCS, trabajos futuros podrían explorar la variación después de la proyección o la inclusión de configuraciones de par roto para mejorar aún más los resultados para fenómenos como el backbending.