Pith. sign in

REVIEW 4 major objections 5 minor 13 references

Charge-exchange reactions with pion and kaon beams in the NA64h experiment at CERN

T0 review · 4 major / 5 minor · reviewed 2026-07-31 · grok-4.5

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

Pith's one-line read Regge fits give charge-exchange cross sections for light neutral mesons to about 10 percent accuracy from 5 to 200 GeV.

desk verdict Usable Regge fits and nuclear tables for NA64h CEX yields, with the nuclear extrapolation carrying more weight than the paper admits. read the letter →

arxiv 2607.28299 v1 pith:BUNNJSML submitted 2026-07-30 hep-ph

Sergei N. Gninenko , Sergey Kuleshov , Valery E. Lyubovitskij , Alexey S. Zhevlakov This is my paper · ORCID
classification hep-ph
keywords charge-exchangereactionsReggeformalismneutralmesonproductionNA64hdarksectorpionbeamskaonnuclearcrosssections
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

The paper supplies a practical set of integral cross-section formulas for producing π⁰, η, η′ and K̄⁰ in charge-exchange reactions of pion and kaon beams on protons and heavier nuclei. These yields are the missing ingredient needed to forecast the sensitivity of dark-sector searches that look for invisible decays of the same mesons at the NA64h experiment. A two-parameter Regge form, fitted once per channel and then scaled by a mildly corrected nuclear factor, reproduces existing data across the full energy range planned for the experiment. The resulting tables and curves therefore become a ready data set for both NA64h projections and any other measurement that relies on the same production mechanism.

What carries the argument

The two-parameter Regge formula σ(s) = σ₀ (s/s₀)^{2α−2} (s₀ = 10 GeV²) together with the nuclear scaling factor Z^{2/3 − 0.15/Z^{2/3}}. The formula converts a handful of measured points into continuous predictions for every beam momentum and every nucleus of interest.

What would settle it

A new high-statistics measurement of any of the four charge-exchange cross sections on a heavy nucleus (for example tungsten or lead) at a beam momentum between 5 and 15 GeV that lies systematically outside the 10 percent error band of the published formulas.

Watch

Extended reading notes

Core claim

Integral charge-exchange cross sections for π⁻p → π⁰n, π⁻p → η(η′)n and K⁻p → K̄⁰n are reproduced to roughly 10 percent by the simple Regge expression σ(s) = σ₀ (s/s₀)^{2α−2} with the four quoted pairs (σ₀, α). Extension from the proton to an arbitrary nucleus is achieved by multiplying by the corrected power Z^{2/3 − 0.15/Z^{2/3}}, which matches the 40 GeV nuclear data and is then used for all other energies and targets.

Load-bearing premise

A single energy-independent intercept α per channel plus the ad-hoc nuclear correction remain valid for every meson, every target from hydrogen to lead, and every momentum down to 5 GeV.

Editorial extensions

If this is right

  • NA64h can convert the tabulated yields directly into projected limits on invisible η, η′ and K̄⁰ decays without further model dependence on production.
  • The same formulas supply the expected meson rates for any other fixed-target experiment that uses pion or kaon beams in the 5–50 GeV window.
  • The corrected Z-power can be tested or refined once new nuclear data at lower energies become available.
  • Kaon-beam charge exchange is placed on the same quantitative footing as the better-measured pion channels.

Reading between the lines

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

  • If the 10 percent accuracy holds at the lowest planned energies, the dominant uncertainty in dark-sector reach will shift from production yield to detector acceptance and background rejection.
  • The mild deviation from pure Z^{2/3} scaling hints that absorption or surface effects are already visible at 40 GeV and may grow at lower momenta, offering a clean target for a dedicated nuclear scan.
  • Once the cross sections are fixed, the same data set can be reused to normalize searches for visible rare decays or for oscillations into sterile states that share the same production vertex.
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

4 major / 5 minor

