Simulation-based design study for the passive shielding of the COSINUS dark matter experiment.

The COSINUS (Cryogenic Observatory for SIgnatures seen in Next-generation Underground Searches) experiment aims at the detection of dark matter-induced recoils in sodium iodide (NaI) crystals operated as scintillating cryogenic calorimeters. The detection of both scintillation light and phonons allows performing an event-by-event signal to background discrimination, thus enhancing the sensitivity of the experiment. The choice of using NaI crystals is motivated by the goal of probing the long-standing DAMA/LIBRA results using the same target material. The construction of the experimental facility is foreseen to start by 2021 at the INFN Gran Sasso National Laboratory (LNGS) in Italy. It consists of a cryostat housing the target crystals shielded from the external radioactivity by a water tank acting, at the same time, as an active veto against cosmic ray-induced events. Taking into account both environmental radioactivity and intrinsic contamination of materials used for cryostat, shielding and infrastructure, we performed a careful background budget estimation. The goal is to evaluate the number of events that could mimic or interfere with signal detection while optimising the geometry of the experimental setup. In this paper we present the results of the detailed Monte Carlo simulations we performed, together with the final design of the setup that minimises the residual amount of background particles reaching the detector volume.

Introduction

According to the Standard Cosmological Model (CDM), Dark Matter (DM) accounts for about % of the total matter density of the Universe. The presence of a dark and non-baryonic matter species is confirmed by a variety of cosmological and astronomical observations. However, a definitive proof of the existence of DM particles either in indirect, direct or accelerator searches, is still missing [1].

Direct searches aim for the detection of DM-induced scattering events in a given target. In this field the DAMA/LIBRA experiment [2], operating 250 kg of high purity tallium-doped sodium iodide (NaI(Tl)) crystals at room temperature, reports a statistically robust (12.9 [3]) annual modulation of the event rate, with period and phase (1 year and 152.5 day) compatible with a halo of DM particles in the milky way [4]. Nevertheless, this result is not confirmed by any other DM experiment. To cross check the DAMA/LIBRA claim in a model-independent way, several experiments using the same target material as DAMA/LIBRA (i.e. NaI crystals) are planned (SABRE [5], PICO-LON [6], COSINUS [7]) or already taking data (COSINE [8], ANAIS [9]). Among them, COSINUS is the only one using a unique detector technology to operate the NaI crystals as scintillating calorimeters (see scheme of detector module in Fig. 1): a double cryogenic read-out channel approach provides the possibility to simultaneously detect the scintillation light and the phonon signal, thus allowing to disentangle e/ events from nuclear recoils on an event-by-event basis. Only below a few keV, the discrimination power decreases, while at higher energies no e/ event will leak into the acceptance region for nuclear recoils. This particle identification technique is a crucial benefit, as most DM models predict DM-nucleus scattering, while e/ events are the dominant background. Comparing to the energy scale probed by scintillation-only experiments, this gives the unique opportunity of testing their results in a largely background-suppressed environment [7]. Moreover, the approach adopted by COSINUS allows to directly set the energy scale via the phonon channel, independent of the light channel and of quenching factor values (see next section). This unique feature is of crucial importance given the recent diverging quenching factor measurements used to interpret scintillation-only experiments (see Ref. [10] and references therein for an introduction to the topic).

Fig. 1

Scheme of a COSINUS detector module consisting of a NaI target crystal coupled to a carrier crystal (phonon detector) and a beaker-shaped light detector. Both channels are read out via a Transition Edge Sensor (TES)

The COSINUS experiment

The largest part of the energy deposited by a particle interacting in a scintillating crystal is converted into lattice vibrations (phonons). These are measurable as heat when operating the detector close to absolute zero (T (mK)). At the same time, a smaller energy fraction of (1–10 %) is emitted in the form of scintillation light. While the amount of scintillation light produced by e/ events strongly differs from that produced by nuclear recoils, the phonon (heat) signal is almost independent of the interacting particle type and can be used to precisely measure the deposited energy. We define the Light Yield (LY) as the ratio between light and phonon signal to quantify the separation between e and nuclear recoil events on an event-by-event basis. By definition, we set the LY of e/ particles to 1 at the energy of the calibration source (typically 122 keV gammas from Co-57). Nuclear recoils and alpha particles feature a lower light yield quantified by particle- and energy-dependent Quenching Factor (QF) values. The QF is defined as the ratio between the scintillation light produced by the energy deposition of a given particle and the amount of light produced by the interaction of e/ particles of the very same energy. Figure 2 shows the simulated event distribution in the LY versus energy plane, assuming a flat electromagnetic background of 1 and an internal K contamination of 600 Bq/kg. For further details the reader is referred to Ref. [7].

