Pith. sign in

REVIEW 4 major objections 4 minor 27 references

Study of Stability and Consistency of EAS Thermal Neutron Detection at ENDA-64

T0 review · 4 major / 4 minor · reviewed 2026-08-07 · deepseek-v4-flash

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

Pith's one-line read ENDA-64, a 64-detector array for air-shower neutrons, is consistent to 6.85% and stable enough for cosmic-ray studies around the knee.

desk verdict First quantitative stability/consistency numbers for ENDA-64; rate-stability results are solid, but the consistency and sand-cube conclusions rest on an unvalidated background sample and a confounded comparison. read the letter →

arxiv 2506.10510 v1 pith:EUSDNGHJ submitted 2025-06-12 hep-ex physics.ins-det

Heng-Yu Zhang , Xin-Hua Ma , Tian-Lu Chen , Shu-Wang Cui , Danzengluobu , Wei Gao , Wen-Chao Gao , Xin-Rui Gao
show 22 more authors
Zi-Ao Gong Hai-Bing Hu Denis Kuleshov Kirill Kurinov Bing-Bing Li Fan-Ping Li Jia-Heng Li Yang Li Hu Liu Mao-Yuan Liu Ye Liu Xi-An Pan Da-Yu Peng Yao-Hui Qi Dong Qu Oleg Shchegolev Yuri Stenkin Tian-Shuang Yang Li-Qiao Yin Hui-Qian Zhang Liang-Wei Zhang Shi-Yuan Zhang
This is my paper
classification hep-exphysics.ins-det PACS 29.40.Mc95.55.Vj96.50.sd
keywords cosmicraysextensiveairshowersthermalneutrondetectiondetectorcalibrationstabilitysandcubeskneeregionENDA-64
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 reports on the calibration, stability, and consistency of ENDA-64, a 64-detector array that records thermal neutrons produced when cosmic-ray air showers strike the ground and surrounding material. The authors want to show that the array is stable and uniform enough to feed cosmic-ray energy-spectrum and composition measurements around the knee, where hadronic signals matter. They find a maximum per-detector calibration inconsistency of 6.85%, a maximum event-rate instability of 4.68%, and a maximum thermal-neutron-rate instability of 11.0%. They also show that clusters without sand cubes under the detectors disagree by up to 18%, and that sand cubes reduce rainwater effects while increasing neutron collection from shower events.

What carries the argument

The central calibration object is the steepness parameter $b$ of the pulse-height distribution, fit by $y=10^{a-bx}$; because $b$ reflects the energy deposited by one captured thermal neutron, equalizing it across detectors by adjusting high voltage is the paper's definition of consistency. Signal/background separation uses the ratio of the triggered-event neutron distribution to the rate-normalized M0 background distribution, with $\lg n > 1.2$ taken as the regime where real EAS events dominate. Stability is measured by fitting an exponential to the time-difference distribution between adjacent events, giving the event rate, and by tracking the thermal neutron rate. The sand-cube effect is assessed by comparing cluster No. 4's distributions across dry and rainy seasons and with and without sand cubes.

What would settle it

Recompute the triggered-to-background ratio using an independent background sample, such as randomly triggered 20 ms gates from non-coincident detector signals, and check whether $\lg n > 1.2$ still marks the EAS-dominated region. If the independent background gives a different threshold, or if the 18% cluster inconsistency changes substantially when the background is subtracted event-by-event, the paper's stability and consistency conclusions would need revision.

Watch

Extended reading notes

Core claim

The central claim is that ENDA-64, after adjusting photomultiplier high voltages, detects single thermal neutrons uniformly enough for physics use: the spread of the pulse-height steepness parameter $b$ within each cluster is at most 6.85% of its mean. Over February–June 2024 the event rate in a cluster is stable to 4.68% and the thermal neutron rate to 11.0%. Comparing the three clusters without sand cubes over August–December 2024, the thermal-neutron distributions agree to within 18% once low-statistics high-count bins are set aside. For cluster No. 4, the data show that sand cubes keep the neutron response stable through the rainy season and that, at neutron counts $\lg n > 1.2$, the sand-cube cluster collects more shower-produced thermal neutrons than the others, which the authors attribute to the 1 m elevation and the increased geometrical factor of the sand target.

