Radiative corrections to two-neutrino double-beta decay

Using heavy-nucleus effective field theory, this paper derives a universal "double-weak Sirlin function" for two-neutrino double-beta decay that significantly distorts electron energy and angular spectra, necessitating the inclusion of these radiative corrections in high-precision Standard Model tests and nuclear structure extractions.

The Gist

Imagine you are trying to listen to a very faint, specific whisper in a crowded, noisy room. That whisper is a rare event in the universe called double-beta decay, where a nucleus changes its identity by spitting out two electrons. Scientists are obsessed with this event because it might hold the key to understanding why the universe is made of matter instead of antimatter.

However, there's a problem. The "crowd" is filled with background noise. In this case, the noise is the standard version of the decay (called two-neutrino double-beta decay), which happens all the time and looks almost exactly like the rare event we are hunting for. To hear the rare whisper, we need to understand the background noise with extreme precision—down to the last decimal point.

This paper is about cleaning up the "noise" so we can hear the "whisper" better. Here is the breakdown of what the scientists did, using some everyday analogies.

1. The Old Map vs. The New GPS

For a long time, physicists calculated the energy of the electrons coming out of this decay using a "rule of thumb." They assumed that because two electrons are coming out, you could just take the math for one electron, double it, and add them together.

Think of it like this: If you are baking a cake and you know how much sugar one egg needs, you might assume two eggs need exactly double the sugar.

• The Old Way: "Two electrons = Two single-electron calculations added together."
• The New Discovery: The authors found that this is wrong. The two electrons don't just sit there; they interact with each other and with the nucleus in a complex dance. It's like realizing that when you put two eggs in a bowl, they don't just sit side-by-side; they mix, change the texture, and affect how the sugar dissolves.

The authors created a new, universal "GPS" (a mathematical function they call the Double-Weak Sirlin Function) that accounts for how the two electrons talk to each other and to the nucleus. They found that the old "double it" method was off by a significant amount.

2. The "Ghost" Photons

In the quantum world, particles are constantly exchanging invisible messengers called photons. Even when no light is actually seen, these "ghost" photons are zipping around, tweaking the energy of the electrons.

The scientists calculated these "ghost" effects for the first time in this specific double-decay scenario.

• The Analogy: Imagine two runners (the electrons) on a track. The old math assumed they were running on a flat, empty track. The new math realizes that the track is actually a bumpy, windy field, and the runners are occasionally bumping into invisible wind gusts (photons) that push them slightly faster or slower.
• The Result: These wind gusts change the shape of the energy curve. The paper shows that this change is huge—it's about the same size as the biggest "structural" changes caused by the nucleus itself.

3. Why This Matters: The "Tuning Fork" Problem

Scientists use these decay measurements to tune their "tuning forks" (theoretical models of the atomic nucleus). If they get the background noise wrong, they will tune the fork incorrectly.

• The Conflict: Recently, experiments (like the CUORE experiment mentioned in the paper) found that their measurements didn't quite match the theoretical predictions. They were confused.
• The Solution: The authors say, "You're confused because you forgot the 'wind gusts' (radiative corrections)!" When you add their new, precise calculation of these ghost photons, the distortion in the data changes.
• The Twist: In some cases, the new "wind gust" correction actually cancels out the other distortions. It's like if you were trying to measure a wobbly table, but you forgot that the floor was also tilting. Once you account for the floor, the table looks perfectly level.

4. The "Bonus" Event: The Flash

The math also predicts a rare side effect: sometimes, instead of just two electrons, the nucleus spits out two electrons and a flash of light (a photon).

• The Analogy: It's like a firework that usually just shoots sparks, but occasionally shoots a spark and a small flash of light.
• The Prediction: The authors calculated how often this happens. It's rare (about 1 in 100 to 1 in 1,000 events), but with the massive new detectors being built, we might finally be able to see these "flash" events. This would be a direct confirmation of their new math.

The Bottom Line

This paper is a massive upgrade to the "instruction manual" for understanding double-beta decay.

1. Stop doubling: You can't just add two single-decay calculations together; you have to calculate the two-electron system as a whole.
2. The noise is loud: The "ghost photon" corrections are just as important as the nuclear structure itself. Ignoring them is like trying to weigh a feather on a scale that hasn't been zeroed out.
3. Re-do the math: Any recent scientific papers that tried to extract nuclear secrets from these decays need to be re-run with this new "Double-Weak Sirlin Function" included.

By fixing these calculations, scientists can finally stop guessing about the background noise and start listening clearly to the rare, universe-changing whispers they are hunting for.

Technical Summary

1. Problem Statement

Two-neutrino double-beta decay ($2\nu\beta\beta$) is a Standard Model (SM) process that serves as both an irreducible background for neutrinoless double-beta decay ($0\nu\beta\beta$) searches and a precision observable for testing nuclear structure and the SM itself. Current and next-generation experiments aim for sub-percent precision in measuring $2\nu\beta\beta$ spectral shapes (electron energy and angular distributions).

However, theoretical descriptions at this precision level are incomplete because they lack a rigorous treatment of electromagnetic radiative corrections specific to the double-weak process.

• Current Approximation: Previous analyses estimated these corrections by simply summing two independent single-$\beta$ decay Sirlin functions.
• The Flaw: This approximation ignores the unique kinematics of $2\nu\beta\beta$, specifically the spectrum of intermediate nuclear excitations and the correlations between the two emitted electrons. It fails to account for photon exchanges between the two electrons and interactions with intermediate nuclear states.

