Pith. sign in

REVIEW 3 major objections 5 minor 19 references

First measurement of flux of the neutron background induced by accelerated neutrinos at the J-PARC facility

T0 review · 3 major / 5 minor · reviewed 2026-07-12 · grok-4.5

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

Pith's one-line read Neutron background flux from the J-PARC neutrino beam is measured at 1.45 × 10^{-7} cm^{-2} s^{-1} POT^{-1}.

desk verdict First usable site-specific neutron flux at J-PARC under an exponential model; solid data chain, but the quoted number is model-dependent by construction. read the letter →

arxiv 2607.02783 v1 pith:HCD4YVVW submitted 2026-07-02 physics.ins-det hep-ex

Hiroshi Ito , Ryo Shibayama , Chisako Ise , Akihiro Minamino , Yusuke Koshio , Masaki Ishitsuka This is my paper · ORCID
classification physics.ins-dethep-ex
keywords neutronfluxJ-PARCneutrinobeambackgroundliquidscintillatorpulse-shapediscriminationNCQEprotonrecoilBGO
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

Neutrino experiments that search for nuclear-recoil signals from neutral-current quasi-elastic scattering are limited by accidental neutrons produced when the beam interacts in surrounding sand and concrete. This paper reports the first direct measurement of that neutron flux at the J-PARC neutrino facility. Using a liquid scintillator that identifies protons by pulse-shape discrimination, a plastic veto against charged particles, and 2.972 × 10^{20} protons on target, the collaboration selected 88 recoil-proton events between 0.98 and 11.60 MeV electron-equivalent. After efficiency and resolution corrections and a fit to Monte Carlo, they extract a flux of [1.45^{+0.22}_{-0.24} (stat.) ± 0.55 (sys.)] × 10^{-7} cm^{-2} s^{-1} POT^{-1} under an exponential energy spectrum. The number supplies a concrete background estimate for future short-baseline recoil detectors at the same site.

What carries the argument

Pulse-shape discrimination in a liquid scintillator, calibrated with a 65 MeV neutron beam, isolates proton recoils from gamma and heavy-ion backgrounds; the observed recoil spectrum is then fitted by re-weighting an isotropic Monte Carlo sample whose kinetic-energy distribution is forced to an exponential form Γ(E_n) = A exp(-E_n/λ).

What would settle it

An independent measurement of the neutron energy spectrum (or angular distribution) at the same location that cannot be described by a single exponential would force a revision of the quoted flux.

Watch

Extended reading notes

Core claim

With 2.972 × 10^{20} protons on target in neutrino mode, 88 neutron-induced proton recoils were observed in a liquid scintillator between 0.98 and 11.60 MeVee after pulse-shape discrimination and charged-particle veto. Comparing the measured spectrum with a re-weighted Monte Carlo simulation yields a neutron flux of [1.45^{+0.22}_{-0.24} (stat.) ± 0.55 (sys.)] × 10^{-7} cm^{-2} s^{-1} POT^{-1} when the kinetic-energy spectrum is taken to be exponential.

Load-bearing premise

The neutron kinetic-energy spectrum is assumed to be a pure exponential whose two free parameters are fitted only to the observed recoil spectrum, and the neutrons are assumed to arrive isotropically.

Editorial extensions

If this is right

  • Future short-baseline NCQE experiments at J-PARC can now estimate their neutron-induced proton-recoil background rate per POT using the reported flux.
  • The measured flux is independent of beam power in the 650–800 kW range, so background scales simply with exposure.
  • The same apparatus can be used to map flux versus detector location and to compare neutrino versus antineutrino modes.
  • Systematic uncertainties are dominated by PSD efficiency (36 %), giving a clear target for improved calibration.