Fig. 2

Simulated COSINUS data in the light yield versus energy plane for a detector reaching the performance goals presented in Ref. [7] and assuming a background level compatible with DAMA/LIBRA: a flat electromagnetic background of 1 and an internal K contamination of 600 Bq/kg responsible for the line at 3 keV. The solid lines mark the upper and lower 90 % boundaries of the e/-band (black) and the nuclear recoil bands for recoils off sodium (blue) and iodine (green), respectively. This plot is based on the energy-dependent quenching factors reported in Ref. [11]. Plot adopted from Ref. [7]

In the baseline design of the COSINUS detector modules (Fig. 1), the NaI target is coupled to a highly sensitive temperature sensor via a thin (1 mm) carrier crystal made of a different material (e.g. ). The interface between carrier and target crystal is made of epoxy resin or silicone oil. The thermometer is a Transition Edge Sensor (TES) consisting of a superconducting tungsten thin film (W-TES, 200 nm) produced via electron-beam evaporation. Over the past 20 years, the CRESST dark matter search [12] has developed and optimised this W-TES technology. The TES is stabilised in its transition between the normal and the superconducting phase. In this configuration the small temperature rise ((K)), due to the energy deposited by a particle interaction in the crystal, is followed by a steep increase in the TES resistance of (m), which is in turn registered by Superconducting Quantum Interference Device (SQUID) amplifiers.

A beaker-shaped silicon light absorber of 40 mm in height and diameter and a wall-thickness of 500 m, encapsulates the NaI crystal. This configuration enhances the light collection efficiency while acting, at the same time, as an active veto against surface -events. The read-out is performed by a TES optimised for the light detector which is directly evaporated onto the silicon surface. Both light and phonon channel TESs are thermally linked to the heat bath. For a summary of the performance obtained during the detector R&D phase, the reader is referred to Refs. [13-17].

The modules described above, composed of both phonon and light detectors, will be arranged in an air-tight copper container to protect NaI crystals, which are highly hygroscopic, from ambient, humid air. The target volume can host up to six layers of containers, each one housing up to twelve modules. Each detector module can be equipped with NaI crystals with a maximum weight of 150 g. A simplified scheme of the overall experimental apparatus is shown in Fig. 3. The space for the target is allocated at the bottom of a cylindrical stainless steel tube, the “dry-well”. This tube is inserted in a tank filled with ultra-pure water and hosts the cryostat. The geometry is conceived to place the detector volume at the center of the tank. For this reason, the dry dilution refrigerator is equipped with a custom-made extension. The water tank will be equipped with photomultiplier tubes (PMTs), thus acting at the same time as passive shielding against ambient radiation and active Cherenkov veto against cosmogenic background. The experimental apparatus will be installed in Hall B of the LNGS underground laboratory.

Fig. 3

Simplified scheme of the cryostat in the dry-well of the water tank. Light blue: stainless steel walls of water tank and dry-well. Orange: copper parts (shielding and cryostat). Solid magenta, dashed black: experimental volume needed for COSINUS-1  (corresponding to a maximum total mass of 3.6 kg of NaI) and -2 (up to a total mass of 10.8 kg of NaI), respectively. Blue: temperature stages of the cryostat. Green: light-tight curtain made of e.g. Tyvek to create a passive layer at the water tank walls. This scheme is based on the favoured solution (Option 4) for the shielding concept as found by MC simulations (see Sect. 5)

