First-Principles Determination of the Proton-Proton Fusion Matrix Element from Lattice QCD

This paper presents a first-principles lattice QCD calculation of the proton-proton fusion matrix element at an unphysical pion mass, demonstrating the feasibility of the approach while highlighting that large uncertainties in two-nucleon scattering parameters currently limit the precision of the extracted low-energy constant $L_{1,A}$.

The Gist

Imagine the Sun as a giant, cosmic power plant. For billions of years, it has been running on a very specific, very slow fuel: proton-proton fusion. This is the process where two tiny hydrogen nuclei (protons) crash into each other and fuse to form a heavier nucleus (deuterium), releasing the energy that keeps us warm and alive.

The problem? This reaction is incredibly rare and happens so slowly that it's almost impossible to measure directly in a lab. Scientists have been trying to calculate exactly how this happens using the laws of physics, but it's like trying to predict the exact path of a single grain of sand in a hurricane.

This paper is a major step forward in solving that puzzle. Here is the story of what the researchers did, explained without the heavy math.

1. The "Digital Sandbox" (Lattice QCD)

To understand how protons fuse, you have to understand the "glue" that holds them together: the Strong Force. This force is governed by a theory called Quantum Chromodynamics (QCD). But QCD is so complex that you can't solve it with a pencil and paper.

Instead, the researchers built a digital sandbox called "Lattice QCD." Imagine the universe isn't smooth, but is actually a giant 3D grid (like a massive chessboard). They placed protons on this grid and used supercomputers to simulate how they interact.

• The Catch: They couldn't simulate the real universe yet. The "glue" (pions) in their simulation was too heavy, like trying to study a feather by weighing a brick. They were working at a "heavy pion mass" (432 MeV), which is not the real world, but it's a necessary stepping stone.

2. The "Ghost" Problem (Excited States)

When you look at a proton in this digital sandbox, it's not just a single, clean object. It's surrounded by a cloud of "noise" or "ghosts" (called excited states). These are temporary, jumpy versions of the proton that make it hard to see the real, calm proton underneath.

• The Old Way: Previous scientists used "pointy" tools to grab the protons. It was like trying to pick up a specific jellyfish in a tank full of them using a single needle. You often grabbed the wrong one or got tangled in the noise.
• The New Way: This team used "bi-local" operators. Imagine instead of a needle, you use a wide, soft net that can gently scoop up the proton from two different spots at once. This net is much better at ignoring the noisy ghosts and grabbing the "real" proton. This allowed them to get a much clearer picture of the energy levels.

3. The "Room Size" Problem (Finite Volume)

Here is the trickiest part. In their digital sandbox, the universe is tiny (a small box). In the real world, the universe is infinite.

When you put two protons in a small box, they bounce off the walls and hit each other more often than they would in an open field. This changes the result of the experiment. It's like trying to measure how fast two people can run a race, but you force them to run inside a tiny elevator. The walls affect their speed.

• The Solution: The researchers used a mathematical "magic spell" (called the Lellouch-Lüscher correction). This spell acts like a translator. It takes the messy results from the tiny, bouncy box and translates them into what would happen in the vast, open universe.
• The Challenge: Because the protons in their simulation were slightly heavier than real ones, they were very sensitive to the size of the box. The "translation" required knowing exactly how the protons bounced off each other, which was hard to measure precisely.

4. The "Two-Handed" Dance (One-body vs. Two-body)

When the protons fuse, they don't just act alone.

• The One-Handed Move: One proton acts on its own (like a solo dancer).
• The Two-Handed Move: The two protons interact with each other while fusing (like a dance partner helping the other spin).

The researchers wanted to measure the "Two-Handed" move. They found that the protons were indeed helping each other, but it was a very subtle effect.

• The Result: They calculated a number (called $L_{1,A}$) that describes how strong this "two-hand" cooperation is.
• The Verdict: Their number was 6.0, with a big "maybe" attached to it (±7.1). It's a bit like saying, "The speed limit is 60 mph, give or take 70." It's a huge range!

