Cosmic neutrino background search experiments as decaying dark matter detectors

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.


I. INTRODUCTION
Particles with lifetimes on cosmological scales that decay to neutrinos arise in a number of new physics contexts, such as lepton number violation and the generation of neutrino masses [1], and can even comprise some or all of the dark matter of the Universe [2][3][4][5]. The present day energy of the neutrinos from these decays can span a wide range, making them accessible at, e.g., existing dark matter direct detection and neutrino experiments [4][5][6].
In this paper we will consider the decays of long-lived particles to neutrinos that carry very little energy today, in the neighborhood of an eV. These neutrinos are slightly more energetic than the standard cosmic neutrino background (CνB)-the thermal relic neutrinos that decoupled from the plasma when the Universe had a temperature of a few MeV. The temperature of the CνB is ∼ 10 −4 eV today, less than the scale of the atmospheric and solar neutrino mass splittings, and therefore at least two mass eigenstates in the CνB are now nonrelativistic; the most energetic of these neutrinos have energy E CνB m ν 0.05 eV. There is an upper limit on the neutrino masses, hence on the energy of the CνB, from observations of the cosmic microwave background (CMB) and baryon acoustic oscillations of m ν < 0.12 eV at 95% confidence level [7]. Terrestrial experiments, in comparison, limit the neutrino mass at 95% confidence level to less than 2.8 eV [8].
The present day number density of CνB neutrinos is very large, with a cosmic average n CνB 330 cm −3 . Despite this enormous density, scattering of CνB neutrinos is highly suppressed since their energies are minuscule and they couple to matter only via the weak interaction. Furthermore, the tiny energies involved make distinguishing from backgrounds very troublesome. Both of these facts mean that detecting the CνB is extremely challenging. The most promising technique to detect the very low energy CνB-neutrino capture on β-decaying nuclei, ν + (A, Z) → e − + (A, Z + 1), where a nucleus, (A, Z), is labeled by its mass number, A, and atomic number, Z-was first proposed by Weinberg in 1962 [9]. This process benefits from the lack of a thresh- * mckeen@triumf.ca old energy and its signature is the production of an electron with energy above the endpoint of natural β-decay, (A, Z) →ν + e − + (A, Z + 1). For a neutrino of mass m ν and energy E ν , the shift above the endpoint is ∼ m ν +E ν . In the case of nonrelativistic CνB neutrinos, this shift is roughly 2m ν . Of course, extremely good resolution on the e − energy is required to resolve this gap given the relatively large rate of β-decay compared to capture of CνB neutrinos.
Recently, the PTOLEMY experiment has proposed [10] to tackle these difficulties using a target of 100 g of tritium (i.e. A = 3, Z = 1 in the expressions above) implanted on a graphene substrate along with MAC-E filtering, radio frequency monitoring, and advanced calorimetry to measure the e − energy. Tritium has a relatively small rate for β-decay, Γ β (17 yr) −1 and implanting it on graphene serves to reduce intrinsic broadening of the e − energy from molecular effects. This could potentially allow for energy resolution as small as 0.1 eV, which is needed to successfully probe the CνB. This paper is organized as follows. In Sec. II, we describe a simple model containing a long-lived particle that decays to neutrinos, and discuss the number density and energy distribution of neutrinos allowed by cosmological observations. We see that decays that take place after recombination, when photons decouple, are the most promising to be probed at experiments sensitive to eVscale neutrinos, such as PTOLEMY. We discuss the signature of this scenario at a PTOLEMY-like experiment in detail in Sec. III, considering both the diffuse and local contributions to the flux, and compare the reach to that of cosmological observations. In Sec. IV, we conclude.