Summary. The manuscript provides a unified phenomenological description of integral charge-exchange (CEX) cross sections for π⁻p → π⁰n, π⁻p → η(η′)n and K⁻p → K̄⁰n, and their extension to nuclear targets, using the two-parameter Regge form σ(s)=σ₀(s/s₀)^{2α−2} with s₀=10 GeV². Parameters σ₀ and α are fitted channel-by-channel to existing hydrogen data over P_Lab ≈ 5–200 GeV (Eqs. 2–4, Tables I–IV); nuclear targets are obtained by the empirical factor Z^{2/3−0.15/Z^{2/3}} (Eqs. 5–6), tuned to 40 GeV π⁰/η data (Tables V–VI). The stated accuracy is ~10%. The results are intended as input for meson-yield estimates in the NA64h dark-sector program at CERN (P_Lab = 5–50 GeV) and for other experiments.

Significance. If the quoted cross sections and their uncertainties are reliable across the full NA64h energy and target matrix, the work supplies a compact, immediately usable data set for sensitivity projections of invisible neutral-meson decays and oscillations. The hydrogen fits themselves are a useful consolidation of scattered older measurements under a single Regge parametrization, with point-by-point residuals that are mostly at the few-to-fifteen-percent level. The nuclear extension and the low-s edge (s∼s₀) are the parts that directly control projected yields; their robustness therefore determines how much the paper advances the experimental program beyond existing compilations.

major comments (4)
  1. [Section II, Eqs. (5)–(6), Tables V–VI] Eqs. (5)–(6) and Tables V–VI: the nuclear factor Z^{2/3−0.15/Z^{2/3}} is fixed solely by the 40 GeV π⁰ and η data set and is then applied unchanged to η′, K̄⁰, every other energy (including the NA64h lower edge P_Lab=5 GeV), and every nucleus from Li to Pb. No independent cross-check, no A-dependent absorption, and no energy dependence of the correction are provided. Because meson yields (and therefore dark-sector sensitivity) scale directly with this factor, the uncontrolled systematic must be quantified or the domain of validity clearly restricted.
  2. [Section II, Eqs. (2)–(3), Tables I–IV] Section II and Figs. 1–4: the same energy-independent intercepts α are used from s=s₀=10 GeV² (P_Lab∼5 GeV) up to several hundred GeV². At the lower edge the Regge asymptotic form is only marginally justified; several low-P_Lab points already show residuals of order 10% (Table IV, P_Lab=5 GeV, R=−10.9%). A quantitative statement of the additional systematic assigned to the s∼s₀ region, or a comparison with alternative low-energy parametrizations, is needed before the formulas can be used for NA64h running at 5–15 GeV.
  3. [Eq. (2), Tables I–VI, Conclusions] The claimed global accuracy “≃10%” is not uniformly supported. In Table II the relative residual R reaches 12–16% at the highest energies; Table IV shows several entries near or above 10%; the nuclear data in Tables V–VI carry experimental uncertainties of 10–30%. The error bands quoted on σ₀ (Eq. 2) appear to be fit uncertainties only. A transparent propagation that folds in data scatter, the nuclear-factor uncertainty, and the low-s extrapolation is required if the numbers are to be used for sensitivity estimates.
  4. [Section II] No information is given on the fit procedure itself (χ² definition, weighting, whether α and σ₀ were fitted simultaneously or sequentially, covariance between parameters). Without this, the quoted ±10% errors and the “predictions” at 50 GeV (Eq. 4) cannot be independently assessed or propagated.
