The projected background for the CUORE experiment

The Cryogenic Underground Observatory for Rare Events (CUORE) is designed to search for neutrinoless double beta decay of 130\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{130}$$\end{document}Te with an array of 988 TeO2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_2$$\end{document} 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 130\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{130}$$\end{document}Te decay half-life of 9 ×\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\times $$\end{document} 1025\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{25}$$\end{document} 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\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{-2}$$\end{document} counts/keV/kg/year. 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 130\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{130}$$\end{document}Te is expected.


Introduction
CUORE is a ton scale experiment [1] with the primary physics goal of searching for neutrinoless double beta (0νββ) decay of 130 Te. Discovery of this phenomenon has been pursued now for several decades [2] 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 TeO 2 crystals [3][4][5][6][7][8][9][10][11][12]. T eO 2 -based bolometers have long been used in 0νββ decay searches because their properties are wellmatched to the requirements of such experiments; they have a e-mail: cuore-spokesperson@lngs.infn.it b Deceased c Presently at: INFN-Laboratori Nazionali di Frascati, Frascati, 00044 Rome, Italy a very low heat capacity and exhibit extremely good energy resolution, while simultaneously serving as both the source of the 0νββ decay and the detector. They also possess low intrinsic background [9,13], and can be operated stably for several years. Moreover the high natural isotopic abundance of the 0νββ decay candidate 130 Te (34.17%) [14] avoids the needs for expensive isotopic enrichment at this stage. MiDBD [3], Cuoricino [6,15] and, more recently, CUORE-0 [12,16] marked important milestones in T eO 2based experiments, with successive improvements to 0νββ 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 TeO 2 or 206 kg of 130 Te) and the expected background event rate in the 0νββ energy region is about 1/6th, ∼10 −2 counts/keV/kg/year. This corresponds to a 90% C.L. sensitivity on the 130 Te decay half-life of 9 × 10 25 years in 5 years of exposure [17].
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).

The CUORE experiment
CUORE will search for the 0νββ decay of 130 Te with a closepacked array of 988 TeO 2 thermal detectors, operated at a temperature of ∼10 mK by means of a custom-made cryostat. The CUORE detectors are arranged in a cylindrical matrix Two copper frames (grey) joined by four columns (green) form the mechanical structure that secures four crystals by means of polytetrafluoroethylen (PTFE) supports (yellow) 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 [18]. Each detector is a 750 g TeO 2 cubic crystal (5 × 5 × 5 cm 3 ) 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 TeO 2 crystals. A neutron-transmutationdoped (NTD) Ge thermistor [19] 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 [20]) 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 [21,22] 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 vesselsacting as thermal shields -which are thermally anchored to the different temperature stages of the cryogenic system and labeled according to their working temperature: 300, 40, 4 K, 600, 50, 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 steel. 1 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 Kouter 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 3 He/ 4 He 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 γ fluxes: a lead shield with ∼25 cm minimal thickness to absorb environmental γ rays, and an outer neutron shield, consisting of 18 cm of polyethylene (to thermalize neutrons) and 2 cm of H 3 BO 3 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 γ rays from the 232 Th decay chain [23].

