---
title: Colorado Underground Research Institute (CURIE)
url: https://www.emergentmind.com/topics/colorado-underground-research-institute-curie
type: topic
---

# Colorado Underground Research Institute (CURIE)

The Colorado Underground Research Institute (CURIE) is a shallow-underground low-background laboratory being developed inside the Edgar Experimental Mine (EEM) in Idaho Springs, Colorado, and operated by the Colorado School of Mines. Its defining technical characteristic is an effective shielding of approximately \(0.415 \pm 0.027\ \mathrm{km.w.e.}\), corresponding to a measured underground muon flux of order \(10^{-1}\ \mathrm{m^{-2}\ s^{-1}}\) and a reduction of the sea-level muon flux by roughly \(700\times\). Within the current literature, CURIE is presented as an accessible, horizontally entered underground environment for detector R\&D, calibration, environmental characterization, and selected low-background physics experiments that can tolerate moderate residual cosmogenic rates, while also serving as a test case for shallow-site muon and muon-induced background modeling [2411.01626; 2510.06150].

## 1. Site, access, and physical configuration

CURIE is located in the EEM near Idaho Springs, Colorado, with the EEM entrance near longitude \(-105.52\) and latitude \(39.74\), at an altitude of \(2401\ \mathrm{m}\) above sea level. The mine has two main tunnels, the Army tunnel and the Miami tunnel, and the Miami tunnel provides the direct, shortest horizontal access to the CURIE sites. Sites 0, 1, and 2 are roughly 400 horizontal meters into the Miami tunnel and are connected internally to the Army tunnel. The vertical height to the surface is approximately 200 m for all sites [2411.01626].

A central operational feature of CURIE is horizontal tunnel access. The published description emphasizes that this enables straightforward logistics, transport, and installation compared with vertical shafts typical of some deep laboratories. Tunnel openings vary from roughly \(1.8\ \mathrm{m} \times 1.8\ \mathrm{m}\) at the narrowest to \(4.5\ \mathrm{m} \times 4.5\ \mathrm{m}\) at the widest, while the Army tunnel has a narrowest access of approximately \(2.4\ \mathrm{m} \times 2.4\ \mathrm{m}\). This geometry is directly relevant to experimental handling, cryogenic delivery, shielding deployment, and staged installation [2411.01626].

The mine is excavated in gneiss, described as analogous to granite but foliated. The reported consequence is excellent rock stability and safety, allowing flexible siting. EEM infrastructure includes 110 V single-phase and 440 V three-phase power, extensive ventilation, compressed air, and water. Site 2 is specifically described as tailored for a cryogenic low-background facility aimed at developing high-resolution, low-threshold detectors and enabling controlled studies of ionizing radiation, electromagnetic disturbances, and vibrations on superconducting circuits, including qubits and cryo-CMOS [2411.01626].

By 2025, two research spaces had become the focus of detailed secondary-background simulations: the Subatomic Particle Hideout (SPH), approximately \(5.2\ \mathrm{m} \times 3.0\ \mathrm{m} \times 3.0\ \mathrm{m}\), brought online in 2024, and Cryolab I (CLI), approximately \(12.2\ \mathrm{m} \times 4.9\ \mathrm{m} \times 4.6\ \mathrm{m}\), under development. The literature characterizes these spaces as representative cavern environments within a non-flat mountain overburden where local geometry and rock composition materially affect underground spectra and source terms [2510.06150].

## 2. Direct muon-flux measurements and overburden characterization

The first quantitative characterization of CURIE’s shielding environment was based on underground muon-flux measurements at Site 0 and Site 1. The detector system used two stacked plastic scintillators with a lead sheet between them, operated in coincidence to suppress gamma backgrounds. The top scintillator \(S1\) measured \(12.5\ \mathrm{cm} \times 14.5\ \mathrm{cm} \times 4.0\ \mathrm{cm}\), the bottom scintillator \(S2\) \(30.0\ \mathrm{cm} \times 48.0\ \mathrm{cm} \times 4.0\ \mathrm{cm}\), and the lead sheet \(18.0\ \mathrm{cm} \times 30.0\ \mathrm{cm} \times 2.0\ \mathrm{cm}\), with each dimension measured to \(\pm 0.5\ \mathrm{cm}\). Each scintillator was coupled to a PMT operated at 800 V bias. The coincidence time window was approximately \(4\ \mu\mathrm{s}\), and the discriminator threshold on the larger panel was approximately 2 MeV [2411.01626].

