{"id":"210f33f9-ba1d-4c4c-abe2-94dcf84f8aea","arxiv_id":"2509.08786","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"First fully microscopic, relativistic calculation of the 2p2h MEC contribution to semi-inclusive neutrino-carbon (νμ, μ−p) cross sections, with T2K flux-folded predictions.","lead":"Physicists have computed, for the first time, the full relativistic two-particle-two-hole contribution to semi-inclusive neutrino-carbon scattering, including pion and Delta meson-exchange currents. This matters because T2K and other neutrino experiments need accurate models of these two-nucleon events to measure neutrino oscillations.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Inclusive validation uses a reduced form-factor set, while the semi-inclusive predictions use the full set; the paper itself says the difference is non-negligible.","rationale":"The reader's verdict is CONDITIONAL and already cites the reduced form-factor set as a secondary concern. My independent stress-test identifies this as the single most load-bearing issue: the validation section validates a different operator set than the one used to produce the central semi-inclusive predictions. This is not just a missing error bar; it is an internal mismatch between the checked object and the claimed result. The reader's weakest_assumption focuses on the RFG+energy-shift model, which is a related but distinct concern: the RFG may be inadequate for semi-inclusive observables, and no 2p2h semi-inclusive data exist to test it. Both concerns point to conditional acceptance, so the verdict is unchanged. The concrete test is a targeted numerical check that would resolve the form-factor mismatch: if the full-form-factor inclusive responses agree with Ref. [31]/[5] within a few percent, the mismatch is benign; if they differ, the central predictions need a caveat or revalidation. I do not see a reason to reject the paper outright, as the formalism and first-calculation claim appear sound and the issue is remediable.","tokens_in":18314,"tokens_out":7562,"duration_ms":92809,"concrete_test":"Recompute the inclusive νμ-12C 2p2h responses of Sec. 4.1 at |q|=400 and 1000 MeV using the full Delta form-factor set (all CiV/A in Eqs. (22)–(23)) instead of only C3V and C5A, and compare with Ref. [31]/[5]. If the CC, CL, or LL responses shift by more than a few percent, then the semi-inclusive results of Sec. 4.2—which use the full set—are not validated by the reported inclusive comparison. In that case, also recompute the flux-folded cross section of Fig. 10 with the reduced set to quantify the expected shift and assess whether the comparison with GENIE is robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central numerical claim rests on the full MEC with all Delta form factors, but the only quantitative validation in Sec. 4.1 and Figs. 3–5 is performed with a restricted current: only C3V and C5A are retained, 'consistently with Ref. [31]'. The text immediately after Fig. 5 states that the omitted subdominant form factors have a non-negligible effect, especially in CC, CL, and LL responses. Yet the semi-inclusive cross section (Eq. 9) and the flux-folded predictions in Figs. 9–10 use the full form-factor set. Therefore the model actually used for the headline predictions is not the model that was validated. The comparison is also circular: it is against the same group's previous reduced-form-factor calculation, not against independent inclusive data or an independent many-body calculation. This leaves the semi-inclusive results without a demonstrated quantitative check, which weakens the claim that this calculation provides a reliable replacement for the inclusive-based GENIE extraction. The novelty of being first to compute this observable is not undermined, but the numerical support for the specific comparison with GENIE is thinner than the word 'validated' implies.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a calculation of the two-particle two-hole (2p2h) contribution to charged-current semi-inclusive neutrino-carbon scattering, (ν_μ, μ^- p), in the relativistic Fermi gas model. The two-body current includes pion-in-flight, seagull, pion-pole, and Δ forward/backward terms. The formalism is adapted from the authors' earlier (e,e'p) study, and the paper validates the inclusive 2p2h response functions against previous results from the same group. It then presents six-fold differential semi-inclusive cross sections, their decomposition into Δ, pionic, interference, and isospin channels, and flux-folded T2K 1μCC0πNp predictions compared with the GENIE/SuSAv2-MEC implementation. The authors claim the first fully microscopic, relativistic, semi-inclusive 2p2h computation, and argue that event generators should move from inclusive-based 2p2h extraction to such a model.","tokens_in":18607,"tokens_out":5732,"duration_ms":67585,"significance":"If the numerical results are reliable, this is a useful and timely contribution: event generators currently implement semi-inclusive 2p2h by making uncontrolled assumptions on top of inclusive calculations. The paper's strengths are the detailed treatment of the ten semi-inclusive response functions, a clear discussion of the leading-proton versus detected-proton observables, and a documented reduction of the nine-dimensional phase space to five dimensions. The claim of novelty appears justified, since previous semi-inclusive 2p2h studies either omitted Δ currents or computed different exclusive observables. However, the quantitative support for the headline predictions is thinner than the word 'validated' suggests: the inclusive check is performed with a reduced form-factor set, while the semi-inclusive results use the full set; and the final comparison is against another model, not against data or an independent calculation.","major_comments":[{"comment":"The inclusive validation is performed with a reduced current: 'only the dominant form factors C3V, C5A, consistently with Ref. [31]'. The text immediately after Fig. 5 states that the subdominant form factors have a non-negligible effect, especially in the CC, CL and LL responses. Yet the semi-inclusive results in Eq. (9) and Figs. 6-10 are computed with all form factors. Thus the model that is validated is not the model used for the headline semi-inclusive predictions. At minimum, the authors should show an inclusive benchmark with the full form-factor set, or quantify the effect of the omitted form factors on the semi-inclusive observables. Without this, the GENIE comparison in Fig. 10 lacks a demonstrated quantitative anchor.","section":"Section 4.1, Figs. 3-5 and text after Fig. 5"},{"comment":"The comparison with GENIE is not a validation. The GENIE/SuSAv2-MEC curve is itself an approximate extraction from an inclusive 2p2h model, supplemented with FSI, and the paper explicitly attributes the observed discrepancies to FSI and to the inclusive-to-exclusive extraction assumptions. The statement that event generators should adopt the present model is therefore a plausible proposal rather than a conclusion supported by the comparison. The authors should either compare their full-model 2p2h contribution with the T2K 1μCC0πNp data (including uncertainties) or with an independent microscopic calculation; otherwise the central practical claim is undersupported.","section":"Section 4.2, Fig. 10"},{"comment":"The energy-shift prescription is ambiguous. The text defines ω~ = ω - E_shift, then states that for 12C the chosen values are pF = 225 MeV, E_shift = 20 MeV, and E_2p2h_shift = 2 E_shift = 40 MeV. It is not stated explicitly whether Eq. (24) is used with E_shift or with E_2p2h_shift when evaluating the 2p2h tensor in Eq. (10). This distinction changes the phase space and the resulting cross sections. Please state explicitly which effective energy transfer enters the 2p2h calculation and justify the factor of two.","section":"Section 3, Eq. (24)"}],"minor_comments":[{"comment":"The caption of Fig. 5 says the comparison is with the computation of Ref. [5], while the main text says the comparison is with responses presented in Ref. [31]; Section 5 again refers to Ref. [5]. Please harmonize the references.","section":"Fig. 5 and Section 5"},{"comment":"Typo: 'As a results' should be 'As a result'.","section":"Section 1"},{"comment":"The theta functions contain a typographical imbalance: θ(|p1| - pF|) should presumably be θ(|p1| - pF). Please correct.","section":"Eq. (10)"},{"comment":"The text says 'To make contact with experimental data', but no T2K data are shown. The current Figs. 9-10 show only the present calculation and the GENIE prediction. Either overlay the T2K data with uncertainties or reword the claim.","section":"Section 4.2"},{"comment":"The novelty claim should be scoped carefully: since Ref. [18] already presents semi-inclusive 2p2h results, the phrase 'first fully microscopic, relativistic, semi-inclusive computation' should explicitly distinguish the inclusion of Δ currents and the relativistic treatment from Ref. [18], to avoid an easily contested priority statement.","section":"Section 4.2, claim of 'first'"}],"recommendation":"major_revision","confidential_remarks":"The central derivation appears careful and the phase-space reduction is a real strength. My main concern is that the numerical model used for the semi-inclusive predictions (full form-factor set) is not the model checked in the inclusive validation (C3V, C5A only), and the text itself concedes the difference is non-negligible. If the authors can supply a full-form-factor inclusive benchmark, or quantify the effect on the semi-inclusive observables, the paper would be substantially strengthened. The GENIE comparison is not an independent benchmark. The novelty of the calculation seems genuine, and I do not see evidence of duplicate publication or missing prior work beyond the scoping point about Ref. [18]."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a genuine step forward. The authors extend their electron-scattering semi-inclusive 2p2h formalism to the weak sector, producing the first relativistic microscopic calculation of semi-inclusive (νμ, μ−p) 2p2h cross sections with Delta meson-exchange currents. That novelty claim holds up against the cited literature. The treatment of the ten response functions, the axial-vector interference, and the isospin sums for a detected proton is careful, and the reduction of the nine-dimensional phase-space integral to five dimensions is documented in enough detail to be reproducible. The pp/pn decomposition and the leading-proton definition study are practically useful for generator tuning. The paper is honest about what the model omits (FSI, correlations, rho meson) and about its computational cost.\n\nThe soft spots are real but not fatal. The inclusive validation in Sec. 4.1 is performed with only the dominant form factors C3V and C5A, deliberately matching Ref. [31]. The semi-inclusive results in Figs. 6–10 use the full set of Delta form factors. The paper explicitly states the subdominant form factors have a non-negligible effect, especially in CC, CL, LL responses. So the model that is validated is not the model that generates the headline predictions. This is not a hidden flaw—it is admitted in the text—but it does mean the quantitative agreement with Ref. [31] cannot be read as a check of the full calculation. The comparison with GENIE in Fig. 10 is suggestive, not a test.\n\nAlso, the RFG+constant-shift model is known to be inadequate for semi-inclusive QE; the authors argue that the two-body current introduces correlations that soften this, but they present no independent test of that premise. Combined with the absence of error bars, the predictions should be treated as the first map of an unmapped observable, not as a final benchmark.\n\nWho is this for: anyone working on neutrino event generators or semi-inclusive measurements at T2K, MINERvA, MicroBooNE, DUNE. They will get a concrete model to compare against and a clear target for improvement.\n\nRecommendation: send it to peer review. A good referee should ask for a validation of the full-form-factor model—for instance, inclusive responses with the full set against independent calculations or data—and for a sensitivity study of the energy shift. But the calculation is new, relevant, and reproducible enough to merit referee time.","headline":"First semi-inclusive 2p2h neutrino cross section with Delta MEC is a real advance, but the validation covers a reduced form-factor version of the model, not the full current used for the predictions.","tokens_in":19116,"tokens_out":2973,"would_cite":true,"duration_ms":32964,"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":"The paper provides the first fully microscopic, relativistic calculation of the two-particle two-hole (2p2h) contribution to semi-inclusive charged-current neutrino-nucleus scattering, specifically the (νμ, μ−p) channel on carbon.","keywords":["two-particle two-hole excitations","semi-inclusive neutrino scattering","meson-exchange currents","relativistic Fermi gas","neutrino-nucleus cross sections","delta resonance","charged-current interactions","T2K"],"falsifier":"A measurement of the T2K 1μCC0πNp cross section as a function of leading proton momentum would settle the matter: if the data follow the inclusive-based Monte Carlo extraction rather than this calculation's peak position and strength, or if the measured pp/pn ratio departs strongly from the predicted value of about 4, the central claim would lose support.","tokens_in":18230,"feed_emoji":"⚛️","tokens_out":4105,"duration_ms":48669,"temperature":0.7,"pith_summary":"The paper extends a relativistic Fermi gas model with meson-exchange currents to compute, for the first time, the full two-particle two-hole (2p2h) contribution to semi-inclusive charged-current neutrino-carbon scattering (νμ, μ−p). It validates the calculation against inclusive neutrino response functions, then produces proton momentum and angle dependent cross sections at fixed energies and folded with the T2K flux. The main message is that a fully microscopic semi-inclusive 2p2h calculation is feasible and differs from the approximations currently used in event generators, both in peak position and strength. A reader should care because semi-inclusive measurements are sensitive to exactly these differences, and a reliable 2p2h model is needed to interpret them.","feed_headline":"First microscopic model of two-nucleon knockout in neutrino scattering","feed_subtitle":"Semi-inclusive proton spectra from two-nucleon currents are computed from first principles; event generators currently rely on inclusive app","key_machinery":"The central object is the semi-inclusive 2p2h nuclear tensor built from a two-body meson-exchange current operator. The current includes pion-in-flight, seagull, pion-pole, and Δ forward/backward contributions, with hadronic and weak form factors. The nuclear tensor integrates over Fermi-gas hole states with Pauli blocking and an energy shift, and the nine-dimensional phase-space integral is reduced analytically to five dimensions. Because azimuthal symmetry is broken in the semi-inclusive process, the cross section is expressed in terms of ten response functions instead of the five used in inclusive scattering.","core_discovery":"The paper claims to be the first to evaluate the 2p2h impact on semi-inclusive neutrino-nucleus scattering from a fully microscopic, relativistic, semi-inclusive computation. Using the relativistic Fermi gas with an energy shift, it builds the two-body nuclear tensor from meson-exchange currents including pion-in-flight, seagull, pion-pole, and Δ forward/backward terms. The model is first checked against inclusive 2p2h response functions, then applied to the semi-inclusive (νμ, μ−p) cross section on 12C. The resulting proton spectra, decomposed into pp and pn channels and into pionic, Δ, and interference contributions, show that the pp channel dominates and the π−Δ interference can change si","pith_inferences":["The constant energy shift fitted to inclusive electron-scattering data may need to be re-fitted once semi-inclusive 2p2h data become available; the formalism allows this without structural changes, and the resulting shift would quantify medium effects in the two-body channel.","The sign-changing π−Δ interference at backward proton angles is a distinctive prediction that dedicated angular measurements could confirm or rule out, providing a sharper test than momentum-integrated cross sections.","If the pp/pn ratio is measured and departs strongly from the predicted value of about 4, that would indicate missing physics—likely short-range correlations—beyond the Fermi-gas basis, as the paper itself suggests when comparing with other models.","The same formalism could be extended to exclusive (νμ, μ−NN) observables, where the two final nucleons are detected together; the leading-proton variable used here is a natural bridge to that more exclusive phase space."],"forward_implications":["If the calculation is correct, event generators should adopt a fully microscopic semi-inclusive 2p2h model rather than extracting exclusive predictions from inclusive results, since the two approaches give visibly different peak positions and strengths.","The predicted dominance of the pp final-state channel, with a pp/pn ratio of about 4, provides a concrete, testable signature in proton-tagged measurements.","The Δ contribution accounts for roughly 80% of the total semi-inclusive 2p2h strength, and the π−Δ interference is non-negligible and can change sign with kinematics, so models omitting either piece are likely to be inaccurate.","The computation is efficient enough for practical use in event generators: for a given four-momentum transfer, the differential cross section can be evaluated to 1% numerical accuracy in about one minute of CPU time.","The model gives a baseline for understanding final-state interactions: since FSI are not included, comparison with data that include FSI should reveal the size of those effects on the proton momentum distribution."],"fun_headline_variants":["First semi-inclusive 2p2h model for neutrino scattering","Two-nucleon knockout: first semi-inclusive ν-nucleus calc","Relativistic MEC model for semi-inclusive two-nucleon knockout","First proton spectra from two-nucleon currents in ν scattering","Semi-inclusive 2p2h: first microscopic neutrino-nucleus result"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The calculation's reliability rests on the assumption that the relativistic Fermi gas with a constant energy shift—fitted to inclusive electron scattering—faithfully represents the initial two-nucleon system for semi-inclusive 2p2h kinematics, even though the same model is known to be inadequate for semi-inclusive one-body knockout.","fun_headline_variants_meta":{"raw":{"variants":["First semi-inclusive 2p2h model for neutrino scattering","Two-nucleon knockout: first semi-inclusive ν-nucleus calc","Relativistic MEC model for semi-inclusive two-nucleon knockout","First proton spectra from two-nucleon currents in ν scattering","Semi-inclusive 2p2h: first microscopic neutrino-nucleus result"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000234,"raw_usage":{"total_tokens":1282,"prompt_tokens":644,"completion_tokens":638,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":388,"completion_tokens_details":{"reasoning_tokens":540}},"tokens_in":388,"tokens_out":638,"duration_ms":7575,"temperature":1.0,"reasoning_tokens":540,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T20:08:26.543246+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A measurement of the T2K 1μCC0πNp cross section as a function of leading proton momentum would settle the matter: if the data follow the inclusive-based Monte Carlo extraction rather than this calculation's peak position and strength, or if the measured pp/pn ratio departs strongly from the predicted value of about 4, the central claim would lose support.","supporting_citations":[],"review_version":1}