CUORE background sources
If it occurs in nature, 0νββ 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 130 Te → 130 Xe +2β − is expected to produce (when both electrons are stopped inside the emitting crystal) a line at Q ββ = 2527.518 ± 0.013 keV, the decay transition energy [24][25][26].
We define the 130 Te 0νββ decay ROI as a 100 keV-wide region in the spectrum centered at Q ββ (2470-2570 keV). A number of sources in addition to 0νββ 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 214 Bi and the 2587 keV Te X-ray escape peak due to the 208 Tl line at ∼2615 keV. The only peak that, given its proximity to Q ββ , cannot be excluded from the ROI is the 60 Co sum peak at 2505 keV, due to simultaneous detection in a single crystal of the two 60 Co γ rays.
The only distinctive signature expected of a 0νββ decay event in CUORE is the energy of the signal and its singlehit characteristic. In the assumption of a uniform distribution of ββ decaying nuclei in the crystal, the probability that both emitted electrons are contained inside the single crystal is ∼88.4% [27]. Therefore, in most cases a 0νββ decay event results in a signal in only one crystal at a time (singlehit 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 signalto-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 [28][29][30][31][32], namely: -long-lived radioactive nuclei such as 40 K, 238 U, and 232 Th (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 60 Co, 137 Cs, and 134 Cs; -cosmogenically-produced radioactive isotopes (e.g. 60 Co caused by fast nucleon interactions in copper and tellurium); Another common source of background is the environmental flux of neutrons, muons and γ rays at the experimental site.
In the case of CUORE, the LNGS average rock overburden of ∼3600 m.w.e. strongly reduces the muon flux [33]. 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 [34][35][36]. 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 [15,37]: 1. multi-Compton events from the 2615 keV γ rays of 208 Tl 2 originating from 232 Th contamination of the cryostat; 2. contamination by 238 U, 232 Th, and their decay products (including the surface implantation of 210 Pb from environmental 222 Rn), on the surface of the TeO 2 crystals; 3. contamination by 238 U, 232 Th, and their decay products, on the surfaces of inert materials in close proximity to the TeO 2 crystals, most likely the copper holder structure.
This result drove most of the strategies adopted for the realization of CUORE [38] 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 γ 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 TeO 2 crystals were processed following a strict radiopurity protocol [13] with a final surface polishing which reduced 232 Th, 238 U, and 210 Pb surface contaminants ∼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 Table 1 List of the background sources contributing to the 0νββ decay ROI of the CUORE-0 experiment (from [42]). The continuum event rate in this region (i.e. excluding the 60 Co sum peak) is 0.058 ± 0.006 counts/keV/kg/year) [27]. 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 [42]  to be the most effective method among those investigated for CUORE [39]. In order to avoid surface recontamination due to airborne radioactivity, which is primarily responsible for 222 Rncontamination, 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 [13,40,41]. To test the effectiveness of these background reduction strategies, the first tower of TeO 2 detectors constructed in the CUORE assembly line [40] was installed in the same cryostat used for the MiDBD and Cuoricino. This tower, named CUORE-0 [38], was also operated in Hall A of LNGS serving both as a technical prototype for CUORE as well as a stand-alone 0νββ decay experiment [12]. 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 [42], as shown in Table 1 (see also Fig. 15 in [42]). The observed background was reduced ∼2.6-fold relative to Cuoricino demonstrating the effectiveness of the cleaning protocols.

CUORE Monte Carlo code
Our simulation is based on the GEANT4 [43] package; its architecture and its implementation of the physics are iden- 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 tical to the one adopted in the simulation successfully used for CUORE-0 background reconstruction [42]. The simulation code generates and propagates primary and secondary particles through the CUORE geometry until they are detected in the TeO 2 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, α particles and heavy ions (nuclear recoils from α 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 TeO 2 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 Sect. 2). Similarly, CuOFE indicates all the OFE copper parts of the experimental setup (mainly the cryostat thermal shields -see Sect. 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 where ρ 0 is the impurity density, x is the distance from the surface into the bulk, and λ 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 [27], we account for particle pulses (i.e. pulses following energy deposits in the crystals) with rise times of ∼0.05 s and two decay time components, a fast decay time of ∼0.2 s and a slower one of ∼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 ±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 ββ of CUORE and realized in CUORE-0 3 [27]). Consistent with what is done in the experiment, we define coincident events as those occurring in multiple bolometers within a ±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 [27], a pileup 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.

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 [44]). The focus was on the most ubiquitous natural contaminants 232 Th and 238 U (with their respective progenies), and on a few cosmogenically activated contaminants: 60 Co, 110 Ag and 110m Ag [45]. All these isotopes emit, in their decay chains, particles that have enough energy to mimic a 0νββ 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 210 Pb in the 238 U chain. Indeed, 222 Rn emanation from any material containing 238 U impurities produces excess concentrations 3 In the CUORE-0 experiment, the projected energy resolution was 5.1 ± 0.3 keV FWHM at Q ββ . of 210 Pb (the only long-lived isotope in the 222 Rn progenies) in the air and dust that, in turn, can contaminate exposed components.
The techniques adopted in the various phases of material screening include γ spectroscopy with heavily-shielded high purity germanium (HPGe) diodes, to investigate bulk contaminations, and α spectroscopy with large-area lowbackground 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 [46,47], Milano-Bicocca [48], Baradello Laboratory [49] and Lawrence Berkeley National Laboratory (LBNL) [50]. Neutron activation analysis was carried out in collaboration with the Laboratory of Applied Nuclear Energy (LENA) in Pavia, Italy [51,52].
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/cm 2 . In fact, in elements sufficiently close to the bolometers α and β particles can also mimic a 0νββ 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 TeO 2 detectors [9,37,53]. These measurements were carried out underground at LNGS using arrays of CUORE-like TeO 2 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.

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 [54]. 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 TeO 2 , 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 ββ , a sizable production cross section, and a relatively long halflife (compared to the time scale of the experiment). Based on the neutron flux at sea-level [55,56], the exposure time of copper and TeO 2 to cosmic rays, and the measured or calculated neutron spallation cross sections [45,[57][58][59][60][61][62][63][64][65], two isotopes fulfill these criteria: 60 Co and 110m Ag. 60 Co (τ 1/2 = 5.27 years, Q = 2.82 MeV) is produced both in copper and tellurium. It β-decays with the simultaneous emission of two γ rays (1.17 MeV and 1.33 MeV). To mimic the energy of a 0νββ event, both the γ rays and the electron must deposit their energy in the same crystal, requiring 60 Co to be either in the TeO 2 crystals or in the copper parts close to the bolometers (CuNOSV ).
110m Ag (τ 1/2 = 250 days) is produced in tellurium. It can either β-decay to the stable isotope 110 Cd (B.R. = 98.7 %, Q = 3.01 MeV), or it can make an isomeric transition to 110 Ag (B.R. = 1.3%, Q = 0.118 MeV), which then further β-decays to 110 Cd (B.R. = 99.7 % and Q = 2.89MeV). 110 Ag has a short half-life (τ 1/2 = 24.6 s) and will therefore be in secular equilibrium with 110m Ag during the expected 5 years of CUORE data-taking. Both 110m Ag and 110 Ag emit numerous photons in their decays. The decay of either isotope can mimic a 0νββ 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 TeO 2 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 TeO 2 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 TeO 2 crystals at the start of CUORE have been esti-mated based on the results of neutron and proton activation measurements published in [65] and [45]. Estimation of activation levels for CuNOSV follows [66] and incorporates the CUORE-0results. 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 110m Ag + 110 Ag (in secular equilibrium with each other) in the T eO 2 [65]; -< 1 nBq/kg of 60 Co in the T eO 2 [65]. This level of contamination is far below both the HPGe sensitivity and the sensitivity achieved in CUORE-0 (30 µBq/kg [42]); -< 35 µBq/kg of 60 Co in the CuNOSV. This limit is in good agreement with the HPGe measurements performed on a few CuNOSV samples (< 25 µBq/kg) and also with the 60 Co activity resulting from the same copper 4 in CUORE-0. Table 3 shows the 238 U and 232 Th 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. 232 Th and 238 U results may come from different assay techniques for the same material).