II. DENSITIES AND DISTRIBUTIONS
The simple model we consider involves only a longlived, real scalar particle J with mass m J . It couples to the light neutrinos ν through an interaction of the form In this expression, ν is a two-component, left-chiral spinor field, and the interaction is gauge invariant if, e.g., it comes from a sterile admixture of the light neutrinos.
This coupling leads to the decay J → νν. Assuming that this is the only J decay mode at tree level, its lifetime is Here and in what follows, we take the neutrino masses to be negligible compared to m J . A natural candidate for J is the majoron associated with the spontaneous breaking of lepton number at a scale f . In this case, the coupling is g = im ν /f which could easily be tiny for f above the TeV scale. We have normalized the J mass on an eV for later convenience. For now, we are agnostic about the flavor structure of the couplings in Eq. (1) but will return to this point in Sec. III. We assume that the J's are produced nonthermally and are nonrelativistic at cosmologically interesting times. 1 We can then simply write down their number density, In this expression, a(t) is the scale factor of the Universe, t 0 = 13.8 Gyr is its age, Ω J is the J energy density in units of the critical energy density of the Universe today, ρ cr,0 = 10.5 h 2 keV/cm 3 , Ω dm = 0.12 h −2 is that of dark matter, and h 0.68 is the Hubble constant today in units of 100 km/s/Mpc. Ω J /Ω dm is the fraction of dark matter that J particles would comprise today if they did not decay, which we normalize here to 5%. In what follows, we will suppress the argument of the scale factor and set its value today to unity, a 0 = a(t 0 ) = 1. We will also require that the J's and their decay products do not grossly disturb the evolution of the Universe from the standard cosmology, with radiation domination at redshifts z = a −1 −1 z eq 3300, followed by matter domination, then more recently vacuum energy domination.
The number density of neutrinos produced in J decays simply follows from Eq. (3), There are essentially three qualitatively different regimes for the J lifetime in terms of its cosmological effects: (i) before neutrino decoupling, τ J < t dec 0.2 s, (ii) after neutrino decoupling and before recombination, t dec < τ J < t rec ∼ 4 × 10 5 yr, and (iii) after recombination, τ J > t rec . Case (i) is not observable since the neutrinos simply thermalize with the plasma. Cases (ii) and (iii) can potentially lead to a nonstandard population of neutrinos today. Crucially, cases (ii) and (iii) affect the observation of the CMB in different ways which we discuss below.
A. Decays before recombination (τJ < trec) In this case, the neutrinos from J decays contribute to the energy density in relativistic species at early times. This is constrained by the observation of the CMB which is record of the Universe at around t rec .
The extra contribution from J decays can be conveniently parametrized by a shift of the effective number of relativistic degrees of freedom, ∆N eff , with where ρ γ andρ ν are the energy densities in photons and neutrinos from J decay, respectively. A nonzero ∆N eff can affect the CMB by changing the expansion rate around the time of last scattering from its standard value. The current 95% confidence level upper limit on ∆N eff from CMB and large scale structure observations is ∆N eff < 0.28 [7]. This upper limit can constrain the number of neutrinos from J decay for τ J < t rec .
To determine this limit, we relate the energy density in neutrinos from J decay to the J number density, Using the fact that photons redshift like radiation, ρ γ ∝ a −4 , we can re-express this equation as where Ω γ = 2.47 × 10 −5 h −2 is the present day contribution of the CMB to the critical density. Assuming that the J's do not come to dominate the energy budget of the Universe at this time, the energy density is dominated by radiation and the scale factor depends on the time as t = t 2 a 2 with t 2 = 7.6 × 10 11 yr. We can then integrate Eq. 7 to find Using Eqs. (4), (5), and (8), we can re-expressñ ν today in terms of ∆N eff , An additional constraint comes from the fact that J particles redshift like matter before their decay so that their energy density can come to exceed that in radiation, causing an early period of matter domination. This would conflict with the usual picture of radiation domination from primordial nucleosynthesis until matterradiation equality at t eq ∼ t rec ; we can use this to limit Ω J for τ J t eq ∼ t rec . The ratio of the energy density in J's to that in radiation during this era is with g * = 3.36. Requiring that this ratio is less than unity gives the constraint Note that this is equivalent to a rather weak limit on the number of relativistic degrees of freedom from CMB observations of ∆N eff < 15. However, it can become an important constraint in a scenario where only a fraction of m J (∼ 0.28/15 2% or less) goes into relativistic degrees of freedom in each J decay. This could be the case, e.g., in neutrino portal dark matter models [11] with a small splitting in the dark sector between a fermion, χ and a scalar φ. These interact with neutrinos via the effective operator φχHL/Λ → (v/Λ)φχν where H and L are the Higgs and lepton doublets respectively. The heavier state in the dark sector, φ for instance, could be long-lived and decay through this operator, φ → χν where the neutrino has energy ∼ m φ − m χ in the φ rest frame.
Lastly, looking at Eq. (9),we might naively think that we can arrange for a number density of neutrinos from J decays before recombination that is much larger than that of the CνB neutrinos. This could be the case if the J's were very light and cold so that their energy density is suppressed while their number density is largethis relies on the J's remaining unthermalized with the CνB neutrinos. This requirement can be used to set an upper limit on the strength of the J-ν interaction or, equivalently, a lower limit on τ J . Production of J's through νν → J happens most readily at T ν ∼ m J and the rate for this is roughly g 2 T ν ∼ 1/τ J . To keep the J's out of equilibrium, we require that this is less than the Hubble rate at T ν ∼ m J which is ∼ m 2 J /M Pl with M Pl 10 19 GeV the Planck mass. This implies that g 2 m J /M Pl or τ J 10 5 yr (eV/m J ) 2 . Using this in Eq. (9) one obtains an upper bound on the present day density of neutrinos from J decays roughly comparable to that in the CνB, O(100 cm −3 ). The energy distribution of such neutrinos today would be indistinguishable from that of the CνB, i.e. they would also be nonrelativistic. Similar situations where nonthermally produced neutrinos, such as the right-chiral component of light Dirac neutrinos, evade constraints on light degrees of freedom and lead to an enhancement of the CνB signal have already been explored in, e.g. [12]. Dark matter that comes into thermal equilibrium with neutrinos after neutrino decoupling but before recombination has been studied extensively in [13].

