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

Bo-Hao Jian; Chuan Liu; Lu-Chang Jin; Xu Feng; Zi-Yu Wang

REVIEW 3 major objections 3 minor 1 cited by

Lattice QCD determines the proton-proton fusion matrix element at pion mass 432 MeV, finding a small two-body current after finite-volume rescattering corrections.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · grok-4.5

2026-07-15 00:08 UTC pith:4LIZGHLB

load-bearing objection First lattice QCD shot at the pp-fusion matrix element with LL rescattering; the raw ME looks clean, but L_{1,A} is still only a naturalness check. the 3 major comments →

arxiv 2603.09554 v3 pith:4LIZGHLB submitted 2026-03-10 hep-lat

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

Zi-Yu Wang , Xu Feng , Bo-Hao Jian , Lu-Chang Jin , Chuan Liu This is my paper

classification hep-lat

keywords lattice QCDproton-proton fusiontwo-body currentsLellouch-Lüscheraxial matrix elementtwo-nucleon systemsstellar nucleosynthesisfinite-volume corrections

verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper reports a first-principles lattice QCD calculation of the matrix element that governs proton-proton fusion, the weak process that begins stellar energy production, performed at a pion mass of roughly 432 MeV. The authors apply Lellouch-Lüscher finite-volume corrections inside a two-nucleon-plus-current framework so that the lattice matrix elements can be converted into infinite-volume quantities that include two-nucleon rescattering. Before those corrections the ratio of the fusion matrix element to the nucleon axial charge is 0.984(10); the small departure from unity signals a non-vanishing two-body current. After the corrections they extract the two-body low-energy constant L_{1,A}=6.0(7.1) fm^{3}, a value still compatible with natural-size phenomenological estimates even though the large uncertainties on the two-nucleon scattering parameters propagate strongly into the final error. The calculation therefore demonstrates both the feasibility of ab-initio lattice studies of weak two-nucleon reactions and the practical obstacles that must be overcome before results at the physical pion mass become competitive.

Core claim

At m_π≈432 MeV the lattice matrix element for proton-proton fusion satisfies ⟨d|J|pp⟩/g_A=0.984(10) before Lellouch-Lüscher rescattering corrections, indicating a small but non-zero two-body current; after the corrections the two-body low-energy constant is L_{1,A}=6.0(7.1) fm^{3}, consistent at the level of naturalness with phenomenological extractions.

What carries the argument

Lellouch-Lüscher finite-volume factors constructed inside a 2+J→2 framework convert the lattice matrix elements into infinite-volume amplitudes while incorporating two-nucleon rescattering; the factors are built from scattering parameters obtained via Lüscher’s formalism on the same ensembles, using bi-local nucleon-nucleon interpolators and a variational analysis over the three lowest momenta to suppress excited-state contamination.

Load-bearing premise

The two-nucleon scattering parameters extracted on the same lattices are accurate enough that their large uncertainties, once inserted into the Lellouch-Lüscher factors, still permit a meaningful determination of the two-body low-energy constant.

What would settle it

A higher-statistics determination of the two-nucleon scattering length and effective range on the same ensembles that shifts L_{1,A} by more than its present 7.1 fm^{3} uncertainty, or a recalculation of the same matrix element that finds ⟨d|J|pp⟩/g_A consistent with exactly 1 after identical finite-volume corrections.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

The same 2+J→2 framework can be repeated at lighter pion masses to remove the dominant systematic associated with the unphysical quark mass.
The measured 1.6 percent deviation from unity quantifies the size of two-body axial currents that must be included in few-body weak processes.
Any future improvement in lattice two-nucleon scattering lengths and effective ranges will translate almost linearly into a tighter error on L_{1,A}.
The calculation supplies a concrete template for other weak two-nucleon reactions that enter stellar modeling or neutrino physics.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

Because scattering-parameter errors dominate the final uncertainty, simultaneous multi-ensemble fits of spectra and matrix elements will be required before physical-point results can compete with phenomenology.
The identical finite-volume machinery can be applied without major redesign to related processes such as neutrino-deuteron breakup or tritium β-decay, where the same two-body axial current appears.
At the physical pion mass the deuteron becomes far more loosely bound, so the finite-volume corrections themselves will change character; the present heavy-mass study mainly validates the numerical strategy rather than the final number.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Desk Editor's Note

First lattice QCD shot at the pp-fusion matrix element with LL rescattering; the raw ME looks clean, but L_{1,A} is still only a naturalness check.

read the letter

The one thing to know is that this is the first lattice attempt to extract the pp-fusion matrix element with Lellouch-Lüscher factors inside a genuine 2+J→2 setup, using bi-local NN operators and a multi-momentum variational basis. At m_π≈432 MeV they report ⟨d|J|pp⟩/g_A=0.984(10) before the LL corrections—small, non-zero two-body current—and after the factors L_{1,A}=6.0(7.1) fm³.

