arxiv: v2 [physics.ins-det] 14 Aug 2017

Transcription

1 Eur. Phys. J. C manuscript No. will be inserte by the etor CUORE Sensitivity to ββ Decay arxiv: v2 [physics.ins-et] 14 Aug 2017 C. Aluino 1, K. Alfonso 2, D. R. Artusa 1,3, F. T. Avignone III 1, O. Azzolini 4, T. I. Banks 5,6, G. Bari 7, J.W. Beeman 8, F. Bellini 9,10, G. Benato 5, A. Bersani 11, M. Biassoni 12,13, A. Branca 14, C. Brofferio 12,13, C. Bucci 3, A. Camacho 4, A. Caminata 11, L. Canonica 15,3, X. G. Cao 16, S. Capelli 12,13, L. Cappelli 3, L. Carbone 13, L. Carani 10, P. Carniti 12,13, N. Casali 9,10, L. Cassina 12,13, D. Chiesa 12,13, N. Chott 1, M. Clemenza 12,13, S. Copello 17,11, C. Cosmelli 9,10, O. Cremonesi 13, R. J. Creswick 1, J. S. Cushman 18, A. D Aabbo 3, I. Dafinei 10, C. J. Davis 18, S. Dell Oro 3,19, M. M. Deninno 7, S. Di Domizio 17,11, M. L. Di Vacri 3,20, A. Drobizhev 5,6, D. Q. Fang 16, M. Faverzani 12,13, G. Fernanes 17,11, E. Ferri 13, F. Ferroni 9,10, E. Fiorini 13,12, M. A. Franceschi 21, S. J. Freeman 6,5,a, B. K. Fujikawa 6, A. Giachero 13, L. Gironi 12,13, A. Giuliani 22, L. Glastone 15, P. Gorla 3, C. Gotti 12,13, T. D. Gutierrez 23, E. E. Haller 8,24, K. Han 25, E. Hansen 15,2, K. M. Heeger 18, R. Hennings-Yeomans 5,6, K. P. Hickerson 2, H. Z. Huang 2, R. Kael 26, G. Keppel 4, Yu. G. Kolomensky 5,6, A. Leer 15, C. Ligi 21, K. E. Lim 18, Y. G. Ma 16, M. Maino 12,13, L. Marini 17,11, M. Martinez 9,10,27, R. H. Maruyama 18, Y. Mei 6, N. Moggi 28,7, S. Morganti 10, P. J. Mosteiro 10, T. Napolitano 21, M. Nastasi 12,13, C. Nones 29, E. B. Norman 30,31, V. Novati 22, A. Nucciotti 12,13, T. O Donnell 32, J. L. Ouellet 15, C. E. Pagliarone 3,33, M. Pallavicini 17,11, V. Palmieri 4, L. Pattavina 3, M. Pavan 12,13, G. Pessina 13, V. Pettinacci 10, G. Piperno 9,10,b, C. Pira 4, S. Pirro 3, S. Pozzi 12,13, E. Previtali 13, C. Rosenfel 1, C. Rusconi 1,3, M. Sakai 2, S. Sangiorgio 30, D. Santone 3,20, B. Schmit 6, J. Schmit 2, N. D. Scielzo 30, V. Singh 5, M. Sisti 12,13, A. R. Smith 6, L. Taffarello 14, M. Tenconi 22, F. Terranova 12,13, C. Tomei 10, S. Trentalange 2, M. Vignati 10, S. L. Wagaarachchi 5,6, B. S. Wang 30,31, H. W. Wang 16, B. Welliver 6, J. Wilson 1, L. A. Winslow 15, T. Wise 18,34, A. Woocraft 35, L. Zanotti 12,13, G. Q. Zhang 16, B. X. Zhu 2, S. Zimmermann 36, S. Zucchelli 37,7 1 Department of Physics an Astronomy, University of South Carolina, Columbia, SC 29208, USA 2 Department of Physics an Astronomy, University of California, Los Angeles, CA 90095, USA 3 INFN Laboratori Nazionali el Gran Sasso, Assergi L Aquila I-67010, Italy 4 INFN Laboratori Nazionali Legnaro, Legnaro Paova I-35020, Italy 5 Department of Physics, University of California, Berkeley, CA 94720, USA 6 Nuclear Science Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA 7 INFN Sezione Bologna, Bologna I-40127, Italy 8 Materials Science Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA 9 Dipartimento Fisica, Sapienza Università Roma, Roma I-00185, Italy 10 INFN Sezione Roma, Roma I-00185, Italy 11 INFN Sezione Genova, Genova I-16146, Italy 12 Dipartimento Fisica, Università Milano-Bicocca, Milano I-20126, Italy 13 INFN Sezione Milano Bicocca, Milano I-20126, Italy 14 INFN Sezione Paova, Paova I-35131, Italy 15 Massachusetts Institute of Technology, Cambrige, MA 02139, USA 16 Shanghai Institute of Applie Physics, Chinese Acaemy of Sciences, Shanghai , China 17 Dipartimento Fisica, Università Genova, Genova I-16146, Italy 18 Department of Physics, Yale University, New Haven, CT 06520, USA 19 INFN Gran Sasso Science Institute, L Aquila I-67100, Italy 20 Dipartimento Scienze Fisiche e Chimiche, Università ell Aquila, L Aquila I-67100, Italy 21 INFN Laboratori Nazionali Frascati, Frascati Roma I-00044, Italy 22 CSNSM, Univ. Paris-Su, CNRS/IN2P3, Universit Paris-Saclay, Orsay, France 23 Physics Department, California Polytechnic State University, San Luis Obispo, CA 93407, USA 24 Department of Materials Science an Engineering, University of California, Berkeley, CA 94720, USA 25 Department of Physics an Astronomy, Shanghai Jiao Tong University, Shanghai , China 26 Physics Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA 27 Laboratorio e Fisica Nuclear y Astroparticulas, Universia e Zaragoza, Zaragoza 50009, Spain 28 Dipartimento Scienze per la Qualità ella Vita, Alma Mater Stuorum Università Bologna, Bologna I-47921, Italy 29 Service e Physique es Particules, CEA / Saclay, Gif-sur-Yvette, France 30 Lawrence Livermore National Laboratory, Livermore, CA 94550, USA 31 Department of Nuclear Engineering, University of California, Berkeley, CA 94720, USA 32 Center for Neutrino Physics, Virginia Polytechnic Institute an State University, Blacksburg, Virginia 24061, USA 33 Dipartimento Ingegneria Civile e Meccanica, Università egli Stu Cassino e el Lazio Merionale, Cassino I-03043, Italy 34 Department of Physics, University of Wisconsin, Mason, WI 53706, USA 35 SUPA, Institute for Astronomy, University of Enburgh, Blackfor Hill, Enburgh EH9 3HJ, UK 36 Engineering Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA 37 Dipartimento Fisica e Astronomia, Alma Mater Stuorum Università Bologna, Bologna I-40127, Italy Receive: ate / Accepte: ate