238 U and 232 Th bulk contamination
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 Sect. 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νββ bolometric experiment could also be produced during reactor exposure was investigated in [67] 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 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 × 10 −9 g/g for 232 Th and 81 × 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 Sect. 6  their shape, and can therefore be efficiently rejected by the standard pulse shape cuts applied by the analysis process.

238 U, 232 Th and 210 Pb 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 λ and ρ 0 . Therefore we proceed by evaluating, for each λ, the contaminant density ρ 0 compatible with the experimental data. The integral of ρ over the contaminated volume, divided by the surface of the sample yields the impurity concentration measured in Bq/cm 2 . After having evaluated the surface impurity concentration for a wide range of depths (λ), 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 µm for T eO 2 and 0.1 µm for any other material. Indeed, these are the most shallow distributions whose effects can be experimentally identified. The maximum depth is ∼10 µm for T eO 2 and CuNOSV (the range of 5 MeV α particles in these materials is ∼10 µm and ∼15 µm, respectively), and 30 µm for PTFE (where the range of 5 MeV α particles is ∼23 µm). Larger depths are almost indistinguishable from bulk contaminations.
While NTDs and PEN were measured through α spectroscopy with SiSB diodes, in the case of the heaters the small size of the sample required a more sensitive technique. The α spectrometer used in this case was an array of two 5 × 5 × 5 cm 3 TeO 2 bolometers operated in the cryogenic facility of Hall C at LNGS. A matrix of heaters, 5 × 5 cm 2 in total area, was oriented towards the crystals and the α induced background was analyzed with the technique illustrated in [5]. In the case of 238 U and 232 Th 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 210 Pb 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 [39] consisted of three small towers, each with 12 TeO 2 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.

