PHYSICAL REVIEW D 103, 103528 (2021)

Cosmological and magnetic fields from violating particle decays

Tanmay Vachaspati 1 and Alexander Vilenkin2 1Physics Department, Arizona State University, Tempe, Arizona 85287, USA 2Institute of Cosmology, Department of Physics and Astronomy, Tufts University, Medford, Massachusetts 02155, USA

(Received 1 February 2021; accepted 30 April 2021; published 20 May 2021)

We estimate the chirality of the cosmological plasma due to parity violating decays of particles, focusing on the example of . The nontrivial chirality is however too small to make a significant contribution to the cosmological magnetic field via the chiral-magnetic effect.

DOI: 10.1103/PhysRevD.103.103528

I. INTRODUCTION particle helicity to the cosmological plasma. The net particle helicity in principle leads to electric currents via The last few decades have seen growing interest in the CME that can generate magnetic helicity. However, cosmic magnetic fields on several fronts [1]. Several ideas accounting only for decays of Standard Model particles, the have been proposed that can generate magnetic fields in net particle helicity is too small to significantly affect cosmology, some of which are directly tied to known cosmological magnetic fields and their helicity. particle physics [2–4] and its possible extensions [1,5–8]. We start by describing the physical effect in some detail The magneto-hydrodynamical (MHD) evolution of cos- in the context of the tau in Sec. II, where we also mological magnetic fields is now understood quite well on estimate the induced electric currents. We find an upper the basis of analytical arguments [9,10] and direct simu- bound to the magnetic helicity that can be generated due to lations [11]. There are claims for an indirect lower bound chiral effects in Sec. III. Our conclusions are summarized on the cosmological magnetic field strength [12–16], and discussed in Sec. IV. though not without debate [17,18], and more direct evidence [19]. Concurrently there are claims of a parity violating signature that can be used to measure the II. CHIRALITY PRODUCTION IN TAU DECAYS magnetic field helicity spectrum [20,21] though with no To illustrate the physics of the effect, in this section we significant detections as yet [22,23]. will discuss the decay of tau leptons in the background of a In parallel to these developments, motivated by heavy- magnetic field and fluid vorticity. Except for small ion collision experiments [24], there has been renewed differences, the physics carries over to the case of decays interest in chiral effects in plasmas, namely the chiral- of other particles. magnetic [25] and chiral-vortical [26] effects (CME and CVE respectively). The CME and CVE have also been applied to the evolution of cosmological and astrophysical A. Particle decay magnetic fields [5,27–33]. In this paper we discuss how Tau leptons decay into (or ) and CME and CVE can effectively arise in standard cosmology with standard particle interactions due to the parity-violat- − − ing decays of heavy leptons and . The basic idea is τ → e þ ντ þ ν¯e ð1Þ that the Standard Model has a number of unstable particles that decay at various cosmological epochs, primarily due to and antitau into and neutrinos the weak interactions. Since the weak interactions violate parity, the decay products are chiral and this provides a net þ þ τ → e þ ν¯τ þ νe: ð2Þ

These decays violate parity since they proceed primarily by Published by the American Physical Society under the terms of the weak interactions. Therefore the tau predominantly the Creative Commons Attribution 4.0 International license. Further distribution of this work must maintain attribution to decays into a relativistic left-handed , while an the author(s) and the published article’s title, journal citation, antitau decays into a relativistic right-handed . Due and DOI. Funded by SCOAP3. to the lepton of the Universe there are more taus

