Determination of nuclear quadrupole moments for $^{25}$Mg, $^{87}$Sr, and $^{135,137}$Ba via configuration-interaction combined with a coupled-cluster approach

Este artículo emplea un enfoque de interacción de configuración más teoría de cúmulos acoplados para calcular los gradientes de campo eléctrico y las constantes de estructura hiperfina de dipolo magnético para los estados de baja energía de Mg, Sr y Ba, permitiendo la determinación precisa de los momentos cuadripolares nucleares para $^{25}$Mg, $^{87}$Sr y $^{135,137}$Ba que revelan discrepancias significativas con los valores adoptados previamente para el estroncio y el bario.

Lo esencial

Imagina el núcleo de un átomo no como una esfera perfecta y redonda, sino como un globo ligeramente aplastado o estirado. Esta forma no es aleatoria; es una "firma" específica de cómo están dispuestos los protones y neutrones en su interior. Los científicos llaman a esta forma momento cuadrupolar nuclear. Conocer la forma exacta de este "globo" es crucial para comprender las reglas fundamentales de la física, desde cómo los átomos se mantienen unidos hasta cómo se comportan los materiales.

Sin embargo, medir esta deformación directamente es increíblemente difícil. Es como intentar adivinar la forma exacta de un globo dentro de una caja sellada y oscura simplemente escuchando el sonido que hace al golpearlo.

El Experimento: Escuchar el "Zumbido" Atómico

En este artículo, el autor, Yong-Bo Tang, actúa como un maestro ingeniero de sonido. Se centra en tres "familias" específicas de átomos: Magnesio (Mg), Estroncio (Sr) y Bario (Ba).

Cuando estos átomos se excitan (como cuando se pulsa la cuerda de una guitarra), emiten un "zumbido" o vibración muy específica llamada estructura hiperfina. Este zumbido es causado por la interacción entre los electrones que giran en el exterior y el núcleo aplastado en el interior.

• La parte medida: Los científicos ya han medido el tono de este zumbido con mucha precisión en un laboratorio.
• El eslabón perdido: Para determinar la forma del núcleo (el momento cuadrupolat) a partir de ese tono, es necesario saber exactamente cómo están dispuestos los electrones alrededor del núcleo. Esta disposición crea un "gradiente de campo eléctrico" (piensa en ello como la pendiente de una colina por la que los electrones ruedan hacia abajo).

El Problema: La Colina es Demasiado Empinada para Calcularla

Calcular la forma de esa "colina" (el gradiente de campo eléctrico) es una pesadilla para las computadoras. Los electrones no se quedan quietos; bailan unos alrededor de otros, empujándose y tirando de sí mismos de formas complejas llamadas correlación electrónica.

• Si ignoras estos bailes, tu cálculo de la colina será erróneo.
• Si intentas calcular cada uno de estos bailes, tu computadora colapsará.

Los intentos previos para calcular esto fueron como tratar de mapear una cordillera usando una foto satelital borrosa. Los resultados eran inconsistentes. Para el Estroncio y el Bario, diferentes estudios dieron respuestas distintas, con algunos resultados variando hasta un 10%.

La Solución: Un Enfoque de "Navaja Suiza" Híbrido

Para resolver esto, Tang desarrolló un nuevo método computacional que combina dos técnicas poderosas:

1. Interacción de Configuración (CI): Esto es como observar cada una de las posibles formas en que los electrones podrían disponerse, una por una. Es exhaustivo pero lento.
2. Clúster Acoplado (CC): Esto es como usar un atajo sofisticado para predecir cómo los electrones se influyen entre sí en grupos. Es rápido pero a veces pierde los detalles minúsculos.

El método de Tang, CI+CC, es lo mejor de ambos mundos. Utiliza el "atajo" para manejar las interacciones grandes y pesadas entre los electrones del núcleo, y luego utiliza el método "exhaustivo" para perfeccionar los detalles de los electrones externos. Es como usar un dron para mapear la forma general de un bosque, y luego enviar un equipo de excursionistas para medir la altura exacta de cada árbol específico.

Los Resultados: Aclarando la Confusión

Utilizando esta "navaja suiza" de alta precisión, Tang calculó los gradientes de campo eléctrico para varios estados de baja energía de Mg, Sr y Ba. Luego combinó sus cálculos con los "zumbidos" experimentales conocidos para determinar las formas nucleares.

Esto es lo que encontró:

