Pith. sign in

REVIEW 3 major objections 5 minor 13 references

Quasielastic Lepton-Nucleus Scattering and the Correlated Fermi Gas Model

T0 review · 3 major / 5 minor · reviewed 2026-08-12 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read 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.

desk verdict 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. read the letter →

arxiv 2412.19810 v2 pith:7TJ5Y3CM submitted 2024-11-27 hep-ph

Sam Carey This is my paper · ORCID
classification hep-ph
keywords CorrelatedFermiGasquasielasticscatteringneutrino-nucleuselectronnuclearformfactorsshort-rangecorrelationsRelativistic
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

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.

What carries the argument

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.

What would settle it

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.

Watch

Extended reading notes

Core claim

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.

Load-bearing premise

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.

Editorial extensions

If this is right

  • 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.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • 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.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

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.

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 (3)
  1. [Section 2, Eq. (2)] 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.
  2. [Section 2, Eq. (2), blocking factor] 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.
  3. [Section 3, Figs. 2-5] 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.
minor comments (5)
  1. [Abstract and Introduction] 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.
  2. [Figure 3] 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.
  3. [Figure 5 caption] 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.
  4. [References] 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.
  5. [Figures 2-5] 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.

Circularity Check

1 steps flagged · score 4.0 of 10

Circularity is limited to the load-bearing self-citation of companion [1] for the analytic CFG formulas; the data comparisons are benchmark checks, not fits.

  1. self citation load bearing [Section 2, after Eq. (2); Abstract; Conclusion]
    "In the CFG model, the nucleons are distributed in regions I, II initially and post-interaction can avail the phase space in regions III, IV and V as depicted in the right side of Figure 1. The total cross section is then calculated by summing the contributions from the six possible transitions. More detailed description of the implementation is discussed in [1]."

    The paper's stated central claim is 'a fully analytic implementation of the Correlated Fermi Gas nuclear model,' but the only concrete content shown is the generic convolution Eq. (2) and a list of six transitions. The actual analytic formulas, the explicit CFG distributions n_i(p) and n_f(p+q), the normalization, the phase-space boundaries, and the λ,c0 values used for the carbon curves are all deferred to [1], a companion paper with overlapping authorship. The implementation and the subsequent 'better agreement' conclusion therefore cannot be derived or checked from this manuscript; they rest entirely on an unshown self-citation. This is a load-bearing citation of the authors' own prior work rather than an independent external result.

full rationale

I did not find a derivation that reduces to its inputs in the sense of Eq. X = Eq. Y by construction, nor any parameter fitted in this paper and then relabeled as a prediction. The CFG model itself is an external published model [2], and the electron-carbon and MiniBooNE data [12,13] are external benchmarks. The paper performs fixed-model comparisons without fitting λ or c0, so the 'better agreement' statement is a consistency/benchmark check rather than an independent prediction. The main circularity concern is structural: the paper's central claim (analytic CFG implementation) is supported only by the sentence 'More detailed description of the implementation is discussed in [1]', where [1] is a companion paper by Bhattacharya, Carey, Cohen, and Paz and thus shares this paper's author. Because the formulas are not displayed, the reader cannot verify that the six transitions exhaust the phase space or that the blocking factor 1−n_f(p+q) is implemented as claimed; correctness collapses if [1] contains an error. I count this as one load-bearing self-citation rather than full circularity, because the comparison figures and the external data/model inputs are present and are not derived from the paper's own fit. A χ² or parameter-uncertainty statement would strengthen the paper, but its absence is a rigor issue, not circularity.

Assumptions & free parameters 4 free parameters · 6 assumptions · 0 invented entities

The paper introduces no new parameters or entities; all model parameters (CFG λ, c0, Fermi momentum, form factor coefficients) are taken from prior literature. The central claim depends on the analytic formulas in the companion paper [1], which are not reproduced here. No invented particles or forces appear.

free parameters (4)
  • CFG parameters λ and c0 = λ ∈ [2.50, 3.00], c0 ∈ [3.21, 5.11] (extremal variation in Fig. 1; central values not stated)
    These parameters define the depleted Fermi sea and high-momentum tail of the CFG model. They come from the CFG literature (Hen et al. [2]) and are not derived or fit in this paper.
  • Fermi momentum pF for carbon = ~225 MeV/c (standard value; not stated in paper)
    Sets the Fermi gas scale in both RFG and the base of CFG. It is taken as an input from nuclear physics literature.
  • Dipole axial mass MA = 1.014 ± 0.014 GeV
    Used in the BBBA dipole axial form factor for neutrino scattering [3]; fitted to neutrino-deuteron data in that reference.
  • Z-expansion form factor coefficients = Multiple coefficients from [4-11]
    Used for BHLT, MBGH, NME22, Mainz22, PNDME, ETMC form factors; all fitted to data or lattice QCD in the cited works.
