Sterile neutrino dark matter catalyzed by a very light dark photon

Gonzalo Alonso-Alvarez´ ∗ and James M. Cline† McGill University, Department of Physics, 3600 University St., Montr´eal,QC H3A2T8 Canada

Sterile neutrinos (νs) that mix with active neutrinos (νa) are interesting dark matter candidates with a rich cosmological and astrophysical phenomenology. In their simplest incarnation, their production is severely constrained by a combination of structure formation observations and X-ray searches. We show that if active neutrinos couple to an oscillating condensate of a very light Lµ Lτ − gauge field, resonant νa-νs oscillations can occur in the early universe, consistent with νs constituting all of the dark matter, while respecting X-ray constraints on νs νaγ decays. Interesting deviations from standard solar and atmospheric neutrino oscillations can→ persist to the present.

I. INTRODUCTION the expansion of the universe eventually shuts off the resonance at late times. 1 The possibility that sterile neutrinos νs constitute the The formation of bosonic condensates is conjectured dark matter (DM) of the universe has attracted steady to be possible in the early universe. The most well-known interest since its inception long ago. Supposing that the examples are those of the QCD axion [30–32] and ul- νs abundance is initially negligible, its relic density can tralight scalar dark matter [33], but an analogous phe- be generated by nonresonant oscillations with an active nomenon is possible for light vector bosons [34]. The con- neutrino species νa though mass mixing [1], with a mixing densate can be described in terms of the classical bosonic −6 angle as small as θ 10 for mνs 100 keV [2]. How- field, which acquires a large primordial vacuum expecta- ever, these values of∼ active-sterile∼ mixing are robustly tion value (VEV) that is initially frozen by Hubble fric- excluded by astrophysical observations. In particular, tion. Once the Hubble scale H becomes smaller than the X-ray searches for the decay νs νaγ, which arises at field’s mass, the field starts oscillating around the mini- → one loop due to the mixing, rule out mνs & 3 keV [3]. mum of its potential, the condensate thereafter diluting This forces nonresonantly produced νs into the warm DM like a pressureless, non-relativistic fluid. regime that is strongly excluded by satellite counts [4–7] There are multiple ways in which a vector condensate and Lyman-α constraints [8–13]. can be produced in a cosmological environment. The These constraints are relaxed in the presence of a large misalignment mechanism, which is most well-known in primordial lepton asymmetry, which can allow for reso- its scalar version [30–32], fails to generate a VEV in the nant oscillations [14–20] in the early universe. But even vector case unless direct couplings to the curvature are in this case the allowed window of parameter space is present [35, 36]. Any vector field is however unavoid- small [21], considering the constraints from Big Bang Nu- ably sourced through quantum fluctuations during infla- cleosynthesis (BBN) on the lepton asymmetry and from tion [36–40], producing a cosmological population that structure formation [22, 23]. Various alternative ways can at late times be effectively described by a nonzero of producing the required νs abundance have been ex- VEV. Finally, an effectively classical vector field config- plored, including decays of Higgs singlets [24–26], sce- uration can originate from the decay of a dark sector narios in which the νs population is colder relative to (pseudo)scalar [41–45]. production through generic oscillations [27], through de- The remainder of this paper is structured as follows. cays of a MeV-GeV scale vector boson [28], or through Firstly, we describe how nonstandard interactions of ac- oscillations modified by neutrino self-interactions [29]. tive neutrinos with a background vector condensate can Here we propose a new mechanism for achieving the resonantly produce νs in Section II. In Section III, the desired νs abundance through resonant oscillations, by arXiv:2107.07524v1 [hep-ph] 15 Jul 2021 identification of this vector with a U(1)L −L gauge bo- coupling the active neutrinos to a very light dark photon, µ τ −12 son is explained, along with the numerous laboratory and with a mass as low as mA0 10 eV. The presence of cosmological constraints on its couplings. Detailed pre- a primordial condensate of such∼ vector boson particles is dictions for the νs relic density are made in Section IV. shown to modify the dispersion relation of the neutrinos Then, a preliminary study is made in Section V of the im- in the early universe. The leading effect is an enhance- pact of the dark photon dynamics at the present time on ment of the active neutrino self-energy, which enables a active neutrino oscillations. A summary and conclusions level crossing with its heavier sterile counterpart and thus are given in Section VI. facilitates resonant oscillations between the two states. The dilution of the vector boson condensate caused by

1 Here, the word condensate should be understood as a coherent ∗ [email protected]; ORCID: 0000-0002-5206-1177 classical field configuration and not necessarily as a Bose-Einstein † [email protected]; ORCID: 0000-0001-7437-4193 condensate in the quantum mechanical sense. 2

II. RESONANT OSCILLATION MECHANISM 103 6 mA0 = 10− eV

] Electroweak 2 5 0 g0A0 = 10− GeV A necessary condition for resonant νa νs oscillations 10

− [GeV

to occur is for the the matter potential of the active neu- | p trinos νa to be increased and cross the sterile neutrino 2 3 Dark photon · 10− ) mass threshold (here and in what follows we always as- 0 Total ˆ A sume that ma ms). For simplicity, here we illustrate · 6  T 10−

the situation for sterile mixing with a single active flavor. −

In practice and due to the smallness of the involved mix- = ~p