2 2 Abstract We report a stuy of the CUORE sensitivity to neutrinoless ouble beta ββ ecay. We use a Bayesian analysis base on a toy Monte Carlo MC approach to extract the exclusion sensitivity to the ββ ecay half-life T at 90% crebility interval CI i.e. the interval containing the true value of T with 90% probability an the 3 σ scovery sensitivity. We consier various backgroun levels an energy resolutions, an escribe the influence of the ata vision in subsets with fferent backgroun levels. If the backgroun level an the energy resolution meet the expectation, CUORE will reach a 90% CI exclusion sensitivity of yr with 3 months, an yr with 5 years of live time. Uner the same contions, the scovery sensitivity after 3 months an 5 years will be yr an yr, respectively. 1 Introuction Neutrinoless ouble beta ecay is a non Stanar Moel process that violates the total lepton number conservation an implies a Majorana neutrino mass component [1, 2]. This ecay is currently being investigate with a variety of ouble beta ecaying isotopes. A recent review can be foun in Ref. [3]. The Cryogenic Unergroun Observatory for Rare Events CUORE [4,5, 6] is an experiment searching for ββ ecay in 130 Te. It is locate at the Laboratori Nazionali el Gran Sasso of INFN, Italy. In CUORE, 988 TeO 2 crystals with natural 130 Te isotopic abunance an a 750 g average mass are operate simultaneously as source an bolometric etector for the ecay. In this way, the ββ ecay signature is a peak at the Q-value of the reaction Q ββ, kev for 130 Te [7,8,9]. Bolometric crystals are characterize by an excellent energy resolution 0.2% Full With at Half Maximum, FWHM an a very low backgroun at Q ββ, which is expecte to be at the 10-2 cts/kev kg yr level in CUORE [10]. The current best limit on ββ ecay in 130 Te comes from a combine analysis of the CUORE-0 [11,12] an Cuoricino ata [13,14]. With a total exposure of 29.6 kg yr, a limit of T > yr 90% CI is obtaine [15] for the ββ ecay half life, T. After the installation of the etector, successfully complete in the summer 2016, CUORE starte the commissioning phase at the beginning of The knowlege of the scovery an exclusion sensitivity to ββ ecay as a function of the measurement live time can be exploite to set the criteria for the unblinng a Decease b Presently at: INFN Laboratori Nazionali Frascati, Frascati Roma I-00044, Italy of the ata an the release of the ββ ecay analysis results. In this work, we ecate our attention to those factors which coul strongly affect the sensitivity, such as the Backgroun Inex BI an the energy resolution at Q ββ. In CUORE, the crystals in the outer part of the array are expecte to show a higher BI than those in the mile [10]. Consiering this an following the strategy alreay implemente by the Gera Collaboration [16,17], we show how the vision of the ata into subsets with fferent BI coul improve the sensitivity. The reporte results are obtaine by means of a Bayesian analysis performe with the Bayesian Analysis Toolkit BAT [18]. The analysis is base on a toy- MC approach. At a cost of a much longer computation time with respect to the use of the mean sensitivity formula [19], this provies the full sensitivity probability stribution an not only its mean value. In Section 2, we review the statistical methos for the parameter estimation, as well as for the extraction of the exclusion an scovery sensitivity. Section 3 escribes the experimental parameters use for the analysis while its technical implementation is summarize in Section 4. Finally, we present the results in Section 5. 2 Statistical Metho The computation of exclusion an scovery sensitivities presente here follows a Bayesian approach: we exploit the Bayes theorem both for parameter estimation an moel comparison. In this work, we use the following notation: H incates both a hypothesis an the corresponng moel; H 0 is the backgroun-only hypothesis, accorng to which the known physics processes are enough to explain the experimental ata. In the present case, we expect the CUORE backgroun to be flat in a 100 kev region aroun Q ββ, except for the presence of a 60 Co summation peak at kev. Therefore, H 0 is implemente as a flat backgroun stribution plus a Gaussian escribing the 60 Co peak. In CUORE-0, this peak was foun to be centere at an energy 1.9 ± 0.7 kev higher than that tabulate in literature [15]. This effect, present also in Cuoricino [14], is a feature of all gamma summation peaks. Hence, we will consier the 60 Co peak to be at kev. H 1 is the backgroun-plus-signal hypothesis, for which some new physics is require to explain the ata. In our case, the physics involve in H 1 contains the

