Gamma-, neutron-, and muon-induced environmental background simulations for 100Mo-based bolometric double-beta decay experiment at Jinping Underground Laboratory

NUCLEAR PHYSICS AND INTERDISCIPLINARY RESEARCH

Gamma-, neutron-, and muon-induced environmental background simulations for ¹⁰⁰Mo-based bolometric double-beta decay experiment at Jinping Underground Laboratory

Wei Chen，

Long Ma ，

Jin-Hui Chen ，

Huan-Zhong Huang

Nuclear Science and Techniques

Vol.34, No.9

Article number 135

Published in print Sep 2023

Available online 20 Sep 2023

DOI：10.1007/s41365-023-01299-9

44806

The sensitivity of an experiment to detect the Majorana neutrino mass via neutrinoless double beta decay (0νββ) strongly depends on the rate of background events that can mimic this decay. One major source of this background is the radioactive emissions from the laboratory environment. In our study, we focused on assessing the background contributions from environmental gamma rays, neutrons, and underground muons to the Jinping bolometric demonstration experiment. This experiment uses an array of lithium molybdate crystal bolometers to probe the potential 0νββ decay of the ¹⁰⁰Mo isotope at the China Jinping Underground Laboratory. We also evaluated the shielding effectiveness of the experimental setup through an attenuation study. Our simulations indicate that the combined background from environmental gamma rays, neutrons, and muons in the relevant ¹⁰⁰Mo 0νββ Q-value region can be reduced to approximately 0.003 cts/kg/keV/yr.

Neutrinoless double beta decayGeant4CUPID

Introduction

Unveiling the properties of neutrinos is fundamental, and searching for Neutrinoless double beta (0νββ) decay stands as a pivotal approach to this goal. Observing the 0νββ decay experimentally can directly substantiate the Majorana nature of neutrinos, implying that neutrinos are their own antiparticles. This would open new frontiers in physics, extending beyond the standard framework of particle physics [1].

Numerous experimental techniques have been employed to detect the rare 0νββ events. While direct evidence of 0νββ remains elusive, strict limits on the half-life of various double beta decay nuclides have been established [2-5]. Among the tools used, crystal bolometers, which operate at ultra-low temperatures, stand out. They are preferred for 0νββ searches due to their high energy resolution, efficient detection, and stable operation [6]. These cryogenic bolometers employ crystals embedded with double-beta decay isotopes to detect the 0νββ decay. The forthcoming CUPID experiment (CUORE with Particle Identification), integrating a high-performance heat-light dual readout to discern alpha particle backgrounds, is projected to enhance the sensitivity of 0νββ detection in the future [7-9].

A critical challenge in 0νββ detection is the potential interference from external radioactivity, as it can significantly hamper the experiment’s potential. Efforts have thus been made to understand and mitigate sources of background interference. Previous research has identified multiple culprits for misleading counts in the 0νββ search energy region of interest (ROI), including cosmogenic activation, natural radioactive contamination from materials, and environmental radioactive emissions in the laboratory [10-14]. Given that cosmogenic activation at sea level notably influences the detectors, experiments are usually conducted in deep mines or under mountains; these locations, shielded by the earth’s crust or mountain rock, offer protection from cosmic radiation [15, 16]. Besides residual cosmogenic background, environmental radioactive interference remains crucial for underground setups. Main sources include gamma rays from rocks, neutrons from surrounding materials, and events induced by underground muons [17, 18]. Passive shields are typically deployed to counter gamma and neutron backgrounds, while muon vetoes are used to effectively tag and exclude underground muon events and their induced activities.

Given that background levels are pivotal for experimental sensitivity, thorough studies on environmental background shielding and subsequent evaluations of any remaining background are crucial during the design phase of the detector system. In this paper, leveraging a Geant4-based Monte Carlo simulation, we delve into the background contributions—due to environmental gamma rays, neutrons, and cosmic muons—to the Jinping CUPID demonstration experiment. This experiment, featuring a 10 kg cryogenic crystal bolometer at the China Jinping Underground Laboratory (CJPL), aims to showcase the advanced technologies intended for the next-generation bolometric experiment targeting the 0νββ decay of ¹⁰⁰Mo. This paper is structured as follows: the simulation setup and methods are discussed in Sect. 2. Findings related to particle attenuation and background are elaborated in Sect. 3. Conclusions and discussions can be found in Sect. 4

Simulation setup

