{"id":"0d94af9a-12a5-4e3c-b81a-2a1f9f3be9da","arxiv_id":"2411.19613","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A non-perturbative N3LO chiral lattice EFT calculation reproduces triton and helion energies and lifetime after fitting two three-nucleon contact constants to triton data.","lead":"This paper computes properties of the triton and helion (hydrogen-3 and helium-3) using a non-perturbative lattice version of chiral effective field theory at next-to-next-to-next-to-leading order. It fits two unknown constants to the triton binding energy and lifetime, then reports the helion energy and charge radii as independent checks.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Truncation-error estimate in Eq. (3) is the weakest link: setting X_N2LO = X_N3LO removes the largest-order term and could understate the error bars that support the agreement claim.","rationale":"The reader identified the truncation-error assumption as the weakest point; I agree. The most load-bearing aspect is not just that X_LO is replaced by X_NLO, but that X_N2LO is set equal to X_N3LO, which zeroes the Q-multiplied final term in Eq. (3). Because the calculation is not a complete N3LO calculation (the 3NF and axial current are only N2LO), the error estimate is more fragile than a standard N3LO truncation estimate. With the Table 1 shifts (NLO to N3LO: about 0.9 MeV in triton energy, about 0.1 fm in radii, about 0.6 yr in half-life), a genuine N2LO value halfway between would make the final term roughly Q times 0.45 MeV, about 0.15 MeV, which is larger than the currently largest truncation term. This would not destroy the central qualitative claim but would enlarge the error bars and weaken 'close to experiment.' A dedicated N2LO computation is the right test. The paper is transparent, the other independent checks are reasonable, and the smearing conclusion is unchanged, so the CONDITIONAL verdict remains appropriate.","tokens_in":7020,"tokens_out":10644,"duration_ms":89112,"concrete_test":"Compute a genuine N2LO result for the three-nucleon observables by using the NLO two-nucleon interaction (switching off the additional N3LO contact terms of Ref. [8]) together with the N2LO 3NF, refitting cE and cD to the same triton energy and half-life, and then inserting X_N2LO into Eq. (3). If the resulting truncation uncertainty for E_3He or R_ch,3H exceeds the quoted error by more than about 50%, the current error bars are understated and the agreement claim needs qualification. A cheaper complement is to recompute the error bars with the physical pion mass in Q (Q=139.6/600) and with Q=1/2 to test sensitivity to the arbitrary scale choice.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that the N3LO simulation gives a helion energy close to experiment and charge radii within a few percent is supported by error bars from Eq. (3). But that formula is only valid for a complete N^nLO calculation. Here the 3NF is kept only at leading (N2LO) order and the axial current only at N2LO; the calculation is not a complete N3LO chiral expansion. The paper sets X_LO = X_NLO and X_N2LO = X_N3LO because no N2LO (or LO) results are computed. The latter assignment is especially dangerous: in Eq. (3) the final term Q |X_N2LO - X_N3LO| (the one with the largest prefactor, Q=1/3) is set to zero. If an actual N2LO result lies between the NLO and N3LO values, this term alone can increase the truncation error by roughly 50% for the triton and helion energies, using the Table 1 NLO-to-N3LO shifts. Since the 'close to experiment' conclusion is drawn within these error bars, underestimating them could change the strength of the central claim. The paper explicitly warns that the estimate 'should be taken with care,' but no alternative systematic uncertainty is provided.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents non-perturbative exact-diagonalization calculations for the three-nucleon system using an N3LO chiral lattice interaction. The three-nucleon contact LECs (cE and cD, or C0,3N and cD) are adjusted to reproduce the triton ground-state energy and beta-decay half-life, with the axial current at N2LO; the helion ground-state energy and the charge radii of triton and helion are then computed as independent checks. The results are compared with experiment, with the NLO version of the same lattice setup, and with recent wave-function-matching lattice calculations. The central finding is that the N3LO calculation yields a helion energy consistent with experiment, charge radii a few percent below experiment, and no significant dependence on the smearing of the three-nucleon contact interaction.","tokens_in":7303,"tokens_out":11326,"duration_ms":92796,"significance":"If taken at face value, the results provide a useful benchmark for non-perturbative chiral lattice EFT in the A=3 sector without a Monte Carlo sign problem. The helion energy and charge radii serve as genuine independent checks because the triton energy and half-life are fit targets, and the paper is transparent about which quantities are fitted. The exact-diagonalization approach and the explicit LEC determination make the calculation reproducible in principle. The quantitative significance is limited by the incompleteness of the N3LO truncation-error estimate and by the single lattice size used for the radii, but the qualitative demonstration of a non-perturbative N3LO three-nucleon simulation is valuable for the lattice-EFT community.","major_comments":[{"comment":"The truncation uncertainty used to set the error bars in Table 1 is not a complete N3LO estimate. With only NLO and N3LO results, setting X_LO = X_NLO and X_N2LO = X_N3LO removes the term with the largest prefactor in Eq. (3), namely Q|X_N2LO - X_N3LO| with Q = 1/3, and the retained 'NLO' calculation still uses the N3LO two-nucleon interaction from Ref. [8] with only the three-nucleon force and the axial current omitted, so X_NLO is not a pure chiral-order result. The tabulated errors should therefore be understood as lower bounds on the truncation error, yet they are used for the comparison with experiment and with Refs. [3,5]. Please provide a conservative alternative (for example, retaining Q|X_NLO - X_N3LO| or performing a sensitivity study with varied Q) or state explicitly in the text and in Table 1 that the quoted errors are not full truncation uncertainties.","section":"Section 3, Eq. (3)"},{"comment":"The fitting criterion of bisecting the LEC interval until the observables are closer to the experimental values than the truncation error bar makes the LEC uncertainties in Eqs. (1)-(2) conditional on the truncation estimate discussed above. Because the acceptance interval is the truncation error, an underestimated truncation error directly shrinks the fitted LEC ranges and the propagated fitting uncertainties. Please propagate a more conservative truncation estimate through the LEC fit, or report the LEC uncertainties separately and state their dependence on the adopted acceptance criterion.","section":"Section 3, LEC fitting"}],"minor_comments":[{"comment":"The statement 'we thus set X_LO = X_NLO and X_N2LO = X_N3LO' should be spelled out by explicitly writing the retained terms in Eq. (3), since a reader may otherwise think the NLO-N2LO difference is also discarded.","section":"Section 3, before Eq. (3)"},{"comment":"The charge radii are computed at a single lattice size L = 9a without the infinite-volume extrapolation used for the energies; a brief estimate of the expected finite-size effect would help assess the few-percent deficit relative to experiment.","section":"Section 3.2"},{"comment":"Marking the fitted rows (triton ground-state energy and triton half-life) with a footnote symbol would make the by-construction agreement visible at a glance.","section":"Table 1"},{"comment":"The statement that the locality effect observed in Ref. [11] likely only occurs in systems with more than three nucleons is speculative; consider softening the wording.","section":"Section 4"},{"comment":"The abstract contains a spacing typo in 'byadjusting' that should be corrected.","section":"Abstract"}],"recommendation":"major_revision","confidential_remarks":"This is a proceedings contribution and the technical calculation is competently executed. The main reason for major revision is that the central quantitative comparison rests on a truncation-error estimate that the authors themselves flag as fragile, and the current error bars in Table 1 are best interpreted as lower bounds. If the authors instead reposition the paper as an explicitly exploratory study with a conservative error budget or a clear statement that the quoted uncertainties are not full truncation uncertainties, a minor revision could be sufficient."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This paper does something genuinely new in a modest way: it shows that a non-perturbative N3LO chiral lattice EFT calculation of the three-nucleon bound states is feasible without wavefunction matching, and it extracts values for the leading 3NF LECs. The helion ground-state energy comes out close to experiment, and the charge radii are a few percent low, with the gap plausibly blamed on omitted exchange currents. The two-nucleon interaction and axial current are inherited from previous work, so the new content is the A=3 computation and the LEC fit -- a legitimate extension rather than a discovery, but a useful one for the subfield.\n\nThe strengths are real: the exact diagonalization and fitting are standard and described cleanly; Table 1 is honest about which quantities are fitted (triton energy and half-life) and which are predictions (helion energy, both charge radii); and the comparison to the wavefunction-matching results is fair. The null result for smearing of the 3N contact is stated without overinterpretation. This is careful work.\n\nThe soft spot is the truncation error in Eq. (3). The formula is designed for a complete N^nLO expansion, but only NLO and N3LO points exist, so the authors set X_LO = X_NLO and X_N2LO = X_N3LO. The second substitution zeroes the term with the largest prefactor, Q|X_N2LO - X_N3LO|, and that term could add roughly 50% to the energy error bars if a real N2LO point lies between NLO and N3LO. The paper explicitly warns that the estimate \"should be taken with care\", and it does -- but the quoted LEC uncertainties and the \"close to experiment\" language lean on those error bars. Since the helion energy is a real prediction and it agrees even under a more conservative error budget, the central claim survives; the problem is more about overstated precision on the LECs. The fitted triton observables are by construction, so they should not be counted as support.\n\nWho is this for? Chiral lattice EFT practitioners and few-nucleon theorists who want a non-perturbative N3LO benchmark or starting LEC values. It deserves refereeing: the computation is clean, the limitations are stated rather than hidden, and the truncation-error issue is addressable in review, not fatal. My advice: send it to review, and ask the authors to add a conservative alternative truncation estimate or an explicit N2LO point, but publish it as an exploratory benchmark.","headline":"A transparent, technically sound A=3 lattice EFT benchmark whose truncation-error estimate is the soft spot -- and the authors already flag it.","tokens_in":7831,"tokens_out":2227,"would_cite":true,"duration_ms":20232,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper demonstrates that a non-perturbative N3LO chiral lattice calculation, with three-nucleon contact couplings fitted to triton binding and half-life, reproduces the helion energy and predicts charge radii a few percent below…","keywords":["chiral effective field theory","nuclear lattice EFT","three-nucleon force","triton half-life","helion charge radius","N3LO","low-energy constants","non-perturbative calculation"],"falsifier":"Compute the same three-nucleon observables at full N2LO order and check whether the N2LO point falls inside the error band estimated from Eq. (3); if it falls outside, the truncation uncertainty is underestimated and the fitted low-energy constants are less constrained than claimed.","tokens_in":6807,"feed_emoji":"⚛️","tokens_out":8236,"duration_ms":59493,"temperature":0.7,"pith_summary":"This paper aims to show that three-nucleon observables can be computed non-perturbatively in chiral effective field theory on a lattice at next-to-next-to-next-to-leading order (N3LO), without expanding the subleading nuclear forces in perturbation theory. After fitting the two three-nucleon contact low-energy constants to the triton ground-state energy and beta-decay half-life, the calculation reproduces the helion ground-state energy within uncertainties and predicts triton and helion charge radii a few percent below the experimental values. The paper also finds that smearing the three-nucleon contact interaction has no significant effect in the three-nucleon system, and that including data from heavier nuclei in the fit does not clearly improve these observables. If the result holds, it demonstrates a reliable non-perturbative route through the higher-order chiral expansion for light nuclei.","feed_headline":"Lattice EFT matches triton and helion without perturbative forces","feed_subtitle":"N3LO chiral calculation reproduces the helion energy and puts charge radii a few percent low.","key_machinery":"The central object is the N3LO lattice Hamiltonian for the three-nucleon system: the previously built N3LO two-nucleon interaction of Ref. [8] together with the leading three-nucleon force $V_{cE}^{(0)}$, $V_{cD}^{(0)}$, and $V_{3N}^{\\rm TPE}$. Because the system is small, the Hamiltonian is diagonalized non-perturbatively: the ground-state energy is obtained by a Lanczos solve on lattices of lengths $L=4a,\\dots,9a$ and then extrapolated to infinite volume with the three-body formula $E(L\\to\\infty)=E_\\infty+E_0 L^{-3/2}\\exp(-L/L_0)$ from Ref. [22]. The carrying objects for the fit are the two low-energy constants $c_D$ and $c_E$ (or $C_{0,3N}$), tuned to the triton half-life and ground-state energy respectively, with the truncation uncertainty of Eq. (3) estimated by setting $X_{\\rm LO}=X_{\\rm NLO}$ and $X_{\\rm N2LO}=X_{\\rm N3LO}$ because only NLO and N3LO results were computed.","core_discovery":"Working on a lattice with spacing $a=1.9733$ fm, the authors solve the few-nucleon Schrödinger equation exactly at N3LO instead of treating subleading forces perturbatively. They take the N3LO two-nucleon lattice interaction from Ref. [8], add the leading three-nucleon force terms $V_{cE}^{(0)}$, $V_{cD}^{(0)}$, and $V_{3N}^{\\rm TPE}$, fix $c_D$ by the triton half-life computed with the N2LO axial current, and fix either $c_E$ or the smeared contact coefficient $C_{0,3N}$ by the triton ground-state energy. The resulting helion energy lands within the experimental error bars, while both charge radii come out a few percent low, which the authors attribute to neglected exchange contributions to the charge density operator. The central quantitative result is that the NLO helion deviation is cured by the N2LO three-nucleon force, and that the choice between unsmeared $c_E$ and smeared $C_{0,3N}$ does not matter at the three-nucleon level.","pith_inferences":["The authors do not test this, but if the charge-radius shortfall is really caused by omitted exchange currents, adding those two-body charge-density operators should shift both radii upward by roughly two to four percent while leaving the ground-state energies nearly unchanged.","By the same logic, running the same Hamiltonian at $A=4$ or $A=6$ would separate the smearing-insensitivity seen here from the multi-nucleon correlations that make the local smeared contact matter in medium-mass binding.","A direct N2LO calculation would provide the missing anchor for the truncation estimate in Eq. (3); if the N2LO point lies outside the band built from NLO and N3LO values, the quoted LEC errors are understated.","Combining this non-perturbative N3LO Hamiltonian with wavefunction matching could push the same fitted interactions to $A\\geq 4$ without a sign problem, giving an independent check on the fitted LECs."],"forward_implications":["The N2LO three-nucleon force removes the NLO discrepancy in the helion ground-state energy, so the triton-helion splitting emerges without fitting to any $A>3$ nucleus.","Triton and helion charge radii are predicted a few percent below experiment, pointing to two-body charge-density (exchange-current) operators as the likely missing input in Eq. (5).","Smearing the three-nucleon contact interaction has no significant effect at $A=3$, so any locality effect seen in the literature likely requires systems with more than three nucleons.","Including bound-state energies of heavier nuclei in the LEC fit, as in the wavefunction-matching approach, does not noticeably improve the three-nucleon charge radii, and a leading-order axial current with neglected truncation uncertainty appears insufficient for the triton half-life.","The constructed N3LO interaction can be applied to nucleon-deuteron scattering through the adiabatic projection method."],"supporting_citations":[{"why":"Supplies the N3LO two-nucleon lattice interaction (24 contact terms, pion exchange, Coulomb) used as the two-body input.","marker":"[8]"},{"why":"Provides the leading three-nucleon force terms $V_{cE}^{(0)}$, $V_{cD}^{(0)}$, and $V_{3N}^{\\rm TPE}$ and the main wavefunction-matching comparison.","marker":"[3]"},{"why":"Gives the N2LO axial current used to compute the Gamow-Teller matrix element in the triton half-life.","marker":"[7]"},{"why":"Sets the pion-nucleon LECs $c_1$, $c_3$, $c_4$ appearing in the two-pion-exchange three-nucleon force.","marker":"[12]"},{"why":"Non-perturbative lattice-EFT triton lifetime calculation used as a comparison for the half-life result.","marker":"[5]"},{"why":"Derives the finite-volume three-body energy extrapolation formula used to take $L$ to infinity.","marker":"[22]"},{"why":"Provides the experimental triton half-life used as the fit target for $c_D$.","marker":"[20]"}],"fun_headline_variants":["Non-perturbative N3LO lattice EFT hits triton and helion","Exact lattice N3LO three-nucleon solution without perturbation","Chiral lattice EFT at N3LO reproduces triton and helion","No-perturbation lattice EFT matches triton and helion energies","Full N3LO lattice chiral EFT nails triton and helion ground states"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The error bars rest on treating the NLO result as a proxy for the missing leading-order result and the N3LO result as a proxy for the missing N2LO result, since only NLO and N3LO computations were actually performed.","fun_headline_variants_meta":{"raw":{"variants":["Non-perturbative N3LO lattice EFT hits triton and helion","Exact lattice N3LO three-nucleon solution without perturbation","Chiral lattice EFT at N3LO reproduces triton and helion","No-perturbation lattice EFT matches triton and helion energies","Full N3LO lattice chiral EFT nails triton and helion ground states"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000302,"raw_usage":{"total_tokens":1756,"prompt_tokens":978,"completion_tokens":778,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":594,"completion_tokens_details":{"reasoning_tokens":675}},"tokens_in":594,"tokens_out":778,"duration_ms":5927,"temperature":1.0,"reasoning_tokens":675,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T05:59:57.634531+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the same three-nucleon observables at full N2LO order and check whether the N2LO point falls inside the error band estimated from Eq. (3); if it falls outside, the truncation uncertainty is underestimated and the fitted low-energy constants are less constrained than claimed.","supporting_citations":[{"cited_title":"Spectrum of three-body bound states in a finite volume","cited_arxiv_id":"1412.4969","evidence_quote":"Derives the finite-volume three-body energy extrapolation formula used to take $L$ to infinity."},{"cited_title":"Simpson,Half-life of tritium and the Gamow-Teller transition rate, Phys","cited_arxiv_id":null,"evidence_quote":"Provides the experimental triton half-life used as the fit target for $c_D$."}],"review_version":1}