Load-bearing premise

The load-bearing assumption is that the minute-long M0 background events, after being rescaled by event rate, represent the background that is actually mixed into the triggered events; if that scaling is wrong, the signal/background separation and the cluster comparisons built on it are biased.

Editorial extensions

If this is right

  • If these stability and consistency numbers hold, ENDA-64 can provide a hadronic, neutron-based channel alongside electromagnetic and muon measurements for cosmic-ray studies around the knee.
  • The quoted instabilities set systematic-error floors for the future energy spectrum and composition analysis; improving high-voltage calibration should reduce the 18% cluster inconsistency.
  • The sand-cube result implies that the array can lower its energy threshold for light-component (proton/helium) knee measurements, because the cubes raise neutron collection at higher neutron counts.
  • The measured inconsistency bounds justify the planned expansion to a larger array, since the per-cluster performance transfers to the larger detector count.
  • The seasonal and moisture sensitivity quantified here will need to be corrected or monitored continuously in the final physics analysis.

Reading between the lines

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

  • The paper leaves implicit that its 11.0% thermal-neutron-rate instability likely includes soil-moisture variation; an explicit correlation of the daily neutron rate with local rainfall or soil moisture would test whether sand cubes suppress this seasonal component.
  • The M0-ratio method for isolating real EAS events could be validated independently by comparing it with a Monte Carlo simulation of the shower neutron signal, which the authors say is planned; if the $\lg n > 1.2$ threshold shifts under simulation, the reported cluster inconsistencies would need revisiting.
  • Because only cluster No. 4 had sand cubes, the with/without-cube comparison conflates the sand-cube effect with any cluster-to-cluster calibration or environmental difference; a swap experiment—moving sand cubes to a different cluster—would isolate the effect.
  • The 18% inconsistency among the three no-cube clusters is stated after excluding high-count bins; if the full distribution is used, the quoted inconsistency could be larger, so future analyses should state the bin range over which the 18% applies.
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 / 4 minor

Summary. The paper reports an instrument-performance study of the ENDA-64 thermal-neutron detector array at LHAASO. The authors calibrate each detector by adjusting PMT high voltages so that the pulse-height distribution for single thermal neutrons is consistent across detectors, and they quantify the residual inconsistency per cluster as the ratio of the standard deviation to the mean of the fitted steepness parameter (maximum 6.85%). They then use time-difference distributions to derive event rates and thermal-neutron rates per cluster, quoting maximum instabilities of 4.68% and 11.0% based on the relative standard deviations of the daily rates. For consistency, they compare thermal-neutron multiplicity distributions among clusters without sand cubes, reporting a maximum difference of 18% after excluding high-multiplicity bins with low statistics. Finally, they compare cluster 4 (with sand cubes) to other clusters and to its own history in dry and rainy seasons, concluding that sand cubes protect the target material from rainwater and enhance collection of neutrons from extensive air showers.

Significance. The stability and calibration numbers (6.85%, 4.68%, 11.0%) are direct, data-driven measurements that are useful for the ENDA-64 collaboration and for any future cosmic-ray physics analysis using this array. The paper makes specific, falsifiable claims that can be checked against the published figures and tables, and the rate-stability analysis is straightforward and reproducible. The main value is the quantitative characterization of detector performance, which is a necessary prerequisite for the planned measurement of the cosmic-ray spectrum around the knee. However, the consistency and sand-cube conclusions are currently supported by comparisons that lack statistical uncertainties and controlled baselines, so the quantitative significance of those specific claims is not yet established.