B. Decays after recombination (τJ > trec)
Nonrelativistic J's that decay to neutrinos after recombination act as a decaying component of the dark matter. This alters the expansion history of the Universe between last scattering and today which can change the precise pattern of CMB angular anisotropies as well as the growth of structure. These effects have been analyzed in detail in Refs. [3,14,15] and even proposed as an explanation of tensions in cosmological data [16]. This limits τ J and the energy density in J's and, consequently, the number density of neutrinos produced in J decay.
For lifetimes short compared to the age of the Universe (but long compared to t rec ), Ref. [14] obtains a bound using CMB observables on the J energy density of Ω J /Ω dm < 0.038 at 95% confidence level. For longer lifetimes, the constraint is roughly Ω J /Ω dm < 0.09 (τ J /15 Gyr). Using these constraints in Eq. (4), the current neutrino number density is then limited to bẽ 12 Gyr τ J 160 Gyr 2.3 (160 Gyr/τ J ) , τ J 160 Gyr.
We observe that neutrino number densities comparable to that in the CνB from J decay are possible for τ J larger or smaller than t rec without being ruled out by cosmological data. We now turn to the question of their energy distribution which greatly impacts their detectability.

C. Energy spectrum
In this section, we find the (cosmic average) energy density of the neutrinos produced in J decays, assuming m ν m J . At the time of their production, each neutrino has energy m J /2. Afterwards, the energy redshifts with the expansion of the Universe so that today it is E ν = (m J /2)a where a is the scale factor at the time of decay (a 0 = 1). Since the J's do not decay all at once, the neutrino spectrum today is not monochromatic, broadened by the expansion of the Universe over the course of the J decays. More quantitatively, the spectrum today is simply related to the change in the comoving number density of neutrinos with scale factor, where the time and Hubble rate are evaluated at a = 2E ν /m J . Assuming that during J decay the Hubble rate depends on the scale factor simply as H = (nt n a n ) −1 (n = 2, 3/2 correspond to radiation-and matter-domination respectively), we can find the "average" energy of the neutrinos today, at which dñ ν /dE ν is maximized, t n is a constant that depends on whether the Universe was radiation-or matter-dominated at the time of decay. Again, assuming the J's and their neutrino decay products do not change the evolution of the Universe from the standard picture, t 2 is given above Eq. (8) and t 3/2 = 1.7 × 10 10 yr which apply to τ J < t eq and τ J > t eq , respectively, where t eq ∼ t rec is the time of matter-radiation equality. Using these expressions the characteristic energy of the neutrinos today is which ignores a slight correction due to the relatively recent transition to vacuum energy domination. Therefore, obtaining relativistic neutrinos today with E ν ∼ eV requires m J 10 keV (m J 10 keV) for τ J t rec (τ J t rec ), as we would expect given the redshift of matter-radiation equality, z eq 3300.
In Fig. 1, we show the cosmic average flux of nonstandard neutrinos as functions of their energy today for m J = 50 keV, τ J = 10 3 yr < t rec and Ω J corresponding to ∆N eff = 0.28 (solid cyan); m J = 1 eV, τ J = 7 × 10 9 yr, Ω J /Ω dm = 0.08 (solid, green); and m J = 3 eV, τ J = 10 11 yr, Ω J = Ω dm (solid, red). We also show the upper limit on the flux from decays prior to recombination such that they contribute ∆N eff = 0.28 (dashed, cyan). To generate these fluxes we used the energy spectrum from Eq. (13) with the Hubble rate calculated for a Universe with contributions to the energy density from baryons, Ω b = 0.05; dark energy, Ω Λ = 0.69; dark matter, Ω dm = 0.26; and (standard) radiation, Ω r = 8.7 × 10 −5 . For τ J > t rec , the dark matter is split up into a decaying component, Ω J , and a nondecaying component, Ω dm − Ω J .
The flux of neutrinos from J → νν, averaged over the Universe, along with CνB, solar ν [17,19], and atmospheric ν [18] neutrinos, as well as antineutrinos from the decay of neutrons and tritons produced in primordial nucleosynthesis [20]. We show the flux for mJ = 50 keV, τJ = 10 3 yr, and, since the decays take place before recombination, ΩJ such that this gives ∆N eff = 0.28 (solid, cyan). In addition, we show the upper limit on any flux from a decay taking place before recombination from the upper limit ∆N eff = 0.28 (dashed, cyan). Also shown are the the fluxes for two benchmark points for decays after trec in Sec. III: mJ = 1 eV, τJ = 7 × 10 9 yr, ΩJ = 0.08 Ω dm (solid, green) and mJ = 3 eV, τJ = 10 11 yr, ΩJ = Ω dm (solid, red).
We also show the results of the recent computations of the low energy thermal component of solar neutrinos [19] and the antineutrino flux from the decay of neutrons and tritons produced in primordial nucleosynthesis [20], which extend into the 0.1-1 eV range but in relatively low numbers.
Because the number density is suppressed with increasing m J , and cosmological limits are strong for small τ J , the most viable scenario to produce an observable number of neutrinos with E ν ∼ eV today involves relatively light J's, of order a few eV, with lifetimes comparable to the age of the Universe. If the light neutrinos are massive, it is possible for early decays, before t rec , of light J's to lead to a population of nonrelativistic neutrinos with present-day energy E ν ∼ m ν . The signal of these neutrinos is indistinguishable from that of the CνB in such a case. Since, as seen in Sec. (II A), such neutrinos' number density can be at most comparable to the CνB, they can only lead to an O(1) increase in the event rate of massive CνB-like neutrinos, similar to scenarios explored in [12]. The standard CνB rate can vary by an O(1) factor depending on whether the neutrinos are Majorana or Dirac [21] as well as gravitational clustering [22] and focusing by the Sun [23]. For these reasons, in what follows we focus on decays of a light J after recombination. Below, we examine the experimental signature of neutrinos from these decays at experiments searching for the CνB such as PTOLEMY.