The effective acceptance was determined with GEANT4 geometric Monte Carlo using geantinos and angular weighting from MUTE’s underground intensity map, yielding
\[
A_{\mathrm{eff}} = 241.15 \pm 23_{\mathrm{sys}} \pm 3.24_{\mathrm{stat}}\ \mathrm{cm^2}.
\]
Efficiency calibration used single-mode operation of the top panel and the relation \(\epsilon = R_{\mathrm{coin}}/R_{\mathrm{single}}\), giving \(\epsilon = 0.823 \pm 0.039\) at Site 0 and \(\epsilon = 0.799 \pm 0.035\) at Site 1. Live times were \(385.5\ \mathrm{h}\) at Site 0 and \(619.7\ \mathrm{h}\) at Site 1 [2411.01626].

Muon and gamma contributions were separated by converting channel number to deposited energy using simulated minimum-ionizing muon energy loss in plastic, \(2\ \mathrm{MeV/cm}\), and fitting the gamma background with
\[
y = a e^{-bx} + c,
\]
where the constant term preserves the muon pedestal below the cutoff and only the exponential term is subtracted. The adjusted total muon counts were \(6786 \pm 60\) for Site 0 and \(10302 \pm 230\) for Site 1. Fluxes were then computed as
\[
\phi_\mu = \frac{N_{\mathrm{tot}}}{T_{\mathrm{live}} \times A_{\mathrm{gen}} \times \epsilon},
\]
giving measured underground muon fluxes of \(0.246 \pm 0.020_{\mathrm{sys}} \pm 0.012_{\mathrm{stat}}\ \mu/\mathrm{m^2/s}\) at Site 0 and \(0.239 \pm 0.025_{\mathrm{sys}} \pm 0.010_{\mathrm{stat}}\ \mu/\mathrm{m^2/s}\) at Site 1 [2411.01626].

These measurements correspond to an approximately \(700\times\) reduction relative to a sea-level reference flux of about \(166.7\ \mathrm{m^{-2}\ s^{-1}}\). The same work reports an average equivalent vertical depth of \(0.415 \pm 0.027\ \mathrm{km.w.e.}\), placing CURIE in the shallow-underground regime. MUTE predicts nearly identical underground muon energy spectra across sites, with mean underground muon energy \(\langle E_\mu\rangle = 98 \pm 2\ \mathrm{GeV}\) [2411.01626].

A notable feature of the measurements is strong angular dependence induced by mountain geometry. At Site 0, a directional comparison between a South-facing orientation \((\theta = 34^\circ, \phi = 270^\circ)\) and a Northwest-facing orientation \((\theta = 34^\circ, \phi = 135^\circ)\) found a measured flux reduction of \(22.8 \pm 2.1\%\), while GEANT4 plus MUTE predicted \(24.0 \pm 2.3\%\). This is direct evidence that topographic anisotropy is experimentally resolvable at CURIE and cannot be collapsed to a single scalar depth without some loss of information [2411.01626].

## 3. Computational pipeline and flat-depth normalization

CURIE has also served as a validation platform for a shallow-site simulation chain combining Daemonflux, MUTE, PROPOSAL, and GEANT4. In the muon-flux study, Daemonflux v0.8.1 generated surface muon fluxes using MCEq cascade equations calibrated to the Global Spline Fit and Data-Driven Model inclusive hadronic data with SIBYLL-2.3c; MUTE v2.0.1 propagated muons through rock using PROPOSAL; and GEANT4 handled detector acceptance and effective-area calculations. The overburden was constructed from USGS DEM data over approximately \(84\ \mathrm{km^2}\) around EEM, with DEM accuracy of roughly 1 m in \(x\)–\(y\) and \(13.6\ \mathrm{cm}\) vertically [2411.01626].

Slant depth was converted to water-equivalent depth using
\[
X(\theta,\phi) = R(\theta,\phi)\times \left(\frac{\rho_r}{\rho_w}\right),
\]
with \(\rho_w = 0.997\ \mathrm{g/cm^3}\) at \(25^\circ\mathrm{C}\). Density models included standard rock \(\rho_{\mathrm{sr}} = 2.650\ \mathrm{g/cm^3}\), a weighted average \(\rho_{\mathrm{avg}} = 2.768\ \mathrm{g/cm^3}\), and quadrant-dependent averages \(Q1 = 2.765\), \(Q2 = 2.802\), \(Q3 = 2.806\), and \(Q4 = 2.752\ \mathrm{g/cm^3}\). Because differences among \(\rho_{\mathrm{avg}}\) and \(\rho_{\mathrm{quad}}\) were below 1% and the difference from standard rock below 3%, the final simulations adopted \(\rho_{\mathrm{avg}}\) and assigned a 3% systematic to inhomogeneities, water saturation, and voids [2411.01626].

