The projected background for the CUORE experiment

TL;DR

This paper presents a detailed background projection for the CUORE experiment, aiming to detect neutrinoless double beta decay of 130Te, by using Monte Carlo simulations informed by material screening and measurements.

Contribution

It provides a comprehensive background model for CUORE based on simulations and material data, guiding the experiment's construction and sensitivity goals.

Findings

Projected background rate in the region of interest is below 10^{-2} counts/keV/kg/y.

Simulation results inform construction strategies to optimize sensitivity.

Background modeling supports reaching the 9×10^{25} year half-life sensitivity.

Abstract

The Cryogenic Underground Observatory for Rare Events (CUORE) is designed to search for neutrinoless double beta decay of 130Te with an array of 988 TeO2 bolometers operating at temperatures around 10 mK. The experiment is currently being commissioned in Hall A of Laboratori Nazionali del Gran Sasso, Italy. The goal of CUORE is to reach a 90\% C.L. exclusion sensitivity on the \tect decay half-life of 9$\times$10$^{25}$ years after 5\,years of data taking. The main issue to be addressed to accomplish this aim is the rate of background events in the region of interest, which must not be higher than 10$^{-2}$\,counts/keV/kg/y. We developed a detailed Monte Carlo simulation, based on results from a campaign of material screening, radioassays, and bolometric measurements, to evaluate the expected background. This was used over the years to guide the construction strategies of the experiment and we use it here to project a background model for CUORE. In this paper we report the results of our study and our expectations for the background rate in the energy region where the peak signature of neutrinoless double beta decay of 130Te is expected.

Tables

Table 1. Table 1: List of the background sources contributing to the 0 ​ ν ​ β ​ β 0 𝜈 𝛽 𝛽 0\nu\beta\beta  decay ROI of the CUORE-0 experiment (from  Q0_2nu ). The continuum event rate in this region (i.e. excluding the 60 Co sum peak) is 0.058 ± 0.006 plus-or-minus 0.058 0.006 0.058\pm 0.006  counts/keV/kg/y)  Q0-analysis . Column (2) reports the fractional contribution of the different sources, as obtained by the CUORE-0 background reconstruction. The Holder is the structure that held the crystals, made from CuNOSV and reproduced identically in CUORE. The Shields represent both the cryostat external and internal shields used in CUORE-0 to shield the detector from environmental and setup radioactivity (see  Q0_2nu for more details).

Component	Fraction [%]
Shields	74.4	$\pm$	1.3
Holder	21.4	$\pm$	0.7
Crystals	2.64	$\pm$	0.14
Muons	1.51	$\pm$	0.06

Table 2. Table 2: Elements of the CUORE experimental setup implemented in the Monte Carlo simulation, with the values of their masses and surface areas. Elements are grouped into two regions: the near region and the far one. Under the OFE category we included also a 3% of NOSV copper parts located in the far region. The near region includes components close to the TeO 2  detectors. Surface contamination is simulated only for elements in the near region, hence surface areas are listed only for them. The surface area indicated for NOSV copper parts only takes into account surfaces directly facing the detectors. The rightmost column lists the short name of the components used in this paper.

Region	Element	Mass [g]	Surface [cm2]	Short name
Near	TeO2 crystals	7.42$\times$105	148200	TeO2
	Si heaters	6.8	158	heaters
	NTD Ge thermistors	42	288	NTDs
	PEN-Cu cables	389	1200	PEN
	Copper wire pads	57	1140	pads
	PTFE supports	5500	29800	PTFE
	NOSV copper parts	9.4$\times$105	278900	CuNOSV
Far	OFE copper parts	7.0$\times$106	-	CuOFE
	Superinsulation layers	17$\times$103	-	SI
	Roman lead shield	4.5$\times$106	-	RomanPb
	Modern lead shield	2.1$\times$106	-	ModernPb
	Stainless steel rods	15.2$\times$103	-	Rods
	300K steel flange	1.9$\times$106	-	300KFlan

Table 3. Table 3: Values and 90% C.L. (95% C.L. for HPGe measurements) upper limits on 232 Th and 238 U bulk contaminations of detector and cryostat materials (the various components are described in Section  2 ), as obtained in the material screening campaign. Activities are expressed in Bq/kg (1 Bq/kg = 246 × \times 10 -9 g/g for 232 Th and 81 × \times 10 -9 g/g for 238 U). The quoted uncertainties are statistical. In the last column, the measurement technique is indicated, as described in the text. We only include the results shown in bold in the background budget evaluation, as discussed in Section  6 .

BULK CONTAMINATIONS OF MATERIALS USED IN FAR AND NEAR REGIONS
Material	232Th	238U	Technique
	[Bq/kg]	[Bq/kg]
TeO2	8.4$\times$10-7	6.7$\times$10-7	bolometric
Glue	8.9$\times$10-4	1.0$\times$10-2	NAA
Au bonding wires	4.1$\times$10-2	1.2$\times$10-2	ICPMS
heaters	3.3$\times$10-4	2.1$\times$10-3	NAA
NTDs	4.1$\times$10-3	1.2$\times$10-2	producer spec.
PEN	1.0$\times$10-3	1.3$\times$10-3	NAA(Th) + HPGe(U)
PTFE	6.1$\times$10-6	2.2$\times$10-5	NAA
CuNOSV	2.0$\times$10-6	6.5$\times$10-5	NAA + HPGe
CuOFE	6.4$\times$10-5	5.4$\times$10-5	HPGe
RomanPb	4.5$\times$10-5	4.6$\times$10-5	HPGe
ModernPb	1.4$\times$10-4	1.4$\times$10-4	NAA
SI	(11$\pm$2)$\times$10-3	2.4$\times$10-3	HPGe
Rods	(4$\pm$2)$\times$10-4	(8$\pm$2)$\times$10-4	HPGe
300KFlan	(4.5$\pm$0.5)$\times$10-3	(1.6$\pm$0.5)$\times$10-3	HPGe

Table 4. Table 4: 90% C.L. upper limits on surface contamination of various CUORE detector components, as obtained in the material screening campaign. To infer the surface contamination from the measured data, different contamination depths were considered (column 1): the limits reported in the table are the ones corresponding to the contamination depth that gives the highest background contribution, as explained in the text. In the last column, the measurement technique is indicated.

SURFACE CONTAMINATIONS OF MATERIALS USED IN THE NEAR REGION
Material	Depth	232Th	238U	210Pb	Technique
	[$\mu$m]	[Bq/cm2]	[Bq/cm2]	[Bq/cm2]
TeO2	0.01-10	1.9$\times$10-9	8.9$\times$10-9	9.8$\times$10-7	bolometric
heaters	0.1-10	3.3$\times$10-6	8.2$\times$10-7	8.2$\times$10-7	bolometric
NTDs	0.1-10	8.0$\times$10-6	5.0$\times$10-6	4.0$\times$10-5	SiSB
PEN	0.1-30	4.0$\times$10-6	5.0$\times$10-6	3.0$\times$10-5	SiSB
PTFE	0.1-30	1.9$\times$10-8	6.8$\times$10-8	-	NAA
CuNOSV	0.1-10	6.8$\times$10-8	6.5$\times$10-8	8.6$\times$10-7	bolometric