2470-0010=2021=103(10)=103528(5) 103528-1 Published by the American Physical Society TANMAY VACHASPATI and ALEXANDER VILENKIN PHYS. REV. D 103, 103528 (2021) than antitaus, and the cosmological plasma gains net left- 1 T2 τ ∼ : ð Þ χ 2 6 handed chirality as taus decay. αT me The decay product electrons are chiral since they are produced by the weak interactions, but chirality is not The excess of antitaus over taus, δnτ, decreases due to preserved for massive particles. Instead, as emphasized in tau decay and is described by the equation Ref. [34] in the context of supernovae and neutron stars, 3 chirality is nearly equal to helicity for ultrarelativistic d 3 a eq ða δnτÞ¼ ðδnτ − δnτÞ: ð7Þ particles, so it is better to think of the final electrons as dt τd being in a definite helicity state. Helicity can only change due to particle interactions. We shall adopt this view in At temperatures below the electroweak phase transition, T ≲ T ∼ 100 τ ≪ t t what follows. EW GeV, we have d , where is the 2 The τ mass is mτ ¼ 1777 MeV and the τ lifetime in its cosmic time. This means that the equilibrium density of −13 rest frame is ττ ¼ 2.9 × 10 s. However, the decaying taus establishes very quickly (compared to the Hubble eq taus are constantly reproduced by reactions inverse to (1), time), and the approximate solution of (7) is δnτ ≈ δnτ . 1 eq (2), so the number density of taus, nτ, remains comparable Inserting (7) in (4) and then using δnτ ≈ δnτ we have to that of until the time a3n d 3 d 3 eq χ −7 ða nχÞ¼− ða δnτ Þ − : ð8Þ tτ ∼ 10 s; ð3Þ dt dt τχ

eq when the cosmic temperature drops to T ∼ mτ. At later With a given δnτ , this equation can be solved in quad- times nτ decreases exponentially. ratures, but we shall instead find an approximate solution. The particle helicity density, nχ, is produced in tau Since we are in the regime where τχ ≪ t, the term on the decays and is dissipated by helicity flipping scatterings and left-hand side can be neglected and we obtain 3 2 due to the chiral . The latter is proportional to α B eq [35], where α ≈ 1=137 is the fine structure constant and B is 3 d δnτ nχ ≈−τχT ; ð9Þ the magnetic field strength, and is much slower than dt T3 helicity flipping scatterings for vanishing or weak magnetic fields. We will ignore the anomalous flipping for now but where we have used aT ≈ const. will discuss it in Sec. III when we consider the effect of Once we determine the equilibrium excess of antitaus eq particle chirality on the generation of magnetic fields. The over taus, denoted by δnτ , we can estimate the chirality evolution of particle helicity density can be described by density of the Universe due to tau decays using (9). the kinetic equation in the relaxation time approximation, B. Equilibrium density 3 3 d a a nχ eq ða3n Þ¼ ðδn − δneqÞ − ; ð Þ The equilibrium density δnτ is given by dt χ τ τ τ τ 4 d χ Z 1 ∞ E − μ E þ μ δneq ¼ dpp2 f τ − f τ ; ð Þ τ 2 10 where 2π 0 T T

þ − x −1 δnτ ¼ nτ − nτ ; ð5Þ wherepffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffifðxÞ¼ðe þ 1Þ is the Fermi distribution, 2 2 E ¼ p þ mτ , and μτ is the chemical potential of τ − þ nτ and nτ are the densities of tau and antitau particles, particles. At T ≫ mτ; μτ we can expand the integrand in eq 2 2 respectively, δnτ is the equilibrium value of δnτ, τd ∼ Eq. (10) in powers of mτ =p and μτ=T. The integrations are ðT=mτÞττ is the decay time of taus (assuming that T>mτ then easily performed and we find −1 and with time dilation taken into account) and τχ is the T ≫ m μ T2 3m2 electron helicity flipping rate. At e, the helicity δneq ≈ τ 1 − τ : ð Þ 2 2 τ 2 2 11 flipping rate is suppressed by a factor me=T compared to 6 2π T the scattering rate α2T and a correction factor of α−1 due to infrared effects [36] (earlier estimates of the flipping rate We assume that the baryon and/or lepton asymmetry of T ≫ T were suppressed by another factor of α [34]), the Universe was generated at EW by some inter- actions beyond the Standard Model, for example by 1Tau particles are also produced and destroyed in scattering pffiffiffiffi τ− þ ν → e− þ ν 2 t ∼ m = NT2 reactions like e τ. We disregard them in what This is easily verified using the relation P , where follows, since they do not change the order of magnitude of the mP is the Planck mass and N is the number of particle species in effect. equilibrium.

103528-2 COSMOLOGICAL CHIRALITY AND MAGNETIC FIELDS FROM … PHYS. REV. D 103, 103528 (2021)