What is actually new is the combination: they measure the finite-volume matrix element, extract the two-nucleon spectrum and scattering parameters on the same ensembles via Lüscher, then fold those into the LL conversion. The pre-correction number is clean and the deviation from unity is stated with a few-percent error, which is useful. They are also frank that the large uncertainties on the 2N scattering lengths and effective ranges “propagate strongly” into the FV factors, so the final LEC error bar exceeds 100 %. That honesty is a plus; they do not oversell the precision.

The soft spot is exactly the one the stress-test flags and the abstract itself admits: once the LL factors are applied, the two-body LEC is only constrained at the level of naturalness. It is compatible with phenomenology, but it does not yet give a quantitatively useful first-principles value. Everything is still at a heavy pion mass, and without the full paper we cannot audit the variational spectrum, the LL derivation, or the error budget. Those are real limitations, not fatal ones for a methods paper.

This is for people who work on nuclear lattice QCD or few-body weak processes. Anyone who needs a controlled number for solar-neutrino or stellar-energy calculations will wait for physical-mass follow-ups, but the technical path is now laid out. I would send it to a serious referee; the calculation is ab initio, the caveats are stated, and the field needs these foundations even when the error bars are still large. Accept for review, expect the referees to push hard on the scattering-parameter systematics and the continuum/physical-mass roadmap.

Referee Report

3 major / 3 minor

Summary. The manuscript reports a lattice QCD calculation of the proton–proton fusion matrix element at m_π ≈ 432 MeV. Bi-local nucleon–nucleon interpolators and a variational analysis with the three lowest momenta are used to control excited-state contamination and to extract the two-nucleon spectrum and scattering parameters via Lüscher’s formalism. Lellouch–Lüscher finite-volume corrections are implemented in a 2+J→2 framework that includes two-nucleon rescattering. Before rescattering corrections the authors obtain ⟨d|J|pp⟩/g_A = 0.984(10); after applying the LL factors they quote L_{1,A} = 6.0(7.1) fm³, which is compatible with phenomenological extractions only at the level of naturalness. The abstract itself emphasizes that large uncertainties in the two-nucleon scattering parameters propagate strongly into the finite-volume corrections, rendering a precise determination of the two-body LEC highly challenging with present lattice inputs.

Significance. A controlled, first-principles extraction of the weak two-nucleon matrix element that initiates stellar hydrogen burning would be a landmark result for nuclear lattice QCD and for solar-model inputs. The work correctly identifies the methodological path (variational spectrum + Lüscher + 2+J→2 Lellouch–Lüscher) and honestly reports that current two-nucleon scattering uncertainties dominate the final error budget. Even a naturalness-level constraint on L_{1,A} at an unphysical pion mass, together with a demonstrated pipeline, constitutes a useful foundation for future calculations nearer the physical point. The explicit acknowledgment of the dominant systematic is itself a strength of the presentation.

major comments (3)

Abstract (result for L_{1,A}): The quoted value L_{1,A}=6.0(7.1) fm³ carries a relative uncertainty larger than 100 %, driven by the propagation of two-nucleon scattering-parameter errors through the Lellouch–Lüscher factors. While the authors state that a precise determination remains highly challenging, the title and the framing as a “first-principles determination” of the fusion matrix element risk overstating the quantitative reach. The manuscript should clarify, already in the abstract and introduction, that the present result constrains only the natural size of the two-body current and does not yet furnish a precision extraction of L_{1,A}.
Abstract (LL conversion): The central claim rests on converting the finite-volume matrix element ⟨d|J|pp⟩/g_A=0.984(10) into the infinite-volume LEC via 2+J→2 Lellouch–Lüscher factors built from the same-ensemble Lüscher scattering parameters. Because only the abstract is available, the derivation of those factors, the variational spectrum, the continuum/volume extrapolations (if any), and the full error budget cannot be inspected. A load-bearing assessment of whether residual systematics are under control relative to the quoted uncertainty is therefore impossible; the full manuscript must supply these details before the conversion can be regarded as validated.
Abstract (pion-mass caveat): All results are obtained at a single, heavy pion mass m_π∼432 MeV. The abstract does not indicate any chiral extrapolation or estimate of the pion-mass dependence of either the one-body-normalized matrix element or L_{1,A}. Given that the two-body current is expected to be sensitive to the long-range pion-exchange contribution, the manuscript should quantify (or at least bound) the expected shift toward the physical point; otherwise the phenomenological comparison remains only qualitative.

minor comments (3)

Abstract: The notation “2+J→2 framework” and the precise definition of the Lellouch–Lüscher factor that multiplies the two-body contribution should be introduced more explicitly for non-specialist readers of the abstract.
Abstract: The phrase “compatible, at the level of naturalness, with phenomenological extractions” would be strengthened by quoting the numerical range of the phenomenological L_{1,A} values being compared against.
Abstract: A brief statement of the lattice spacing(s), volume(s), and number of gauge configurations used would help the reader gauge the statistical and discretization context of the quoted 0.984(10) result.

