The PHENIX experiment at the Relativistic Heavy Ion Collider has measured the differential cross section of $\phi$(1020)-meson production at forward rapidity in $p$$+$$p$ collisions at $\sqrt{s}=$510 GeV via the dimuon decay channel. The partial cross section in the rapidity and $p_T$ ranges $1.2<|y|<2.2$ and $2<p_T<7$ GeV/$c$ is $\sigma_\phi=[2.28 \pm 0.09\,{\rm (stat)} \pm 0.14\,{\rm (syst)} \pm 0.27\, {\rm (norm)}] \times 10^{-2}$~mb. The energy dependence of $\sigma_\phi$ ($1.2<|y|<2.2, \; 2<p_T<5$ GeV/$c$) is studied using the PHENIX measurements at $\sqrt{s}=$200 and 510 GeV and the Large-Hadron-Collider measurements at $\sqrt{s}=$2.76 and 7 TeV. The experimental results are compared to various event generator predictions ({\sc pythia6, pythia8, phojet, ampt, epos3,} and {\sc epos-lhc}).

We present a detailed measurement of charged two-pion correlation functions in 0%-30% centrality $\sqrt{s_{_{NN}}}=200$ GeV Au$+$Au collisions by the PHENIX experiment at the Relativistic Heavy Ion Collider. The data are well described by Bose-Einstein correlation functions stemming from L\'evy-stable source distributions. Using a fine transverse momentum binning, we extract the correlation strength parameter $\lambda$, the L\'evy index of stability $\alpha$ and the L\'evy length scale parameter $R$ as a function of average transverse mass of the pair $m_T$. We find that the positively and the negatively charged pion pairs yield consistent results, and their correlation functions are represented, within uncertainties, by the same L\'evy-stable source functions. The $\lambda(m_T)$ measurements indicate a decrease of the strength of the correlations at low $m_T$. The L\'evy length scale parameter $R(m_T)$ decreases with increasing $m_T$, following a hydrodynamically predicted type of scaling behavior. The values of the L\'evy index of stability $\alpha$ are found to be significantly lower than the Gaussian case of $\alpha=2$, but also significantly larger than the conjectured value that may characterize the critical point of a second-order quark-hadron phase transition.

Production of $\pi^0$ and $\eta$ mesons has been measured at midrapidity in Cu$+$Au collisions at $\sqrt{s_{_{NN}}}$=200 GeV. Measurements were performed in $\pi^0(\eta)\rightarrow\gamma\gamma$ decay channel in the 1(2)-20 GeV/$c$ transverse momentum range. A strong suppression is observed for $\pi^0$ and $\eta$ meson production at high transverse momentum in central Cu$+$Au collisions relative to the $p$$+$$p$ results scaled by the number of nucleon-nucleon collisions. In central collisions the suppression is similar to Au$+$Au with comparable nuclear overlap. The $\eta/\pi^0$ ratio measured as a function of transverse momentum is consistent with $m_T$-scaling parameterization down to $p_T=$2 GeV/$c$, its asymptotic value is constant and consistent with Au$+$Au and $p$$+$$p$ and does not show any significant dependence on collision centrality. Similar results were obtained in hadron-hadron, hadron-nucleus, and nucleus-nucleus collisions as well as in $e^+e^-$ collisions in a range of collision energies $\sqrt{s_{_{NN}}}=$3--1800 GeV. This suggests that the quark-gluon-plasma medium produced in Cu$+$Cu collisions either does not affect the jet fragmentation into light mesons or it affects the $\pi^0$ and $\eta$ the same way.

We have measured direct photons for $p_T<5~$GeV/$c$ in minimum bias and 0\%--40\% most central events at midrapidity for Cu$+$Cu collisions at $\sqrt{s_{_{NN}}}=200$ GeV. The $e^{+}e^{-}$ contribution from quasi-real direct virtual photons has been determined as an excess over the known hadronic contributions in the $e^{+}e^{-}$ mass distribution. A clear enhancement of photons over the binary scaled $p$$+$$p$ fit is observed for $p_T<4$ GeV/$c$ in Cu$+$Cu data. The $p_T$ spectra are consistent with the Au$+$Au data covering a similar number of participants. The inverse slopes of the exponential fits to the excess after subtraction of the $p$$+$$p$ baseline are 285$\pm$53(stat)$\pm$57(syst)~MeV/$c$ and 333$\pm$72(stat)$\pm$45(syst)~MeV/$c$ for minimum bias and 0\%--40\% most central events, respectively. The rapidity density, $dN/dy$, of photons demonstrates the same power law as a function of $dN_{\rm ch}/d\eta$ observed in Au$+$Au at the same collision energy.

