Pith. sign in

REVIEW 2 major objections 4 minor 28 references

Accelerator neutrinos can measure the low-density potentials that decide when hyperons appear inside neutron stars, with a forecast precision of about 6 MeV on the Lambda potential.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

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

T0 review · grok-4.5

2026-07-13 04:21 UTC pith:4DZ2KN3D

load-bearing objection Clean, candid forecast that neutrino FSI can pin low-density U_Λ to ~6 MeV; U_Σ is systematics-limited and M_max is not claimed from the data. the 2 major comments →

arxiv 2607.09255 v1 pith:4DZ2KN3D submitted 2026-07-10 nucl-th astro-ph.HEhep-ph

Accelerator neutrinos as a probe of in-medium hyperon potentials

Jaroslaw Nowak This is my paper
classification nucl-th astro-ph.HEhep-ph
keywords accelerator neutrinoshyperon potentialsfinal-state interactionsneutron starshyperon puzzleliquid argonSBNDDUNE
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

Neutron stars above two solar masses require a stiff equation of state, yet hyperons that form at high density soften it and threaten to undercut the observed masses. The key inputs are the in-medium potentials felt by Lambda and Sigma hyperons. This paper shows that charged-current neutrino and antineutrino interactions produce those hyperons inside argon nuclei, so their final-state interactions leave measurable imprints on trapped fractions, escape momenta and a clean Sigma-plus tag. A Monte Carlo study of SBND and DUNE beams finds that these observables respond monotonically to the potentials. After detector reconstruction and marginalisation over the low-density density slope, the projected uncertainty on the Lambda potential is roughly 6 MeV, with systematics independent of hypernuclear spectroscopy; the Sigma potential is limited by poorly known hyperon-nucleon cross sections. The same measured depths can be fed directly into a relativistic mean-field equation of state to predict neutron-star mass-radius curves and tidal deformabilities, anchoring the low-density side of the hyperon puzzle with a new terrestrial handle.

Core claim

Charged-current (anti)neutrino interactions create Lambda and Sigma hyperons inside the nucleus; their subsequent final-state interactions therefore encode the in-medium potentials U_Lambda and U_Sigma. In the StrangeMC simulation the trapped-Lambda fraction, escaping hyperon momenta and a kaon-vetoed Sigma-plus tag respond monotonically to those potentials at SBND and DUNE. A detector-level Fisher forecast, marginalised over the low-density slope, yields delta U_Lambda approximately 6 MeV with systematics distinct from hypernuclear data, while U_Sigma remains limited by hyperon-nucleon cross-section uncertainties.

What carries the argument

The StrangeMC multi-layer Monte Carlo that produces hyperons via Cabibbo-suppressed and associated-production channels, then propagates them through an intranuclear cascade under the two-parameter potential form U_Y(rho) = U_Y(rho_0)(rho/rho_0)^gamma plus a high-density turn-over term, converting the exit-energy shift into the trapped fraction and momentum observables used for the forecast.

Load-bearing premise

The calculation treats a simple energy-threshold criterion (outgoing energy below the free mass) as a faithful enough proxy for real hypernucleus capture that the relative trapped-Lambda fraction stays a clean, monotonic measure of the Lambda potential.

What would settle it

A dedicated antineutrino reverse-horn run at SBND that reconstructs the trapped-Lambda fraction and mean escaping momenta and finds no monotonic dependence on the assumed potential depths, or a measured bias larger than the projected 6 MeV once the transport prescription is varied.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • The low-density anchors U_Lambda(rho_0) and U_Sigma(rho_0) become measurable with accelerator neutrinos at the few-MeV level, independent of hypernuclear spectroscopy.
  • Those measured depths can be inserted directly into a relativistic mean-field equation of state to generate mass-radius and tidal-deformability sequences that confront NICER and gravitational-wave data.
  • A kaon-vetoed Sigma-plus sample of a few hundred to a few thousand events provides a pure final-state-interaction tag of the potential difference that controls Sigma-minus onset.
  • An unplanned SBND reverse-horn antineutrino run becomes scientifically motivated because its soft quasi-elastic sample cleanly separates U_Lambda from U_Sigma.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • If the 6 MeV Lambda precision is realised, tension between hypernuclear and neutrino extractions would become a quantitative diagnostic of transport modelling rather than an uncontrolled systematic.
  • The same cascade framework could be re-used for other light nuclei or for proton-beam hyperon production to test whether the trapped-fraction observable is nucleus-independent.
  • Once external hyperon-nucleon cross-section data improve, the presently systematics-limited Sigma handle could become competitive, tightening the Sigma-minus onset density that most strongly affects the maximum mass.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

