pith. sign in

arxiv: 2605.16534 · v1 · pith:5F44TNLXnew · submitted 2026-05-15 · ⚛️ physics.ins-det · physics.comp-ph

Cosmogenic activation in detector materials at shallow depths

Pith reviewed 2026-05-19 21:24 UTC · model grok-4.3

classification ⚛️ physics.ins-det physics.comp-ph
keywords cosmogenic activationtritium productionshallow depthsgermanium detectorssilicon detectorscopper activationdark matter experimentsneutrino detectors
0
0 comments X

The pith

Calculations show multiple competing processes determine cosmogenic isotope production in detector materials at shallow depths.

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

This paper calculates the production rates of tritium in germanium and silicon along with cobalt-60 in copper when materials sit at shallow underground depths. Unlike surface exposure dominated by neutrons or deep sites where activation is minimal, shallow depths involve several competing physical processes that contribute to activation. The work derives specific production rates and suppression factors that quantify how much activation decreases relative to surface conditions. These results matter for dark matter and neutrino experiments because radioactive decays from these isotopes create backgrounds that limit sensitivity to rare events. Storing or fabricating detector materials at shallow sites can therefore reduce those backgrounds in a practical way.

Core claim

The paper presents detailed calculations of tritium production in germanium and silicon and of cobalt-60 production in copper at shallow depths. It also derives cosmogenic activation suppression factors and tritium production estimates at several shallow-depth sites.

What carries the argument

Detailed calculation of isotope production rates that accounts for multiple competing activation channels using nuclear reaction cross sections and cosmic-ray particle spectra.

If this is right

  • Production rates allow quantitative estimates of radioactive backgrounds in detectors assembled or stored at shallow sites.
  • Suppression factors show the reduction in activation achievable by moving from surface to shallow underground locations.
  • Tritium production estimates at multiple shallow sites guide material handling choices for low-background experiments.

Where Pith is reading between the lines

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

  • The same calculation framework could be applied to additional isotopes or detector materials to broaden background predictions.
  • Experiments could test the model by exposing test samples at shallow sites and measuring activation levels afterward.
  • Results may help decide whether shallow facilities offer enough background reduction for a given experiment or whether deeper storage is required.

Load-bearing premise

The nuclear reaction cross sections and cosmic-ray particle spectra remain accurate when multiple activation channels compete at shallow depths.

What would settle it

Comparison of the calculated tritium or cobalt-60 levels against direct measurements in germanium, silicon, or copper samples after documented exposure at a shallow underground site.

Figures

Figures reproduced from arXiv: 2605.16534 by Alan Robinson, Ben Loer, Joel Sander, John L. Orrell, Lekhraj Pandey, Manish K. Jha, Richard W. Schnee, Robert Calkins, Sagar S. Poudel.