(B − L)-violating leptoquark decays. This asymmetry was where h ¼hA · Bi is the magnetic helicity. Assuming that then redistributed between the Standard Model leptons and the initial magnetic helicity and the final particle helicity quarks by sphaleron processes. Let us denote the lepton-to- vanish, we get T ≪ T η ratio at EW by L. Generically we expect the α η m2 T chemical potentials of light baryons and leptons to be of the h ¼ n ∼ L τ n ð Þ μ =T ∼ η ∼ η η ∼ 10−9 max 4π χ 2 m γ 16 same order, and τ L B [37,38], where B 4πme P is the observed baryon-to-photon ratio. However, observa- h tions allow for ηL to be much larger, with an upper bound where we have used (14). We compare max to magnetic −2 derived from big bang nucleosynthesis of ηL ≲ 10 [39]. helicity that could be induced due to [3,4], T In the high-temperature regime, when is large compared η μ =T ≈ B −5 −3 −45 2 to all relevant particle masses, we have τ const, with a hB ∼ nγ ∼ 10 cm ∼ 10 G Mpc ð17Þ 2 2 α mass correction Oðm =T Þ [40]. Then Eq. (11) becomes

δneq m2 where we have used the known cosmic τ ≈ Cη − Kη τ ; ð Þ density and are using natural units. Then T3 L L T2 12 η αm2 T where C and K are Oð1Þ numerical constants.3 The mass h ∼ L τ h ∼ 10−6h ∼ 10−51 2 ð Þ max η 4πm2 m B B G Mpc 18 correction term in (12) can be qualitatively understood as B e P follows. As the temperature decreases, it becomes ener- T ∼ 100 getically favorable to transfer the conserved τ-lepton where we have used GeV in the numerical η ∼ 10−2 number from τ particles to τ neutrinos. The excess τ-lepton estimate and the maximum allowed value L . ∼175 number is also decreased as a result [40]. Even if the decay of top quarks with mass GeV ∼1 Substituting Eq. (12) in (9) we obtain to down quarks with mass MeV is considered, h ∼ 10−3h ∼ 10−48 2 max B G Mpc. 2 _ nχ ≈−3KηLτχmτ T: ð13Þ Comparing to observations, even with the most conservative lower bound of 10−19 G on Mpc scales, we T_ ¼ −T=2t t ∼ m =T2 τ get an estimate for the observed helicity ∼10−38 G2 Mpc With , P and χ from Eq. (6), this gives (omitting numerical factors) which is 10 orders of magnitude larger than the most optimistic estimate above. Alternately, we can estimate the 2 ηLmτ T magnetic field strength on Mpc scales implied by the n ∼ n ; ð Þ −24 χ 2 m γ 14 helicity and this gives 10 G as compared to the observed αme P lower bound of 10−19 G. 3 where nγ ∼ T is the photon number density. This estimate was derived assuming T ≫ mτ, but it still IV. CONCLUSIONS applies at T ∼ mτ. Reactions (1), (2) remain in equilibrium We have shown that the decays of certain Standard when T drops well below mτ. In this regime, δnτ and nχ Model particles can lead to a chiral cosmological plasma decrease exponentially. around the epoch of the electroweak phase transition. The Similar formulas can be written down for the decay of final result for the chiral asymmetry due to tau-lepton other unstable particles. The largest helicity is injected by decays is given in (14). However, the asymmetry is sup- the decay of the heaviest particle into the lightest particle −2 pressed by the lepton to photon ratio (ηL ≲ 10 ) and the and at the highest temperature. effect on magnetic field helicity generation is very weak as we have shown in Sec. III. Nonetheless it may be of interest III. MAGNETIC HELICITY that the cosmological plasma was chiral at the earliest As noted in Ref. [32], the maximum magnetic helicity moments even within the Standard Model of particle that can be obtained from particle helicity can be derived physics. A chiral medium may also arise in extensions from the equation, which can be written as a of the Standard Model that include parity violating conservation law, decays of very heavy particles to light Standard Model T at temperatures greater than EW. Such decays 4π would result in stronger magnetic fields with larger nχ þ h ¼ constant; ð15Þ α magnetic helicity.

3This estimate assumes that taus are the heaviest particles ACKNOWLEDGMENTS present in equilibrium at temperature T. If a heavier particle is present in equilibrium, it too will contribute to the mass We thank the participants of the Nordita workshop on correction and may change the estimate. “Quantum Anomalies and Chiral Magnetic Phenomena,”