• Magnesio (25Mg): El resultado fue una coincidencia perfecta con experimentos previos. Es como sintonizar una radio y encontrar la estación con total claridad. La forma calculada concuerda con lo que se encontró utilizando experimentos de "rayos X muónicos" (otra forma de medición de alta tecnología).
• Estroncio (87Sr): Aquí, la trama se complica. El resultado de Tang sugiere que el núcleo está aproximadamente un 10% más aplastado que el valor actualmente aceptado en los libros de texto. El valor antiguo provenía de observar un ión de Estroncio (un átomo que perdió un electrón), mientras que Tang observó el átomo neutro. La diferencia sugiere que el cálculo antiguo pudo haber pasado por alto algunos bailes electrónicos sutiles.
• Bario (135,137Ba): De manera similar al Estroncio, los resultados de Tang para el Bario son aproximadamente un 4% diferentes de los valores actualmente aceptados derivados de iones de Bario.

La Conclusión

El artículo concluye que, si bien el método funciona maravillosamente para el Magnesio, existe una discrepancia significativa para el Estroncio y el Bario en comparación con los valores de "estándar de oro" actualmente utilizados por los científicos.

Tang sugiere que la diferencia podría deberse a que los cálculos del "estándar de oro" actuales pasaron por alto un tipo específico de interacción electrónica llamada excitación triple (donde tres electrones interactúan simultáneamente). Así como un coro suena diferente si tres cantantes armonizan de una manera que nadie predijo, estas interacciones triples podrían estar desplazando el "tono" de la forma de un átomo.

En resumen: El autor construyó un mejor modelo computacional para medir la forma de los núcleos atómicos. Para el Magnesio, el modelo confirmó lo que ya sabíamos. Para el Estroncio y el Bario, el modelo sugiere que las mediciones "oficiales" actuales podrían estar ligeramente erradas, insinuando que necesitamos mirar más de cerca cómo interactúan tres electrones para obtener la verdadera forma de estos núcleos atómicos.

Resumen técnico

Resumen Técnico: Determinación de los Momentos Cuadrupolares Nucleares para $^{25}\text{Mg}$, $^{87}\text{Sr}$, y $^{135,137}\text{Ba}$

Planteamiento del Problema
El momento cuadrupolar eléctrico nuclear ($Q$) es un parámetro fundamental que describe la simetría de la estructura nuclear, con implicaciones significativas en la física atómica, nuclear y de la materia condensada. Si bien las constantes de estructura hiperfina experimentales ($B$) pueden medirse con alta precisión, la extracción de valores precisos de $Q$ requiere cálculos teóricos igualmente precisos de los gradientes de campo eléctrico ($q$). La literatura existente revela discrepancias notables en los valores recomendados para los momentos cuadrupolares nucleares de $^{87}\text{Sr}$ y $^{137}\text{Ba}$, particularmente al comparar los valores derivados de átomos neutros frente a los de iones. Por ejemplo, los valores derivados del átomo de Sr neutro son aproximadamente un 8% mayores que los recomendados a partir de datos de iones, y existen inconsistencias similares para el Ba. Además, los enfoques teóricos estándar suelen tener dificultades para contabilizar los complejos efectos de correlación electrónica necesarios para cálculos de gradiente de alta precisión en sistemas pesados de alcalinotérreos.

Metodología
Los autores emplean un enfoque híbrido relativista que combina el método de Interacción de Configuraciónación (CI) con el método de Clúster Acoplado (CC), denominado método RCI+all-order. Este marco está diseñado para contabilizar exhaustivamente las correlaciones electrónicas núcleo-núcleo, núcleo-valencia y valencia-valencia, las cuales son decisivas para el cálculo preciso del gradiente de campo eléctrico.

Componentes metodológicos clave incluyen:

• Hamiltoniano e Interacción: La interacción hiperfina se modela utilizando operadores de tensores esféricos para interacciones de dipolo magnético ($k=1$) y cuadrupolo eléctrico ($k=2$). El núcleo se trata como una esfera uniformemente magnetizada y cargada.
• Esquema de Cálculo Híbrido: El método utiliza un conjunto de base finito de funciones de tipo Gaussiano de temperatura uniforme para expandir las funciones de onda de Dirac. Los potenciales de correlación se derivan mediante cinco estrategias computacionales distintas para evaluar la incertidumbre:
  1. CI + Teoría de Perturbación de Muchos Cuerpos de segundo orden (MBPT(2)).
  2. CI + Clúster Acoplado Linealizado de Singulares y Dobles (LCCSD).
  3. CI + Clúster Acoplado Completo de Singulares y Dobles (CCSD).
  4. CI + LCCSD con parámetros de reescalado (CI+LCCSDs).
  5. CI + CCSD con parámetros de reescalado (CI+CCSDs).