3 3 backgroun processes as well as ββ ecay. The latter is moele as a Gaussian peak at Q ββ. E represents the ata. It is a list of N energy bins centere at the energy E i an containing n i event counts. The energy range is [2470; 2570] kev. This is the same range use for the CUORE-0 ββ ecay analysis [15], an is boune by the possible presence of peaks from 214 Bi at kev an 208 Tl X- ray escape at 2585 kev [15]. While an unbinne fit allows to fully exploit the information containe in the ata, it can result in a long computation time for large ata samples. Given an energy resolution of 5 kev FWHM an using a 1 kev bin with, the ±3 sigma range of a Gaussian peak is containe in 12.7 bins. With the 1 kev binning choice, the loss of information with respect to the unbinne fit is negligible. Γ is the parameter escribing the ββ ecay rate for H 1 : Γ = ln 2 T. 1 θ is the list of nuisance parameters escribing the backgroun processes in both H 0 an H 1 ; Ω is the parameter space for the parameters θ. 2.1 Parameter Estimation We perform the parameter estimation for a moel H through the Bayes theorem, which yiels the probability stribution for the parameters base on the measure ata, uner the assumption that the moel H is correct. In the ββ ecay analysis, we are intereste in the measurement of Γ for the hypothesis H 1. The probability stribution for the parameter set Γ, θ is: P Γ, θ E, H1 = Ω 0 P E Γ, θ, H 1 πγ πθ P E. 2 Γ, θ, H 1 πγ πθ θ Γ The numerator contains the contional probability P E Γ, θ, H 1 of finng the measure ata E given the moel H 1 for a set of parameters Γ, θ, times the prior probability π for each of the consiere parameters. The prior probability has to be chosen accorng to the knowlege available before the analysis of the current ata. For instance, the prior for the number of signal counts Γ might be base on the half-life limits reporte by previous experiments while the prior for the backgroun level in the region of interest ROI coul be set base on the extrapolation of the backgroun measure outsie the ROI. The enominator represents the overall probability to obtain the ata E given the hypothesis H 1 an all possible parameter combinations, P E H 1. The posterior probability stribution for Γ is obtaine via marginalization, i.e. integrating P Γ, θ E, H 1 over all nuisance parameters θ: P Γ H 1, E = Ω P Γ, θ E, H 1 θ. 3 For each moel H, the probability of the ata given the moel an the parameters has to be efine. For a fixe set of experimental ata, this correspons to the likelihoo function [20]. Diving the ata into N subsets with inex characterize by fferent backgroun levels, an consiering a binne energy spectrum with N bins an a number n of events in the bin i of the subset spectrum, the likelihoo function is expresse by the prouct of a Poisson term for each bin : P E Γ, θ, H = L E Γ, θ, H = N =1 N i=1 e λ λ n n!, 4 where λ is the expectation value for the bin. The best-fit is efine as the set of parameter values Γ, θ for which the likelihoo is at its global maximum. In the practical case, we perform the maximization on the log-likelihoo ln L E Γ, θ, H = N =1 N i=1 λ + ln λ n, 5 where the ative terms ln n! are roppe from the calculation. The fference between H 0 an H 1 is manifeste in the formulation of λ. As mentione above, we parametrize H 0 with a flat stribution over the consiere energy range, i.e. [2470; 2570] kev: f bkg E = 1 E max E min 6 plus a Gaussian stribution for the 60 Co peak: 1 f Co E = exp E µ Co 2 2π σ 2σ 2. 7 The expecte backgroun counts in the bin correspons to the integral of f bkg E in the bin times the total number of backgroun counts M bkg for the subset : λ bkg = E max E min M bkg f bkg EE 8

