Elastic deuteron-deuteron scattering within Nuclear Lattice Effective Field Theory

Resumen Técnico: Dispersión Elástica de Deuterón-Deuterón dentro de la Teoría de Campo Efectivo de Red Nuclear

Planteamiento del Problema
Las predicciones teóricas precisas para la Nucleosíntesis del Big Bang (BBN) requieren un conocimiento preciso de las tasas de reacción dominantes, específicamente las reacciones de deuterón-deuterón $d(d, n)^3\text{He}$ y $d(d, p)^3\text{H}$. Las discrepancias actuales en las predicciones de la abundancia primordial de deuterio ("anomalía del deuterón") dependen de la fuente de las tasas de reacción utilizadas (por ejemplo, NACRE II frente a PRIMAT). Para resolver esto, es necesaria una realización ab initio de la física de dispersión subyacente. El problema físico de cuatro nucleones en el sector $J^\pi = 2^+$ es un problema multicanal que involucra canales de $d+d$ y canales de reordenamiento ($p+^3\text{H}$, $n+^3\text{He}$). Cálculos previos utilizando métodos como el Método de Grupo Resonante (RGM) o las ecuaciones de Yakubovskii han arrojado diferencias significativas en los resultados, particularmente en el canal $S=2$ ($^5S_2$). Este artículo aborda el problema de dispersión de canal único $d+d$, $^5S_2$, como un paso fundamental hacia un cálculo de canales acoplados completo dentro de la Teoría de Campo Efectivo de Red Nuclear (NLEFT).

Metodología
Los autores emplean NLEFT, una formulación de red de la Teoría de Campo Efectivo Quiral (ChEFT), donde los nucleones se colocan en una red cúbica con condiciones de contorno periódicas. El cálculo utiliza interacciones quirales en el Orden Siguiente a Siguiente al Orden Principal (N3LO).

• Acoplamiento de la Función de Onda: Para manejar eficientemente el Hamiltoniano N3LO de alta fidelidad, los autores aplican una transformación unitaria de rango finito. Esto mapea el complejo Hamiltoniano N3LO a un Hamiltoniano más simple ($H_S$) a distancias cortas ($r < R$) preservando la física original a distancias largas. El Hamiltoniano simple incluye energía cinética, una interacción de contacto de corto alcance y el intercambio de un pión. La corrección ($\Delta H$) se trata de forma perturbativa, asegurando una convergencia rápida y suprimiendo las oscilaciones de signo de Monte Carlo.
• Método de Proyección Adiabática: Para aumentar la eficiencia computacional, los estados iniciales de dos cústers se construyen a partir de dos cústers de deuterón (cada uno con proyección de espín $m_d = +1$), seleccionando de forma única el canal de espín-quinteto. El movimiento relativo se proyecta sobre armónicos esféricos ($\ell=0, m_\ell=0$) para aislar el canal $^5S_2$. El sistema evoluciona en tiempo euclidiano utilizando una matriz de transferencia para filtrar los estados de menor energía.
• Regularización de la Matriz de Norma: Un desafío técnico crítico surge porque la base de cústers radiales desarrolla autovalores de la matriz de norma pequeños en tiempos de proyección euclidea largos, lo que provoca un mal condicionamiento. Los autores investigan y comparan dos procedimientos de estabilización:
  1. Regularización de Tikhonov: Añadir un parámetro de regularización a la matriz de norma para evitar la amplificación de los modos propios mal resueltos.
  2. Proyección de Modos de Norma: Proyectar el Hamiltoniano sobre los autovectores de la matriz de norma con comportamientos físicos bien definidos (autovalores por encima de un corte).
• Extracción de Desplazamientos de Fase: Los autores utilizan una estrategia de potencial auxiliar. Se añade un potencial auxiliar real cerca del límite radial para desplazar los autovalores de energía discretos de forma continua. Al ajustar las funciones de onda resultantes a formas asintóticas modificadas por Coulomb, extraen los desplazamientos de fase corregidos por Coulomb.
• Análisis de Rango Efectivo: Los desplazamientos de fase extraídos se ajustan a una expansión de rango efectivo modificada por Coulomb para determinar la longitud de dispersión ($^5a_{dd}$) y el rango efectivo ($^5r_{dd}$).

Resultos Clave

• Consistencia de la Estabilización: Tanto la regularización de Tikhonov como la proyección de modos de norma producen desplazamientos de fase corregidos por Coulomb consistentes dentro de las incertidumbres estadísticas y numéricas.
• Comportamiento del Desplazamiento de Fase: Los desplazamientos de fase calculados son más negativos a medida que aumenta el momento en comparación con cálculos previos (incluyendo RGM, Yakubovskii y métodos variacionales de estado ligado). Esto indica una repulsión efectiva más fuerte en el canal $^5S_2$.
• Parámetros de Dispersión: El análisis de rango efectivo modificado por Coulomb arroja:
  • Longitud de dispersión: $^5a_{dd} = (12.96 \pm 0.26)$ fm.
  • Rango efectivo: $^5r_{dd} = (3.62 \pm 0.79)$ fm.
    La longitud de dispersión es sustancialmente mayor que los valores reportados en estudios previos (por ejemplo, $7.8 \pm 0.3$ fm en la Ref. [19]).
• Verificaciones Sistemáticas: Los resultados son consistentes con un cálculo en una red más grande ($L=21.1$ fm), lo que sugiere que los efectos de volumen finito son insignificantes, aunque las incertidumbres estadísticas en esa verificación específica fueron grandes. Se encontró que las interacciones de tres nucleones tienen solo un efecto pequeño en los desplazamientos de fase.

Significado y Reivindicaciones
El artículo establece el primer referente de red nuclear para la dispersión de deuterón-deuterón en el canal $^5S_2$. Los autores declaran explícitamente que estos resultados proporcionan una base para futuros cálculos de canales acoplados de reacciones inducidas por deuterón relevantes para la BBN. No pretenden haber resuelto el problema completo de canales acoplados ni determinado la matriz de dispersión completa para el sector $J^\pi=2^+$; más bien, han restringido la base asintótica al canal único de $d+d$ $^5S_2$. El trabajo demuestra la viabilidad de tratar este canal con N3LO en NLEFT e identifica los requisitos numéricos (específicamente respecto a la estabilización de la matriz de norma) para extender el enfoque a canales de reacción acoplados más complejos. Los autores señalan que el origen de las discrepancias con la literatura previa no puede identificarse de forma inequívoca, pero probablemente se deba a diferencias en las interacciones nucleares, los espacios de canales y las aproximaciones numéricas.