2 major / 4 minor

Summary. The paper proposes that charged-current (anti)neutrino interactions on argon produce Λ and Σ hyperons inside the nucleus, so that final-state interactions encode the in-medium potentials U_Y(ρ) that control hyperon onset in neutron-star matter. Using the StrangeMC cascade with a turn-over potential (Eq. 1) that cleanly separates the low-density pair (U_Y(ρ_0), γ) from the high-density parameters, the authors show that the trapped-Λ fraction, escaping hyperon momenta and a kaon-vetoed FSI-Σ^+ tag respond monotonically to U_Λ and U_Σ at SBND and DUNE (Fig. 1). A detector-level Fisher forecast, after marginalising over γ, yields δU_Λ ≈ 6 MeV with systematics distinct from hypernuclear spectroscopy; U_Σ remains limited by YN cross-section uncertainties. The same potentials are inserted into a GM1 RMF EOS to produce mass–radius and tidal sequences, and a joint Bayesian fit with external priors is used to illustrate the mapping onto M_max.

Significance. If the forecast holds, accelerator neutrinos supply a fourth, independent terrestrial handle on the low-density hyperon potentials that anchor hyperonic equations of state. The work is unusually candid about the hierarchy of systematics (U_Λ robust at the ~10 MeV level after γ marginalisation; U_Σ systematics-limited by YN scattering) and about the fact that the neutrino data constrain only the sub-saturation anchor, not the high-density slope or M_max itself. The end-to-end chain from reconstructed observables through a standard RMF to multi-messenger quantities, together with the explicit quantification of transport and FSI biases, makes the proposal falsifiable and complementary to hypernuclear, Σ-atom and heavy-ion constraints.

major comments (2)
  1. Simulation and Projected reach sections: the central U_Λ observable is the relative trapped-Λ fraction, defined by the transport-level energy-threshold criterion E_out ≤ m_Y. The paper itself reports that switching from an exit-energy shift to force-integrated gradient transport already biases the extracted U_Λ by -6 MeV—comparable to the γ-marginalised statistical reach of ~5.6 MeV. Because this proxy is load-bearing for the claim that the trapped fraction is a clean, monotonic handle, the forecast should either (i) adopt the more conservative of the two prescriptions as baseline and quote the difference as a systematic, or (ii) demonstrate that the relative (not absolute) trapped fraction remains stable under a broader set of capture models (e.g., simple optical-potential capture probabilities).
  2. Projected reach and joint-fit paragraphs: the Fisher matrix is built from a schematic liquid-argon response (V0 tag, momentum smearing, tracking thresholds). While the paper states robustness to detector model and exposure, no quantitative variation of reconstruction efficiency, background contamination of the kaon-vetoed Σ^+ sample, or Λ o pπ^- purity is shown. A short table or appendix quantifying how δU_Λ and δU_Σ degrade under ±20 % efficiency or realistic background fractions would make the projected 6 MeV and 3 MeV figures more credible for experimental planning.
minor comments (4)
  1. Fig. 1 caption and surrounding text: the slight rise of ⟨p_Λ⟩ with deeper U_Λ is correctly attributed to survivor selection, but a one-sentence clarification that the trapped fraction (not the momentum) is the primary U_Λ handle would prevent misreading of the blue curve.
  2. Eq. (1) and the EOS section: the high-density parameters (c_Y, β) are stated to vanish at ρ_0, yet the joint fit still varies c_Λ. A brief remark that the neutrino likelihood is independent of c_Λ (as already noted later) would make the separation of scales fully transparent at first appearance.
  3. References: the companion paper [21] is cited for all technical details of StrangeMC; a short public note or repository link (even if the full code remains on request) would improve reproducibility for readers who wish to re-run the Fisher forecast.
  4. Typographical consistency: “Σ− atoms” versus “Σ- atoms” and occasional missing spaces around MeV units appear in a few places; a light copy-edit pass would remove them.