Contamination from CUORE-0 analysis
In [42], 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. Tables 5 and 6 summarize the results obtained for the impurity concentrations in materials that belong to the same production batches of CUORE, i.e. the TeO 2 crystals and the NOSV copper (that in CUORE-0 analysis was taken as representative of the whole detector holder structuresee later). For the bulk contamination levels, the sensitivity improvements obtained by the CUORE-0 analysis are minor (compare the limits reported for T eO 2 in Tables 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 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 [42]. 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   Tables 5 and 6 with their Index numbering 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 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 [42]. 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 T eO 2  Holder (see Tables 5, 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 and 3.9 MeV (i.e. the interval dominated by degraded α contributions from surface contaminants) is compared to the rate expected from a 210 Pb 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 232 Th and 238 U 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).

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 instrumen-tal 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) 238 U, 232 Th and their progenies in T eO 2 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. (near region) cosmogenically activated isotopes in T eO 2 and CuNOSV, using the activities discussed in Sect. 5.2; 3. (far region) 238 U and 232 Th in CuOFE, RomanPb, Mod-ernPb, 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 ( MC i ) is evaluated. MC i 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 RO I ) and the number of simulated decays (N i decays ): 5

MC i
In the case of sources with uncorrelated activities (points 2 and 3 in the previous list), individual contributions to the BI are given by: where A i is the activity of the source i, is the ROI width (100 keV) and M is the TeO 2 mass. The factor MC i × 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 [27]: the probability for the source i to produce an event in the ROI ( MC i ), the probability that this event is triggered and properly reconstructed ( tr ), the probability that this event is not accidentally in coincidence with an unrelated event in a different bolometer ( acc ) and the probability that it passes the pulse shape cuts ( P S A ). 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 [27]. tr and P S A depend on data acquisition and analysis, and we expect them to be similar to CUORE-0 ones [27], tr = 98.5% and P S A = 93.7%. acc depends on the probability of accidental coincidences and is evaluated to be 99% (as expected for a ±5 ms coincidence window and an average signal rate of 1 mHz/bolometer). The product of acc , tr and P S A is source independent (we assume here all efficiencies to be energy independent), and we will refer to it as instrumental efficiency, instr = 91.4%.
The BIs computed following this procedure have a statistical uncertainty that is obtained propagating the statistical uncertainty of A i and a systematic error that derives both from A i (material activity measurements have systematics of the order of ±5%) and from MC i (±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. Figure 3 shows that we have two groups of sources, T eO 2 contaminants and Holder contaminants, with strong correlations within each group. Thus we define two cumulative Background Indexes, BI T eO 2 and BI Holder , 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 MC i × 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 T eO 2 and BI Holder 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, BI Holder 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 Sect. 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 BI Holder .  Therefore, the final systematic errors are ±10% for BI T eO 2 and +25% −10% for BI Holder . Table 8 and Fig. 6 summarize the final results: in both we report also the contribution expected from external sources ("External") as derived in [33].
The total projected BI in the ROI of CUORE is equal to [1.02 ± 0.03(stat) +0. 23 −0.10 (syst)] × 10 −2 counts/keV/kg/year. 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 [42]).
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. TeO 2 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νββ half-life of 130 Te, 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. Figure 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 60 Co 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: 130 Te 2νββ decay and 40 K in T eO 2 and CuNOSV bulk, all fixed at the CUORE-0 results [42]. Due to the lack of information on 40 K activities of other materials, the 1460 keV 40 K 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 60 Co at 2505 keV (discussed earlier) and one from 208 Tl at 2615 keV.

Conclusion
In this paper the different sources contributing to the CUORE background figure of merit are analyzed in detail. The study CuNOSV: cosmogenic activation (   Fig. 7 Anticoincidence spectrum resulting from the sources analyzed in this paper (histograms corresponding to the different sources are stacked one over the other, the dotted contour shows the total predicted spectrum). Only sources with activities incompatible with zero are considered, with the exception of the 60 Co contamination in CuNOSV (reported at its 90% C.L. upper limit). Due to their expected large con-tribution to the shape of the spectrum, the 40 K in T eO 2 and in CuNOSV as well as the 130 Te 2νββ decay contributions are also included. The region around the ROI (i.e. the interval 2470-2570 keV) is shown in the right panel. Two peaks are visible, they are: the 60 Co line at 2505 keV line and the 208 Tl line at 2615 keV (see text for more details) 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 [42] as input for the activity of all the common sources (mainly T eO 2 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/year. 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 ± 0.03(stat) +0. 23 −0.10 (syst)] × 10 −2 counts/keV/kg/year. 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.