• Correcciones: Los elementos de la matriz de transición incluyen la Aproximación de Fase Aleatoria (RPA), efectos de Brueckner del núcleo, radiación estructural y correcciones de normalización (denominadas colectivamente correcciones de alto orden o HO), así como interacciones de dos partículas (TP) hasta el segundo orden.
• Extracción de Q: El momento cuadrupolar nuclear se extrae utilizando la relación $Q = B / (234.9648867 \cdot q)$, donde $B$ es la constante de estructura hiperfina medida experimentalmente y $q$ es el gradiente de campo eléctrico calculado teóricamente.

Contribuciones Clave y Resultados
El estudio calcula sistemáticamente los gradientes de campo eléctrico y las constantes de estructura hiperfina de dipolo magnético ($A$) para los estados de baja energía de $^{25}\text{Mg}$, $^{87}\text{Sr}$ y $^{135,137}\text{Ba}$ neutros. Los estados específicos analizados son:

• Mg: $3s3p \ ^3P_{1,2}$
• Sr: $5s4d \ ^1D_2$ y $5s5p \ ^3P_{1,2}$
• Ba: $6s5d \ ^3D_{1,2,3}$, $6s5d \ ^1D_2$, y $6s6p \ ^1P_1$

Energía y Constantes Hiperfinas:
Las energías calculadas y las constantes de dipolo magnético ($A$) muestran un excelente acuerdo con los datos experimentales (NIST) y otros resultados teóricos. Se identifica al método CI+LCCSD como el más fiable, con los valores finales recomendados para $A$ coincidiendo con las mediciones experimentales dentro de sus incertidumbres (típicamente dentro del 2% para la mayoría de los estados). Esto valida el enfoque híbrido CI+CC para sistemas de alcalinotérreos.

Momentos Cuadrupolares Nucleares ($Q$):
Combinando los gradientes calculados con los valores experimentales de $B$, los autores determinan los siguientes momentos cuadrupolares nucleares:

• $^{25}\text{Mg}$: $0.203(2)$ b. Este resultado está en perfecto acuerdo con el experimento de rayos X mesónicos ($0.201(3)$ b) y con las determinaciones previas basadas en estructura hiperfina.
• $^{87}\text{Sr}$: $0.336(4)$ b. Este valor es aproximadamente un 10% mayor que el valor adoptado actualmente de $0.305(2)$ b (derivado del estado $4d_{5/2}$ de Sr$^+$), pero se alinea bien con los valores derivados del estado $5p_{3/2}$ de Sr$^+$ y los estados de Sr neutro ($5s5p \ ^3P_{1,2}$).
• $^{137}\text{Ba}$: $0.246(4)$ b. Es aproximadamente un 4% mayor que el valor recomendado de $0.236(3)$ b (derivado de los estados $5d_{3/2,5/2}$ de Ba$^+$). Sin embargo, concuerda excelentemente con los valores derivados del estado $6p_{3/2}$ de Ba$^+$.
• $^{135}\text{Ba}$: $0.161(2)$ b.

Significancia y Reivindicaciones
El artículo sostiene que el enfoque de interacción de configuración más clúster acoplado proporciona un marco robusto y preciso para determinar los momentos cuadrupolares nucleares, particularmente para núcleos inestables o pesados donde la medición directa es difícil.

Los autores destacan una discrepancia significativa entre los valores derivados de átomos neutros y los de iones para Sr y Ba. Específicamente, los resultados para $^{87}\text{Sr}$ y $^{137}\text{Ba}$ derivados de átomos neutros en este trabajo difieren en un ~10% y ~4%, respectivamente, de los valores actualmente adoptados en la literatura (que dependen de datos de iones). Los autores sugieren que estas discrepancias pueden provenir de la omisión de las contribuciones de excitaciones triples en los cálculos de clúster acoplado estándar, señalando que trabajos recientes sobre iones indican que las excitaciones triples contribuyen con un 1–2% a las constantes hiperfinas.

El estudio concluye que, si bien sus resultados para $^{25}\text{Mg}$ son consistentes con los datos establecidos, los hallazgos para $^{87}\text{Sr}$ y $^{137}\text{Ba}$ desafían los valores actualmente adoptados derivados de estados iónicos. Los autores afirman que sus extracciones basadas en átomos neutros son fiables, respaldadas por la alta precisión de sus cálculos de la constante de dipolo magnético, pero reconocen que es necesario un mayor trabajo teórico que incluya excitaciones triples e interacciones de tres cuerpos para resolver plenamente las diferencias entre las determinaciones basadas en átomos neutros e iones.