Lattice Calculation of Short-Range Contributions to Neutrinoless Double-Beta Decay $\pi^-\to\pi^+ ee$ at Physical Pion Mass

Este artículo presenta un cálculo de QCD en retículo de los elementos de matriz $\pi^- \to \pi^+ ee$ a la masa física del pión utilizando fermiones de pared de dominio, empleando un método novedoso para restar efectos de alrededor del mundo y una renormalización no perturbativa para reducir las incertidumbres y conciliar las discrepancias en estudios previos de contribuciones de corto alcance al doble decaimiento beta sin neutrinos.

Lo esencial

El panorama general: Resolver un misterio cósmico

Imagina que el universo es un rompecabezas gigante, y una de las piezas faltantes es la naturaleza de los neutrinos. Estas son partículas diminutas, parecidas a fantasmas, que atraviesan todo a gran velocidad. Los científicos quieren saber: ¿Son fermiones "Dirac" (como los electrones normales, donde una partícula y su antipartícula son distintas) o fermiones "Majorana" (donde una partícula es su propia antipartícula)?

La única manera de resolver este misterio es observar un evento muy raro llamado desintegración doble beta sin neutrinos. Es como ver dos átomos cambiar espontáneamente en átomos diferentes y escupir dos electrones, pero sin liberar ningún neutrino. Si vemos que esto sucede, prueba que los neutrinos son sus propias antipartículas.

El problema: Una señal ruidosa

Para predecir si ocurrirá este evento, los físicos tienen que realizar cálculos matemáticos complejos. Dividen el cálculo en dos partes:

1. La parte de larga distancia: Como un susurro que viaja a través de una habitación.
2. La parte de corta distancia: Como un grito que ocurre justo al lado de tu oído.

Este artículo se centra en la parte de corta distancia. Específicamente, están calculando cómo dos piones (partículas formadas por quarks) interactúan para producir dos electrones. Piensa en esto como medir el "volumen" de ese grito.

El conflicto: Dos equipos diferentes de científicos habían intentado medir previamente este "volumen" utilizando supercomputadoras (llamadas QCD de red). Sin embargo, sus resultados discrepaban en un factor de dos. Era como si un equipo dijera que el grito era de 60 decibelios y el otro dijera que era de 120 decibelios. Esta enorme discrepancia hacía difícil confiar en las predicciones para el misterio de los neutrinos.

La solución: Una nueva forma de limpiar los datos

Los autores de este artículo decidieron realizar su propio experimento para zanjar la cuestión. Utilizaron una supercomputadora masiva para simular el mundo subatómico. Pero se enfrentaron a un problema técnico específico: "Efectos de dar la vuelta al mundo".

La analogía: Imagina que estás grabando una conversación en una habitación pequeña y con eco, con una pared circular. Si aplaudes, el sonido viaja hacia adelante, golpea la pared, da la vuelta a la habitación y regresa a ti desde atrás. En la simulación por computadora, la "habitación" es la cuadrícula del espacio-tiempo. Debido a que la cuadrícula es finita, las partículas pueden viajar todo el camino alrededor del bucle e interferir con la medición, creando un confuso "eco" que arruina los datos.

La innovación: Los métodos anteriores intentaban adivinar cómo cancelar estos ecos. Este equipo inventó un nuevo método de sustracción.

• En lugar de adivinar, aislaron la señal del "eco" directamente de los datos.
• Calcularon exactamente qué tan fuerte era el eco y lo restaron de la señal principal.
• El resultado: El "ruido" desapareció, dejando una señal limpia y estable (un "plateau") en la que podían confiar.

La verificación: Comprobando la regla

Para asegurarse de que su nuevo método no estaba defectuoso, verificaron su trabajo contra un estándar conocido. Calcularon un valor específico (llamado "parámetro de bolsa") que había sido medido por otros equipos antes.

• Su resultado coincidió perfectamente con el estándar confiable.
• Cuando compararon su resultado con el equipo que tenía la diferencia de "factor de dos", descubrieron que sus números eran exactamente el doble de los números del otro equipo.
• La conclusión: Resulta que el otro equipo probablemente utilizó una "regla" ligeramente diferente (convención de normalización) para sus mediciones. Una vez que se tiene en cuenta esa diferencia, los puntos de datos en realidad se alinean muy bien. El método de los autores confirma que su cálculo es correcto y resuelve la confusión.

El resultado final

El equipo calculó con éxito la contribución de "corto alcance" al proceso de desintegración doble beta sin neutrinos con mucha mayor precisión que antes.