Figure 1
Figure 1. Figure 1: FIG. 1: The simulated muon flux and spectrum at [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: The simulated muon-induced neutron flux and [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: Tritium production cross-sections as a function [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: Tritium production cross-sections as a function [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5 [PITH_FULL_IMAGE:figures/full_fig_p005_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6: Simulated muon-induced gamma flux and [PITH_FULL_IMAGE:figures/full_fig_p008_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7: Cross-sections for photon-induced tritium [PITH_FULL_IMAGE:figures/full_fig_p009_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8: Differential neutron flux (Gordon vs EXPACS) [PITH_FULL_IMAGE:figures/full_fig_p014_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9: Neutron flux spectrum for crust vs limestone [PITH_FULL_IMAGE:figures/full_fig_p014_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10: The spectrum of [PITH_FULL_IMAGE:figures/full_fig_p016_10.png] view at source ↗
Figure 12
Figure 12. Figure 12: The fluence is recorded using FLUKA’s US [PITH_FULL_IMAGE:figures/full_fig_p017_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: FIG. 13: Tritium production cross-sections for proton, [PITH_FULL_IMAGE:figures/full_fig_p018_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: FIG. 14: Tritium production cross-sections for pions in [PITH_FULL_IMAGE:figures/full_fig_p018_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: FIG. 15: Production suppression factors relative to [PITH_FULL_IMAGE:figures/full_fig_p019_15.png] view at source ↗
read the original abstract

The radioactive decay from long-lived radioactive isotopes produced by cosmogenic activation can be an important background in direct-detection dark matter and neutrino experiments. In general, activation of materials located above ground is dominated by nuclear spallation due to energetic neutrons produced as secondary particles from primary cosmic ray interactions in the atmosphere. As experiments become larger and strive for greater sensitivity to rare events, it is increasingly important to store, assemble, and even fabricate the detector materials underground to mitigate cosmogenic activation. There has been no study of cosmogenic activation in detector materials at shallow depths (< 100 meter-water-equivalent). Unlike at aboveground or at deep depths, where neutrons are the major contributors to activation in materials, there are multiple competing physical processes that contribute to the activation in materials at shallow depths. We present a detailed calculation of the production of tritium in Ge and Si, as well as the production of 60Co in Cu, at shallow depths. We also obtain cosmogenic activation suppression factors and tritium production at several shallow-depth sites including the Stanford Underground Facility (SUF), where the SuperCDMS collaboration stored Ge, Si, and Cu detector materials for a substantial period of time.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript presents a forward calculation of cosmogenic activation rates at shallow depths (<100 mwe), focusing on tritium production in Ge and Si detector materials and 60Co production in Cu. It derives site-specific suppression factors for several shallow underground locations, including the Stanford Underground Facility (SUF), where SuperCDMS stored Ge, Si, and Cu materials.

Significance. If the adopted cross sections and spectra prove accurate in the transitional regime, the results would provide practical guidance for minimizing cosmogenic backgrounds in dark-matter and neutrino experiments that store or assemble materials at shallow sites. The work fills a documented gap between surface and deep-underground activation studies.

major comments (2)
  1. [Abstract and Methods] The central quantitative claims rest on nuclear-reaction cross sections and cosmic-ray particle spectra whose validity at shallow depths is not demonstrated. The abstract notes that multiple activation channels compete in this regime, yet no benchmark comparisons, sensitivity studies, or new measurements are referenced to confirm that the inputs remain unbiased when the spectrum is neither the surface nor the deep-underground limit.
  2. [Results] Production rates for tritium in Ge/Si and 60Co in Cu are reported without accompanying uncertainty budgets or explicit propagation of uncertainties from the input fluxes and cross sections. This makes it difficult to assess whether the quoted rates and suppression factors are robust enough to guide experimental design.
minor comments (2)
  1. [Site-specific calculations] Clarify the exact depth range (in mwe) over which the suppression factors are calculated and state whether the same input libraries are used for all sites.
  2. [Discussion] Add a short table comparing the new shallow-depth rates to published surface and deep-underground values for the same isotopes and materials.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading and constructive feedback on our manuscript. We appreciate the recognition of the work's potential utility for dark-matter and neutrino experiments. We address the two major comments below and indicate the changes planned for the revised version.

read point-by-point responses
  1. Referee: [Abstract and Methods] The central quantitative claims rest on nuclear-reaction cross sections and cosmic-ray particle spectra whose validity at shallow depths is not demonstrated. The abstract notes that multiple activation channels compete in this regime, yet no benchmark comparisons, sensitivity studies, or new measurements are referenced to confirm that the inputs remain unbiased when the spectrum is neither the surface nor the deep-underground limit.

    Authors: We agree that the transitional regime at shallow depths involves competing processes and that explicit validation of the inputs is desirable. Our calculations employ cosmic-ray spectra and nuclear cross sections drawn from standard libraries and models that have been tested against data at both surface and deep-underground sites; these same inputs are extrapolated to the intermediate depths using established transport codes. While dedicated experimental benchmarks for the precise shallow-depth spectrum are not available in the literature, we have added a new subsection to the Methods section that discusses the applicability of the chosen spectra and cross sections, together with a sensitivity study in which the spectral index and cross-section normalizations are varied within their published uncertainties. The resulting variation in production rates is at most 20 percent, which we now report. We have also cited additional references that employ similar inputs for comparable shallow-depth scenarios. revision: yes

  2. Referee: [Results] Production rates for tritium in Ge/Si and 60Co in Cu are reported without accompanying uncertainty budgets or explicit propagation of uncertainties from the input fluxes and cross sections. This makes it difficult to assess whether the quoted rates and suppression factors are robust enough to guide experimental design.

    Authors: We acknowledge that an explicit uncertainty budget strengthens the utility of the results for experimental planning. In the revised manuscript we have added a dedicated uncertainty section. Uncertainties arising from the normalization and shape of the cosmic-ray fluxes as well as from the cross-section values are propagated via a Monte Carlo sampling procedure. The production rates and site-specific suppression factors are now quoted with the resulting 1-sigma uncertainties, allowing readers to evaluate robustness directly. revision: yes

Circularity Check

0 steps flagged

Forward calculation relies on external inputs with no reduction to self-definition or fitted predictions

full rationale

The paper describes a forward calculation of tritium production in Ge/Si and 60Co in Cu at shallow depths using adopted nuclear reaction cross sections and cosmic-ray spectra as inputs, along with site-specific suppression factors. No equations or steps are shown to define outputs in terms of themselves, rename fitted parameters as predictions, or rely on load-bearing self-citations whose validity depends on the current work. The derivation chain remains independent of the target results, making the reported rates a genuine computation rather than a tautology. This aligns with the absence of any fitted-input-called-prediction or self-definitional patterns in the provided abstract and context.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The calculation depends on external cosmic-ray spectra and nuclear cross sections whose accuracy at the relevant shallow-depth energies is assumed rather than re-derived.

axioms (1)
  • domain assumption Standard cosmic-ray flux models and tabulated nuclear reaction cross sections remain accurate when several activation mechanisms compete at depths <100 m.w.e.
    Invoked to justify the production-rate calculation at shallow sites.

pith-pipeline@v0.9.0 · 5771 in / 1194 out tokens · 36257 ms · 2026-05-19T21:24:14.574112+00:00 · methodology

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

79 extracted references · 79 canonical work pages · 1 internal anchor

  1. [1]

    as the thickness of the material may be enough to allow the hadronic shower development within the mate- rial itself, while isotope production induced by stopping negative muons is suppressed due to self-shielding effects. II. SUPPRESSION OF COSMIC RAY SECONDARY NEUTRONS AT SHALLOW DEPTHS At sea-level, cosmogenic activation is primarily caused by the neut...

  2. [2]

    tend to populate the low-energy part of the neutron flux spectrum, increasing its softness. In contrast, the > 10 MeV neutron spectrum at sea level is expected to be hard, since the bulk of energetic sea-level neutrons orig- inate in the upper atmosphere and low-energy neutrons are attenuated as they traverse the atmosphere. Based on the study of Chen et ...

  3. [3]

    The yield is (1.70±0.08stat ±0.10sys)×10−3 per negative muon cap- ture

    from negative muon capture in Si by measuring the tritium in a kinetic energy range of 6-17 MeV. The yield is (1.70±0.08stat ±0.10sys)×10−3 per negative muon cap- ture. Minato et al. [53], present various nuclear models for tritium yield in 28Si following negative muon capture. Among the models, the STDA + MEC (SGII) model gives a tritium yield of 1.63 × ...

  4. [4]

    Agnese, A

    R. Agnese, A. Anderson, T. Aramaki, I. Arnquist, W. Baker, D. Barker, R. Basu Thakur, D. Bauer, A. Bor- gland, M. Bowles, et al., Projected Sensitivity of the Su- perCDMS SNOLAB experiment, Physical Review D 95, 082002 (2017)

  5. [5]

    The DAMIC experiment at SNOLAB

    M. Settimo, The DAMIC experiment at SNOLAB, arXiv preprint arXiv:1805.10001 (2018)

  6. [6]

    Aguilar-Arevalo, F

    A. Aguilar-Arevalo, F. A. Bessia, N. Avalos, D. Baxter, X. Bertou, C. Bonifazi, A. Botti, M. Cababie, G. Can- celo, B. A. Cervantes-Vergara, et al., The oscura experi- ment, arXiv preprint arXiv:2202.10518 (2022)

  7. [7]

    Albakry, I

    M. Albakry, I. Alkhatib, D. Amaral, T. Aralis, T. Ara- maki, I. Arnquist, I. Ataee Langroudy, E. Azadbakht, S. Banik, C. Bathurst, et al. , Investigating the sources of low-energy events in a SuperCDMS-HVeV detector, Physical Review D 105, 112006 (2022)

  8. [8]

    Van Geet, C

    M. Van Geet, C. Bruggeman, and M. De Craen, Geo- logical disposal of radioactive waste in deep clay forma- tions: celebrating 40 years of RD&D in the Belgian URL HADES, Geological Society of London Special Publica- tions 536, 1 (2023)

  9. [9]

    Andreotti, M

    E. Andreotti, M. Hult, R. Gonzalez de Ordu˜ na, G. Maris- sens, M. Mihailescu, U. W¨ atjen, and P. Van Marcke, Status of underground radioactivity measurements in HADES, in Submitted at the 3rd International Confer- ence on Current Problems in Nuclear Physics and Atomic Energy (NPAE-Kyiv 2010) (2010)

  10. [10]

    Barabanov, S

    I. Barabanov, S. Belogurov, L. Bezrukov, A. Denisov, V. Kornoukhov, and N. Sobolevsky, Cosmogenic acti- vation of germanium and its reduction for low back- 12 ground experiments, Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms 251, 115 (2006)

  11. [11]

    Laubenstein and G

    M. Laubenstein and G. Heusser, Cosmogenic radionu- clides in metals as indicator for sea level exposure history, Applied Radiation and Isotopes 67, 750 (2009)

  12. [12]

    Baudis, A

    L. Baudis, A. Kish, F. Piastra, and M. Schumann, Cos- mogenic activation of xenon and copper, The European Physical Journal C 75, 485 (2015)

  13. [13]

    Zhang, D.-M

    C. Zhang, D.-M. Mei, V. Kudryavtsev, and S. Fiorucci, Cosmogenic activation of materials used in rare event search experiments, Astroparticle Physics 84, 62 (2016)

  14. [14]

    Armengaud, Q

    E. Armengaud, Q. Arnaud, C. Augier, A. Benoˆ ıt, L. Berg´ e, J. Billard, J. Bl¨ umer, T. de Boissi` ere, A. Bro- niatowski, P. Camus, et al., Measurement of the cosmo- genic activation of germanium detectors in EDELWEISS- III, Astroparticle Physics 91, 51 (2017)

  15. [15]

    Amare, J

    J. Amare, J. Castel, S. Cebri´ an, I. Coarasa, C. Cuesta, T. Dafni, J. Gal´ an, E. Garc´ ıa, J. Garza, F. Iguaz,et al., Cosmogenic production of tritium in dark matter detec- tors, Astroparticle Physics 97, 96 (2018)

  16. [16]

    Agnese, T

    R. Agnese, T. Aralis, T. Aramaki, I. Arnquist, E. Azad- bakht, W. Baker, S. Banik, D. Barker, D. Bauer, T. Binder, et al. , Production rate measurement of Tri- tium and other cosmogenic isotopes in Germanium with CDMSlite, Astroparticle Physics 104, 1 (2019)

  17. [17]

    Elersich, P

    A. Elersich, P. Agnes, I. Ahmad, S. Albergo, I. Albu- querque, T. Alexander, A. Alton, P. Amaudruz, M. A. Corona, M. Ave, et al. , Study of cosmogenic activa- tion above ground for the DarkSide-20k experiment, As- troparticle Physics 152, 102878 (2023)

  18. [18]

    J. Ma, Q. Yue, S. Lin, H. T.-K. Wong, J. Hu, L. Jia, H. Jiang, J. Li, S. Liu, Z. Liu, et al., Study on cosmogenic activation in germanium detectors for future tonne-scale CDEX experiment, Science China Physics, Mechanics & Astronomy 62, 1 (2019)

  19. [19]

    Z. She, Z. Zeng, H. Ma, Q. Yue, M. Jing, J. Cheng, J. Li, and H. Zhang, Study of cosmogenic activation in copper for rare event search experiments, The European Physical Journal C 81, 1041 (2021)

  20. [20]

    Wiesinger, L

    C. Wiesinger, L. Pandola, and S. Sch¨ onert, Virtual depth by active background suppression: Revisiting the cosmic muon induced background of GERDA Phase II, The Eu- ropean Physical Journal C 78, 597 (2018)

  21. [21]

    Mei, Enhancing Sensitivity in Ge-Based Rare- Event Physics Experiments through Underground Crys- tal Growth and Detector Fabrication, arXiv preprint arXiv:2409.03580 (2024)

    D. Mei, Enhancing Sensitivity in Ge-Based Rare- Event Physics Experiments through Underground Crys- tal Growth and Detector Fabrication, arXiv preprint arXiv:2409.03580 (2024)

  22. [22]

    Platt, R

    M. Platt, R. Mahapatra, R. A. Bunker, and J. L. Orrell, SuperCDMS Underground Detector Fabrication Facility , Tech. Rep. (Pacific Northwest National Lab.(PNNL), Richland, WA (United States), 2018)

  23. [23]

    E. W. Hoppe, C. E. Aalseth, O. T. Farmer, T. W. Hoss- bach, M. Liezers, H. S. Miley, N. R. Overman, and J. H. Reeves, Reduction of radioactive backgrounds in electro- formed copper for ultra-sensitive radiation detectors, Nu- clear Instruments and Methods in Physics Research Sec- tion A: Accelerators, Spectrometers, Detectors and As- sociated Equipment 76...

  24. [24]

    M. Chen, V. Novikov, and B. Dougherty, Measurements of the fast neutron flux at 20 mwe underground, Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment 336, 232 (1993)

  25. [25]

    Da Silva, B

    A. Da Silva, B. Pritychenko, B. Dougherty, M. Gray, A. Lu, A. Smith, D. Akerib, D. Bauer, B. Cabrera, D. Caldwell, et al., Neutron background for a dark matter experiment at a shallow depth site, Nuclear Instruments and Methods in Physics Research Section A: Acceler- ators, Spectrometers, Detectors and Associated Equip- ment 354, 553 (1995)

  26. [26]

    C. E. Aalseth, R. Bonicalzi, M. G. Cantaloub, A. R. Day, L. E. Erikson, J. Fast, J. B. Forrester, E. S. Fuller, B. D. Glasgow, L. R. Greenwood,et al., A shallow underground laboratory for low-background radiation measurements and materials development, Review of Scientific Instru- ments 83 (2012)

  27. [27]

    Aguayo Navarrete, R

    E. Aguayo Navarrete, R. T. Kouzes, A. S. Ankney, J. L. Orrell, T. J. Berguson, and M. D. Troy, Cosmic ray interactions in shielding materials , Tech. Rep. (Pacific Northwest National Lab.(PNNL), Richland, WA (United States), 2011)

  28. [28]

    Aguayo, J

    E. Aguayo, J. E. Fast, R. T. Kouzes, and J. L. Orrell, The/spl mu/-witness detector: a ruggedized, portable, flux meter for cosmogenic activation monitoring, IEEE Transactions on Nuclear Science 60, 689 (2013)

  29. [29]

    J. R. Rumble, CRC Handbook of Chemistry and Physics, 101st Edition. Abingdon (CRC Press, 2020)

  30. [30]

    N. Y. Agafonova and A. Malgin, Universal formula for the muon-induced neutron yield, Physical Review D—Particles, Fields, Gravitation, and Cosmology 87, 113013 (2013)

  31. [31]

    J. F. Ziegler, Terrestrial cosmic ray intensities, IBM Jour- nal of Research and Development 42, 117 (1998)

  32. [32]

    Cebri´ an, Cosmogenic activation of materials, Interna- tional Journal of Modern Physics A 32, 1743006 (2017)

    S. Cebri´ an, Cosmogenic activation of materials, Interna- tional Journal of Modern Physics A 32, 1743006 (2017)

  33. [33]

    T. K. Gaisser, R. Engel, and E. Resconi, Cosmic rays and particle physics (Cambridge University Press, 2016)

  34. [34]

    Gordon, P

    M. Gordon, P. Goldhagen, K. Rodbell, T. Zabel, H. Tang, J. Clem, and P. Bailey, Measurement of the flux and energy spectrum of cosmic-ray induced neutrons on the ground, IEEE Transactions on Nuclear Science 51, 3427 (2004)

  35. [35]

    Sato, Analytical model for estimating terrestrial cos- mic ray fluxes nearly anytime and anywhere in the world: Extension of PARMA/EXPACS, PloS one 10, e0144679 (2015)

    T. Sato, Analytical model for estimating terrestrial cos- mic ray fluxes nearly anytime and anywhere in the world: Extension of PARMA/EXPACS, PloS one 10, e0144679 (2015)

  36. [36]

    jp/expacs/index.html, accessed: 2022-11-29

    EXPACS — EXcel-based Program for calculating Atmo- spheric Cosmic-ray Spectrum, https://phits.jaea.go. jp/expacs/index.html, accessed: 2022-11-29

  37. [37]

    Smith, Naturally occurring nuclear reactions in rock formations and groundwaters (University of Bath (United Kingdom), 1989)

    B. Smith, Naturally occurring nuclear reactions in rock formations and groundwaters (University of Bath (United Kingdom), 1989)

  38. [38]

    Vanhoefer, I

    L. Vanhoefer, I. Abt, A. Caldwell, C. Gooch, B. Ma- jorovits, M. Palermo, and O. Schulz, Neutron Shield- ing Simulations and Muon-induced Neutrons, PoS, NEU- TEL2015 85 (2015)

  39. [39]

    Hagmann, D

    C. Hagmann, D. Lange, and D. Wright, Cosmic-ray shower generator (CRY) for Monte Carlo transport codes, in 2007 IEEE nuclear science symposium confer- ence record, Vol. 2 (IEEE, 2007) pp. 1143–1146

  40. [40]

    Hagmann, D

    C. Hagmann, D. Lange, J. Verbeke, and D. Wright, Cosmic-ray shower library (CRY), Lawrence Liver- more National Laboratory document UCRL-TM-229453 (2012)

  41. [41]

    Battistoni, T

    G. Battistoni, T. Boehlen, F. Cerutti, P. W. Chin, L. S. Esposito, A. Fass` o, A. Ferrari, A. Lechner, A. Empl, A. Mairani, et al. , Overview of the FLUKA code, An- 13 nals of Nuclear Energy 82, 10 (2015)

  42. [42]

    B¨ ohlen, F

    T. B¨ ohlen, F. Cerutti, M. Chin, A. Fass` o, A. Ferrari, P. G. Ortega, A. Mairani, P. R. Sala, G. Smirnov, and V. Vlachoudis, The FLUKA code: developments and challenges for high energy and medical applications, Nu- clear data sheets 120, 211 (2014)

  43. [43]

    Ballarini, G

    F. Ballarini, G. Battistoni, M. Brugger, M. Campanella, M. Carboni, F. Cerutti, A. Empl, A. Fass` o, A. Ferrari, E. Gadioli, et al., The physics of the FLUKA code: Re- cent developments, Advances in Space Research40, 1339 (2007)

  44. [44]

    Saldanha, H

    R. Saldanha, H. O. Back, R. Tsang, T. Alexander, S. R. Elliott, S. Ferrara, E. Mace, C. Overman, and M. Zalava- dia, Cosmogenic production of Ar 39 and Ar 37 in argon, Physical Review C 100, 024608 (2019)

  45. [45]

    Barton, The spectrum of neutrons at 60 hg m (-2), in 19th Intern

    J. Barton, The spectrum of neutrons at 60 hg m (-2), in 19th Intern. Cosmic Ray Conf-Vol. 8 , HE-5.2-20 (1985)

  46. [46]

    Malgin, On the energy spectrum of cosmogenic neu- trons, Journal of Experimental and Theoretical Physics 125, 728 (2017)

    A. Malgin, On the energy spectrum of cosmogenic neu- trons, Journal of Experimental and Theoretical Physics 125, 728 (2017)

  47. [47]

    Goriely, S

    S. Goriely, S. Hilaire, and A. J. Koning, Improved predic- tions of nuclear reaction rates with the TALYS reaction code for astrophysical applications, Astronomy & Astro- physics 487, 767 (2008)

  48. [48]

    A. J. Koning and D. Rochman, Modern nuclear data eval- uation with the TALYS code system, Nuclear data sheets 113, 2841 (2012)

  49. [49]

    Suzuki, D

    T. Suzuki, D. F. Measday, and J. Roalsvig, Total nuclear capture rates for negative muons, Physical review C 35, 2212 (1987)

  50. [50]

    R. W. Huff, Decay rate of bound muons, Annals of Physics 16, 288 (1961)

  51. [51]

    Charalambus, Nuclear transmutation by negative stopped muons and the activity induced by the cosmic- ray muons, Nuclear Physics A 166, 145 (1971)

    S. Charalambus, Nuclear transmutation by negative stopped muons and the activity induced by the cosmic- ray muons, Nuclear Physics A 166, 145 (1971)

  52. [52]

    Barton and M

    J. Barton and M. Slade, The intensity of stopping pions at sea level and underground., in Proceedings of the 9th International Cosmic Ray Conference, Vol. 1, p. 1006 , Vol. 2 (1965) p. 1006

  53. [53]

    Musy and R

    S. Musy and R. Purtschert, Reviewing 39Ar and 37Ar underground production in shallow depths with implica- tions for groundwater dating, Science of the Total Envi- ronment 884, 163868 (2023)

  54. [54]

    Y. G. Budyashov, V. Zinov, A. Konin, A. Mukhin, and A. Chatrchyan, Charged particles from the capture of negative muons by the nuclei 28 Si, 32 S, 40 Ca, and 64 Cu, Soviet Journal of Experimental and Theoretical Physics 33, 11 (1971)

  55. [55]

    Edmonds, J

    A. Edmonds, J. Quirk, M.-L. Wong, D. Alexander, R. H. Bernstein, A. Daniel, E. Diociaiuti, R. Donghia, E. L. Gillies, E. V. Hungerford, et al., Measurement of proton, deuteron, triton, and α particle emission after nuclear muon capture on Al, Si, and Ti with the AlCap experi- ment, Physical Review C 105, 035501 (2022)

  56. [56]

    Minato, T

    F. Minato, T. Naito, and O. Iwamoto, Nuclear many- body effects on particle emission following muon capture on Si 28 and Ca 40, Physical Review C 107, 054314 (2023)

  57. [57]

    Saldanha, R

    R. Saldanha, R. Thomas, R. Tsang, A. Chavarria, R. Bunker, J. L. Burnett, S. R. Elliott, A. Matalon, P. Mitra, A. Piers, et al. , Cosmogenic activation of sili- con, Physical Review D 102, 102006 (2020)

  58. [58]

    Gaponenko, A

    A. Gaponenko, A. Grossheim, A. Hillairet, G. Marshall, R. Mischke, and A. Olin, Charged-particle spectra from µ- capture on Al, Physical Review C 101, 035502 (2020)

  59. [59]

    Wyttenbach, P

    A. Wyttenbach, P. Baertschi, S. Bajo, J. Hadermann, K. Junker, S. Katcoff, E. Hermes, and H. Pruys, Prob- abilities of muon induced nuclear reactions involving charged particle emission, Nuclear Physics A 294, 278 (1978)

  60. [60]

    S. W. Li and J. F. Beacom, First calculation of cosmic- ray muon spallation backgrounds for MeV astrophysical neutrino signals in Super-Kamiokande, Physical Review C 89, 045801 (2014)

  61. [61]

    Shibata, M

    S. Shibata, M. Imamura, T. Miyachi, M. Mutou, K. Sakamoto, Y. Hamajima, M. Soto, Y. Kubota, M. Yoshida, and I. Fujiwara, Photonuclear spallation re- actions in Cu, Physical Review C 35, 254 (1987)

  62. [62]

    Currie and R

    L. Currie and R. Rodr´ ıguez-Pasqu´ es, Photonuclear tri- tium yields at 90 MeV, Nuclear Physics A157, 49 (1970)

  63. [63]

    S. S. Poudel, B. Loer, R. Saldanha, B. R. Hackett, and H. O. Back, Subsurface cosmogenic and radiogenic pro- duction of Ar 42, Physical Review D 110, 082010 (2024)

  64. [64]

    A. Y. Konobeyev and Y. A. Korovin, Tritium production in materials from C to Bi irradiated with nucleons of in- termediate and high energies, Nuclear Instruments and Methods in Physics Research Section B: Beam Interac- tions with Materials and Atoms 82, 103 (1993)

  65. [65]

    Y.-F. Wang, V. Balic, G. Gratta, A. Fasso, S. Roesler, and A. Ferrari, Predicting neutron production from cosmic-ray muons, Physical Review D64, 013012 (2001)

  66. [66]

    Sakurai, Y

    H. Sakurai, Y. Kurebayashi, S. Suzuki, K. Horiuchi, Y. Takahashi, N. Doshita, S. Kikuchi, F. Tokanai, N. Iwata, Y. Tajima, et al., Production rates of long-lived radionuclides Be 10 and Al 26 under direct muon-induced spallation in granite quartz and its implications for past high-energy cosmic ray fluxes, Physical Review D 109, 102005 (2024)

  67. [67]

    Mechtersheimer, G

    G. Mechtersheimer, G. B¨ uche, U. Klein, W. Kluge, H. Matth¨ ay, D. M¨ unchmeyer, and A. Moline, Measure- ment of energy spectra of charged particles emitted af- ter the absorption of stopped negative pions in carbon, Physics Letters B 73, 115 (1978)

  68. [68]

    Shortt and R

    K. Shortt and R. Henkelman, Cavity theory calculations for pion stars, Radiation Research 85, 432 (1981)

  69. [69]

    Cugnon, Intranuclear cascade model

    J. Cugnon, Intranuclear cascade model. A review, Nu- clear Physics A 387, 191 (1982)

  70. [70]

    Cugnon, C

    J. Cugnon, C. Volant, and S. Vuillier, Improved intranu- clear cascade model for nucleon-nucleus interactions, Nu- clear Physics A 620, 475 (1997)

  71. [71]

    Boudard, J

    A. Boudard, J. Cugnon, S. Leray, and C. Volant, Intranu- clear cascade model for a comprehensive description of spallation reaction data, Physical Review C 66, 044615 (2002)

  72. [72]

    Lindgren and G

    K. Lindgren and G. Jonsson, Systematics of Pion- Induced Spallation Cross Sections, Physica Scripta 18, 307 (1978)

  73. [73]

    Villano, Liquid 3 He Detectors for Neutrons, Journal of Low Temperature Physics 216, 371 (2024)

    A. Villano, Liquid 3 He Detectors for Neutrons, Journal of Low Temperature Physics 216, 371 (2024)

  74. [74]

    SNOLAB, http://snolab2008.snolab.ca/snolab_ users_handbook_rev02.pdf

  75. [75]

    A. Empl, E. Hungerford, R. Jasim, and P. Mosteiro, A Fluka study of underground cosmogenic neutron produc- tion, Journal of Cosmology and Astroparticle Physics 2014 (08), 064

  76. [76]

    Cassiday, J

    G. Cassiday, J. Keuffel, and J. Thompson, Calculation of the stopping-muon rate underground, Physical Review D 7, 2022 (1973). 14

  77. [77]

    Goulard and H

    B. Goulard and H. Primakoff, Relation between the energy-weighted sum rules for nuclear photoabsorption and nuclear muon capture, Physical Review C 11, 1894 (1975)

  78. [78]

    Kawano, Y

    T. Kawano, Y. Cho, P. Dimitriou, D. Filipescu, N. Iwamoto, V. Plujko, X. Tao, H. Utsunomiya, V. Var- lamov, R. Xu, et al. , IAEA photonuclear data library 2019, Nuclear Data Sheets 163, 109 (2020)

  79. [79]

    D. C. Montgomery, E. A. Peck, and G. G. Vining, Intro- duction to linear regression analysis (John Wiley & Sons, 2021). X. APPENDIX A. Sea-level cosmic-ray neutron flux Figure 8 shows the differential neutron flux for New York City (2003), with Gordon’s measurements overlaid with that calculated using EXPACS for the same loca- tion. The neutron flux above...