Table 5. Table 5: Values and 90% C.L. upper limits on bulk contaminations in the TeO 2  crystals and the Holder based on the background model of CUORE-0  Q0_2nu . Contaminants are identified as follows: when the progenitor is indicated all the chain is assumed in secular equilibrium, in all other cases single isotopes or sub-chains of 232 Th and 238 U are considered. Each contaminant is assigned a unique index for future reference.

BULK CONTAMINATION of TeO2 and Holder
Material	Source	Index	Activity [Bq/kg]
TeO2	210Po	1	(2.39$\pm$0.11)$\times$10-6
	210Pb	2	(1.37$\pm$0.19)$\times$10-6
	232Th only	3	(7$\pm$3)$\times$10-8
	228Ra-208Pb	4	3.5$\times$10-8
	238U–230Th	5	7.5$\times$10-9
	230Th only	6	(2.8$\pm$0.3)$\times$10-7
	226Ra–210Pb	7	7$\times$10-9
Holder	232Th	18	2.1$\times$10-6
	238U	19	1.2$\times$10-5

Table 6. Table 6: 90% C.L. upper limits and values for surface contamination of the TeO 2  crystals and the Holder based on the background model of CUORE-0  Q0_2nu . Contaminants are identified as follows: when the progenitor is indicated all the chain is assumed in secular equilibrium, in all other cases single isotopes or sub-chains of 232 Th and 238 U are considered. Each contaminant is assigned a unique index for future reference.

SURFACE CONTAMINATION of TeO2 and Holder
Material	Contamination	Index	Activity [Bq/cm2]
TeO2	232Th only .01$\mu$m	8	(3$\pm$1)$\times$10-10
	228Ra-208Pb .01$\mu$m	9	(2.32$\pm$0.12)$\times$10-9
	238U-230Th .01$\mu$m	10	(2.07$\pm$0.11)$\times$10-9
	230Th only .01$\mu$m	11	(1.15$\pm$0.14)$\times$10-9
	226Ra-210Pb .01$\mu$m	12	(3.14$\pm$0.10)$\times$10-9
	210Pb .001$\mu$m	13	(6.02$\pm$0.08)$\times$10-8
	210Pb 1$\mu$m	14	(8.6$\pm$0.8)$\times$10-9
	210Pb 10$\mu$m	15	2.7$\times$10-9
	232Th 10$\mu$m	16	(7.8$\pm$1.4)$\times$10-10
	238U 10$\mu$m	17	3.3 $\times$10-11
Holder	210Pb .01$\mu$m	20	(2.9$\pm$0.4)$\times$10-8
	210Pb .1$\mu$m	21	(4.3$\pm$0.5)$\times$10-8
	210Pb 10$\mu$m	22	1.9 $\times$10-8
	232Th 10$\mu$m	23	(5.0$\pm$1.7)$\times$10-9
	238U 10$\mu$m	24	(1.38$\pm$0.16)$\times$10-8

Table 7. Table 7: 90% C.L. upper limits for surface contamination of the PTFE supports based on the background model of CUORE-0  Q0_2nu with the additional information coming from Fig.  4 .

SURFACE CONTAMINATION of PTFE
Material	Contamination	Activity [Bq/cm2]
PTFE	210Pb 10$\mu$m	7.6 $\times$10-8
	232Th 10$\mu$m	1.5 $\times$10-8
	238U 10$\mu$m	1.9 $\times$10-8

Table 8. Table 8: Main background indexes expected for the various components of CUORE (limits are 90% C.L.).

Region	Source	Background Index
		[counts/keV/kg/y]
External	Environmental $\gamma$	3.6$\times$10-4
	Environmental n	(8$\pm$6)$\times$10-6
	Environmental $\mu$	(1.0$\pm$0.2)$\times$10-4
Far	Rods and 300KFlan: natural radioactivity	(7$\pm$3)$\times$10-5
	SI: natural radioactivity	(2.2$\pm$0.4)$\times$10-4
	ModernPb: natural radioactivity	1.8$\times$10-4
	RomanPb: natural radioactivity	2.9$\times$10-3
	CuOFE: natural radioactivity	1.7$\times$10-3
Near	TeO2 : cosmogenic activation	6.7$\times$10-5
	CuNOSV: cosmogenic activation	2.2$\times$10-6
	CuNOSV: natural radioactivity	(8.7$\pm$0.3)$\times$10-3
	TeO2 : natural radioactivity	(1.2$\pm$0.1)$\times$10-3

Table 9. Table 9: Expected background indexes (statistical errors only) for the different layers of the CUORE tower array. The number of crystals in each layer is reported.

Layer	Number of crystals	Background Index
		[counts/keV/kg/y]
0	528	(0.82$\pm$0.02)$\times$10-2
1	272	(1.17$\pm$0.04)$\times$10-2
2	164	(1.36$\pm$0.04)$\times$10-2
3	24	(1.78$\pm$0.07)$\times$10-2

Equations

$$\rho ={\rho }_0\times \exp (-x/\lambda ),$$

$${ϵ}_i^{MC}=N_{ROI}^i/N_{decays}^i.$$

$${\text{BI}}_i=\frac{A_i\times ({ϵ}_i^{MC}\times {ϵ}^{instr})}{\Delta \times M}$$

Full text

††journal: Eur. Phys. J. C

1 Introduction

CUORE is a ton scale experiment Q_NIMA with the primary physics goal of searching for neutrinoless double beta ($0\nu\beta\beta$) decay of 130Te. Discovery of this phenomenon has been pursued now for several decades DBDReview and, if observed, it would provide crucial evidence for lepton number violation as well as open the door to physics models seeking to explain the matter-antimatter asymmetry in the universe.

Installed in the underground Hall A of Laboratori Nazionali del Gran Sasso (LNGS), Italy, CUORE is currently in the commissioning phase. The experiment is the result of a long standing activity focused on the optimization of single particle thermal detectors (bolometers) based on TeO2 crystals MiDBD ; rad2005 ; rad2006 ; qino-2011 ; qino-120Te ; qino-exc ; ccvr ; Q-axions ; Q-DM ; Q0_PRL . TeO2 -based bolometers have long been used in $0\nu\beta\beta$ decay searches because their properties are well-matched to the requirements of such experiments; they have a very low heat capacity and exhibit extremely good energy resolution, while simultaneously serving as both the source of the $0\nu\beta\beta$ decay and the detector. They also possess low intrinsic background ccvr ; JCG312 , and can be operated stably for several years. Moreover the high natural isotopic abundance of the $0\nu\beta\beta$ decay candidate 130Te (34.17%) isotopic_abundance avoids the needs for expensive isotopic enrichment at this stage.

MiDBD MiDBD , Cuoricino qino-2011 ; qino-2008 and, more recently, CUORE-0 Q0_PRL ; Q0-InitialPerformances marked important milestones in TeO2–based experiments, with successive improvements to $0\nu\beta\beta$ decay sensitivity obtained through better energy resolution, increases in detector mass, and reduction of background. CUORE is the latest step in this series of evolution: with respect to CUORE-0, the active mass is 19 times higher (742 kg of TeO2 or 206 kg of 130Te) and the expected background event rate in the $0\nu\beta\beta$ energy region is about 1/6th, $\sim$10*-2* counts/keV/kg/y. This corresponds to a 90% C.L. sensitivity on the 130Te decay half-life of 9$\times$1025 years in 5 years of exposure BanksKe .