103528-3 TANMAY VACHASPATI and ALEXANDER VILENKIN PHYS. REV. D 103, 103528 (2021) especially Axel Brandenburg and Kohei Kamada for feed- Department of Energy, Office of High Energy Physics, back. We also thank Matt Baumgart, Andrew Long, under Award No. DE-SC0019470 at Arizona State Cecilia Lunardini, Igor Shovkovy, and Tracy Slatyer for University. A. V. is supported by the National Science discussions. The work of T. V. is supported by the U.S. Foundation Grant No. 1820872.

[1] T. Vachaspati, Progress on cosmological magnetic fields, extragalactic gamma-ray background, Astrophys. J. 752,22 arXiv:2010.10525. (2012). [2] T. Vachaspati, Magnetic fields from cosmological phase [18] T. C. Arlen, V. V. Vassiliev, T. Weisgarber, S. P. Wakely, and transitions. Phys. Lett. B 265, 258 (1991). S. Yusef Shafi, Intergalactic magnetic fields and gamma ray [3] J. M. Cornwall, Speculations on primordial magnetic hel- observations of extreme TeV blazars, Astrophys. J. 796,18 icity, Phys. Rev. D 56, 6146 (1997). (2014). [4] T. Vachaspati, Estimate of the Primordial Magnetic Field [19] W. Chen, J. H. Buckley, and F. Ferrer, Search for GeV? Ray Helicity, Phys. Rev. Lett. 87, 251302 (2001). Pair Halos around Low Redshift Blazars, Phys. Rev. Lett. [5] M. Joyce and M. E. Shaposhnikov, Primordial Magnetic 115, 211103 (2015). Fields, Right-Handed Electrons, and the Abelian Anomaly, [20] H. Tashiro, W. Chen, F. Ferrer, and T. Vachaspati, Search for Phys. Rev. Lett. 79, 1193 (1997). CP violating signature of intergalactic magnetic helicity in [6] M. McNeil Forbes and A. R. Zhitnitsky, Primordial Galactic the gamma ray sky, Mon. Not. R. Astron. Soc. 445, L41 Magnetic Fields from Domain Walls at the QCD Phase (2014). Transition, Phys. Rev. Lett. 85, 5268 (2000). [21] W. Chen, B. D. Chowdhury, F. Ferrer, H. Tashiro, and T. [7] T. Stevens, M. B. Johnson, L. S. Kisslinger, and E. M. Vachaspati, Intergalactic magnetic field spectra from diffuse Henley, Non-Abelian Higgs model of magnetic field gen- gamma rays, Mon. Not. R. Astron. Soc. 450, 3371 (2015). eration during a cosmological first-order electroweak phase [22] J. Asplund, G. Johannesson, and A. Brandenburg, On the transition, Phys. Rev. D 85, 063003 (2012). measurement of handedness in Fermi Large Area Telescope [8] F. Miniati, G. Gregori, B. Reville, and S. Sarkar, Axion- data, Astrophys. J. 898, 124 (2020). Driven Cosmic Magnetogenesis during the QCD Crossover, [23] M. Kachelriess and B. C. Martinez, Searching for primor- Phys. Rev. Lett. 121, 021301 (2018). dial helical magnetic fields, Phys. Rev. D 102, 083001 [9] R. Banerjee and K. Jedamzik, Evolution of cosmic magnetic (2020). fields: From the very early universe, to recombination, to the [24] D. E. Kharzeev, The chiral magnetohydrodynamics of QCD present, Phys. Rev. D 70, 123003 (2004). fluid at RHIC and LHC, J. Phys. G 38, 124061 (2011). [10] K. Jedamzik and G. Sigl, The evolution of the large-scale [25] A. Vilenkin, Equilibrium parity violating current in a tail of primordial magnetic fields, Phys. Rev. D 83, 103005 magnetic field, Phys. Rev. D 22, 3080 (1980). (2011). [26] A. Vilenkin, Macroscopic parity violating effects: [11] A. Brandenburg, T. Kahniashvili, S. Mandal, A. R. Pol, fluxes from rotating black holes and in rotating thermal A. G. Tevzadze, and T. Vachaspati, Evolution of hydro- radiation, Phys. Rev. D 20, 1807 (1979). magnetic turbulence from the electroweak phase transition, [27] A. Boyarsky, J. Frohlich, and O. Ruchayskiy, Self-Consis- Phys. Rev. D 96, 123528 (2017). tent Evolution of Magnetic Fields and Chiral Asymmetry in [12] A. Neronov and I. Vovk, Evidence for strong extragalactic the Early Universe, Phys. Rev. Lett. 108, 031301 (2012). magnetic fields from Fermi observations of TeV blazars, [28] H. Tashiro, T. Vachaspati, and A. Vilenkin, Chiral effects Science 328, 73 (2010). and cosmic magnetic fields, Phys. Rev. D 86, 105033 [13] W. Essey, S. Ando, and A. Kusenko, Determination of (2012). intergalactic magnetic fields from gamma ray data, Astro- [29] M. Dvornikov and V. B. Semikoz, Lepton asymmetry part. Phys. 35, 135 (2011). growth in the symmetric phase of an electroweak plasma [14] K. Dolag, M. Kachelriess, S. Ostapchenko, and R. Tomas, with hypermagnetic fields versus its washing out by Lower limit on the strength and filling factor of extragalactic sphalerons, Phys. Rev. D 87, 025023 (2013). magnetic fields, Astrophys. J. 727, L4 (2011). [30] M. Dvornikov and V. B. Semikoz, Instability of magnetic [15] M. Wood, J. Biteau, R. Caputo, M. Di Mauro, and M. fields in electroweak plasma driven by neutrino asymme- Meyer, Preliminary results of the Fermi high-latitude ex- tries, J. Cosmol. Astropart. Phys. 05 (2014) 002. tended source catalog, Proc. Sci., ICRC2017 (2018) 647. [31] M. Dvornikov and V. B. Semikoz, Influence of the turbulent [16] M. Ackermann et al., The search for spatial extension in motion on the chiral magnetic effect in the early universe, high-latitude sources detected by the Fermi large area Phys. Rev. D 95, 043538 (2017). telescope, Astrophys. J. Suppl. 237, 32 (2018). [32] A. Brandenburg, J. Schober, I. Rogachevskii, T. Kahniashvili, [17] A. E. Broderick, P. Chang, and C. Pfrommer, The cosmo- A. Boyarsky, J. Frohlich, O. Ruchayskiy, and N. Kleeorin, logical impact of luminous TeV blazars I: Implications of The turbulent chiral-magnetic cascade in the early universe, plasma instabilities for the intergalactic magnetic field and Astrophys. J. 845, L21 (2017).

