Ab initio calculations of nuclear charge radii across and beyond ${}^{132}$Sn: Putting chiral EFT nuclear interactions to the test

Este estudio utiliza cálculos *ab initio* de tipo Bogoliubov coupled cluster para demostrar que, aunque las interacciones nucleares de la teoría efectiva de campo quiral (χEFT) reproducen parcialmente los radios de carga en la cadena de isótopos de estaño, ninguna de las interacciones probadas logra capturar simultáneamente todas las características clave observadas, lo que subraya la necesidad de nuevas mediciones experimentales y cálculos de mayor precisión para refinar estos modelos teóricos.

Lo esencial

Imagina que el núcleo de un átomo es como una ciudad en constante expansión. En el centro viven los protones (cargados positivamente) y los neutrones (sin carga), todos apretujados y bailando juntos bajo la influencia de una fuerza invisible llamada "fuerza fuerte".

El tamaño de esta ciudad se mide por su "radio de carga". Los científicos quieren predecir exactamente cuán grande es esta ciudad para diferentes versiones de la misma ciudad (isótopos), especialmente en la familia del Estao (Sn), que es como una larga fila de hermanos desde el más pequeño hasta el más grande.

Aquí está lo que hicieron los autores de este estudio, explicado de forma sencilla:

1. El Reto: Predecir el tamaño de la ciudad

Los científicos tienen una "receta" teórica muy avanzada llamada Teoría de Campo Efectivo Quiral (χEFT). Es como un manual de instrucciones muy detallado que dice cómo deben comportarse los protones y neutrones entre sí.

Sin embargo, hasta ahora, estas recetas tenían un problema:

• A veces decían que la ciudad era demasiado pequeña.
• A veces no podían explicar por qué la ciudad crecía de golpe en ciertos momentos (como cuando se añaden ciertos vecinos).

2. La Prueba: La "Ciudad de Estao"

Los investigadores usaron superordenadores para simular la familia del Estao, desde el isótopo 96 hasta el 150. Se centraron en dos momentos clave:

• El "bache" en 100Sn: Un punto donde la ciudad cambia de tamaño drásticamente.
• El "codo" en 132Sn: Un punto famoso donde la ciudad da un salto de tamaño inesperado (como si de repente se pusiera un anillo de grasa extra).

Usaron tres versiones diferentes de su "receta" (tres Hamiltonianos) para ver cuál funcionaba mejor:

1. La receta vieja (1.8/2.0 EM): Predecía ciudades demasiado pequeñas.
2. La receta ajustada (ΔNNLOGO): Mejoraba un poco, pero seguía fallando en los saltos de tamaño.
3. La receta nueva y afinada (1.8/2.0 EM7.5): ¡Esta fue la ganadora! Logró predecir el tamaño correcto en el punto de 132Sn, ¡pero con un truco!

3. El Problema: La solución "correcta" por la razón "incorrecta"

Aquí viene la parte divertida y preocupante. La receta que funcionó mejor para el punto 132Sn (la 1.8/2.0 EM7.5) lo hizo porque imaginó mal la estructura de la ciudad más allá de ese punto.

• La analogía de la escalera: Imagina que la ciudad crece subiendo escaleras. La receta correcta debería decir que, al llegar al piso 82 (N=82), los vecinos suben a un piso muy alto y ancho (una órbita llamada 1h9/2), lo que hace que la ciudad se expanda mucho (el "codo" en 132Sn).
• El error: Esa misma receta predijo que, al llegar al piso 92 (N=92), los vecinos bajarían a un piso estrecho y pequeño. Esto crearía un "codo invertido" (la ciudad se encogería un poco), algo que no debería pasar.

Básicamente, la receta adivinó el tamaño correcto en un punto por casualidad, pero su mapa del futuro (más allá de 132Sn) estaba equivocado.

4. ¿Qué aprendemos de esto?

El estudio nos dice que:

• Nuestras recetas aún no son perfectas: Ninguna de las tres versiones probadas pudo explicar todo el comportamiento de la ciudad del Estao a la vez.
• Necesitamos más datos: Para arreglar la receta, necesitamos medir el tamaño de las ciudades del Estao que son más raras y pesadas (más allá de 134Sn) y las que son más ligeras (hacia 100Sn). Es como si necesitáramos ver más fotos de la ciudad para entender sus reglas de crecimiento.
• La computación debe mejorar: Los cálculos actuales son muy buenos, pero para ser perfectos, necesitan incluir "correcciones de triplets" (como si contáramos no solo parejas de vecinos, sino también grupos de tres que interactúan) y correcciones relativistas.