PHENIX has measured the azimuthal correlations of muon pairs from charm and bottom semi-leptonic decays in $p$+$p$ collisions at $\sqrt{s}=200$ GeV, using a novel analysis technique utilizing both unlike- and like-sign muon pairs to separate charm, bottom and Drell-Yan contributions. The dimuon measurements combined with the previous electron-muon and dielectron measurements span a wide range in rapidity, and are well described by PYTHIA Tune A. Through a Bayesian analysis based on PYTHIA Tune A, we show that leading order pair creation is the dominant $(76\%\pm^{14}_{19}\%)$ contribution for $b\bar{b}$ production, whereas the data favor the scenario in which next-to-leading-order processes dominate $c\bar{c}$ production. The small contribution of next-to-leading-order processes in $b\bar{b}$ production at the collision energies of the Relativistic Heavy Ion Collider contrasts with the case at Large-Hadron-Collider energies, where next-to-leading-order processes are expected to dominate.

The PHENIX collaboration presents first measurements of low-momentum ($0.4<p_T<3$ GeV/$c$) direct-photon yields from Au$+$Au collisions at $\sqrt{s_{_{NN}}}$=39 and 62.4 GeV. For both beam energies the direct-photon yields are substantially enhanced with respect to expectations from prompt processes, similar to the yields observed in Au$+$Au collisions at $\sqrt{s_{_{NN}}}$=200. Analyzing the photon yield as a function of the experimental observable $dN_{\rm ch}/d\eta$ reveals that the low-momentum ($>$1\,GeV/$c$) direct-photon yield $dN_{\gamma}^{\rm dir}/d\eta$ is a smooth function of $dN_{\rm ch}/d\eta$ and can be well described as proportional to $(dN_{\rm ch}/d\eta)^\alpha$ with $\alpha{\sim}$1.25. This new scaling behavior holds for a wide range of beam energies at the Relativistic Heavy Ion Collider and Large Hadron Collider, for centrality selected samples, as well as for different, $A$$+$$A$ collision systems. At a given beam energy the scaling also holds for high $p_T$ ($>5$\,GeV/$c$) but when results from different collision energies are compared, an additional $\sqrt{s_{_{NN}}}$-dependent multiplicative factor is needed to describe the integrated-direct-photon yield.

The experimental study of the collisions of heavy nuclei at relativistic energies has established the properties of the quark-gluon plasma (QGP), a state of hot, dense nuclear matter in which quarks and gluons are not bound into hadrons. In this state, matter behaves as a nearly inviscid fluid that efficiently translates initial spatial anisotropies into correlated momentum anisotropies among the produced particles, producing a common velocity field pattern known as collective flow. In recent years, comparable momentum anisotropies have been measured in small-system proton-proton ($p$$+$$p$) and proton-nucleus ($p$$+$$A$) collisions, despite expectations that the volume and lifetime of the medium produced would be too small to form a QGP. Here, we report on the observation of elliptic and triangular flow patterns of charged particles produced in proton-gold ($p$$+$Au), deuteron-gold ($d$$+$Au), and helium-gold ($^3$He$+$Au) collisions at a nucleon-nucleon center-of-mass energy $\sqrt{s_{_{NN}}}$~=~200 GeV. The unique combination of three distinct initial geometries and two flow patterns provides unprecedented model discrimination. Hydrodynamical models, which include the formation of a short-lived QGP droplet, provide a simultaneous description of these measurements.