2. Methodology

The authors employ Heavy-Nucleus Effective Field Theory (EFT) to perform the first explicit calculation of ultrasoft radiative corrections for double-weak processes.

• Framework: The EFT includes fields for the initial ($0^+$), final ($0^+$), and intermediate ($1^+$) nuclear states. It integrates out nucleon degrees of freedom, organizing the expansion in the ratio of lepton energy to the Fermi momentum ($E_e/k_F$).
• Lagrangian: The leading-order Lagrangian includes kinetic terms for the nuclei and weak decay vertices coupling the nuclei to leptons (electrons and neutrinos).
• Diagrammatic Calculation: The authors calculate all virtual and real photon diagrams at order $\mathcal{O}(\alpha)$ (Fig. 1 in the paper). Crucially, this includes topologies absent in single-$\beta$ decay:
  • Photons attached to the intermediate nuclear state.
  • Photon exchange between the two outgoing electrons.
• Renormalization: Ultraviolet (UV) divergences are handled via dimensional regularization and renormalization of the vector coupling $g_V$. Infrared (IR) divergences cancel between virtual loops and real photon emission (bremsstrahlung).
• Expansion: The results are derived for generic nuclear excitation energies ($\omega_n$) and expanded in the limit $\omega_n \gg \epsilon$ (where $\epsilon$ is the lepton energy sum), allowing for a universal expression independent of specific nuclear matrix elements (NMEs) at leading order.

3. Key Contributions

The paper introduces several novel theoretical constructs and calculations:

1. The "Double-Weak Sirlin Function":

   • The authors derive a universal radiative correction factor for $2\nu\beta\beta$, analogous to the Sirlin function in single-$\beta$ decay.
   • Unlike the single-$\beta$ Sirlin function (which depends only on electron energy and $Q$-value), this new function depends on both electron energies and their relative angle ($y_{12}$).
   • It explicitly accounts for the intermediate nuclear excitation spectrum.

2. Explicit Calculation of Unique Topologies:

   • The work quantifies corrections arising from photon exchange between the two electrons and photon coupling to the intermediate nucleus, which were previously uncalculated.

3. Separation of Universal and Nuclear-Structure Dependent Terms:

   • The authors separate the radiative corrections into a universal part (independent of NMEs at leading order) and a nuclear-structure dependent part.
   • They demonstrate that the nuclear-structure dependent component is negligible at current experimental sensitivities, while the universal part is significant.

4. Key Results

• Deviation from Single-$\beta$ Approximation:

  • The "double-weak Sirlin function" differs significantly from the naive sum of two single-$\beta$ Sirlin functions across the entire phase space.
  • The naive approximation overestimates the total radiative corrections. For example, in $^{76}\text{Ge}$, the full calculation yields a different spectral shape compared to the sum-of-two approximation.

• Magnitude of Distortions:

  • The radiative corrections induce distortions in the electron energy spectrum and angular distributions that are comparable in size to the leading nuclear-structure correction parameterized by the ratio of NMEs, $\xi_{31} = M^{(-3)}_{GT} / M^{(-1)}_{GT}$.
  • In some isotopes (e.g., $^{76}\text{Ge}$), the radiative distortion almost completely cancels the distortion induced by $\xi_{31}$ over a wide energy range. In others (e.g., $^{100}\text{Mo}$), the effects nearly cancel in the angular distribution.

• Impact on $\xi_{31}$ Extraction:

  • Because radiative effects mimic or cancel nuclear structure effects, extracting $\xi_{31}$ from precision $2\nu\beta\beta$ data without including these corrections leads to biased results.
  • The authors suggest that recent extractions of $\xi_{31}$ (e.g., from CUORE) that did not include these specific double-weak corrections may need revision.

• Radiative Decay Branching Ratio ($2\nu\beta\beta + \gamma$):

  • The calculation predicts the branching ratio for the radiative mode where a photon is emitted.
  • For a photon energy cut $x_\gamma = E_\gamma/Q \sim 0.01$, the branching ratio is $\mathcal{O}(10^{-2} - 10^{-3})$. This suggests the mode is observable in high-statistics next-generation experiments, though high-energy photons are rare.

5. Significance and Implications

• Nuclear Structure Constraints: The results imply that future high-precision measurements of $2\nu\beta\beta$ cannot be used to constrain nuclear matrix elements (crucial for $0\nu\beta\beta$ predictions) unless the double-weak radiative corrections are consistently included in the analysis.
• Standard Model Tests: The inclusion of these corrections is essential for using $2\nu\beta\beta$ as a precision probe of electroweak dynamics in nuclei and for searching for physics beyond the Standard Model.
• Experimental Guidance: The paper provides the theoretical framework and code (via a Python Jupyter notebook) for experimental collaborations (like LEGEND, nEXO, CUPID) to incorporate these corrections into their spectral shape analyses. It also motivates dedicated searches for the $2\nu\beta\beta + \gamma$ channel.
• Theoretical Advancement: This work establishes a rigorous EFT framework for double-weak processes, paving the way for calculating higher-order corrections (e.g., $\mathcal{O}(\alpha^2 Z)$) and other subleading effects.

In conclusion, the paper demonstrates that the "double-weak Sirlin function" is a critical, non-negligible component of $2\nu\beta\beta$ phenomenology, necessitating a revision of how nuclear structure parameters are extracted from experimental data.