The COSINUS experiment is planning a phased approach. In the first phase, called COSINUS-1, up to two layers of detectors (i.e. 24 modules) will be installed. The goal is to provide, exploiting the particle identification capability of our detectors, a cross-check of the DAMA/LIBRA result independent of the dark matter halo (astrophysics) and of the interaction mechanism (particle physics) apart from demanding that dark matter interacts with nuclei and not with electrons in the target. As discussed in detail in Ref. [18], assuming a nuclear recoil threshold on the phonon signal of about 1 keV and the radiopurity level and performance discussed in Ref. [7], the exposure needed to fulfil the COSINUS-1 objective is of the order of 1000 . Depending on the final mass of the single NaI crystal, that will depend on the final phonon sensor configuration (currently under optimisation), the goal can be achieved running about 20 crystals, 50(100) g each, for 3(1.5) years.

In the second phase, called COSINUS-2, the search for the annual modulation signal is planned by increasing the target mass up to the maximum capacity (72 modules in total).

The reach of the physics goal of the project imposes the need of an accurate background budget estimation taking into account both environmental and intrinsic radiation. The objective of the setup is to efficiently shield the former with appropriate materials, while reducing the amount of the latter originating from the unavoidable radioactive contamination contained in both shielding and infrastructure.

Background sources and shielding concept

COSINUS is the unique NaI experiment featuring a discrimination capability between e-particles and nuclear recoil events on an event-by-event basis. In the search for DM particles scattering off atomic nuclei, neutrons interacting in the target crystals can mimic a DM signal and thus have to be considered as the most dangerous type of background. Nevertheless, gamma radioactivity has to be properly taken into account especially in the very low energy region where the discrimination power is weaker (see Fig. 2).

Neutrons, in underground sites, can be produced in reactions initiated both by natural radioactivity (radiogenic neutrons) and cosmic rays (cosmogenic neutrons). Radiogenic neutrons are produced via spontaneous fission and interactions originating from natural radioactive contaminants (i.e. U and Th) in rock, water and concrete surrounding the apparatus as well as in the materials used in the experimental setup. Cosmogenic neutrons are instead produced either in muon-induced spallation processes or by secondary particles generated in muon-induced cascades, i.e. through photon- or hadron-induced spallation or disintegration reactions [19]. While the radiogenic neutron energy is in the fast range (up to 10 MeV), cosmogenic neutrons have energies extending up to few GeV.

Ambient gamma rays originate from long-lived natural radioisotopes K, U and Th present in the rock. Moreover, cosmic muons can induce gamma radiation via bremsstrahlung both in the rock and in high-Z materials possibly used as passive shielding.

Table 1

The five different options for our shielding configuration, featuring different thicknesses of water, Pb, Cu and PE. See Fig. 4 for a schematic view of the examined configurations

Option	Tank	Water	Dry-well	Inner shielding			Cryostat	Top shielding
	(steel)	radius	(stainless steel)	(Pb)	(Cu)	(PE)	(Cu)	(Pb)	(Cu)	(PE)
	(cm)	(cm)	(cm)	(cm)	(cm)	(cm)	(cm)	(cm)	(cm)	(cm)
1	1.5	150	0.4	10	15	10	0.8	10	15	10
2	1.5	200	0.4	0	15	10	0.8	0	40	10
3	1.5	200	0.4	0	15	0	0.8	0	40	0
4	1.5	300	0.4	0	8	0	0.8	0	30	0
5	1.5	300	0.4	0	0	0	0.8	0	40	0

Fig. 4

The five different experimental setups considered for the MC simulation. The Cu cryostat containing the detector volume (cyan box), which will house the NaI detectors, is inserted in a thin stainless steel structure (dry-well) hosting the shielding layers. The detector volume is positioned at the center of a water-filled stainless steel tank. Details about the thickness of the various layers for the five setup options can be found in Table 1