PHENIX reports differential cross sections of $\mu\mu$ pairs from semileptonic heavy-flavor decays and the Drell-Yan production mechanism measured in $p$$+$$p$ collisions at $\sqrt{s}=200$ GeV at forward and backward rapidity ($1.2<|\eta|<2.2$). The $\mu\mu$ pairs from $c\bar{c}$, $b\bar{b}$, and Drell-Yan are separated using a simultaneous fit to unlike- and like-sign muon pair spectra in mass and $p_T$. The azimuthal opening angle correlation between the muons from $c\bar{c}$ and $b\bar{b}$ decays and the pair-$p_T$ distributions are compared to distributions generated using {\sc pythia} and {\sc powheg} models, which both include next-to-leading order processes. The measured distributions for pairs from $c\bar{c}$ are consistent with {\sc pythia} calculations. The $c\bar{c}$ data presents narrower azimuthal correlations and softer $p_T$ distributions compared to distributions generated from {\sc powheg}. The $b\bar{b}$ data are well described by both models. The extrapolated total cross section for bottom production is $3.75{\pm}0.24({\rm stat}){\pm}^{0.35}_{0.50}({\rm syst}){\pm}0.45({\rm global})$[$\mu$b], which is consistent with previous measurements at the Relativistic Heavy Ion Collider in the same system at the same collision energy, and is approximately a factor of two higher than the central value calculated with theoretical models. The measured Drell-Yan cross section is in good agreement with next-to-leading-order quantum-chromodynamics calculations.

Dihadron and isolated direct photon-hadron angular correlations are measured in $p$$+$$p$ collisions at $\sqrt{s}=200$ GeV. The correlations are sensitive to nonperturbative initial-state and final-state transverse momentum $k_T$ and $j_T$ in the azimuthal nearly back-to-back region $\Delta\phi\sim\pi$. In this region, transverse-momentum-dependent evolution can be studied when several different hard scales are measured. To have sensitivity to small transverse momentum scales, nonperturbative momentum widths of $p_{\rm out}$, the out-of-plane transverse momentum component perpendicular to the trigger particle, are measured. These widths are used to investigate possible effects from transverse-momentum-dependent factorization breaking. When accounting for the longitudinal momentum fraction of the away-side hadron with respect to the near-side trigger particle, the widths are found to increase with the hard scale; this is qualitatively similar to the observed behavior in Drell-Yan and semi-inclusive deep-inelastic scattering interactions. The momentum widths are also studied as a function of center-of-mass energy by comparing to previous measurements at $\sqrt{s}=510$ GeV. The nonperturbative jet widths also appear to increase with $\sqrt{s}$ at a similar $x_T$, which is qualitatively consistent to similar measurements in Drell-Yan interactions. To quantify the magnitude of any transverse-momentum-dependent factorization breaking effects, calculations will need to be performed to compare to these measurements.

We report the transverse single-spin asymmetries of $J/\psi$ production at forward and backward rapidity, $1.2<|y|<2.2$, as a function of $J/\psi$ transverse momentum ($p_T$) and Feynman-$x$ ($x_F$). The data analyzed were recorded by the PHENIX experiment at the Relativistic Heavy Ion Collider in 2015 from $p$$+$$p$, $p$$+$Al, and $p$$+$Au collisions with transversely polarized proton beams at $\sqrt{s_{_{NN}}}=200$ GeV. At this collision energy, single-spin asymmetries for heavy-flavor particle production of $p$$+$$p$ collisions provide access to the spin-dependent gluon distribution and higher-twist correlation functions inside the nucleon, such as the gluon Qiu-Sterman and trigluon correlation functions. Proton+nucleus collisions offer an excellent opportunity to study nuclear effects on the correlation functions. The data indicate negative asymmetries at the two-standard-deviation level in the $p$$+$Au data for $p_T<2$ GeV/$c$ at both forward and backward rapidity, while in $p$$+$$p$ and $p$$+$Al collisions the asymmetries are consistent with zero within the range of experimental uncertainties.

We present measurements of elliptic and triangular azimuthal anisotropy of charged particles detected at forward rapidity $1<|\eta|<3$ in Au$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV, as a function of centrality. The multiparticle cumulant technique is used to obtain the elliptic flow coefficients $v_2\{2\}$, $v_2\{4\}$, $v_2\{6\}$, and $v_2\{8\}$, and triangular flow coefficients $v_3\{2\}$ and $v_3\{4\}$. Using the small-variance limit, we estimate the mean and variance of the event-by-event $v_2$ distribution from $v_2\{2\}$ and $v_2\{4\}$. In a complementary analysis, we also use a folding procedure to study the distributions of $v_2$ and $v_3$ directly, extracting both the mean and variance. Implications for initial geometrical fluctuations and their translation into the final state momentum distributions are discussed.