Circularity Check

0 steps flagged

No significant circularity: ab initio lattice matrix element converted via standard LL/Lüscher formalism; large 2N uncertainties are a precision limitation, not a definitional loop.

full rationale

Only the abstract is available. The claimed result is a direct lattice QCD measurement of the finite-volume pp-fusion matrix element at m_π ≃ 432 MeV, followed by conversion to an infinite-volume two-body LEC via the standard Lellouch-Lüscher (2+J o2) and Lüscher formalisms. The abstract reports ⟥d|J|pp⟩/g_A = 0.984(10) before rescattering corrections and L_{1,A} = 6.0(7.1) fm^{3} after them. The same ensembles supply both the matrix element and the 2N scattering parameters that enter the LL factors; this is ordinary lattice practice, not a self-definitional or fitted-input circularity. The large propagated uncertainty is explicitly acknowledged and does not force the central number by construction. No uniqueness theorem, ansatz smuggled via self-citation, or renaming of a known empirical pattern is invoked. Self-citation risk is low on the abstract. The derivation is therefore self-contained against external benchmarks at the level of the available text; score 0 is the honest finding. (Precision/naturalness caveats belong to correctness risk, not circularity.)

Axiom & Free-Parameter Ledger

3 free parameters · 3 axioms · 0 invented entities

Abstract-only review: free parameters and axioms are those standard to unphysical-mass lattice QCD plus the LL mapping assumptions. No new particles or forces are invented; L_{1,A} is an EFT LEC being determined, not postulated. Dominant free inputs are the lattice pion mass, volume, and the 2N scattering lengths/effective ranges that enter the LL factor.

free parameters (3)

m_π (lattice pion mass) = ~432 MeV
Calculation performed at m_π~432 MeV, far above physical; all matrix elements and scattering parameters are at this mass and must later be extrapolated.
2N scattering parameters (a, r, ...) = not quoted in abstract
Extracted via Lüscher formalism on the same ensembles and then fed into LL factors; their large uncertainties dominate the final L_{1,A} error.
lattice spacing / volume set
Standard lattice discretization and finite-volume parameters that set the scale and the FV spectrum; not given numerically in the abstract.

axioms (3)

domain assumption Lellouch-Lüscher finite-volume formalism correctly maps the 2+J→2 matrix element including two-nucleon rescattering at the simulated mass and volumes.
Central conversion step from FV lattice data to infinite-volume matrix element and L_{1,A}; validity assumed for the kinematics used.
domain assumption Bi-local NN interpolators plus a three-momentum variational basis sufficiently suppress excited-state contamination for the ground-state matrix element.
Stated as the method that enables clean spectrum and matrix-element extraction; residual contamination is a systematic.
standard math Standard continuum QCD and electroweak axial current at the quark level, with lattice discretization artifacts under control or negligible for the quoted precision.
Background of any lattice QCD weak matrix element calculation.

pith-pipeline@v1.1.0-grok45 · 6227 in / 2812 out tokens · 30724 ms · 2026-07-15T00:08:10.168284+00:00 · methodology

0 comments

read the original abstract

Proton-proton fusion is the fundamental weak reaction initiating stellar energy production, and a first-principles determination of its matrix element remains a long-standing goal of nuclear theory. We present a lattice QCD calculation of the pp fusion matrix element at m_pi~432 MeV. We implement Lellouch-Luscher (LL) finite-volume (FV) corrections within a 2+J->2 framework, accounting for two-nucleon (2N) rescattering, to relate FV matrix elements to infinite-volume counterparts. Excited-state contamination is suppressed using bi-local nucleon-nucleon interpolating operators and a variational analysis with three lowest momenta. This enables determination of 2N energy spectrum and scattering parameters via Luscher's FV formalism. Before including rescattering effects in the LL factor, we obtain <d|J|pp>/g_A = 0.984(10), where g_A is the axial charge. The deviation from unity indicates a small nonvanishing 2-body current contribution. Our analysis shows that rescattering effects in LL factors substantially modify the 2-body contribution, while large uncertainties in 2N scattering parameters propagate strongly into FV corrections. Thus, precise determination of the 2-body low-energy constant L_{1,A} remains highly challenging with current lattice inputs. Despite the large uncertainty, L_{1,A}=6.0(7.1) fm^3 is compatible, at the level of naturalness, with phenomenological extractions. This work demonstrates feasibility and intrinsic challenges of ab initio lattice QCD calculations of weak 2N reactions, and establishes a foundation for future studies at or near the physical pion mass.

discussion (0)

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

Analytic decomposition of two-body electroweak processes with left-hand cuts
hep-lat 2026-07 accept novelty 5.0

An on-shell decomposition of 2+J o2 electroweak amplitudes isolates OPE poles, logs, and triangle singularities, leaving only smooth short-distance functions.