The construction of CUORE required us to overcome a number of challenges, foremost among them was the design of an apparatus simultaneously meeting the stringent cryogenic and radioactivity constraints. This necessitated meticulous material selection, optimization of production and handling protocols, development and validation of special cleaning procedures, and improvements in sensitivity of radioactive assay techniques. A detailed Monte Carlo simulation was used to evaluate the effects of contamination in the different elements of the apparatus. This provided guidelines for radiopurity requirements and indicated where design modifications could mitigate effects from contamination that was unavoidable. In this paper, we discuss the final steps of this work, presenting the background budget of the CUORE experiment, namely the evaluation of the various contributions to the background event rate in the energy region of interest (ROI).

2 The CUORE experiment

CUORE will search for the $0\nu\beta\beta$ decay of 130Te with a close–packed array of 988 TeO2 thermal detectors, operated at a temperature of $\sim$10 mK by means of a custom–made cryostat. The CUORE detectors are arranged in a cylindrical matrix of 19 vertical towers, each one composed of 13 planes of 4 detectors modules supported by a copper frame. The rigidity of the structure is secured by four copper columns that connect each plane to the next one, as shown in Fig. 1. All the copper components of the detector tower are made of NOSV copper, a special copper alloy suitable for cryogenic use produced by Aurubis aurubis .

Each detector is a 750 g TeO2 cubic crystal ($5\times 5\times 5$ cm3) secured inside the copper frame by PTFE supports. These supports are designed to reduce the amount of mechanical stress on the crystals when cooled to cryogenic temperatures, and are the only mechanical parts in contact with the TeO2 crystals. A neutron-transmutation-doped (NTD) Ge thermistor 95Haller is glued with Araldit Rapid Epoxy onto the surface of each crystal to serve as a thermometer. The temperature increase produced by energy deposition due to particle interaction is converted into a voltage signal, which in turn is read out by a specially designed electronics chain. Electrical contacts to the detectors are obtained by means of Au bonding wires. These connect the thermistors to special flexible flat tapes (PEN-Cu cables PENcalbes ) that run vertically along the tower, housed in NOSV copper wire trays. Both the Au bonding wires and the PTFE supports act as thermal links between the crystal–thermometer system (i.e. the bolometer) and the heat bath (i.e. the dilution refrigerator), to restore the bolometers to the operating temperature after each particle interaction. Finally, a small Si heater 98SiHeater ; 12SiHeater is glued onto each crystal and is used to generate reference thermal pulses for the off-line correction of temperature drifts.

The CUORE detector array is operated in vacuum inside a custom-made cryogenic apparatus that complies with very stringent requirements regarding the lowest temperature reached, the mechanical vibration levels, the stability and reliability over long periods, and the radiopurity of all the materials in use. The complete system is shown in Fig. 2.

The cryostat is made up of six cylindrical nested copper vessels—acting as thermal shields—which are thermally anchored to the different temperature stages of the cryogenic system and labeled according to their working temperature: 300 K, 40 K, 4 K, 600 mK, 50 mK, and 10 mK thermal shields. Each vessel is closed on the top by a thick plate; all plates are made of copper but the 300 K one, that is made out of stainless steel111The 300 K vessel is not shown in Fig. 2. Its upper part (about 1/6 of its length) is also made of stainless steel.. All the copper vessels and covers but the 10 mK one are made of OFE copper (C10100 type), an oxygen-free copper alloy suitable for e-beam welding and high-vacuum sealing. The 300 K—Outer Vacuum Chamber (OVC)—and the 4 K—Inner Vacuum Chamber (IVC)—thermal shields are vacuum–tight. Also, the 40 K and the 4 K vessels are surrounded by several layers of super-insulation foil for cryogenic needs (not shown in Fig. 2). The apparatus is a cryogen free cryostat, and the various temperature stages are maintained by a complex cooling system which includes a set of pulse tube cryocoolers and a custom made high-power 3He/4He dilution refrigerator.

The detector towers hang from a thick copper disk—the Tower Support Plate (TSP)—placed inside and thermalized to the 10 mK thermal shield. The TSP is mechanically decoupled from the rest of the experimental setup in order to minimize the transmission of mechanical vibrations to the detectors. Both the TSP and the 10 mK shield are made of NOSV copper. Moreover they are covered, on the side facing the detector array, by thin NOSV copper tiles. This design allowed all the copper parts that face the detector, including the tower skeleton, to be cleaned with the same procedure specially tuned to minimize surface radioactive contaminants (which was not possible for the large and massive parts such as the thermal shield or the TSP).

Two cold lead shields are used to protect the array from the background contributions coming from the cryogenic apparatus. One is a 30 cm thick disk of modern lead together with 6.4 cm of NOSV copper placed just above the TSP inside the 10 mK shield, and thermally anchored to the 50 mK shield, to protect the detector from radioactivity coming from above. The other is a 6 cm layer of ancient Roman lead, thermally anchored to the 4 K vessel, which shields the detector on the side and on the bottom. The Roman lead is mechanically supported by stainless steel bars. Both the modern and Roman lead shields are thermalized by NOSV copper spacers. Outside the cryostat, at room temperature, two shields surround the setup for the abatement of environmental neutron and $\gamma$ fluxes: a lead shield with $\sim$25 cm minimal thickness to absorb environmental $\gamma$ rays, and an outer neutron shield, consisting of 18 cm of polyethylene (to thermalize neutrons) and 2 cm of H3BO3 powder (to capture thermal neutrons).

For the periodic energy calibration of the detector, twelve thoriated tungsten wires are deployed into the cryostat; six are guided into NOSV copper tubes that are placed in fixed positions within the 10 mK volume, and six are guided to the outside of the 50 mK shield. This allows for a near-uniform irradiation of all the crystals with $\gamma$ rays from the 232Th decay chain 17Jeremy .

3 CUORE background sources

If it occurs in nature, $0\nu\beta\beta$ decay can be detected by measuring the energy of the two emitted electrons. Since bolometers are calorimetric detectors, the events in which both electrons are stopped inside the emitting crystal produce a monochromatic peak in the energy spectrum. Particularly, the process 130Te $\rightarrow$ 130Xe $+2\beta^{-}$ is expected to produce (when both electrons are stopped inside the emitting crystal) a line at $Q_{\beta\beta}$ = 2527.518$\pm$ 0.013 keV, the decay transition energy QValueFrank ; QValueNick ; QValueRaman .

We define the 130Te $0\nu\beta\beta$ decay ROI as a 100 keV–wide region in the spectrum centered at $Q_{\beta\beta}$ (2470–2570 keV). A number of sources in addition to $0\nu\beta\beta$ decay can produce events in this energy region, including natural and artificial environmental radioactivity, cosmogenically activated isotopes, and cosmic rays. In most cases these sources produce a continuum in the ROI, that is explicitly chosen to exclude both the 2448 keV line of 214Bi and the 2587 keV Te X-ray escape peak due to the 208Tl line at $\sim$2615 keV. The only peak that, given its proximity to $Q_{\beta\beta}$, cannot be excluded from the ROI is the 60Co sum peak at 2505 keV, due to simultaneous detection in a single crystal of the two 60Co $\gamma$ rays.