minor comments (5)
  1. [Introduction and Eq. (5)] Inconsistent notation for the nuclear factor: the Introduction writes “Z^{2/3}−0.15Z^{−2/3}” (additive), while Eq. (5) writes Z^{2/3−0.15/Z^{2/3}} (modified exponent). The tables are consistent only with the exponent form; the text must be aligned.
  2. [Throughout] Typos and grammar: “at nuclei target”, “for P_Lab running”, “CER reactions”, “π− +p→η[→2γ]+n” spacing, “10 Gev”, “Fit(eff. Reggeon)” in figure labels, repeated “H, Li” in the nuclei list of Figs. 5–8.
  3. [Figs. 5–8] Figures 5–8 show only the central curves; the stated 10% error bands are mentioned in the caption text but not drawn, making visual assessment of the nuclear extrapolation difficult.
  4. [Tables II–III, Eqs. (2)–(4)] The paper never states whether the η and η′ cross sections are corrected for the 2γ branching fractions or are the observed 2γ rates; the notation σ^{η;2γ} suggests the latter, but this should be explicit for users who need total meson yields.
  5. [References] Reference list is adequate for the data but omits more recent global Regge/CEX analyses that could serve as external benchmarks (beyond the single JPAC citation).

Circularity Check

3 steps flagged · score 4.0 of 10

Regge σ₀, α and the nuclear exponent are fitted to the same CEX data the paper then tabulates as ‘theory’/‘predictions’; honest phenomenology, not a closed definitional loop.

  1. fitted input called prediction [Sec. II, Eqs. (1)–(4) and Tables I–IV]
    "The parameters σ₀=σ(s₀) and α are fitted from data. ... The final expressions for the fitted integral cross sections including errors ∼10% read [Eq. 2]. ... Our predictions for the PLab=50 GeV/s are σπ−H=(8.09±0.38)µb, σηH=(1.04±0.07)µb, ..."

    σ₀ and α are fixed by the same hydrogen energy scans that the tables then report as theory-vs-data agreement. The 50 GeV ‘predictions’ are direct evaluations of those fitted functions inside the fitted PLab window (data exist from 5–200 GeV), so they are statistically forced by the fit rather than independent forecasts.

  2. fitted input called prediction [Sec. II, Eqs. (5)–(6) and Tables V–VI]
    "Extension to arbitrary nuclei is to multiply with factor Z^{2/3−0.15/Z^{2/3}}. ... we modified the power of the nuclei charge Z as Z^{2/3}→Z^{2/3−0.15/Z^{2/3}} to describe data for heavier nuclei."

    The −0.15/Z^{2/3} correction is chosen expressly so that σ_Z matches the 40 GeV nuclear π⁰ and η measurements. Tables V–VI then present that same tuned formula against those measurements as validation; the reported agreement is by construction of the exponent adjustment.

1 more flagged steps
  1. self citation load bearing [Sec. II, paragraph on nuclear extension; Ref. [2]]
    "Note in Ref. [2] we found that the above formula is slightly modified. In particular, we modified the power of the nuclei charge Z as Z^{2/3}→Z^{2/3−0.15/Z^{2/3}} to describe data for heavier nuclei."

    The non-standard nuclear factor that carries the multi-target results is justified only by the authors’ prior paper [2] (overlapping author list), not by an independent external derivation. It is not the sole support for the hydrogen fits, so the step is secondary rather than fully load-bearing for the whole claim.

full rationale

The paper is a two-parameter Regge phenomenology plus an empirical nuclear rescaling. σ₀ and α in Eqs. (2)–(3) are explicitly fitted to the hydrogen CEX data that Tables I–IV then compare to ‘theory’; the PLab=50 GeV numbers labeled ‘predictions’ (Eq. 4) are simply the same fitted curves evaluated inside the fitted energy range. The nuclear factor Z^{2/3−0.15/Z^{2/3}} (Eq. 5) is introduced solely ‘to describe data for heavier nuclei’ and is checked only against the 40 GeV π⁰/η points of Tables V–VI that motivated the −0.15 correction (via self-citation to the authors’ Ref. [2]). That is ordinary curve-fitting presented with mild predictive language, not a self-definitional identity: the inputs are external measured cross sections, the functional form is the standard Regge ansatz, and extrapolations to η′, K̄⁰, other energies, and other nuclei are genuine uses of the fit rather than tautologies. No uniqueness theorem or load-bearing external ‘derivation’ is smuggled in. Score 4 reflects partial fitted-input-called-prediction and non-load-bearing self-citation without elevating a data-driven parametrization to circular derivation.

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