Circularity Check

0 steps flagged

No significant circularity: neutrino sensitivity forecast is independent of the EOS mapping; M_max posterior is explicitly driven by external high-density prior, not by construction from neutrino inputs.

full rationale

The paper is a Monte-Carlo sensitivity/forecast study. StrangeMC generates hyperon observables (trapped-Λ fraction, escaping momenta, kaon-vetoed Σ+ tag) as a function of input potentials U_Y(ρ0) and γ; a Fisher matrix then projects the statistical reach on those same parameters after detector reconstruction. That chain is not circular: the observables are simulated responses, not fitted and then re-predicted. The identical functional form (Eq. 1) is deliberately reused when the extracted U_Y(ρ0) is inserted into a GM1 RMF to set the scalar couplings g_σY; this is a consistent mapping, not a reduction of a prediction to its input. The paper states repeatedly that the resulting M_max posterior (2.21 +0.04/-0.15 M_⊙) is delivered by the external heavy-ion c_Λ prior after marginalisation, not by the neutrino likelihood, which has no sensitivity above ρ0. Companion-paper self-citations supply simulation details but do not supply an unverified uniqueness theorem or force the central claim; the response curves and Fisher ellipses are shown in the present text. No self-definitional loop, no fitted-input-called-prediction, and no ansatz smuggled as a theorem appear. The only mild self-reference is ordinary method-paper citation, scored at 1.

Axiom & Free-Parameter Ledger

5 free parameters · 4 axioms · 2 invented entities

The central forecast rests on a phenomenological potential form, a transport-level trapping proxy, calibrated elementary cross sections, a schematic detector response, and a standard RMF EOS whose high-density behaviour is set by external priors. No new fundamental entities are introduced; the free parameters are the usual potential depths and slopes plus the YN cross-section uncertainties that dominate the U_Σ budget.

free parameters (5)
  • U_Λ(ρ0), U_Σ(ρ0)
    Baseline depths (-28 MeV, +30 MeV) taken from hypernuclear and Σ-atom phenomenology; the forecast measures deviations around these values.
  • γ (low-density slope)
    Exponent in the potential form; data give only δγ≈0.8, so it is marginalized and contributes the dominant statistical degradation of δU_Λ from 0.3 to ~6 MeV.
  • c_Y, β (high-density parameters)
    Control supra-saturation stiffness; fixed by external heavy-ion priors or left free; neutrino data have no sensitivity.
  • YN elastic and conversion cross sections
    O(50 %) uncertainty produces ~150 MeV bias on extracted U_Σ; treated as the leading systematic.
  • GM1 RMF couplings and SU(6) vector ratios
    Standard mean-field parameter set used to map U_Y into M_max and tidal deformability; choice of GM1 vs GM3 already spans ~0.3 M_☉.
axioms (4)
  • ad hoc to paper Hyperon potential of the turn-over form U_Y(ρ)=U_Y(ρ0)(ρ/ρ0)^γ + c_Y[(ρ/ρ0)^β-(ρ/ρ0)^γ] cleanly separates low- and high-density regimes.
    Eq. (1); chosen so that the second term vanishes at saturation and the neutrino data probe only (U_Y(ρ0),γ).
  • domain assumption An energy-threshold criterion E_out ä m_Y is a sufficient transport-level proxy for hypernucleus capture; the relative trapped fraction remains monotonic in U_Λ.
    Simulation section; absolute formation rates are not claimed, but the relative observable is treated as robust.
  • domain assumption Elementary ΔS=0,1 production cross sections calibrated to existing data and MicroBooNE are adequate for the forecast.
    Companion paper and citations [13,17-20]; associated-production model spread of factor 3-5 is acknowledged.
  • domain assumption Standard nonlinear σ-ω-ρ RMF with GM1 parameters and SU(6) vector couplings correctly maps measured U_Y(ρ0) into the high-density EOS.
    From the potential to the maximum mass section; the paper notes that density-dependent or chiral-EFT frameworks would widen the M_max band further.