We utilized the Geant4 simulation toolkit (version 4.10.5) to model radiation shielding and assess background contributions [19]. Geant4 is renowned for offering an extensive array of particle-material interaction cross sections, allowing precise characterization of particle production and transport properties [20-22]. For our study, we employed the "Shielding" physics list, which is tailored for simulations associated with deep underground rare event experiments. This physics list incorporates a series of interaction models, each detailing specific physical processes. For neutrons with energies up to 20 MeV, we leveraged high-precision models and cross-sections to depict processes such as elastic and inelastic scattering, capture, and fission. To study muon-induced backgrounds, we also applied the "QGSP_BERT" physics list as a supplementary check, given its suitability for cosmic ray-related processes. Event generation for gamma rays and neutrons was based on the in-situ spectral measurements taken in the CJPL experimental hall, as depicted in Fig. 1. We assumed isotropic radiation for both gamma and neutron sources. The initiation points for these environmental background sources were set in line with the hall’s dimensions (14 m in height and 12 m in width). Moreover, the muon simulation factored in experimentally measured zenith and azimuthal angle distributions.

Fig. 1

(Color online) Energy spectra of CJPL environmental gamma (left) and neutron (right). The sampled Monte Carlo spectra (solid lines) are shown in comparison with the results obtained by the in-situ measurements (dashed lines) [23, 28].

Figure 2 shows a schematic view of the experimental setup, including the detector crystal array, background shields surrounding the detector, and cooling system used to maintain low temperatures. The configurations of the detector and shielding systems are implemented in the simulation code. The CUPID-CJPL demonstrator comprised an array of thirty-six lithium molybdate (Li₂MoO₄) or LMO cryogenic bolometers, having a cumulative crystal mass of approximately 10 kg. The LMO crystals not only serve as the source for the 0νββ decay but also function as heat absorbers. This array is organized into nine layers, with each layer containing four bolometers. Each of these bolometers is cooled to an ultra-low temperature of about 10 mK using a dilution refrigerator. The individual bolometer module integrates a 4.5 cm×4.5 cm×4.5 cm LMO cubic crystal and is equipped with a heat-light dual-phase readout. For these modules, 280 g LMO crystals were specifically chosen. The cooling apparatus, or the cryostat, consists of three concentric vessels dedicated to thermal shielding. To complement these thermal barriers, radiation shields were strategically positioned both inside and surrounding the cryostat to counteract radiation emanating from the cryogenic system. Within this arrangement, a substantial copper shield (having a thickness of 12 cm and a diameter of 50 cm) was situated above the detector tower to mitigate external radioactive interference. To further fortify the system, passive shielding materials were employed. This includes a 12 cm thick layer of copper, designed to diminish environmental gamma rays, and a 20 cm thick layer of 5%-borated polyethylene (PE) tailored to intercept environmental neutrons. The salient features of this experimental setup are detailed in Table 1.

Key elements of the experimental setup.

Element	Material	Mass (kg)	Surface (cm²)
Crystal array	LMO	10.0	4374.0
Thermal shield	Cu-OFE	61.7	23850.9
Top Cu shield	Cu-OFE	209.7	4736.9
Lateral Cu shield	Cu-OF	2990.4	60880.0
Lateral PE shield	PE	1343.9	114300.0

Fig. 2

(Color online) Schematic view of the cryogenic system, shielding system (a), and crystal detector array (b) implemented in the Monte Carlo simulation.

Results and discussion

3.1

Attenuation study

Numerous experimental studies have established that materials with high atomic mass numbers, such as lead and copper, excel at attenuating gamma rays. Notably, copper typically possesses lower inherent radioactivity than lead. As a result, high-purity oxygen-free copper has emerged as a preferred material for crafting shielding systems in sensitive underground experiments related to rare events.

For the Jinping bolometric demonstration, the design of the background shielding system was informed by a Monte Carlo simulation. Utilizing the Geant4 toolkit, we simulated particle interactions with materials to quantitatively assess the attenuation capacities of the shielding against gamma rays and neutrons.

Figure 3 shows the gamma-ray absorption fraction as a function of the gamma energy for different copper shielding thicknesses. The incident gamma energy was sampled based on the experimentally measured CJPL in situ gamma spectra (Fig 1). Our results indicate that the gamma absorption at different energies improves significantly by increasing the thickness of the material. In addition to high-energy environmental gamma rays, fast neutrons with energies above 1 MeV are serious background sources for rare underground event experiments. Borated polyethylene (BPE) is an effective material for fast neutron shielding [25, 24]. The results of the neutron-shielding performance are shown in Fig. 4 compared with the gamma results. The black and blue curves represent the gamma and neutron absorptions of copper and 5% BPE, respectively. It was observed that the neutron shielding was improved by increasing the thickness of BPE. A reduction in gamma and neutron fluxes by a factor of ~10 can be realized by deploying a 12-cm thick copper shield and an 8-cm 5% BPE shield, respectively. Thus, for the CJPL bolometric experiment, with 12-cm thick copper and 20-cm thick 5%-borated polyethylene external shields, we expect a reduction exceeding 10-fold in the environmental gamma and neutron fluxes within the detector system.