4 4 where E min an E max are the left an right margins of the energy bin, respectively. Consiering bins of size δe an expressing M bkg as function of the backgroun inex BI, of the total mass m an of the measurement live time t, we obtain: λ bkg = BI m t δe. 9 Similarly, the expectation value for the 60 Co stribution on the bin is: λ Co = E max E min M Co exp E µ Co 2 2π σ 2σ 2 E, 10 where M Co is the total number of 60 Co events for the subset an can be reefine as function of the 60 Co event rate, R Co: M Co = R Co m t. 11 The total expectation value λ for H 0 is then: λ = λ bkg + λ Co. 12 In the case of H 1 an ational expectation value for ββ ecay is require: λ = E max E min M exp E Q ββ 2 2π σ 2σ 2 E. 13 The number of ββ ecay events in the subset is: M = Γ N A m A f 130 ε tot m t, 14 where N A is the Avogaro number, m a an f 130 are the molar mass an the isotopic abunance of 130 Te an ε tot is the total efficiency, i.e. the prouct of the containment efficiency ε MC obtaine with MC simulations an the instrumental efficiency ε instr. 2.2 Exclusion Sensitivity We compute the exclusion sensitivity by means of the 90% CI limit. This is efine as the value of T corresponng to the 90% quantile of the posterior Γ stribution: T 90% CI = T : ln 2/T 0 P Γ H1, E Γ = An example of posterior probability for Γ an the relative 90% CI limit is shown in Fig. 1, top. Flat prior stributions are use for all parameters, as escribe in Sec. 3. In the Bayesian approach, the limit is a statement regarng the true value of the consiere physical quantity. In our case, a 90% CI limit on T is to be interprete as the value above which, given the current knowlege, the true value of T lies with 90% probability. This ffers from a frequentist 90% C.L. limit, which is a statement regarng the possible results of the repetition of ientical measurements an shoul be interprete as the value above which the best-fit value of T woul lie in the 90% of the imaginary ientical experiments. In orer to extract the exclusion sensitivity, we generate a set of N toy-mc spectra accorng to the backgroun-only moel, H 0. We then fit spectra with the backgroun-plus-signal moel, H 1, an obtain the T an 90% CI stribution Fig. 1, bottom. Its me- ˆT 90% CI is referre as the mean sensitivity. For a real experiment, the experimental T limit P Γ H1, E Number of toy-mc exp ln 2 T 90% c.i ˆT 3000 Γ [yr -1 ] 1 σ 2 σ 3 σ T 90% c.i. [yr] T Fig. 1: Top: marginalize probability stribution for Γ relative to one toy-mc experiment. The 90% CI limit on is incate. Bottom: stribution of the 90% CI limits on T for 105 toy-mc experiments. The re vertical line correspons to the mean sensitivity, while the three colors epict the 1, 2 an 3 σ quantiles of the stribution. Both plots are obtaine with 1 yr live time, a 10-2 cts/kev kg yr BI an 5 kev energy resolution.

5 5 is expecte to be above/below ˆT 90% CI with 50% probability. Alternatively, one can consier the moe of the stribution, which correspons to the most probable T limit. The exact proceure for the computation of the exclusion sensitivity is the following: for each subset, we generate a ranom number of backgroun events N bkg accorng to a Poisson stribution with mean λ bkg ; for each subset, we generate N bkg events with an energy ranomly stribute accorng to f bkg E; we repeate the proceure for the 60 Co contribution; we fit the toy-mc spectrum with the H 1 moel Eq. 2, an marginalize the likelihoo with respect to the parameters BI an R Co Eq. 3; we extract the 90% CI limit on T ; we repeat the algorithm for N toy-mc experiments, an buil the stribution of T 90% CI. 2.3 Discovery Sensitivity The scovery sensitivity provies information on the require strength of the signal amplitue for claiming that the known processes alone are not sufficient to properly escribe the experimental ata. It is compute on the basis of the comparison between the backgrounonly an the backgroun-plus-signal moels. A metho for the calculation of the Bayesian scovery sensitivity was introuce in Ref. [21]. We report it here for completeness. In our case, we assume that H 0 an H 1 are a complete set of moels, for which: P H 0 E + P H 1 E = The application of the Bayes theorem to the moels H 0 an H 1 yiels: P H 0 E = P E H0 πh0 P E P H 1 E = P E H1 πh1 P E. 17 In this case, the numerator contains the probability of measuring the ata E given the moel H: P E H0 = P E θ, H0 πθ θ Ω P E H1 = P E Γ 18, θ, H 1 Ω 0 πγ πθ θ Γ, while the prior probabilities for the moels H 0 an H 1 can be chosen as 0.5 so that neither moel is favore. The enominator of Eq. 17 is the sum probability of obtaining the ata E given either the moel H 0 or H 1 : P E = P E H0 πh0 + P E H1 πh1. 19 At this point we nee to efine a criterion for claiming the scovery of new physics. Our choice is to quote the 3 σ mean scovery sensitivity, i.e. the value of for which the posterior probability of the backgroun-only moel H 0 given the ata is smaller than in 50% of the possible experiments. In other wors: T ˆT 3σ = T : P H 0 E < for N/2 experiments. 20 The etaile proceure for the etermination of the scovery sensitivity is: we prouce a toy-mc spectrum accorng to the H 1 moel with an arbitrary value of T ; we fit the spectrum with both H 0 an H 1 ; we compute P H 0 E; we repeat the proceure for N toy-mc spectra using the same T ; we repeat the routine with fferent values of T until the contion of Eq. 20 is satisfie. The iteration is implemente using the bisection metho until a precision is obtaine on the mean P H 0 E. 3 Experimental Parameters The fit parameters of the H 1 moel are BI, R Co an Γ, while only the first two are present for H 0. If the ata are vie in subsets, fferent BI an R Co fit parameter are consiere for each subset. On the contrary, the inverse ββ half-life is common to all subsets. Prior to the assembly of the CUORE crystal towers, we performe a screening survey of the employe materials [22,23,24,25,26,27,28,29]. From these measurements, either a non-zero activity was obtaine, or a 90% confience level C.L. upper limit was set. Ationally, the raoactive contamination of the crystals an holers was also obtaine from the CUORE-0 backgroun moel [30]. We evelope a full MC simulation of CUORE [10], an we use the results of the screening measurements an of the CUORE-0 backgroun moel for the normalization of the simulate spectra. We then compute the BI at Q ββ using the output of the simulations. In the present stuy, we consier only those