We have measured the cross section and single spin asymmetries from forward $W^{\pm}\rightarrow\mu^{\pm}\nu$ production in longitudinally polarized $p$$+$$p$ collisions at $\sqrt{s}=510$ GeV using the PHENIX detector at the Relativistic Heavy Ion Collider. The cross sections are consistent with previous measurements at this collision energy, while the most forward and backward longitudinal single spin asymmetries provide new insights into the sea quark helicities in the proton. The charge of the W bosons provides a natural flavor separation of the participating partons.

The CUORE experiment, a ton-scale cryogenic bolometer array, recently began operation at the Laboratori Nazionali del Gran Sasso in Italy. The array represents a significant advancement in this technology, and in this work we apply it for the first time to a high-sensitivity search for a lepton-number--violating process: $^{130}$Te neutrinoless double-beta decay. Examining a total TeO$_2$ exposure of 86.3 kg$\cdot$yr, characterized by an effective energy resolution of (7.7 $\pm$ 0.5) keV FWHM and a background in the region of interest of (0.014 $\pm$ 0.002) counts/(keV$\cdot$kg$\cdot$yr), we find no evidence for neutrinoless double-beta decay. The median statistical sensitivity of this search is $7.0\times10^{24}$ yr. Including systematic uncertainties, we place a lower limit on the decay half-life of $T^{0\nu}_{1/2}$($^{130}$Te) > $1.3\times 10^{25}$ yr (90% C.L.). Combining this result with those of two earlier experiments, Cuoricino and CUORE-0, we find $T^{0\nu}_{1/2}$($^{130}$Te) > $1.5\times 10^{25}$ yr (90% C.L.), which is the most stringent limit to date on this decay. Interpreting this result as a limit on the effective Majorana neutrino mass, we find $m_{\beta\beta}<(110 - 520)$ meV, where the range reflects the nuclear matrix element estimates employed.

We report in situ neutron background measurements at the Kuo-Sheng Reactor Neutrino Laboratory (KSNL) by a hybrid neutron detector (HND) with a data size of 33.8 days under identical shielding configurations as during the neutrino physics data taking. The HND consists of BC-501A liquid and BC-702 phosphor powder scintillation neutron detectors, which is sensitive to both fast and thermal neutrons, respectively. Neutron-induced events for the two channels are identified and differentiated by pulse shape analysis, such that background of both are simultaneously measured. The fast neutron fluxes are derived by an iterative unfolding algorithm. Neutron induced background in the germanium detector under the same fluxes, both due to cosmic-rays and ambient radioactivity, are derived and compared with the measurements. The results are valuable to background understanding of the neutrino data at the KSNL. In particular, neutron-induced background events due to ambient radioactivity as well as from reactor operation are negligible compared to intrinsic cosmogenic activity and ambient $\gamma$-activity. The detector concept and analysis procedures are applicable to neutron background characterization in similar rare-event experiments.

We present measurements of azimuthal correlations of charged hadron pairs in $\sqrt{s_{_{NN}}}=200$ GeV Au$+$Au collisions after subtracting an underlying event using a model that includes higher-order azimuthal anisotropy $v_2$, $v_3$, and $v_4$. After subtraction, the away-side ($\Delta\phi\sim\pi)$ of the highest transverse-momentum trigger ($p_T>4$ GeV/$c$) correlations is suppressed compared to that of correlations measured in $p$$+$$p$ collisions. At the lowest associated particle $p_T$, the away-side shape and yield are modified. These observations are consistent with the scenario of radiative-jet energy loss. For the lowest-$p_T$ trigger correlations, an away-side yield exists and we explore the dependence of the shape of the away-side within the context of an underlying-event model. Correlations are also studied differentially versus event-plane angle $\Psi_n$. The angular correlations show an asymmetry when selecting the sign of the trigger-particle azimuthal angle with respect to the $\Psi_2$ event plane. This asymmetry and the measured suppression of the pair yield out of plane is consistent with a path-length-dependent energy loss. No $\Psi_3$ dependence can be resolved within experimental uncertainties.

