Investigating Role of Electron Correlation Effects via Triple Excitations for Precise Evaluation of Energies and Hyperfine Structure Constants in $^{23}$Na

Resumen Técnico: El papel de las excitaciones triples en la estructura hiperfina del $^{23}$Na

Planteamiento del problema
La determinación precisa de las constantes de la estructura hiperfina (HFSC) es crítica para comprender la interacción entre la correlación electrónica y los efectos relativistas en la región nuclear. Aunque el sodio ($^{23}$Na) es un átomo alcalino relativamente ligero con una estructura simple de un solo electrón de valencia, los cálculos previos de todos los órdenes de la teoría de muchos cuerpos relativista para sus constantes de estructura hiperfina de dipolo magnético ($A_{hf}$) han mostrado discrepancias notables con los datos experimentales. Específicamente, los cálculos anteriores de Relativistic Coupled-Cluster (RCC) utilizando el método de Singles and Doubles (RCCSD) arrojaron valores que eran ligeramente mayores o sustancialmente menores que los resultados experimentales, dependiendo de la inclusión de términos no lineales. Además, se encontró previamente que la inclusión perturbativa de excitaciones triples en el $^{23}$Na casi cancelaba las contribuciones no lineales de RCCSD, pero el papel individual de los efectos de correlación de orden superior y su interacción con las correcciones relativistas (Breit, QED y Bohr–Weisskopf) permanecía incompletamente comprendido. Este estudio aborda la necesidad de un entendimiento cuantitativo de estas contribuciones para establecer la confiabilidad de los métodos de muchos cuerpos para el $^{23}$Na y para proporcionar un referente para sistemas alcalinos más pesados.

Metodología
Los autores emplearon la teoría de relativistic coupled-cluster (RCC) con la inclusión explícita de excitaciones triples para calcular los potenciales de ionización (IP) y las constantes de estructura hiperfina ($A_{hf}$ y la constante de cuadrupolo eléctrico $B_{hf}$) para 11 estados de baja energía del $^{23}$Na. El marco computacional incluyó:

• Hamiltoniano y correcciones: Se utilizó el Hamiltoniano de Dirac-Coulomb (DC) como base. Para una alta precisión, se incorporaron correcciones de la interacción de Breit, efectos de la Electrodinámica Cuántica (QED) (específicamente el potencial de Uehling y Wichmann-Kroll para la polarización del vacío, y el enfoque de potencial radiativo de Flambaum-Ginges para la autoenergía) y el efecto Bohr–Weisskopf (distribución de magnetización nuclear).
• Enfoques de muchos cuerpos: Para evaluar sistemáticamente la importancia de diferentes contribuciones físicas, los autores compararon los resultados a través de múltiples niveles de aproximación:
  • Dirac-Hartree-Fock (DHF).
  • Teoría de Muchos Cuerpos Relativista (RMBPT) hasta el tercer orden.
  • Aproximación de Fase Aleatoria (RPA) y sus extensiones con orbitales Brückner (BO), radiación estructural (SR) y normalización (Nm).
  • Métodos de Relativistic Coupled-Cluster: RCCSD (Singles y Doubles), RCCSDTv (Singles, Doubles y triples de valencia) y el RCCSDT completo (Singles, Doubles y triples completos).
• Conjuntos de bases: Los cálculos utilizaron Orbitales de Tipo Gaussiano (GTO) con hasta 40 funciones por momento angular orbital ($l$) hasta simetría $i$ ($l=6$), asegurando una descripción equilibrada de las regiones cercanas al núcleo y asintóticas.

Contribuciones clave y resultados

1. Potenciales de Ionización (IPs):
   Los IPs calculados mostraron que los efectos de correlación electrónica aumentan sistemáticamente los valores desde el nivel DHF hasta el nivel RCCSDTv. Sin embargo, la inclusión de excitaciones triples de núcleo completas en el método RCCSDT condujo a una ligera reducción, llevando los valores teóricos a un excelente acuerdo con los datos experimentales de NIST (discrepancias < 0.04%). Las correcciones relativistas y de QED tuvieron un impacto marginal en los IPs.

2. Constantes de dipolo magnético ($A_{hf}$):

   • Sensibilidad a la correlación: A diferencia de los IPs, las constantes $A_{hf}$ exhiben una sensibilidad extrema a la correlación electrónica debido a la localización del operador en la región nuclear. Los valores de DHF subestimaron el $A_{hf}$ del estado fundamental en casi un 30%.
   • Papel de las excitaciones triples: El método RCCSD subestimó consistentemente los valores experimentales, mientras que el RCCSDTv (que incluye solo triples de valencia) los sobreestimó. El método RCCSDT completo, que incluye excitaciones triples de núcleo, llevó los resultados a un estrecho acuerdo con el experimento. El estudio demuestra que las excitaciones triples de núcleo cancelan parcialmente las contribuciones de los triples de valencia.
   • Efectos directos vs. indirectos: Un análisis detallado de los términos RCC reveló que la influencia primaria de las excitaciones triples es indirecta. Aunque los términos explícitos que involucran operadores de excitación triple ($T_3, S_{3v}$) son pequeños, la inclusión de triples en las ecuaciones RCC modifica de manera autoconsistente las amplitudes de excitación simple y doble ($T_1, S_{1v}, S_{2v}$), las cuales dominan el resultado final. Esto subraya la necesidad de un tratamiento autoconsistente de los triples en lugar de uno perturbativo.
   • Correcciones nucleares y relativistas: Aunque las contribuciones individuales de Breit, QED y BW son pequeñas, su efecto combinado es comparable a las correcciones de correlación electrónica de orden superior. La inclusión de estas correcciones fue esencial para reconciliar los resultados de RCCSDT con las observaciones experimentales.

3. Constantes de cuadrupolo eléctrico ($B_{hf}$):
   Los efectos de correlación electrónica fueron menos pronunciados para $B_{hf}$ que para $A_{hf}$. Se encontró que los efectos de polarización del núcleo (CP) son el mecanismo de correlación dominante. De manera similar a $A_{hf}$, la inclusión de excitaciones triples (RCCSDT) redujo ligeramente los valores en comparación con RCCSDtv, indicando una cancelación parcial entre las contribuciones de los triples de núcleo y de valencia.

Significado y afirmaciones
El artículo afirma que lograr una precisión sub-porcentual en los cálculos de la estructura hiperfina para átomos alcalinos requiere un tratamiento riguroso de la correlación electrónica, específicamente la inclusión autoconsistente de excitaciones triples. El estudio demuestra que:

• Las correcciones de orden inferior, tanto relativistas como nucleares (BW), desempeñan papeles comparables a los efectos de correlación electrónica y son indispensables para la precisión.
• La interacción entre las excitaciones triples de núcleo y de valencia es crítica; omitir los triples de núcleo o tratar los triples de forma perturbativa conduce a errores sistemáticos.
• El alto grado de acuerdo entre las predicciones teóricas basadas en RCCSDT y los datos experimentales valida la confiabilidad de las funciones de onda atómicas generadas.
• Este trabajo sirve como guía para comprender el papel de las excitaciones triples en sistemas alcalinos más pesados, donde se espera que los efectos de correlación sean aún más pronunciados.

Los autores concluyen que su comparación sistemática de métodos de orden inferior contra el enfoque de todos los órdenes de RCCSDT proporciona un marco robusto para aislar y cuantificar diferentes clases de efectos de correlación electrónica, esenciales para desarrollar teorías de muchos cuerpos confiables aplicables a una amplia gama de sistemas cuánticos.