The central numerical claims rest entirely on the Regge power-law ansatz, four pairs of fitted (σ₀, α), a fixed s₀=10 GeV², and one empirically adjusted nuclear exponent. No new dynamical entity is introduced; everything else is standard high-energy phenomenology or published cross-section data.

free parameters (4)
  • σ₀(π⁰), σ₀(η), σ₀(η′), σ₀(K̄⁰) = 109±13, 29±2, 1.05±0.10, 145±25 μb
    Normalization of each channel’s integral cross section at s=s₀; fixed by fit to proton data.
  • α_π, α_η, α_η′, α_K = 0.405, 0.260, 0.200, 0.248
    Effective Regge intercepts controlling the energy slope; fixed by fit to the same proton data.
  • nuclear exponent correction 0.15 = 0.15
    Ad-hoc shift of the naive Z^{2/3} scaling, tuned so that the 40 GeV nuclear data on heavier targets are reproduced.
  • s₀ = 10 GeV²
    Reference scale in the Regge factor (s/s₀)^{2α−2}; chosen by hand as 10 GeV² to mark the start of the fitted regime.
assumptions (3)
  • domain assumption Integral CEX cross sections at large s and small |t| obey the single-trajectory Regge form σ(s)=σ₀(s/s₀)^{2α−2} with energy-independent α.
    Stated in Sec. II as the working framework; standard but not derived here, and applied down to s≈s₀.
  • ad hoc to paper Nuclear A-dependence factors as Z^{2/3−0.15/Z^{2/3}} independent of meson species and beam energy.
    Introduced in Sec. II to correct pure Z^{2/3}; justified only by 40 GeV π⁰/η data on a handful of nuclei.
  • domain assumption Published integral cross-section data sets (Refs. [8–10]) are accurate enough that a global ~10% fit residual is meaningful.
    All fit quality claims rest on taking those measurements at face value without re-assessment of systematics.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Charge-exchange reactions with pion and kaon beams in the NA64h experiment at CERN." pith.science (2026). https://pith.science/paper/BUNNJSML

@misc{pith2026260728299,
  author       = {Pith},
  title        = {Pith review of: Charge-exchange reactions with pion and kaon beams in the NA64h experiment at CERN},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/BUNNJSML}},
  note         = {Machine review of arXiv:2607.28299}
}
abstract

We discuss a new approach to search for dark sector coupled to quarks in the invisible decays and oscillations of light neutral mesons produced in the charge-exchange (CEX) reactions with pion and kaon beams $\pi^-(K^-)+ (A, Z) \to \pi^0, \eta, \eta' (\bar K^0) + (A,Z-1)$ at nuclei target $(A,Z)$ in the NA64h experiment at CERN for $P_{\rm Lab}$ running in the interval from 5 to 50 GeV. For estimating of the projection sensitivity for the proposed searches the knowledge of the meson yield is crucial. This work is dedicated to the accurate evaluation of CEX cross sections for the wide range of energies and target nuclei which could be also useful for other experiments.

Figures

Figures reproduced from arXiv: 2607.28299 by the authors.

