{"id":"3d03d912-925b-4104-9417-347e487fb73e","arxiv_id":"2412.19810","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"An analytic Correlated Fermi Gas implementation is compared with the Relativistic Fermi Gas model, showing that flux-averaged neutrino-carbon cross sections barely distinguish the two nuclear models but respond visibly to the axial form factor.","lead":"This paper presents an analytic version of the Correlated Fermi Gas model for quasielastic lepton-nucleus scattering and compares it with the standard Relativistic Fermi Gas model on carbon data. The comparison suggests that flux-averaged neutrino cross sections are more sensitive to the nucleon form factor than to the nuclear model, which is useful for neutrino oscillation experiments.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central claim is unverifiable here: the analytic CFG formulas, final-state Pauli-blocking choice, and the λ,c0 values used for carbon are all deferred to companion [1], so the six-transition enumeration and the 'better agreement' conclusion cannot be checked from this paper.","rationale":"The reader's weakest assumption matches my reading. The paper's novelty and central claim is the analytic CFG implementation, yet Section 2 gives no explicit CFG distribution, no transition integrals, and no normalization; it explicitly defers to [1]. Therefore the strongest claim cannot be checked from this manuscript alone. The most specific technical soft spot is the final-state phase space: Eq. (2) includes the Pauli-blocking factor [1−n_f(p+q)], but n_f is never defined. If n_f is taken to be the same CFG distribution, transitions into region II are not fully blocked and should contribute, making the six listed transitions incomplete. If n_f is instead a step function, that is an unstated modeling choice that changes the physics. Either way, the comparison is not reproducible from the text. The 'better agreement' conclusion is also underquantified: no χ², no parameter uncertainties, and no stated values of λ and c0 for the CFG curves, so even if the formulas are correct, the secondary claim is suggestive rather than established. I do not see an internal contradiction that forces rejection; the companion paper may resolve every issue. The appropriate response is therefore to keep the conditional verdict, requiring the authors to make the explicit formulas and quantitative comparisons available.","tokens_in":5802,"tokens_out":10074,"duration_ms":100196,"concrete_test":"Fix the Fig. 2 left kinematics (E_i=480 MeV, θ=60°, carbon) and independently evaluate Eq. (2) by numerical quadrature using the CFG momentum distribution of Ref. [2], with the specific λ and c0 values stated in companion [1]. Compare this direct integral to the sum of the displayed six analytic transitions. If the sums disagree, or if transitions ending in region II contribute non-negligibly when 1−n_f(p+q)>0, the six-transition enumeration is incomplete and the analytic implementation is not validated.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 2 displays only the generic convolution formula, Eq. (2), and says that the total CFG cross section is the sum of six transitions (I→III, I→IV, I→V, II→III, II→IV, II→V), with the actual implementation 'discussed in [1]'. No explicit form of n_i(p), n_f(p+q), the normalization N, the region boundaries, or the λ and c0 values used for the carbon curves is given. The unavoidable consequence is that the paper's stated central claim—a fully analytic CFG implementation—cannot be derived or reproduced from the manuscript alone. The most concrete technical ambiguity is the final-state blocking factor 1−n_f(p+q) in Eq. (2): if n_f is the same CFG distribution, transitions into region II (I→II and II→II) are kinematically allowed and the six listed transitions may not exhaust the phase space; if n_f is instead a step function with full blocking below some cutoff, that choice is not stated. Either way, the figures cannot be regenerated. The secondary claim that CFG shows better agreement with electron-carbon data is made by visual inspection of Fig. 2 left, with no χ², no parameter uncertainties, and no statement of which λ,c0 pair produced the curve. The correctness of the central claim therefore rests entirely on an external companion that this paper does not expose.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports an analytic implementation of the Correlated Fermi Gas (CFG) nuclear model for quasielastic lepton-nucleus scattering, building on a companion paper [1]. It compares CFG with the Relativistic Fermi Gas (RFG) model for electron-carbon and neutrino-carbon data, using several vector and axial form-factor parametrizations. The stated central claim is that the implementation is fully analytic, and the secondary claim is that CFG shows better agreement with electron-carbon data for specific kinematics. The comparisons are made with fixed nuclear model or fixed form factor, and the paper concludes that nuclear-model differences are visible in electron scattering but partially washed out in flux-averaged neutrino scattering.","tokens_in":6124,"tokens_out":2919,"duration_ms":27079,"significance":"If the analytic CFG implementation is correct and reproducible, it would provide a fast closed-form nuclear model that incorporates short-range correlations, useful for benchmarking and for separating nuclear-model uncertainties from form-factor uncertainties in quasielastic lepton-nucleus scattering. The paper makes a useful qualitative comparison of available form-factor parametrizations and of CFG versus RFG, and it does not fit any parameters, instead using published inputs. However, the significance is currently limited because the central derivation is not present in this manuscript: the reader cannot verify or reproduce the analytic implementation from the text and figures, and the claimed better agreement is not quantified.","major_comments":[{"comment":"The central claim of a fully analytic implementation is not supported within this manuscript. Equation (2) gives the generic convolution formula, but the CFG-specific ingredients are not shown: the explicit forms of n_i(p) and n_f(p+q), the normalization N, the phase-space boundaries for the six transitions, and the values of lambda and c0 used for the carbon curves are all deferred to the companion paper [1]. The text only states that the total cross section is the sum of six transitions and that 'More detailed description of the implementation is discussed in [1].' Consequently, the analytic expressions cannot be derived or checked from the present paper, and the six-transition enumeration as the complete CFG phase space remains an assertion rather than a demonstrated result. Please include the explicit analytic formulas and parameter values, either in the text or in an appendix.","section":"Section 2, Eq. (2)"},{"comment":"The final-state Pauli blocking factor [1 - n_f(p+q)] is ambiguous for the CFG model. If n_f is the same CFG distribution that includes region II, then transitions into region II (I to II and II to II) are kinematically allowed, and the six listed transitions (I to III, I to IV, I to V, II to III, II to IV, II to V) would not exhaust the phase space. If n_f is instead a step function that fully blocks final states below a certain momentum, that choice must be stated explicitly, since it changes the kinematics and the claimed analytic result. Please specify the precise form of n_f used in Eq. (2) and justify why the six listed transitions are complete.","section":"Section 2, Eq. (2), blocking factor"},{"comment":"The secondary claim that 'The CFG model shows better agreement with data for specific kinematics in electron scattering' is made by visual inspection only. No error bars or uncertainty bands are shown for the model curves, no chi-squared or other goodness-of-fit quantity is reported, and the text does not state which value of (lambda, c0) produced the CFG curves in each figure. Because the CFG parameters were not fitted in this work but are taken from the companion paper, the reader cannot assess whether the agreement is statistically meaningful or merely a consequence of the chosen parameter values. Please provide quantitative comparison metrics and state the model parameters used in each figure.","section":"Section 3, Figs. 2-5"}],"minor_comments":[{"comment":"The phrasing 'Here we present a fully analytic implementation' overstates what this paper alone delivers, since the detailed implementation is in the companion paper [1]. Please rephrase to make clear that the present contribution summarizes and applies that implementation.","section":"Abstract and Introduction"},{"comment":"The right panel of Figure 3 has the vertical axis labeled 'd2σ/dΩdω (103 b/sr-GeV)', while the left panel uses '103 nb/sr-GeV'. The units in the right panel appear to be a typo, as the numerical magnitudes are the same; please correct to nb/sr-GeV or clarify.","section":"Figure 3"},{"comment":"The figure caption lists 'BBBA' among axial form factor parametrizations, but BBBA is primarily a vector form factor parametrization. Please clarify how the BBBA set is used in the neutrino-scattering calculation, or rename the label to indicate the full set of form factors used.","section":"Figure 5 caption"},{"comment":"Reference [1] is cited only by arXiv number. Provide the title and, if available, the journal reference or a DOI, so that readers can locate the companion paper and the specific equations used for the CFG implementation.","section":"References"},{"comment":"The Q2 axis in Figures 2 and 3 appears as a secondary axis with a label 'Q2 (10-3 GeV2)', but the relationship between omega and Q2 is not explained in the text or captions. Please clarify the kinematic mapping or remove the secondary axis if it is redundant.","section":"Figures 2-5"}],"recommendation":"major_revision","confidential_remarks":"This is a proceedings-style contribution whose central technical content lives in the companion paper [1]. For a journal publication, the manuscript is not self-contained: the analytic formulas, the CFG phase-space derivation, and the parameter values used are not disclosed, and the data comparisons lack quantitative metrics. The main results are therefore not verifiable from this paper alone. If the companion paper is accepted and the authors add the necessary technical details, the work could become a useful concise application note; in its current form, it reads as an extended abstract. The editor may also wish to consider whether a proceedings contribution with this level of technical detail meets the journal's standards for original research."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Sam's NuFact proceedings. The short version: the analytic CFG implementation is presented by reference — the derivation sits in the authors' companion [1], and this paper runs the comparisons. If you want the closed-form CFG formulas, they are not here. Read this for the benchmarks.\n\nWhat is genuinely useful: the clean separation of nuclear-model and form-factor effects. The MiniBooNE comparison makes a point worth keeping — after flux averaging, RFG and CFG are nearly indistinguishable, while the choice of axial form factor gives a visible spread. That is the right practical message for neutrino experiments, and it is stated without overclaiming. At fixed electron kinematics, the CFG-vs-RFG comparison nicely shows the redistribution of strength from the quasielastic peak into the tail, and the transition decomposition in Fig. 2 is pedagogically clear. The axial form-factor list is current and relevant, and the citation practice is honest. The electron comparison is a benchmark rather than a prediction, since the CFG parameters came from electron data, and the paper does not pretend otherwise.\n\nThe soft spots, in proportion. The stress-test concern is real: Eq. (2) has the factor 1 − n_f(p+q), and the paper never says what n_f is. If n_f is the CFG distribution, final states in the tail are only partially blocked and the six listed transitions do not obviously exhaust the phase space — I→II and II→II would contribute. If n_f is a step function with full blocking below p_F, that is a choice that should be stated. Either way the six-transition enumeration cannot be checked from this text and the figures cannot be regenerated. Pointing to [1] is legitimate for a proceedings, but the ambiguity is cheap to resolve in-text, and a referee should ask for it. Second, \"better agreement\" with the Barreau data is visual: no chi-squared, no error bars, and no statement of which (lambda, c0) pair produced the curves in Figs. 2-5. Since the claim is hedged to \"specific kinematics,\" I call that a minor gap, not a fatal one. Third, small typos: Fig. 3 right says \"b\" where it should be \"nb,\" and Fig. 5's caption has \"MINER νA.\"\n\nBottom line: thin as a standalone paper, but the benchmarks are useful, the message is sensible, and I see no load-bearing error. The derivation belongs in [1]; this paper deserves a serious referee who asks for the blocking factor, the parameter choices, and a quantitative agreement measure — then it is fine as a proceedings.","headline":"Useful benchmark comparisons in a proceedings paper whose core CFG derivation lives in the companion [1]; the blocking-factor ambiguity flagged in the stress test is real, and the referee should pin it down along with a quantitative agreement measure.","tokens_in":6630,"tokens_out":6615,"would_cite":false,"duration_ms":59103,"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":"This paper presents a fully analytic implementation of the Correlated Fermi Gas nuclear model for quasielastic lepton-nucleus scattering, and uses it to separate nuclear-model from form-factor uncertainties.","keywords":["Correlated Fermi Gas","quasielastic scattering","neutrino-nucleus scattering","electron scattering","nuclear form factors","short-range correlations","Relativistic Fermi Gas"],"falsifier":"An independent numerical Monte Carlo integration of the nuclear tensor with the same CFG momentum distribution, the same model parameters, and the same BBBA or BHLT form factors should reproduce the dashed CFG curves for 480 MeV electrons at 60 degrees; any residual difference beyond the numerical error would show that the analytic implementation, or the six-transition decomposition, is missing part of the phase space.","tokens_in":5606,"feed_emoji":"⚛️","tokens_out":5229,"duration_ms":44198,"temperature":0.7,"pith_summary":"This paper presents a fully analytic implementation of the Correlated Fermi Gas (CFG) nuclear model for quasielastic electron- and neutrino-nucleus scattering. The implementation decomposes the cross section into six transitions between depleted Fermi-sea and high-momentum tail regions, replacing the standard Relativistic Fermi Gas (RFG) phase space. Using it, the paper separates form-factor effects from nuclear-model effects against electron-carbon and neutrino-carbon data. The central finding is that CFG redistributes cross-section strength from the quasielastic peak into the tail and agrees better with electron-carbon data than RFG for the kinematics shown, while flux-averaged neutrino cross sections change little with the nuclear model.","feed_headline":"Fully analytic Fermi gas model sharpens quasielastic cross sections","feed_subtitle":"Six-transition implementation separates nuclear-model from form-factor effects in electron and neutrino scattering.","key_machinery":"The machinery is the six-transition phase-space decomposition of the CFG model: initial nucleons live in the depleted Fermi-sea region I and high-momentum tail region II, and after the interaction can populate final regions III, IV and V, so the total cross section is the sum of I to III, I to IV, I to V, II to III, II to IV and II to V. The nuclear tensor carries the physics with the occupation factors $n_i(p)[1-n_f(p+q)]$ and the energy-conserving delta function, and the hadron tensor is built from the chosen vector (BBBA or BHLT) or axial form factors. The parameters $\\lambda$ and $c_0$ set the tail depletion, and a companion paper contains the detailed formulas.","core_discovery":"The central claim is that a fully analytic CFG implementation exists and works: cross sections from the nuclear-tensor expression are obtained by summing six closed-form transitions among the five momentum regions defined in the paper. This yields concrete comparisons: at electron energy 480 MeV and scattering angle 60 degrees, the CFG curves shift strength away from the quasielastic peak and improve agreement with the electron-carbon data relative to RFG; for flux-averaged neutrino-carbon scattering, the nuclear-model dependence almost disappears, while the choice of axial form factor produces a visible spread. The paper's intended lesson is that CFG is a viable fast analytic benchmark that separates the two main systematic uncertainties in lepton-nucleus quasielastic scattering.","pith_inferences":["A natural next test is to apply the same six-transition implementation to oxygen or argon targets; if the model parameters can be determined for those nuclei, the analytic form makes target-systematics scans nearly cost-free.","The near-invisibility of CFG versus RFG in flux-averaged neutrino data implies that differential observables, such as recoil nucleon momentum or energy-transfer bins, will be required to pin down correlation effects, an implication the paper hints at with its semi-inclusive outlook.","Because the axial-form-factor spread is larger than the nuclear-model spread at the studied neutrino kinematics, oscillation analyses that rely on this cross section may need to marginalize over axial form factors before claiming sensitivity to nuclear models."],"forward_implications":["CFG curves can be produced in closed form, so quasielastic cross-section estimates no longer require per-event numerical sampling of the correlated momentum distribution.","In electron-carbon scattering, moving strength from the peak to the tail is the feature that brings CFG closer to data, so tail-sensitive kinematics are where nuclear-model effects show up.","For flux-averaged neutrino-carbon measurements, nuclear-model differences are too small to resolve, and axial form-factor parametrization is the dominant visible spread.","Because the two effects are separated with one model fixed at a time, comparisons of this kind give a benchmark for assigning which systematic uncertainty dominates in a given observable."],"supporting_citations":[{"why":"Companion paper that supplies the detailed analytic six-transition formulas on which this implementation rests.","marker":"[1]"},{"why":"Original CFG model definition with the depleted Fermi sea and high-momentum tail regions.","marker":"[2]"},{"why":"Provides the BBBA dipole parametrization used for the vector and axial form factors.","marker":"[3]"},{"why":"Provides the BHLT z-expansion vector form factor used for the electron-scattering comparisons.","marker":"[4]"},{"why":"Provides the MBGH axial form factor used in the neutrino-scattering comparisons.","marker":"[5]"},{"why":"Provides the Mainz 22 axial form factor used in the neutrino-scattering comparisons.","marker":"[8]"},{"why":"Electron-carbon scattering data used as the benchmark that CFG is compared against.","marker":"[12]"},{"why":"Flux-averaged neutrino-carbon scattering data used to compare nuclear models and axial form factors.","marker":"[13]"}],"fun_headline_variants":["Analytic CFG model separates nuclear and form-factor effects","Fully analytic Fermi gas model isolates scattering uncertainties","Analytic CFG sharpens quasielastic lepton-nucleus calculations","New analytic CFG benchmark decouples nuclear and form-factor uncertainties"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the analytic six-transition formulas, which this paper does not display but attributes to the companion paper, are correct and that those six transitions exhaust the CFG phase space; if either fails, the plotted CFG curves are not what the paper claims them to be.","fun_headline_variants_meta":{"raw":{"variants":["Analytic CFG model separates nuclear and form-factor effects","Fully analytic Fermi gas model isolates scattering uncertainties","Analytic CFG sharpens quasielastic lepton-nucleus calculations","New analytic CFG benchmark decouples nuclear and form-factor uncertainties"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00067,"raw_usage":{"total_tokens":2968,"prompt_tokens":771,"completion_tokens":2197,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":387,"completion_tokens_details":{"reasoning_tokens":2126}},"tokens_in":387,"tokens_out":2197,"duration_ms":14616,"temperature":1.0,"reasoning_tokens":2126,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T11:33:48.807557+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"An independent numerical Monte Carlo integration of the nuclear tensor with the same CFG momentum distribution, the same model parameters, and the same BBBA or BHLT form factors should reproduce the dashed CFG curves for 480 MeV electrons at 60 degrees; any residual difference beyond the numerical error would show that the analytic implementation, or the six-transition decomposition, is missing part of the phase space.","supporting_citations":[{"cited_title":"Quasielastic Lepton-Nucleus Scattering and the Correlated Fermi Gas Model","cited_arxiv_id":"2405.05342","evidence_quote":"Companion paper that supplies the detailed analytic six-transition formulas on which this implementation rests."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Original CFG model definition with the depleted Fermi sea and high-momentum tail regions."},{"cited_title":"Bradford, A","cited_arxiv_id":null,"evidence_quote":"Provides the BBBA dipole parametrization used for the vector and axial form factors."},{"cited_title":"Borah, R.J","cited_arxiv_id":null,"evidence_quote":"Provides the BHLT z-expansion vector form factor used for the electron-scattering comparisons."},{"cited_title":"Meyer, M","cited_arxiv_id":null,"evidence_quote":"Provides the MBGH axial form factor used in the neutrino-scattering comparisons."},{"cited_title":"Djukanovic et al","cited_arxiv_id":null,"evidence_quote":"Provides the Mainz 22 axial form factor used in the neutrino-scattering comparisons."},{"cited_title":"Barreau et al","cited_arxiv_id":null,"evidence_quote":"Electron-carbon scattering data used as the benchmark that CFG is compared against."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Flux-averaged neutrino-carbon scattering data used to compare nuclear models and axial form factors."}],"review_version":1}