major comments (4)
  1. [Section 3.2, Fig. 4] The definition of the real-EAS-dominated region (lg(n)>1.2) relies on normalizing M0 background events to triggered events using event rates. M0 events are 20 ms windows triggered by software once per minute, whereas triggered events are 20 ms gates opened by an EAS charged-particle trigger. These two types of windows can have different dead-time, afterpulse, and acceptance conditions, and it is not demonstrated that the M0 distribution faithfully represents the background inside triggered events. Because this threshold is used in Section 3.4 to interpret the sand-cube ratios, this unvalidated background treatment is load-bearing for the sand-cube conclusions.
  2. [Section 3.3, Fig. 7] The maximum inconsistency of 18% between the three clusters without sand cubes is stated after 'disregarding the higher neutron counts, which have low statistics and consequently large statistical errors,' but no quantitative definition of that cut is given, and the ratio panel in Fig. 7 has no error bars. Without a stated multiplicity range and statistical uncertainties, the 18% value cannot be reproduced or distinguished from a statistical fluctuation. Please specify the fitted lg(n) range, the binning, and the statistical uncertainty on the maximum deviation.
  3. [Section 3.4, Fig. 8] The rain-protection claim for sand cubes is based on a confounded comparison. The 'rainy season without sand cubes' period (Aug 16-Oct 14, 2023) is compared against a 'dry season with sand cubes' baseline (Feb 9-Jun 3, 2024), so differences between these two curves can be due to seasonal or inter-annual variations in soil moisture, temperature, or other environmental factors, rather than to the presence of the sand cubes alone. A controlled comparison would require the same cluster with and without sand cubes in the same season, or a quantitative measurement of soil moisture and rain during the two periods.
  4. [Section 3.4, Fig. 9] The comparison between cluster 4 (with sand cubes) and the average of the other three clusters (without sand cubes) shows a ratio that rises above 1 for lg(n)>1.2, but the text attributes this rise to the cluster being 1 m above the ground, not to the sand cubes themselves. This attribution contradicts the abstract and conclusion statement that 'the sand cubes help the cluster to increase collection of neutrons generated by EAS events.' Since the sand mass, the detector elevation, and the observation period are all different between the compared samples, the EAS-neutron enhancement cannot be uniquely assigned to the sand cubes without a simulation or a dedicated test.
minor comments (4)
  1. [Section 3.2] The text states that for lg(n)>1 the ratio R is greater than 1.2, but later uses lg(n)>1.2 as the real-EAS threshold; please make the threshold consistent and specify how it is chosen.
  2. [General figures] The ratio panels in Figs. 4, 7, 8, and 9 do not show statistical error bars. Adding error bars would allow the reader to assess whether features such as the 18% inconsistency and the sand-cube excess are significant.
  3. [Tables 1-3] The captions do not define M, S, and R; please state explicitly that M is the mean, S the standard deviation, and R the ratio S/M expressed as a percentage.
  4. [Section 3.4] The three periods in Fig. 8 are from different years (2023 and 2024), but this is only stated in the figure caption; the text should note the non-overlapping observation windows and the potential for long-term drift or seasonal effects.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: all reported quantities are direct data summaries or post-calibration residuals; no derivation reduces to its own inputs.

full rationale

This is an experimental characterization paper, not a derivation paper. The stability and consistency numbers are direct statistical summaries (mean M, standard deviation S, and ratio R = S/M) of measured quantities: the fitted steepness parameter b of pulse-height distributions, the event rate lambda fitted from time-difference distributions, and the thermal neutron rate. None of these results is used as an input to define itself: the 6.85% calibration inconsistency is the residual spread of b after high-voltage adjustment, not a parameter fitted to produce that spread; the 4.68% and 11.0% instabilities are spreads of measured rates over time; and the 18% cross-cluster inconsistency is a direct comparison of measured neutron-count distributions. The lg(n) > 1.2 region is empirically identified from triggered-to-M0 background ratios and then used consistently when interpreting sand-cube effects, but this is an operational selection, not a self-definitional reduction. Self-citations ([11, 23, 26]) are used only for trigger description and prior sand-cube installation background; they are not load-bearing uniqueness arguments and do not force any result. The acknowledged confounds in the sand-cube comparison (elevation, season, observation period) are correctness risks, not circularity. Therefore the central claims are self-contained measurements with no circular derivation.

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