( 9 ing angles, the more general results for mixing with mul- a 10− V tiple flavors can be easily deduced by adding the individ- | 2 2 ms = (10 keV) ual single-flavor contributions. The Schr¨odingerequation 12 10− 2 1 0 1 2 describing the oscillations is governed by the Hamiltonian 10 10 10 10− 10− T [GeV] 2 2 1 maa mas Va 0 Has = p + 2 2 + (1) 2p mas mss 0 0 Figure 1. Absolute value of Va, the matter potential for active     neutrinos, times 2p, twice the neutrino proper momentum, as in the (νa, νs) flavor basis and in the relativistic limit. a function of the temperature of the universe. The green line Here, p is the neutrino momentum and Va is the corre- shows the EW potential in Eq. (2), the purple line corresponds sponding matter potential. Resonant flavor transitions to the one induced by the dark photon condensate, Eq. (5), 2 2 2 and the orange line shows the sum of both contributions. occur when 2pVa +maa = mss ∼= ms, indicating the need Solid (dashed) lines indicate where Va or its components are for a large positive contribution to Va which does not 2 exist within the standard model (SM). In fact, the elec- positive (negative). Resonances occur when 2pVa crosses ms, the square of the sterile neutrino mass (drawn as a red line troweak (EW) matter potential in the absence of lepton for an exemplary value of ms = 10 keV). For this figure, pa- asymmetries, which we denote by VW , is negative and −6 0 0 −5 rameters were fixed at mA0 = 10 eV and g A0 = 10 GeV, given by with the neutrino momentum taken to be antiparallel to the ˆ0 0 vector VEV, ~p = T A . A0 denotes the initial amplitude of 14π 2 2 2 4 2 4 the dark photon field− before it starts oscillating at H m. VW = (3 s )s G T p C T p, (2) ∼ −45α − W W F ≡ − W ∼ at temperatures T below the electroweak phase transition temperature TEW but higher than the mass of the cor- responding charged lepton [46]. Here, sW = sin θW is the can directly read off the new contribution to the matter sine of the Weinberg angle and we have taken the neu- potential, which becomes Va = VW + δVa, with trino to be approximately on shell (E ∼= p) to remove the energy-dependence from the general expression for the 1 0 0 2 0 0 self-energy. To accurately treat the system at all tem- δVa = (g A ) 2g ~p A~ . (5) 2p − · peratures and momenta, rather than using the limiting   expression for T TEW in (2), we take the general form  The total matter potential for an exemplary choice of of VW described in appendix A in the following numerical model parameters is shown in Fig. 1. The time evolution treatment. of the dark photon field is obtained from its classical A positive contribution to V can arise by coupling the a equation of motion active neutrino to a new bosonic field with a nonvanishing VEV. SU(2)L gauge invariance motivates us to consider ¨0 ˙ 0 2 0 0 A + 3H(t)A + mA0 A = 0, (6) a vector boson Xµ with the gauge coupling

0 0 which is solved by = L¯a i∂/ g X/ ma La , (3) L − −   0 0 J1/4(mA0 t) where La is the lepton doublet. Assuming that the VEV A (t) = 1.08 A0 1/4 , (7) ∼ (m 0 t) is spacelike and using the gauge freedom to set X0 = 0, A the dispersion relation in the absence of mixingh andi in where J1/4 is a Bessel function of fractional order, and the relativistic limit is 0 we denote by A0 the amplitude of the field at early times, 0 1 0 0 2 0 ~0 2 when A is still frozen by Hubble damping. As can be E = p + VW + (g A ) 2g ~p A + ma , (4) 2 ∼ 2p − · seen in Fig. 1, the condition for resonance, Va = ms/2p,   can occur when the the dark photon contribution almost 0 0 i.e. where A = X /a is magnitude of the physical vector exactly cancels the electroweak one, , δVa ∼= VW , field, a is the cosmological scale factor, and p is the neu- or when the vector potential oscillates at temperatures− trino proper momentum [35]. From this expression we H(T ) . mA0 . 3

6 III. COUPLING THE DARK PHOTON TO 10−

NEUTRINOS 8 10− Neff > 0.5 0 g 10 10− To couple A0 to active neutrinos in an anomaly-free Stars 12 (Higgs) 10 (LISA) way, the simplest phenomenologically viable choice is to − 13 GeV SN1987a 10 Superradiance Superradiance gauge Lµ Lτ . Although Le Lτ and Le Lµ are also 14 SSB ' − − − 10− decay (CMB)T anomaly free, there are very stringent constraints on fifth gauge coupling ν τ 16 0 L 10− forces, such as that mediated by a light A coupling to −

µ (St¨uckelberg)

0 L 18 electrons [47], which would rule out couplings g down to GeV 10− 13 (1) −25 −20 10 the 10 -10 level. In the Lµ Lτ case, fifth forces U 20 NS binaries − 10− ' between the muons present in neutron stars can alter Λ UV  = g0/70 22 the orbital period and the GW emission of neutron star 10− 20 18 16 14 12 10 8 6 4 2 0 2 4 10− 10− 10− 10− 10− 10− 10− 10− 10− 10− 10 10 10 binaries [48], leading to strong bounds, which however Dark photon mass m [eV] −10 A0 only apply to very light (mA0 . 10 eV) mediators. We therefore focus on Lµ Lτ gauge bosons with masses − Figure 2. Constraints on a U(1)Lµ−Lτ gauge boson in the m 0 1 eV in what follows. A . mass vs. gauge coupling plane, including the loop-induced ki- 0 The U(1)Lµ−Lτ gauge symmetry forbids mixing be- netic mixing with  = g /70 [54]. Bounds are shown from fifth tween the active neutrinos in the absence of spontaneous forces on neutron star binaries [48] (purple), modifications of symmetry breaking. It is therefore necessary to intro- Neff during BBN [55–57] (blue), observations of supernova duce extra scalars charged under the symmetry that can SN1987A [58, 59] (green), imprints of neutrino decay in the generate the observed neutrino mixing pattern after ac- CMB [60–64] (red, see text for explanation of the different quiring suitable VEVs. Many proposals in this direc- shades), stellar emission [65–68] (brown), and black hole su- tion have been put forward involving extra Higgs dou- perradiance [69–71] (yellow, the lighter band labelled “LISA” is a forecast and not a current exclusion). Above the dashed blets [49], soft-breaking terms [50] and/or right-handed black lines, the temperature of spontaneous symmetry break- neutrinos that are also assumed to carry Lµ Lτ charge − ing (for a Higgsed dark photon) or the UV cutoff of the theory [51–53]. The common element of all these models is the (for a St¨uckelberg mass) are lower than 1013 GeV, which may presence of U(1)Lµ−Lτ symmetry-breaking scalars, which interfere with the production of the dark photon during high- inevitably contribute to the Lµ Lτ dark photon mass. scale inflation. −

+ − 0 2 µ µ γA leads to an increase in Neff [56]. At smaller III.1. Constraints on the Lµ Lτ gauge coupling ↔ − mA0 , the enhanced interaction of the longitudinal com- ponent of the vector with neutrinos makes νν A0A0 → A variety of astrophysical constraints on g0 are shown the most efficient process to produce dark photons in in Fig. 2. First, the duration of the observed neutrino the early universe [57], yielding stronger limits with de- burst from supernova 1987A constrains the amount of creasing mA0 , shown by the blue line in Fig. 2. Al- energy that could have been carried away by A0 from though these early-universe limits are currently weaker ± the supernova. The abundance of µ and νµ in the than those from SN1987A, they have the potential to sig- proto-neutron star allows a sufficiently strongly coupled nificantly improve, in light of the expected sensitivity of Lµ Lτ gauge boson to be copiously produced and de- next-generation Cosmic Microwave Background (CMB) − plete energy from the exploding star. The production of experiments to deviations in Neff at the 0.03 level [72]. 0 0 A is dominated by semi-Compton scattering for mA0 The stronger coupling of the longitudinal A polar- 10−4 eV [59] and by neutrino pair-coalescence for lighter∼ ization also enhances the rate of active neutrino decays 0 masses [58]. The constraint strengthens with decreasing νi νj + A for mA0 mν , which can distort the →  i mA0 (green region of Fig. 2) because the coupling between CMB angular and frequency spectra [60, 61]. A reeval- neutrinos and the longitudinal dark photon mode scales uation of this effect in light of Planck 2018 data [62] 0 9 3 as g mν /mA0 at energies E mA0 . This is a consequence obtained the limit τνi & 10 s (mνi /50 meV) , strongly  of the gauge current nonconservation in the presence of constraining Lµ Lτ models [63]. This was recently dis- − spontaneous symmetry breaking, as is necessary to allow puted in Ref. [64], which found a much weaker constraint 5 6 5 for ν ν oscillations, and can be understood using the τν (4 10 4 10 ) s (mν /50 meV) . In view µ τ i & × − × × i Goldstone− boson equivalence theorem [57]. of this discrepancy, we conservatively impose the weaker Similar conditions occur in the early universe, where neutrino interactions with the dark photon can alter the effective number of degrees of freedom (Neff ) and 2 This is only possible if the universe reheats above the muon 0 spoil the light-element abundance predictions of Big threshold (TR > mµ); otherwise A production relies on the in- −6 Bang Nucleosynthesis. For masses mA0 & 4 10 eV teraction with neutrinos [55], which yields bounds that are about and couplings g0 4 10−9, thermalization× through 3 orders of magnitude weaker for high-mass dark photons. & × 4 constraint, corresponding to the dark red region in Fig. 2, As we will see in the following sections, the field values while the lighter shading depicts the stronger limit. corresponding to these inflationary-sourced abundances Beyond the direct limits on g0, loop processes involving are largely compatible with the ones that favor sterile µ and τ leptons induce kinetic mixing between A0 and γ of neutrino dark matter production. strength  g0/70 [54]. This allows A0 to couple to elec- trons, giving∼ additional constraints from stellar emission of dark photons [65–68] that become relevant for larger III.2. Higgs versus St¨uckelberg mass masses mA0 & 0.1 eV. Finally, observable consequences of the superradiance of spinning black holes triggered by For U(1) theories, it is consistent to explain a gauge bo- a light vector can rule out the existence of light bosons son mass using either the Higgs or the St¨uckelberg mech- in a given mass window [69–71], independently of their anism. The generation of the observed neutrino mixing nongravitational interactions.3 Fig. 2 displays the limit 2 matrix implies that at least part of mA0 should arise from from Ref. [71] in dark yellow, while the lighter yellow re- the Higgs mechanism, but it is logically possible that gion indicates the range of mA0 that could be probed by both mechanisms contribute, and the latter could domi- future GW observatories like LISA [75]. nate. Although the low-energy phenomenology is equiva- Fig. 2 shows that the CMB constraints on decaying lent for either case, important differences arise when the 0 neutrinos require g to be very small for dark photons UV completions of the model are considered. in the mass region of interest. However, even such small 0 0 For the case of a Higgsed vector, the manifestly gauge- couplings are easily compatible with an initial g A0 that invariant Lagrangian can be written as can copiously produce sterile neutrinos by the mechanism that we propose. To characterize the allowed range of 2 2 0 0 0 0 2 2 mφ g A0 values, we note that the energy density of the vector Higgs = ∂µ + ig Aµ φ + λ φ , (11) condensate A0 redshifts like dark matter once the field L | | − 2λ ! h i
starts oscillating at a temperature Tosc. The relic density in A0 is related to its mass and the initial field value where φ is the Higgs field with mass mφ and quartic self- h i 0 by [35] coupling λ. At temperatures below TSSB mφ/g , spon- taneous symmetry breaking occurs and∼ the scalar field 0 2 ΩA0 A0 mA0 acquires a VEV, vφ = mφ/√2λ. This in turn generates 5 (T ) , (8) 0 osc 12 0 ΩDM ' F 10 GeV eV a mass for the gauge boson, mA = g vφ. Note that vφ is r 0   completely determined by g and mA0 , but freedom still 3/4 where (Tosc) = [g?(Tosc)/3.36] /[g?S(Tosc)/3.91] is an remains to vary the Higgs mass and self-coupling as long F (1) function encoding the dependence on the number of as the ratio mφ/√2λ is kept fixed. O relativistic degrees of freedom. For a given value of mA0 , As an example, for the benchmark values mA0 = −6 0 −14 the condition that the relic density of dark photons does 10 eV and g . 10 , motivated by the discus- not exceed a given fraction of ΩDM sets a limit on the sion in the next section, one finds vφ & 100 MeV and 0 13 size of the initial A0. Combining this with the previously TSSB & √2λ 10 GeV. For an (1) quartic coupling, described limits on g0 constrains the dark photon· would already beO massive during high- 0 scale inflation, as is necessary to generate the initial A0 1/2 3/4 0 0 1 ΩA0 mA0 during that epoch. This is also the case for any Hig- g A0 < 0.5 GeV −6 , ΩDM 10 eV gsed dark photon with mass and gauge coupling below (Tosc)   F   (9) the corresponding dashed black line in Fig. 2, which in- p −10 −2 which is valid for 10 eV . mA0 . 10 eV, for which cludes most of the viable parameter space. We thus con- the CMB constraint on ν decays dominates the limits on clude that a dark photon mass generated via the Higgs the size of g0 (c.f. Fig. 2). mechanism is largely compatible with the requirements Several mechanisms have been proposed to seed the ini- of the sterile neutrino production dynamics under con- tial dark photon abundance. Perhaps the simplest and sideration. most general one is the gravitational production of lon- The other possibility is the St¨uckelberg mechanism, gitudinal modes during inflation [37], which yields a relic which can be implemented via the Lagrangian density of dark photons that depends on the Hubble rate 1 2 0 0 2 during inflation as St¨uckelberg = f ∂µθ + g A . (12) L 2 θ µ 2 ΩA0 mA0 HI
−6 14 . (10) Gauge invariance is recovered by demanding that the ΩDM ' 6 10 eV 10 GeV r ×   St¨uckelberg field θ shifts appropriately under a gauge transformation. In this case, the dark photon mass 0 mA0 = g fθ is proportional to the decay constant fθ. 3 Interactions of the dark photon with the plasma around astro- The main difference with the Higgs case is that in the physical black holes may lead to a weakening or nullification of St¨uckelberg model the mA0 0 limit lies at an infinite the constraints [73, 74]. distance in field space, in a→ supersymmetric context. It 5 was argued in [76] that this leads to restrictions in the take into account the decoherence effects due to elastic dark photon parameters based on quantum gravity con- scatterings, we perform the time average jectures. In particular, the Swampland Distance Conjec- −Γat ture [77] postulates that traversing large distances in field Pa→s = Γa dt e Pa→s space necessarily causes a tower of states to become light h i Z 2 4 and invalidates the effective theory under consideration. 1 θ0ms = 2 2 2 2 4 .(17) In our case, this is conjectured to occur when approach- ∼ 2 (pVa m /2) + (pΓa/2) + θ m − s 0 s ing the m 0 0 limit, motivating the introduction of a A Note that elastic scatterings suppresses the oscillation UV cutoff of→ the effective theory, probability if their rate Γa is larger than the oscillation frequency ω. It is therefore of crucial importance to cor- mA0 MP rectly estimate the rate of active neutrino scattering in ΛUV . . (13) s g0 the early universe plasma. While the approximation (15) for Γa is convenient for analytic estimates, for our numer- To maintain predictivity for the inflationary production ical analysis we use the more accurate results that have of dark photons, we should demand this cutoff to lie been tabulated in Refs. [85, 86]. We find that using the above the scale of inflation. For high-scale inflation, this more precise Γa evaluation can correct the sterile neu- is only possible if the gauge coupling lies below the lower trino yield by up to an order of magnitude. dashed black line in Fig. 2. Although this argument has The relative abundance of sterile versus active neutri- a number of possible loopholes,4 it tends to disfavor the nos of momentum p at a given temperature is commonly St¨uckelberg mechanism as the origin of the dark photon denoted by R(p, T ). This quantity satisfies the Boltz- 5 mass in the parameter space of interest for νs dark matter mann equation dR/dt = Γs(1 R), whose solution is production. − ∞ dT 0 Γ R(p, T ) = 1 exp s . (18) − − H(T 0)T 0  ZT  IV. STERILE NEUTRINO DARK MATTER The total abundance of sterile neutrinos relative to active ones is obtained by the phase space integration To compute the abundance of resonantly produced sterile neutrinos, including matter effects, a rigorous ap- d 3p R(p, T )f(p/T ) Y = , (19) proach is to use the density matrix form of the Boltz- s/a d 3p f(p/T ) R mann equations [79–81]. However, this is computation- ally demanding and accurate results can be found though where f(p/T ) denotes theR Fermi-Dirac distribution func- a simpler, intuitive procedure [82–84], based on the pro- tion. Once the relative abundance Ys/a is known, the duction rate of νs through νa-νs oscillations. This rate number density of sterile neutrinos can be written as can be written as ns = Ys/anνa , where the active neutrino abundance is −3 nνa = (3/77)s, and s = 2891 cm is the current entropy Γs = Γa Pa→s (14) density. The fractional contribution to the cosmic energy · h i density is thus Ωs = msns/ρcrit. and is proportional to the rate of elastic νa scattering below the weak scale [46], IV.1. Evaluation of the sterile neutrino abundance 7π Γ G2 T 4p C2T 4p, (15) a ∼= F a 24 ≡ In general, it is numerically challenging to evaluate the abundance (19), but several simplifications help. First, and to the time-dependent oscillation probability Pa→s = 2 2 since we are interested in small abundances, we can lin- sin (ωt/2) sin (2θm). Here, the matter mixing angle is earize the exponential in (18) and then reverse the order 2 of the p and T integrals in (19). After this, the prob- 2 4θ0 sin (2θm) = , (16) ∼ 2 2 2 ability (17) can be integrated exactly over the angle α 4θ + (2pVa/m 1) 0 s − between A~ and ~p to give 2 2 where θ0 = mas/ms is the vacuum mixing angle and 1 2 ∼ 2 2 2θ0 −1 2ab ω = Va + (ma ms)/2p is the oscillation frequency. To d cos α Pa→s = tan 2 2 2 ,(20) − −1 h i ab a b + c Z  − 

4 5 The conjectures referred to above have not yet been proven, and This expression assumes that T & 5 MeV so that the active neu- it may be possible to model-build them away [78]. Moreover, trinos are still in equilibrium and total neutrino number is not inflation may occur at a lower scale or the production of dark conserved. I practice, since we are interested in solutions where photons may take place at a lower energies through a different R  1, this caveat is numerically unimportant even at lower process. temperatures. 6 where that a 2θ0 and therefore Ys/a θ0. This expectation is confirmed' by the numerical results∝ described below. 2 0 0 2 pΓa g A p The dependencies of the sterile neutrino yield on the a = 4θ0 + 2 , b = 2 2 , s ms ms, other parameters—the νs mass ms, the dark photon mass   0 0 0 2 mA , and the product of the gauge coupling and the ini- g A pVW 0 0 c = 1 2 . (21) tial field value g A0—are more complex and their under- − m − m2  s  s standing requires a full numerical study, to which we now turn. With this, Ys/a and the resulting cosmological abundance can be written as 2 2MP θ0 0 0 IV.2. Numerical results Y = Is(ms, mA0 , g A ), (22) s/a 3 1.66 ζ(3) 0 · 2 2 ms The results of the evaluating Eq. (24) numerically, be- Ωsh = 3.6 θ MP Is , (23) ∼ 0 keV fore and after doing the integration over T , are shown in   Figs. 3 and 4. Fig. 3 illustrates the temperature evolution where the 2d integral of the differential sterile neutrino abundance, demon- 2 strating that bulk of the production occurs at temper- q f(q)Γa −1 2ab Is = dT dq tan (24) atures T 1 GeV. At this time the Hubble rate has g (T ) T 3ab a2 b2 + c2  Z ?  −  already fallen below 0.6 mA0 , after which the vector field is no longer frozen∼ but is instead oscillating around must in generalp be computed numerically. We can how- the minimum of its potential at A0 = 0 [36]. The tem- ever make use of the narrow width approximation to gain perature at which this transition occurs is some further analytic understanding. In terms of the angle between ~p and A~0, the resonance 2 0 0 0.6 mA0 Mp occurs near p cos α = ms/ g A , and is narrow if Tosc , (27) | | ≡ s 1.66 √g∗ ∆(p cos α) 2 2 2 = (pΓa/m ) + 4θ 1. (25) p cos α s 0  and for T < Tosc, the amplitude redshifts as q 3/2 Since we are only interested in θ0 1, the condition 0 0 T  2 A A0 . (28) for narrow resonance is simply pΓa ms. Given that ' Tosc the phase space integral introduced above is dominated   3 −6 by p T , this implies GF T ms, hence T 1 GeV The benchmark value m 0 = 10 eV of the dark pho- ∼ . . A for ms = 10 keV. This is satisfied for the parameters of ton mass chosen in Fig. 3 corresponds to Tosc 22 GeV, ∼ ' interest here, and allows one to approximate the integral at which point g∗ = 87.2 [87]. The sterile neutrino pro- Eq. (20) by duction thereafter occurs at the subsequent narrow reso- nances when the dark photon field value passes through 1 2 θ0 zero and the resonant condition b c is satisfied. d cos θ Pa→s ∼= 2π Θ( b c ) . (26) | | ≥ | | −1 h i ab | | − | | These events are responsible for the bulk of the sterile Z | | neutrino production and are visible as sharp vertical fea- The q-integral in Eq. (24) can then be carried out an- tures in Fig. 3. alytically in terms of polylog functions, but the task The basic properties of Figs. 3 and 4 can be understood of solving for the limits of integration (the roots of in terms of the conversion probability given in Eq. (17). b = c ) becomes too numerically demanding once the Production of νs is most efficient when the denominator vector| | | field| unfreezes and starts oscillating. This eventu- is minimized, which occurs when the resonance condition 2 ally makes this procedure more costly than the straight pVa = ms/2 is satisfied, and when the full width of the 2 2 4 1/2 two-dimensional numerical integration. resonance ((pΓa/2) + θ0ms) is small. Since the aver- Nevertheless, the expression (26) is helpful to under- age momentum is p 3T , the latter condition is fulfilled ∼ stand the dependence of the sterile neutrino yield on the for some optimal temperature T∗ given by vacuum mixing angle θ0. The width of the resonance is 1 determined by a competition between the vacuum mixing 2 6 3 2 θ0ms T∗ Γa(T∗) θ0 m = T∗ . (29) angle and the active elastic scattering rate, as can be seen 2 ∼ s ⇒ ∼ C2 in Eq. (17). This fact singles out two distinct regimes in  a  the calculation. For resonances occurring at high T the For the fiducial parameters shown in Fig. 3, this temper- 2 scattering rate is large, so that typically 2θ0 pΓa/ms. ature is T∗ 0.06 GeV. It is important to note that the 2  ' In this case, a pΓa/ms and Eq. (26) implies that the estimate (29) is only reliable if resonant conversions are ' 2 abundance of sterile neutrinos scales as Ys/a θ0. On occurring at T T∗. This is only the case if (i) the dark the other hand, the width of resonances occurring∝ at low photon field has∼ already begun oscillating at that tem- T is controlled by the vacuum mixing angle, implying perature (i.e. if Tosc > T∗) and if (ii) the amplitude of the 7

100 105 m = 1.0 10 6 eV A0 − 5 4 2 11 DM × g0A0 = 10− GeV 10 θ0 = 10− 2 12 Ω 2 θ0 = 1.0 10− 10 3 2 12 − × θ = 10− >

10 0 1% 0 ms = 10 keV 7 10% 3 g0A0 = 10− GeV 2 ms = 100 keV A

g0A0 = 10− GeV 10 Ω 10 4 0 − 1 1 g0A0 = 10− GeV 10 s/a

T 0 m = 10 keV

DM s

Y 10

d 6

10− Ω d

/ 1 s 10 T − Ω 8 2 ms = 1 keV 10− 10− 3 10− 10 4 10− 10− 5 10− 12 10− 6 10 4 10 3 10 2 10 1 100 101 10− − − − − 10 4 10 3 10 2 10 1 100 101 102 103 104 105 T [GeV] − − − − 6 3/4 g0A0 10− eV X 5 ≡ 10− GeV mA0 Figure 3. Differential sterile neutrino abundance Ys/a pro-     duced as a function of temperature, for several initial values 0 0 of g A , and fixed mA0 , θ0, and ms as indicated in the figure. Figure 4. Integrated sterile neutrino relic density Ωs as a func- The vertical jumps in the curves are due to fast oscillations 0 0 −3/4 tion of the product X g A0 mA0 , for several values of ms in the dark photon field. 2 ∝ · and θ0. In the area shaded in grey the dark photon overcon- tributes to the dark matter energy density, the dashed and dotted lines indicate a minimum 10% and 1% contribution 0 oscillations around T∗ is still large enough such that the given the constraints on g and mA0 shown in Fig. 2. sterile neutrino mass threshold can be overcome, which 0 0 2 0 requires that g A (T∗) & ms/T∗. When A becomes too small to satisfy the resonance condition, sterile neutrino the DM. production falls quite abruptly, as Fig. 3 shows in the 2 0 0 To project our findings onto the ms-θ0 plane, we low-T regions. For small values of g A0, like the one cor- responding to the blue curve in Fig. 3, resonant conver- marginalize over the dependence of the mechanism on sions halt too early to efficiently produce a large sterile the dark photon parameters, determining at given val- ues of ms and θ0 the maximum νs abundance that can neutrino abundance. 0 0 0 be generated by varying g A and mA0 . The result is Since νs production is determined by A at the time of 0 resonance, the dependence on dark photon parameters is shown in Fig. 5, in which the white region indicates where our mechanism can successfully generate νs constituting only through the combination 0 all of the DM. In the grey areas, the νs + A combined 3/4 g0A0 10−6 eV abundance falls above or below the observed DM den- X 0 g0A0 . (30) −5 T >Tosc sity. The region excluded by X-ray searches [3] is shown ≡ 10 GeV m 0 ∝ |    A  in red, while the purple dashed line depicts the poten- Hence, in Fig. 4 we choose to show the integrated νs tial reach of the proposed Athena mission [88]. In the abundance as a function of this quantity. There, one can orange region, sterile neutrinos are too light to be bound see that Ωs is maximized at some intermediate value of in the DM halos of dwarf spheroidal galaxies [89], and X, which depends on ms through the resonance condi- in the green-dotted one they are likely to be inconsistent 0 0 2 with observations of the Lyman-α forest [90], although a tion g A ms/T . One can understand the trends in Fig. 4 leading∼ to the peaks: at large X the yield is a de- reevaluation of the constraints using the distorted νs mo- creasing function of X, since larger g0A0 produces a nar- mentum spectrum in our mechanism would be necessary rower resonance, while at small X, the resonance condi- in order to establish a robust limit. tion stops being fulfilled before the optimal temperature It is worth noting that the contour along which T∗ is reached, resulting in less νs production. the nonresonant Dodelson-Widrow (DW) mechanism [1] The nonmonotonic nature of the curves in Fig. 4 im- would give the observed relic density, shown in Fig. 5, plies that there can be two values of X consistent with does not quantitatively agree with previous evaluations νs constituting the observed relic density. For sufficiently in the literature (see, e.g., Fig. 14 in [3]). We have re- small X, resonant production no longer occurs, and in- computed it as a by-product of our analysis, since it cor- stead the nonresonant Dodelson-Widrow mechanism [1] responds to the g0 0 limit of our model, finding 10 → ∼ dominates. This corresponds to the flat regions in Fig. times larger abundance in the high-ms region relative to 4, which are insensitive to the properties of A0. For very previous analyses of DW production [3]. The discrep- large X, the A0 itself starts to constitute a significant ancy comes from the use of the simplified approximation part of the dark matter, leading to the ruled-out grey (15) for the active neutrino damping rate Γa in previous part of the figure, as well as intermediate regions where literature, instead of the more quantitative results from A0 contributes a subdominant yet significant fraction of Ref. [85, 86] that we have employed. 8

10 8 − Nonresonant production Phase 10 space 10− Ωs + ΩA0 > ΩDM 2 0 Athena θ 12 10− Ly α Viable 10 14 − νs DM X rays 16 10−

Vacuum mixing angle Ωs < ΩDM 18 10−

20 10− 0.5 1 2 5 10 20 50 100 200 500 Sterile neutrino mass mν [keV]

2 Figure 5. Region of the ms-θ0 plane where sterile neutrino dark matter can be produced via dark photon-induced resonant oscillations, corresponding to the unshaded white area. In the areas shaded in grey, the mechanism either underproduces or overproduces dark matter as explained in the text. The blue line depicts the line where non-resonant oscillations yield the observed dark matter abundance in the absence of the Lµ Lτ gauge interactions. We also show constraints coming from X-ray searches (in red) and prospects for the Athena telescope− (purple dashed line), as well as the area excluded by phase-space considerations in dwarf galaxies (in orange) and the one disfavored by Lyman-α observations (green dots).

Fig. 5 clearly showcases that the new mechanism of We will present an exhaustive study in the future [94]; 0 A -induced resonant production of νs opens the sterile here we make a preliminary investigation to show that neutrino DM parameter space dramatically. This is par- there is no conflict with the allowed region in Fig. 5. ticularly relevant because previously untested regions of To study the effect of the new terms in Eq. (5) in ter- it will soon be probed as X-ray searches and structure restrial neutrino oscillations, we need the present value 0 formation observables continue to grow in sensitivity. Fi- A of the amplitude of the vector in the solar neighbor- nally, we note that the mechanism can also accommodate hood, which is related to its initial value by a sterile neutrino DM explanation of the 3.5 keV X-ray −10 3/4 0 0 line [91, 92], marked by the yellow star in Fig. 5. 0 0 −13 10 eV g A0 g A 1.3 10 eV (31) ' × m 0 1 keV  A    due to the redshifting of the energy density of A0 before V. ACTIVE NEUTRINO OSCILLATIONS structure formation, and the subsequent concentration av 5 of DM in the galaxy, ρDM/ρDM 3 10 . The value of 0 0 ∼ × The framework under study can also be tested through g A is therefore much smaller than the typical neutrino its potential impact on laboratory measurements of neu- energy in oscillation experiments, so that the (g0A0)2, trino oscillations, particularly at long baseline experi- term in Eq. (5) can be neglected, and the Hamiltonian ments sensitive to solar (νe νµ) or atmospheric (νµ takes the form → → ντ ) conversions [93]. This can happen if the new con- 2 2 1 mµ mµτ 0 0 1 0 tributions to the neutrino matter potential δVa from the Hµτ = p + + g A~ pˆ − , (32) 0 2p m2 m2 · 0 1 oscillating A field are still significant at the present time. µτ τ !   9

10 2 m = 10 eV, g A = 500 keV 0 A0 − 0 0 small momenta p . ∆m /4mA , a distortion of the shape 1.0 of the survival probability curve becomes evident, while

) Vacuum oscillation µ

