Neutrons produced by cosmic ray muons constitute a significant background for underground experiments investigating neutron oscillations, neutrinoless double beta decay, dark matter, and other rare event signals. This work benchmarks measured neutron yields and neutron multiplicities--with a focus on data from the Daya Bay Reactor Neutrino Experiment--against comprehensive simulations using three GEANT4 hadronic physics lists. These simulations are further refined via a TALYS-based adjustment of hadronic cross sections. For the BERT-based models, the adjustment reduces the discrepancy in the total neutron yield from about 20% to approximately 6%, while for the BIC-based models it improves the agreement from roughly 13% to the sub-percent level (~0.3%), indicating a markedly better consistency of the BIC-based models with the experimental data. Nevertheless, a clear tension persists: simulations systematically underproduce single-neutron events while overproducing multi-neutron events. The study establishes a reproducible benchmark for cosmogenic neutron modeling and proposes a targeted refinement strategy--including channel-specific reweighting and intranuclear cascade parameter tuning--to guide future model development.
@article{arxiv.2601.11980,
title = {Comprehensive study of cosmogenic neutron production in large liquid scintillator detectors},
author = {Yijian Jiang and Jie Cheng and Haoqi Lu and Yaoguang Wang},
journal= {arXiv preprint arXiv:2601.11980},
year = {2026}
}