En resumen

Los científicos usaron superordenadores para probar si sus teorías sobre cómo crecen los núcleos atómicos eran correctas. Descubrieron que, aunque algunas teorías adivinan bien el tamaño en un punto específico, fallan al predecir cómo se comportarán en el futuro.

La moraleja: Tenemos un mapa muy bueno, pero necesitamos más exploradores (experimentos) y mejores herramientas de cálculo para dibujar el mapa completo de la "ciudad nuclear" y entender las reglas exactas que gobiernan el universo.

Resumen técnico

Título: Cálculos ab initio de radios de carga nuclear a través y más allá de $^{132}$Sn: Poniendo a prueba las interacciones nucleares de la Teoría del Campo Efectivo Quiral (χEFT)

1. El Problema

La descripción precisa de los radios de carga nuclear representa uno de los desafíos más difíciles para las teorías de ab initio (de primeros principios) basadas en la Teoría del Campo Efectivo Quiral (χEFT).

• Desafíos previos: En sistemas más ligeros, las interacciones nucleares estándar tienden a subestimar sistemáticamente los radios de carga absolutos y no logran reproducir comportamientos complejos, como el comportamiento parabólico invertido entre $^{40}$Ca y $^{48}$Ca o el aumento marcado hacia $^{52}$Ca.
• El caso del Estaño (Sn): La cadena isotópica del Estaño ($Z=50$) es un laboratorio ideal para probar estas interacciones debido a su extensión (desde $^{96}$Sn hasta $^{150}$Sn) y a la presencia de números mágicos ($N=50$ y $N=82$). Recientemente, mediciones experimentales han revelado un "codo" (kink) pronunciado en el radio de carga al cruzar el isótopo doblemente mágico $^{132}$Sn (específicamente en $^{134}$Sn), así como datos en isótopos deficientes en neutrones.
• La incógnita: Ninguna de las interacciones χEFT ajustadas ("fine-tuned") disponibles hasta la fecha ha logrado reproducir simultáneamente todas las características clave (energías de enlace, radios absolutos y la evolución de los desplazamientos isotópicos) en esta región de masa pesada.

2. Metodología

Los autores realizaron cálculos ab initio utilizando el marco de la Teoría de Clúster Acoplado de Bogoliubov (BCC):

• Nivel de aproximación: Se trabajó a nivel de singles y doubles (BCCSD), truncando el operador de clúster $T$ y el operador de desexcitación $\Lambda$ en sus componentes de un y dos cuasipartículas.
• Hamiltonianos: Se emplearon tres interacciones nucleares basadas en χEFT de última generación, todas ajustadas a diferentes propiedades:
  1. 1.8/2.0 (EM): Ajustada para energías de enlace, pero subestima los radios.
  2. $\Delta$NNLOGO: Incluye grados de libertad de delta ($\Delta$) y se ajusta a propiedades de saturación de materia nuclear para mejorar los radios absolutos.
  3. 1.8/2.0 (EM7.5): Recientemente propuesta, con una constante de baja energía de cuerpo tres ($c_D$) reajustada para reproducir la energía y el radio de $^{16}$O, mejorando drásticamente los radios en núcleos cerrados.
• Espacio de base: Se utilizó una base de oscilador armónico esférico con 13 capas principales ($e_{max}=12$) y frecuencias óptimas ($\hbar\omega = 12$ MeV).
• Operadores: Se calcularon los radios cuadráticos medios de carga ($R_{ch}^2$) incluyendo correcciones de centro de masa, tamaño del protón, tamaño del neutrón, corrección spin-órbita y el término relativista de Darwin-Foldy.
• Validación: Los resultados se compararon con datos experimentales de espectroscopía láser (ISOLDE), cálculos VS-IMSRG(2) y, en algunos casos, con cálculos de teoría de campo efectivo en red (NLEFT).

3. Contribuciones Clave