ν for large p these disappear, except for an increase in the 0.8 → oscillation length relative to the standard vacuum oscil-

µ p > 1 GeV ν ( p = 0.01 GeV lation result, which has the same effect as a decrease in P 2 0.6 ∆m . This effect can be understood by using second-order 0.4 time-dependent perturbation theory to approximately solve Eq. (33). One finds an imaginary secular term grow- ing linearly in time that can be resummed to produce a 0.2 survival probability

µ shift δH in the effective Hamiltonian, appearing in the ν p = 0.1 GeV evolution operator exp(i(H0 + δH)t). The effect of δH 0.0 0 1 2 3 4 5 6 is precisely the same as a multiplicative correction to the ∆m2 2 Dimensionless time tˆ = 23 t vacuum mixing parameter ∆m , 4p 0 0 2 2 2 g A sin 2θ Figure 6. Muon neutrino survival probability in atmospheric ∆m ∆m 1 , (35) → − 2m 0 neutrino oscillations as a function of the dimensionless time  A  ! parameter tˆ defined in eq. (34), modified by interactions with an oscillating dark photon background. For this illustrative written here for a general mixing angle θ. By demanding 0 0 −10 2 example, we fix g A0 = 500 keV and mA0 = 10 eV. Curves that the fractional change in ∆m does not exceed the are shown for neutrino momenta p > 1 GeV (green), p = experimental error (which may be too stringent since for 0.1 GeV (red), and p = 0.01 GeV (orange), as well as the high p this effect is degenerate with the vacuum neutrino standard prediction when g0A0 = 0 (blue) as a reference. squared mass difference), and using Eq. (31), one can derive the limit

m 0 7/3 g0A0 55 keV A . (36) where the vector is oscillating at the present time as 0 . 10−10 eV ~0 ~0   A (t) = A cos(mA0 t), andp ˆ is a unit vector. This sys- tem can be| solved| numerically and compared with stan- Although this bound is strong for small dark photon 0 0 masses, it degrades steeply with increasing mA0 and dard oscillations, leading to constraints on g A and mA0 0 becomes subleading to the constraint in Eq. (36) for to avoid too big a deviation. 0 −8 1/3 m 4 10 (ΩA0 /ΩDM) eV. As the production It is convenient to rescale to a dimensionless time vari- A & × able tˆ = (∆m2/4p)t for each momentum p, where ∆m2 of the νs relic density can always be mediated by dark 2 −3 2 photons heavier than this threshold, this constraint does is the atmospheric (∆m23 2.5 10 eV ) or solar 2 −5 2 ' × not further restrict any of the parameter space shown in (∆m12 7.4 10 eV ) mass splitting. With this, the Schr¨odingerequation' × takes the form Fig. 5. For lighter dark photons, neutrino oscillations constitute a promising way to indirectly test the sterile neutrino DM paradigm presented in this paper. Since d ψνµ (cos 2θ + h) sin 2θ ψνµ i = − , the experimental error in ∆m2 is relatively smaller than dtˆ ψν ∼ sin 2θ cos 2θ + h ψν 23  τ     τ  2 (33) that for ∆m12, we find that solar neutrino oscillations with give a somewhat less stringent limit than (36). 0 0 p g A cos α 4p m 0 h(t) = 4 | | cos A tˆ ∆m2 ∆m2 VI. SUMMARY AND CONCLUSIONS   h0 cos ωt,ˆ (34) ≡ In this work, we have shown how the existence of a ~0 where α is the angle between ~p and A . In the tˆvariable, cosmological dark photon condensate can catalyze the all momentum modes have the same oscillation frequency production of keV sterile neutrino dark matter in the in the absence of the new physics effects from h(t). early universe. If active neutrinos couple to the dark As an example, we consider the case of atmospheric photon, which is assumed to be very light and weakly oscillations and approximate the vacuum νµ-ντ mixing coupled, their dispersion relation is modified in such a angle as θ = 45◦ in Eq. (33) for simplicity. Here we also way that resonant conversions into a sterile flavor can set α = 0; a more detailed study will average over the take place. For concreteness, we have implemented the directions of the neutrinos relative to A~0, slightly weak- coupling between the dark photon and the neutrinos by ening the ensuing constraints. gauging the anomaly-free combination Lµ Lτ , although Results for a set of parameters for which the new other possibilities may be viable. − physics effect would be appreciable in current atmo- Our main finding is that sterile neutrinos with masses spheric oscillation experiments are shown in Fig. 6. For in the few to 500 keV range and θ2 10−18 ∼ 0 ∼ − 10

−10 10 mixings, produced through the resonant oscilla- relation is parametrized as ω = k bL with tions enabled by the dark photon, may constitute all − −2 ¯ ¯ the dark matter of the universe. The mass and interac- cw B(0, mZ ) + 2B(mf , mW ), broken phase bL = παw 2 ¯ tion strength of the dark photon can span a large range, (2 + 1/cw)B(0, 0). sym. phase −12 0 −10  10 eV m 0 1 eV, and g 10 . (A1) . A . . ¯ The new resonant mechanism is much more effective Here, cw = cos θW and the function B depends on the than its nonresonant counterpart, commonly known as neutrino energy ω and momentum p, as well as T , mZ Dodelson-Widrow production [1]. This can be seen in and mf , the mass of the charged lepton µ or τ. Treating Fig. 4, where the nonresonant yield corresponds to the the thermal contribution as a perturbation, we can set low plateau of the curves at the left of the figure, and in ω = p, the unperturbed on-shell relation. This simplifies ∼ ¯ the form of B(0, mZ ) to Fig. 5, where the favored region for νs DM production extends to much lower mixings than the blue line repre- 1 dk m2 k pk2 ¯ Z + senting the DW prediction. Compared to another pop- B(0, mZ ) = L (k) 4 nb(z) p2 8π2 2  2 −  ular mechanism for resonant production of sterile neu- Z "  z z  trino dark matter [14], our proposal has the appeal of 2 mZ + not requiring the existence of a large primordial lepton + L (k) 4pk nf (k) , (A2) 2 1 − asymmetry.   # The small dark photon mass can in principle be gener- where nf and nb are the Fermi-Dirac or Bose-Einstein ated by either the Higgs or the St¨uckelberg mechanism. distribution functions,  = k2 + m2 , and However, if we assume that the dark photon condensate z Z 2 2 was produced gravitationally during high-scale inflation, [m + 2p(k +pz)][m + 2p(k z)] L+(k) = ln Z Z − , the St¨uckelberg option is disfavored by the Swampland 2 [m2 2p(k +  )][m2 2p(k  )]  Z z Z z  Distance Conjecture [76]. 2 − − − + mZ + 4pk We have also discussed how the existence of a Lµ Lτ L (p) = ln . (A3) − 1 m2 4pk dark photon condensate in the universe can impact the  Z −  oscillation pattern of µ and τ neutrinos. This may have This way, B becomes a function of p only, as we have interesting implications for long baseline or atmospheric assumed in the text. A smooth crossover at the EWPT neutrino oscillation experiments. We have only presented 2 2 implies that mZ,W (T ) ∼= mZ,W (1 T /TEW ), which fur- a preliminary study of this effect here, and we will analyze ther ensures a smooth dependence− on T near the EWPT. it in more detail in an upcoming publication [94]. At high temperatures, we can further simplify bL by Finally, it is worth noting that for some combinations taking mf ∼= 0 in eq. (A1), which allows us to write of parameters the dark photon itself can make a signifi- cant contribution to the dark matter abundance. In this −mf /T ¯ −2 ¯ bL ∼= παw 2e B(0, mW (T )) + cw B(0, mZ (T )) case, the DM would be an admixture of dark photons and  (A4) sterile neutrinos. This scenario could open additional pa- to describe either phase, with the factor e−mf /T sup- rameter space by relaxing the constraints on the sterile pressing the charged-current contribution for T < m . neutrino mass and mixing arising from X-ray searches f and structure formation observables, and perhaps even help reconcile the DM explanation of the 3.5 keV X- T = 100 GeV T = 99.99 GeV ray line with recent observations that seem to challenge T = 99.9 GeV 2 T = 99 GeV it [95]. We leave this interesting possibility for future T = 95 GeV T = 90 GeV investigation. T = 80 GeV 0 T = 70 GeV T = 50 GeV T = 30 GeV /T | T = 10 GeV L Acknowledgment. We thank M. Escudero, -2 T = √10 GeV T = 1 GeV S. Hannestad, K. Kainulainen, and M. Laine for help- | b T = 0.1 GeV ful correspondence, and J. Jaeckel for thoughtful com- 10 -4 ments on the manuscript. This work was supported log by NSERC (Natural Sciences and Engineering Research -6 Council, Canada). GA is supported by the McGill Space Institute through a McGill Trottier Chair Astrophysics -8 Postdoctoral Fellowship. -4 -2 0 2 4 log10 k / T Appendix A: Neutrino self-energy Figure 7. bL versus µ-neutrino momentum for a series of temperatures| as| shown in the legend, taking the EWPT tem- We can deduce the thermal contribution to the µ or perature to be TEW = 100 GeV. The dashed curves are the τ neutrino self-energy from eqs. (4.48,4.52) of ref. [96]. results of a fit as discussed in the text. In the relativistic limit, the correction to the dispersion 11

The neutrino thermal self-energy is thus obtained simply and eventually disappears. as VW = bL. At temperatures T TEW, the inte- gral (A2) can− be performed analytically, leading to We find that the curves have a nearly universal shape, such that any one of them can be obtained from another by shifting it in logarithmic space. This saves computa- 14π 2 −mf /T 2 4 VW (T TEW) GF sw 2e + cw p T . tional effort, since then it is possible to obtain bL /T as  ' −45α | |   (A5) a function of p/T at any T , from its functional form at a

In Fig. 7 we display bL /T versus p/T for a series of fixed reference temperature, in terms of the two functions | | temperatures. In agreement with eq. (2), bL is negative of T that describe how the curve is shifted horizontally at low p and T < TEW = 100 GeV, while it is everywhere and vertically in the plane of log bL /T and log p/T . 10 | | 10 positive for T > TEW . For T < TEW , bL changes sign These fits to the actual values are shown as dotted curves at some critical value of p, which decreases as T TEW in Fig. 7. →

[1] S. Dodelson and L. M. Widrow, “Sterile-neutrinos as BOSS and XQ-100,” JCAP 06 (2017) 047, dark matter,” Phys. Rev. Lett. 72 (1994) 17–20, arXiv:1702.03314 [astro-ph.CO]. arXiv:hep-ph/9303287. [13] A. Garzilli, A. Magalich, T. Theuns, C. S. Frenk, [2] K. N. Abazajian, “Sterile neutrinos in cosmology,” C. Weniger, O. Ruchayskiy, and A. Boyarsky, “The Phys. Rept. 711-712 (2017) 1–28, arXiv:1705.01837 Lyman-α forest as a diagnostic of the nature of the dark [hep-ph]. matter,” Mon. Not. Roy. Astron. Soc. 489 no. 3, (2019) [3] A. Boyarsky, M. Drewes, T. Lasserre, S. Mertens, and 3456–3471, arXiv:1809.06585 [astro-ph.CO]. O. Ruchayskiy, “Sterile neutrino Dark Matter,” Prog. [14] X.-D. Shi and G. M. Fuller, “A New dark matter Part. Nucl. Phys. 104 (2019) 1–45, arXiv:1807.07938 candidate: Nonthermal sterile neutrinos,” Phys. Rev. [hep-ph]. Lett. 82 (1999) 2832–2835, arXiv:astro-ph/9810076. [4] S. Horiuchi, P. J. Humphrey, J. Onorbe, K. N. [15] A. D. Dolgov, S. H. Hansen, S. Pastor, S. T. Petcov, Abazajian, M. Kaplinghat, and S. Garrison-Kimmel, G. G. Raffelt, and D. V. Semikoz, “Cosmological “Sterile neutrino dark matter bounds from galaxies of bounds on neutrino degeneracy improved by flavor the Local Group,” Phys. Rev. D 89 no. 2, (2014) oscillations,” Nucl. Phys. B 632 (2002) 363–382, 025017, arXiv:1311.0282 [astro-ph.CO]. arXiv:hep-ph/0201287. [5] M. R. Lovell, S. Bose, A. Boyarsky, S. Cole, C. S. [16] P. D. Serpico and G. G. Raffelt, “Lepton asymmetry Frenk, V. Gonzalez-Perez, R. Kennedy, O. Ruchayskiy, and primordial nucleosynthesis in the era of precision and A. Smith, “Satellite galaxies in semi-analytic cosmology,” Phys. Rev. D 71 (2005) 127301, models of galaxy formation with sterile neutrino dark arXiv:astro-ph/0506162. matter,” Mon. Not. Roy. Astron. Soc. 461 no. 1, (2016) [17] M. Laine and M. Shaposhnikov, “Sterile neutrino dark 60–72, arXiv:1511.04078 [astro-ph.CO]. matter as a consequence of nuMSM-induced lepton [6] A. Schneider, “Astrophysical constraints on resonantly asymmetry,” JCAP 06 (2008) 031, arXiv:0804.4543 produced sterile neutrino dark matter,” JCAP 04 [hep-ph]. (2016) 059, arXiv:1601.07553 [astro-ph.CO]. [18] A. Boyarsky, O. Ruchayskiy, and M. Shaposhnikov, [7] J. F. Cherry and S. Horiuchi, “Closing in on Resonantly “The Role of sterile neutrinos in cosmology and Produced Sterile Neutrino Dark Matter,” Phys. Rev. D astrophysics,” Ann. Rev. Nucl. Part. Sci. 59 (2009) 95 no. 8, (2017) 083015, arXiv:1701.07874 [hep-ph]. 191–214, arXiv:0901.0011 [hep-ph]. [8] A. Boyarsky, J. Lesgourgues, O. Ruchayskiy, and [19] L. Canetti, M. Drewes, T. Frossard, and M. Viel, “Lyman-alpha constraints on warm and on M. Shaposhnikov, “Dark Matter, Baryogenesis and warm-plus-cold dark matter models,” JCAP 05 (2009) Neutrino Oscillations from Right Handed Neutrinos,” 012, arXiv:0812.0010 [astro-ph]. Phys. Rev. D 87 (2013) 093006, arXiv:1208.4607 [9] M. Viel, G. D. Becker, J. S. Bolton, and M. G. [hep-ph]. Haehnelt, “Warm dark matter as a solution to the small [20] T. Venumadhav, F.-Y. Cyr-Racine, K. N. Abazajian, scale crisis: New constraints from high redshift and C. M. Hirata, “Sterile neutrino dark matter: Weak Lyman-α forest data,” Phys. Rev. D 88 (2013) 043502, interactions in the strong coupling epoch,” Phys. Rev. arXiv:1306.2314 [astro-ph.CO]. D 94 no. 4, (2016) 043515, arXiv:1507.06655 [10] A. Garzilli, A. Boyarsky, and O. Ruchayskiy, “Cutoff in [astro-ph.CO]. the Lyman alpha forest power spectrum: warm IGM [21] K. Perez, K. C. Y. Ng, J. F. Beacom, C. Hersh, or warm dark\ matter?,” Phys. Lett. B 773 (2017) S. Horiuchi, and R. Krivonos, “Almost closing the 258–264, arXiv:1510.07006 [astro-ph.CO]. νMSM sterile neutrino dark matter window with [11] J. Baur, N. Palanque-Delabrouille, C. Y`eche, NuSTAR,” Phys. Rev. D 95 no. 12, (2017) 123002, C. Magneville, and M. Viel, “Lyman-alpha Forests cool arXiv:1609.00667 [astro-ph.HE]. Warm Dark Matter,” JCAP 08 (2016) 012, [22] K. Abazajian, G. M. Fuller, and M. Patel, “Sterile arXiv:1512.01981 [astro-ph.CO]. neutrino hot, warm, and cold dark matter,” Phys. Rev. [12] C. Y`eche, N. Palanque-Delabrouille, J. Baur, and D 64 (2001) 023501, arXiv:astro-ph/0101524. H. du Mas des Bourboux, “Constraints on neutrino [23] F. L. Bezrukov and M. Shaposhnikov, “Searching for masses from Lyman-alpha forest power spectrum with dark matter sterile neutrino in laboratory,” Phys. Rev. 12

D 75 (2007) 053005, arXiv:hep-ph/0611352. arXiv:1810.07196 [hep-ph]. [24] A. Kusenko, “Sterile neutrinos, dark matter, and the [43] M. Bastero-Gil, J. Santiago, L. Ubaldi, and pulsar velocities in models with a Higgs singlet,” Phys. R. Vega-Morales, “Vector dark matter production at Rev. Lett. 97 (2006) 241301, arXiv:hep-ph/0609081. the end of inflation,” JCAP 04 (2019) 015, [25] K. Petraki and A. Kusenko, “Dark-matter sterile arXiv:1810.07208 [hep-ph]. neutrinos in models with a gauge singlet in the Higgs [44] A. J. Long and L.-T. Wang, “Dark Photon Dark Matter sector,” Phys. Rev. D 77 (2008) 065014, from a Network of Cosmic Strings,” Phys. Rev. D 99 arXiv:0711.4646 [hep-ph]. no. 6, (2019) 063529, arXiv:1901.03312 [hep-ph]. [26] M. Drewes et al., “A White Paper on keV Sterile [45] M. Bastero-Gil, J. Santiago, L. Ubaldi, and Neutrino Dark Matter,” JCAP 01 (2017) 025, R. Vega-Morales, “Dark photon dark matter from a arXiv:1602.04816 [hep-ph]. rolling inflaton,” arXiv:2103.12145 [hep-ph]. [27] A. Kusenko, F. Takahashi, and T. T. Yanagida, “Dark [46] D. Notzold and G. Raffelt, “Neutrino Dispersion at Matter from Split Seesaw,” Phys. Lett. B 693 (2010) Finite Temperature and Density,” Nucl. Phys. B 307 144–148, arXiv:1006.1731 [hep-ph]. (1988) 924–936. [28] B. Shuve and I. Yavin, “Dark matter progenitor: Light [47] M. B. Wise and Y. Zhang, “Lepton Flavorful Fifth vector boson decay into sterile neutrinos,” Phys. Rev. D Force and Depth-dependent Neutrino Matter 89 no. 11, (2014) 113004, arXiv:1403.2727 [hep-ph]. Interactions,” JHEP 06 (2018) 053, arXiv:1803.00591 [29] A. De Gouvˆea,M. Sen, W. Tangarife, and Y. Zhang, [hep-ph]. “Dodelson-Widrow Mechanism in the Presence of [48] J. A. Dror, R. Laha, and T. Opferkuch, “Probing Self-Interacting Neutrinos,” Phys. Rev. Lett. 124 no. 8, muonic forces with neutron star binaries,” Phys. Rev. D (2020) 081802, arXiv:1910.04901 [hep-ph]. 102 no. 2, (2020) 023005, arXiv:1909.12845 [hep-ph]. [30] J. Preskill, M. B. Wise, and F. Wilczek, “Cosmology of [49] E. Ma, D. P. Roy, and S. Roy, “Gauged Lµ Lτ with the Invisible Axion,” Phys. Lett. B 120 (1983) 127–132. large muon anomalous magnetic moment and− the [31] L. Abbott and P. Sikivie, “A Cosmological Bound on bimaximal mixing of neutrinos,” Phys. Lett. B 525 the Invisible Axion,” Phys. Lett. B 120 (1983) 133–136. (2002) 101–106, arXiv:hep-ph/0110146. [32] M. Dine and W. Fischler, “The Not So Harmless [50] S. Choubey and W. Rodejohann, “A Flavor symmetry Axion,” Phys. Lett. B 120 (1983) 137–141. for quasi-degenerate neutrinos: Lµ Lτ ,” Eur. Phys. J. [33] L. Hui, J. P. Ostriker, S. Tremaine, and E. Witten, C 40 (2005) 259–268, arXiv:hep-ph/0411190− . “Ultralight scalars as cosmological dark matter,” Phys. [51] J. Heeck and W. Rodejohann, “Gauged Lµ Lτ Rev. D 95 no. 4, (2017) 043541, arXiv:1610.08297 Symmetry at the Electroweak Scale,” Phys.− Rev. D 84 [astro-ph.CO]. (2011) 075007, arXiv:1107.5238 [hep-ph]. [34] A. E. Nelson and J. Scholtz, “Dark Light, Dark Matter [52] K. Asai, K. Hamaguchi, and N. Nagata, “Predictions and the Misalignment Mechanism,” Phys. Rev. D 84 for the neutrino parameters in the minimal gauged

(2011) 103501, arXiv:1105.2812 [hep-ph]. U(1)Lµ−Lτ model,” Eur. Phys. J. C 77 no. 11, (2017) [35] P. Arias, D. Cadamuro, M. Goodsell, J. Jaeckel, 763, arXiv:1705.00419 [hep-ph]. J. Redondo, and A. Ringwald, “WISPy Cold Dark [53] T. Araki, K. Asai, J. Sato, and T. Shimomura, “Low

Matter,” JCAP 06 (2012) 013, arXiv:1201.5902 scale seesaw models for low scale U(1)Lµ−Lτ [hep-ph]. symmetry,” Phys. Rev. D 100 no. 9, (2019) 095012, [36] G. Alonso-Alvarez,´ T. Hugle, and J. Jaeckel, arXiv:1909.08827 [hep-ph]. “Misalignment & Co.: (Pseudo-)scalar and vector dark [54] A. Kamada and H.-B. Yu, “Coherent Propagation of matter with curvature couplings,” JCAP 02 (2020) 014, PeV Neutrinos and the Dip in the Neutrino Spectrum arXiv:1905.09836 [hep-ph]. at IceCube,” Phys. Rev. D 92 no. 11, (2015) 113004, [37] P. W. Graham, J. Mardon, and S. Rajendran, “Vector arXiv:1504.00711 [hep-ph]. Dark Matter from Inflationary Fluctuations,” Phys. [55] G.-y. Huang, T. Ohlsson, and S. Zhou, “Observational Rev. D 93 no. 10, (2016) 103520, arXiv:1504.02102 Constraints on Secret Neutrino Interactions from Big [hep-ph]. Bang Nucleosynthesis,” Phys. Rev. D 97 no. 7, (2018) [38] Y. Ema, K. Nakayama, and Y. Tang, “Production of 075009, arXiv:1712.04792 [hep-ph]. Purely Gravitational Dark Matter: The Case of [56] M. Escudero, D. Hooper, G. Krnjaic, and M. Pierre, Fermion and Vector Boson,” JHEP 07 (2019) 060, “Cosmology with A Very Light Lµ Lτ Gauge Boson,” − arXiv:1903.10973 [hep-ph]. JHEP 03 (2019) 071, arXiv:1901.02010 [hep-ph]. [39] A. Ahmed, B. Grzadkowski, and A. Socha, [57] J. A. Dror, “Discovering leptonic forces using “Gravitational production of vector dark matter,” nonconserved currents,” Phys. Rev. D 101 no. 9, (2020) arXiv:2005.01766 [hep-ph]. 095013, arXiv:2004.04750 [hep-ph]. [40] E. W. Kolb and A. J. Long, “Completely Dark Photons [58] Y. Farzan, “Bounds on the coupling of the Majoron to from Gravitational Particle Production During light neutrinos from supernova cooling,” Phys. Rev. D Inflation,” arXiv:2009.03828 [astro-ph.CO]. 67 (2003) 073015, arXiv:hep-ph/0211375. [41] P. Agrawal, N. Kitajima, M. Reece, T. Sekiguchi, and [59] D. Croon, G. Elor, R. K. Leane, and S. D. McDermott, F. Takahashi, “Relic Abundance of Dark Photon Dark “Supernova Muons: New Constraints on Z’ Bosons, Matter,” Phys. Lett. B 801 (2020) 135136, Axions, and ALPs,” arXiv:2006.13942 [hep-ph]. arXiv:1810.07188 [hep-ph]. [60] S. Hannestad, “Structure formation with strongly [42] R. T. Co, A. Pierce, Z. Zhang, and Y. Zhao, “Dark interacting neutrinos - Implications for the cosmological Photon Dark Matter Produced by Axion Oscillations,” neutrino mass bound,” JCAP 02 (2005) 011, Phys. Rev. D 99 no. 7, (2019) 075002, arXiv:astro-ph/0411475. [61] S. Hannestad and G. Raffelt, “Constraining invisible 13

neutrino decays with the cosmic microwave [79] K. Enqvist, K. Kainulainen, and J. Maalampi, background,” Phys. Rev. D 72 (2005) 103514, “Refraction and Oscillations of Neutrinos in the Early arXiv:hep-ph/0509278. Universe,” Nucl. Phys. B 349 (1991) 754–790. [62] M. Escudero and M. Fairbairn, “Cosmological [80] G. Raffelt, G. Sigl, and L. Stodolsky, “NonAbelian Constraints on Invisible Neutrino Decays Revisited,” Boltzmann equation for mixing and decoherence,” Phys. Phys. Rev. D 100 no. 10, (2019) 103531, Rev. Lett. 70 (1993) 2363–2366, arXiv:1907.05425 [hep-ph]. arXiv:hep-ph/9209276. [Erratum: Phys.Rev.Lett. 98, [63] M. Escudero, J. Lopez-Pavon, N. Rius, and S. Sandner, 069902 (2007)]. “Relaxing Cosmological Neutrino Mass Bounds with [81] G. Sigl and G. Raffelt, “General kinetic description of Unstable Neutrinos,” JHEP 12 (2020) 119, relativistic mixed neutrinos,” Nucl. Phys. B 406 (1993) arXiv:2007.04994 [hep-ph]. 423–451. [64] G. Barenboim, J. Z. Chen, S. Hannestad, I. M. [82] K. Kainulainen, “Light Singlet Neutrinos and the Oldengott, T. Tram, and Y. Y. Y. Wong, “Invisible Primordial Nucleosynthesis,” Phys. Lett. B 244 (1990) neutrino decay in precision cosmology,” 191–195. arXiv:2011.01502 [astro-ph.CO]. [83] K. S. Babu and I. Z. Rothstein, “Relaxing [65] H. An, M. Pospelov, and J. Pradler, “New stellar nucleosynthesis bounds on sterile-neutrinos,” Phys. constraints on dark photons,” Phys. Lett. B 725 (2013) Lett. B 275 (1992) 112–118. 190–195, arXiv:1302.3884 [hep-ph]. [84] J. M. Cline, “Constraints on almost Dirac neutrinos [66] H. An, M. Pospelov, J. Pradler, and A. Ritz, “Direct from neutrino - anti-neutrino oscillations,” Phys. Rev. Detection Constraints on Dark Photon Dark Matter,” Lett. 68 (1992) 3137–3140. Phys. Lett. B 747 (2015) 331–338, arXiv:1412.8378 [85] T. Asaka, M. Laine, and M. Shaposhnikov, “Lightest [hep-ph]. sterile neutrino abundance within the nuMSM,” JHEP [67] E. Hardy and R. Lasenby, “Stellar cooling bounds on 01 (2007) 091, arXiv:hep-ph/0612182. http://www. new light particles: plasma mixing effects,” JHEP 02 laine.itp.unibe.ch/neutrino-rate/index.html. (2017) 033, arXiv:1611.05852 [hep-ph]. [Erratum: JHEP 02, 028 (2015)]. [68] H. An, M. Pospelov, J. Pradler, and A. Ritz, “New [86] J. Ghiglieri and M. Laine, “Neutrino dynamics below limits on dark photons from solar emission and keV the electroweak crossover,” JCAP 07 (2016) 015, scale dark matter,” Phys. Rev. D 102 (2020) 115022, arXiv:1605.07720 [hep-ph]. http://www.laine.itp. arXiv:2006.13929 [hep-ph]. unibe.ch/production-midT/Gamma.dat. [69] M. Baryakhtar, R. Lasenby, and M. Teo, “Black Hole [87] M. Laine and M. Meyer, “Standard Model Superradiance Signatures of Ultralight Vectors,” Phys. thermodynamics across the electroweak crossover,” Rev. D 96 no. 3, (2017) 035019, arXiv:1704.05081 JCAP 07 (2015) 035, arXiv:1503.04935 [hep-ph]. [hep-ph]. [88] A. Neronov and D. Malyshev, “Toward a full test of the [70] V. Cardoso, P. Pani, and T.-T. Yu, “Superradiance in νMSM sterile neutrino dark matter model with rotating stars and pulsar-timing constraints on dark Athena,” Phys. Rev. D 93 no. 6, (2016) 063518, photons,” Phys. Rev. D 95 no. 12, (2017) 124056, arXiv:1509.02758 [astro-ph.HE]. arXiv:1704.06151 [gr-qc]. [89] J. Alvey, N. Sabti, V. Tiki, D. Blas, K. Bondarenko, [71] V. Cardoso, O. J. C. Dias, G. S. Hartnett, A. Boyarsky, M. Escudero, M. Fairbairn, M. Orkney, M. Middleton, P. Pani, and J. E. Santos, “Constraining and J. I. Read, “New constraints on the mass of the mass of dark photons and axion-like particles fermionic dark matter from dwarf spheroidal galaxies,” through black-hole superradiance,” JCAP 03 (2018) Mon. Not. Roy. Astron. Soc. 501 no. 1, (2021) 043, arXiv:1801.01420 [gr-qc]. 1188–1201, arXiv:2010.03572 [hep-ph]. [72] CMB-S4 Collaboration, K. N. Abazajian et al., [90] J. Baur, N. Palanque-Delabrouille, C. Yeche, “CMB-S4 Science Book, First Edition,” A. Boyarsky, O. Ruchayskiy, E. Armengaud, and arXiv:1610.02743 [astro-ph.CO]. J. Lesgourgues, “Constraints from Ly-α forests on [73] D. Blas and S. J. Witte, “Quenching Mechanisms of non-thermal dark matter including resonantly-produced Photon Superradiance,” Phys. Rev. D 102 (2020) sterile neutrinos,” JCAP 12 (2017) 013, 123018, arXiv:2009.10075 [hep-ph]. arXiv:1706.03118 [astro-ph.CO]. [74] E. Cannizzaro, A. Caputo, L. Sberna, and P. Pani, [91] E. Bulbul, M. Markevitch, A. Foster, R. K. Smith, “Plasma-photon interaction in curved spacetime I: M. Loewenstein, and S. W. Randall, “Detection of An formalism and quasibound states around nonspinning Unidentified Emission Line in the Stacked X-ray black holes,” arXiv:2012.05114 [gr-qc]. spectrum of Galaxy Clusters,” Astrophys. J. 789 (2014) [75] P. Amaro-Seoane et al., “Laser Interferometer Space 13, arXiv:1402.2301 [astro-ph.CO]. Antenna,” arXiv e-prints (Feb., 2017) arXiv:1702.00786, [92] A. Boyarsky, O. Ruchayskiy, D. Iakubovskyi, and arXiv:1702.00786 [astro-ph.IM]. J. Franse, “Unidentified Line in X-Ray Spectra of the [76] M. Reece, “Photon Masses in the Landscape and the Andromeda Galaxy and Perseus Galaxy Cluster,” Phys. Swampland,” JHEP 07 (2019) 181, arXiv:1808.09966 Rev. Lett. 113 (2014) 251301, arXiv:1402.4119 [hep-th]. [astro-ph.CO]. [77] H. Ooguri and C. Vafa, “On the Geometry of the String [93] Super-Kamiokande Collaboration, Y. Fukuda et al., Landscape and the Swampland,” Nucl. Phys. B 766 “Evidence for oscillation of atmospheric neutrinos,” (2007) 21–33, arXiv:hep-th/0605264. Phys. Rev. Lett. 81 (1998) 1562–1567, [78] N. Craig and I. Garcia Garcia, “Rescuing Massive arXiv:hep-ex/9807003. Photons from the Swampland,” JHEP 11 (2018) 067, [94] G. Alonso-Alvarez,´ K. Bleau, and J. M. Cline, arXiv:1810.05647 [hep-th]. 14

“Distortion of νµ-ντ oscillations by a dark photon (2020) 1465–1467, arXiv:1812.06976 [astro-ph.CO]. field,” in preparation (2021) . [96] C. Quimbay and S. Vargas-Castrillon, “Fermionic [95] C. Dessert, N. L. Rodd, and B. R. Safdi, “The dark dispersion relations in the standard model at finite matter interpretation of the 3.5-keV line is inconsistent temperature,” Nucl. Phys. B 451 (1995) 265–304, with blank-sky observations,” Science 367 no. 6485, arXiv:hep-ph/9504410.