{"id":"e4768a5a-20b7-4476-90e7-061ad9f0fc7f","arxiv_id":"2412.18636","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"SuSAv2 model predictions are compared with NOvA and MicroBooNE inclusive charged-current neutrino cross sections, with good overall agreement for NOvA electron neutrinos and systematic underestimation for muon neutrinos and argon at higher energies.","lead":"This paper tests the SuSAv2 nuclear model against charged-current neutrino cross sections measured by NOvA and MicroBooNE on carbon- and argon-based detectors. It finds the model broadly matches NOvA electron-neutrino data but underestimates some muon-neutrino and MicroBooNE results, pointing to missing strength in high-energy inelastic channels.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Underestimation attributed to missing inelastic strength is not supported by the paper's own channel decomposition: MicroBooNE discrepancies occur where QE dominates, and NOvA νμ shows overprediction from DIS.","rationale":"The reader's weakest_assumption focuses on the superscaling factorization. That is a legitimate foundational concern, but the paper has previously validated SuSAv2 against electron and neutrino data on carbon, and the present data comparison is itself a test of the full model. The more directly load-bearing issue for the central claim is the causal attribution of the residual discrepancy to 'missing strength in the inelastic channels.' The paper's own figures and text undermine this attribution: MicroBooNE discrepancies appear in the QE-dominated total cross section, and NOvA νμ overpredicts from DIS. Without a systematic residual analysis or alternate-channel scaling test, the conclusion is not uniquely supported. This does not invalidate the model comparison—the χ2 values are useful—but it means the paper's central interpretive claim should be treated as conditional, which matches the reader's CONDITIONAL verdict. I therefore keep the verdict unchanged.","tokens_in":13867,"tokens_out":5722,"duration_ms":55641,"concrete_test":"Perform a residual analysis on the published figures: for each kinematic bin, compute the ratio (data - model)/model and the inelastic fraction (RES+SoftDIS+TrueDIS)/total from the model curves. Test whether the residual correlates positively with the inelastic fraction across bins. Then, using the NOvA νe data and covariance, refit the model with two free scale factors, s_QE for QE+MEC and s_inel for RES+SoftDIS+TrueDIS. If the best fit has s_QE ≈ 1 and s_inel > 1, the inelastic attribution is supported; if s_QE deviates or the correlation is absent, the attribution fails. This uses only the data and model already in the paper.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that the observed underestimation 'can be ascribed to some missing strength in the inelastic channels' is not established by the presented comparison and is partly contradicted by the paper's own channel decomposition. For MicroBooNE, the flux peaks at 0.8 GeV and the text states QE is the dominant channel; the total cross section is underestimated (Fig. 3 left), so missing inelastic strength cannot explain the discrepancy—the conclusion itself concedes this ('not expected to explain these discrepancies'). For NOvA νμ, the model overpredicts data at high muon kinetic energies, attributed to DIS (Fig. 2), the opposite sign of a 'missing strength'. The NOvA νe forward-angle bin has DIS contributions below 15% and QE/RES each ~35%, so the residual there is not obviously inelastic-dominated. Additionally, no theoretical uncertainty bands are given, and the quoted χ2 values lack degrees of freedom/covariance details, so 'reproduce well' is not quantified. Thus the causal attribution to inelastic channels is the weakest link: it requires the channel decomposition to be exact and the discrepancy to vanish when inelastic strength is increased, neither of which is demonstrated.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper confronts the SuSAv2 model, supplemented by RFG-based 2p2h-MEC, DCC resonances, and Bodek-Ritchie-based soft/true DIS contributions, with published CC-inclusive neutrino scattering data from NOvA (electron and muon neutrino double-differential cross sections on a mixed carbon-rich target) and MicroBooNE (muon neutrino total, single-differential, and double-differential cross sections on argon). The central positive claim is that the model reproduces the shape and magnitude of most measured distributions, with a residual underestimation at very forward angles and high lepton energies that the authors attribute to missing strength in the inelastic channels. The paper also reports chi-square values, including a shape-only factor for NOvA muon neutrinos, and discusses the use of a smearing matrix for the MicroBooNE comparison.","tokens_in":14134,"tokens_out":5526,"duration_ms":56344,"significance":"If the comparison is taken at face value, the paper provides a useful, wide-energy-range validation of the SuSAv2 framework on two different nuclear targets, including argon, and it produces a channel-by-channel decomposition that can inform neutrino event generators. A genuine strength is that no parameter is fitted to the NOvA or MicroBooNE data under comparison, so the exercise is a forward prediction rather than a fit. The use of published experimental data, the transparent listing of the reaction channels in Table I, and the explicit handling of the MicroBooNE smearing matrix are also assets. However, the quantitative support for the paper's main interpretive claim is weakened by the absence of theoretical uncertainty bands, the lack of statistical details (degrees of freedom, covariance definition), and several places where the channel attribution is either internally inconsistent or contradicted by the displayed decomposition.","major_comments":[{"comment":"The attribution of the MicroBooNE underestimation to missing inelastic strength is internally inconsistent. In Section III B the text says that in all panels of Fig. 3 a similar underestimation 'can be ascribed to some missing strength in the inelastic channels,' yet the same section and the conclusions state that at the MicroBooNE flux peak (around 0.8 GeV) the quasielastic channel dominates and that 'the possible lack of strength presented by our model in the inelastic channel is not expected to explain these discrepancies.' Since the underpredicted total and single-differential cross sections in Fig. 3 are dominated by QE at these kinematics, the paper must either remove this attribution for the MicroBooNE data of Ref. [56] or provide a quantitative channel-level decomposition, after smearing, showing that inelastic contributions are actually substantial in the discrepant kinematic region.","section":"Section III B and Section IV"},{"comment":"The central claim that the observed underestimation can be ascribed to missing inelastic strength is partly contradicted by the NOvA muon-neutrino results. In the discussion of Fig. 2 the paper states that at high muon kinetic energies the model overpredicts the data 'mostly due to the deep inelastic scattering contribution,' which is the opposite sign of a missing inelastic contribution. Moreover, for the NOvA electron-neutrino forward-angle bins the text reports that the DIS contributions fall below 15% and that QE and RES are each about 35% of the total, so the residual underprediction there is not obviously dominated by inelastic channels. The conclusion should be rephrased to specify, for each data set and kinematic region, which channel is over- or underpredicted and with what sign; a single blanket attribution to inelastic strength is not supported.","section":"Section III A (NOvA)"},{"comment":"The quantitative support for 'reproduce well' is not auditable as presented. The manuscript quotes chi-square values (e.g., 14.4 in Fig. 1, 3.0 for the NOvA nu-mu shape factor, 740.1 and 289.5 for MicroBooNE) without stating the number of data points, the degrees of freedom, or the precise covariance matrix used. No theoretical uncertainty bands are shown on any prediction, so differences between the model and data cannot be separated from model uncertainty. To support the paper's comparative claims, the authors should provide, for each fit, the number of points, the chi-square per degree of freedom, and a clear description of how the covariance matrices and the smearing regularization bias enter the calculation.","section":"Section III (Figs. 1-9)"},{"comment":"The superscaling factorization is assumed to hold across the entire energy spectrum, 'from quasielastic to deep inelastic scattering,' and is then applied to argon targets at NOvA and MicroBooNE kinematics. Since the SuSAv2 scaling function fSuSAv2 is extracted from electron-scattering data, largely on lighter nuclei and in the quasielastic region, the transfer to argon and to DIS kinematics is a substantial extrapolation. The paper should either provide validation for this extrapolation (for example, a comparison with inclusive electron-argon scattering data or a sensitivity study varying fSuSAv2) or explicitly list it as a limitation that weakens the conclusion about missing inelastic strength.","section":"Section II, Eq. (3)"}],"minor_comments":[{"comment":"The axis labels should be made consistent and explicit: Fig. 1 uses E_e while the text discusses electron energy, and Fig. 2 uses E_mu with the text referring to muon kinetic energy; the authors should state whether the plotted quantity is kinetic or total energy.","section":"Figures 1 and 2"},{"comment":"In the middle and right panels, the horizontal-axis labels (E_mu and omega) should be defined in the caption, since the text interchangeably refers to 'muon energy' and 'transferred energy.'","section":"Fig. 3"},{"comment":"The notation 'chi2_shape = chi2 * Fshape, where Fshape is the only-shape factor' is unclear; Fshape should be defined in the text, and the value 3.0 should be explained as a multiplicative factor or as a separate shape-only chi-square.","section":"Fig. 2 legend"},{"comment":"The sentence noting that the underestimation 'is not observed in a previous work [38]' is ambiguous because the following sentences discuss the difference between MicroBooNE data sets [3] and [56]; the relation between these two measurements should be clarified.","section":"Section III B"},{"comment":"The entry for SoftDIS ('SuSAv2 inelastic - SuSAv2-DCC') should be expanded to state that both terms are evaluated with the same kinematic limits and the same single-nucleon structure functions, so the difference isolates the non-resonant contribution cleanly.","section":"Table I"}],"recommendation":"major_revision","confidential_remarks":"The central interpretational claim needs to be reworked: Section III B attributes the MicroBooNE discrepancy to inelastic strength while Section IV concedes that this cannot explain the discrepancy. This inconsistency, together with the missing statistical details, makes the manuscript not ready in its current form. The paper is a phenomenological comparison rather than a new formalism, and its contribution would be strengthened by a more careful, channel-by-channel, kinematics-resolved statement of where the model succeeds and where it fails."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThis is a straightforward extension of the SuSAv2 model to the NOvA and MicroBooNE CC-inclusive data. It is not conceptually new—the framework, the DCC resonance input, and the Bodek-Ritchie DIS structure functions are all prior work by the same group—but it is a legitimate and honest validation step. The good news: the comparison is internally consistent, the smearing-matrix procedure for MicroBooNE is documented, and the paper gives chi-square values for the main comparisons. I found no circularity: no parameters are fitted to the NOvA or MicroBooNE data under comparison.\n\nThe soft spots are in the interpretation, not the mechanics. The repeated claim that underestimation \"can be ascribed to some missing strength in the inelastic channels\" does not hold up against the paper's own channel decomposition. For MicroBooNE, the flux peaks at 0.8 GeV, QE dominates, and the total cross section is underestimated—missing inelastic strength cannot fill that gap. The conclusions section itself concedes this, saying the inelastic lack is \"not expected to explain these discrepancies.\" For NOvA muon neutrino, the model overpredicts at high muon energies, attributed to DIS—the opposite sign of a missing-strength problem. So the central causal attribution is shaky. What remains is a legitimate documentation of where SuSAv2 works and where it does not.\n\nThe other issue is quantitative: there are no theoretical uncertainty bands on the curves, and the quoted chi-squares lack degrees of freedom or covariance details. That makes \"reproduce well\" hard to audit, and the comparison across models is suggestive rather than rigorous.\n\nNone of this is fatal. The paper is a useful data point for people building generators or comparing nuclear models against these measurements. It deserves a serious referee, but the referee should push for a revision that either drops or substantially qualifies the inelastic attribution, adds some error estimate on the theory side, and reports chi-square properly.\n\nWho is this for: neutrino-nucleus modelers, oscillation-experiment MC developers, and anyone tracking how well SuSAv2 transfers to argon. I would cite it if I worked in that space.\n\nRecommendation: send to peer review with revisions expected.","headline":"A careful SuSAv2 comparison with NOvA and MicroBooNE data, but the central attribution of discrepancies to missing inelastic strength is not supported by the paper's own channel decomposition.","tokens_in":14638,"tokens_out":2550,"would_cite":true,"duration_ms":23489,"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 claims that the SuSAv2 framework, which factorizes neutrino-nucleus scattering into single-nucleon responses and nuclear scaling functions, reproduces the shape and value of NOvA and MicroBooNE charged-current inclusive cross…","keywords":["neutrino-nucleus scattering","charged-current inclusive cross sections","SuSAv2 superscaling model","quasielastic and inelastic channels","meson exchange currents","deep inelastic scattering","NOvA","MicroBooNE"],"falsifier":"Take a high-statistics inclusive electron-argon scattering measurement at $Q^2$ near 1 GeV$^2$/c$^2$ and invariant mass $W$ between 1.4 and 2.1 GeV and compare it with the SuSAv2 inelastic channel sum: if the same forward-angle, high-energy deficit appears, the missing strength lies in the nuclear scaling functions rather than in the neutrino-specific inelastic models.","tokens_in":13688,"feed_emoji":"⚛️","tokens_out":8519,"duration_ms":79555,"temperature":0.7,"pith_summary":"This paper claims that the SuSAv2 model, built on superscaling and relativistic mean field theory, can describe the charged-current inclusive neutrino cross sections recently measured by NOvA (carbon-rich target) and MicroBooNE (argon) across neutrino energies from tens of MeV to about 20 GeV. The model splits the cross section into five channels—quasielastic, two-particle two-hole meson-exchange currents, resonances, soft deep inelastic, and true deep inelastic—and the comparison locates a systematic shortfall at very forward angles and high lepton energies, which the authors ascribe to missing strength in the inelastic channels. If the picture is right, the same factorization can be used to improve event generators in the few-GeV regime, and the deficit gives a concrete place to look for missing nuclear or single-nucleon inelastic contributions.","feed_headline":"Five-channel neutrino model reproduces NOvA and MicroBooNE data","feed_subtitle":"SuSAv2's channel decomposition fits inclusive cross sections; the shortfall sits in forward, high-energy inelastic channels.","key_machinery":"The central object is the SuSAv2 scaling-function factorization: each reaction channel is written as a single-nucleon response multiplied by a nuclear scaling function $f(\\psi)$ extracted from relativistic mean field calculations and electron-scattering data, with the inelastic response integrated over the reduced invariant mass $\\mu_X$ using a generalized scaling variable $\\psi_X$ and single-nucleon inelastic structure functions. The resonance channel uses the Dynamical Coupled-Channels (DCC) model from the Osaka group, the DIS part uses the Bodek-Ritchie parameterization split into SoftDIS (below $W_X = 2.1$ GeV) and TrueDIS (above), and MEC comes from relativistic Fermi gas calculations. This machinery lets each discrepancy be assigned to a named channel, which is what carries the argument that the deficit is inelastic rather than quasielastic.","core_discovery":"The paper's central claim is that the SuSAv2 channel decomposition—QE, MEC, RES, SoftDIS, TrueDIS—reproduces the shape and normalization of the recent NOvA and MicroBooNE charged-current inclusive measurements well enough to expose a pattern: underestimation at very forward angles and high lepton energies, which the authors attribute to missing strength in the inelastic channels. For NOvA electron neutrinos the overall agreement is good and the model's $\\chi^2$ is below that of NuWro, GiBUU and GENIE; for NOvA muon neutrinos the model underestimates the peak region and overpredicts at high muon energies, with the shape-only $\\chi^2$ factor of 3 tracing to the DIS description. For MicroBooNE, the model underpredicts the 2022 total and single-differential cross sections on argon, while the 2024 three-dimensional data agree for neutrino energies below 1.6 GeV and are underpredicted above, where inelastic channels matter most. The authors conclude that the discrepancies point to the inelastic description rather than to the QE or MEC pieces, and that new inelastic ingredients are needed.","pith_inferences":["The paper's attribution of the MicroBooNE deficit to inelastic channels is not fully tested by its own comparison, because at average energies below 1 GeV the QE channel dominates; a CC0π measurement on argon at the same flux would isolate whether the QE scaling function itself transfers to argon.","One way to sharpen the missing-inelastic-strength claim would be to repeat the comparison with electron scattering on argon: the same forward-angle deficit in $(e,e')$ data would point to the nuclear scaling functions rather than neutrino-specific structure functions.","The negative cross-section bins produced by the MicroBooNE smearing matrix mean part of the apparent discrepancy could be an unfolding artifact; a forward-folded comparison without the regularization would test this."],"forward_implications":["If the SuSAv2 channel decomposition is right, event generators that adopt it will inherit a concrete high-energy shortfall: forward-angle, high-lepton-energy inclusive events will be underpredicted until the inelastic single-nucleon structure functions are augmented.","The NOvA electron-neutrino agreement implies the factorization works over a broad energy range on a mixed carbon/chlorine target, supporting its use in near-detector oscillation analyses.","The MicroBooNE pattern implies that below 1.6 GeV the QE-dominated model is adequate on argon, so future argon measurements should focus on the inelastic threshold region to find the missing strength.","The shape-only $\\chi^2$ factor of 3 for NOvA muon neutrinos, if confirmed, means the normalization can be tuned by 2p2h adjustments but the angular shape requires changes in the DIS treatment."],"supporting_citations":[{"why":"Provides the NOvA electron-neutrino CC-inclusive double-differential data that the model is compared with.","marker":"[26]"},{"why":"Provides the NOvA muon-neutrino CC-inclusive data where the model underestimates the peak and shows a shape-only chi-squared factor of 3.","marker":"[58]"},{"why":"Provides the MicroBooNE argon total and single-differential cross sections that the model underpredicts.","marker":"[56]"},{"why":"Provides the MicroBooNE three-dimensional differential data that agree below 1.6 GeV and show the deficit at higher energies.","marker":"[57]"},{"why":"Establishes the SuSAv2-DCC treatment of resonances and DIS against electron, T2K, and MINERvA data, the baseline this work extends.","marker":"[38]"},{"why":"Supplies the Dynamical Coupled-Channels single-nucleon resonance structure functions used for the RES channel.","marker":"[39–41]"},{"why":"Supplies the relativistic Fermi gas two-particle two-hole MEC calculation used for the MEC channel.","marker":"[24,25]"},{"why":"Supplies the RMF-based scaling functions that define the SuSAv2 quasielastic and inelastic responses.","marker":"[19,20]"}],"fun_headline_variants":["SuSAv2 fits neutrino data, but inelastic channels lag","Inclusive neutrino data: SuSAv2 good except forward inelastic","SuSAv2 model exposes inelastic shortfall in neutrino cross sections","NOvA/MicroBooNE: SuSAv2 undercounts forward inelastic events"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The factorization is assumed to hold across the entire energy spectrum, from quasielastic to deep inelastic scattering, and for argon as well as carbon; if the scaling function extracted from quasielastic electron-carbon data does not transfer to argon or to high-energy inelastic kinematics, the five-channel decomposition loses its basis.","fun_headline_variants_meta":{"raw":{"variants":["SuSAv2 fits neutrino data, but inelastic channels lag","Inclusive neutrino data: SuSAv2 good except forward inelastic","SuSAv2 model exposes inelastic shortfall in neutrino cross sections","NOvA/MicroBooNE: SuSAv2 undercounts forward inelastic events"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000495,"raw_usage":{"total_tokens":2404,"prompt_tokens":897,"completion_tokens":1507,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":513,"completion_tokens_details":{"reasoning_tokens":1423}},"tokens_in":513,"tokens_out":1507,"duration_ms":11835,"temperature":1.0,"reasoning_tokens":1423,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T05:32:15.349581+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a high-statistics inclusive electron-argon scattering measurement at $Q^2$ near 1 GeV$^2$/c$^2$ and invariant mass $W$ between 1.4 and 2.1 GeV and compare it with the SuSAv2 inelastic channel sum: if the same forward-angle, high-energy deficit appears, the missing strength lies in the nuclear scaling functions rather than in the neutrino-specific inelastic models.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the NOvA electron-neutrino CC-inclusive double-differential data that the model is compared with."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the NOvA muon-neutrino CC-inclusive data where the model underestimates the peak and shows a shape-only chi-squared factor of 3."},{"cited_title":"Measurement of three-dimensional inclu- sive muon-neutrino charged-current cross sections on argon with the microboone detector,","cited_arxiv_id":null,"evidence_quote":"Provides the MicroBooNE three-dimensional differential data that agree below 1.6 GeV and show the deficit at higher energies."},{"cited_title":"Gonzalez-Rosa, G","cited_arxiv_id":null,"evidence_quote":"Establishes the SuSAv2-DCC treatment of resonances and DIS against electron, T2K, and MINERvA data, the baseline this work extends."}],"review_version":1}