invented entities (2)
  • StrangeMC multi-layer Monte Carlo no independent evidence
    purpose: Generate bound-nucleon initial state, ΔS=0,1 production, and hyperon intranuclear cascade with the potential of Eq. (1) for SBND and DUNE beams.
    Described as forked from LUNAR and fully specified only in the companion paper; no independent public validation yet.
  • Kaon-vetoed FSI-Σ+ tag no independent evidence
    purpose: Isolate pure final-state-interaction Σ+ events that the weak vertex cannot produce (ΔS=ΔQ forbids them), providing a low-background probe of U_Σ-U_Λ.
    Kinematically allowed by charge exchange; the paper estimates 1-1.4 % of antineutrino QE hyperons. Conceptually new as a neutrino observable, though the underlying FSI physics is standard.

pith-pipeline@v1.1.0-grok45 · 12647 in / 3841 out tokens · 28825 ms · 2026-07-13T04:21:15.398784+00:00 · methodology

0 comments
read the original abstract

Charged-current (anti)neutrino interactions create $\Lambda$ and $\Sigma$ hyperons \emph{inside} the nucleus, making hyperon final-state interactions a terrestrial probe of the in-medium potentials that govern hyperon onset in neutron stars. In the StrangeMC simulation the trapped-$\Lambda$ fraction, escaping hyperon momenta, and a kaon-vetoed $\Sigma^+$ tag respond monotonically to $U_\Lambda$ and $U_\Sigma$ at SBND and DUNE. Marginalised over the low-density slope, a forecast gives $\delta U_\Lambda\approx6\MeV$ with systematics distinct from hypernuclear data; $U_\Sigma$ is limited by hyperon--nucleon cross sections.

Figures

Figures reproduced from arXiv: 2607.09255 by Jaroslaw Nowak.