III. OBSERVING THE NEUTRINOS FROM DARK MATTER DECAY
A. Neutrino capture from diffuse J → νν We are interested in the detection prospects for neutrinos with E ν ∼ eV today which could show up in experiments searching for direct evidence of the CνB. The most promising technique to detect these neutrinos is neutrino capture on β-decaying nuclei. As a benchmark detector setup, we focus on the recently proposed experiment PTOLEMY [10], which hopes to use a tritium target of mass M T = 100 g and to achieve an energy resolution of about 0.1 eV (full width at half maximum) on recoiling electrons near the endpoint. This could probe neutrinos with energies as low as 0.1 eV.
The cross section for an electron neutrino with E ν keV to capture on tritium is [21] σ = 3.83 × 10 −45 cm 2 .
Using this with Eq. (4) leads to a capture rate on tritium for neutrinos from J decays, averaged over the Universe, of R cos = 4.52 yr In this expression f νe is the fraction of these neutrinos that are of electron flavor. While we have thus far not specified the flavor content of the neutrinos that J couples to, we will discuss this briefly here. Note that if J is indeed the majoron [1] it couples to the neutrino mass eigenstates with strength proportional to their masses, and thus the flavor content of these neutrinos does not oscillate [4]. First, consider the case that the neutrinos are light and not degenerate, m ν 0.1 eV. If the mass hierarchy is normal then J → ν 3 ν 3 is the dominant mode and f νe 0.03. If instead it is inverted, Γ J→ν1ν1 Γ J→ν2ν2 and f νe 1/2. If the neutrinos are relatively heavy and degenerate with m ν 0.1 eV the rate into all mass eigenstates is comparable and f νe 1/3. Different scenarios where J couples to states that are not mass eigenstates lead to flavor oscillations and can give different values of f νe .
Since we are focusing on τ J > t rec , we can use the same limit on Ω J as a function of τ J from [14] that lead to the upper limit onñ ν in Eq. (12) to limit their capture rate.
Doing so gives R cos 3.44 yr 12 Gyr τ J 160 Gyr 2.3 (160 Gyr/τ J ) , τ J 160 Gyr, which, as expected, is comparable to the rates expected from the CνB. However, as we see from Eq. (15), the neutrinos' energies can be around an eV today without requiring the neutrino masses to be large.