The only distinctive signature expected of a $0\nu\beta\beta$ decay event in CUORE is the energy of the signal and its single-hit characteristic. In the assumption of a uniform distribution of $\beta\beta$ decaying nuclei in the crystal, the probability that both emitted electrons are contained inside the single crystal is $\sim$88.4% Q0-analysis . Therefore, in most cases a $0\nu\beta\beta$ decay event results in a signal in only one crystal at a time (single-hit event). On the contrary a large fraction of background events occur in more than one detector at a time (multi-hit event). The rejection of multi-hit events improves the signal-to-noise ratio in the ROI by reducing the background by a larger amount than the signal, and is therefore a powerful active background rejection tool in the CUORE experiment.

The most common background sources for experiments located in underground laboratories are the radioactive contaminants of the construction materials  bkg_EXO ; bkg_EDELWEISS ; bkg_GERDA ; bkg_KZ ; bkg_MAJORANA , namely:

• •

long–lived radioactive nuclei such as 40K, 238U, and 232Th (the latter two are progenitors of radioactive decay chains);

• •

anthropogenic radioactive isotopes (mainly produced by industrial processes or by human–induced nuclear activity in the atmosphere) such as 60Co, 137Cs, and 134Cs;

• •

cosmogenically–produced radioactive isotopes (e.g. 60Co caused by fast nucleon interactions in copper and tellurium);

Another common source of background is the environmental flux of neutrons, muons and $\gamma$ rays at the experimental site. In the case of CUORE, the LNGS average rock overburden of $\sim 3600$ m.w.e. strongly reduces the muon flux CuoreExternal . Moreover, the simultaneous effects of the outer neutron and lead shields, as well as the off-line rejection of multi–hit events, are expected to reduce the background contribution of this source in the ROI to a negligible level, compared to background from radioactive contaminants, as discussed in MuFlux1 ; MuFlux2 ; MuFlux3 . Therefore, in this paper we focus on radioactive contaminants, since they drive the sensitivity of CUORE. A similar condition held also for CUORE precursors. In particular, three dominant sources were identified as contributors to the event rate recorded in the ROI of the MiDBD and Cuoricino experiments qino-2008 ; EPJ2 :

multi-Compton events from the 2615 keV $\gamma$ rays of 208Tl222This is the only environmental gamma line with a B.R. $>$ 1 % and an energy higher than $Q_{\beta\beta}$. originating from 232Th contamination of the cryostat; 2. 2.

contamination by 238U, 232Th, and their decay products (including the surface implantation of 210Pb from environmental 222Rn), on the surface of the TeO2 crystals; 3. 3.

contamination by 238U, 232Th, and their decay products, on the surfaces of inert materials in close proximity to the TeO2 crystals, most likely the copper holder structure.

This result drove most of the strategies adopted for the realization of CUORE Q0-detector including: the minimization of both the amount of material used in the support structure and the space between the crystals (the latter improves the efficiency of background rejection through the * multi-hit* cut); the stringent control of material contamination; and the design of highly efficient shielding from radioactivity. A completely new design was developed to achieve at least a 20-fold reduction of the $\gamma$ ray induced background (source 1 in the previous list) which dominated in the previous experiments. Special cleaning processes were developed to reduce sources 2 and 3. In particular all the TeO2 crystals were processed following a strict radiopurity protocol JCG312 with a final surface polishing which reduced 232Th, 238U, and 210Pb surface contaminants $\sim$6 fold relative to what previously achieved. Copper is by far the most abundant material in the ambient surfaces directly facing the detectors. Therefore, a special cleaning protocol was optimized for all NOSV copper parts constituting the detector holder (frames, columns, wire trays, etc.) as well as for the thin copper plates which surround the whole detector array. The procedure was developed at Laboratori Nazionali di Legnaro, Italy, and includes a sequence of different steps of surface treatments—Tumbling, Electropolishing, Chemical etching, and Magnetron plasma cleaning. This particular methodology, called TECM cleaning, was shown to be the most effective method among those investigated for CUORE TTT . In order to avoid surface recontamination due to airborne radioactivity, which is primarily responsible for 222Rncontamination, particular care was devoted to the production, storage and assembly of all the detector components. These were handled and stored in nitrogen atmosphere and the 19 towers, once assembled, were continuously flushed with nitrogen. The final installation of the CUORE towers in the cryogenic system took place in a dedicated clean-room environment, with continuously circulating radon-depleted atmosphere and constant monitor of the radon level. Details can be found in JCG312 ; 14AssemblyLine ; LNGS15 .

To test the effectiveness of these background reduction strategies, the first tower of TeO2 detectors constructed in the CUORE assembly line 14AssemblyLine was installed in the same cryostat used for the MiDBD and Cuoricino. This tower, named CUORE-0 Q0-detector , was also operated in Hall A of LNGS serving both as a technical prototype for CUORE as well as a stand-alone $0\nu\beta\beta$ decay experiment Q0_PRL . Thanks to the high quality of its data, a very detailed background reconstruction was possible, disentangling the major sources contributing to the ROI event rate. While no change was expected in the background contribution of the cryostat shields because the new tower was operated in the older cryostat, a reduction of the contribution of the other two sources was expected from the improved background abatement protocols. CUORE-0 results proved to be in good agreement with expectations Q0_2nu , as shown in Table 1 (see also Fig. 15 in Q0_2nu ). The observed background was reduced $\sim$2.6-fold relative to Cuoricino demonstrating the effectiveness of the cleaning protocols.

4 CUORE Monte Carlo code

Our simulation is based on the GEANT4 Agostinelli2003250 package; its architecture and its implementation of the physics are identical to the one adopted in the simulation successfully used for CUORE-0 background reconstruction  Q0_2nu .

The simulation code generates and propagates primary and secondary particles through the CUORE geometry until they are detected in the TeO2 crystals. It outputs the energy and time of the energy depositions; the time is used to properly take into account correlations in nuclear decay chains.

We then take the output from GEANT4 and apply a detector response function and other readout features.

Our code is based on the 4.9.6.p03 version of the GEANT4 simulation toolkit, and employees the Livermore Physics List. It includes the propagation of photons, electrons, $\alpha$ particles and heavy ions (nuclear recoils from $\alpha$ emission), as well as neutrons and muons. All primary particles and secondaries are propagated down to keV energies, with a tracking cut optimization inserted in the different detector volumes to balance simulation accuracy and speed (e.g. the tracking cut in lead is set to 1 cm, while in copper to 1 mm). The generation of nuclear transitions is based on an ad hoc implementation of the G4RadioactiveDecay database, performing a concatenation of the tabulated single isotopes decays to correctly simulate chains of radioactive decays in secular equilibrium taking into account their specific time structures.