This is an experimental performance paper; no new physical entities or hand-tuned parameters are introduced. The only fitted quantities are the parameters of empirical distributions (steepness b of the pulse-height spectrum, event rate lambda from the time-difference exponential), which are measured from the data rather than imposed. The analysis relies on a few domain assumptions about background subtraction and signal identification.

assumptions (4)
  • domain assumption The middle part of the pulse height distribution fitted by Eq. 2 corresponds to single thermal neutron captures in 10B.
    Location: Section 3.1. This assumption makes the spread of the fitted steepness parameter b into a measure of inter-detector calibration consistency; if this region includes multi-neutron or fast-neutron background, the reported 6.85% inconsistency would not truly reflect single-neutron detection.
  • domain assumption Event times follow an exponential distribution, so event rate can be extracted from the fit to Eq. 3.
    Location: Section 3.3. Used to convert time-difference histograms into daily event rates; any non-Poisson arrival structure would bias the stability estimates.
  • domain assumption M0 background events, normalized by event rate, represent the background contribution within triggered events.
    Location: Section 3.2. This normalization underpins the signal/background classification (lg(n)>1.2) used in all distribution comparisons; no independent validation is provided.
  • domain assumption The threshold lg(n)>1.2 separates real EAS events from background.
    Location: Section 3.2 and Figure 4. The ratio crossing 1.2 is treated as the signal region, but the same threshold is used for all clusters and seasons without a simulation-based justification.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Study of Stability and Consistency of EAS Thermal Neutron Detection at ENDA-64." pith.science (2026). https://pith.science/paper/EUSDNGHJ

@misc{pith2026250610510,
  author       = {Pith},
  title        = {Pith review of: Study of Stability and Consistency of EAS Thermal Neutron Detection at ENDA-64},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/EUSDNGHJ}},
  note         = {Machine review of arXiv:2506.10510}
}
read the original abstract

Introduction:Electron-Neutron Detector Array (ENDA) is designed to measure thermal neutrons produced by hadronic interactions between cosmic ray extensive air showers (EAS) and the surrounding environment as well as electrons around the cores of EAS. ENDA is located within Large High Altitude Air Shower Observatory (LHAASO). ENDA was expanded from an initial 16 detectors to 64 detectors in April 2023, so called ENDA-64, and has been running alongside LHAASO. The stability and consistency of neutron detection are crucial for laying a solid foundation for subsequent data analysis and physical results. Methods:We obtain the stability by studying variations of event rate and thermal neutron rate in each cluster and the consistency by comparing distribution of number of thermal neutrons between clusters. Additionally, we investigate the specific influences of the rainy and dry seasons, as well as the presence or absence of sand cubes under the detectors, to examine the environmental factors affecting neutron measurement performance. Results:The calibration results indicate good consistency in thermal neutron detection across the clusters, with the maximum inconsistency of 6.85%. The maximum instability of event rate and thermal neutron rate over time are 4.68% and 11.0% respectively. The maximum inconsistency between the clusters without the sand cubes is 18%. The use of sand cubes is effective in protecting the target material from rainwater, and the sand cubes help the cluster to increase collection of neutrons generated by EAS events.

Figures

Figures reproduced from arXiv: 2506.10510 by the authors.