Reading between the lines

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

  • If the true spectrum is harder or softer than the fitted exponential, the absolute flux (defined as the integral Aλ) will shift even if the recoil rate in the 1–12 MeVee window remains fixed.
  • A compact BGO array with the same liquid-scintillator tag could convert this background measurement into a simultaneous NCQE cross-section measurement once the neutron component is subtracted.
  • The large PSD systematic suggests that a dual-scintillator or capture-gated detector would reduce the dominant uncertainty for the next generation of flux measurements.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

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 the first measurement of the accelerator-neutrino-induced neutron background flux at the J-PARC neutrino facility (NM building B2 floor). Using a plastic-scintillator veto, BGO crystals, and a BC-501A liquid scintillator with pulse-shape discrimination, the authors select 88 proton-recoil events in the 0.98–11.60 MeVee window from a neutrino-mode exposure of 2.972 × 10^20 POT. After beam-on/off subtraction, charged-particle veto, and PSD cuts calibrated at RCNP, they re-weight a Geant4 neutron sample with an assumed exponential spectrum Γ(E_n) = A exp(-E_n/λ) (Eq. 1) and extract the integrated flux φ = Aλ = [1.45^{+0.22}_{-0.24} (stat.) ± 0.55 (sys.)] × 10^{-7} cm^{-2} s^{-1} POT^{-1}. Three successive setups yield consistent rates; the dominant systematic is the PSD selection (36 %). The result is intended as input for future short-baseline NCQE measurements.

Significance. A site-specific neutron-flux measurement at J-PARC is of practical value for any short-baseline NCQE or nuclear-recoil experiment planned in the NM hall. The experimental chain—energy calibration with tagged γ sources, 65 MeV neutron-beam PSD calibration at RCNP, beam-on/off subtraction, and multi-setup consistency—is carefully documented and constitutes a solid first data point. The large PSD systematic and the model dependence of the quoted flux limit the precision, but the measurement remains useful under the stated assumptions and will help guide detector placement and shielding design.

major comments (3)
  1. Sec. 3.4.3, Eq. (1) and Sec. 5.2: the reported flux is defined as φ = Aλ, the integral of a pure exponential from 0 to ∞. Neutrons are generated flat in 1–10 000 MeV and re-weighted; the observed recoil window (0.98–11.60 MeVee) maps, after quenching, to only a limited band of true E_n. Consequently A and λ are strongly correlated (Fig. 7) and the high-energy tail that still contributes to the integral is essentially unconstrained. The paper should either (i) quote the flux only inside a well-defined energy window that is actually constrained by the data, or (ii) demonstrate that alternative spectral shapes (power-law, broken exponential, etc.) yield fluxes consistent within the quoted systematics. Without such a check the single number cannot be treated as a spectrum-independent flux.
  2. Sec. 3.4.4: an isotropic angular distribution is assumed because production points and wall geometry are poorly known. Because the detectors are small and sit near walls, a strongly forward- or wall-directed component would change both the effective solid angle and the recoil spectrum. A brief sensitivity study (e.g., cos θ distributions weighted toward the beam axis or the nearest wall) should be added so that readers can judge how much of the 0.55 systematic already covers this uncertainty.
  3. Sec. 5.3: the PSD systematic is quoted as 36 % and dominates the total error. The text states that the PSD cut boundaries were varied, but does not specify the variation range or the resulting change in selected event count. Given that the cut is taken from a 65 MeV mono-energetic beam while the in-situ spectrum is continuous and softer, a more transparent propagation (or an alternative PSD metric) is needed before the 36 % figure can be accepted as complete.
minor comments (5)
  1. Table 1 and abstract: total POT is written both as 2.972 × 10^20 and as 29.72 × 10^19; standardize to a single notation.
  2. Fig. 3 caption: “root temperature” should be “room temperature”.
  3. Throughout: “Liqid scintillator” / “fl ame” / “P ARC” spacing errors should be corrected.
  4. Sec. 4.2: the pulse-separation limit of 0.4 µs is stated without a reference or measurement; a short justification would help.
  5. Appendix A, Table A.1: σ_0 is negative; a brief remark that this is an empirical offset (not a physical resolution) would avoid confusion.

Circularity Check

1 steps flagged · score 2.0 of 10

