Electromagnetic sum rules for 22O from coupled-cluster theory

Resumen Técnico: Reglas de Suma Electromagnética para $^{22}$O mediante la Teoría de Clúster Acoplado

Planteamiento del Problema
El artículo aborda el desafío de describir las funciones de respuesta electrodébiles en núcleos de masa media, centráose específicamente en el isótopo rico en neutrones $^{22}$O. Si bien los métodos ab initio han descrito con éxito las propiedades del estado fundamental, el cálculo de las funciones de respuesta que involucran el espectro continuo sigue siendo difícil debido a la complejidad de manejar estados no ligados. Los autores pretenden determinar la polarizabilidad dieléctrica ($\alpha_D$) para $^{22}$O, una cantidad vinculada a la regla de suma $m^{-1}$ de la función de respuesta, y comparar estas predicciones teóricas con los datos experimentales disponibles derivados de mediciones de secciones eficaces de fotoneutrones.

Metodología
El estudio emplea el enfoque de la Transformada Integral de Lorentz (LIT) combinado con el método de Clúster Acoplado (LIT-CC). Este marco evita el cálculo directo de los estados continuos mediante la transformación de la función de respuesta dependiente de la energía $R(\omega)$ en un problema de tipo estado ligado, a través de la convolución con un núcleo lorentziano $L(\sigma, \Gamma)$.

Componentes metodológicos clave:

• Marco de Muchos Cuerpos: El estado fundamental $|\Psi_0\rangle$ se parametriza utilizando el ansatz exponencial $|\Psi_0\rangle = e^{\hat{T}}|\Phi_0\rangle$, donde $|\Phi_0\rangle$ es un estado de referencia de Hartree-Fock esférico. El operador de clúster $\hat{T}$ se expande para incluir excitaciones de un solo y dobles (CCSD) y las principales excitaciones triples (CCSDT-1).
• Solución LIT: La LIT se calcula como la norma de los pseudoestados auxiliares $|\Psi_R(z)\rangle$ y $\langle\Psi_L(z^*)|$, que son soluciones de ecuaciones de tipo Schrödinger que involucran un Hamiltoniano transformado por similitud $\bar{H}$ y un operador de excitación $\bar{\Theta}$. Estos estados se resuelven utilizando un ansatz de Ecuación de Movimiento (EOM) dentro del espacio de configuración $np$-$nh$, utilizando el algoritmo de Lanczos no simétrico.
• Interacciones: Los cálculos utilizan dos Hamiltonianos de la teoría de campo efectivo quiral que contienen tanto fuerzas de dos como de tres nucleones: NNLOsat y $\Delta$NNL$O_{G}$(394). También se realiza una comparación con la interacción Entem-Machleidt EM(500) (sin fuerzas de tres nucleones) utilizando coeficientes de Lanczos publicados previamente.
• Esquemas de Truncamiento: Los autores comparan los resultados entre dos niveles de truncamiento: "D" (CCSD para el estado fundamental y el operador de transición) y "T-1" (CCSDT-1 para el estado fundamental, con el operador de transición aproximado mediante las amplitudes de CCSD).
• Reglas de Suma: En lugar de invertir la LIT para reconstruir $R(\omega)$, los autores calculan directamente los momentos de la función de respuesta. La polarizabilidad dieléctrica $\alpha_D$ está relacionada con el momento $m^{-1}$.

Resultados Clave

• Distribución de la Fuerza Dipolar: Ambas interacciones quirales predicen una concentración de la fuerza dipolar cerca de los 10 MeV, cerca de la energía de separación de neutrones (6.9 MeV), lo cual es consistente con el modo de baja energía observado experimentalmente. Una segunda concentración aparece en 20–25 MeV, correspondiente a la Resonancia Dipolar Gigante (GDR).
• Dependencia del Truncamiento: Al comparar los esquemas "D" y "T-1", se observa que la inclusión de correlaciones triples (T-1) desplaza ligeramente los picos de la GDR hacia energías más altas y reduce las reglas de suma dipolares. Sin embargo, la distribución general permanece similar, lo que sugiere que los resultados son robustos frente al truncamiento de muchos cuerpos.
• Valores de Polarizabilidad:
  • Para el rango de energía hasta el límite experimental ($\epsilon = 24.9$ MeV), las predicciones teóricas ($\alpha_D \approx 0.56 - 0.61$ fm$^3$) muestran un buen acuerdo con el valor experimental de $0.41(13)$ fm$^3$.
  • Al extrapolar a energía infinita ($\epsilon \to \infty$), los valores teóricos de $\alpha_D$ aumentan a aproximadamente $0.83 - 0.86$ fm$^3$. Los autores atribuyen la discrepancia entre el valor teórico completo y el valor experimental (que se detiene en 24.9 MeV) a que los datos experimentales omiten componentes de emisión de partículas cargadas (por ejemplo, emisión de protones) que contribuyen a la GDR.
• Dependencia de la Interacción: Los resultados de NNLOsat y $\Delta$NNL$O_{G}$(394) son consistentes entre sí y con cálculos previos. Se encuentra que la inclusión de fuerzas de tres nucleones y correlaciones de muchos cuerpos de orden superior (triples) mejora la descripción de la distribución de fuerza de baja energía en comparación con los cálculos que utilizan únicamente fuerzas de dos nucleones (EM(500)).

Significancia y Reivindicaciones
El artículo sostiene que el método LIT-CC proporciona un marco ab initio robusto para estudiar las reglas de suma electromagnética en isótopos ricos en neutrones como $^{22}$O. La principal significancia radica en demostrar que:

1. El método puede predecir con éxito la polarizabilidad dieléctrica y reproducir la respuesta dipolar de baja energía observada en experimentos.
2. La inclusión de fuerzas de tres nucleones y de correlaciones de muchos cuerpos de orden superior (triples) es esencial para una descripción precisa de la distribución de la fuerza.
3. El acuerdo entre la teoría y el experimento en la región de baja energía valida el enfoque LIT-CC para núcleos de capa abierta, a pesar de los desafíos planteados por los datos experimentales incompletos (rango de energía limitado y canales de decaimiento no resueltos).

Los autores concluyen que el trabajo en curso se centrará en la reconstrucción de la función de respuesta completa $R(\omega)$ mediante la inversión de la LIT para permitir comparaciones directas con las secciones eficaces experimentales tanto en núcleos de capa cerrada como de capa abierta.