Why is this a big deal if the number is so uncertain?

You might think, "If the answer is so fuzzy, why publish?"

1. It's the First Time: This is the first time anyone has tried to calculate this specific fusion process from the very bottom up (from the quarks and gluons) without just guessing based on experiments. It's like building a bridge from scratch to see if it holds, rather than just copying an old blueprint.
2. The Method Works: They proved that their "soft net" (bi-local operators) and their "translation spell" (finite-volume corrections) work. They showed that the "two-hand" dance is real and significant.
3. The Future: The big uncertainty comes from the fact that their simulation used "heavy" protons. The next step is to run the simulation with the real weight of protons. Once they do that, the "noise" will go down, and the answer will become sharp and precise.

The Bottom Line

This paper is a proof of concept. The researchers successfully built a digital model of the Sun's most important reaction. They didn't get the perfect final answer yet because the simulation was still a bit "heavy" and the math was incredibly difficult, but they proved the path forward is clear.

They showed us that with better computers and better "nets" to catch the protons, we will soon be able to calculate exactly how the Sun shines, purely from the laws of physics, without needing to look at the Sun itself.

Technical Summary

Here is a detailed technical summary of the paper "First-Principles Determination of the Proton-Proton Fusion Matrix Element from Lattice QCD" by Wang et al.

1. Problem Statement

Proton-proton fusion ($pp \to d e^+ \nu_e$) is the initiating reaction of the proton-proton chain, governing stellar energy production and evolution. While the cross-section is crucial for the Standard Solar Model and neutrino physics, a first-principles determination of its matrix element from Quantum Chromodynamics (QCD) has remained a long-standing challenge.

Previous lattice QCD attempts (e.g., by the NPLQCD Collaboration) were limited to unphysical pion masses ($m_\pi \sim 806$ MeV) where two-nucleon systems form deeply bound states, allowing for simplified finite-volume corrections. However, at lighter, more physical pion masses, the system may be shallowly bound or in a scattering regime, where finite-volume effects and two-nucleon rescattering become substantial. Furthermore, extracting matrix elements for two-nucleon systems is hindered by the "signal-to-noise" problem and severe excited-state contamination.

2. Methodology

The authors present a lattice QCD calculation at a pion mass of $m_\pi \approx 432$ MeV using $2+1$ flavor domain-wall fermions. The study employs a systematic framework to address the specific challenges of weak two-nucleon reactions:

• Interpolating Operators:

  • To suppress excited-state contamination, the authors utilize bi-local (Point-Split, PS) interpolating operators where nucleons are spatially separated, rather than local (hexa-quark) operators.
  • They implement a variational analysis using a $3 \times 3$correlation matrix constructed from operators with the three lowest momenta ($p=0, p=(2\pi/L)\hat{e}_x, p=(2\pi/L)(\hat{e}_x+\hat{e}_y)$).
  • Algorithmic Optimization: A specialized contraction algorithm is developed for bi-local sources. By treating the sink as a product of three-quark blocks and expanding the source fully, they utilize quark permutation symmetries to reduce computational redundancy (achieving $\sim 31\%$ efficiency for two-point functions).

• Weak Matrix Element Extraction:

  • The calculation uses the sequential-source propagator (background-field) technique combined with the summation method to enhance the signal-to-noise ratio.
  • Isospin symmetry is exploited via the Wigner-Eckart theorem to relate flavor-changing weak currents (required for $pp \to d$) to flavor-neutral currents, allowing the reuse of optimized contraction algorithms.
  • The matrix element is extracted from the ratio of three-point to two-point correlation functions, $R_3(t)$, using a two-state fit to isolate the ground state and subtract residual excited-state contamination.

