Shell model description of the $N=82$ isotonic chain with a new effective interaction

Resumen Técnico: Descripción del Modelo de Capas de la Cadena Isotónica de N = 82 con una Nueva Interacción Efectiva

Planteamiento del Problema
La estructura nuclear de la cadena isotónica de $N=82$, caracterizada por protones de valencia en la capa $g_{7/2}d_{5/2}s_{1/2}d_{3/2}h_{11/2}$ por encima del núcleo $^{132}\text{Sn}$, ha sido objeto de extensos estudios. Si bien interacciones efectivas previas (por ejemplo, JJ56PNA, KH5082, CW5082) han proporcionado información valiosa, estas exhiben discrepancias sistemáticas al compararse con los datos experimentales. Notablemente, estas interacciones tienden a sobreestimar las energías de excitación de los estados $3^-$ y no logran predecir correctamente los espines y paridades de los estados fundamentales para los isotonos de $A$ impar más pesados ($Z > 64$). Además, se requieren predicciones teóricas fiables para núcleos ricos en protones cerca y más allá de la línea de gotera de protones (por ejemplo, $^{155}\text{Ta}$, $^{156}\text{W}$), donde los datos experimentales son escasos o inexistentes.

Metodología
Los autores emplean el modelo de capas de interacción de configuración completa utilizando el código BIGSTICK dentro de un espacio de modelos que comprende todos los orbitales de protones entre $Z=50$ y $Z=82$. La innovación central radica en el desarrollo de una nueva interacción efectiva de alta calidad optimizada mediante el Análisis de Componentes Principales (PCA).

1. Espacio de Parámetros: El Hamiltoniano incluye 5 energías de partícula única (SPEs) y 160 elementos de matriz de dos cuerpos (TBMEs) protón-protón, sumando un total de 165 parámetros de entrada.
2. Estrategia de Optimización: En lugar de ajustar todos los parámetros directamente, los autores utilizan PCA para truncar el espacio de parámetros. Realizan un ajuste de mínimos cuadrados a 204 puntos de datos experimentales, que incluyen:
   • 21 energías de enlace.
   • 101 energías de excitación en núcleos even-even (centrándose en estados yrast y estados no-yrast específicos como $2^+_2, 3^-$, etc.).
   • 82 energías de excitación en núcleos de $A$ impar.
3. Implementación de PCA: La matriz de error del ajuste se diagonaliza. El espacio de parámetros se trunca a los 30 componentes principales (combinaciones lineales de los parámetros originales) asociados con los autovalores más pequeños (incertidumbres más pequeñas). Este enfoque mitiga el sobreajuste (overfitting) mientras retiene los parámetros más sensibles.
4. Restricciones: Se aplican factores de escala de masa estándar a las TBMEs. Los efectos Coulomb se incorporan implícitamente a través del procedimiento de ajuste. Se utilizan cargas efectivas ($e_\pi = 1.5e$) y factores de amortiguación ($g_s = 5.586 \times 0.7$) para las propiedades electromagnéticas.

Contribuciones Clave y Resultados

• Optimización de la Interacción: La nueva interacción logra una desviación cuadrática media (RMSD) de 102 keV a través de los 204 puntos de datos. Esto representa una mejora significativa respecto a las interacciones previas, que a menudo presentaban desviaciones de escalas de MeV para las energías de enlace.
• Energías de Enlace y de Separación: El modelo reproduce las energías de enlace para $^{134}\text{Te}$ hasta $^{151}\text{Tm}$ con un RMSD de 59 keV. Predice con éxito los excesos de masa atómica y las energías de separación para núcleos más allá del alcance experimental actual, incluyendo $^{156}\text{W}$ a través de $^{159}\text{Ir}$. Los cálculos sugieren que $^{153}\text{Lu}$ está más allá de la línea de gotera de protones, mientras que $^{154}\text{Hf}$ permanece ligado. Para $^{155}\text{Ta}$, el modelo predice la emisión de un solo protón, y para $^{156}\text{W}$, predice una energía de separación de un protón positiva pero una energía de separación de dos protones negativa, lo cual es consistente con la supresión observada de la desintegración de dos protones.
• Espectros de Baja Energía:
  • Núcleos Even-Even: El modelo reproduce la tendencia sistemática de las energías de excitación $2^+_1$, identificando correctamente el pico en $^{146}\text{Gd}$ ($Z=64$) y el máximo local en $^{140}\text{Ce}$, confirmando el cierre de subcapa de $Z=64$. Resuelve las discrepancias previas al predecir con precisión las energías de excitación $3^-$, particularmente el mínimo observado en $^{146}\text{Gd}$.
  • Núcleos de A Impar: La interacción reproduce correctamente los espines y paridades de los estados fundamentales en toda la cadena, incluyendo los desafiantes estados fundamentales $1/2^+$ de $^{147}\text{Tb}$ y $3/2^+$ de $^{155}\text{Ta}$ y $^{157}\text{Re}$, que anteriormente eran difíciles de describir.
