Sterile neutrino dark matter catalyzed by a very light dark photon Gonzalo Alonso-Alvarez´ ∗ and James M. Cliney 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 y [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.