• Finite-Volume Formalism ($2+J \to 2$):

  • The study implements the Lellouch-Lüscher formalism for $2+J \to 2$ processes to relate finite-volume lattice matrix elements to infinite-volume physical quantities.
  • This involves calculating the Lellouch-Lüscher factor, which depends on the derivative of the scattering phase shift. This factor accounts for the distortion of the two-nucleon wavefunction in the finite volume due to rescattering.
  • The infinite-volume amplitude is decomposed into one-body (single-nucleon) and two-body (contact interaction) contributions within Pionless Effective Field Theory ($\not\pi$EFT).

3. Key Contributions

• **First Application of $2+J \to 2$Formalism:** This is the first lattice QCD study to explicitly incorporate two-nucleon rescattering effects into the Lellouch-Lüscher factor for proton-proton fusion at$m_\pi \approx 432$ MeV.
• Operator Optimization: The work demonstrates that bi-local interpolators combined with variational analysis significantly suppress excited-state contamination compared to traditional local (hexa-quark) operators, which suffer from "pseudo-plateau" issues.
• Algorithmic Efficiency: The development of a tailored contraction algorithm for bi-local dibaryon operators reduces the computational cost of Wick contractions, making such calculations feasible.
• Systematic Uncertainty Analysis: The paper provides a rigorous cross-validation of energy shifts and matrix elements using ratio methods, two-state fits, and variational methods, highlighting the dangers of residual excited-state contamination in precision nuclear calculations.

4. Key Results

• Energy Spectrum & Scattering Parameters:
  • The ground states of the deuteron ($^3S_1$) and diproton ($^1S_0$) are overwhelmingly dominated by zero-momentum components ($>99\%$ overlap).
  • Extracted energy shifts: $\Delta_{^1S_0} = -19.4(7.9)$ MeV and $\Delta_{^3S_1} = -25.3(7.6)$ MeV.
  • Scattering parameters (via Effective Range Expansion) show large uncertainties but are consistent with a shallow bound state or attractive scattering state at this pion mass.
• Matrix Element Extraction:
  • The raw matrix element (normalized by the nucleon axial charge $g_A$) using bi-local operators is:
    $$\frac{|\langle d | J | pp \rangle|}{g_A} = 0.984(10)$$
  • The deviation from unity ($-0.016(10)$) indicates a small but non-vanishing contribution from two-body currents.
• Low-Energy Constant ($L_{1,A}$):
  • After applying the Lellouch-Lüscher finite-volume corrections (which can enhance the two-body contribution by a factor of $\sim 2$), the scale-independent constant is determined as $\tilde{L}_{1,A} = -0.11(12) \text{ fm}^{-1}$.
  • The renormalization-scale dependent constant at $\mu = m_\pi$ is:
    $$L_{1,A} = 6.0(7.1) \text{ fm}^3$$
  • Uncertainty Source: The large uncertainty is primarily driven by the imprecise determination of two-nucleon scattering parameters (scattering length and effective range) from the lattice, which propagate strongly into the finite-volume corrections.

5. Significance and Conclusion

• Validation of Feasibility: The work proves that ab initio lattice QCD calculations of weak two-nucleon reactions are feasible, even at unphysical pion masses, provided that advanced operator construction and finite-volume formalisms are used.
• Consistency with Phenomenology: Despite large statistical and systematic uncertainties, the extracted value of $L_{1,A}$ is consistent in magnitude and order with existing phenomenological extractions from experimental data (e.g., neutrino-deuteron scattering, tritium beta decay).
• Future Outlook: The study identifies the precision of two-nucleon scattering parameters as the bottleneck for reducing uncertainties in $L_{1,A}$. It suggests that future calculations at or near the physical pion mass would be advantageous, as experimental scattering lengths could be used as inputs, bypassing the need to infer them from noisy lattice energy shifts.
• Broader Impact: The methodology established here provides a robust framework for studying other weak two-nucleon processes, such as neutrinoless double-beta decay and muon capture, which are critical for testing fundamental physics beyond the Standard Model.