Flux is the integral of a two-parameter exponential ansatz fitted to the same recoil spectrum; ordinary parametric extraction under an explicit model assumption, not a self-referential derivation.

  1. fitted input called prediction [Sec. 3.4.3 Eq. (1) and Sec. 5.2]
    "We assumed the neutron energy spectrum to be an exponential function of the neutron kinetic energy En, Γ(En)=A exp(-En/λ) ... The neutron flux is estimated as φ=Aλ, which corresponds to the integral of Eq. (1) over the neutron energy En from zero to infinity. ... the combined neutron flux is determined to be (1.45^{+0.22}_{-0.24})×10^{-7} cm^{-2} s^{-1} POT^{-1}"

    A and λ are free parameters fitted by re-weighting the simulated recoil spectrum to the same data that define the measurement; the quoted flux is then exactly the product Aλ by definition of the integral of the ansatz. The numerical result is therefore the fitted normalization of the assumed shape, not an independent integral of a measured rate. Mild and fully disclosed; standard for parametric flux extraction.

full rationale

The paper is an experimental flux measurement, not a first-principles derivation or independent prediction. Neutrons are generated flat in energy, re-weighted by free parameters A and λ of an assumed exponential Γ(E_n)=A exp(-E_n/λ), and the parameters are fitted to the observed proton-recoil spectrum (0.98–11.60 MeVee) via χ²; the reported flux is then defined as φ≡Aλ (the integral of that ansatz). This is transparent model-dependent spectral fitting, not a circular claim that forces a result by construction or renames a tautology as a prediction. Quenching factors, energy resolution, and PSD cuts come from external calibrations (RCNP beam test and literature). No uniqueness theorem, self-citation chain, or ansatz smuggled via prior author work is load-bearing for the central number. The large A–λ correlation and unconstrained high-energy tail are model-dependence issues, not circularity. Score 2 reflects only the mild definitional step φ=Aλ after the fit; the measurement remains self-contained under the stated assumptions.

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

The central flux number rests on two free spectral parameters fitted to the data, on literature quenching factors, on an isotropic-direction assumption, and on standard detector-response models. No new physical entities are postulated.

free parameters (2)
  • A (exponential amplitude) = 1.0 × 10^5 (best-fit combined)
    Normalization of the assumed neutron kinetic-energy spectrum; fitted jointly with λ to the three recoil spectra.
  • λ (exponential decay constant) = 1.45 × 10^2 MeV (best-fit combined)
    Characteristic energy scale of the assumed neutron spectrum; fitted jointly with A; flux defined as Aλ.
assumptions (4)
  • ad hoc to paper Neutron kinetic-energy spectrum is exponential: Γ(E_n) = A exp(-E_n/λ)
    Stated in Sec. 3.4.3; no independent measurement of the spectrum shape at the site is available.
  • ad hoc to paper Neutron directions are isotropic
    Stated in Sec. 3.4.4 because production points and wall geometry are poorly known.
  • domain assumption Quenching factors for protons and carbon recoils in BC-501A follow the literature curves of Refs. 17 and 18
    Adopted without re-measurement for the present liquid scintillator (Sec. 3.4.2).
  • domain assumption Geant4.11.1 correctly models neutron scattering and energy deposition inside the detectors once the spectrum is re-weighted
    Standard Monte-Carlo tool used throughout Sec. 3.4.

how reviews work

0 comments
Cite this review

Pith. "Pith review of First measurement of flux of the neutron background induced by accelerated neutrinos at the J-PARC facility." pith.science (2026). https://pith.science/paper/HCD4YVVW

@misc{pith2026260702783,
  author       = {Pith},
  title        = {Pith review of: First measurement of flux of the neutron background induced by accelerated neutrinos at the J-PARC facility},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/HCD4YVVW}},
  note         = {Machine review of arXiv:2607.02783}
}
abstract

