Study of neutron production for 360 GeV cosmic muons

Basic Information

• Paper ID: 2409.04169
• Title: Study of neutron production for 360 GeV cosmic muons
• Authors: Xinshun Zhang, Jinjing Li, Shaomin Chen, et al. (JNE Collaboration)
• Classification: hep-ex (High Energy Physics - Experiment)
• Publication Date: December 12, 2024
• Paper Link: https://arxiv.org/abs/2409.04169

Abstract

The China Jinping Underground Laboratory (CJPL), characterized by its 2400-meter vertical rock overburden, serves as an excellent location for studying solar, terrestrial, and supernova neutrinos. Its unparalleled depth provides an opportunity to study high-energy cosmic ray muons with an average energy of approximately 360 GeV. Based on 1178 days of data collected since 2017 by the Jinping Neutrino Experiment (JNE) 1-ton prototype neutrino detector, this paper presents a detailed study of muon-related backgrounds. For the first time, we provide a detailed discussion of the significant impact of muon secondary particle leakage due to the detector's finite size on neutron yield measurements. Through analysis of 493 cosmic ray muon candidate events and 13.6±5.7 cosmogenic neutron candidate events, combined with comprehensive assessment of detection efficiency and uncertainties, we obtain a muon flux of (3.56±0.16stat.±0.10syst.)×10^-10 cm^-2s^-1 and a cosmogenic neutron yield in LAB liquid scintillator of (3.37±1.41stat.±0.31syst.)×10^-4 μ^-1g^-1cm^2.

Research Background and Motivation

1. Core Problem: This study aims to precisely measure the flux and cosmogenic neutron yield of high-energy cosmic ray muons at CJPL, one of the world's deepest underground laboratories.
2. Significance:
   • CJPL's 2400-meter vertical rock overburden results in muons reaching the underground facility with extremely high average energy (~360 GeV)
   • Cosmogenic neutrons constitute an important background for rare event detection such as solar neutrinos and supernova relic neutrinos
   • β-decay of cosmogenic isotopes (such as ^11C, ^9Li, etc.) is difficult to distinguish from neutrino-induced electron recoil
3. Existing Limitations:
   • Previous studies utilized only partial datasets with limited precision
   • Lack of detailed investigation into finite-size detector effects
   • Absence of neutron yield measurement data for 360 GeV high-energy muons
4. Research Motivation:
   • Provide fundamental data for low-background measurements in the JNE experiment
   • Refine empirical formulas for neutron yield from high-energy muons
   • Investigate finite-size effects in small-scale detectors

Core Contributions

1. Provided the first measurement of neutron yield in liquid scintillator at the world's highest average muon energy (360 GeV)
2. First detailed investigation of finite-size detector effects on cosmogenic product yield measurements
3. Precise measurement of CJPL muon flux using complete dataset, improving precision by approximately 25%
4. Established a complete GEANT4 simulation framework incorporating topographical features
5. Proposed finite-size correction factor fsize, providing reference for other small-scale detector experiments

Methodology Details

Experimental Apparatus Design

1-ton Prototype Detector Structure:

• Central target volume: spherical acrylic container with 645 mm diameter
• Liquid scintillator: LAB solvent + PPO fluorophore + bis-MSB wavelength shifter
• Photomultiplier tubes: 30 Hamamatsu R5912 8-inch PMTs
• Shielding system: pure water + 5 cm lead wall + stainless steel tank

Simulation Framework

Dual-layer Simulation System:

1. Topography-related Simulation:
   • Primary muon generation using MCEq and modified Gaisser formula
   • Mountain topography modeling based on NASA SRTM3 dataset
   • Delaunay triangulation method for mountain surface modeling
2. Detector-related Simulation:
   • Complete detector geometry
   • Three GEANT4 physics lists: QGSP-BIC-HP, QGSP-BERT-HP, FTFP-BERT-HP
   • Waveform saturation effect simulation

Data Analysis Methods

Muon Event Selection:

• Energy threshold: >60 MeV
• Background suppression: rmax < 0.15 (flash suppression), npeak < 40 (electronic noise suppression)
• Direction reconstruction: template-matching based time-of-arrival analysis

Neutron Event Selection:

• Muon-neutron coincidence: time window 20, 1020 μs
• Neutron candidate energy: 0, 4 MeV
• Signal extraction: two-dimensional unbinned likelihood fit (energy + time)

Efficiency and Correction Factors

Muon Detection Efficiency: ε = εg × εd + εs = 1.7%

• Geometric efficiency εg = 2.1%
• Detection efficiency εd = 82.2%
• Shower efficiency εs = 0.3‰

Neutron Yield Correction: Nn = Nobs × (1-fnon-target)×(1-fshower)×fspill / εnt×εt×fsize

Key correction factors:

• fnon-target = 0.24 (non-target correction)
• fshower = 0.21 (muon shower correction)
• fspill = 0.61 (neutron spillout correction)
• fsize = 0.63 (finite-size correction, first detailed investigation in this work)

Experimental Setup

