{"id":"1d8c24b8-f731-48d1-ae84-27466ff9b8ce","arxiv_id":"2607.26958","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"A spectroscopic reparametrization of N4LO contact three-nucleon forces lets nine of thirteen LECs be determined from elastic nucleon-deuteron scattering via an RBF emulator.","lead":"A new spectroscopic basis reorganizes short-range three-nucleon contact forces by angular momentum and parity, so each constant maps to specific nucleon-deuteron partial waves. With that basis the authors build a cheap emulator and show nine of thirteen low-energy constants can be fixed from elastic Nd data.","discovery_kind":"new_method","skeptic_critique":{"model":"grok-4.5","headline":"Progressive fits already show order-one/sign-flipping LEC shifts that exceed quoted statistical errors, so “reliably determined” is not secured even inside the incomplete model.","rationale":"The reader correctly flags incompleteness of the 3NF (no N3LO+ long/intermediate range, no Coulomb, angle cut) as the main reason the fit is only exploratory and scores CONDITIONAL. That concern is real and author-acknowledged (Secs. III A, III C, IV). The more immediate load-bearing soft spot for the strongest claim as written—“9 of 13 … reliably determined”—is already visible inside the published progressive fits and χ²: LEC values move by more than their quoted errors and sometimes flip sign when channels are added, s6 is unnatural, and χ²/datum=4.44 means the statistical errors understate uncertainty. This is the same reliability issue the reader ties to future completion of the Hamiltonian, but it does not require waiting for N3LO pieces; it is testable on the existing Table V setup via a naturalness-prior re-fit. I therefore agree on CONDITIONAL and on the spirit of the weakest assumption, but locate the sharpest present failure mode in the internal progressive-fit instability rather than only in the prospective N3LO restoration. No change of verdict is warranted: the spectroscopic basis, PWD appendices, emulator QA (App. D), and held-out predictions remain a solid methods contribution; the over-claim is confined to the strength of the word “reliably” and the precise count nine.","tokens_in":38497,"tokens_out":835,"duration_ms":47169,"concrete_test":"Re-minimize the full-fit χ² of Sec. III C / Table V last column under a hard naturalness prior |si|≤2 (and optionally rescale the covariance by √(χ²/datum)≈2.1). Report which of s1,s3,s5–s11 remain interior and nonzero at >2σ and the new χ²/datum. If s6 (and possibly s9) hit the boundary and/or several others lose significance, the “9 of 13 reliably determined” claim does not hold even inside the present Hamiltonian.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim (abstract; Sec. III C; Table V) is that 9 of 13 spectroscopic LECs are reliably fixed by elastic Nd data. Within the authors’ own incomplete Hamiltonian, that claim rests on the full-fit column of Table V and on the statement that LECs are stable and only weakly correlated once s2=s4=0. But the progressive columns of the same table already show shifts larger than the quoted 1σ errors, including sign flips: s6 goes from −6.461±0.487 (5/2−+3/2−+1/2−) to +6.294±0.472 (full fit); s5 flips −0.505→+1.464; s10 flips −1.343→+0.823 when the 1/2− sector is opened. The authors themselves note (Sec. III C) that changes from simplest to fullest fit typically exceed statistical errors and point to unquantified systematics, and that χ²/datum=4.44 (vs 19.26 with contacts off). With χ²/N≃4.4 the covariance used for “statistical errors” is mis-scaled if the model is incomplete, and s6≃6 lies well outside the natural window |si|≲2 they adopt from Table I and NN experience. Thus “reliably determined” is not established by the internal evidence even before N3LO long-range 3NFs are restored; the spectroscopic basis and emulator remain useful, but the count “9 of 13” and the reliability language over-reach what Table V and the χ² support.","agreement_with_reader":"partial"},"referee_report":{"model":"grok-4.5","summary":"The paper introduces a spectroscopic basis for the 13 subleading contact three-nucleon operators of chiral EFT at N4LO, reparametrizing the standard Ei LECs into Si combinations that each act in a single JP(T) channel (or a single Nd partial-wave/mixing sector). Explicit partial-wave decompositions, triton expectation-value benchmarks, and the inverse transformation are supplied. Using this basis the authors map the sensitivity of elastic Nd analyzing powers, build a low-dimensional RBF emulator of the Nd transfer matrix, and perform progressive exploratory fits of cD, cE and the isospin-1/2 spectroscopic LECs to selected pd data at 10, 70 and 135 MeV. They conclude that nine of the thirteen LECs (cD, cE, s1, s3, s5–s11) can be constrained by elastic Nd scattering, while s2 and s4 are largely redundant for elastic observables and s12, s13 live only in T=3/2.","tokens_in":38946,"tokens_out":1569,"duration_ms":38903,"significance":"A transparent spectroscopic organization of the N4LO contact 3NF is a genuine and reusable contribution: it clarifies which LECs control which Nd waves, reduces parameter-space dimension enough for a simple RBF emulator, and supplies analytic partial-wave formulae and triton benchmarks that the community can adopt. The emulator validation (App. D) and held-out predictions at 200 MeV and for Cij are concrete strengths. If the determinability claim survives once longer-range 3NF pieces are restored, the work would materially ease the LEC-fitting bottleneck for precision chiral 3NFs. Even as an exploratory study the basis and methodology are immediately useful to LENPIC-style analyses.","major_comments":[{"comment":"Abstract and Sec. III C / Table V: the central claim that “9 of 13 LECs can be reliably determined” is stronger than the internal evidence supports. Progressive columns of Table V show order-one and sign-flipping shifts that exceed the quoted 1σ errors (e.g. s6: −6.461±0.487 → +6.294±0.472; s5: −0.505 → +1.464; s10: −1.343 → +0.823 when the 1/2− sector is opened). The authors themselves note that changes typically exceed statistical errors and that systematics are unquantified, and the full-fit χ²/datum = 4.44 implies the covariance used for “statistical errors” is mis-scaled under an incomplete model. The language should be softened to “constrainable / determinable within the present incomplete Hamiltonian,” with an explicit caveat that values will shift once N3LO/N4LO long- and intermediate-range 3NFs are restored (already flagged in Secs. III A, III C, IV).","section":"Abstract; Sec. III C; Table V"},{"comment":"Sec. II and Table V full fit: the adopted natural window |si| ≲ 2 (motivated by Table I and NN experience) is violated by s6 ≈ 6 (and, in intermediate fits, |s9| ≈ 3.5–4). While the authors correctly state that naturalness cannot be judged until longer-range 3NFs are included, the abstract and summary still present the nine LECs as “reliably” fixed of natural size. Either the naturalness discussion should be moved into the main claim language, or the fit should be repeated with a soft naturalness prior / penalty so that the reader can see how much of the χ² improvement is purchased by unnaturally large si.","section":"Sec. II; Sec. III C; Table V; Sec. IV"},{"comment":"Sec. III C data selection and Eq. (7): fits omit Coulomb, restrict θcm ∈ [60°, 160°], drop 200 MeV and Cij from the χ², and retain a database that has not been cleaned for mutual inconsistencies (the KVI vs RIKEN 135 MeV cross-section tension is discussed but not resolved by a 3σ-style rejection). With χ²/N ≈ 4.4 it is unclear how much of the remaining discrepancy is model incompleteness versus outlier data. A short robustness check—refitting after a transparent outlier cut or after inflating experimental errors to force χ²/N ≈ 1—would make the “9 of 13” count more credible even inside the exploratory setup.","section":"Sec. III C; Table III; Eq. (7)"}],"minor_comments":[{"comment":"Table I: the text notes that ⟨VE5,6,9…13⟩ differ from Ref. [26]; a one-sentence remark on whether the discrepancy is traced to antisymmetrization, regulator, or a typographical issue in [26] would help readers who rely on both papers.","section":"Table I; Sec. II"},{"comment":"Figs. 1–8 lower panels: the fixed-angle slices are useful, but the chosen angle (open square) is not always near the maximum sensitivity; a brief note on the selection criterion would avoid the impression of cherry-picking.","section":"Sec. III A; Figs. 1–8"},{"comment":"Eq. (5) and the discussion of s2 redundancy: the Born-level argument that s2 is absorbable into E and s3 is clear; stating explicitly that s2 was fixed to zero only after verifying Δ(χ²/datum) ≲ 1% for |s2| = 2 would make the exclusion fully reproducible.","section":"Sec. II; Sec. III C"},{"comment":"Appendix A is long and valuable; a compact machine-readable supplement (or a note that the expressions are available on request / in a repository) would increase uptake.","section":"Appendix A"},{"comment":"Minor typography: “APPLICA TIONS”, “SCA TTERING” (Sec. III heading) and occasional missing spaces around math in the arXiv text should be cleaned in production.","section":"Sec. III"},{"comment":"Fig. 16 correlation matrix: the strongest correlation (cD–s11 ≈ 0.71) is mentioned in the text; adding the numerical value in the caption would help skimmers.","section":"Fig. 16; Sec. III C"}],"recommendation":"minor_revision","confidential_remarks":"The spectroscopic-basis contribution is solid and the paper is a natural fit for a nuclear-theory journal. The main risk is over-selling the fit results in the abstract; once the reliability language is tempered to match the exploratory, incomplete-Hamiltonian setting, the manuscript is close to acceptable. I do not see a load-bearing calculational error. The skeptic’s concern about progressive LEC instability is real but is already partially acknowledged by the authors and is fixable by wording and a short robustness note rather than new physics."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The real advance here is the spectroscopic reparametrization of the N4LO contact 3NF (Eq. 3, Table II) plus the channel-by-channel RBF emulator that makes multi-LEC Nd fits cheap. That is not just cosmetics: Si map cleanly onto 2S/2P/4P waves and mixing angles, the inverse map is given, and Appendix A supplies independent partial-wave formulas with triton benchmarks that correct some prior Ei expectation values. For anyone fitting 3NFs in the LENPIC orbit, this is practical infrastructure.\n\nThe sensitivity plots and progressive fits are done carefully. Emulator QA (App. D) shows interpolation error orders of magnitude below the 3NF signal. Held-out 200 MeV analyzing powers and Cij look reasonable. s2 and s4 are rightly dropped for elastic data; s12,s13 are T=3/2 and out of scope. Nine LECs do move the χ² and improve Ay and several tensors relative to NN + N2LO 2π only.\n\nThe soft spot is the reliability language, not the method. Table V already shows order-one and sign-flipping shifts when sectors are opened (s6 −6.5 → +6.3, s5 and s10 flip), larger than the quoted statistical errors. χ²/datum ends at 4.44; s6 sits well outside the |si|≲2 natural window they advertise. The Hamiltonian is incomplete—no N3LO/N4LO long-range 3NF, no Coulomb, angle cuts—and the authors say so and call the study exploratory. So “reliably determined” oversells what the internal evidence supports; the LECs will move again when the longer-range pieces come in. That does not kill the paper; it just means the abstract’s count should be read as “constrainable in this setup,” not “fixed for the full EFT.”\n\nMath, citations, and few-body methodology look solid. Self-citation is normal for this group’s ongoing program. This is for people who actually fit 3NFs or build Nd emulators—not a general nuclear-structure audience.\n\nI would send it to peer review. Engage with the basis and the emulator; treat the LEC values as provisional until N3LO long-range 3NF is restored.","headline":"The spectroscopic basis and emulator are genuinely useful; the “9 of 13 reliably determined” claim is softer than the abstract suggests once you look at Table V’s progressive shifts.","tokens_in":39610,"tokens_out":589,"would_cite":true,"duration_ms":18934,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"Rewriting short-range three-nucleon forces in a spectroscopic basis lets nine of thirteen low-energy constants be fixed from elastic nucleon-deuteron data.","keywords":["three-nucleon force","chiral EFT","spectroscopic basis","nucleon-deuteron scattering","low-energy constants","contact interactions","radial basis function emulator"],"falsifier":"Repeat the same elastic Nd fits after the symmetry-preserving N3LO (and N4LO long-range) three-nucleon contributions are included; if several of the nine Si then become unconstrained or wildly unnatural, the claim fails.","tokens_in":39346,"feed_emoji":"⚛️","tokens_out":976,"duration_ms":25877,"temperature":0.7,"pith_summary":"Three-nucleon short-range forces at next-to-next-to-next-to-leading order involve thirteen low-energy constants that are hard to pin down because each operator mixes many partial waves. This paper rewrites those operators as a spectroscopic basis in which each new constant controls only one total-angular-momentum and parity channel in nucleon-deuteron scattering. The reduced parameter space makes a simple radial-basis-function emulator practical, so the constants can be fit directly to elastic cross sections and analyzing powers. Exploratory fits show that nine of the thirteen constants are reliably determined by existing elastic data; two are largely redundant for elastic observables and two live only in the isospin-3/2 sector. The rewrite therefore turns an intractable multiparameter problem into a transparent, channel-by-channel determination that can guide the construction of precision three-nucleon forces.","feed_headline":"Nine of 13 three-nucleon constants fixed from Nd data","feed_subtitle":"A spectroscopic rewrite turns a 13-parameter mess into channel-by-channel fits that elastic scattering can actually constrain.","key_machinery":"The spectroscopic basis: an invertible linear map from the original Ei operators onto new constants Si that each act in only one JP partial-wave sector of Nd scattering (Table II and Eqs. 3–5), collapsing the fit space enough for radial-basis-function interpolation.","core_discovery":"Once the thirteen subleading contact three-nucleon operators are linearly transformed into the spectroscopic basis S1…S13, each Si contributes to a single JP channel (or vanishes) in elastic nucleon-deuteron scattering. With that basis the authors build an RBF emulator and demonstrate that nine of the thirteen dimensionless coefficients can be extracted from elastic Nd data at 10, 70 and 135 MeV, while s2 and s4 remain essentially unconstrained by elastic observables and s12, s13 are invisible in the T=1/2 system.","pith_inferences":["The same spectroscopic map should make eigenvector-continuation or Woodbury emulators for three-nucleon scattering dramatically cheaper, because each channel now depends on only a handful of LECs.","Once s12 and s13 are fixed from A=4 or neutron-matter calculations, the full contact 3NF can be frozen and used as a controlled input for medium-mass structure, testing whether the remaining discrepancies are truly long-range.","The opposite impact of the 3/2− LECs on nucleon versus deuteron Ay suggests a clean experimental program that measures both analyzing powers at the same energy to isolate mixing angles."],"forward_implications":["Elastic Nd data alone can fix the isospin-1/2 short-range 3NF sector once the spectroscopic rewrite is used.","The Ay puzzle at low energy is at least half-resolved by a natural-size s11 together with the 3/2− LECs.","The N2LO contact cD shifts negative once N4LO contacts are free, moving closer to the value preferred by tritium beta decay.","s2 and s4 must be constrained by breakup or four-body data; s12 and s13 require systems with T=3/2 components.","A complete N3LO Nd analysis can reuse the same spectroscopic emulator after the enhanced linear combinations of Si and Fi are restored."],"fun_headline_variants":["Spectroscopic basis pins nine of 13 three-nucleon constants from Nd data","Nine of 13 subleading 3N constants fixed by elastic Nd scattering","New basis lets elastic Nd data determine nine 3N low-energy constants","Channel-by-channel fits constrain nine short-range three-nucleon forces","Elastic nucleon-deuteron data constrains nine of thirteen 3N contacts"],"cache_read_input_tokens":32896,"weakest_assumption_plain":"That the nine constants remain determinable and of roughly natural size after the still-missing longer-range three-nucleon pieces beyond N2LO are restored to the Hamiltonian.","fun_headline_variants_meta":{"raw":{"variants":["Spectroscopic basis pins nine of 13 three-nucleon constants from Nd data","Nine of 13 subleading 3N constants fixed by elastic Nd scattering","New basis lets elastic Nd data determine nine 3N low-energy constants","Channel-by-channel fits constrain nine short-range three-nucleon forces","Elastic nucleon-deuteron data constrains nine of thirteen 3N contacts"]},"model":"grok-4.5","effort":"low","cost_usd":0.004795,"raw_usage":{"total_tokens":1339,"prompt_tokens":701,"num_sources_used":0,"completion_tokens":107,"cost_in_usd_ticks":47948000,"prompt_tokens_details":{"text_tokens":701,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":531,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":701,"tokens_out":107,"duration_ms":9808,"temperature":1.0,"reasoning_tokens":531,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-30T15:49:26.940955+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"Repeat the same elastic Nd fits after the symmetry-preserving N3LO (and N4LO long-range) three-nucleon contributions are included; if several of the nine Si then become unconstrained or wildly unnatural, the claim fails.","supporting_citations":[],"review_version":1}