In the simulation the geometries of the detector, of the cryostat, and of the internal and external shields are reproduced (see Fig. 2), namely: the TeO2 crystals, the copper structure holding the array (i.e. frames and columns), the PTFE supports, the wire trays, the NTD Ge thermistors, the calibration source guiding tubes, the lead support steel bars, the various thermal shields and other cryostat parts, the internal and external lead shields, and the external polyethylene shield. In Table 2 we list the masses, surface areas, and materials corresponding to each element implemented in the Monte Carlo simulation. CuNOSV, henceforth, includes all the NOSV copper components of the detector, i.e. frames, columns, wire–trays, calibration source tubes, 10 mK shield, TSP, all the thin plates covering the TSP and the 10 mK shield, and the two NOSV copper disks enclosing the internal modern lead shield (see Section 2). Similarly, CuOFE indicates all the OFE copper parts of the experimental setup (mainly the cryostat thermal shields—see Section 2).

In our convention, elements are geometrical volumes made of the same material that have a similar production history and therefore similar contamination levels. Each element of the simulation can be studied as an active source whose bulk and/or surface radioactive pollutants are simulated independently in order to evaluate the effect of its contamination, particularly in the ROI. We group elements into two regions (near and far) according to their proximity to the bolometers. This is useful when discussing surface contaminations, as it will become clear in the following.

The distribution of impurities in the bulk of the different materials is assumed to be spatially uniform, while impurities on the surfaces are modeled according to the diffusion process with an exponential density profile

[TABLE]

where $\rho_{0}$ is the impurity density, $x$ is the distance from the surface into the bulk, and $\lambda$ is the mean penetration depth of the impurity.

For reproducing the experimental data from the Monte Carlo simulation, we also model the detector time and energy response. Assuming a detector response similar to that of CUORE-0 detectors Q0-analysis , we account for particle pulses (i.e. pulses following energy deposits in the crystals) with rise times of $\sim 0.05$ s and two decay time components, a fast decay time of $\sim 0.2$ s and a slower one of $\sim 1.5$ s. We account for the timing resolution of each crystal by summing energy depositions that occur in the same crystal within a time window of $\pm$ 5 ms (pulses with longer time distances can be distinguished). Once the simulated events are correctly correlated in time, the resulting energy depositions are smeared with a Gaussian energy response function assuming an average energy resolution of 5 keV FWHM (i.e. the design energy resolution at $Q_{\beta\beta}$ of CUORE and realized in CUORE-0333In the CUORE-0 experiment, the projected energy resolution was $5.1\pm 0.3$ keV FWHM at $Q_{\beta\beta}$.Q0-analysis ). Consistent with what is done in the experiment, we define coincident events as those occurring in multiple bolometers within a $\pm$ 5 ms window. A multiplicity parameter is added to each event, given by the number of bolometers in coincidence for that event. Single-hit events are identified by multiplicity = 1. Finally, as done on CUORE-0 real data Q0-analysis , a pile-up cut of 7.1 s around each event is implemented, rejecting pulses occurring in the same bolometer within 3.1 s before or 4 s after the selected event, and dead–time is accounted for.

5 Material assay

An extensive screening campaign has preceded the selection and procurement of all materials used in the construction of CUORE (for details on material assay techniques commonly adopted in this field see bkgLTR ). The focus was on the most ubiquitous natural contaminants 232Th and 238U (with their respective progenies), and on a few cosmogenically activated contaminants: 60Co, 110Ag and 110mAg EttoreCosmogenic . All these isotopes emit, in their decay chains, particles that have enough energy to mimic a $0\nu\beta\beta$ decay. Even if the mentioned isotopes do not include all the natural or comogenic radionuclei, all other isotopes are expected to yield minor, if not negligible, contributions in the ROI.

The selection of all the CUORE components comprised the certification of the material activity, both for bulk and surface contaminations. Surface contamination can occur during the processes of machining and cleaning, or during exposure to contaminated air. In addition, the effect of breaks in the secular equilibrium of a radioactive chain can be quite relevant for surface contaminations, as it is the case for 210Pb in the 238U chain. Indeed, 222Rn emanation from any material containing 238U impurities produces excess concentrations of 210Pb (the only long-lived isotope in the 222Rn progenies) in the air and dust that, in turn, can contaminate exposed components.

The techniques adopted in the various phases of material screening include $\gamma$ spectroscopy with heavily-shielded High Purity Germanium (HPGe) diodes, to investigate bulk contaminations, and $\alpha$ spectroscopy with large-area low-background Silicon Surface Barrier (SiSB) diodes, to screen for surface contaminations. Additionaly, Inductively Coupled Plasma Mass Spectrometry (ICPMS) and Neutron Activation Analysis (NAA), both particularly suitable for small samples, were also used. Alpha and gamma spectroscopy and ICPMS analysis were performed at LNGS 04Laubenstein ; Heusser-gempi , Milano-Bicocca 13Sala , Baradello Laboratory 04Brofferio and Lawrence Berkeley National Laboratory (LBNL) 13Thomas . Neutron activation analysis was carried out in collaboration with the Laboratory of Applied Nuclear Energy (LENA) in Pavia, Italy NAA1 ; NAA2 .

All materials constituting the CUORE experimental setup, as well as those in contact with the detector components during the production and cleaning phases, have thus been carefully selected according to their bulk contamination levels. Then, for materials in the far region (see Table 2), careful cleaning and storage with standard techniques are enough to ensure that surface contributions are negligible. This was routinely confirmed with SiSB detectors. On the other hand, in the near region the radiopurity requirement for material surfaces can be as strict as a few nBq/cm2. In fact, in elements sufficiently close to the bolometers $\alpha$ and $\beta$ particles can also mimic a $0\nu\beta\beta$ event. For these particles the detection efficiency is much higher (since the absorption probability is nearly 1), and even contaminants with a small activity can produce a relevant background rate in this case. In critical cases when the required sensitivity was not achievable with SiSB detectors, NAA and ICPMS have been exploited to certify material surface contamination.

Finally, in the few cases where all the above techniques failed in reaching the required sensitivity, bulk and/or surface contamination levels were determined through particle spectroscopy with TeO2 detectors ccvr ; EPJ2 ; EPJ1 . These measurements were carried out underground at LNGS using arrays of CUORE-like TeO2 bolometers. One of the most remarkable results obtained in this case is the evaluation of the radioactive sources contaminating the CUORE-0 detector with an unprecedented precision, as will be discussed later.

In the next sections the evaluation of contaminant activities in all the materials used for the CUORE construction will be presented. The evaluation of bulk contaminants will be discussed first, followed by surface ones. The CUORE-0 results, both for bulk and surface contaminations, will instead be discussed in two dedicated sections.

5.1 Cosmogenic activation

Cosmogenic activation is a well-known mechanism for the production of radioactive nuclei in the bulk of materials, mainly through spallation processes. At sea level, cosmic rays are comprised of charged pions, protons, electrons, neutrons, and muons, with relative flux intensities of roughly 1:13:340:480:1420 RadioProtection . Neutrons are clearly the dominant source of cosmogenic activation above ground. At the CUORE experimental site at LNGS, the cosmic ray flux is decreased to an almost negligible level (only muons survive the rock overburden), and hence cosmogenic activation is drastically reduced. Therefore, in this paper we analyze only activation of materials before their underground storage.