The general strategy in reducing the flux of ambient and radiogenic neutrons and gamma rays is based on the use of successive layers of passive absorbers. An outermost layer of a low-A material (typically polyethylene (PE) or water) acts as an efficient moderator for ambient neutrons. To reduce the environmental gamma flux, an appropriate thickness of a high-Z material (usually lead) is needed. A layer of copper is then used to mitigate the gamma radiation emitted in the decay of the long-lived Pb isotope contained in lead. Finally, to deal with radiogenic neutrons produced in the shielding materials, a further layer of a low-A material is often used as the innermost shield close to the detector volume. In addition to these layers, the natural shielding provided by the rock overburden above the LNGS underground halls (3600 m.w.e. [20]), where the COSINUS apparatus will be located, assures a reduction of the muon flux by a factor of 10 with respect to the surface [21]. Nonetheless, especially in high-Z materials like Pb used as gamma ray absorbers, muon-induced spallation and subsequent hadronic showers result in a further emission of neutrons. A muon veto to tag such events is therefore crucial. The use of a tank filled with ultra-pure water and instrumented with PMTs satisfies both, the requirement for an external passive shield absorbing ambient radiation and the need for an active veto. Moreover, ultra-pure water has a lower contamination level than PE and a proper tuning of the water thickness allows reducing the amount of high-Z material layers, thus finally minimising the sources of radiogenic as well as those of muon-induced backgrounds.

Simulation

To estimate the overall number of particles reaching the detector volume in our setup, we performed Monte Carlo (MC) simulations based on Geant4 (v10.2.3) [22-24], SOURCES4C [25] with extended alpha energies and higher-A target isotope lists and MUSUN [26].

Starting from the basic configuration sketched in Fig. 3 and detailed in Sect. 2, we simulated five simplified experimental layouts as shown in Fig. 4. The cryostat, represented by a single Cu layer with a thickness of 0.8 cm, containing the detector volume (cyan box), is inserted in a thin (0.4 cm) steel structure, the dry-well, which also hosts the inner passive shielding layers. They may consist of one or two layers of high-Z materials (Pb and Cu, or Cu only) and potentially a further layer of polyethylene (PE). We moreover considered distinct thicknesses of the water layer (between 150 and 300 cm) while, for the stainless steel tank, we assumed an average wall thickness of 1.5 cm for all the options under study. The thicknesses of all shielding layers used in the five configurations depicted in Fig. 4 are listed in detail in Table 1.

For each material and its relative thickness used in the different design options shown in Fig. 4, we compute the radiogenic neutron and gamma fluxes, assuming the intrinsic radioactive contamination values quoted in literature, as reported in Table 2. The intrinsic neutron and gamma yields are used as a benchmark and compared to the attenuation of the environmental fluxes offered by the various shielding configurations to optimise the thicknesses of the different layers. The goal of the optimisation is to reduce the ambient background to a negligible level with respect to the unavoidable intrinsic radiogenic background.

Table 2

Activity of the main contaminants in the materials considered for the COSINUS setup. Secular equilibrium is assumed for the decay chains. The rightmost column states the neutron yield due to and spontaneous fission reactions attained using the SOURCES4C code [25]

Material	^238U	^235U	^232Th	Ref.	^40K	^60Co	^137Cs	Ref.	Neutron yield
	(mBq/kg)	(mBq/kg)	(mBq/kg)		(mBq/kg)	(mBq/kg)	(mBq/kg)		(cm^-3 s^-1)
Stainless steel	<0.2	–^∗	<0.1	[30]	<5.2	1.9	<0.6	[31]	3.0·10^-12
Pb	<0.01	–^∗	<0.07	[30]	–	–	–	–	1.2·10^-13
Cu	<0.065	–^∗	<0.002	[32]	<0.34	0.21	<0.03	[31]	6.6·10^-13
PE	<3.8	<0.37	<0.14	[33]	0.7	<0.1	0.06	[33]	9.4·10^-12

no measured value given, natural abundance assumed

The background rates will nonetheless depend on the final detector layout (e.g. size and number of modules) which is currently under definition. A detailed background model will therefore be the topic of future work. For the preferred shielding option, however, we will briefly discuss the effect of detector size and segmentation at the end of this section.

Besides using passive shields, an active muon tagging system will be crucial to reduce the dominating neutron background contribution presented in Sect. 4.3. For this purpose, we will use the water tank as an active Cherenkov veto. Its tagging efficiency will depend on the size of the tank and on the number, distribution, trigger condition and performance of PMTs used to measure the Cherenkov light.

Ambient neutrons