Figure 1
Figure 1. (a) Aerial view of ENDA-64 within the LHAASO Array. ENDs and four clusters of ENDA and electron detectors (EDs) and muon detectors (MDs) of LHAASO are marked. (b) ENDs of cluster #4 with sand cubes. The detector mainly composed of scintillator, scintillation light collecting cone, Photomultiplier Tube (PMT), the front-end electronics (FEE), and black housing. The scintillator is made of ZnS(Ag) and B2O3 alloy, depos… view at source ↗
Figure 2
Figure 2. Schematic diagram of the END. 1 - the high-voltage input port, 2 - DIU connected to the 8th dynode of the PMT, 3 - IU connected to the 5th dynode of the PMT, 4 - a black tank for the detector housing, 5 - the PMT fixed holder, 6 - the PMT, 7 - the scintillation light collecting cone, and 8 - the scintillator. [27]. 4/13 [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Pulse height distribution. Black line shows the measurement from one detector; blue line represents the fitting function (Eq. 2) [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figures from the paper (6 more)
Figure 4
Figure 4. Figure 4: Thermal neutron distributions. Upper panel: [PITH_FULL_IMAGE:figures/full_fig_p006_4.png]
Figure 5
Figure 5. Figure 5: Time difference distribution. Black line shows time differences between adjacent events in one cluster during one day; red line represents the fitting function (Eq. 3) [PITH_FULL_IMAGE:figures/full_fig_p007_5.png]
Figure 6
Figure 6. Figure 6: Event rate variation. Measured in one cluster from February 9 to June 3, 2024 (MJD: Modified Julian Day). 3.3 Instability and inconsistency It is important to estimate the systematic uncertainty of the array by measuring instability and inconsistency using real data. E…
Figure 7
Figure 7. Figure 7: Thermal neutron consistency analysis. Upper panel: [PITH_FULL_IMAGE:figures/full_fig_p009_7.png]
Figure 8
Figure 8. Figure 8: Thermal neutron distributions with/without sand cubes. Upper [PITH_FULL_IMAGE:figures/full_fig_p010_8.png]
Figure 9
Figure 9. Figure 9: Thermal neutron comparison with sand cubes. Upper Panel: [PITH_FULL_IMAGE:figures/full_fig_p011_9.png]

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

27 extracted references · 27 canonical work pages

  1. [1]

    Antoni and et al

    T. Antoni and et al. Kascade measurements of energy spectra for elemental groups of cosmic rays: Results and open problems.Astroparticle Physics, 24(1):1–25, 2005

    work page 2005

  2. [2]

    Amenomori and et al

    M. Amenomori and et al. Are protons still dominant at the knee of the cosmic-ray energy spectrum.Physics Letters B, 632(1):58–64, 2006

    work page 2006

  3. [3]

    Bartoli and et al

    B. Bartoli and et al. Knee of the cosmic hydrogen and helium spectrum below 1 pev measured by argo-ybj and a cherenkov telescope of lhaaso.Phys. Rev. D, 92:092005, Nov 2015

    work page 2015

  4. [4]

    Measurements of all-particle energy spectrum and mean logarithmic mass of cosmic rays from 0.3 to 30 pev with lhaaso-km2a.Phys

    Zhen Cao and et al. Measurements of all-particle energy spectrum and mean logarithmic mass of cosmic rays from 0.3 to 30 pev with lhaaso-km2a.Phys. Rev. Lett., 132:131002, Mar 2024

    work page 2024

  5. [5]

    H¨ orandel

    J¨ org R. H¨ orandel. Models of the knee in the energy spectrum of cosmic rays. Astroparticle Physics, 21(3):241–265, 2004

    work page 2004

  6. [6]

    Cosmic ray energy spectrum from measurements of air showers.Frontiers of Physics, 8, 04 2013

    Thomas Gaisser and et al. Cosmic ray energy spectrum from measurements of air showers.Frontiers of Physics, 8, 04 2013

    work page 2013

  7. [7]

    S. I. Nikolsky. The Cause of the EAS Spectrum Break. InInternational Cosmic Ray Conference, volume 6 ofInternational Cosmic Ray Conference, page 105, January 1997

    work page 1997

  8. [8]

    On special features of the longitudinal development of extensive air showers and on the spectrum of cosmic rays.Physics of Atomic Nuclei, 71:98–110, 01 2011

    Yuri Stenkin. On special features of the longitudinal development of extensive air showers and on the spectrum of cosmic rays.Physics of Atomic Nuclei, 71:98–110, 01 2011

    work page 2011

Show all 27 references
  1. [9]

    Chapter 1 lhaaso instruments and detector technology *

    Xin-Hua Ma and et al. Chapter 1 lhaaso instruments and detector technology *. Chinese Physics C, 46(3):030001, mar 2022

    2022

  2. [10]

    Electron–neutron detector array (enda).Physics of Atomic Nuclei, 84:941–946, 11 2021

    Bing-Bing Li and et al. Electron–neutron detector array (enda).Physics of Atomic Nuclei, 84:941–946, 11 2021

    2021

  3. [11]

    Progress of electron–neutron detector array (enda)

    Da-Yu Peng and et al. Progress of electron–neutron detector array (enda). Physics of Atomic Nuclei, 86:1056–1062, 2023

    2023

  4. [12]

    Research on the knee region of cosmic ray by using a novel type of electron–neutron detector array.Frontiers of Physics, 19:44200, 2024

    Bing-Bing Li and et al. Research on the knee region of cosmic ray by using a novel type of electron–neutron detector array.Frontiers of Physics, 19:44200, 2024

    2024

  5. [13]

    Stenkin and et al

    Yuri V. Stenkin and et al. On the Neutron Bursts Origin.Modern Physics Letters A, 17(26):1745–1751, jan 2002. 12/13

    2002

  6. [14]

    Yuri V. Stenkin. On the prisma project.Nuclear Physics B - Proceedings Supplements, 196:293–296, 2009

    2009

  7. [15]

    Thermal neutron flux produced by eas at various altitudes

    Yuri Tenkin and et al. Thermal neutron flux produced by eas at various altitudes. Chinese Physics C, 37(1):015001, jan 2013

    2013

  8. [16]

    Yu. V. Stenkin and et al. Neutrons in extensive air showers.Physics of Atomic Nuclei, 70:1088–1099, 06 2007

    2007

  9. [17]

    Yu.V. Stenkin. Thermal neutrons in eas: a new dimension in eas study.Nuclear Physics B - Proceedings Supplements, 175-176:326–329, 2008

    2008

  10. [18]

    D. D. Djappuev and et al. Compact multicomponent array for EAS study (MULTICOM). InInternational Cosmic Ray Conference, page 822, 2001

    2001

  11. [19]

    Gromushkin and et al

    D. Gromushkin and et al. The protoprisma array for eas study: first results. Journal of Physics Conference Series, 409:2044–, 02 2013

    2013

  12. [20]

    Bartoli and et al

    B. Bartoli and et al. Detection of thermal neutrons with the prisma-ybj array in extensive air showers selected by the argo-ybj experiment.Astroparticle Physics, 81:49–60, 2016

    2016

  13. [21]

    Seasonal and lunar month periods observed in natural neutron flux at high altitude.Pure and Applied Geophysics, 174:2763–2771, 07 2017

    Yuri Stenkin and et al. Seasonal and lunar month periods observed in natural neutron flux at high altitude.Pure and Applied Geophysics, 174:2763–2771, 07 2017

    2017

  14. [22]

    Response of the environmental thermal neutron flux to earthquakes.Journal of Environmental Radioactivity, 208-209:105981, 2019

    Yuri Stenkin and et al. Response of the environmental thermal neutron flux to earthquakes.Journal of Environmental Radioactivity, 208-209:105981, 2019

    2019

  15. [23]

    Li and et al

    B.-B. Li and et al. Eas thermal neutron detection with the prisma-lhaaso-16 experiment.Journal of Instrumentation, 12(12):P12028, dec 2017

    2017

  16. [24]

    Performance of the thermal neutron detector array in yangbajing, tibet for cosmic ray eas detection.Astrophysics and Space Science, 365:123, 07 2020

    Mao-Yuan Liu and et al. Performance of the thermal neutron detector array in yangbajing, tibet for cosmic ray eas detection.Astrophysics and Space Science, 365:123, 07 2020

    2020

  17. [25]

    Yang and et al

    F. Yang and et al. Correlation between thermal neutrons and soil moisture measured by enda.Journal of Instrumentation, 18(05):P05020, may 2023

    2023

  18. [26]

    Influence of soil environment on performance of eas electron–neutron detector array.Astrophysics and Space Science, 367:75, 08 2022

    Di-Xuan Xiao and et al. Influence of soil environment on performance of eas electron–neutron detector array.Astrophysics and Space Science, 367:75, 08 2022

    2022

  19. [27]

    Readout electronics for thermal neutron detector array

    Chunhui Dong and et al. Readout electronics for thermal neutron detector array. Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, 946:162639, 2019. 13/13

    2019

Pith tools

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