Figure 2
Figure 2. Figure 2: FIG. 2. GM1 RMF predictions with the hyperon couplings [PITH_FULL_IMAGE:figures/full_fig_p002_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Projected 1 [PITH_FULL_IMAGE:figures/full_fig_p003_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. The propagated posterior on the neutron-star maxi [PITH_FULL_IMAGE:figures/full_fig_p003_4.png] view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

28 extracted references · 24 linked inside Pith

  1. [1]

    P. B. Demorest, T. Pennucci, S. M. Ransom, M. S. E. Roberts, and J. W. T. Hessels, Nature467, 1081 (2010), arXiv:1010.5788

    Pith/arXiv arXiv 2010

  2. [2]

    Antoniadiset al., Science340, 1233232 (2013), arXiv:1304.6875

    J. Antoniadiset al., Science340, 1233232 (2013), arXiv:1304.6875

    Pith/arXiv arXiv 2013

  3. [3]

    Fonsecaet al., Astrophys

    E. Fonsecaet al., Astrophys. J. Lett.915, L12 (2021), arXiv:2104.00880

    Pith/arXiv arXiv 2021

  4. [4]

    Chatterjee and I

    D. Chatterjee and I. Vida˜ na, Eur. Phys. J. A52, 29 (2016), arXiv:1510.06306

    Pith/arXiv arXiv 2016

  5. [5]

    Tolos and L

    L. Tolos and L. Fabbietti, Prog. Part. Nucl. Phys.112, 103770 (2020), arXiv:2002.09223

    Pith/arXiv arXiv 2020

  6. [6]

    G. F. Burgio, H.-J. Schulze, I. Vida˜ na, and J.-B. Wei, Prog. Part. Nucl. Phys.120, 103879 (2021), arXiv:2105.03747

    Pith/arXiv arXiv 2021

  7. [7]

    A. Gal, E. V. Hungerford, and D. J. Millener, Rev. Mod. Phys.88, 035004 (2016), arXiv:1605.00557

    Pith/arXiv arXiv 2016

  8. [8]

    P. K. Sahaet al., Phys. Rev. C70, 044613 (2004), arXiv:nucl-ex/0405031

    Pith/arXiv arXiv 2004

  9. [9]

    Bednarek, P

    I. Bednarek, P. Haensel, J. L. Zdunik, M. Bejger, and R. Ma´ nka, Astron. Astrophys.543, A157 (2012), arXiv:1111.6942

    Pith/arXiv arXiv 2012

  10. [10]

    Ohnishi, S

    A. Ohnishi, S. Jinno, Y. Murase, and Y. Nara, EPJ Web Conf. (2022), arXiv:2210.17202

    Pith/arXiv arXiv 2022

  11. [11]

    Y. Nara, S. Jinno, Y. Murase, and A. Ohnishi, Phys. Rev. C (2022), arXiv:2208.01297

    Pith/arXiv arXiv 2022

  12. [12]

    Kohno, Y

    M. Kohno, Y. Nara,et al., arXiv preprint (2025), arXiv:2508.19560

    Pith/arXiv arXiv 2025

  13. [13]

    S. K. Singh and M. J. Vicente Vacas, Phys. Rev. D74, 053009 (2006), arXiv:hep-ph/0606235

    Pith/arXiv arXiv 2006

  14. [14]

    Thorpe, J

    C. Thorpe, J. Nowak, K. Niewczas, J. T. Sobczyk, and C. Juszczak, Phys. Rev. C104, 035502 (2021), arXiv:2010.12361

    Pith/arXiv arXiv 2021

  15. [15]

    Benitez Galan, L

    A. Benitez Galan, L. Alvarez-Ruso, M. Rafi Alam, I. Ruiz Simo, and M. J. Vicente Vacas, Phys. Rev. D 109, 033001 (2024), arXiv:2305.17004

    Pith/arXiv arXiv 2024

  16. [16]

    MicroBooNE Collaboration, Phys. Rev. Lett.130, 231802 (2023)

    2023

  17. [17]

    Rafi Alam, I

    M. Rafi Alam, I. Ruiz Simo, M. Sajjad Athar, and M. J. Vicente Vacas, Phys. Rev. D82, 033001 (2010), arXiv:1004.5484

    Pith/arXiv arXiv 2010

  18. [18]

    M. R. Alam, M. S. Athar, S. Chauhan, and S. K. Singh, Int. J. Mod. Phys. E25, 1650010 (2016), arXiv:1303.5924 [hep-ph]

    Pith/arXiv arXiv 2016

  19. [19]

    Fatima, M

    A. Fatima, M. Sajjad Athar, and S. K. Singh, Phys. Rev. D (2025), arXiv:2507.20754 [hep-ph]

    Pith/arXiv arXiv 2025

  20. [20]

    MicroBooNE Collaboration, Phys. Rev. Lett.135, 251804 (2025), arXiv:2503.00291 [hep-ex]

    arXiv 2025

  21. [21]

    J. A. Nowak (2026), companion paper, submitted to Phys. Rev. D

    2026

  22. [22]

    N. K. Glendenning and S. A. Moszkowski, Phys. Rev. Lett.67, 2414 (1991)

    1991

  23. [23]

    Hinderer, Astrophys

    T. Hinderer, Astrophys. J.677, 1216 (2008), arXiv:0711.2420

    Pith/arXiv arXiv 2008

  24. [24]

    LIGO Scientific Collaboration and Virgo Collaboration, Phys. Rev. Lett.119, 161101 (2017), arXiv:1710.05832

    Pith/arXiv arXiv 2017

  25. [25]

    LIGO Scientific Collaboration and Virgo Collaboration, Phys. Rev. Lett.121, 161101 (2018), arXiv:1805.11581

    Pith/arXiv arXiv 2018

  26. [26]

    LIGO Scientific Collaboration and Virgo Collaboration, Phys. Rev. X9, 011001 (2019), arXiv:1805.11579

    Pith/arXiv arXiv 2019

  27. [27]

    M. C. Milleret al., Astrophys. J. Lett.918, L28 (2021), arXiv:2105.06979

    Pith/arXiv arXiv 2021

  28. [28]

    Nowak, LUNAR: a Monte Carlo generator for bound- nucleon decay in liquid argon (2026), arXiv:2606.30872 [hep-ph]

    J. Nowak, LUNAR: a Monte Carlo generator for bound- nucleon decay in liquid argon (2026), arXiv:2606.30872 [hep-ph]

    Pith/arXiv arXiv 2026