Fig. 3

(Color online) Simulated results of the gamma absorption for varying thicknesses of the copper brick.

Fig. 4

(Color online) Simulated results of the attenuation coefficient as a function of copper shielding thickness for gamma and borated polyethylene thickness for neutron. The input energy spectra for gamma and neutron are shown in Fig. 1.

3.2

Environmental gamma induced background

Different background sources have the potential to simulate a 0νββ signal within the decay-energy ROI of 0νββ. A predominant source of such radiation is environmental gamma rays. For underground experiments, these gamma rays largely stem from the decay of radioactive nuclei and processes induced by high-energy neutrons. Moreover, the lab’s infrastructure and its experimental apparatus can introduce additional gamma background. In our research, we delved into evaluating the contribution of environmental gamma background to a ¹⁰⁰Mo 0νββ detection using a bolometric method.

The gamma radioactivity of the CJPL environment has been extensively profiled in prior studies through in-situ gamma-ray measurements employing a high-purity germanium spectrometer [26-30]. The naturally occurring gamma spectrum spans up to 3 MeV; events exceeding this energy are typically the result of interactions between neutrons, muons, and rocks. Given that the isotope ¹⁰⁰Mo boasts a relatively high Qββ value (approximately 3034.40(17)keV), we designated the ROI to range from 3000–3060 keV. This decision factored in the energy resolution of the crystal bolometer when assessing the background. The approach to parameterizing the energy resolution is as follows: $σ (E) = \sqrt{0.49 + {(0.058 \times \sqrt{E})}^{2}},$ (1) where E denotes the measured energy in keV based on the CUPID-Mo study [31].

The gamma spectrum measured using a high-purity germanium (HPGe) detector in CJPL Phase II was used as the input for our Geant4 simulation [26]. In this study, given that we focused on the ¹⁰⁰Mo 0νββ ROI, only high-energy gamma-rays (Eγ>3000 keV) were considered. The integral flux of the high-energy gamma rays was calculated as 2.51 ± 0.07 × 10^-6 cm^-2s^-1 based on the experimentally measured results.

In bolometric experiments, the array structure of the detector system allows us to apply the anti-coincidence cut criterion. This helps in rejecting high-energy external background events that pass through several crystal modules. Figure 5 illustrates the concept of “multiplicity”. An event that deposits energy in a single crystal module is defined as having “multiplicity=1”. Conversely, if the energy is deposited across multiple crystal modules, it’s classified as a “multiplicity>1” event. By selecting the “multiplicity=1” events, we can significantly filter out high-energy external gamma rays, as well as neutrons and muons, that traverse multiple crystals.

Fig. 5

(Color online) Definition of the multiplicity.

Figure 6 presents the simulated results for the gamma-induced background. The results for different cutting criteria are shown for comparison. It is observed that by requiring “multiplicity=1”, the environmental gammas contribute approximately 1.3 × 10^-3 cts/keV/kg/yr background in the ¹⁰⁰Mo 0νββ ROI (3000-3060 keV). The uncertainty of the result is estimated by considering the statistical uncertainties of the flux measurements and selection of PhysicsList (“Shielding” vs “QGSP_BERT_HP”). The relative difference was found to be within 10%.

Fig. 6

(Color online)Simulated results of the gamma-induced background.

3.3

Environmental neutron induced background

High-energy neutrons present another significant background for 0νββ experiments. In deep underground laboratories, neutron fluxes below 10 MeV primarily arise from spontaneous fission (such as from U-238) and (α,n) reactions, which involve interactions between alphas from natural radioactive sources and light target nuclei within the rock. Beyond production by the inherent radioactivity of materials, neutrons can also be produced by cosmic-ray muons and their secondary particles through nuclear reactions [32, 33]. The energy of these muon-induced neutrons can reach up to several gigaelectronvolts.

Systematic measurements of the fluxes of low-energy thermal neutrons and high-energy fast neutrons have been performed in deep underground laboratories [23, 34-36]. In our simulation study, we evaluated the neutron-induced background based on the in situ fast neutron spectrum measured at the CJPL Hall-A [37]. Monte Carlo neutrons with kinetic energies exceeding 1 MeV were sampled based on the experimental flux spectrum results. The total flux was approximately (1.51 ± 0.03) × 10^-7 cm^-2s^-1 for E_n>1 MeV without the polyethylene shielding. The results for the residual neutron-induced background with the copper and PE shielding setups are shown in Fig. 7. We observed that the environmental neutrons contribute approximately 1.6 × 10^-3 cts/keV/kg/yr in the ¹⁰⁰Mo 0νββ ROI with "multiplicity=1" cut applied. It should be noted that there are currently no published experimental data on CJPL environmental neutron flux exceeding 20 MeV. Thus, high-energy neutrons (E_n>20 MeV) were not considered in our simulations. Nevertheless, according to existing simulation studies performed for LNGS, the integral flux exceeding 10 MeV (mainly from muon-induced processes) is three orders of magnitude lower than that below 10 MeV (mainly from radioactivity-induced processes) [38, 41]. For the CJPL, owing to the larger depth, the flux of cosmic-muon-induced neutrons was even lower. The relative background contribution from E_n>20 MeV high energy neutron was expected to be low.