• Evaluación exhaustiva: Se realizó el primer estudio sistemático de radios de carga en toda la cadena de isótopos de Estaño ($96 \le A \le 150$) utilizando BCCSD con tres Hamiltonianos χEFT modernos.
• Análisis de la sensibilidad: Se demostró que la región más allá de $N=82$ (más allá de $^{132}$Sn) actúa como un "campo de juego" único donde la sensibilidad de los resultados a los detalles del Hamiltoniano se amplifica, permitiendo discriminar entre diferentes interacciones.
• Descomposición de desplazamientos: Se analizó la estructura de los desplazamientos isotópicos, separando sus componentes lineales y parabólicas, y se identificó cómo las subcapas y la estructura de capas de partículas individuales influyen en estos patrones.
• Identificación de limitaciones: Se reveló que, aunque algunas interacciones pueden reproducir un "codo" específico, lo hacen a menudo debido a una estructura de capas incorrecta, lo que sugiere que la concordancia con el experimento puede ser fortuita.

4. Resultados Principales

• Radios Absolutos:
  • El Hamiltoniano 1.8/2.0 (EM) subestima significativamente los radios (error RMS ~5.6%).
  • Los Hamiltonianos $\Delta$NNLOGO y 1.8/2.0 (EM7.5) mejoran drásticamente la predicción de radios absolutos (reduciendo el error a ~1.2% y ~0.7% respectivamente), corrigiendo la deficiencia del modelo de campo medio.
• Desplazamientos Isotópicos y el "Codo" en $^{132}$Sn:
  • Ninguna interacción captura perfectamente todas las características.
  • El Hamiltoniano 1.8/2.0 (EM7.5) reproduce con gran precisión el aumento experimental del radio entre $^{132}$Sn y $^{134}$Sn. Sin embargo, esto se debe a que predice incorrectamente que los neutrones ocupan la orbital $1h_{9/2}$ justo después de $N=82$. Esta ocupación incorrecta genera un aumento artificialmente grande del radio.
  • Los Hamiltonianos 1.8/2.0 (EM) y $\Delta$NNLOGO predicen correctamente la ocupación de la orbital $2f_{7/2}$ (coherente con el espín experimental de $^{133}$Sn), pero subestiman la magnitud del aumento del radio en $^{134}$Sn.
• Predicción en $^{142}$Sn:
  • El Hamiltoniano 1.8/2.0 (EM7.5) predice un "codo invertido" en $^{142}$Sn, correlacionado con un aumento no físico en la energía de separación de dos neutrones ($S_{2n}$). Esto pone en duda la validez de su éxito previo en $^{134}$Sn.
  • Los otros dos Hamiltonianos muestran una variación suave en esta región, sin el codo invertido.
• Componente Parabólica: Entre $N=50$ y $N=82$, la componente cuasi-parabólica de los desplazamientos isotópicos se reproduce bien por los Hamiltonianos 1.8/2.0 (EM) y $\Delta$NNLOGO, pero se subestima significativamente por el 1.8/2.0 (EM7.5).

5. Significado y Conclusiones

• Desafío para la Teoría: El estudio concluye que lograr una descripción consistente de las energías de enlace y los radios de carga más allá de $^{132}$Sn sigue siendo un desafío mayor para los cálculos ab initio y los Hamiltonianos χEFT.
• Falsabilidad de las Interacciones: La región más allá de $N=82$ es crítica para validar las interacciones nucleares. El hecho de que el Hamiltoniano 1.8/2.0 (EM7.5) reproduzca el dato de $^{134}$Sn por la "razón incorrecta" (estructura de capas errónea) subraya la necesidad de medir isótopos más allá de $^{134}$Sn y hacia $^{100}$Sn para discriminar entre modelos.
• Futuro de los Cálculos: Para alcanzar precisiones subporcentuales y resolver las discrepancias actuales, se requiere:
  1. Incluir correcciones de triples ($T_3$) en la función de onda (actualmente omitidas en BCCSD).
  2. Incorporar correcciones de densidad de carga de dos cuerpos en el operador.
  3. Realizar nuevas mediciones experimentales de desplazamientos isotópicos en isótopos extremos.
  4. Considerar correlaciones colectivas verdaderas que son difíciles de capturar con métodos de expansión de referencia única.

En resumen, el trabajo demuestra que la cadena de isótopos de Estaño actúa como un test riguroso que expone las limitaciones actuales de las interacciones χEFT, señalando la necesidad de refinamientos teóricos (triples, operadores de dos cuerpos) y nuevas mediciones experimentales para avanzar en la comprensión de la fuerza nuclear fuerte en sistemas pesados.