6 6 Table 1: Input parameters for the prouction of toy-mc spectra. BI [cts/kev kg yr] R Co [cts/kg yr] 1.02 ± 0.03stat syst backgroun contributions for which a non-zero activity is obtaine from the available measurements. The largest backgroun consists of α particles emitte by U an Th surface contaminations of the copper structure holng the crystals. Ationally, we consier a 60 Co contribution normalize to the 90% C.L. limit from the screening measurement. In this sense, the effect of a 60 Co backgroun on the CUORE sensitivity is to be hel as an upper limit. Given the 60 Co importance especially in case of sub-optimal energy resolution, we preferre to maintain a conservative approach in this regar. In the generation of the toy-mc spectra, we take into account the 60 Co half life 5.27 yr, an set the start of ata taking to January The parameter values use for the prouction of the toy-mc are reporte in Tab. 1. The quote uncertainty on the BI comes from the CUORE MC simulations [10]. We prouce the toy-mc spectra using the best-fit value of the BI. In a secon time, we repeat the analysis after increasing an ecreasing the BI by an amount equivalent to its statistical an systematic uncertainties combine in quarature. After running the fit on the entire crystal array as if it were a unique etector, we consiere the possibility of ving the ata grouping the crystals with a similar BI. Namely, being the backgroun at Q ββ ominate by surface α contamination of the copper structure, the crystals facing a larger copper surface are expecte to have a larger BI. This effect was alreay observe in CUORE-0, where the crystals in the uppermost an lowermost levels, which ha 3 sies facing the copper shiel, were characterize by a larger backgroun than those in all other levels, which were expose to coppper only on 2 sies. Consiering the CUORE geometry, the crystals can be vie in 4 subsets with fferent numbers of expose faces. Corresponngly, they are characterize by fferent BI, as reporte in Tab. 2. A major ingreent of a Bayesian analysis is the choice of the priors. In the present case, we use a flat prior for all parameters. In particular, the prior stribution for Γ is flat between zero an a value large enough to contain > % of its posterior stribution. This correspons to the most conservative choice. Any other reasonable prior, e.g. a scale invariant prior on Γ, woul yiel a stronger limit. A fferent prior choice base on the real characteristic of the experimental spectra might be more appropriate for BI an R Co in the analysis of the CUORE ata. For the time being the lack of ata prevents the use of informative priors. As a cross-check, we performe the analysis using the BI an 60 Co rate uncertainties obtaine by the backgroun buget as the σ of a Gaussian prior. No significant fference was foun on the sensitivity ban because the Poisson fluctuations of the generate number of backgroun an 60 Co events are ominant for the extraction of the Γ posterior probability stribution. Tab. 3 lists the constant quantities present in the formulation of H 0 an H 1. All of them are fixe, with the exception of the live time t an the FWHM of the ββ ecay an 60 Co Gaussian peaks. We perform the analysis with a FWHM of 5 an 10 kev, corresponng to a σ of 2.12 an 4.25 kev, respectively. Regarng the efficiency, while in the toy-mc prouction the BI an R Co are multiplie by the instrumental efficiency 1, in the fit the total efficiency is use. This is the prouct of the containment an instrumental efficiency. Also in this case, we use the same value as for CUORE-0, i.e. 81.3% [15]. We evaluate the exclusion an scovery sensitivities for fferent live times, with t ranging from 0.1 to 5 yr an using logarithmically increasing values: t i = 1.05 t i 1. Table 2: Crystal subsets with fferent expecte α backgroun in CUORE. The values of BI an R Co are taken from [10]. Subset Free Number of BI [cts/ R Co [cts Name Sies crystals kev kg yr] /kg yr] Inner Mile Mile Outer Table 3: Constants use in H 0 an H 1. Constant Symbol Value Detector Mass m kg Avogaro number N A mol 1 Molar mass m A g/mol Live Time t yr Efficiency ε tot 81.3% 130 Te abunance f ββ Q-value Q ββ kev 60 Co peak position µ Co kev Energy resolution FWHM 5, 10 kev 1 The containment efficiency is alreay encompasse in BI an R Co [10].