This article reports the measurement of the neutron background flux accidentally induced by the neutrino beam at the J-PARC facility in Japan. In particular, neutrino-nucleus neutral-current quasi-elastic (NCQE) scattering events, $\nu~(\bar{\nu})+\rm{N}\to\nu~(\bar{\nu})+\rm{N'}$, can be obscured by a massive background of neutron-nucleus scattering events. The neutrons are produced within the materials (such as sand and concrete) located between the beam dump and the experimental area where the detectors are placed. We measured the accidental neutron events at the J-PARC neutrino facility using BGO and liquid scintillation detectors, with a plastic scintillation detector serving as an active veto counter. Based on a neutrino-mode data set of $2.972 \times 10^{20}$ POT, we observed 88 neutron-induced recoil proton events within an electron-equivalent recoil energy range of 0.98-11.60 $\rm MeV_{ee}$, selected via pulse-shape discrimination in the liquid scintillation detector. Accounting for the detection efficiency and resolution, and through comparison with simulations, the neutron flux was determined to be$[1.45^{+0.22}_{-0.24}~\rm{(stat.)} \pm 0.55~\rm{(sys.)}] \times 10^{-7}~\rm{cm}^{-2} \rm{s}^{-1}\rm{POT}^{-1}$,assuming an exponential neutron energy spectrum. This result will contribute to the evaluation of neutron backgrounds for neutrino experiments at the J-PARC facility.

Figures

Figures reproduced from arXiv: 2607.02783 by the authors.

Figure 1
Figure 1. Front-end data acquisition system overview. In peri [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Setup for periods 1, 2, and 3. The plastic scintillati [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Monitored root temperature in period-1, -2, and -3 (t [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figures from the paper (3 more)
Figure 4
Figure 4. Figure 4: Hit time distribution in the plastic scintillation d [PITH_FULL_IMAGE:figures/full_fig_p005_4.png]
Figure 5
Figure 5. Figure 5: Distribution of deposited energy versus the PSD valu [PITH_FULL_IMAGE:figures/full_fig_p005_5.png]
Figure 6
Figure 6. Figure 6: Observed deposited energy spectra in electron recoi [PITH_FULL_IMAGE:figures/full_fig_p006_6.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

19 extracted references

  1. [1]

    G. T. Garvey, W . C. Louis, D. H. White, Determination of proton strange form factors from νp elastic scattering, Phys. Rev. C 48 (1993) 761–765

    1993

  2. [2]

    Golan, et al., Extraction of axial mass and strangenes s values from the MiniBooNE neutral current elastic cross section measurement, Phys

    T. Golan, et al., Extraction of axial mass and strangenes s values from the MiniBooNE neutral current elastic cross section measurement, Phys. Rev. C 88 (2013) 024612

    2013

  3. [3]

    Abe, et al., First measurement of the strange axial cou - pling constant using neutral-current quasielastic intera c- tions of atmospheric neutrinos at KamLAND, Phys

    S. Abe, et al., First measurement of the strange axial cou - pling constant using neutral-current quasielastic intera c- tions of atmospheric neutrinos at KamLAND, Phys. Rev. D 107 (2023) 072006

    2023

  4. [4]

    Harada, et al., Search for Astrophysical Elec- tron Antineutrinos in Super-Kamiokande with 0.01% Gadolinium-loaded Water, ApJL 951 (2) (2023) L27

    M. Harada, et al., Search for Astrophysical Elec- tron Antineutrinos in Super-Kamiokande with 0.01% Gadolinium-loaded Water, ApJL 951 (2) (2023) L27

    2023

  5. [5]

    Honda, et al., Improvement of low energy atmospheric neutrino flux calculation using the JAM nuclear interac- tion model, Phys

    M. Honda, et al., Improvement of low energy atmospheric neutrino flux calculation using the JAM nuclear interac- tion model, Phys. Rev. D 83, 123001 (2011)

    2011

  6. [6]

    Honda, et al., Calculation of atmospheric neutrino flu x using the interaction model calibrated with atmospheric muon data, Phys

    M. Honda, et al., Calculation of atmospheric neutrino flu x using the interaction model calibrated with atmospheric muon data, Phys. Rev. D 75, 043006 (2007)

    2007

  7. [7]

    Hayato, A NEUTRINO INTERACTION SIMU- LA TIONPROGRAM LIBRARY NEUT, Acta

    Y . Hayato, A NEUTRINO INTERACTION SIMU- LA TIONPROGRAM LIBRARY NEUT, Acta. Physica Polonica B, 40, 9, 2477 (2009)

    2009

  8. [8]

    Hayato, L

    Y . Hayato, L. Pickering, The NEUT neutrino interac- tion simulation program library, Eur. Phys. J. Spec. Top. (2021) 230: 4469–4481

    2021

Show all 19 references
  1. [9]

    A. M. Ankowski, et al., Analysis of γ-Ray Production in Neutral-Current Neutrino-Oxygen Interactions at En- ergies above 200 MeV, Phys. Rev. Lett. 108, 052505

  2. [10]

    Abe, et al., Measurement of neutrino and antineutrin o neutral-current quasielasticlike interactions on oxygen by detecting nuclear deexcitation γrays, Phys

    K. Abe, et al., Measurement of neutrino and antineutrin o neutral-current quasielasticlike interactions on oxygen by detecting nuclear deexcitation γrays, Phys. Rev. D 100 (2019) 112009

    2019

  3. [11]

    Ommura, et al., A method to measure the quenching factor for recoil energy of oxygen in bismuth germanium oxide scintillators, Journal of Instrumentation 18 (04) (2023) T04006

    Y . Ommura, et al., A method to measure the quenching factor for recoil energy of oxygen in bismuth germanium oxide scintillators, Journal of Instrumentation 18 (04) (2023) T04006

    2023

  4. [12]

    Ito, et al., Measurement of neutrino-induced neu- tron and γ-ray background for BGO-based detectors at J-P ARC, Proposals for Nuclear and Particle Physics Ex- periments at J-P ARC

    H. Ito, et al., Measurement of neutrino-induced neu- tron and γ-ray background for BGO-based detectors at J-P ARC, Proposals for Nuclear and Particle Physics Ex- periments at J-P ARC

  5. [13]

    Abe, et al., T2K neutrino flux prediction, Phys

    K. Abe, et al., T2K neutrino flux prediction, Phys. Rev. D 87 (2013) 012001

    2013

  6. [14]

    Ovsiannikova, et al., The new experiment W AGASCI for water to hydrocarbon neutrino cross section measure- ment using the J-P ARC beam, J

    T. Ovsiannikova, et al., The new experiment W AGASCI for water to hydrocarbon neutrino cross section measure- ment using the J-P ARC beam, J. Phys.: Conf. Ser. 675 (2016) 012030

    2016

  7. [15]

    Fukuda, et al., First neutrino event detection with n u- clear emulsion at J-P ARC neutrino beamline, Theor

    T. Fukuda, et al., First neutrino event detection with n u- clear emulsion at J-P ARC neutrino beamline, Theor. Exp. Phys. 063C02

  8. [16]

    Fukuda, et al., First demonstration of an emulsion multi-stage shifter for accelerator neutrino experiments in J-P ARC T60, Theor

    T. Fukuda, et al., First demonstration of an emulsion multi-stage shifter for accelerator neutrino experiments in J-P ARC T60, Theor. Exp. Phys. 063H02

  9. [17]

    Y oshida, et al., Light output response of KamLAND liquid scintillator for protons and 12C nuclei, Nucl

    S. Y oshida, et al., Light output response of KamLAND liquid scintillator for protons and 12C nuclei, Nucl. Instr. Meth. 622 (3) (2010) 574–582

    2010

  10. [18]

    A we, et al., Measurement of proton quenching in a plastic scintillator detector, JINST 16 P02035 (2021)

    C. A we, et al., Measurement of proton quenching in a plastic scintillator detector, JINST 16 P02035 (2021)

    2021

  11. [19]

    Ito, et al., Analyzing the neutron and γ-ray emis- sion properties of an americium–beryllium tagged neutron source, Nucl

    H. Ito, et al., Analyzing the neutron and γ-ray emis- sion properties of an americium–beryllium tagged neutron source, Nucl. Instr. Meth. A 1057 (2023) 168701. 9

    2023

Pith tools

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