Air gaps in mountain trajectories were corrected using a path-length algorithm that scans repeated \(R\) values along \((\theta,\phi)\) columns over the full \(2\pi\) steradian. At \(1^\circ\) angular sampling, the correction increased the average flux by 1–2%, and a \(+1.5\%\) correction was adopted. With full-scale simulations of approximately 30 million propagated muons, the predicted total underground muon fluxes were \(0.227 \pm 0.023\ \mathrm{m^{-2}\ s^{-1}}\) for Site 0, \(0.217 \pm 0.022\ \mathrm{m^{-2}\ s^{-1}}\) for Site 1, and \(0.259 \pm 0.026\ \mathrm{m^{-2}\ s^{-1}}\) for Site 2, in agreement with measured values at roughly the \(1\sigma\) level [2411.01626].

A distinct contribution of the 2024 work is a new depth–intensity relationship for converting non-flat mountain overburden to equivalent flat depth:
\[
\Phi_\mu(H) = A \times \left(\frac{H_0}{H}\right)^n \times e^{-H/H_0},
\]
with \(A = 0.044\ \mu/\mathrm{m^2/s}\), \(n = 2.145\), and \(H_0 = 1.094\ \mathrm{km.w.e.}\). Using this relation, the measured equivalent depths were \(0.413 \pm 0.018\ \mathrm{km.w.e.}\) for Site 0 and \(0.418 \pm 0.020\ \mathrm{km.w.e.}\) for Site 1; the simulated values were \(0.426 \pm 0.014\ \mathrm{km.w.e.}\), \(0.434 \pm 0.014\ \mathrm{km.w.e.}\), and \(0.405 \pm 0.013\ \mathrm{km.w.e.}\) for Sites 0, 1, and 2, respectively. For the measured fluxes, Mitrica et al. would yield \(0.546\)–\(0.552\ \mathrm{km.w.e.}\), while Mei–Hime would yield \(0.306\)–\(0.314\ \mathrm{km.w.e.}\), so the new relation is presented as better aligned with world data at shallow depths [2411.01626].

The same study identifies this as the first successful validation of the Daemonflux+MUTE+PROPOSAL framework against measurements at depths below \(1\ \mathrm{km.w.e.}\). In practice, this establishes CURIE not only as an experimental site but also as a benchmark environment for shallow-overburden transport and normalization methods [2411.01626].

## 4. Muon-induced secondary backgrounds

A subsequent end-to-end study modeled muon-induced secondaries for SPH and CLI with a two-stage mute–GEANT4 framework. Stage 1 used mute v2.0.1 to propagate atmospheric muons through the mountain profile and produce angle-dependent underground differential muon fluxes \(\Phi_\mu(E_\mu,\theta,\phi)\). Stage 2 used Geant4 v11-02-01 to transport these muons through a simplified, site-specific rock-and-cavern geometry and to simulate production and transport of electromagnetic, hadronic, and neutrino secondaries. The laboratory was embedded in a homogeneous hemispherical rock volume with a minimum 4 m thickness from hemisphere boundary to cavern wall to ensure shower equilibrium, and secondaries were tracked to their first entry through the rock–cavern boundary [2510.06150].

The rock models were built from SEM analyses of local samples. For SPH, the density was \(2.69 \pm 0.06\ \mathrm{g\ cm^{-3}}\) with mass fractions, in percent, of Si 36.22, O 49.84, Al 5.89, K 1.65, Na 3.47, Mg 1.08, Fe 1.24, and Other 0.61; \(\langle Z\rangle = 11.06\), \(\langle A\rangle = 22.27\), and \(\langle Z^2/A\rangle = 5.49\). For CLI, the density was \(2.71 \pm 0.07\ \mathrm{g\ cm^{-3}}\) with Si 36.01, O 49.71, Al 5.47, K 2.40, Na 2.42, Mg 1.84, Fe 1.17, and Other 0.98; \(\langle Z\rangle = 11.13\), \(\langle A\rangle = 22.41\), and \(\langle Z^2/A\rangle = 5.53\). The paper explicitly links these quantities to ionization, spallation, pair production, and bremsstrahlung rates [2510.06150].