7 7 4 Fit Proceure We perform the analysis with the software BAT v DEV [21], which internally uses CUBA [31] v4.2 for the integration of multi-mensional probabilities an the Metropolis-Hastings algorithm [32] for the fit. The computation time epens on the number of samples rawn from the consiere probability stribution. For the exclusion sensitivity, we raw 10 5 likelihoo samples for every toy-mc experiment, while, ue to the higher computational cost, we use only 10 3 for the scovery sensitivity. For every combination of live time, BI an energy resolution, we run toy-mc experiments for the exclusion scovery sensitivity stuy. In the case of the scovery sensitivity, we chose the number of toy-mc experiments as the minimum for which a 2% relative precision was achievable on the mean sensitivity. For the exclusion sensitivity, it was possible to increase both the number of toy-mc experiments an iterations, with a systematic uncertainty on the mean sensitivity at the per mil level. 5 Results an Discussion 5.1 Exclusion Sensitivity The stributions of 90% CI limit as a function of live time with no ata subvision are shown in Fig. 2. For all BI values an all live times, the FWHM of 5 kev yiels a 45% higher sensitivity with respect to a 10 kev resolution. The mean sensitivity after 3 months an 5 years of ata collection in the two consiere cases are reporte in Tab. 4. The epenence of the mean sensitivity on live time is typical of a backgroun-ominate experiments: namely, CUORE expects about one event every four ays in a ±3σ region aroun Q ββ. The results in Tab. 4 show also the importance of energy resolution an suggest to put a strong effort in its optimization. As a cross check, we compare the sensitivity just obtaine with that provie by the analytical metho presente in [19] an shown in ark green in Fig. 2. The analytical metho yiels a slightly higher sentitivity for short live times, while the two techniques agree when the ata sample is bigger. We justify this with the fact that the uncertainty on the number of backgroun counts obtaine with the Bayesian fit is slightly larger than the corresponng Poisson uncertainty assume in the analytical approach [33], hence the limit on T is systematically weaker 2. The effect becomes less an less strong with in- 2 See the scussion of the pulls for a more etaile explanation. creasing ata samples, i.e. with growing live time. With a resolution of 5 kev, the fference goes from 8% after 3 months to < 0.1% after 5 years, while for a 10 kev FWHM the fference is 6% after 3 months an 4% after 5 years. One remark has to be one concerning the values reporte in [19]: there we gave a 90% C.I. exclusion sensitivity of yr with 5 yr of live time. This is 5% higher than the result presente here an is explaine by the use of a fferent total efficiency: 87.4% in [19] an 81.3% in this work. We then extract the exclusion sensitivity after ving the crystals in 4 subsets, as escribe in Sec. 3. The mean exclusion sensitivity values after 3 months an 5 years of ata collection with one an 4 subsets are reporte in Tab. 4. The vision in subsets yiels only a small improvement at the percent level in mean sensitivity. Base on this results only, one woul conclue that ving the ata into subsets with fferent BI is not worth the effort. This conclusion is not always true, an strongly relies on the exposure an BI of the consiere subsets. As an example, we repeate a toy analysis assuming a BI of 10-2 cts/kev kg yr, an with two subsets of equal exposure an BI cts/kev kg yr an cts/kev kg yr, respectively. In this case, the vision of the ata in to two subsets yiels a 10% improvement after 5 yr of ata taking. Hence, the ata subvision is a viable option for the final analysis, whose gain strongly epens on the experimental BI of each channel. Similarly, we expect the CUORE bolometers to have fferent energy resolutions; in CUORE-0, these range from 3 kev to 20 kev FWHM [34]. In the real CUORE analysis a further splitting of the ata can be one by grouping the channels with similar FWHM, or by keeping every channels separate. At the present stage it is not possible to make reliable prections for the FWHM stribution among the crystals, so we assume an average value of 5 or 10 kev throughout the whole work. Ieally, the final CUORE ββ ecay analysis shoul be performe keeping the spectra collecte by each crystal separate, ationally to the usual vision of the ata into ata sets comprise by two calibration runs [15]. Assuming an average frequency of one calibration per month, the total number of energy spectra woul be Assuming a fferent but stationary BI for each crystal, an using the same 60 Co rate for all crystals, the fit moel woul have 10 3 parameters. This represents a major obstacle for any existing implementation of the Metropolis-Hastings or Gibbs sampling algorithm. A possible way to aress the problem might be the use of fferent algorithms, e.g. neste sampling [35, 36], or a partial analytical solution of the likelihoo maximization.