Figure 1
Figure 1. FIG. 1: Fit of the integral cross section of the [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2: Fit of the integral cross section of the [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3: Fit of the integral cross section of the [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figures from the paper (5 more)
Figure 4
Figure 4. Figure 4: FIG. 4: Fit of the integral cross section of the [PITH_FULL_IMAGE:figures/full_fig_p005_4.png]
Figure 5
Figure 5. Figure 5: FIG. 5: Results for the integral cross section of the [PITH_FULL_IMAGE:figures/full_fig_p006_5.png]
Figure 6
Figure 6. Figure 6: FIG. 6: Results for the integral cross section of the [PITH_FULL_IMAGE:figures/full_fig_p006_6.png]
Figure 7
Figure 7. Figure 7: FIG. 7: Results for the integral cross section of the [PITH_FULL_IMAGE:figures/full_fig_p007_7.png]
Figure 8
Figure 8. Figure 8: FIG. 8: Results for the integral cross section of the [PITH_FULL_IMAGE:figures/full_fig_p007_8.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

13 extracted references · 1 linked inside Pith

  1. [1]

    Yu. M. Andreevet al.(NA64 Collaboration), Phys. Rev. Lett.133, 121803 (2024)

    2024

  2. [2]

    S. N. Gninenko, D. V. Kirpichnikov, S. Kuleshov, V. E. Lyubovitskij, and A. S. Zhevlakov, Phys. Rev. D109, 075021 (2024)

    2024

  3. [3]

    S. N. Gninenko, Phys. Rev. D91, 015004 (2015)

    2015

  4. [4]

    S. N. Gninenko and N. V. Krasnikov, Phys. Rev. D92, 034009 (2015); Mod. Phys. Lett. A31, 1650142 (2016)

    2015

  5. [5]

    S. N. Gninenko, D. V. Kirpichnikov, N. V. Krasnikov, S. Kuleshov, V. E. Lyubovitskij, and A. S. Zhevlakov, JHEP06, 235 (2025)

    2025

  6. [6]

    S. N. Gninenko, N. V. Krasnikov and V. A. Matveev, Natural Sci. Rev.1, 5 (2024),

    2024

  7. [7]

    S. N. Gninenko and N. V. Krasnikov, Phys. Rev. D111, 11 (2025)

    2025

  8. [8]

    A. V. Barneset al., Phys. Rev. Lett.37, 76 (1976); O. I. Dahlet al., Phys. Rev. Lett.37, 80 (1976); W. D. Apelet al.(Serpukhov-CERN Collaboration), Phys. Lett. B83, 131 (1979); Nucl. Phys. B152, 1 (1979); Nucl. Phys. B154, 189 (1979); V. I. Belousovet al., Phys. Lett. B43, 76 (1973); V. D. Apokinet al., Sov. J. Nucl. Phys.46, 644 (1987); S. V. Donskovet a...

    1976

Show all 13 references
  1. [9]

    V. D. Apokinet al., Sov. J. Nucl. Phys.35, 219 (1982)

    1982

  2. [10]

    Astburyet al., Phys

    P. Astburyet al., Phys. Lett.23, 396 (1966); G. W. Brandenburget al., Phys. Rev. D15, 617 (1977); W. D. Apelet al. (Serpukhov-CERN Collaboration), Nucl. Phys. B129, 275 (1977); F. G. Binonet al.(Serpukhov-Brussels-Annecy(LAPP) and Serpukhov-CERN Collaborations), Sov. J. Nucl. ...

    1966

  3. [11]

    I. I. Strakovsky, EPJ Web Conf.130, 06005 (2016); M. Amaryanet al.(KLF Collaboration), [arXiv:2008.08215 [nucl-ex]]; M. J. Amaryanet al.(KLF Collaboration), Mod. Phys. Lett. A39, 2450063 (2024)

    2016 arXiv

  4. [12]

    Regge, Nuovo Cim.14, 951 (1959); V

    T. Regge, Nuovo Cim.14, 951 (1959); V. N. Gribov, Sov. Phys. JETP26, 414(1968); P. D. B. Collins,An Introduction to Regge Theory and High-Energy Physics, Cambridge Univ. Press, 2009

    1959

  5. [13]

    Nyset al.(JPAC Collaboration), Phys

    J. Nyset al.(JPAC Collaboration), Phys. Rev. D98, 034020 (2018). 5 Fit(eff. Reggeon) 10 20 30 40 50 0 20 40 60 80 100 120 PLab (GeV) σ[p π-  n π0] (μ b) Fit(eff. Reggeon) 50 100 150 200 0 20 40 60 80 100 120 PLab (GeV) σ[p π-  n π0] (μ b) FIG. 1: Fit of the integral cross se...

    2018

Pith tools

Reviewed July 31, 2026 · model on record in the stance chip above.