Rather than sampling muons from a single mean underground energy, the simulations used the full angle-dependent \(\Phi_\mu(E_\mu,\theta,\phi)\). For SPH, \(\langle E_\mu\rangle \approx 98\ \mathrm{GeV}\), but replacing the angle-dependent distribution by a single mean energy elevated secondary fluxes by about 7.6% for neutrons. The authors therefore treat angle-dependent sampling as necessary for accurate background prediction in complex topography [2510.06150].

| Research space | Muon-induced neutron flux at rock–cavern boundary | Muon-induced \(\gamma\)-ray flux at rock–cavern boundary |
|---|---:|---:|
| SPH | \((3.78 \pm 0.61_{\mathrm{sys}})\times10^{-3}\ \mathrm{m^{-2}\ s^{-1}}\) | \((5.54 \pm 0.91_{\mathrm{sys}})\times10^{-1}\ \mathrm{m^{-2}\ s^{-1}}\) |
| CLI | \((3.97 \pm 0.65_{\mathrm{sys}})\times10^{-3}\ \mathrm{m^{-2}\ s^{-1}}\) | \((6.51 \pm 1.06_{\mathrm{sys}})\times10^{-1}\ \mathrm{m^{-2}\ s^{-1}}\) |

For context, the primary muon fluxes at the boundary were \((2.39 \pm 0.25_{\mathrm{sys}})\times10^{-1}\ \mathrm{m^{-2}\ s^{-1}}\) for SPH and \((2.59 \pm 0.26_{\mathrm{sys}})\times10^{-1}\ \mathrm{m^{-2}\ s^{-1}}\) for CLI. Neutron thresholds were also reported. At SPH, the neutron fluxes above 1 MeV, 10 MeV, 100 MeV, and 1000 MeV were \((1.62 \pm 0.26_{\mathrm{sys}})\times10^{-3}\), \((8.32 \pm 1.36_{\mathrm{sys}})\times10^{-4}\), \((3.44 \pm 0.56_{\mathrm{sys}})\times10^{-4}\), and \((1.49 \pm 0.24_{\mathrm{sys}})\times10^{-5}\ \mathrm{m^{-2}\ s^{-1}}\), respectively; at CLI the corresponding values were \((1.79 \pm 0.29_{\mathrm{sys}})\times10^{-3}\), \((9.39 \pm 1.54_{\mathrm{sys}})\times10^{-4}\), \((3.91 \pm 0.64_{\mathrm{sys}})\times10^{-4}\), and \((1.67 \pm 0.27_{\mathrm{sys}})\times10^{-5}\ \mathrm{m^{-2}\ s^{-1}}\) [2510.06150].

The predicted neutron spectrum has four components: a thermal peak around \(2.5\times10^{-8}\ \mathrm{MeV}\), an epithermal/fast region, evaporation neutrons roughly \(0.1\)–\(10\ \mathrm{MeV}\), and a spallation component peaking near \(100\ \mathrm{MeV}\) with a high-energy tail to multiple GeV. An absorption-like feature just above \(10^{-3}\ \mathrm{MeV}\) is described as plausibly linked to resonances in \({}^{23}\mathrm{Na}\) and \({}^{39}\mathrm{K}\) in the local rock. The study also reports that the electromagnetic component dominates the total muon-induced flux at CURIE depth, while at approximately GeV energies the proton flux approaches and then exceeds the neutron flux [2510.06150].

Uncertainties are explicitly systematic-dominated: \(\pm 5\%\) from correlated rock density and muon-energy effects, \(\pm 10\%\) from the Sternheimer parameter approximation, \(\pm 10\%\) from live-time normalization, \(\pm 6.5\%\) from laboratory wall-area uncertainty, and \(\pm 0.5\%\) from Geant4 physics-list choice. Muon bundles are neglected on the basis of shallow-depth measurements indicating single-muon fractions \(\gtrsim 97\%\) in the 100–300 m.w.e. range [2510.06150].

## 5. Experimental scope and infrastructure requirements

The published CURIE studies consistently frame the facility as suitable for detector R\&D, low-threshold detector development, calibration, and environmental studies on superconducting circuits where moderate muon rates are acceptable. The literature also states that experiments using passive shielding and/or active muon veto systems can mitigate residual muon-induced backgrounds such as spallation neutrons. Horizontal access, stable rock, and existing infrastructure are presented as expanding the range of deployable instruments relative to limited U.S. shallow facilities [2411.01626; 2510.06150].

