Simulation of Muon-induced Backgrounds for the Colorado Underground Research Institute (CURIE)

Abstract: We present a comprehensive Monte Carlo simulation of muon-induced backgrounds for the Colorado Underground Research Institute (CURIE), a shallow-underground facility with $\approx 415$~m.w.e. overburden. Using coupled \textsc{mute} and \textsc{geant4} frameworks, we characterize the production and transport of muon-induced secondaries through site-specific rock compositions and geometries, establishing a proof-of-concept for high-precision, end-to-end simulations. Our simulations employ angular-dependent muon energy distributions, which improve secondary flux accuracy. For the Subatomic Particle Hideout and Cryolab I research spaces, we predict total muon-induced neutron fluxes of $(3.78 \pm 0.61_{\text{sys}}) \times 10<sup>{-3}$~m$<sup>{-2}$s$<sup>{-1}$ and $(3.97 \pm 0.65_{\text{sys}}) \times 10<sup>{-3}$~m$<sup>{-2}$s$<sup>{-1}$, respectively, consistent with empirical depth parameterizations. The simulated neutron energy spectra exhibit the expected thermal, epithermal, evaporation, and spallation components extending to GeV energies. Electromagnetic backgrounds are expected to dominate the total flux, with $\gamma$-ray components of $(5.54 \pm 0.91_{\text{sys}}) \times 10<sup>{-1}$~m$<sup>{-2}$s$<sup>{-1}$ and $(6.51 \pm 1.06_{\text{sys}}) \times 10<sup>{-1}$~m$<sup>{-2}$s$<sup>{-1}$ for the Subatomic Particle Hideout and Cryolab I facilities, respectively. These results provide quantitative background predictions for experimental design and sensitivity projections at shallow- and deep-underground facilities. They further demonstrate that local geology and overburden geometry influence muon-induced secondary yields and energy spectra, emphasizing the need for site-specific simulations for accurate underground background characterization. Therefore, the simulation framework has been made publicly available for the broader low-background physics community to enable meaningful inter-facility comparisons.