In CUORE, the most abundant materials are TeO2, copper, and lead. Among these, only Te and Cu isotopes exhibit large enough cross sections for radioisotope production via cosmogenic activation. In order to contribute to the ROI, cosmogenic nuclei must have a transition energy greater than $Q_{\beta\beta}$, a sizable production cross section, and a relatively long half-life (compared to the time scale of the experiment). Based on the neutron flux at sea-level 04Gordon ; cosmicflux , the exposure time of copper and TeO2 to cosmic rays, and the measured or calculated neutron spallation cross sections act1 ; act2 ; act3 ; act4 ; act5 ; act6 ; act7 ; act8 ; BarbaraCosmogenic ; EttoreCosmogenic , two isotopes fulfill these criteria: 60Co and 110mAg.

60Co ($\tau_{1/2}=5.27$ years, Q = 2.82 MeV) is produced both in copper and tellurium. It $\beta$-decays with the simultaneous emission of two $\gamma$ rays (1.17 MeV and 1.33 MeV). To mimic the energy of a $0\nu\beta\beta$ event, both the $\gamma$ rays and the electron must deposit their energy in the same crystal, requiring 60Co to be either in the TeO2 crystals or in the copper parts close to the bolometers (CuNOSV).

110mAg ($\tau_{1/2}=250$ days) is produced in tellurium. It can either $\beta$-decay to the stable isotope 110Cd (B.R. = $98.7\,\%$, Q = 3.01 MeV), or it can make an isomeric transition to 110Ag (B.R. = $1.3\%$, Q = 0.118 MeV), which then further $\beta$-decays to 110Cd (B.R. = $99.7\,\%$ and Q = 2.89 MeV). 110Ag has a short half-life ($\tau_{1/2}=24.6$ s) and will therefore be in secular equilibrium with 110mAg during the expected 5 years of CUORE data-taking. Both 110mAg and 110Ag emit numerous photons in their decays. The decay of either isotope can mimic a $0\nu\beta\beta$ event when the total energy deposited in a single crystal by the photons and electrons is within the ROI.

To minimize the activation levels of TeO2 and copper (both NOSV and OFE), their exposure to cosmic rays was reduced as much as possible. Approximately 3 months elapsed from crystal growth to underground storage at LNGS for the TeO2 crystals, followed by an average cooling time of 4 years before use. 4 months elapsed from casting of the raw materials to underground storage for all the machined and cleaned CuNOSV components. Requirements were less stringent for OFE copper since this is used in the far region. The activation levels of TeO2 crystals at the start of CUORE have been estimated based on the results of neutron and proton activation measurements published in BarbaraCosmogenic and EttoreCosmogenic . Estimation of activation levels for CuNOSV follows SusanaCosmogenic and incorporates the CUORE-0 results. The estimates of the activation levels are given as conservative upper limits, evaluated approximately one year prior to the beginning of CUORE data taking. These are:

• •

$<20$ nBq/kg of 110mAg + 110Ag (in secular equilibrium with each other) in the *TeO2 * BarbaraCosmogenic ;

• •

$<1$ nBq/kg of 60Co in the *TeO2 * BarbaraCosmogenic . This level of contamination is far below both the HPGe sensitivity and the sensitivity achieved in CUORE-0 (30 $\mu$Bq/kg Q0_2nu );

• •

$<35~{}\mu$Bq/kg of 60Co in the CuNOSV. This limit is in good agreement with the HPGe measurements performed on a few CuNOSV samples ($<25~{}\mu$Bq/kg) and also with the 60Co activity resulting from the same copper444Same copper here refers to copper belonging to the same batch and having the same production history, with the only difference being a longer undergound storage period in the case of the CUORE copper. in CUORE-0.

5.2 238U and 232Th bulk contamination

Table 3 shows the 238U and 232Th bulk activities of the different CUORE materials as obtained in the radioactive assay campaign. The detection efficiency in each measurement is determined using a GEANT4-based Monte Carlo simulation that reproduces the detector geometry and distribution of contaminants in the sample. For each isotope we report the result (or upper limit) of the most sensitive method employed, even if more than one technique may have been used for the same material (i.e. 232Th and 238U results may come from different assay techniques for the same material).

In the case of the bulk contamination of CuNOSV, we report in Table 3 the limits obtained with direct measurements and in Table 5 the limits obtained in the CUORE-0 analysis (the latter are discussed in the next section). As discussed in Section 6, the results obtained with the CUORE-0 detector are used for the CUORE background budget evaluation.

The bulk contamination of the NTDs here refers to the impurity concentration certified by the manufacturer for undoped wafers. The doping is done in a nuclear reactor producing a large number of short-lived radioactive isotopes. The hypothesis that long-lived isotopes potentially dangerous for a $0\nu\beta\beta$ bolometric experiment could also be produced during reactor exposure was investigated in articoloNTD and rejected. After the doping process, an ohmic contact is created on the Ge surface. This operation could cause contamination of the thermistor, that is analyzed as a surface contribution. It is anyway worth noting that signals originating from nuclear decays occurring in the thermistor volume are deformed in their shape, and can therefore be efficiently rejected by the standard pulse shape cuts applied by the analysis process.

5.3 238U, 232Th and 210Pb surface contamination

In Table 4, we report the most sensitive upper limits obtained for the surface activities of the materials used in the near region.

As it was done for the bulk contamination, a Monte Carlo simulation was used to determine, in each measurement, the detection efficiency of the surface contamination. While for bulk contamination the free parameter is only the bulk activity of the sample (always considered to be uniformly distributed in its volume), in the case of surface contamination the impurity distribution is described according to Eq. 1 using two parameters. In all the measurements discussed in this section, experimental data are not enough to measure both $\lambda$ and $\rho_{0}$. Therefore we proceed by evaluating, for each $\lambda$, the contaminant density $\rho_{0}$ compatible with the experimental data. The integral of $\rho$ over the contaminated volume, divided by the surface of the sample yields the impurity concentration measured in Bq/cm2. After having evaluated the surface impurity concentration for a wide range of depths ($\lambda$), we chose the one producing the highest background contribution, and quoted this limit in Table 4. The minimum depth considered in this analysis is 0.001 $\mu$m for *TeO2 * and 0.1 $\mu$m for any other material. Indeed, these are the most shallow distributions whose effects can be experimentally identified. The maximum depth is $\sim$10 $\mu$m for *TeO2 * and CuNOSV (the range of 5 MeV $\alpha$ particles in these materials is $\sim$10 $\mu$m and $\sim$15 $\mu$m, respectively), and 30 $\mu$m for PTFE (where the range of 5 MeV $\alpha$ particles is $\sim$23 $\mu$m). Larger depths are almost indistinguishable from bulk contaminations.

While NTDs and PEN were measured through $\alpha$ spectroscopy with SiSB diodes, in the case of the heaters the small size of the sample required a more sensitive technique. The $\alpha$ spectrometer used in this case was an array of two 5$\times$5$\times$5 cm3 TeO2 bolometers operated in the cryogenic facility of Hall C at LNGS. A matrix of heaters, 5$\times$5 cm2 in total area, was oriented towards the crystals and the $\alpha$ induced background was analyzed with the technique illustrated in rad2006 . In the case of 238U and 232Th impurities in PTFE, the best upper limits are obtained with the NAA technique (analyzing the results as if all measured contaminants were contained in a surface layer). Giving information only on the progenitor, such technique is insensitive to 210Pb contamination. A bolometric detector array (TTT, for Three Tower Test) was used to compare different copper surface treatments and to analyze the CuNOSV contamination level. The TTT detector TTT consisted of three small towers, each with 12 TeO2 detectors and enclosed inside its own copper box. The limits reported in Table 4 are those obtained with the TECM cleaning, wich showed the lowest background level and was therefore chosen as the baseline cleaning protocol for all the CuNOSV pieces of the CUORE detector.