The energy spectrum of ambient neutrons in the LNGS underground halls, taken from Ref. [27], is shown in Fig. 5a. This spectrum was inferred from a detailed study revisiting the neutron flux at LNGS and is well in agreement with measurements. The integral flux in the range 1−500 keV is , while the integrated flux above 500 keV is .

Fig. 5

The plots show spectra of background components external to our experimental setup and the respective fractions of surviving particles as a function of the thickness of the water layer. The rather low-energy ambient neutrons are much more efficiently shielded than comparably high-energy muon-induced neutrons. The efficiency for shielding ambient gamma rays with water lies between the other two values

We performed a Geant4-based simulation, propagating neutrons, generated according to the above spectrum, through a large water cuboid. We conservatively considered neutrons entering perpendicularly to the water surface. The fraction of surviving ambient neutrons as a function of the thickness of the water layer is reported in Fig. 5d. An attenuation factor of () can be obtained for a thickness of less than 1 m (1.5 m). The flux of ambient neutrons reaching the detector volume, in any of the considered shielding configurations, will thus be lower than . Projecting this number to the surface of the detector region would translate to less than year entering the detector volume. The ambient neutron flux is hence reduced to a negligible level with respect to the potential background sources presented in the subsequent sections. Moreover, surviving neutrons are efficiently thermalised, so that they do not contribute to the background via nuclear recoils.

To cope with particles impinging the target volume perpendicularly from the top, through the neck of the dry-well where no water is present, we foresee to deploy a layer of a low-Z material around the extension of the cryostat (1.5 m) between the first two cooling stages and the dilution unit (see Fig. 3). For details on the cryostat design the reader is referred to Ref. [29]. This layer will either be composed of PE or could consist of thin PE tanks filled with ultra-pure water and shaped to fit the space around the cryostat prolongation. The latter solution would assure a lower intrinsic radioactive contamination level and facilitate the mechanical handling when lifting the cryostat out of the dry-well.

Radiogenic neutrons

For the intrinsic activity of the materials employed in the various shielding layers and in the infrastructure of the COSINUS apparatus, we assumed the upper limits reported by other rare event searches listed in Table 2. To estimate the flux and the energy spectrum of neutrons originating in the various materials via spontaneous fission and (,n) reactions, we performed a simulation using the SOURCES4C code [25] and assuming secular equilibrium for U and Th decay chains. The integrated radiogenic neutron yield is reported, for each material, in the last column of Table 2, while the energy spectra are shown in Fig. 6. Systematic uncertainties on the neutron flux obtained from SOURCES are discussed in Sect. 4.6. The results are used as input for the Geant4 simulation.

Fig. 6

Energy spectrum of radiogenic neutrons obtained using the SOURCES4C code [25] for PE, Pb and Cu shielding layers. The double peak in the energy distribution of PE is due to the presence of C, which has a high cross section for interactions. All spectra are normalized to unit integral

The result of the simulation is reported in Table 3: for each shielding material within a given design option, we evaluate the number of neutrons reaching the target volume. The sum of all contributions is reported in the last column of the table.

Table 3

Radiogenic neutron background contributions originating from distinct materials used in the various design options. Stated values include 1 statistical uncertainties

Design option	Neutron source	Events with neutrons	Total
		entering detector volume
		(year^-1)	(year^-1)
Option 1	PE	8.86±0.01	9.17±0.01
	Cu	(2.96±0.03)·10^-1
	Pb	(7.41±0.46)·10^-3
	Steel	(3.62±0.34)·10^-3
Option 2	PE	8.91±0.01	9.18±0.01
	Cu	(2.67±0.03)·10^-1
	Steel	(2.06±0.24)·10^-3
Option 3	Cu	2.16±0.01	2.17±0.01
	Steel	(1.19±0.06)·10^-2
Option 4	Cu	(8.44±0.05)·10^-1	(9.31±0.07)·10^-1
	Steel	(8.72±0.16)·10^-2
Option 5	Cu (cryostat)	(2.09±0.02)·10^-1	(4.22±0.05)·10^-1
	Steel	(2.13±0.03)·10^-1