assumptions (6)
  • domain assumption Impulse approximation: the cross-section factorizes into lepton tensor and nuclear tensor with single-nucleon knockout
    Eq. (1)-(2) assume the lepton interacts with one nucleon and the rest of the nucleus is a spectator.
  • standard math Energy-momentum conservation at the nucleon level is enforced by an on-shell delta function
    Eq. (2) includes δ(ϵp - ϵ'_{p+q} + q0).
  • domain assumption The CFG momentum distribution is fully described by a depleted Fermi sea plus a high-momentum tail with parameters λ and c0
    Section 2 and Fig. 1 right; this is the defining model from [2], not derived here.
  • domain assumption The six initial-to-final phase-space transitions exhaust the CFG phase space
    Section 2 says the total cross section is the sum of six transitions, with no proof of completeness.
  • domain assumption External datasets (Barreau [12], MiniBooNE [13]) are reliable and have been properly corrected for backgrounds
    Section 3 compares directly to these data without re-analysis.
  • domain assumption The cited form factor parametrizations [3-11] are valid at the kinematics considered
    All form factors are taken as external inputs.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Quasielastic Lepton-Nucleus Scattering and the Correlated Fermi Gas Model." pith.science (2026). https://pith.science/paper/7TJ5Y3CM

@misc{pith2026241219810,
  author       = {Pith},
  title        = {Pith review of: Quasielastic Lepton-Nucleus Scattering and the Correlated Fermi Gas Model},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/7TJ5Y3CM}},
  note         = {Machine review of arXiv:2412.19810}
}
read the original abstract

The study of neutrino-nucleus scattering processes is important for the new generation neutrino experiments for better understanding of the neutrino oscillation phenomenon. A significant source of uncertainty in the cross-section comes from limitations in our knowledge of nucleon and nuclear effects in the scattering process. Here we present a fully analytic implementation of the Correlated Fermi Gas nuclear model for the quasi-elastic lepton nucleus scattering. The implementation is then used to separately compare form factor effects and nuclear model effects, specifically the Relativistic Fermi Gas nuclear model, for both electron-carbon and neutrino-carbon scattering data.

Figures

Figures reproduced from arXiv: 2412.19810 by the authors.

Figure 1
Figure 1. Left: the comparison between the different axial form factor parameterizations: Dipole (BBBA) using MA = 1.014 ± 0.014 GeV [3], MBGH [5], RQCD 20 [6], NME 22 [7], Mainz 22 [8], MINERvA [9], PNDME 23 [10] and ETMC 23 [11] Right: The distribution of nucleons in the CFG model with the different lines corresponding to the extremal variation of λ and c0. 3 Results First, we compare cross-section between the two nuclear m… view at source ↗
Figure 2
Figure 2. Left: A comparison of the RFG (red) and CFG (blue) nuclear model with the Carbon scattering data. Right: The contribution to the CFG implementation from various transitions. The incoming electron energy is 480 MeV and scattering off at 60◦ . In [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. A comparison of the BHLT (teal) and BBBA (purple) vector form factor with the Carbon scattering [PITH_FULL_IMAGE:figures/full_fig_p003_3.png] view at source ↗
Figures from the paper (2 more)
Figure 4
Figure 4. Figure 4: A comparison of RFG (red) and CFG (blue) models of flux averaged neutrino - carbon scattering [PITH_FULL_IMAGE:figures/full_fig_p004_4.png]
Figure 5
Figure 5. Figure 5: A comparison of BBBA, MBGH, NME22, Mainz22, MINER [PITH_FULL_IMAGE:figures/full_fig_p004_5.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

13 extracted references · 13 canonical work pages

  1. [1]

    Quasielastic Lepton-Nucleus Scattering and the Correlated Fermi Gas Model

    B. Bhattacharya, S. Carey, E. O. Cohen, and G. Paz. arXiv: 2405.05342 [hep-ph]

    work page Pith review arXiv

  2. [2]

    O. Hen, B. A Li, W. J. Guo, L. B. Weinstein, and E. Piasetzky. Phys. Rev. C, 91(2): 025803 (2015)

    work page 2015

  3. [3]

    Bradford, A

    R. Bradford, A. Bodek, H. S. Budd, and J. Arrington. Nucl. Phys. B Proc. Suppl., 159: 127–132 (2006)

    work page 2006

  4. [4]

    Borah, R.J

    K. Borah, R.J. Hill, G. Lee, and O. Tomalak. Phys. Rev. D, 102(7): 074012 (2020)

    work page 2020

  5. [5]

    Meyer, M

    A.S. Meyer, M. Betancourt, R. Gran, and R.J. Hill. Phys. Rev. D, 93(11): 113015 (2016)

    work page 2016

  6. [6]

    Bali et al

    G.S. Bali et al. JHEP, 05: 126 (2020)

    work page 2020

  7. [7]

    Park et al

    S. Park et al. Phys. Rev. D, 105(5): 054505 (2022)

    work page 2022

  8. [8]

    Djukanovic et al

    D. Djukanovic et al. Phys. Rev. D, 106(7): 074503 (2022)

    work page 2022

Show all 13 references
  1. [9]

    Cai et al

    T. Cai et al. Nature, 614(7946): 48–53 (2023)

    2023

  2. [10]

    Jang et al

    Y-C. Jang et al. Phys. Rev. D, 109(1): 014503 (2024)

    2024

  3. [11]

    Alexandrou et al

    C. Alexandrou et al. Phys. Rev. D, 109(3): 034503 (2024)

    2024

  4. [12]

    Barreau et al

    P. Barreau et al. Nucl. Phys. A, 402(3): 515–540 (1983)

    1983

  5. [13]

    A. A. Aguilar-Arevalo et al. Phys. Rev. D, 81: 092005 (2010). 5

    2010

Pith tools

Reviewed August 12, 2026 · model on record in the stance chip above.