8 8 Exposure [kg yr] BI: 1.02 ± cts/kev kg yr FWHM: 5 kev Exposure [kg yr] BI: 1.02 ± cts/kev kg yr FWHM: 10 kev 90% c.i. limit [yr] % c.i. limit [yr] T CUORE-0 + CUORICINO Mean 1 σ 2 σ 3 σ T CUORE-0 + CUORICINO Mean 1 σ 2 σ 3 σ Live time [yr] Live time [yr] Fig. 2: 90% CI exclusion sensitivity for a 5 kev left an a 10 kev FWHM right. The re crosses correspon to the mean sensitivity for a BI of 10-2 cts/kev kg yr, while the smaller black crosses correspon to the mean sensitivity obtaine after shifting the BI up an own by an amount equivalent to its statistical an systematic uncertainty summe in quarature. The uncertainty on BI yiels a +4% 9% effect on the mean sensitivity. The fferent colore areas epict the ranges containing the 68.3, 95.5 an 99.7% of the toy-mc experiments. They are compute for each live time value separately as escribe in Fig. 1, bottom. We also show the sensitivity compute as in [19] in ark green. The horizontal green line at yr correspons to the limit obtaine with CUORE-0 an Cuoricino [15]. Table 4: Mean exclusion sensitivity for fferent energy resolutions an fferent subset numbers. FWHM N ˆT at 0.25 yr ˆT at 5 yr [kev] [yr] [yr] We perform two further cross-checks in orer to investigate the relative importance of the flat backgroun an the 60 Co peak. In the first scenario we set the BI to zero, an o the same for the 60 Co rate in the secon one. In both cases, the ata are not vie into subsets, an resolutions of 5 an 10 kev are consiere. With no flat backgroun an a 5 kev resolution, no 60 Co event leaks in the ±3σ region aroun Q ββ even after 5 yr of measurement. As a consequence, the 90% CI limits are stribute on a very narrow ban, an the mean sensitivity reaches yr after 5 yr of ata collection. On the contrary, if we assume a 10 kev FWHM, some 60 Co events fall in the ββ ecay ROI from the very beginning of the ata taking. This results in a strong asymmetry of the sensitivity ban. In the secon crosscheck, we keep the BI at cts/kev kg yr, but an λ only Eqs. 9 an 13. This is also confirme by the computation of the sensitivity for the optimistic scenario without the 1.9 kev shift of the 60 Co peak use in the stanar case. We test the fit correctness an bias computing the pulls, i.e. the normalize resiuals, of the number of counts assigne to each of the fit components. Denoting with N bkg an N Co the number of generate backgroun an 60 Co events, respectively, an with M bkg an M Co the corresponng number of reconstructe events, the pulls are efine as: set the 60 Co rate to zero. In both cases, the fference with respect to the stanar scenario is below 1%. We can conclue that the 60 Co peak with an initial rate of cts/kg yr is not worrisome for a resolution of up to 10 kev, an that the lower sensitivity obtaine with 10 kev FWHM with respect to the 5 kev case is ascribable to the relative amplitue of λ bkg r bkgco = M bkgco N bkgco σ M bkgco, 21 where σ M bkgco is the statistical uncertainty on M bkgco given by the fit. For an unbiase fit, the stribution of the pulls is expecte to be Gaussian with a unitary root mean square RMS. In the case of exclusion sensitivity, we