• Propiedades Electromagnéticas:
  • Valores B(E2): El modelo concuerda generalmente bien con los valores experimentales de $B(E2)$. Sin embargo, se observan desviaciones para transiciones específicas en $^{138}\text{Ba}$ y $^{139}\text{La}$, lo que sugiere una posible inversión en el orden de los estados $9/2^+$ calculados en $^{139}\text{La}$.
  • Momentos Magnéticos: Los momentos dipolares magnéticos ($\mu$) calculados concuerdan bien con el experimento (discrepancias $< 0.5 \mu_N$).
  • Fortalezas B(E3): El modelo subestima las fuerzas $B(E3; 3^- \to 0^+)$ por un factor de aproximadamente dos. Los autores atribuyen esto a las contribuciones faltantes de las excitaciones de neutrones a través de la brecha de $N=82$, señalando que aumentar la carga efectiva a 2.0 pone los resultados en concordancia.
• Análisis Estructural:
  • Brechas de Subcapa: Las energías de partícula única de campo medio (MSPEs) revelan una brecha de subcapa de 1.788 MeV en $Z=64$ (entre $1d_{5/2}$ y $0h_{11/2}$) y de 1.062 MeV en $Z=58$ (entre $0g_{7/2}$ y $1d_{5/2}$).
  • Senioridad: El estudio analiza la estructura de senioridad de los estados yrast. Los estados $6^+_1$ en isotonos más ligeros están dominados por configuraciones $(0g_{7/2})^2$ o $0g_{7/2}1d_{5/2}$, mientras que en isotonos más pesados ($Z > 64$) corresponden a $(0h_{11/2})^2$.
  • Isomerismo: El modelo predice varios estados excitados de muy baja energía en núcleos de $A$ impar (por ejemplo, un estado $1/2^+$ en $^{155}\text{Ta}$ a 4 keV), lo que sugiere un potencial isomerismo.

Significado y Reivindicaciones
El artículo sostiene que la interacción recién optimizada proporciona una descripción sistemática y fiable de la cadena isotónica de $N=82$, resolviendo las discrepancias de larga data en las energías de excitación y los espines de los estados fundamentales que afectaron a los modelos anteriores. Al utilizar PCA para optimizar la interacción contra un gran conjunto de datos, los autores demuestran que una descripción de una sola capa principal con una interacción optimizada es suficiente para capturar características estructurales complejas, como el cierre de la subcapa de $Z=64$ y la evolución de los patrones de senioridad, sin necesidad de invocar excitaciones a gran escala a través de las capas principales.

Los autores enfatizan que esta interacción sirve como una base sólida para estudios futuros de núcleos cercanos a $N=82$. Citan explícitamente su aplicación (y uso continuo) en:

1. Estudios de la evolución de la forma nuclear y transiciones de fase cuántica en núcleos de tierras raras.
2. Investigaciones sistemáticas de los isotonos even-even de $N=80$.
3. La primera medición espectroscópica de $^{150}\text{Yb}$.

El trabajo concluye proporcionando la interacción optimizada como material suplementario para facilitar la investigación posterior en esta región de masa.