103528-4 COSMOLOGICAL CHIRALITY AND MAGNETIC FIELDS FROM … PHYS. REV. D 103, 103528 (2021)

[33] Y. Masada, K. Kotake, T. Takiwaki, and N. Yamamoto, [37] V. A. Kuzmin, V. A. Rubakov, and M. E. Shaposhnikov, Chiral magnetohydrodynamic turbulence in core-collapse Anomalous electroweak baryon number nonconservation supernovae, Phys. Rev. D 98, 083018 (2018). and GUT mechanism for baryogenesis, Phys. Lett. B 191, [34] D. Grabowska, D. B. Kaplan, and S. Reddy, Role of the 171 (1987). electron mass in damping chiral plasma instability in super- [38] J. A. Harvey and M. S. Turner, Cosmological baryon and novae and neutron stars, Phys. Rev. D 91, 085035 (2015). in the presence of electroweak [35] D. G. Figueroa and M. Shaposhnikov, Anomalous non- number violation, Phys. Rev. D 42, 3344 (1990). conservation of fermion/chiral number in Abelian gauge [39] D. J. Schwarz and M. Stuke, Does the CMB prefer a theories at finite temperature, J. High Energy Phys. 04 leptonic Universe? New J. Phys. 15, 033021 (2013). (2018) 026. [40] A. I. Bochkarev, S. Yu.Khlebnikov, and M. E. Shaposhnikov, [36] A. Boyarsky, V. Cheianov, O. Ruchayskiy, and O. Sobol, Sphalerons and baryogenesis: Electroweak CP violation at Evolution of the Primordial Axial Charge across Cosmic high temperatures, Nucl. Phys. B329, 493 (1990). Times, Phys. Rev. Lett. 126, 021801 (2021).

103528-5