9 9 3σ [yr] ˆT Exposure [kg yr] BI: 1.02 ± cts/kev kg yr CUORE-0 + CUORICINO 5 kev putation cost, we compute the 3 σ scovery sensitivity for a FWHM of 5 an 10 kev with no crystal subvision. As shown in Fig. 3, with a 5 kev energy resolution CUORE has a 3 σ scovery sensitivity superior to the limit obtaine from the combine analysis of Cuore-0 an Cuoricino ata [15] after less than one month of operation, an reaches yr with 5 yr of live time. Also in this case, the pulls are characterize by an RMS smaller than expecte, but no bias is present ue to the use of H 1 for both the generation an the fit of the toy-mc spectra. 10 kev Live time [yr] Fig. 3: 3σ scovery sensitivity with a BI of cts/kev kg yr an an FWHM of 5 kev. obtain r bkg = 0.2 ± 0.4 an r Co = 0.1 ± 0.5 for all live times. The fact that the pull stributions are slightly shifte incates the presence of a bias. Its origin lies in the Bayesian nature of the fit an is that all fit contributions are constraine to be greater than zero. We perform a cross-check, by extenng the range of all parameters BI, R Co an Γ to negative values. Uner this contion, the bias sappears. In ation to this, an explanation is neee for the small RMS of the pull stributions. This is mainly ue to two effects: first, the toy-mc spectra are generate using H 0, while the fit is performe using H 1 ; secon, the statistical uncertainties on all parameters are larger than the Poisson uncertainty on the number of generate events. To confirm the first statement, we repeat the fit using H 0 instea of H 1 an we obtain pulls with zero mean an an RMS 0.8, which is closer to the expecte value. Finally, we compare the parameter uncertainty obtaine from the fit with the Poisson uncertainty for the equivalent number of counts, an we fin that the fference is of O20%. 5.2 Discovery Sensitivity The extraction of the scovery sensitivity involves fits with the backgroun-only an the backgroun-plus-signal moels. Moreover, two multi-mensional integrations have to be performe for each toy-mc spectrum, an a loop over the ββ ecay half-life has to be one until the contion of Eq. 20 is met. Due to the high com- 6 Conclusion an Outlook We implemente a toy-mc metho for the computation of the exclusion an scovery sensitivity of CUORE using a Bayesian analysis. We have highlighte the influence of the BI an energy resolution on the exclusion sensitivity, showing how the achievement of the expecte 5 kev FWHM is esirable. Ationally, we have shown how the vision of the ata into subsets with fferent BI coul yiel an improvement in exclusion sensitivity. Once the CUORE ata collection starts an the experimental parameters are available, the sensitivity stuy can be repeate in a more etaile way. As an example, non-gaussian spectral shapes for the ββ ecay an 60 Co peaks can be use, an the systematics of the energy reconstruction can be inclue. Acknowlegments The CUORE Collaboration thanks the rectors an staff of the Laboratori Nazionali el Gran Sasso an the technical staff of our laboratories. CUORE is supporte by The Istituto Nazionale Fisica Nucleare INFN; The National Science Founation uner Grant Nos. NSF-PHY , NSF-PHY , NSF- PHY , NSF-PHY , NSF-PHY , NSF-PHY , NSF-PHY , NSF-PHY , an NSF-PHY ; The Alfre P. Sloan Founation; The University of Wisconsin Founation; Yale University; The US Department of Energy DOE Office of Science uner Contract Nos. DE-AC02-05CH1-1231, DE-AC52-07NA27344, an DE-SC ; The DOE Office of Science, Office of Nuclear Physics uner Contract Nos. DE-FG02-08ER41551 an DE-FG03-00ER41138; The National Energy Research Scientific Computing Center NERSC.

10 10 References 1. J. Schechter an J. W. F. Valle, Phys. Rev. D25, M. Duerr, M. Linner, an A. Merle, JHEP 06, S. Dell Oro et al., Av. High Energy Phys. 2016, D. R. Artusa et al., [CUORE Collaboration], Av. High Energy Phys. 2015, C. Arnabol et al., [CUORE Collaboration], Nucl. Instrum. Meth. A518, C. Arnabol et al., [CUORE Collaboration], Astropart. Phys. 20, M. Reshaw et al., Phys. Rev. Lett. 102, N. D. Scielzo et al., Phys. Rev. C80, S. Rahaman et al., Phys. Lett. B703, C. Aluino et al., [CUORE Collaboration], Eur. Phys. J. C 2017, x.oi.org/ /epjc/ s C. Aluino et al., [CUORE Collaboration], JINST 11, P D. R. Artusa et al., [CUORE Collaboration], Eur. Phys. J. C74, C. Arnabol et al., [Cuoricino Collaboration], Phys. Rev. C78, E. Anreotti et al., [Cuoricino Collaboration], Astropart. Phys. 34, K. Alfonso et al., [CUORE Collaboration], Phys. Rev. Lett. 115, M. Agostini et al., GERDA, Phys. Rev. Lett. 111, M. Agostini et al., Nature 544, A. Calwell et al., Comput. Phys. Commun. 180, F. Alessanria et al., [CUORE Collaboration], 2011, arxiv: v F. James, Statistical Methos in Experimental Physics, 2 e. Worl Scientific, A. Calwell an K. Kroninger, Phys. Rev. D74, F. Alessanria et al., [CUORE Collaboration], Astropart. Phys. 35, A. F. Barghouty et al., Nucl. Instrum. Meth. B295, B. S. Wang et al., Phys. Rev. C92, F. Alessanria et al., [CUORE Collaboration], Astropart. Phys. 45, E. Anreotti et al., [Cuoricino Collaboration], Astropart. Phys. 34, F. Bellini et al., Astropart. Phys. 33, E. Anreotti et al., JINST 4, P A. Giachero, Characterization of cryogenic bolometers an ata acquisition system for the CUORE experiment, PhD thesis, Università egli Stu Genova, C. Aluino et al., [CUORE Collaboration], Eur. Phys. J. C77, T. Hahn, Comput. Phys. Commun. 168, D. D. L. Minh an D. L. P. Minh, Communications in Statistics - Simulation an Computation 44, G. Cowan, K. Cranmer, E. Gross, an O. Vitells, 2011, arxiv: C. Aluino et al., [CUORE Collaboration], Phys. Rev. C93, F. Feroz et al., Mon. Not. Roy. Astron. Soc. 398, W. J. Hanley et al., Mon. Not. Roy. Astron. Soc. 450, L