• Eliminaron los "ecos" (efectos de dar la vuelta al mundo) que estaban arruinando los datos anteriores.
• Utilizaron dos "lentes" matemáticas diferentes (esquemas de renormalización) para asegurar que sus matemáticas eran sólidas.
• Proporcionaron un número definitivo y de alta precisión que ayuda a los físicos a predecir si eventualmente veremos esta desintegración rara en experimentos de la vida real.

En resumen: Construyeron un mejor microscopio, limpiaron el ruido estático y confirmaron que el desacuerdo anterior era simplemente una cuestión de usar diferentes cintas métricas. Ahora, la comunidad científica tiene un número confiable para ayudar a resolver el misterio de la masa de los neutrinos.

Resumen técnico

A continuación se presenta un resumen técnico detallado del artículo "Cálculo en retículo de las contribuciones de corto alcance a la desintegración doble beta sin neutrinos $\pi^-\to\pi^+ee$ a la masa física del pión" (CERN-TH-2025-132).

1. Planteamiento del Problema

La desintegración doble beta sin neutrinos ($0\nu\beta\beta$) es un proceso crítico para determinar si los neutrinos son fermiones de Majorana y para sondear la violación del número leptónico (LNV). La tasa de desintegración depende de los elementos de matriz nucleares (NME). En la Teoría de Campo Efectivo del Modelo Estándar (SMEFT), las contribuciones de corto alcance a la $0\nu\beta\beta$ surgen de operadores de dimensión 9. Estas contribuciones están mediadas por el proceso $\pi^- \to \pi^+ ee$ a través del intercambio de piones entre nucleones.

El Problema Central: Los cálculos previos de QCD en retículo de los elementos de matriz de corto alcance para $\pi^- \to \pi^+ ee$ (específicamente las Refs. [30] y [31]) reportaron resultados significativamente discrepantes. Para el parámetro de bolsa $B_\pi^{MS}$ a $\mu=3$ GeV, la Ref. [30] reportó $\approx 0.42$, mientras que la Ref. [31] reportó $\approx 0.20$—una discrepancia de aproximadamente un factor de dos. Esta inconsistencia obstaculiza predicciones teóricas fiables para las tasas de desintegración $0\nu\beta\beta$. Además, los estudios anteriores adolecieron de:

• Incertidumbres Sistemáticas: Efectos significativos de "alrededor del mundo" (ATW) debido a la propagación hacia atrás del pión en retículos temporales periódicos, lo que arruina la extracción de elementos de matriz del estado fundamental.
• Extrapolación Quiral: Algunos estudios dependieron de masas de pión no físicas, requiriendo una extrapolación al punto físico, introduciendo una dependencia adicional del modelo.

2. Metodología

Los autores realizaron un cálculo de QCD en retículo de alta precisión utilizando ensambles de Fermiones de Pared de Dominio (DWF) generados por la Colaboración RBC/UKQCD.

A. Configuración del Retículo

• Ensamble: Se utilizaron ensambles $N_f = 2+1$, incluyendo dos con masas físicas de pión (48I y 64I) con diferentes espaciados de retículo ($a^{-1} \approx 1.73$ y $2.36$ GeV) para permitir una extrapolación directa al continuo. También se utilizaron ensambles de masa no física para verificaciones cruzadas.
• Operadores: Se calcularon elementos de matriz para cinco operadores de cuatro quarks de dimensión 9 ($O_1, O_2, O_3, O'_1, O'_2$) relevantes para la transición $\pi^- \to \pi^+ ee$.

B. Mitigación de Efectos "Alrededor del Mundo" (ATW)

Una innovación mayor de este trabajo es un nuevo método para manejar los efectos ATW, que surgen cuando los piones se propagan a través del límite temporal periódico, contaminando las funciones de correlación.

• Análisis: Los autores analizaron la estructura temporal de las funciones de correlación de tres puntos, identificando contribuciones del diagrama deseado (A) y de los diagramas ATW (B, C, D).
• Nuevo Método de Sustracción: En lugar de depender de métodos de razón (R1 o R2) que solo suprimen parcialmente los efectos ATW, los autores propusieron una técnica de sustracción directa.
  • Identificaron regiones temporales donde la contribución ATW de un solo pión cruzando el límite (Diagrama B) domina.
  • Reconstruyeron la contribución ATW directamente a partir de los datos del retículo utilizando separaciones temporales específicas ($t_{\pi\pi}, t_\pi$).
  • Fórmula: El elemento de matriz corregido se obtiene restando el término ATW reconstruido:
    $$O^{(sub)}_i(t) = N_\pi^{-2} e^{2m_\pi t} \left[ C^i_3(t,t) - 2 C^i_3(t_{\pi\pi}, T - t_{\pi\pi} - t_\pi) e^{-m_\pi(T-t_\pi-2t_{\pi\pi})} \right]$$
  • Este método restauró con éxito mesetas estables en las gráficas de masa efectiva, reduciendo significativamente los errores sistemáticos en comparación con los métodos de razón anteriores.