B. Dark matter decay in the galaxy
Thus far, we have discussed only the contribution to a neutrino capture signal from diffuse J decays averaged over the entire Universe. Since the J's act as dark matter today if τ J is comparable to or larger than the age of the Universe, their local density in the Milky Way can be greatly enhanced over the cosmic average. Using a Navarro-Frenck-White dark matter profile [24], ρ(r) ∝ r −1 (1 + r/r s ) −2 with scale radius r s = 24 kpc and local dark matter density of 0.3 GeV/cm 3 along with the capture cross section on tritium in Eq. (16), we find a rate for capture from neutrinos produced in the Milky Way of This signal is essentially monochromatic with E ν = m J /2. For τ J 5 Gyr, this rate is negligible but becomes important as τ J is increased, so that for τ J 100 Gyr it is roughly comparable to that from cosmically averaged J decays, R cos . In Fig. 2, we show contours for a capture rate on tritium of 1 yr −1 (100 g) −1 [10 yr −1 (100 g) −1 ] as solid [dotted] curves in the Ω J /Ω dm vs. τ J parameter space, for (from right to left) m J = 1, 3 eV. We have chosen f νe = 1/2, appropriate if J were a majoron decaying to light neutrinos with an inverted mass hierarchy. Note that these are raw signal rates and do not impose a cut to remove background from β-decay. In addition, we shade the 2σ exclusion region on a decaying dark matter component from the analysis of CMB data in Ref. [14]. We also point out the two benchmark points in this parameter space whose fluxes were shown in Fig. 1: m J = 1 eV, τ J = 7 Gyr, Ω J /Ω dm = 0.08 as a green club and m J = 3 eV, τ J = 10 11 yr, Ω J /Ω dm = 1 as a red spade.
Contours showing capture rates on tritium of 1 yr −1 (100 g) −1 [solid] and 10 yr −1 (100 g) −1 [dotted] in the space of τJ and the J energy density in units of the total dark matter energy density, ΩJ /Ω dm (J acts as a decaying dark matter component). mJ has been chosen to be, from right to left, 1 (green) and 3 (red) eV. In all cases we have taken the electron neutrino fraction of the J decays to be fν e = 1/2. The 2σ exclusion from CMB data in [14] is the shaded gray region. We also show the locations of the two benchmark points which are plotted in Figs. 1 and 3: mJ = 1 eV, τJ = 7 Gyr, ΩJ /Ω dm = 0.08 (green club) and mJ = 3 eV, τJ = 10 11 yr, ΩJ /Ω dm = 1 (red spade).