Fig. 7

(Color online) Simulated results of the neutron-induced background.

3.4

Underground cosmic muon induced background

Cosmogenic backgrounds play a pivotal role in rare-event search experiments [39]. At Earth’s surface, cosmic rays predominantly consist of high-energy muons [40]. The significant rock overburden of CJPL (equivalent to 6700 m of water) dramatically reduces the cosmic muon flux underground—by about eight orders of magnitude compared to the flux at sea level. As a result, the cosmogenic background doesn’t pose a significant threat to the Jinping bolometric experiments. Nevertheless, for the heightened sensitivity objectives of upcoming 0νββ experiments, a thorough evaluation of the residual muon-induced background is imperative.

To estimate the muon-induced background level, it’s essential to understand both the energy and angular distributions of underground muons in a given experimental cavern. The muon energy was derived using a parametric function [41]: $\frac{d N}{d E_{μ}} = C e^{- b h (γ_{μ} - 1)} {(E_{μ} + ε_{μ} (1 - e^{- b h}))}^{- γ_{μ}},$ (2) where b=0.4/km⋅w⋅e, γμ=3.77, and Eμ=693 GeV were used, as suggested.

Figure 8 shows the underground muon spectrum from a Geant4-based simulation study and compares it with the Mei-Hime parameterization results [42]. The parametric function accurately described the Geant4 simulation results. For the angular distribution, additional information regarding rock overburden is required to determine the flux in a specific direction. In this study, the zenith and azimuthal angle distributions of the underground muons were sampled according to JNE measurements [43]. The experimental distributions can be described using the parametric function discussed in Ref. [44]. The total underground muon flux used as the input for the simulation was 3.53 × 10^-10 cm^-2s^-1 based on the measurement at CJPL Phase II.

Fig. 8

(Color online) Energy spectrum of the CJPL muons. Results from the parametric function are shown in comparison with that from the Monte Carlo simulation.

The simulated results of the muon-induced background contributions are shown in Fig. 9. We observed that the background contribution of muon events, including muon-induced secondaries, corresponds to a low level of 2.2 × 10^-5 cts/keV/kg/yr in the ¹⁰⁰Mo 0νββ ROI. It was also observed that the multiplicity anti-coincidence cut significantly suppressed the muon-induced background by an order of magnitude.

Fig. 9

(Color online)Simulated results of the muon-induced background contribution.

Conclusion

We conducted a simulation study evaluating the background shielding for the Jinping ¹⁰⁰Mo-based bolometric demonstration experiment aimed at searching for neutrinoless double-beta decay. Having established the shielding arrangement through an attenuation study, we assessed the environmental gamma, neutron, and muon-induced backgrounds. Our findings indicate that within the ¹⁰⁰Mo 0νββ ROI, gamma rays and neutrons contribute to background levels of approximately 1.3 × 10^-3 cts/keV/kg/yr and 1.6 × 10^-3 cts/keV/kg/yr, respectively. The underground cosmic ray muon’s contribution was minor, only reaching 2.2 × 10^-5 cts/keV/kg/yr in the ROI when the anti-coincidence cut was applied. Additionally, given the bolometer’s slow data acquisition rate, potential pile-up events, notably coincident low-energy gamma rays, could also add to the discernible experimental background. It’s worth noting that our simulations were done for a demonstration experiment without supplementary passive background shielding. In CJPL Phase II, a meticulously designed combined shielding apparatus (comprising 20 cm of lead and 100 cm of PE) was implemented for radioactive background protection. We anticipate a significant reduction in the combined contributions from environmental gamma and neutron backgrounds to levels below 10^-3 cts/keV/kg/y, aligning with our objectives for this demonstration experiment. For the 0νββ experiment, background plays a pivotal role, directly influencing experimental sensitivity. Therefore, curbing background contribution emerges as a critical facet of bolometric experimental studies and advancements. Regarding the environmental background, precise in-situ measurements of both neutron and gamma spectra, extending into high-energy regions, are paramount for designing shielding systems tailored to future high-sensitivity 0νββ research endeavors at CJPL.

References

S.M. Bilenky, C. Giunti,

Neutrinoless double-beta decay: A probe of physics beyond the Standard Model

. Int. J. Mod. Phys. A 30, 1530001 (2015). doi: 10.1142/S0217751X1530001X