5.4 Contamination from CUORE-0 analysis

In Q0_2nu , the sources contributing to the CUORE-0 event rate were reconstructed by fitting a number of Monte Carlo simulations to the measured spectra. The fit was performed with a Bayesian approach which allows to exploit any previous knowledge on material contamination by defining priors on source activities. For many sources, a sensitivity on contaminant concentration better than that achieved with standard techniques was obtained, along with a more efficient disentanglement of contamination species and a detailed study of secular equilibrium violations in radioactive chains.

Table 5 and Table 6 summarize the results obtained for the impurity concentrations in materials that belong to the same production batches of CUORE, i.e. the TeO2 crystals and the NOSV copper (that in CUORE-0 analysis was taken as representative of the whole detector holder structure —see later). For the bulk contamination levels, the sensitivity improvements obtained by the CUORE-0 analysis are minor (compare the limits reported for TeO2 in Table 3 and 5 or those reported for CuNOSV in Table 3 and for Holder in Table 5). On the contrary, for surface contaminations the improvements are significant: the CUORE-0 analysis allowed to derive a model for the impurity density profiles meanwhile achieving a high sensitivity in the determination of their activity. Moreover, since the procedure adopted in CUORE-0 background reconstruction was to simultaneously fit all the simulated spectra to the measured one, it was possible to find the best evaluation for each contaminant concentration as well as the correlation factor among different sources (see Fig. 3).

A separate discussion applies to the small parts used to build the CUORE-0 detector — i.e. NTD thermistors, PEN cables, Si heaters, Au bonding wires, glue and PTFE supports — also belonging to the same production batch of CUORE. In fact, in CUORE-0 analysis these small mass/surface components exhibit spectra that are completely degenerate with those of the NOSV copper detector structure; therefore only the latter was considered in the model under the name of Holder (see Table 5 and Table 6), chosen to underline that the resulting activities also include the small parts contribution. This contribution is expected to be negligible in the CUORE-0 background, as confirmed by Fig. 4. In these plots the experimental rate of the different planes of the CUORE-0 detector in the energy region between 2.7 MeV and 3.9 MeV (i.e. the interval dominated by degraded $\alpha$ contributions from surface contaminants) is compared to the rate expected from a 210Pb  contamination in the NOSV copper detector structure or in the PTFE supports (chosen as representative of all the small parts, because of their bigger mass/surface) —similar plots are obtained for 232Th  and 238U  contaminants. The degeneracy between these two sources can be broken by examining the detector plane dependence of the counting rate. This is illustrated in Fig. 4. The detectors in the upper and lower planes of the CUORE-0 tower (1+13 floors in Fig. 4) face a larger copper surface of the 10 mK shield, while seeing a slightly smaller surface of PTFE and other small parts compared to the middle floors of CUORE-0. The CUORE-0 data shown in Fig. 4 suggest that the contribution to the overall background in the 2.7-3.9 MeV region from PTFE supports is compatible with zero. The resulting PTFE supports surface contaminations are reported in Table 7 (to be compared to the corresponding ones in Table 4).

6 Background budget

In this section we discuss contributions from material radioactivity to the CUORE background. From now on, the event rate in the ROI will be denoted as the Background Index, BI, and its estimation will include also the instrumental efficiency (discussed later), in order to properly represent the experimental event rate expected in CUORE.

Following the discussion of the previous sections, the sources that give sizable contribution to the BI are:

(near region) 238U, 232Th and their progenies in TeO2 and Holder —the latter as representative of both the CuNOSV and the small parts. For both bulk and surface activities of these sources we use CUORE-0 results as reported in Tables 5 and 6 along with the respective correlation factors; 2. 2.

(near region) cosmogenically activated isotopes in TeO2 and CuNOSV, using the activities discussed in Section 5.2; 3. 3.

(far region) 238U and 232Th in CuOFE, RomanPb, ModernPb, SI, Rods, 300KFlan (refer to Table 2). Activities are those reported in Table 3.

The BI is evaluated as follows. For each source $i$ (i.e. a definite detector element contaminated by a specific contaminant) we run a MC simulation and project the simulated spectra for the 988 bolometers. As we plan to do with the experimental data, we enforce an anticoincidence cut (i.e. we reject multi-hit events). The 988 spectra are finally summed and the source efficiency ($\epsilon_{i}^{MC}$) is evaluated. $\epsilon_{i}^{MC}$ is defined as the probability for the source $i$ to produce an event in the ROI of the sum spectrum; it is therefore given by the ratio between the number of single-hit events in the ROI ($N^{i}_{ROI}$) and the number of simulated decays ($N^{i}_{decays}$)555In the case of 60Co cosmogenic activation in CuNOSV, the 2505 keV peak contribution is subtracted.:

[TABLE]

In the case of sources with uncorrelated activities (points 2 and 3 in the previous list), individual contributions to the BI are given by:

[TABLE]

where $A_{i}$ is the activity of the source $i$, $\Delta$ is the ROI width (100 keV) and $M$ is the TeO2 mass. The factor $\epsilon_{i}^{MC}\times\epsilon^{instr}$ measures the probability of observing a contribution to the measured ROI event rate caused by source $i$. It is a product of conditional probabilities Q0-analysis : the probability for the source $i$ to produce an event in the ROI ($\epsilon_{i}^{MC}$), the probability that this event is triggered and properly reconstructed ($\epsilon^{tr}$), the probability that this event is not accidentally in coincidence with an unrelated event in a different bolometer ($\epsilon^{acc}$) and the probability that it passes the pulse shape cuts ($\epsilon^{PSA}$). These cuts are used to remove events that are either non-signal-like or are in some way not handled well by the data processing software Q0-analysis .

$\epsilon^{tr}$ and $\epsilon^{PSA}$ depend on data acquisition and analysis, and we expect them to be similar to CUORE-0 ones Q0-analysis , $\epsilon^{tr}$= 98.5% and $\epsilon^{PSA}$= 93.7%. $\epsilon^{acc}$ depends on the probability of accidental coincidences and is evaluated to be 99% (as expected for a $\pm$ 5 ms coincidence window and an average signal rate of 1 mHz/bolometer). The product of $\epsilon^{acc}$, $\epsilon^{tr}$ and $\epsilon^{PSA}$ is source independent (we assume here all efficiencies to be energy independent), and we will refer to it as instrumental efficiency, $\epsilon^{instr}$= 91.4%.

The BIs computed following this procedure have a statistical uncertainty that is obtained propagating the statistical uncertainty of Ai and a systematic error that derives both from Ai (material activity measurements have systematics of the order of $\pm$5%) and from $\epsilon_{i}^{MC}$ ($\pm$5%, due to approximations in the description of geometry and detector response to particle interaction).