C. Signature
Just as is the case with the CνB, detecting the neutrinos from J decay is challenging because the signal of an electron above the β-decay endpoint can be swamped by the enormous rate from β-decay due to the finite e − energy resolution. To mitigate this background, PTOLEMY aims for roughly 0.1 eV resolution on the electron energy to be able to observe the CνB for neutrinos with m ν 0.1 eV. For reference, Dirac (Majorana) CνB neutrinos would lead to a rate of 4 (8) yr −1 for 100 g of tritium [21]. We see that J → νν decays at late times can lead to comparable event rates with E ν above threshold and, unlike the CνB, this can even be the case if m ν < 0.1 eV.
To illustrate the signal, we show the e − kinetic energy distribution for capture on tritium in Fig. 3, assuming a resolution of 0.1 eV as proposed by PTOLEMY, for two benchmark points currently allowed by cosmological data, m J = 1 eV, τ J = 7 × 10 9 yr, Ω J /Ω dm = 0.08 (solid, green) and m J = 3 eV, τ J = 10 11 yr, Ω J /Ω dm = 1 (solid, red) taking f νe = 1/2 as above. These benchmark points are also marked in Fig. 2 and their cosmically averaged fluxes are shown in Fig. 1. Also shown in Fig. 3 is the background from the standard tritium β-decay with m ν = 0 as a solid gray line. The rates are shown as functions of the e − kinetic energy K e minus the m ν = 0 endpoint energy K end . The m J = 1 eV and 3 eV benchmarks here give rates of 4.8 and 7.0 yr −1 (100 g) −1 for K e − K end > 0.2 eV, respectively, compared to 0.66 yr −1 (100 g) −1 from the β-decay background. As a comparison, we also plot the total rate in the case of a standard CνB signal with a m ν = 0.15 eV Majorana neutrino as a dashed yellow curve.
Note that, for the diffuse components of the signals, we used, as in Fig. 1, the energy spectrum from Eq. (13), Ω b = 0.05, Ω Λ = 0.69, Ω dm = 0.26 with the dark matter split up into a decaying and nondecaying component, Ω J and Ω dm − Ω J , respectively. The contribution from J decays within the Milky Way gives bumps with a width determined by the detector resolution. These two sources of neutrinos lead to the interesting signature of a peak at K e −K end = m J /2 from local J decays and a shoulder extending from m J /2 down to the endpoint from redshifted diffuse decays, with roughly comparable signal strength in each. The size of the signal from local J decays is subject to uncertainty coming from our ignorance of the precise dark matter distribution of the Milky Way as well as its overall normalization; observation of a nonzero signal could of course help shed light on this issue.

IV. CONCLUSIONS
This paper has explored a new physics scenario that searches for the CνB can impact: light dark matter that decays to neutrinos. While the event rates at a PTOLEMY-like experiment in this scenario are not enormous and the range of m J accessible is not large, over a spread of several eV, a comparable signal strength to those expected from the CνB is possible. Furthermore, the signal described here is distinct from that of the CνB, and can potentially be observable even if the e − energy resolution is poorer than what is needed for CνB detection given the current upper bound on m ν .
While a robust statistical analysis of possible signals is beyond the scope of this study, if the electron energy resolution can be kept to the roughly 0.1 eV level, excesses beyond the β-decay background could potentially be seen. Seeing such a signal would of course be tremendously exciting, opening up a world of new physics implications. Even providing an upper bound on the rate, however, would constrain part of the parameter space of dark matter decaying to neutrinos that is not currently probed by cosmological observations. This would be useful information, especially in light of nagging tensions in the cosmological data, such as the discrepancy in the extraction of the Hubble constant between CMB and local observations and the σ 8 problem. This signature also does not suffer from the same degeneracies that affect the extraction of conclusions from cosmological observations. Furthermore, in the event of a positive signal, the two distinct sources of neutrinos, from J decays averaged over the Universe and those in the Milky Way could allow for the dark matter distribution of the Milky Way to be studied in more detail.
We mention here the possibility that similar models with a higher multiplicity of neutrinos in the decay of dark matter χ, e.g. χ → 2J → 4ν, could lead to larger event rates. However, this would come at the cost of diluting the energy of the neutrinos today, potentially masking the distinguishing spectral features in the capture rate. Other models, more complicated than the simple one described here, with potentially more involved cosmological histories could also lead to larger enhancements.
While we have not discussed them in detail, other proposals to search for the CνB, such as using laser interferometers [25], that do not involve looking for electrons above the β-decay endpoint could also be impacted by the decay of dark matter to neutrinos. While potential signals in such a setup are not likely to lead to large effects, assessing the sensitivity of new techniques in general to the signal described in this paper would be interesting. Moreover, other interesting features to consider would be the effects of gravitational clustering of the neutrinos from J decay and potential anisotropies in the signal from decays in the Milky Way, along the lines explored in [26].