Efficient calculation of two-neutrino double-beta-decay nuclear matrix elements

Resumen Técnico: Cálculo Eficiente de los Elementos de Matriz Nucleares de la Doble Beta-Desintegración de Dos Neutrinos

Planteamiento del Problema
Los elementos de matriz nucleares (NME, por sus siglas en inglés) fiables son fundamentales para interpretar los experimentos de desintegración beta-doble ($\beta\beta$) y para vincular las vidas medias medidas con la física fundamental de la interacción débil. Mientras que el modo de dos neutrinos ($2\nu\beta\beta$) es un proceso del Modelo Estándar observado en varios núcleos, el modo sin neutrinos ($0\nu\beta\beta$) sigue siendo la firma principal de la violación de la conservación del número leptónico. Los análisis recientes de modelos de capas mediante métodos estadísticos indican una fuerte correlación entre los NME de $2\nu\beta\beta$ y $0\nu\beta\beta$, lo que sugiere que el cálculo preciso de los primeros puede restringir los segundos.

Sin embargo, calcular los NME de $2\nu\beta\beta$ ($M_{2\nu}$) presenta un desafío computacional significativo. El enfoque estándar requiere una sumatoria explícita sobre todos los estados $1^+$ en el núcleo impar-par. En espacios de modelos grandes (por ejemplo, núcleos pesados o espacios de valencia expandidos como $sd-pf$), el número de estados contribuyentes se vuelve prohibitivo, y la suma acumulativa puede converger de forma lenta o no monotónica. La diagonalización completa del Hamiltoniano de muchos cuerpos para obtener estos estados es a menudo computacionalmente inviable, particularmente para los experimentos de próxima generación que apuntan a sensibilidades de hasta $10^{28}$ años.

Metodología
El artículo presenta una implementación mejorada del método de la función de fuerza basado en iteraciones de Lanczos para calcular $M_{2\nu}$ sin requerir la diagonalización completa. El enfoque sigue el formalismo de la Ref. [24] y utiliza una prescripción numérica detallada en la Ref. [25].

1. Iteraciones de Lanczos: El método inicia las iteraciones de Lanczos utilizando estados "puerta" de Gamow-Teller (GT) construidos aplicando el operador GT a los estados fundamentales inicial ($0^+_i$) o final ($0^+_f$).
2. Matriz Tridiagonal: En lugar de almacenar los autovectores completos del Hamiltoniano de muchos cuerpos, el algoritmo construye una matriz de Lanczos tridiagonal mucho más pequeña. El NME se aproxima sumando los autovalores y autovectores de esta pequeña matriz.
3. Coeficientes: Los coeficientes de expansión ($g^\pm_m$) se derivan de los vectores de Lanczos y de los denominadores de energía asociados con los estados intermedios.
4. Eficiencia: Este enfoque evita los costes de almacenamiento y procesamiento de la diagonalización del Hamiltoniano completo. Rastrea la convergencia de forma iterativa utilizando los propios vectores de Lanczos, ofreciendo un incremento de eficiencia computacional de aproximadamente un factor de dos en comparación con implementaciones anteriores de la función de fuerza.
5. Extensión a Órdenes Superiores: El marco se extiende para calcular NME de órdenes superiores ($M_{2\nu-3}$ y $M_{2\nu-5}$) requeridos para tratamientos de factores de espacio de fase expandidos mediante series de Taylor y observables relacionados de captura electrónica de dos neutrinos. Estos términos involucran potencias más altas del denominador de energía, lo que los hace más sensibles a los estados de baja energía.

Contribuciones Clave y Resultados
El autor aplica este método a varios isótopos experimentalmente importantes ($^{48}\text{Ca}$, $^{76}\text{Ge}$, $^{82}\text{Se}$, $^{128,130}\text{Te}$, y $^{136}\text{Xe}$) utilizando Hamiltonianos efectivos estándar (por ejemplo, GXPF1A, KB3G, JUN45, GCN2850, GCN5082, SVD) en sus correspondientes espacios de valencia.

• Evaluación de Referencia (Benchmarking): El método es rigurosamente probado comparando las "sumas acumuladas" obtenidas de la expansión de Lanczos contra sumas explícitas sobre cientos o miles de estados $1^+$ intermedios.
  • Para $^{48}\text{Ca}$ (capa-$pf$), el resultado de 25 iteraciones de Lanczos reproduce la suma directa sobre 1000 estados, capturando el ascenso rápido, los cambios en las regiones de fuerza y la meseta final estable.
  • Para núcleos más pesados como $^{82}\text{Se}$ y $^{136}\text{Xe}$, donde la diagonalización directa de miles de estados es difícil, el método de Lanczos reproduce el comportamiento de convergencia de las sumas directas (por ejemplo, sobre 2000 estados para $^{136}\text{Xe}$).
  • El método captura con éxito las contribuciones de la región de la resonancia gigante de GT y maneja las cancelaciones no triviales típicas de los elementos de matriz de $2\nu\beta\beta$.
• Términos de Orden Superior: Las pruebas de referencia para $M_{2\nu-3}$ y $M_{2\nu-5}$ en $^{48}\text{Ca}$ muestran una excelente concordancia con las sumas explícitas. Esto confirma que los mismos coeficientes de Lanczos utilizados para el NME principal pueden generalizarse a potencias inversas más altas del denominador de energía.
• Factores de Atenuación (Quenching): El estudio determina los factores de atenuación GT óptimos ($q_{opt}$) para varios Hamiltonianos con el fin de reproducir los elementos de matriz empíricos. Los resultados resaltan que estos factores son calibraciones fenomenológicas dependientes de la interacción específica y el espacio de valencia, más que constantes universales.

Significado y Reivindicaciones
El artículo afirma que el método de la función de fuerza de Lanczos propuesto proporciona una solución práctica y eficiente para evaluar $M_{2\nu}$ en espacios de modelos donde la suma directa es computacionalmente costosa o imposible.

• Precisión: El método preserva la precisión de la suma explícita donde es posible realizar comparaciones de referencia, incluso en casos que involucran cancelaciones significativas y contribuciones de estados de alta energía.
• Utilidad para $0\nu\beta\beta$: Dada la fuerte correlación estadística entre los NME de $2\nu\beta\beta$ y $0\nu\beta\beta$, los cálculos fiables y eficientes de $2\nu\beta\beta$ pueden ayudar a identificar Hamiltonianos efectivos y renormalizaciones de operadores que sean favorables para las predicciones de $0\nu\beta\beta$.
• Escalabilidad: La eficiencia del método permite la exploración de espacios de modelos ampliados (por ejemplo, incluyendo excitaciones de $2\hbar\omega$) que anteriormente eran inaccesibles mediante la diagonalización directa.
• Perspectiva Futura: El autor sugiere que este marco puede servir como un banco de pruebas complementario a gran escala para los cálculos ab initio emergentes, proporcionando muestras estadísticas y restricciones fenomenológicas derivadas de las vidas medias medidas.

Este trabajo no propone nuevas firmas experimentales ni pretende resolver el problema de la $0\nu\beta\beta$ directamente, sino que posiciona el cálculo eficiente de los NME de $2\nu\beta\beta$ como una herramienta necesaria para reducir las incertididades teóricas en la búsqueda más amplia de la desintegración sin neutrinos.