In the case of the activities extrapolated from CUORE-0 analysis (sources at point 1 in the previous list), the determination of the background index follows a slightly different procedure. Fig. 3 shows that we have two groups of sources,* TeO2* contaminants and Holder contaminants, with strong correlations within each group. Thus we define two cumulative Background Indexes, BI${}_{TeO_{2}}$ and BIHolder, and we associate to both of them a probability density function (pdf), obtained as follows. The joint pdf describing contaminant activities, as derived by the CUORE-0 background fit reconstruction, is sampled. For each sampled point, the activity of the $i^{th}$ source is weighted by its respective CUORE efficiency $\epsilon_{i}^{MC}\times\epsilon^{instr}$ and the sum of the activities induced in the ROI by sources belonging to the same group is calculated, thus taking into account the correlations. The pdf’s associated with BI${}_{TeO_{2}}$ and BIHolder are shown in Fig. 5. In both cases, bulk contamination gives a minor contribution (1–2 orders of magnitude lower) with respect to surface one, which represents therefore the dominant source. The mean of each pdf is taken as the estimation of the BI and the RMS is used to compute its statistical uncertainty.

In the adopted procedure, BIHolder is computed under the hypothesis of neglecting the small parts contribution to the ROI rate. This assumption, though completely acceptable in the case of CUORE-0 analysis (see discussion in Section 5.4), may not be thoroughly valid in CUORE because of the different scaling of the geometrical efficiencies of CuNOSV and small parts going from CUORE-0 to CUORE. As a matter of fact, the various towers of the CUORE detector are facing the same amount of small parts surface but different amounts of CuNOSV surface, depending on their position in the array with respect to the CuNOSV thermal shield of the cryostat. Though small parts contribution is compatible with zero in the case of CUORE-0 model, their introduction in CUORE model adds an additional +15% systematic error to BIHolder. Therefore, the final systematic errors are $\pm$10% for BI${}_{TeO_{2}}$ and ${}_{-10\%}^{+25\%}$ for BIHolder.

Table 8 and Fig. 6 summarize the final results: in both we report also the contribution expected from external sources (“External”) as derived in CuoreExternal .

The total projected BI in the ROI of CUORE is equal to [1.02$\pm$0.03(stat)${}^{+0.23}_{-0.10}$(syst)]$\times$10*-2* counts/keV/kg/y. This number is obtained by summing all the contributions incompatible with zero listed in Table 8 (i.e. we exclude sources for which we have an upper limit). The BI is by far dominated by the Holder contribution, mainly ascribed to degraded alpha particles from surface contaminants (see Q0_2nu ).

The CUORE predicted BI can be compared to the BIs measured by precursor experiments (in units of counts/keV/kg/y):

• $\bullet$

$0.058\pm 0.004$ (stat.) $\pm$ 0.002 (syst.) ($\epsilon^{instr}$= 92.0%) in CUORE-0;

• $\bullet$

$0.153\pm 0.006$  ($\epsilon^{instr}$= 94.7%) in Cuoricino.

Finally, we report in Table  9 the BIs in different subgroups of detectors (“layers”), grouped by the number of crystal sides directly facing the CuNOSV thermal shield. Layer-0 contains those crystals with no sides facing the CuNOSV shield (i.e. TeO2 crystal fully shielded by other detectors), Layer-1 includes crystals with only one side facing the CuNOSV shield, and so on. This grouping reflects the fact that contaminations on the CuNOSV surface are the most significant contribution to the background in the ROI. As a result, we expect different layers to have different projected sensitivities to the $0\nu\beta\beta$ half-life of 130Te, resulting in a global improvement in the final CUORE sensitivity.

The sources analyzed in this study are expected to be the major contributors to CUORE background index. Therefore, the results of this work can be used for a tentative prediction of the CUORE spectrum over a wide energy range. Fig. 7 shows the anticoincidence spectrum recorded by the 988 bolometer array predicted by this analysis. Among the contaminants analyzed in this paper, we summed only those that have activities incompatible with zero (i.e. those used for the BI evaluation), with the exception of 60Co in the CuNOSV, which is fixed at its 90% C.L. upper limit. Indeed, even if it does not contribute to the BI, it is likely that the 2505 keV peak will appear in the ROI as it did for CUORE-0, although its activity in CUORE has a large uncertainty. Three additional sources are considered because of their high activity: 130Te $2\nu\beta\beta$ decay and 40K in TeO2 and CuNOSV bulk, all fixed at the CUORE-0 results Q0_2nu . Due to the lack of information on 40K activities of other materials, the 1460 keV 40K peak could be underestimated.

The right panel of Fig. 7 shows the predicted spectrum in the region around the ROI. Two peaks are visible, one from 60Co at 2505 keV (discussed earlier) and one from 208Tl at 2615 keV.

7 Conclusion

In this paper the different sources contributing to the CUORE background figure of merit are analyzed in detail. The study exploit the available data on bulk and surface contamination of CUORE construction materials and is based on the final design of the experiment.

All possible background sources are investigated by means of accurate dedicated measurements, and using the most sensitive available techniques with realistic experimental setup. A specially developed simulation code allows us to calculate the effects of the various background sources and to compare them with the required experimental sensitivity.

A dedicated R&D program (for material selection or improvement of the measurement sensitivity) was initiated each time an experimental result was not consistent with the corresponding design sensitivity.

The equivalence of the CUORE and the CUORE-0 detector tower structure and materials allowed us to fully exploit the results of the CUORE-0 background model reference fit Q0_2nu as input for the activity of all the common sources (mainly TeO2 and CuNOSV bulk and surface contaminations).

Results of the background model developed in this work are reported in Table 8. None of the bulk contaminations appear worrisome for the CUORE target background rate of 10*-2* counts/keV/kg/y. Despite being more serious, surface contaminations of detector components (the detectors themselves and the materials directly surrounding them), are nevertheless within the CUORE goals for the background counts. The final estimate for the BI is [1.02$\pm$0.03(stat)${}^{+0.23}_{-0.10}$(syst)]$\times$10*-2* counts/keV/kg/y.

In conclusion, as shown in Fig. 6, the work and processes involving material selection and surface cleaning were successful in yielding a result consistent with the background budget goal of the CUORE experiment.

8 Acknowledgments

The CUORE Collaboration thanks the directors and staff of the Laboratori Nazionali del Gran Sasso and the technical staff of our laboratories. This work was supported by the Istituto Nazionale di Fisica Nucleare (INFN); the National Science Foundation under Grant Nos. NSF-PHY-0605119, NSF-PHY-0500337, NSF-PHY-0855314, NSF-PHY-0902171, NSF-PHY-0969852, NSF-PHY-1307204, NSF-PHY-1314881, NSF-PHY-1401832, and NSF-PHY-1404205; the Alfred P. Sloan Foundation; the University of Wisconsin Foundation; and Yale University. This material is also based upon work supported by the US Department of Energy (DOE) Office of Science under Contract Nos. DE-AC02-05CH11231, DE-AC52-07NA27344, and DE-SC0012654; and by the DOE Office of Science, Office of Nuclear Physics under Contract Nos. DE-FG02-08ER41551 and DE-FG03-00ER41138. This research used resources of the National Energy Research Scientific Computing Center (NERSC).