As expected, less employed material leads to a lower number of radiogenic neutrons reaching the detector volume. Moreover, PE, despite its small mass fraction, results in the largest contribution to the total radiogenic neutron flux, (10 events/year), due to the higher contamination levels with respect to other shielding materials (see Table 2) and higher (,n) yield. Radiogenic neutrons originating from the intrinsic radioactivity of target crystals do not contribute to the overall background. Neutrons produced in (,n) reactions inside the crystal can indeed be tagged by detecting the preceding particle having a typical energy in the order of a few MeV, thus inducing a saturated signal. A clearly saturated signal is also obtained for intrinsic neutrons derived from spontaneous fission, as they are emitted together with nuclear fragments with energies of the order of 100 MeV.

Cosmogenic neutrons

As mentioned in Sect. 3, muons reaching the underground halls at LNGS are suppressed by a factor of with respect to sea level thanks to the shielding provided by an average vertical rock overburden of 3600 m.w.e. [20, 21]. Nevertheless, surviving muons produce neutrons through spallation processes in the rock or in high-Z materials surrounding the detector. To estimate the cosmogenic contribution to the overall background, we performed a simulation of the cosmic muon flux reaching the underground site by using the MUSUN code [26]. The results obtained are compatible with values found in literature [19]: the mean energy of muons reaching the underground halls is 270 GeV, the average zenith angle is 0.67 rad and both the distribution of the zenith and azimuthal angles are compatible with the profile of the Gran Sasso mountain (see Fig. 7).

Fig. 7

Spectrum of the (a) energy, (b) zenith angle and (c) azimuthal angle of muons reaching the LNGS underground laboratory, obtained using the MUSUN [26] simulation code

As output of the MUSUN simulation we obtained an integral flux of year, through the surface of a cuboid with dimension m: position, direction and energy of each simulated muon are stored and used as input for the Geant4 simulation. Muons are first propagated through a thickness of at least 2.5 m of rock (to properly take into account the interaction with the rock surrounding the apparatus), and then through the five different options of our experimental setup listed in Table 1. We assumed the water tank positioned at a fixed distance of 0.5 m from one of the walls of the underground hall. Figure 5c shows the energy distribution of cosmogenic neutrons produced in rock and shielding. As expected, the spectrum extends to much higher values with respect to the radiogenic one (see Fig. 6 for comparison).

For each design option in Fig. 4, we simulated between 15 and 30 million muons, roughly corresponding to an exposure between 6.5 and 13 years. The number of events in which cosmogenic neutrons, produced either in the rock or in the experimental setup, enter the detector volume is reported in Table 5. The reason for quoting the number of events rather than the number of neutrons is that in a single muon-induced event many neutrons, as part of a hadronic shower, may enter the target volume at the same time. The various design options exhibit comparable performance, since a notable fraction of neutrons is produced in the rock surrounding the setup. Nonetheless, it is worth noting that the use of a Pb layer in the shielding results in a higher amount of muon-induced neutrons, as the neutron yield scales with the atomic mass number as A [19], while a PE layer provides an additional shielding against these neutrons. Compared to the radiogenic background, however, the cosmogenic contribution is one or two orders of magnitude higher for all the experimental options considered. To not limit our sensitivity we thus plan the use of an active Cherenkov muon veto by instrumenting the water tank with PMTs.

Table 5

List of all the background contributions in every shielding option considered in the Geant4 simulations. The numbers are attained using the simplified geometrical setup detailed in Sect. 3 and represent the number of events per year, in which at least one particle of the corresponding background source enters the detector volume. Stated values include 1 statistical uncertainties

We started dedicated simulations to evaluate the necessary number of PMTs and their arrangement in the tank. The final results will be the subject of a future work. Preliminary analyses already show that a veto system equipped with about 30 PMTs reduces the cosmogenic neutron background by roughly two orders of magnitude, thus reaching the level of the estimated radiogenic one.

Radiogenic gammas

The radiogenic gamma flux originates from the natural decay chains as well as from additional contaminants like K, Co or Cs. To evaluate this contribution we performed a dedicated Geant4 simulation using the contamination levels stated in Table 2. The estimated radiogenic gamma rates are reported in Table 4. The numbers show that PE and steel layers lead to high contributions, if no additional shield is placed closer to the cryostat. Different contaminants may play important roles in different materials. While the largest fraction of events which have its source in PE or Cu originate from the U decay chain, the dominating contribution coming from steel originates from its Co contamination.

