REVIEW 3 major objections 5 minor 41 references
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.
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 →
T0 review · deepseek-v4-flash
2026-08-04 20:08 UTC pith:L7DAOKSR
load-bearing objection 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. the 3 major comments →
Two-particle two-hole excitations in semi-inclusive neutrino-nucleus scattering
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
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
What carries the argument
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.
Load-bearing premise
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.
What would settle it
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.
If this is right
- 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.
Where Pith is reading between the lines
- 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.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
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.
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 (3)
- [Section 4.1, Figs. 3-5 and text after Fig. 5] 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 4.2, Fig. 10] 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 3, Eq. (24)] 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.
minor comments (5)
- [Fig. 5 and Section 5] 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 1] Typo: 'As a results' should be 'As a result'.
- [Eq. (10)] The theta functions contain a typographical imbalance: θ(|p1| - pF|) should presumably be θ(|p1| - pF). Please correct.
- [Section 4.2] 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 4.2, claim of 'first'] 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.
Circularity Check
No significant circularity: the semi-inclusive 2p2h cross section is computed directly from the MEC matrix elements and RFG phase space; no fitted input or self-cited result forces the headline prediction.
full rationale
The derivation chain is self-contained. The semi-inclusive tensor W_{A(N),2p2h} is computed from Eq. (10) using the RFG phase space and the MEC matrix elements of Eq. (12), with the current operators, couplings, and form factors taken from Refs. [27,28]. The only model parameters, p_F and E_shift, are fixed by electron-scattering data and the de Forest prescription, not by the neutrino semi-inclusive observable being predicted. No semi-inclusive datum is fitted, and the GENIE comparison in Fig. 10 is an external, non-fitted check. The inclusive validation in Sec. 4.1 is an internal consistency check: Eq. (26) defines the inclusive tensor as the integral of the semi-inclusive tensor, so comparing with the authors' earlier inclusive results in Refs. [5]/[31] verifies the phase-space integration and computation but does not enter as an input into the semi-inclusive calculation. This is self-citation, but it is not load-bearing for the central derivation. The one flagged limitation is in Sec. 4.1 after Fig. 5: "In these results the ∆ interaction vertex is described using only the dominant form factors, to be consistent with Ref. [31]. Note, however, that in the other results presented in this work all the ∆ form factors ... are included. The effect of the subdominant form factors is not negligible..." That is a validation gap for the full-current predictions, not a circular step: the full MEC prediction is still computed from the model rather than extracted from the reduced-form-factor benchmark. No self-definitional, fitted-input-called-prediction, or uniqueness-imported-from-authors circularity is present.
Axiom & Free-Parameter Ledger
free parameters (2)
- Fermi momentum p_F =
225 MeV
- Energy shift E_shift =
20 MeV
axioms (4)
- domain assumption The nuclear tensor is approximated as the sum of 1p1h and 2p2h contributions, omitting many-body currents beyond two-body and medium corrections to self-energy.
- domain assumption The relativistic Fermi gas with a constant energy shift provides an adequate framework for 2p2h semi-inclusive responses.
- domain assumption The MEC is given by the chiral Lagrangian of Ref. [27] with the form factors and Delta propagator of Ref. [28].
- standard math The de Forest prescription is applied: form factors are evaluated at the exact omega, while other terms in the hadronic tensor use the energy-shifted omega.
read the original abstract
A calculation of the two-particle two-hole contribution to the semi-inclusive $(\nu_\mu,\mu^-N)$ cross section on carbon is performed in the framework of the relativistic Fermi gas model. The process is driven by meson-exchange currents encompassing all contributions involving the exchange of a single pion and the excitation of a $\Delta$ resonance. The calculation is validated against the inclusive $(\nu_\mu,\mu^-)$ response functions already existing in the literature, and then extended for the first time to the semi-inclusive channel. Results are presented both at fixed neutrino energy and folded with the T2K neutrino flux.
Figures
Reference graph
Works this paper leans on
-
[1]
INTRODUCTION Understanding and accurately describing the properties of neutrinos is one of the primary and most challenging goals of the neutrino physics community, nowadays especially focused in the precise determination of the PMNS oscillation parameters. Particular attention is paid to the CP violating Dirac phase, the main purpose of the current long-...
2009
-
[2]
[24] for semi-inclusive electron-nucleus scattering is extended to describe the charged-current (CC) semi-inclusive neutrino-nucleus interactions
FORMALISM In this Section, the formalism developed in Ref. [24] for semi-inclusive electron-nucleus scattering is extended to describe the charged-current (CC) semi-inclusive neutrino-nucleus interactions. In a semi-inclusive process one or more hadrons in the final state are detected in coin- cidence with the outgoing lepton. We focus here on the specifi...
-
[3]
effective
MODEL To compute the 2p2hhadronic tensor, we adopt as wave-function basis the Relativistic Fermi Gas (RFG) model, which treats the nucleus as a system of non-interacting nucleons described by Dirac spinors, correlated only by the Pauli exclusion principle. It is worth mentioning that, in the QE channel, the RFG model is inadequate to provide semi-inclusiv...
-
[4]
parallel kinematics
RESUL TS In this section we present the numerical results of the calculation illustrated above. Before showing the semi-inclusive 2p2hcross sections, we validate our calculation by evaluating the inclusive weak 2p2hresponse functions and comparing the results with the ones published in Ref. [31]. The values of the parameters used in this work are listed i...
-
[5]
Planes Complementarios de I+D+i
CONCLUSIONS Motivated by the promising results presented in [24] for the 2p2hcontribution to the electromagnetic (e, e′p) process, in this work we have extended the calculation to the weak sector and made predictions for neutrino-carbon semi-inclusive (ν µ, µ−p) observables. The weak model has been first validated by comparing with the inclusive results o...
2022
-
[6]
A. A. Aguilar-Arevalo et al. (MiniBooNE), Phys. Rev. D81, 092005 (2010), 1002.2680
Pith/arXiv arXiv 2010
-
[7]
M. Martini, M. Ericson, G. Chanfray, and J. Marteau, Phys. Rev. C80, 065501 (2009), 0910.2622
Pith/arXiv arXiv 2009
-
[8]
M. Martini, M. Ericson, G. Chanfray, and J. Marteau, Phys. Rev. C81, 045502 (2010), 1002.4538
Pith/arXiv arXiv 2010
-
[9]
J. Nieves, I. Ruiz Simo, and M. J. Vicente Vacas, Phys. Rev. C83, 045501 (2011), 1102.2777
Pith/arXiv arXiv 2011
-
[10]
I. Ruiz Simo, J. E. Amaro, M. B. Barbaro, A. De Pace, J. A. Caballero, and T. W. Donnelly, J. Phys. G44, 065105 (2017), 1604.08423
Pith/arXiv arXiv 2017
-
[11]
T. Van Cuyck, N. Jachowicz, R. Gonz´ alez-Jim´ enez, J. Ryckebusch, and N. Van Dessel, Phys. Rev. C95, 054611 (2017), 1702.06402
Pith/arXiv arXiv 2017
-
[12]
N. Rocco, C. Barbieri, O. Benhar, A. De Pace, and A. Lovato, Phys. Rev. C99, 025502 (2019), 1810.07647
Pith/arXiv arXiv 2019
-
[13]
A. Lovato, J. Carlson, S. Gandolfi, N. Rocco, and R. Schiavilla, Phys. Rev. X10, 031068 (2020), 2003.07710
Pith/arXiv arXiv 2020
- [14]
-
[15]
X. G. Lu et al. (MINERvA), Phys. Rev. Lett.121, 022504 (2018), 1805.05486
Pith/arXiv arXiv 2018
- [16]
-
[17]
P. Abratenko et al. (MicroBooNE), Phys. Rev. Lett.125, 201803 (2020), 2006.00108
Pith/arXiv arXiv 2020
-
[18]
P. Abratenko et al. (MicroBooNE), Phys. Rev. D102, 112013 (2020), 2010.02390
Pith/arXiv arXiv 2020
-
[19]
P. Abratenko et al. (MicroBooNE), Phys. Rev. Lett.131, 101802 (2023), 2301.03706
Pith/arXiv arXiv 2023
-
[20]
O. Moreno, T. W. Donnelly, J. W. Van Orden, and W. P. Ford, Phys. Rev. D90, 013014 34 (2014), 1406.4494
Pith/arXiv arXiv 2014
-
[21]
J. M. Franco-Patino, R. Gonz´ alez-Jim´ enez, S. Dolan, M. B. Barbaro, J. A. Caballero, G. D. Megias, and J. M. Udias, Phys. Rev. D106, 113005 (2022), 2207.02086
Pith/arXiv arXiv 2022
-
[22]
J. M. Franco-Patino, S. Dolan, R. Gonz´ alez-Jim´ enez, M. B. Barbaro, J. A. Caballero, and G. D. Megias, Phys. Rev. D109, 013004 (2024), 2304.01916
Pith/arXiv arXiv 2024
-
[23]
T. Van Cuyck, N. Jachowicz, R. Gonz´ alez-Jim´ enez, M. Martini, V. Pandey, J. Ryckebusch, and N. Van Dessel, Phys. Rev. C94, 024611 (2016), 1606.00273
Pith/arXiv arXiv 2016
-
[24]
V. L. Martinez-Consentino, A. M. Cantizani, and J. E. Amaro, Phys. Rev. C109, 015502 (2024), 2310.12642
Pith/arXiv arXiv 2024
-
[25]
V. L. Martinez-Consentino and J. E. Amaro, Symmetry16(2024), ISSN 2073-8994, URL https://www.mdpi.com/2073-8994/16/2/247
2024
-
[26]
S. Dolan, G. D. Megias, and S. Bolognesi, Phys. Rev. D101, 033003 (2020), 1905.08556
Pith/arXiv arXiv 2020
-
[27]
Y. Hayato and L. Pickering, Eur. Phys. J. ST230, 4469 (2021), 2106.15809
Pith/arXiv arXiv 2021
-
[28]
Khachatryan et al
M. Khachatryan et al. (CLAS, e4v), Nature599, 565 (2021)
2021
-
[29]
V. Belocchi, M. B. Barbaro, A. De Pace, and M. Martini, Phys. Rev. C109, 065502 (2024), 2401.13640
Pith/arXiv arXiv 2024
-
[30]
P. R. Casale, J. E. Amaro, V. Belocchi, M. B. Barbaro, A. De Pace, and M. Martini, Phys. Rev. C112, 024603 (2025), 2503.08391
Pith/arXiv arXiv 2025
-
[31]
P. R. Casale, J. E. Amaro, V. Belocchi, M. B. Barbaro, and M. Martini (2025), 2507.20246
arXiv 2025
-
[32]
Scherer and M
S. Scherer and M. R. Schindler, A Primer for Chiral Perturbation Theory (2012)
2012
-
[33]
E. Hernandez, J. Nieves, and M. Valverde, Phys. Rev. D76, 033005 (2007), hep-ph/0701149
Pith/arXiv arXiv 2007
-
[34]
De Forest, Nuclear Physics A392, 232 (1983), ISSN 0375-9474, URLhttps://www
T. De Forest, Nuclear Physics A392, 232 (1983), ISSN 0375-9474, URLhttps://www. sciencedirect.com/science/article/pii/0375947483901240
arXiv 1983
-
[35]
C. Maieron, T. W. Donnelly, and I. Sick, Phys. Rev. C65, 025502 (2002), nucl-th/0109032
Pith/arXiv arXiv 2002
-
[36]
I. Ruiz Simo, J. E. Amaro, M. B. Barbaro, A. De Pace, J. A. Caballero, G. D. Megias, and T. W. Donnelly, Phys. Lett. B762, 124 (2016), 1607.08451
Pith/arXiv arXiv 2016
-
[37]
Belocchi, Ph.D
V. Belocchi, Ph.D. thesis, Universit` a degli Studi di Torino, Italy (2025)
2025
-
[38]
Niewczas, Ph.D
K. Niewczas, Ph.D. thesis, U. Gent (main) (2023)
2023
-
[39]
J. M. Franco-Patino, J. Gonzalez-Rosa, J. A. Caballero, and M. B. Barbaro, Phys. Rev. C 102, 064626 (2020), 2010.14937
Pith/arXiv arXiv 2020
-
[40]
M. J. Dekker, P. J. Brussaard, and J. A. Tjon, Phys. Rev. C49, 2650 (1994)
1994
-
[41]
E. E. van Faassen and J. A. Tjon, Phys. Rev. C28, 2354 (1983). 35
1983
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.