The Snowmass 2021 topical report on underground facilities for the Cosmic Frontier does not explicitly mention CURIE, Henderson Mine, Edgar Mine, or any Colorado-based underground facility proposals or sites. It does, however, articulate facility requirements for next-generation underground dark-matter programs that are directly applicable to the concept of a Colorado-based laboratory [2212.00868].

For noble-liquid experiments, the report identifies next-generation detector scales of 50–300 tons and states that the detector must be in an experiment hall that is large both in floor space and ceiling height, with room for external active veto and shielding. It further specifies significant continuous cooling power, underground liquid-nitrogen capacity extrapolating from XENON1T’s \(10\ \mathrm{m^3}\) tank to a next-generation scale of order \(100\ \mathrm{m^3}\), and on-site cryogenic distillation columns. The cited examples include a 5.5 m purification column releasing some 10s of kW of heat into the cavern and a 3.8 m radon distillation column requiring a few kW of cooling power. The report also calls for large underground low-radon clean rooms with staging areas, cleanliness better than ISO-6, radon concentration better than \(100\ \mathrm{mBq/m^3}\), and substantial vertical space and cranes for heavy lifts [2212.00868].

For cryogenic bolometers, the report states that a vertical space of at least 4 m is highly beneficial for \({}^3\mathrm{He}/{}^4\mathrm{He}\) dilution refrigeration, that environmental vibrations should be below \(10^{-7}\ g\sqrt{\mathrm{Hz}}\) at all frequencies, and that superconducting electronics such as SQUIDs require quiet electromagnetic conditions and minimized DC magnetic fields. It also emphasizes the value of underground R\&D space for developmental and R\&D-scale efforts because long-duration events and pile-up complicate above-ground characterization [2212.00868].

For other technologies, including skipper CCDs, point-contact Ge detectors, low-pressure gas directional detectors such as CYGNUS, scintillator experiments, and superheated water detectors, the report states that present facility needs are relatively compact in footprint and would generally be met if the environmental controls required by noble-liquid and cryogenic-bolometer experiments are available. This suggests that a Colorado laboratory with the environmental conditions measured and modeled at CURIE could support compact early-stage directional or low-threshold programs while remaining especially relevant to R\&D-scale deployments [2212.00868].

## 6. Role, limitations, and broader context

CURIE occupies an intermediate position between surface or minimally shielded installations and deep underground laboratories. The measured reduction in muon flux, approximately \(700\times\), exceeds the values cited for PNNL at approximately 30 m.w.e. (\(\approx 6\times\)) and Fermilab shallow sites at approximately 225–300 m.w.e. (\(\approx 200\times\)–\(400\times\)). At the same time, CURIE remains a shallow-underground facility with residual muon intensities of order \(10^{-1}\ \mathrm{m^{-2}\ s^{-1}}\), not a deep laboratory in the SNOLAB or SURF sense [2411.01626].

Several limitations follow directly from this depth regime. Directional muon-flux variations due to mountain topography require orientation-aware shielding and veto strategies. Muon-induced neutrons and spallation products in rock and shielding require continued characterization and mitigation. The 2025 simulation further shows that local geology and overburden geometry alter yields and energy spectra, so simple depth-based extrapolation is insufficient for high-fidelity background budgeting [2411.01626; 2510.06150].

A common misconception would be to treat all \(0.4\ \mathrm{km.w.e.}\)-class sites as interchangeable once total muon flux is matched. The CURIE studies argue against that simplification: they report measurable azimuthal anisotropy, site-specific rock-composition effects, and a 5–8% bias in secondary fluxes when angle-dependent muon energies are replaced by a single \(\langle E_\mu\rangle\). The Snowmass report also cautions, in a different context, that depth requirements for future low-threshold cryogenic-bolometer experiments remain difficult to assess until the relevant R\&D matures [2411.01626; 2510.06150; 2212.00868].

Within the U.S. underground ecosystem, the Snowmass report concludes that current host laboratories are mostly full and that there is a clear need for expanded underground space tailored to needs beyond the late 2020s. A plausible implication is that a facility with CURIE’s measured accessibility, infrastructure, and environmental characterization could augment national capacity for shallow- to intermediate-depth low-background work, especially for technology development, environmental studies, and experiments that can combine passive shielding with active vetoing. In that sense, CURIE is best understood not as a substitute for deep laboratories, but as a quantified and increasingly well-characterized platform for the segment of low-background research that benefits from substantial muon attenuation without the logistical constraints of deep-lab access [2212.00868; 2411.01626].

Source: https://www.emergentmind.com/topics/colorado-underground-research-institute-curie