Table 4

Radiogenic gammas, originating from distinct materials used in the various design options, entering the detector volume. Stated values include 1 statistical uncertainties

Design option	Gamma source	Events with gammas	Total
		entering detector volume
		(year^-1)	(year^-1)
Option 1	PE	(5.41±0.13)·10^6	(5.68±0.14)·10^6
	Cu	(2.74±0.11)·10^5
	Pb	(8.90±0.89)·10^2
	Steel	(2.70±0.69)·10^1
Option 2	PE	(5.41±0.13)·10^6	(5.68±0.14)·10^6
	Cu	(2.74±0.11)·10^5
	Steel	(1.16±0.10)·10^3
Option 3	Cu	(4.06±0.12)·10^5	(4.08±0.13)·10^5
	Steel	(2.36±0.15)·10^3
Option 4	Cu	(3.97±0.11)·10^5	(4.46±0.14)·10^5
	Steel	(4.89±0.25)·10^4
Option 5	Cu (cryostat)	(1.51±0.02)·10^5	(2.09±0.04)·10^6
	Steel	(1.94±0.04)·10^6

Ambient gammas

For the simulation of the ambient gamma radioactivity in the LNGS underground halls, we adopted the spectrum shown in Fig. 5b taken from Ref. [28]: the integrated gamma flux is 0.23 cm s. The gamma spectrum has been used as input for the Geant4 simulation computing the number of gamma rays transmitted through the water tank. The attenuation factor as a function of water thickness is reported in Fig. 5d: an attenuation factor of the order of () is obtained for 200 cm (300 cm) of water thickness. The use of a high-Z material (Pb and/or Cu) can further reduce the residual flux. The resulting particle rates reaching the detector volume are reported in Table 5. In all shielding options, the thicknesses of the various layers enable a reduction of the ambient gamma flux down to, or below the level of the estimated radiogenic gamma background (see Sect. 4.4).

Systematic uncertainties

Besides statistical uncertainties in the Monte Carlo evaluations and the uncertainties related to the assumed literature values of contamination levels, we provide, in the following, an estimation of systematic uncertainties in the physics models of employed simulation codes.

SOURCES4C carries an uncertainty in the amount of radiogenic neutrons produced via (,n) reactions. While spontaneous fission is well-described, (,n) cross-sections are based on tabulated databases, which may comprise differing values. In Ref. [34], various codes for calculating -induced neutron yields are compared and SOURCES4C is found to agree to experimental data within 20 % in most of the cases.

MUSUN was tested against data of various experiments and found to be consistent with their results [26]. Particularly, the code was validated against measured data of the LVD experiment [35], which is situated in hall A of LNGS, thus yielding a negligible systematic uncertainty in the computed amount, energy spectrum and angular distribution of muons reaching the LNGS underground halls.

Uncertainties in Geant4 mainly persist in the muon-induced neutron yield and in data-driven physics models. In Ref. [36], measurement and simulation of the muon-induced neutron production in a deep underground laboratory are checked against each other. The setting is comparable to the LNGS, so that the quoted agreement within 25 % may be assumed as an estimate on the systematic error. For the neutron transport below 20 MeV, we employ the data-driven high precision neutron transport model. A comparison in Ref. [37] between Geant4 and MCNPX [38] showed an agreement better than 20 %, which we can assume to be an upper limit on the systematic error on low-energy neutron propagation.

Discussion of the results

In Table 5, the rate of particles reaching the detector volume in the considered shielding options is listed. Looking at these results and considering the discussions in the previous sections, we can draw the following conclusions.

Neutrons are the most dangerous background particles. They can mimic a dark matter signal by inducing nuclear recoils in the detector. Since the cosmogenic neutron background is the dominant contribution, the use of an active Cherenkov veto is mandatory to reduce the rate of cosmogenic events to or below the level of the radiogenic one. The use of a large instrumented water tank volume can provide such a muon veto and allows, at the same time, to minimise the amount of further inner shielding layers (Pb or Cu), which in turn leads to a reduction of both the radiogenic neutron and gamma fluxes.

