{"id":"b1c1dabd-49cf-4ebd-aca9-30e9ca264f04","arxiv_id":"1908.08603","paper_version":1,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"low","formal_verification":"none","parameter_count":3,"one_line_summary":"A topical review of neutral-current quasielastic neutrino-nucleus scattering, comparing relativistic models and highlighting the role of strange form factors and multi-nucleon effects.","lead":"Neutral-current neutrino-nucleus scattering is a key source of systematic uncertainty for oscillation experiments, and this paper surveys the theoretical models used to describe it. It explains how nuclear effects and the nucleon's strange-quark content enter the cross sections, and how current models compare with data from MiniBooNE and BNL.","discovery_kind":"review","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The review's central 'impossibility' claim for a single axial mass in CCQE and NCE rests on MiniBooNE comparisons that mix true and reconstructed Q2, different efficiency corrections, and mostly omit non-QE backgrounds; the claim may be analysis-dependent.","rationale":"The reader's weakest assumption concerned selection bias and correct reproduction of the numerical examples. My concern is more specific but related: the data comparisons used to draw the main conclusion are not all subjected to a common detector-response, efficiency, and background treatment. This matters because the MiniBooNE NCE measurement is made in terms of the total reconstructed nucleon kinetic energy T, and Q2_QE is inferred via equation (34). Without consistent folding, the fitted axial mass and the comparison among models can shift enough to change the qualitative conclusion. The review is otherwise careful and comprehensive: it includes many caveats, acknowledges the RFG exception and the role of 2p-2h and MEC contributions, and points to the NuWro result where a standard MA becomes acceptable when np-nh and backgrounds are included. The central claim is therefore presented with more nuance than the reader's strongest_claim suggests, but the evidence for 'general impossibility' is still not fully convincing. A conditional acceptance is appropriate: the article should explicitly state that the conclusion depends on the choice of observable definition, the inclusion of backgrounds, and the detector-response treatment, and it should recommend that future comparisons use a common reconstruction-based framework.","tokens_in":49230,"tokens_out":4267,"duration_ms":45634,"concrete_test":"Take the spectral-function model used for figure 20 (Ref. [69]), compute its double-differential cross section, fold it with the MiniBooNE flux, apply the published efficiency correction functions of equation (32), and use the migration matrices from the MiniBooNE NCE data release (Ref. [220]) to convert true nucleon kinetic energies to reconstructed T and Q2_QE. Then add the irreducible pion-production background in the same way as the NuWro analysis (Ref. [227]) and refit MA to the MiniBooNE NCE differential cross section. If a single value of MA around 1.3 GeV can describe both the fitted CCQE and NCE data within this same model, the Section 6 'impossibility' claim is not supported; if the inconsistency persists after this common treatment, the claim is corroborated.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The review's principal conclusion in Section 6 is that 'it is generally impossible to describe, within the same model, both CCQE and NCE data using the same value of MA'. This conclusion is supported mainly by comparing model calculations with MiniBooNE NCE data in Section 5.2 (e.g., figures 20-34). The load-bearing condition is that these comparisons are apples-to-apples in their treatment of the measured observable. They are not. Equation (32) defines the CH2 cross section with efficiency correction functions Ci, and equation (34) defines Q2_QE = 2mT, where T is the total reconstructed kinetic energy of all detected nucleons. Some curves are plotted versus true Q2 while others are plotted versus reconstructed Q2_QE or T, and only some calculations include detector smearing via migration matrices, as in Ref. [67] and the NuWro analysis of Ref. [227]. The MiniBooNE T-based measurement sums over all final-state nucleons, whereas many model curves are single-nucleon knockout predictions without the pion-production background and without multi-nucleon final states. The review itself notes in Section 6 that 'it is not quite clear how the multinucleon component shows up in the experimental data.' A concrete illustration of the fragility is given by the NuWro analysis in figure 32: once np-nh contributions and backgrounds are included, the fitted MA becomes consistent with the world average (MA=1.10 GeV), in contrast to the larger values extracted by simpler IA models. Thus, the claimed incompatibility of CCQE and NCE axial masses may be an artifact of inconsistent observable definitions and missing background/model components, not a robust physical conclusion.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"Giusti and Ivanov review the theoretical description of neutral-current (NC) quasielastic neutrino-nucleus scattering in the kinematic regime of accelerator neutrino experiments. They present the general NC/CC formalism and the single-nucleon weak current (Section 2), then survey the main model families: relativistic impulse approximation and optical-potential approaches, spectral-function models, superscaling (SuSA), RPA and np-nh treatments, and meson-exchange currents (Section 3). They discuss the strange nucleon form factors and proposed ratios/asymmetries for extracting them (Section 4), and compare model predictions with BNL E734 and MiniBooNE neutral-current elastic data, paying attention to flux averaging and to the difference between true and reconstructed Q2 (Section 5). The concluding remarks emphasize that one-body IA models without multi-nucleon channels underpredict data unless an enhanced axial mass is used, that no single axial mass describes both CCQE and NCE data within such models, and that consistent treatment of multi-nucleon and two-body current contributions is needed for DUNE and Hyper-K.","tokens_in":49515,"tokens_out":9778,"duration_ms":97515,"significance":"The manuscript is a comprehensive and generally balanced topical review. Its strengths are the systematic organization of a large and technically diverse literature, the explicit caveats attached to model families (notably the warning that the RGF treatment is designed for inclusive observables and may include channels absent in semi-inclusive NCE measurements), the reproducible presentation of the flux-averaging equations, and the extensive reference list. It does not present new derivations or machine-checked results; its conclusions are literature-based judgments, but that is appropriate for the review format. If the axial-mass conclusion is qualified as proposed below, the review will serve as a useful reference for both specialists and non-specialists and as a basis for systematic studies in DUNE and Hyper-K.","major_comments":[{"comment":"The concluding statement that 'it is generally impossible to describe, within the same model, both CCQE and NCE data using the same value of MA' is stronger than the evidence assembled in Section 5.2 and should be qualified. The comparisons underlying this conclusion are not uniform in the measured observable: Eq. (32) introduces per-channel efficiency functions C_i for the CH2 cross section, Eq. (34) defines the reconstructed variable Q2_QE = 2mT with T the sum of all final-state nucleon kinetic energies, and the figures mix true and reconstructed quantities (e.g., Fig. 24 shows true Q2 in the left panel and reconstructed energy in the right panel), with detector smearing applied only in some calculations (Refs. [67], [227]). In addition, the NuWro analysis in Fig. 32, which includes np-nh contributions and detector response, yields MA = 1.10^{+0.13}_{-0.15} GeV, consistent with the deuterium world average; this is a same-model counterexample if the claim is not restricted to IA one-nucleon-knockout models. I recommend rewriting the bullet to say explicitly that the impossibility claim applies to one-body-current IA models without multi-nucleon/background channels, and to acknowledge that a model-independent no-go conclusion is not established because no single analysis treats all models with identical unfolding and efficiency corrections.","section":"Section 6, concluding remarks (fourth bullet)"},{"comment":"The comparison figures would be much easier to interpret if each panel explicitly stated whether the plotted quantity is true Q2, reconstructed Q2_QE = 2mT (Eq. (34)), or nucleon kinetic energy T, and whether efficiency corrections and migration matrices were applied. The current text mixes these definitions: Fig. 20 uses Eq. (32) with efficiency functions C_i, Fig. 22 uses Q2_QE, Fig. 24 shows true energy (left) and reconstructed energy (right), and Fig. 28 compares several models without a unified folding prescription. As a result, the reader cannot determine how much of the spread between model and data is due to physics rather than to the treatment of detector response. A short summary table or an explicit sentence per figure would materially strengthen the evidentiary basis of the Section 6 conclusions.","section":"Section 5.2, Eq. (32) and Figs. 20-28"}],"minor_comments":[{"comment":"Equation (32) as printed contains a duplicated term: the carbon-neutron contribution 3/7 C_{nu n,C} dsigma_{nu n,C}/dQ2 appears twice, and the bound-proton contribution is absent. The intended decomposition for CH2 should contain one 3/7 term for carbon protons and one 3/7 term for carbon neutrons.","section":"Section 5.2, Eq. (32)"},{"comment":"The sentence listing the strange form factors reads 'F s 1 , F s 1 , and G s A'; the second entry should be F s 2 (the strange vector magnetic form factor).","section":"Section 2.5, text before Eq. (20)"},{"comment":"The phrase 'An example in shown in figure 4' contains a typo; 'in' should be 'is'.","section":"Section 3.3, p. 14"},{"comment":"The phrase 'nuclei wih A > 4' contains a typo; 'wih' should be 'with'.","section":"Section 3.6, p. 26"},{"comment":"The caption's error-bar sentence 'in the nu(nu) case' should read 'in the nu(bar nu) case'; the same typo appears in the text describing the antineutrino panel.","section":"Section 5.2, Fig. 28 caption"},{"comment":"The word 'understimanted' should be 'underestimated'.","section":"Section 5.2, p. 46"}],"recommendation":"major_revision","confidential_remarks":"To the editor: the manuscript is a competent and useful topical review that fits the journal's scope. The only substantive concern is the overbroad axial-mass conclusion; I believe a revision that qualifies the claim and clarifies the observable definitions will resolve it. I do not see a need for new calculations or additional data comparison. The self-citation level is high but appropriate for a review in which the authors' groups have made central contributions to the subject."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Carlotta - quick take on Giusti and Ivanov's NC review. This is a solid, genuinely useful review. It covers the NC QE model landscape more thoroughly than most existing reviews, which focus on CC, and it does a good job explaining the difference between inclusive CC and semi-inclusive NC observables. The discussion of strange form factors and the ratios proposed to isolate them is clear and practical. The paper is honest about model limitations, including the RGF caveat (it may include channels absent from semi-inclusive NCE) and the unresolved question of how multinucleon strength appears in the T-based MiniBooNE observable.\n\nThe main soft spot is the concluding claim that it is 'generally impossible' to describe CCQE and NCE data with the same axial mass. That statement is too categorical. The MiniBooNE comparisons it rests on mix true and reconstructed Q2, apply different efficiency corrections, and most of the quoted model curves are one-nucleon knockout without np-nh or pion backgrounds. More importantly, the review itself cites the NuWro analysis in figure 32, where adding np-nh brings MA back to 1.10 GeV, consistent with the world average. So the 'impossibility' applies to simple IA models, not to models with two-body currents. The authors do mention these caveats in places, but the bullet in Section 6 doesn't integrate them, and a reader could easily take away a stronger conclusion than the evidence supports.\n\nSecondary issue: many of the numerical examples come from the authors' own papers. That's not circular, since these are review examples, but the selection could bias the impression of which models work best. Independent comparisons would strengthen it.\n\nThat said, this is a review, not a new result, and on those terms it is accurate and fair. The field needs this kind of map, and it deserves a serious referee. I'd suggest accepting it, with a request that the authors soften the axial-mass conclusion to explicitly distinguish IA-only models from models with multinucleon contributions, and acknowledge the observable-definition dependence more prominently.\n\nFor your reading group: worth a look for people working in neutrino cross sections, though it's long. I'd cite it if I needed a reference for the NC model landscape. Serious thinker: yes.","headline":"Useful, balanced NC QE scattering review, but the 'impossible axial mass' conclusion is more conditional than the abstract/conclusions suggest once you factor in the NuWro np-nh fit and MiniBooNE's reconstructed-Q2 caveats.","tokens_in":50081,"tokens_out":3200,"would_cite":true,"duration_ms":34123,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The review concludes that nuclear effects, not a larger axial mass, explain neutrino-nucleus data, and that multi-nucleon and two-body current contributions must be included explicitly.","keywords":["neutrino-nucleus scattering","neutral current","quasielastic scattering","axial mass","strange form factors","two-body currents","final-state interactions","nuclear effects"],"falsifier":"A single high-statistics experiment measuring CCQE and NCE on the same nucleus with a well understood flux, analyzed with one model that includes two-body currents and the world-average $M_A$, would settle it: success would confirm the review's diagnosis, and failure would show the missing physics is elsewhere.","tokens_in":48989,"feed_emoji":"⚛️","tokens_out":6965,"duration_ms":72638,"temperature":0.7,"pith_summary":"This review synthesizes the theory of neutral-current quasielastic neutrino-nucleus scattering in the hundreds-of-MeV to few-GeV range, the regime relevant for next-generation long-baseline experiments. It argues that nuclear effects, not the nucleon's axial form factor, are the leading systematic uncertainty, and it identifies a central tension: models built on the impulse approximation with one-nucleon knockout underpredict both charged-current and neutral-current elastic cross sections unless the axial mass $M_A$ is artificially raised. The review's key conclusion is that no single value of $M_A$ can describe charged-current and neutral-current data within the same one-body model, so multi-nucleon excitations and two-body meson-exchange currents must be included explicitly. It also shows that ratios of cross sections, especially proton-to-neutron and neutrino-to-antineutrino ratios, suppress nuclear-model uncertainties and are the practical route to the strange axial form factor. If this picture is right, future oscillation analyses cannot rely on tuning $M_A$; they need consistent reaction models covering all contributing channels.","feed_headline":"Neutrino data need multi-nucleon effects, not a bigger axial mass","feed_subtitle":"A theory review shows charged- and neutral-current scattering can only be reconciled with consistent two-body currents.","key_machinery":"The machinery is the contraction of the leptonic tensor with a hadronic tensor $W^{\\mu\\nu}$, built from the single-nucleon weak current whose axial form factor is a dipole in $Q^2$ with mass $M_A$ and whose isoscalar component can carry strange form factors. Around this current, the review compares competing descriptions of the nucleus: relativistic impulse approximation with plane or distorted waves, spectral functions with short-range correlations, the superscaling approach built from electron-scattering data, RPA treatments with np-nh excitations, and Green's-function formulations of final-state interactions. The load-bearing observables are flux-averaged cross sections and, above all, ratios (proton-to-neutron, neutrino-to-antineutrino, NC-to-CC) that cancel nuclear-model uncertainties.","core_discovery":"The review's central claim is that the apparent discrepancy between the axial mass extracted from deuterium data ($M_A\\simeq 1.03$ GeV) and the larger values needed to fit carbon-target data is not a property of the nucleon but a symptom of missing nuclear physics. One-nucleon knockout models, even with sophisticated final-state interactions, underestimate the measured neutral-current elastic cross sections; models that add collective RPA correlations, multi-nucleon (np-nh) excitations, and two-body meson-exchange currents can describe the same data with the standard axial mass. The review further claims that charged-current and neutral-current data cannot both be fit with a single $M_A$ inside one-body models, and that cross-section ratios are the observables that best isolate the strange axial form factor of the nucleon.","pith_inferences":["If the axial-mass discrepancy is caused by missing two-body currents, then effective $M_A$ values extracted from CCQE and NCE should converge once the same 2p-2h model is used; this convergence is a quantitative test that future combined fits can perform.","The review's RGF results suggest final-state inelasticity can mimic two-body current strength; a model that tracks explicit final-state channels could separate the two and remove a possible double-counting.","For argon targets, building the spectral function from a dedicated $(e,e'p)$ measurement would test whether short-range correlations in argon differ enough from carbon to change NC neutrino event rates in next-generation detectors."],"forward_implications":["Oscillation analyses must treat the 'quasielastic' sample as a mix of one-nucleon and multi-nucleon events; otherwise the reconstructed neutrino energy is biased.","Adding the same two-body current to both CCQE and NCE predictions should remove the need for separate effective axial masses, making the two channels consistent.","Cross-section ratios such as $R(p/n)$ and the $\\nu$-$\\bar\\nu$ asymmetry become primary tools for strange form factors because nuclear effects largely cancel.","Models that pass the electron-scattering test are necessary but not sufficient; neutrino data impose additional constraints from the axial channel.","New detectors with better final-state resolution will constrain hadronic observables and flux, reducing the model dependence highlighted here."],"supporting_citations":[{"why":"Supplies the standard deuterium-based axial mass $M_A=1.026\\pm0.021$ GeV that the review treats as the baseline.","marker":"[39]"},{"why":"Supplies the charged-current data whose fit motivated the larger effective axial mass and the controversy.","marker":"[40]"},{"why":"Shows that a spectral-function model fitting the CCQE data with $M_A=1.6$ GeV fails to reproduce NCE data, motivating the inconsistency claim.","marker":"[68]"},{"why":"Compares spectral-function predictions with NCE $Q^2$ data to demonstrate that no single axial mass describes both channels.","marker":"[69]"},{"why":"Gives RMF and SuSA fits to NCE data, extracting effective axial masses and strangeness values used in the comparison.","marker":"[125]"},{"why":"Provides RPA-plus-multinucleon predictions that reproduce NCE data with the standard axial mass, supporting the multi-nucleon explanation.","marker":"[164]"},{"why":"Provides the neutrino NCE dataset used in most model-data comparisons and in the strange axial extraction.","marker":"[207]"},{"why":"Provides the antineutrino NCE dataset used for the $\\bar\\nu/\\nu$ ratio comparison.","marker":"[213]"},{"why":"Simultaneously fits $M_A$ and $\\Delta s$ in a Monte Carlo generator and finds the standard $M_A$ once np-nh contributions are included.","marker":"[227]"}],"fun_headline_variants":["Nuclear effects, not axial mass, explain neutrino data","Multi-nucleon effects rescue standard axial mass","Theorists: neutrino scattering needs two-body currents","Neutrino puzzles solved by nuclear correlations","Beyond one-body: nuclear theory fixes neutrino data"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The conclusions assume that the selected numerical examples and data comparisons, many from the authors' prior work, are representative of the broader literature and are folded with correct flux normalizations and detector efficiencies; if the neutral-current sample or the model inputs are not representative, the claimed axial-mass discrepancy could be an artifact.","fun_headline_variants_meta":{"raw":{"variants":["Nuclear effects, not axial mass, explain neutrino data","Multi-nucleon effects rescue standard axial mass","Theorists: neutrino scattering needs two-body currents","Neutrino puzzles solved by nuclear correlations","Beyond one-body: nuclear theory fixes neutrino data"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000642,"raw_usage":{"total_tokens":2929,"prompt_tokens":894,"completion_tokens":2035,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":510,"completion_tokens_details":{"reasoning_tokens":1963}},"tokens_in":510,"tokens_out":2035,"duration_ms":12991,"temperature":1.0,"reasoning_tokens":1963,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:34:16.393404+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A single high-statistics experiment measuring CCQE and NCE on the same nucleus with a well understood flux, analyzed with one model that includes two-body currents and the world-average $M_A$, would settle it: success would confirm the review's diagnosis, and failure would show the missing physics is elsewhere.","supporting_citations":[{"cited_title":"Martini, M","cited_arxiv_id":null,"evidence_quote":"Provides RPA-plus-multinucleon predictions that reproduce NCE data with the standard axial mass, supporting the multi-nucleon explanation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the neutrino NCE dataset used in most model-data comparisons and in the strange axial extraction."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the antineutrino NCE dataset used for the $\\bar\\nu/\\nu$ ratio comparison."},{"cited_title":"Golan, K","cited_arxiv_id":null,"evidence_quote":"Simultaneously fits $M_A$ and $\\Delta s$ in a Monte Carlo generator and finds the standard $M_A$ once np-nh contributions are included."}],"review_version":1}