{"id":"3b9bfc8d-5bdb-468b-a6a1-21529c8191a0","arxiv_id":"2603.09554","paper_version":3,"verdict":"CONDITIONAL","confidence":"LOW","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"At m_π≈432 MeV, lattice QCD gives ⟨d|J|pp⟩/g_A=0.984(10) before rescattering LL corrections and L_{1,A}=6.0(7.1) fm³ after, with large FV-propagated errors.","lead":"Lattice QCD is used to compute the proton-proton fusion matrix element at a pion mass of about 432 MeV, including finite-volume corrections for two-nucleon rescattering. The result is a first-principles handle on a key stellar weak reaction and on the two-body axial low-energy constant, with large remaining uncertainties.","discovery_kind":"first_principles","skeptic_critique":{"model":"grok-4.5","headline":"Large 2N scattering uncertainties dominate LL factors, so L_{1,A}=6.0(7.1) fm^{3} only constrains naturalness, not a precise two-body current.","rationale":"The reader correctly isolates the dominant soft spot: the propagation of 2N scattering-parameter errors through the LL factors. The abstract already flags this limitation and presents L_{1,A} only as a naturalness-level check, so the central claim as worded survives, but remains conditional on improved scattering inputs (and ultimately physical-mass ensembles). No stronger internal inconsistency appears in the available material; the CONDITIONAL/LOW-confidence verdict is therefore left unchanged.","tokens_in":2138,"tokens_out":447,"duration_ms":21181,"concrete_test":"From the full paper extract the numerical 2N scattering lengths/effective ranges (with covariances) used for the LL factors; recompute the LL correction and the resulting L_{1,A} by Monte-Carlo sampling of those parameters within their quoted errors. If the reconstructed central value and uncertainty match 6.0(7.1) and the error remains >100 % relative, the concern is confirmed.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim rests on converting the reported FV matrix element into the infinite-volume two-body LEC L_{1,A} via Lellouch-Lüscher factors in the 2+J→2 framework. Those factors are built from the same-ensemble 2N scattering parameters obtained with Lüscher’s method. The abstract itself states that the large uncertainties on those parameters “propagate strongly into FV corrections,” producing the final result L_{1,A}=6.0(7.1) fm^{3} whose relative error exceeds 100 %. Consequently the number is compatible with phenomenological extractions only at the level of naturalness; a quantitatively useful first-principles determination of the two-body current is not yet achieved. This is the least secure link between the lattice matrix element and the claimed LEC.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.5","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.","tokens_in":2358,"tokens_out":1099,"duration_ms":17869,"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":[{"comment":"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}.","section":null},{"comment":"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.","section":null},{"comment":"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.","section":null}],"minor_comments":[{"comment":"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.","section":null},{"comment":"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.","section":null},{"comment":"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.","section":null}],"recommendation":"uncertain","confidential_remarks":"Only the abstract was available for this review; a definitive recommendation (accept / minor / major revision) cannot be issued until the full manuscript, including the derivation of the LL factors, the variational spectrum, the scattering-parameter extraction, and the complete error budget, can be examined. On the basis of the abstract alone the scientific narrative is careful and the dominant limitation is openly stated, but the quantitative claim in the title appears stronger than the reported uncertainty on L_{1,A} supports. I recommend that the editor obtain the full text before proceeding."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","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³.\n\nWhat 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.\n\nThe 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.\n\nThis 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.","headline":"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.","tokens_in":2985,"tokens_out":559,"would_cite":false,"duration_ms":8406,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"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.","keywords":["lattice QCD","proton-proton fusion","two-body currents","Lellouch-Lüscher","axial matrix element","two-nucleon systems","stellar nucleosynthesis","finite-volume corrections"],"falsifier":"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.","tokens_in":3036,"feed_emoji":"⚛️","tokens_out":1043,"duration_ms":18545,"temperature":0.7,"pith_summary":"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.","feed_headline":"Lattice QCD finds small two-body current in pp fusion","feed_subtitle":"At heavy pion mass the matrix element is 0.984 before rescattering; after corrections L1A matches natural size","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"fun_headline_variants":["Lattice QCD finds 0.984 pp-fusion matrix element at heavy pion mass","Small two-body current in lattice QCD proton-proton fusion","pp fusion matrix element near unity before LL rescattering","Lattice QCD extracts L1A of 6.0 fm3 with large uncertainty","Ab initio lattice QCD probes weak two-nucleon current in fusion"],"cache_read_input_tokens":128,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Lattice QCD finds 0.984 pp-fusion matrix element at heavy pion mass","Small two-body current in lattice QCD proton-proton fusion","pp fusion matrix element near unity before LL rescattering","Lattice QCD extracts L1A of 6.0 fm3 with large uncertainty","Ab initio lattice QCD probes weak two-nucleon current in fusion"]},"model":"grok-4.5","effort":"low","cost_usd":0.003924,"raw_usage":{"total_tokens":1318,"prompt_tokens":900,"num_sources_used":0,"completion_tokens":95,"cost_in_usd_ticks":39240000,"prompt_tokens_details":{"text_tokens":900,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":323,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":900,"tokens_out":95,"duration_ms":3332,"temperature":1.0,"reasoning_tokens":323,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-15T00:08:10.168284+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"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.","supporting_citations":[],"review_version":1}