Understanding of particle physics will be incomplete without knowing the true nature i.e. Dirac or Majorana, and absolute effective mass of neutrinos. This will lead us into the world of beyond the standard model physics. In the current era of experimental technology it is probably only neutrinoless double beta decay (0$\nu\beta\beta$), which has the ability to clarify the nature and provide the absolute effective mass of neutrinos. Observation of 0$\nu\beta\beta$ implies not only neutrinos are Majorana particles but also the lepton number violation. There are already intense activities going on over the world-wide for the experimental searches of 0$\nu\beta\beta$. In these experiments, some are currently in their running phase and several experiments have been projected. With the same motivation as the other experiments, the TIN-TIN (The INdia$^{'}$s TIN detector) experiment has been projected by the scientific community of India. This experiment, in coming future, will be housed in the approved underground laboratory of INO (India based Neutrino Observatory), to search 0$\nu\beta\beta$ using enriched $^{124}$Sn isotope as the target mass. In the present paper we have calculated the required sensitivity for $^{124}$Sn isotope in terms of exposure ($\beta\beta{isotope}$ mass $\times$ time) and background rate (counts/ton-year-keV) at the projected FWHM energy resolution. Such that, after achieving the required sensitivity the TIN-TIN can really see the events of 0$\nu\beta\beta$ decay, if it exists in nature.

CUORE - the Cryogenic Underground Observatory for Rare Events - is an experiment searching for the neutrinoless double-beta ($0\nu\beta\beta$) decay of $^{130}$Te with an array of 988 TeO$_2$ crystals operated as bolometers at $\sim$10 mK in a large dilution refrigerator. With this detector, we aim for a $^{130}$Te $0\nu\beta\beta$ decay half-life sensitivity of $9\times10^{25}$ y with 5 y of live time, and a background index of $\lesssim 10^{-2}$ counts/keV/kg/y. Making an effort to maintain radiopurity by minimizing the bolometers' exposure to radon gas during their installation in the cryostat, we perform all operations inside a dedicated cleanroom environment with a controlled radon-reduced atmosphere. In this paper, we discuss the design and performance of the CUORE Radon Abatement System and cleanroom, as well as a system to monitor the radon level in real time.

TeO2 bolometers have been used for many years to search for neutrinoless double beta decay in 130-Te. CUORE, a tonne-scale TeO2 detector array, recently published the most sensitive limit on the half-life, $T_{1/2}^{0\nu} > 1.5 \times 10^{25}\,$yr, which corresponds to an upper bound of $140-400$~meV on the effective Majorana mass of the neutrino. While it makes CUORE a world-leading experiment looking for neutrinoless double beta decay, it is not the only study that CUORE will contribute to in the field of nuclear and particle physics. As already done over the years with many small-scale experiments, CUORE will investigate both rare decays (such as the two-neutrino double beta decay of 130-Te and the hypothesized electron capture in 123-Te), and rare processes (e.g., dark matter and axion interactions). This paper describes some of the achievements of past experiments that used TeO2 bolometers, and perspectives for CUORE.

Recently, multiparticle-correlation measurements of relativistic $p/d/^3$He$+$Au, $p$$+$Pb, and even $p$$+$$p$ collisions have shown surprising collective signatures. Here we present beam-energy-scan measurements of 2-, 4-, and 6-particle angular correlations in $d$$+$Au collisions at $\sqrt{s_{_{NN}}}$=200, 62.4, 39, and 19.6 GeV. We also present measurements of 2- and 4-particle angular correlations in $p$$+$Au collisions at $\sqrt{s_{_{NN}}}$=200 GeV. We find the 4-particle cumulant to be real-valued for $d$$+$Au collisions at all four energies. We also find that the 4-particle cumulant in $p$$+$Au has the opposite sign as that in $d$$+$Au, indicating that the correlations are geometrical and therefore collective in origin. Further we find that the 6-particle cumulant agrees with the 4-particle cumulant in $d$$+$Au collisions at 200 GeV, which indicates that nonflow effects are subdominant and provides strong evidence of collectivity.

During 2015 the Relativistic Heavy Ion Collider (RHIC) provided collisions of transversely polarized protons with Au and Al nuclei for the first time, enabling the exploration of transverse-single-spin asymmetries with heavy nuclei. Large single-spin asymmetries in very forward neutron production have been previously observed in transversely polarized $p$$+$$p$ collisions at RHIC, and the existing theoretical framework that was successful in describing the single-spin asymmetry in $p$$+$$p$ collisions predicts only a moderate atomic-mass-number ($A$) dependence. In contrast, the asymmetries observed at RHIC in $p$$+$$A$ collisions showed a surprisingly strong $A$ dependence in inclusive forward neutron production. The observed asymmetry in $p$$+$Al collisions is much smaller, while the asymmetry in $p$$+$Au collisions is a factor of three larger in absolute value and of opposite sign. The interplay of different neutron production mechanisms is discussed as a possible explanation of the observed $A$ dependence.

We have performed a search for neutrinoless $\beta^{+}\hspace{-0.2em}EC$ decay of $^{120}$Te using the final CUORE-0 data release. We describe a new analysis method for the simultaneous fit of signatures with different event topology, and of data subsets with different signal efficiency, obtaining a limit on the half-life of the decay of $T_{1/2}>1.6\cdot10^{21}$ yr at $90\%$ CI. Combining this with results from Cuoricino, a predecessor experiment, we obtain the strongest limit to date, corresponding to $T_{1/2}>2.7\cdot10^{21}$ yr at $90\%$ CI.

The observation of neutrinoless double-beta decay (0${\nu}{\beta}{\beta}$) would show that lepton number is violated, reveal that neutrinos are Majorana particles, and provide information on neutrino mass. A discovery-capable experiment covering the inverted ordering region, with effective Majorana neutrino masses of 15 - 50 meV, will require a tonne-scale experiment with excellent energy resolution and extremely low backgrounds, at the level of $\sim$0.1 count /(FWHM$\cdot$t$\cdot$yr) in the region of the signal. The current generation $^{76}$Ge experiments GERDA and the MAJORANA DEMONSTRATOR utilizing high purity Germanium detectors with an intrinsic energy resolution of 0.12%, have achieved the lowest backgrounds by over an order of magnitude in the 0${\nu}{\beta}{\beta}$ signal region of all 0${\nu}{\beta}{\beta}$ experiments. Building on this success, the LEGEND collaboration has been formed to pursue a tonne-scale $^{76}$Ge experiment. The collaboration aims to develop a phased 0${\nu}{\beta}{\beta}$ experimental program with discovery potential at a half-life approaching or at $10^{28}$ years, using existing resources as appropriate to expedite physics results.

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.

CUORE is a tonne-scale cryogenic detector operating at the Laboratori Nazionali del Gran Sasso (LNGS) that uses tellurium dioxide bolometers to search for neutrinoless double-beta decay of $^{130}$Te. CUORE is also suitable to search for low energy rare events such as solar axions or WIMP scattering, thanks to its ultra-low background and large target mass. However, to conduct such sensitive searches requires improving the energy threshold to 10 keV. In this paper, we describe the analysis techniques developed for the low energy analysis of CUORE-like detectors, using the data acquired from November 2013 to March 2015 by CUORE-0, a single-tower prototype designed to validate the assembly procedure and new cleaning techniques of CUORE. We explain the energy threshold optimization, continuous monitoring of the trigger efficiency, data and event selection, and energy calibration at low energies in detail. We also present the low energy background spectrum of CUORE-0 below 60keV. Finally, we report the sensitivity of CUORE to WIMP annual modulation using the CUORE-0 energy threshold and background, as well as an estimate of the uncertainty on the nuclear quenching factor from nuclear recoils in CUORE-0.

We report a study of the CUORE sensitivity to neutrinoless double beta ($0\nu\beta\beta$) decay. We used a Bayesian analysis based on a toy Monte Carlo (MC) approach to extract the exclusion sensitivity to the $0\nu\beta\beta$ decay half-life ($T_{1/2}^{0\nu}$) at $90\%$ credibility interval (CI) -- i.e. the interval containing the true value of $T_{1/2}^{0\nu}$ with $90\%$ probability -- and the $3 \sigma$ discovery sensitivity. We consider various background levels and energy resolutions, and describe the influence of the data division in subsets with different background levels. If the background level and the energy resolution meet the expectation, CUORE will reach a $90\%$ CI exclusion sensitivity of $2\cdot10^{25}$ yr with $3$ months, and $9\cdot10^{25}$ yr with $5$ years of live time. Under the same conditions, the discovery sensitivity after $3$ months and $5$ years will be $7\cdot10^{24}$ yr and $4\cdot10^{25}$ yr, respectively.