PHYSICAL REVIEW D 100, 015028 (2019)
Cosmic neutrino background search experiments as decaying dark matter detectors
David McKeen* TRIUMF, 4004 Wesbrook Mall, Vancouver, British Columbia V6T 2A3, Canada
(Received 1 April 2019; published 18 July 2019)
We investigate the possibility that particles that are long-lived on cosmological scales, making up part or all of the dark matter, decay to neutrinos that have present-day energies around an eV. The neutrinos from these decays can potentially be visible at experiments that hope to directly observe the cosmic neutrino background through neutrino capture on tritium, such as PTOLEMY. In the context of a simple model that can realize such decays, we discuss the allowed signatures at a PTOLEMY-like experiment given current cosmological constraints.
DOI: 10.1103/PhysRevD.100.015028
I. INTRODUCTION enormous density, scattering of CνB neutrinos is highly suppressed since their energies are minuscule and they Particles with lifetimes on cosmological scales that couple to matter only via the weak interaction. Further- decay to neutrinos arise in a number of new physics more, the tiny energies involved make distinguishing from contexts, such as lepton number violation and the gener- backgrounds very troublesome. Both of these facts mean ation of neutrino masses [1], and can even comprise some that detecting the CνB is extremely challenging. The most or all of the dark matter of the Universe [2–5]. The present- promising technique to detect the very low energy CνB— day energy of the neutrinos from these decays can span a − neutrino capture on β-decaying nuclei, ν þðA; ZÞ → e þ wide range, making them accessible at, e.g., existing dark ðA; Z þ 1Þ, where a nucleus, ðA; ZÞ, is labeled by its mass matter direct detection and neutrino experiments [4–6]. number, A, and atomic number, Z—was first proposed by In this paper we will consider the decays of long-lived Weinberg in 1962 [9] and more recently studied in [10,11]. particles to neutrinos that carry very little energy today, in This process benefits from the lack of a threshold energy the neighborhood of an eV. These neutrinos are slightly and its signature is the production of an electron with more energetic than the standard cosmic neutrino back- energy above the end point of natural β decay, ðA; ZÞ → ground (CνB)—the thermal relic neutrinos that decoupled − ν¯ þ e þðA; Z þ 1Þ. For a neutrino of mass mν and energy from the plasma when the Universe had a temperature of a Eν, the shift above the end point is ∼mν þ Eν. In the case of few MeV. The temperature of the CνBis∼10−4 eV today, nonrelativistic CνB neutrinos, this shift is roughly 2mν.Of less than the scale of the atmospheric and solar neutrino course, extremely good resolution on the e− energy is mass splittings, and therefore at least two mass eigenstates required to resolve this gap given the relatively large rate ν in the C B are now nonrelativistic; the most energetic of of β decay compared to capture of CνB neutrinos. ≃ ≳ 0 05 these neutrinos have energy ECνB mν . eV. There is Recently, the PTOLEMY experiment has proposed [12] an upper limit on the neutrino masses, hence on the energy to tackle these difficulties using a target of 100 g of tritium ν of the C B, from observations of the cosmic microwave (i.e., A ¼ 3, Z ¼ 1 in the expressions above) implanted on Pbackground (CMB) and baryon acoustic oscillations of a graphene substrate along with MAC-E (Magnetic 0 12 mν< . eV at 95% confidence level [7]. Terrestrial Adiabatic Collimation with Electrostatic) filtering, radio experiments, in comparison, limit the neutrino mass at frequency monitoring, and advanced calorimetry to mea- 95% confidence level to less than 2.8 eV [8]. sure the e− energy. Tritium has a relatively small rate for β The present-day number density of CνB neutrinos is very −1 decay, Γβ ≃ ð17 yrÞ and implanting it on graphene serves ≃ 330 −3 large, with a cosmic average nCνB cm . Despite this to reduce intrinsic broadening of the e− energy from molecular effects. This could potentially allow for energy *[email protected] resolution as small as 0.1 eV, which is needed to success- fully probe the CνB. Published by the American Physical Society under the terms of This paper is organized as follows. In Sec. II, we describe the Creative Commons Attribution 4.0 International license. Further distribution of this work must maintain attribution to a simple model containing a long-lived particle that the author(s) and the published article’s title, journal citation, decays to neutrinos, and discuss the number density and and DOI. Funded by SCOAP3. energy distribution of neutrinos allowed by cosmological
2470-0010=2019=100(1)=015028(8) 015028-1 Published by the American Physical Society DAVID MCKEEN PHYS. REV. D 100, 015028 (2019) observations. We see that decays that take place after matter, and h ≃ 0.68 is the Hubble constant today in units 100 Ω Ω recombination, when photons decouple, are the most prom- of km=s=Mpc. J= dm is the fraction of dark matter ising to be probed at experiments sensitive to eV-scale that J particles would comprise today if they did not decay, neutrinos, such as PTOLEMY. We discuss the signature of which we normalize here to 5%. In what follows, we will this scenario at a PTOLEMY-like experiment in detail in suppress the argument of the scale factor and set its value Sec. III, considering both the diffuse and local contributions today to unity, a0 ¼ aðt0Þ¼1. We will also require that to the flux, and compare the reach to that of cosmological the J’s and their decay products do not grossly disturb observations. In Sec. IV, we conclude. the evolution of the Universe from the standard cosmology, with radiation domination at redshifts z ¼ a−1 − 1≳ ≃ 3300 II. DENSITIES AND DISTRIBUTIONS zeq , followed by matter domination, then more recently vacuum energy domination. The simple model we consider involves only a long- The number density of neutrinos produced in J decays lived, real scalar particle J with mass mJ. It couples to the simply follows from Eq. (3), light neutrinos ν through an interaction of the form − τ 2Ω ρ 1 − t= J g ˜ ð Þ¼ J cr;0 e L ¼ − νν þ ð Þ nν t 3 int J H:c: 1 m a 2 J − τ 130 Ω Ω 1 − t= J J= dm eV e In this expression, ν is a two-component, left-chiral spinor ¼ ; ð4Þ cm3 0.05 m a3 field, and the interaction is gauge invariant if, e.g., it comes J from a sterile admixture of the light neutrinos. This → νν where we use a tilde to distinguish this population from the coupling leads to the decay J . Assuming that this standard neutrinos. is the only J decay mode at tree level, its lifetime is There are essentially three qualitatively different regimes 32π 10−15 2 for the J lifetime in terms of its cosmological effects: 9 eV τ ≃ 0 2 τ ¼ ¼ 2 × 10 yr : ð2Þ (i) before neutrino decoupling, J 015028-2 COSMIC NEUTRINO BACKGROUND SEARCH EXPERIMENTS AS … PHYS. REV. D 100, 015028 (2019) dρ˜ν mJnJ (∼0.28=15 ≃ 2% or less) goes into relativistic d.o.f. in þ 4Hρ˜ν ¼ : ð6Þ dt τJ each J decay. This could be the case, e.g., in neutrino portal dark matter models [13] with a small splitting in the dark −4 Using the fact that photons redshift like radiation, ργ ∝ a , sector between a fermion, χ and a scalar ϕ. These interact we can reexpress this equation as with neutrinos via the effective operator ϕχ¯HL=Λ → ðv=ΛÞϕχν¯ where H and L are the Higgs and lepton ðρ˜ ρ Þ d ν= γ mJnJ ΩJ a − τ doublets respectively. The heavier state in the dark sector, ¼ ¼ t= J ð Þ e ; 7 ϕ for instance, could be long-lived and decay through this dt ργτJ Ωγ τJ operator, ϕ → χν where the neutrino has energy ∼mϕ − mχ −5 −2 in the ϕ rest frame. where Ωγ ¼ 2.47 × 10 h is the present-day contribution of the CMB to the critical density. Assuming that the J’sdo Lastly, looking at Eq. (9), we might naively think that we J not come to dominate the energy budget of the Universe at can arrange for a number density of neutrinos from decays before recombination that is much larger than that of this time, the energy density is dominated by radiation and 2 the CνB neutrinos. This could be the case if the J’s were the scale factor depends on the time as t ¼ t2a with ¼ 7 6 1011 very light and cold so that their energy density is sup- t2 . × yr. We can then integrate Eq. (7) to find pressed while their number density is large—this relies on pffiffiffi rffiffiffiffiffi ’ ν ρ˜ π Ω τ the J s remaining unthermalized with the C B neutrinos. ν ¼ J J This requirement can be used to set an upper limit on the ρ 2 Ω γ t≫τ γ t2 strength of the J − ν interaction or, equivalently, a lower J rffiffiffiffiffiffiffiffiffiffiffiffiffi Ω τ limit on τJ. Production of J’s through νν → J happens most ¼ 0.15 J J : ð8Þ ∼ Ω 103 readily at Tν mJ and the rate for this is roughly dm yr 2 g Tν ∼ 1=τJ. To keep the J’s out of equilibrium, we require that this is less than the Hubble rate at Tν ∼ m which is Using Eqs. (4), (5), and (8), we can reexpress n˜ ν today in J ∼ 2 ≃1019 terms of ΔN , mJ=MPl with MPl GeV the Planck mass. This eff 2 ≲ τ ≳ 105 ð Þ2 implies that g mJ=MPl or J yr eV=mJ . sffiffiffiffiffiffiffiffiffiffiffiffiffi 1 0 103 Δ 103 Using this in Eq. (9) one obtains an upper bound on the . × Neff eV yr n˜ νðt0Þ¼ : ð9Þ present-day density of neutrinos from J decays roughly 3 0 28 τ −3 cm . mJ J comparable to that in the CνB, Oð100 cm Þ. The energy distribution of such neutrinos today would be indistin- An additional constraint comes from the fact that J guishable from that of the CνB; i.e., they would also be particles redshift like matter before their decay so that their nonrelativistic. energy density can come to exceed that in radiation, Similar situations where nonthermally produced neutri- causing an early period of matter domination. This would nos, such as the right-chiral component of light Dirac conflict with the usual picture of radiation domination from neutrinos, evade constraints on light d.o.f. and lead to an primordial nucleosynthesis until matter-radiation equality enhancement of the CνB signal have already been explored ∼ Ω τ ≲ ∼ at teq trec; we can use this to limit J for J teq trec. in, e.g., [14]. Dark matter that comes into thermal equi- The ratio of the energy density in J’s to that in radiation librium with neutrinos after neutrino decoupling but before during this era is recombination has been studied extensively in [15]. ρJ 2 mJnJ 2 ΩJ − τ τ t ¼ ¼ t= J B. Decays after recombination ( J > rec) ρ ρ Ω ae r g γ g γ ’ rffiffiffiffi Nonrelativistic J s that decay to neutrinos after recombi- 2 ΩJ t − τ nation act as a decaying component of the dark matter. This ¼ e t= J ; ð10Þ g Ωγ t2 alters the expansion history of the Universe between last scattering and today which can change the precise pattern with g ¼ 3.36. Requiring that this ratio is less than unity of CMB angular anisotropies as well as the growth of gives the constraint structure. These effects have been analyzed in detail in Refs. [3,16,17] and even proposed as an explanation of sffiffiffiffiffiffiffiffiffiffiffiffiffi tensions in cosmological data [18]. This limits τ and the Ω 103 J J 23 yr ð Þ energy density in J’s and, consequently, the number density Ω < τ : 11 dm J of neutrinos produced in J decay. For lifetimes short compared to the age of the Universe Note that this is equivalent to a rather weak limit on the (but long compared to trec), Ref. [16] obtains a bound using Ω Ω number of relativistic d.o.f. from CMB observations CMB observables on the J energy density of J= dm < Δ 15 of Neff < . However, it can become an important 0.038 at 95% confidence level. For longer lifetimes, the Ω Ω 0 09ðτ 15 Þ constraint in a scenario where only a fraction of mJ constraint is roughly J= dm < . J= Gyr . Using 015028-3 DAVID MCKEEN PHYS. REV. D 100, 015028 (2019) these constraints in Eq. (4), the current neutrino number these expressions the characteristic energy of the neutrinos density is then limited to be today is ( 95 pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi eV av 1 10−5 τ 103 τ ≲ n˜ νðt0Þ ≲ Eν × J= yr; J trec 3 m ¼ ð15Þ cm J m ½0 04ðτ Þ2=3 0 5 τ ≳ 8 J min . J=Gyr ; . ; J trec; > 1 τ ≲ 12 <> ; J Gyr 2 3ð1 − 2τ Þ 12 ≲ τ ≲ 160 which ignores a slight correction due to the relatively × > . t0= J ; Gyr J Gyr :> recent transition to vacuum energy domination. Therefore, 2.3ð160 Gyr=τ Þ; τ ≳ 160 Gyr: J J obtaining relativistic neutrinos today with Eν ∼ eV requires ≳ 10 ≲ 10 τ ≲ τ ≳ ð12Þ mJ keV (mJ keV) for J trec ( J trec), as we would expect given the redshift of matter-radiation equal- ≃ 3300 We observe that neutrino number densities comparable ity, zeq . to that in the CνB from J decay are possible for τJ larger or In Fig. 1, we show the cosmic average flux of non- smaller than trec without being ruled out by cosmological standard neutrinos as functions of their energy today for ¼ 50 τ ¼ 103 Ω data. We now turn to the question of their energy distri- mJ keV, J yr