C. Renormalización

• Esquema: Se realizó una renormalización no perturbativa utilizando el método RI/SMOM (Roma-Southampton) con cinemática no excepcional.
• Esquemas: Se emplearon dos esquemas para estimar los errores de truncamiento perturbativo: $(\gamma_\mu, \gamma_\mu)$ y $(\not{q}, \not{q})$.
• Emparejamiento: Los resultados se convirtieron al esquema $\overline{MS}$ a $\mu = 3$ GeV utilizando coeficientes de emparejamiento de un bucle. La diferencia entre los dos esquemas RI/SMOM se utilizó para estimar la incertidumbre sistemática derivada del truncamiento perturbativo.

3. Resultados Clave

A. Elementos de Matriz Desnudos

El estudio extrajo elementos de matriz desnudos $\langle \pi^+ | O_i | \pi^- \rangle$ a la masa física del pión.

• Verificación Cruzada con Parámetros de Bolsa: Los autores verificaron su normalización comparando el elemento de matriz de $O_3$ (relacionado con el parámetro de bolsa del mesón neutro $B_K$) con resultados revisados por FLAG anteriores. Los resultados coincidieron perfectamente con la Ref. [38] (que utilizó la misma convención de normalización) pero fueron aproximadamente el doble de grandes que los resultados en la Ref. [31].
• Discrepancia de Normalización: La relación de los resultados de los autores con la Ref. [31] es consistentemente $\approx 2.0$ para los cinco operadores. Esto sugiere un error global de factor de normalización en la Ref. [31], en lugar de una discrepancia física.

B. Elementos de Matriz Renormalizados y Parámetro de Bolsa

Los resultados finales en el esquema $\overline{MS}$ ($\mu = 3$ GeV) son:

• Parámetro de Bolsa: $B_\pi^{MS}(3 \text{ GeV}) = 0.3769(24)_{stat}(76)_{PT}(1)_L$.
  • Este valor es consistente con la Ref. [30] ($0.421$) dentro de las incertidumbres, pero significativamente diferente de la Ref. [31] ($0.197$).
  • La relación con la Ref. [31] es $1.91(18)$, confirmando el problema de normalización de un factor de 2.
• Extrapolación al Continuo: Utilizando los ensambles de masa física (48I y 64I), los autores realizaron una extrapolación directa al continuo, eliminando las incertidumbres asociadas con la extrapolación quiral.
• Reducción de Incertidumbre: Las incertidumbres totales en este trabajo se redujeron significativamente en comparación con estudios anteriores debido a:
  1. Cálculo directo a la masa física del pión.
  2. Precisión estadística mejorada.
  3. Eliminación efectiva de errores sistemáticos ATW.
  4. Estimación robusta de los errores de truncamiento perturbativo.

4. Significado y Conclusión

• Resolución de la Discrepancia: Este trabajo proporciona una verificación cruzada independiente que respalda fuertemente la normalización utilizada en la Ref. [30] y la revisión FLAG, mientras identifica un probable error de normalización en la Ref. [31]. La discrepancia de un factor de 2 se resuelve como un problema de convención de normalización en lugar de un desacuerdo fundamental de física.
• Avance Metodológico: El método de sustracción propuesto para los efectos ATW es una técnica robusta que puede aplicarse a otros cálculos en retículo que involucren mesones ligeros en volúmenes periódicos, mejorando la fiabilidad de los futuros cálculos de elementos de matriz $0\nu\beta\beta$.
• Impacto en la Física: Al proporcionar elementos de matriz precisos, extrapolados al continuo y a masas físicas de pión, este trabajo reduce las incertidumbres de QCD no perturbativas en la predicción teórica de las tasas de desintegración $0\nu\beta\beta$. Esto es crucial para interpretar los resultados experimentales de las próximas generaciones de experimentos (por ejemplo, LEGEND, nEXO, CUPID) para determinar la naturaleza de la masa del neutrino y la escala de nueva física.

En resumen, este artículo establece un nuevo referente para las contribuciones de corto alcance a la desintegración doble beta sin neutrinos, resolviendo las inconsistencias anteriores mediante una metodología rigurosa, simulaciones con masas físicas y un tratamiento novedoso de los artefactos de volumen finito.