Dataset

• Data Collection Period: July 31, 2017 to September 2, 2023
• Effective Data Acquisition Time: 1178 days
• Operating Phases: 4 different trigger condition phases
• Data Quality Control: 4.9% of problematic data excluded

Calibration System

• Gain Calibration: "RollingGain" method based on PMT dark noise
• Time Calibration: utilizing natural radioactive decay products, precision ±5 ns
• Energy Calibration: based on ^208Tl 2.61 MeV γ peak, systematic error 5%

Trigger Condition Evolution

Experimental Results

Muon Flux Measurement

Measurement Result: ϕμ = (3.56±0.16stat.±0.10syst.)×10^-10 cm^-2s^-1

Systematic Error Analysis:

• Energy calibration: ±1.6%
• Acrylic container radius: ±1.5%
• Lead shielding thickness: ±0.2%
• Total systematic error: ±2.7%

Cosmogenic Neutron Yield

Measurement Result: Yn = (3.37±1.41stat.±0.31syst.)×10^-4 μ^-1g^-1cm^2

Comparison with Theoretical Predictions:

Finite-Size Effect Investigation

Key Findings:

• Secondary particle leakage in small detectors significantly affects neutron yield measurements
• Correction factor fsize = 0.63, indicating approximately 37% of produced neutrons are missed
• Effect varies with detector size, muon energy, and target material

Particle Contribution Analysis:

• High-energy γ rays: photonuclear reactions
• Protons, neutrons, π±: inelastic scattering
• π^-: thermal capture reactions

Global Fitting Results

Updated Empirical Formula Parameters: Yn = a × Eμ^b

• a = (3.9±0.7)×10^-6 μ^-1g^-1cm^2
• b = 0.77±0.03

Related Work

Underground Laboratory Muon Measurements

This study contrasts with measurement results from the following laboratories:

• Daya Bay: 70 GeV average energy
• KamLAND: ~200 GeV
• Borexino: ~270 GeV
• Aberdeen Tunnel: ~160 GeV

History of Neutron Yield Measurements

From Hertenberger (1995) to the present study, covering the complete energy range from low to 360 GeV, providing the highest energy point constraint for empirical formulas.

Conclusions and Discussion

Main Conclusions

1. Successfully measured neutron yield at the world's highest average muon energy
2. Verified applicability of empirical formulas in the high-energy region
3. First quantification of finite-size effects in small detectors
4. Provided precise background parameters for future CJPL experiments

Limitations

1. Statistical Precision Limited: small number of neutron candidate events (13.6±5.7)
2. Model Dependence: uncertainties in physical process simulation
3. Geometric Simplification: incomplete consideration of mountain interior structures (e.g., tunnels)
4. Energy Range: applicable only at the 360 GeV energy point

Future Directions

1. Larger-scale Detectors: reduce statistical errors and finite-size effects
2. Multi-energy-point Measurements: refine empirical formulas
3. Other Cosmogenic Isotopes: yield studies of ^9Li, ^8He, etc.
4. Refined Simulation: incorporation of more complex mountain structures

In-Depth Evaluation

Strengths

1. Methodological Innovation: first systematic investigation of finite-size effects, proposing fsize correction factor
2. Experimental Completeness: utilization of complete 1178-day dataset with multiple systematic error assessments
3. Result Convincingness: consistency with multiple theoretical predictions, filling high-energy region gap
4. Technical Advancement: complete topographical simulation, precise detector response modeling

Weaknesses

1. Limited Statistical Sample: small number of neutron events affects precision
2. Model Uncertainties: reliability of GEANT4 physical processes in extreme high-energy region
3. Geometric Idealization: actual mountain structure more complex than simulation
4. Background Complexity: potential systematic bias in accidental coincidence background treatment

Impact

1. Academic Contribution: provides important background parameters for underground high-energy physics experiments
2. Practical Value: directly serves future neutrino experiments at CJPL
3. Methodological Significance: finite-size correction method generalizable to other small-scale detectors
4. International Impact: fills global data gap in high-energy muon neutron yield measurements

Applicable Scenarios

1. Underground Neutrino Experiments: background assessment for solar and supernova neutrino detection
2. Dark Matter Search: cosmogenic background estimation
3. Small Detector Design: finite-size effect correction
4. Monte Carlo Simulation Verification: high-energy physics process model validation

References

This paper cites 43 important references, covering:

• Underground laboratory muon measurements (Daya Bay, KamLAND, Borexino, etc.)
• Cosmogenic neutron theory (Mei & Hime, 2006; Li & Beacom, 2015)
• Simulation tools (GEANT4, FLUKA, MCEq)
• Neutrino experiment background studies (SNO, Super-Kamiokande, etc.)

Overall Assessment: This is a high-quality experimental physics paper with important contributions in both methodology and measurement precision. The systematic investigation of finite-size effects is particularly valuable methodologically, providing important reference for small-scale detector experiments. The measurement results establish a solid foundation for future CJPL experiments while enriching the global database of high-energy muon neutron yield data.