View metadata, citation and similar papers at core.ac.uk brought to you by CORE
provided by Caltech Authors - Main PHYSICAL REVIEW D 94, 043515 (2016) Sterile neutrino dark matter: Weak interactions in the strong coupling epoch
Tejaswi Venumadhav,1,5 Francis-Yan Cyr-Racine,1,2,6 Kevork N. Abazajian,3 and Christopher M. Hirata4 1California Institute of Technology, Mail Code 350-17, Pasadena, California 91125, USA 2Jet Propulsion Laboratory, California Institute of Technology, Pasadena, California 91109, USA 3Center for Cosmology, Department of Physics and Astronomy, University of California, Irvine, Irvine, California 92697, USA 4Center for Cosmology and Astroparticle Physics (CCAPP), The Ohio State University, 191 West Woodruff Lane, Columbus, Ohio 43210, USA 5School of Natural Sciences, Institute for Advanced Study, Einstein Drive, Princeton, New Jersy 08540, USA 6Department of Physics, Harvard University, Cambridge, Massachusetts 02138, USA (Received 14 August 2015; published 15 August 2016) We perform a detailed study of the weak interactions of standard model neutrinos with the primordial plasma and their effect on the resonant production of sterile neutrino dark matter. Motivated by issues in cosmological structure formation on small scales, and reported x-ray signals that could be due to sterile neutrino decay, we consider 7 keV-scale sterile neutrinos. Oscillation-driven production of such sterile neutrinos occurs at temperatures T ≳ 100 MeV, where we study two significant effects of weakly charged species in the primordial plasma: (1) the redistribution of an input lepton asymmetry; (2) the opacity for active neutrinos. We calculate the redistribution analytically above and below the quark-hadron transition, and match with lattice QCD calculations through the transition. We estimate opacities due to tree-level processes involving leptons and quarks above the quark-hadron transition, and the most important mesons below the transition. We report final sterile neutrino dark matter phase space densities that are significantly influenced by these effects, and yet relatively robust to remaining uncertainties in the nature of the quark- hadron transition. We also provide transfer functions for cosmological density fluctuations with cutoffs at k ≃ 10h Mpc−1, that are relevant to galactic structure formation.
DOI: 10.1103/PhysRevD.94.043515
I. INTRODUCTION plasma [10–14], or (ii) that sterile neutrinos are massive enough to form the inferred population of dark matter (DM) Deep in the radiation dominated epoch of the Universe, in the Universe (see e.g. Ref. [15]). Sterile neutrinos with the three neutrinos present in the standard model (SM) of masses in the keV range act as DM in the CMB era, but are particle physics [1] make up a significant population of relativistic in the BBN era, when they do not significantly relativistic species within the primeval cosmic plasma. We impact the expansion rate due to their negligible energy have strong evidence of their existence at these early density (compared to the Fermi-Dirac value). epochs from probes of the primordial Universe such as Early works in this direction studied right-handed sterile the cosmic microwave background (CMB) (probing tem- ≈ 0 1 − 100 ∼ 0 25 neutrinos with masses ms . keV, produced by perature T . eV) [2], and the synthesis of light the oscillation of left-handed SM neutrinos [16–20]. The elements during the epoch of big bang nucleosynthesis mixing angle between the SM and sterile neutrinos is fixed (BBN) which depends on the neutron-to-proton ratio set by the present day DM abundance. In the original ∼ 1 5 at Tdec . MeV [3], the temperature of weak neutrino Dodelson-Widrow scenario [17], sterile neutrinos are decoupling. Above this temperature, SM neutrinos interact produced with a momentum distribution reflecting that with species that carry weak charge, through which they of the active neutrino species, and thus constitute “warm” remain coupled to the primordial plasma [4]. DM [21–24]. However, small-scale structure formation There is a long history of speculation about additional [25–32] and x-ray observations [33–37] appear to be in neutrino species (see Ref. [5] for a recent review). Owing to significant conflict with the fiducial Dodelson-Widrow the precise measurement of the invisible decay width of the scenario, hence prompting the search for alternative sterile SM Z boson [1], any extra neutrino species must be “sterile” neutrino production mechanisms [18,38–55]. (i.e. electroweak singlets) [6]. Furthermore, precise mea- In this paper, we examine in detail the resonant production surements of the CMB [2,7] and of the primeval abundance of sterile neutrinos in the presence of a small primordial of light elements [8,9] strongly constrain the presence of lepton asymmetry. Originally proposed by Shi and Fuller extra relativistic species in the early Universe. These [38], this production mechanism makes use of a small lepton constraints indicate that (i) unlike SM neutrinos, light sterile asymmetry to modify the plasma’s interaction with SM neutrinos never fully thermalize with the rest of the cosmic neutrinos in such a manner as to resonantly produce sterile
2470-0010=2016=94(4)=043515(28) 043515-1 © 2016 American Physical Society TEJASWI VENUMADHAV et al. PHYSICAL REVIEW D 94, 043515 (2016) neutrinos at particular momenta [18,41,56]. This generically redistribution modifies the dynamics of the resonant sterile results in a “colder” DM distribution which improves neutrino production, resulting in a modified final PSD. consistency with models of cosmological structure forma- Secondly, we incorporate several new elements to the tion [57–65], while requiring a modest primordial lepton calculations of the neutrino opacity (i.e. the imaginary part asymmetry, which is relatively poorly constrained [66–70]. of the self-energy) at temperatures 10 MeV ≤ T ≤ 10 GeV. Sterile and active neutrino mixing, which is needed for Accurate neutrino opacities are needed since they basically the former’s production, also leads to their decay [71,72]. control the production rate of sterile neutrinos through cosmic For typical values of the sterile neutrino mass this predicts epochs. Early works on neutrino interactions in the early an x-ray flux from the DM distribution in the low redshift Universe [4,18,84] assumed that neutrinos largely scattered off Universe [18,33]. This has been the subject of much recent relativistic particles and thus scaled their cross sections with interest, due to hints of an excess flux at ∼3.5 keV in the center-of-mass (CM) energy. In addition, these calcula- stacked x-ray spectra of several galaxy clusters [73] and in tions also neglected the effects of particle statistics. Under Γ observations of M31, the Milky Way, and Perseus [74,75]. these two simplifying assumptions, the opacity ðEνα Þ for an There is currently an active debate on the existence, input neutrino of energy Eνα is of the form significance and interpretation of this excess [76–83].In Γ λ 2 4 the present work, we use this tentative signal as a ðEνα Þ¼ ðTÞGFT Eνα ; ð2Þ motivation to study in detail the physics of sterile neutrino λ production in the early Universe, but the machinery we where GF is the Fermi coupling constant, and ðTÞ develop is more generally applicable to the broader is a constant that depends on the number and type of available parameter space of the Shi-Fuller mechanism. relativistic species in the cosmic plasma. References [20,85] We present here an updated calculation of resonantly subsequently developed a framework to include particle produced sterile neutrinos and relax several simplifications masses, loop corrections, and particle statistics in the neutrino that had been adopted previously in the literature. opacity calculation. In the present work, we add previously Furthermore, we leverage recent advances in our under- neglected contributions to the opacity such as two- and three- standing of the quark-hadron transition in order to include a body fusion reactions, and also use chiral perturbation theory more realistic treatment of the strongly interacting sector. to compute the hadronic contribution to the opacity below the Our motivation is twofold: (a) improve the treatment of quark-hadron transition. We find both quantitative and quali- lepton asymmetry, which is a crucial beyond-SM ingredient tative modifications to the form of Eq. (2).Whereverwe − in the mechanism, and (b) provide realistic sterile neutrino present matrix elements, we use the þþþmetric signature. phase space densities (PSDs) and transfer functions for Thirdly, we fold the asymmetry redistribution and matter fluctuations, which are starting points for studying opacity calculations into the sterile neutrino production cosmological implications on small scales. Our improve- computation, and provide updated PSDs for the range of ments to the sterile neutrino production calculation can parameters relevant to the x-ray excess. As part of this broadly be classified in three categories. process, we carefully review and correct the numerical Firstly, we study how the cosmic plasma reprocesses a implementation of the sterile neutrino production used in primordial lepton asymmetry. For models that can explain Ref. [63]. Our sterile neutrino production code is publicly ‑ the above x-ray excess, the majority of sterile neutrinos are available at https://github.com/ntveem/sterile dm.We produced at temperatures above 100 MeV [56]. At these finally use the updated sterile neutrino PSDs in a standard temperatures, there is a significant population of either cosmological Boltzmann code [86] and provide new dark quarks or mesons, depending on whether the temperature is matter transfer functions. above or below the quark-hadron transition. Since these We organize the paper such that the beginning sections hadronic species are coupled to neutrinos and charged deal with SM physics, while the later ones apply their leptons through weak processes, the establishment of results to sterile neutrino DM. We introduce the production chemical equilibrium among the different constituents of mechanism in Sec. II. We then study the asymmetry the cosmic plasma will automatically redistribute the redistribution in Sec. III and active neutrino opacities in primordial asymmetry within the charged and neutral Sec. IV. Finally, we apply these results to sterile neutrino leptons, and lead to a charge asymmetry within the production in Sec. V, evaluate transfer functions for matter hadronic sector. An illustrative example is the reaction fluctuations in Sec. VI, and finish with a discussion of our assumptions and uncertainties in Sec. VII. We collect þ þ νμ þ μ ⇌π ; ð1Þ technical details into the appendixes. which can redistribute an initial neutrino asymmetry into II. OVERVIEW OF RESONANT STERILE asymmetries of the charged leptons and hadrons (keeping NEUTRINO PRODUCTION the net baryon number constant). At lower temperatures, the asymmetry is redistributed to a lesser degree between the In this section, we briefly review the resonant production leptonic flavors. As we discuss in the body of the paper, this of sterile neutrinos in the early Universe. We first present
043515-2 STERILE NEUTRINO DARK MATTER: WEAK … PHYSICAL REVIEW D 94, 043515 (2016) the specific scenario that we consider in this work, and then discuss the Boltzmann formalism used to compute the out- of-equilibrium production of sterile neutrinos. We finally discuss how the presence of the thermal bath and lepton asymmetry change the neutrino self-energy and govern the sterile neutrino production. We refer the reader to Refs. [18,38,56] for more details.
A. Assumptions In our current study, we focus on the following scenario. (1) We consider an extra sterile neutrino species, νs, that is massive compared to the SM neutrinos νe=μ=τ, which we take to be effectively massless. The propagation (light/heavy) and interaction (active/ sterile) eigenstates are related by a unitary trans- formation, the most general version of which is a FIG. 1. Sterile neutrino DM parameter space: shaded regions are 4 × 4 matrix. We assume that the sterile neutrino consistent with the x-ray signal at 1,2 and 3σ. The blue contours mixes with only one of the SM ones, which we take show the best determined parameters from the analysis of Ref. [73], to be the muon neutrino, i.e., who use the stacked x-ray spectrum of 73 clusters obtained with the MOS detectors on XMM Newton. Statistically consistent Ψ θ θ Ψ νμ cos sin 0 signals are found in their core-removed Perseus spectrum, and ¼ : ð3Þ M31 [74]. The lines show constraints at the 90% level from Ψν − sin θ cos θ Ψ s ms Chandra observations of M31 (H14) [87], stacked dwarf galaxies (M14) [88], and Suzaku observations of Perseus (T15) [89]. Stars The fields on the left- and right-hand sides are mark the models that we study in the body of the paper. interaction and mass (ms) eigenstates, respectively, and θ is the active-sterile mixing angle. The choice of a muon neutrino is arbitrary, and reflects the with constraints from Chandra observations of M31 [87], choice of previous work [18,63]. stacked dwarf galaxies [88], and Suzaku observations of (2) We assign a nonzero lepton asymmetry to the Perseus [89]. The stars show a range of mixing angles at a primordial plasma. In the general case, each SM specific value of ms, and mark models that we study in flavor has its own asymmetry, but we assume a Secs. V and VI. nonzero value only for the mu flavor (i.e. the one that For all these models, the bulk of the sterile neutrinos are mixes with the sterile neutrino): produced at temperatures well below the masses of the weak gauge bosons (∼80 GeV), but above weak decou- Δˆ Δˆ − ≡ Lˆ δ Lˆ ∼ 1 5 nνα þ nα α ¼ αμ μ; ð4Þ pling at T . MeV [56]. Active-active neutrino oscil- lations in the primordial plasma are suppressed at these ˆ where the dimensionless asymmetry ΔnA in species temperatures [90]; hence it is consistent to assign individual A is the temperature-scaled difference between asymmetries in Eq. (4) and neglect electron and tau the particle and antiparticle densities, Δnˆ A≡ neutrino mixing in Eq. (3). − 3 δ ðnA nAÞ=T , and αμ is the Kronecker delta. In general, entropy-scaled asymmetries are preferable, B. Boltzmann formalism since they are conserved through epochs of annihi- In its full generality, out-of-equilibrium sterile neutrino lation. However, the definition used in Eq. (4) production (via oscillations) is best described by the simplifies the process of comparison with lattice evolution of the two-state density matrix of the neutrinos QCD calculations in Sec. III. We fix by hand the mu in the active-sterile (interaction) basis [91–94]. lepton asymmetry at high temperatures to produce For the parameter range in Fig. 1, most sterile neutrinos Ω 2 0 1188 the canonical DM density, DMh ¼ . in the are produced above temperatures T ≳ 100 MeV. At these current epoch [2]. temperatures, the two-state system is collision dominated; In the rest of the paper, we use a hat to indicate i.e. the interaction contribution dominates the vacuum temperature scaled quantities. We choose to study the oscillations. In this regime, the evolution of the density parameter space of interest for resonantly produced sterile matrix separates out and yields a quasiclassical Boltzmann neutrino DM consistent with the recent x-ray signal. transport equation for the diagonal terms, which are the 2 2θ Figure 1 shows a section of the ms and sin plane with PSDs of the active and sterile components [95–97]. The contours for the unidentified lines of Refs. [73,74], along Boltzmann equation for the sterile neutrino PSD is
043515-3 TEJASWI VENUMADHAV et al. PHYSICAL REVIEW D 94, 043515 (2016) Z X 3 3 ∂ ∂ d pa d pi 4 4 fν ðp; tÞ − Hp fν ðp; tÞ¼ ð2πÞ δ ðp þ p þ − p − Þ ∂ s ∂ s 3 3 a i t p ð2πÞ 2Ea ð2πÞ 2Ei νxþaþ →iþ X 1 2 × hP ðνμ → ν ; p; tÞið1 − fν Þ jMj → ν f ð1 ∓ f Þð1 − fν Þ 2 m s s iþ aþ μþ i a μ X − ν → ν ; 1 − M 2 1 ∓ hPmð s μ p; tÞifνs ð fνμ Þ j jνμþaþ →iþ fa ð fiÞ : ð5Þ
We can write an analogous equation for the antineu- quantum-damped and collisionally driven sterile neutrino trinos. Here, the fðpÞ are PSDs for particles with three- production is [18] momentum p and energy E, and H is the Hubble expansion rate. The right-hand side sums over all reactions that ∂ ∂ fν ðp; tÞ − Hp fν ðp; tÞ Pconsume or produce a muon neutrino. The symbol ∂t s ∂p s jMj2 denotes the squared and spin-summed matrix Γ νμ ðpÞ element for the reaction, and the multiplicative factors of ≈ hP ðνμ ↔ ν ; p; tÞi½fν ðp; tÞ − fν ðp; tÞ ; 2 m s μ s ð1 ∓ fÞ implement Pauli blocking/Bose enhancement respectively. The factor of 1=2 accounts for the fact that ð8Þ only one (i.e. the muon neutrino) state in the two-state – with a related equation for antineutrinos. There are subtleties system interacts [91 93]. The Pm are active-sterile oscil- lation probabilities in matter, which depend on the vacuum with the effects of quantum damping in the case of resonance mixing angle θ, and are modified by interactions with the [98], but tests with the full density matrix formalism find that medium. The latter are parametrized by the neutrino self- the quasiclassical treatment is appropriate [56]. energy [84], and the “quantum damping” rate for active neutrinos. In terms of these quantities, the oscillation C. Asymmetry and thermal potentials probabilities are [96,97] We now expand on the origins of the asymmetry and thermal potentials appearing in Eq. (6). These potentials ν ↔ ν ; 1 2 Δ2 22θ hPmð μ s p; tÞi¼ð = Þ ðpÞsin encapsulate the self-energy of propagating active neutrinos × fΔ2ðpÞsin22θ þ D2ðpÞ due to interactions with the plasma. Under the conditions we are interested in, there are three contributions to the 2 −1 þ½ΔðpÞ cos 2θ − VL − VthðpÞ g : neutrino self energy: (a) an imaginary part proportional to ð6Þ the net neutrino opacity, (b) a real part due to finite weak gauge boson masses (Vth), and (c) a real part proportional to L We have introduced the symbol ΔðpÞ for the vacuum asymmetries in weakly interacting particles (V ). We Δ ≡ δ 2 2 follow the treatment in Ref. [84], and recast it in terms oscillation rate, ðpÞ mνμ;ν = p, and split the neutrino s of the quantities that we compute later. self-energy into the lepton asymmetry potential VL, and the th Figure 2 shows lowest-order contributions to active thermal potential V [the asymmetry contribution enters neutrinos’ self-energy. Thick red lines are thermal propa- with the opposite sign in the version of Eq. (6) for gators of weakly charged species in the background antineutrinos]. The quantity DðpÞ is the quantum damping plasma. There are two corrections—bubbles and tadpoles, rate, and equals half the net interaction rate of active neutrinos [the factor of half enters for the same reason as it does in Eq. (5)]. The net interaction rate for a muon neutrino is Z X 3 3 Γ d pa d pi νμ ðpÞ¼ 3 3 ð2πÞ 2Ea ð2πÞ 2Ei νxþaþ →iþ × ð2πÞ4δ4ðp þ p þ − p − Þ X a i M 2 1 ∓ × j jνμþaþ →iþ fa ð fiÞ ð7Þ
We simplify the phase space integrals in Eq. (5) by using FIG. 2. Lowest-order contributions to a propagating active detailed balance to equate the forward and backward neutrino’s self-energy. Red lines are thermal propagators. In − reaction rates. The resulting Boltzmann equation for (a), f is any species with weak charge. In (b), f ¼ να, α .
043515-4 STERILE NEUTRINO DARK MATTER: WEAK … PHYSICAL REVIEW D 94, 043515 (2016) θ shown in Figs. 2(a) and 2(b) respectively. The background In the above equations, W is the weak mixing angle, and fermion is a lepton of the same flavor in the former, and any MZ=W are the masses of the weak gauge bosons. The δ ρ ρ weakly charged species in the latter. symbol αβ is a Kronecker delta; the quantities α and να A massless active neutrino’s “dressed” propagator is are net energy densities of charged and neutral leptons, respectively; and Δn and Δn are densities of the baryon −1 − 1 − γ 2 B Q Gνα ðpνα Þ¼pνα bνα ðpνα Þuð 5Þ= ; ð9aÞ number B and electric charge Q, respectively. The standard model baryon number asymmetry is small compared to the ð0Þ ð1Þω ω − bνα ðpνα Þ¼bνα þ bνα να ; να ¼ pνα · u: ð9bÞ lepton asymmetry of interest [3]; hence it can be set to zero for the purposes of this calculation. ’ Here, pνα and u are the neutrino and plasma s four- According to the assumptions in the first part of this γμ section, the plasma starts out with a net lepton asymmetry momenta, v is shorthand for vμ, and bνα is the left- handed neutrino’s self-energy. Equation (9b) divides this in the mu flavor. As we showed in Sec. I, this asymmetry is self-energy into two contributions that affect the particle redistributed between muons and muon neutrinos. and antiparticle poles of Eq. (9a) differently. Figure 3 Moreover leptons of other flavors acquire asymmetries illustrates their association with asymmetry and thermal that respect Eq. (4), and the strong fluid acquires a net Δ potentials: electric charge density nQ to maintain overall neutrality. Equation (12a) shows how the asymmetry potential L ð0Þ depends on the redistributed asymmetries, which we study Vν ¼ bν ; ð10Þ α α in the ensuing section.
th ð1Þ Vν ðEν Þ¼bν Eν : ð11Þ α α α α III. REDISTRIBUTION OF AN INPUT ASYMMETRY The authors of Ref. [84] computed these terms by summing over all species in Fig. 2. Both kinds of diagrams contribute Weak processes couple leptonic and hadronic degrees of to the asymmetry potential, while only bubble diagrams freedom in the primordial plasma. In this section, we study 1 contribute to the thermal potential. We write the answer in this coupling’s effect on lepton asymmetries. We define terms of the leptons’ asymmetries, and the densities of the relevant susceptibilities in Sec. III A, and compute them strong fluid’s conserved quantities: over a range of temperatures in Sec. III B. pffiffiffi X 1 L 2 A. Definitions and parametrization Vν ¼ 2G δαβ − þ 2sin θ Δnβ− α F 2 W β∈fe;μ;τg Let us consider the primordial plasma at temperatures X 1 above the quark-hadron transition temperature, TQCD. The þ ð1 þ δαβÞΔnν − Δn following reactions couple leptons of different flavors, and β 2 B β∈fe;μ;τg the quark and lepton sectors: 1 − 2 2θ Δ − − þð sin WÞ nQ ; ð12aÞ να þ β ⇌νβ þ α ; ð13aÞ
pffiffiffi ν αþ⇌ 8 2 ρ ρ α þ a þ b; ð13bÞ th GF να α Vν ðEν Þ¼− þ Eν : ð12bÞ α α 3 2 2 α MZ MW where a and b are quarks with charges of þ2=3 and −1=3 respectively. Free quarks no longer exist at temperatures below TQCD, and the reactions in Eq. (13b) transition to ones involving mesons, like Eq. (1). In principle, we could study the effect of all these reactions on input asymmetries, but it is a daunting task, one that is further complicated by the quark-hadron transition. The following consideration of the relevant time
1During this preparation of this manuscript, we became aware of Ref. [99], which points out the relevance of this effect to sterile neutrino production, and estimates it under the simplifying Stefan-Boltzmann approximation for free quarks, along with a modified number of colors as a phenomenological correction for FIG. 3. Matter potentials for massless neutrinos in the plasma’s the effects of confinement. Wherever possible, we will try to rest frame: filled and unfilled circles are poles at finite and zero compare our results to those of Ref. [99], which we will refer to as temperature, respectively. See Sec. IV for the imaginary shift. “GL15” in the rest of the paper.
043515-5 TEJASWI VENUMADHAV et al. PHYSICAL REVIEW D 94, 043515 (2016)
μˆ − μˆ − − μˆ 0 α ∈ μ τ scales suggests a solution. In the radiation dominated να α Q ¼ ; fe; ; g: ð16Þ ≈ 2 105 −1 1=2 2 era, the Hubble rate is H × s g ðT=GeVÞ . μˆ At temperatures above the quark-hadron transition, the Here Q is the chemical potential for adding a unit of Γ ≃ 2 5≈ rates of reactions in (13) are ðTÞ GFT electric charge. The conserved quantities’ densities are 2 × 1014 s−1ðT=GeVÞ5, while the relevant rates below related to their chemical potentials by their susceptibilities: the transition are those of pion decays. The most important ! ! ! þ þ channel for the latter is the muonic decay, π → μ þ νμ, Δˆ χˆQ χˆQB μˆ nQ 2 11 Q which is faster than the Hubble rate (Γπþ→μþ ν ¼ : ð17Þ þ μ ¼ Δˆ BQ B μˆ nB χˆ χˆ B 3.8 × 107 s−1). Thus, a significant number of the reactions 11 2 in Eq. (13) are faster than the Hubble rate.2 Equation (17), along with net charge and baryon number This has two primary consequences. Firstly, high reaction conservation, yields the constraint equations rates enforce kinetic equilibrium; i.e. all active species’ PSDs approach the Fermi-Dirac or Bose-Einstein forms. Secondly, Δˆ χˆBQμˆ χˆBμˆ ≈ 0 forward and backward reactions are in chemical equilibrium, nB ¼ 11 Q þ 2 B ; ð18Þ one effect of which is to equate the chemical potentials for X Q QB 0 χˆ μˆ χˆ μˆ − Δˆ − both sides (the Saha equation). However, it has another ¼ 2 Q þ 11 B nα implication—the plasma’s complicated internal dynamics α∈fe;μ;τg X can be abstracted into a few parameters or susceptibilities that Q QB ¼ χˆ μˆ þ χˆ μˆ − χˆα− μˆ α− : ð19Þ completely specify its response to small external “forces,” or 2 Q 11 B α∈ e;μ;τ in this case, input asymmetries. All that remains is to compute f g the susceptibilities relevant to our problem. Equations (15), (16), (18) and (19) are eight linear We now define a few useful quantities and notation. equations for eight unknowns. The resulting asymmetries Given any conserved quantity F, the symbol μ denotes its F (obtained via their chemical potentials) are the “redistrib- chemical potential. The asymmetry Δnˆ A, in a particle A,isa uted” input lepton asymmetries Lα. function of its chemical potential μˆ ≡ μ =T. The quan- A A We symbolically represent the solutions as tities Δnˆ A and μˆ A are small, and in the linearized limit, related by X ∂μˆ μˆ A Lˆ A ¼ ˆ α; ð20Þ ∂Lα Δnˆ A ¼ χˆAμˆ A; ð14Þ α∈fe;μ;τg 2 ∂μˆ ∂Lˆ where χˆA ≡ χA=T is the number-density susceptibility. where the coefficients ð A= αÞ depend on the suscep- The lepton asymmetries in the three flavors are tibilities of both the leptons and the strong fluid. We also express the redistributed asymmetries as ˆ Lα ¼ Δnˆ α− þ Δnˆ ν α X ∂Δˆ X ∂μˆ nA ˆ A ˆ ¼ χˆα− μˆ α− þ χˆν μˆ ν ; Δnˆ ¼ Lα ¼ χˆ Lα: ð21Þ α α A ∂Lˆ A ∂Lˆ α∈fe;μ;τg α α∈fe;μ;τg α α ∈ fe; μ; τg: ð15Þ At the temperatures of interest, the lepton susceptibilities The strong fluid is described by the densities of its are essentially given by the free particle, or Stefan- conserved quantities: the charge and baryon-number den- Δˆ Δˆ 3 Boltzmann, formula: sities nQ and nB, respectively. The chemical equilib- Z qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi rium of the reactions in Eq. (13) implies g ∞ χˆ ˆ − A ˆ ˆ 2 ˆ 0 ˆ 2 ˆ 2 AðmAÞ¼ 2 dpp nFð p þ mAÞ: ð22Þ π 0 2 þ þ The electronic channel for the pion decay, π → e þ νe is 3 −1 helicity suppressed (Γπþ→ þ ν ¼ 4.7 × 10 s ) and of the order e þ e In this equation, g and mˆ ≡ m =T are the spin ≃ 50 A A A of the Hubble rate at temperatures T MeV; hence one might nˆ 0 x worry that leptons with electronic flavor depart from equilibrium. degeneracy and mass respectively, and Fð Þ¼ This is resolved by the observation that they are coupled to ðd=dxÞf1=½expðxÞþ1 g is the derivative of the Fermi- muonic species by other nonhelicity-suppressed, and conse- Dirac distribution. The strongly interacting fluid’s suscep- þ ν ↔ μþ ν quently faster, reactions such as e þ e þ μ and tibilities are considerably more complicated, especially þ þ ¯ μ ↔ e þ νe þ νμ. 3 near the quark-hadron transition. We evaluate them using We do not follow the strangeness, S, since it is not conserved a number of techniques: perturbative quantum chromody- in weak reactions. Above the transitions, it disappears at the Γ ≃ 2 2 5 ≈ 1013 −1 5 namics (QCD) at high temperatures, matching to lattice Cabbibo-suppressed rate S jVusj GFT s ðT=GeVÞ , while below the transition the relevant rate is the kaon inverse QCD results near the transition, and a hadron resonance gas Γ 8 1 107 −1 lifetime, K ¼ . × s . (HRG) approximation at low temperatures.
043515-6 STERILE NEUTRINO DARK MATTER: WEAK … PHYSICAL REVIEW D 94, 043515 (2016) B. Susceptibilities of the strongly interacting plasma g4 ~2 2 s αMS g3 ¼ gs þ 7 ; ð27Þ In this section, we compute the strongly interacting ð4πÞ2 E plasma’s susceptibilities to baryon number and electric αMS charge fluctuations. The susceptibilities are the following and the functions En are given by derivatives of the QCD pressure, π2 XNf MS 2 2 α d 4C F1 mˆ ; μˆ ; ∂ pˆ E1 ¼ A 45 þ A ð i iÞ ð28Þ χˆX QCD 1 2 ¼ 2 ; ð23Þ i¼ ∂μˆ μˆ 0 X X¼ XNf 1 MS 2 2 2 α −d F2 mˆ ; μˆ 1 6F2 mˆ ; μˆ ∂ pˆ E2 ¼ A 6 ð i iÞ½ þ ð i iÞ χˆXY QCD 1 11 ¼ ; ð24Þ i¼ ∂μˆ ∂μˆ μˆ μˆ 0 2 X Y X; Y ¼ ˆ μ mi 3 2 ˆ 2 μˆ þ 2 ln þ F2ðmi ; iÞ 4π mi where X; Y ∈ fB; Qg; μˆ X ≡ μX=T is the chemical potential ˆ d C of the conserved charge X; and the pressure pQCD is given − 2 ˆ ˆ 2 μˆ − A A miF4ðmi ; iÞ 144 ; ð29Þ by the logarithm of the QCD partition function ZQCD, γ 1 22C μe E 1 pQCD αMS A pˆ ≡ ¼ ln Z ðV;T;μ ; μ Þ; ð25Þ E7 ¼ ln þ QCD T4 VT3 QCD Q B 3 4πT 22 2 XNf μ − 2 ˆ 2 μˆ where V is the volume. In Eq. (24), the off-diagonal term ln þ F3ðmi ; iÞ : ð30Þ XY 3 mi χˆ11 encodes the correlation between the fluctuations of i¼1 conserved charges X and Y. Note that the susceptibilities in ≡ 2 − 1 ≡ Eqs. (23) and (24) are dimensionless. Here, dA Nc and CA Nc stand for the gauge-group The sterile neutrino production calculation carried in the constants for the adjoint and fundamental representation of present work requires knowledge of these susceptibilities SUðNcÞ, respectively. In this work, we adopt the standard 3 over a broad range of temperatures, both above and below value of Nc ¼ . In the above, Nf is the number of quark μ the quark-hadron transition. At very high temperatures flavors, is the energy scale at which the masses and the T ≫ T , the QCD pressure can be computed using a coupling constant are evaluated (not to be confused with QCD γ standard perturbative approach, while at intermediate tem- the chemical potentials), and E is the Euler-Mascheroni … peratures T ≃ T , perturbative techniques become inad- constant. The functions F1; ;F4 are given in QCD Appendix A. equate; we must rely on lattice calculations (see e.g. We also need a prescription for the running of the [100,101]) to compute susceptibilities through the quark- coupling constant g and of the quark masses with the hadron transition. At low temperatures T