Radiogenic background simulations show that the PE layer as innermost shielding is disadvantageous. Although the gamma background can be efficiently discriminated by using the double read-out channel approach, part of its low-energy fraction may potentially leak into the region of interest for dark matter search. In addition, the total rate in the detectors affects the dead-time of the experiment. Thus, also the minimisation of the gamma background is of importance.

From these considerations, we conclude that options number 4 or 5 in Table 1 best fulfill our requirements. Option 4 is the safer and thus preferred configuration, as it definitely leads to a reduced gamma rate and is less sensitive to the intrinsic radioactivity level of the stainless steel of the dry-well, which is yet to be determined.

It is worth mentioning that the amount of radiogenic background was computed assuming the contamination levels found in literature and reported in Table 2. We are currently performing an intensive screening campaign on all the materials that will be employed for the construction of the apparatus taking the values assumed in the simulations as a benchmark for selection. At the time of writing, we got preliminary results about the screening measurements of the different copper samples that will be used for the manufacturing of the cryostat, detector holders and shielding. The bulk contamination levels we found are in good agreement with or even lower than those assumed in Table 2. For the stainless steel of the dry-well we are currently performing screening measurements for material selection.

The actual background rate measured in our detectors will depend on the size, number, position and housing of the modules. By taking the example of shielding option 4, we briefly discuss the effect of using a segmented target made of an array of small crystals compared to the use of a single crystal of the same integrated mass. For this purpose, we ran two identical simulations: one using a single 2.4 kg NaI cube, the other using 24 smaller 100 g cubes, distributed over a larger volume. Looking at both cosmogenic and radiogenic neutrons, the segmented arrangement leads to an enhanced total rate for the same integrated detector mass, because of the larger surface to volume ratio. On the other hand, the segmentation allows for an anti-coincidence cut, removing events in which two or more detector modules trigger simultaneously. This cut is more efficient with respect to cosmogenic than radiogenic background, because in the first case more particles with higher energies enter the detector volume on average. In our simulation, we applied the anti-coincidence cut assuming a detection threshold of 1 keV. In both arrangements the results yield similar nuclear recoil rates due to cosmogenic neutrons within statistical uncertainties. The obtained rates are of the order of cts kgyear. On the other hand, the radiogenic neutron background is systematically slightly higher in the segmented detector setup, with a rate of cts kgyear compared to cts kgyear. A more detailed study of the total estimated background rate will be the topic of a future publication once the final geometry of the detector volume will be assessed. However, the preliminary number of muon-induced neutron events already demonstrates the necessity to use a Cherenkov muon veto for COSINUS.

Summary and conclusion

In this work, we performed a conceptual design study of the shielding for the COSINUS experiment. We took into account the major background sources external to the NaI detector volume as well as a simplified geometry and compared different shielding configurations. The aim was to minimise the amount of ambient and cosmogenic background particles reaching the detector volume, while using the unavoidable intrinsic radiogenic background generated in the employed material as a benchmark.

This study was used as the basis of the design of the COSINUS apparatus whose construction started October 25, 2021.

Following the layout of option 4 in Table 1, the experimental setup will consist of a water tank with diameter and height of 7 m, providing a shielding of at least 3 m of water all around the detector volume where the NaI crystals will be inserted. The diameter of the dry-well has been set in order to host the 8 cm layer of Cu around the cryostat, extending in height to line up with the inner Cu shield inside the cryostat, as shown in Fig. 3. In this option no lead or polyethylene layers are foreseen. To suppress the dominant cosmogenic neutron contribution, the water tank will be instrumented with PMTs to serve as an active Cherenkov veto.

A publication describing the performance and the geometry of the muon veto is currently in preparation.

A dedicated screening campaign of all materials that will be employed in the construction of the setup and shielding, as well as of those used for detector assembly, is ongoing. The screening results will be used for material selection and as an input to a simulation considering a detailed geometrical description of the entire setup and detector modules for the background modeling in the region of interest for the Dark Matter search.