Temperature Dependence of the Neutron Lifespan Cheng Tao Yanga, Jeremiah Birrellb, Johann Rafelskia aDepartment of Physics, The University of Arizona, Tucson, Arizona 85721, USA bDepartment of Mathematics and Statistics, University of Massachusetts Amherst, Amherst, MA 01003, USA (Dated: May 18, 2018) Current precision big bang nucleosynthesis (BBN) studies motivate us to revisit the neutron lifespan in the plasma medium of the early universe. The mechanism we explore is the Fermi-blocking of decay electrons and neutrinos by plasma. As result, neutrons live longer and we find a significant 6.4% modification of neutron abundance in the BBN era arising from in plasma modification of the neutron lifetime. This effect can influence the final abundances of the light elements in BBN. I. INTRODUCTION and the neutron-proton mass difference (in medium). f The values considered are in the range Xn = 0:15 ∼ 0:17 Element production during BBN is influenced by sev- [1]. A dynamical approach to neutron freezeout is nec- essary to fully understand Xf , we hope to return to this eral parameters, e.g. baryon to photon ratio ηb, num- n challenge in the near future. ber of neutrino species Nν , and neutron to proton ratio, as controlled by both the dynamics of neutron freeze- Following freeze-out the neutron is subject to natural decay and normally the neutron lifetime in vacuum τ 0 out at temperature Tf ≈ 0:8MeV and neutron lifetime. n Since about 200 seconds pass between neutron freeze- is used c.f. Eq. (2) to calculate the neutron abundance out and BBN at T ≈ 0:07MeV, the neutron lifetime resulting in the `desired' value Xn(TBBN ) ≈ 0:13 when is one of the important parameter controlling BBN el- BBN starts. To improve precision a dynamically evolving ement yields [1]. However, the neutron decay in cosmic neutron yield needs to be studied and for this purpose we plasma medium has seemingly never been considered be- explore here the neutron decay which occurs in medium, fore. Therefore, here we study the Fermi suppression ef- not vacuum. This leads to neutron lifespan dependence fects which lengthen the neutron lifetime in the early uni- on temperature of the cosmic medium as the decay occurs verse. The medium dependence of particle decay was rec- for a particle emerged in plasma of electron/positron, ognized by Kuznetsova et al [2], we will use their method neutrino/antineutrino, (and protons). to study the plasma effects on neutron decay in the cos- Two physical effects of the medium influence the neu- mic medium. tron lifetime in the early universe noticeably: After freeze-out, the neutron abundance is subject to • Fermi suppression factors from the medium: Dur- natural decay ing the temperature range Tf > T > TBBN , elec- trons and neutrinos in the background plasma can n −! p + e + νe ; (1) reduce the neutron decay rate by Fermi suppression where the current experimental neutron lifetime quoted to the neutron decay rate. Furthermore, the neu- 0 by the Particle Data Group [3] is τn = 880:2±1:0 sec. The trino background can still provide the suppression latest measurement for neutron lifetime by using magne- after electron/positron pair annihilation becomes togravitational trap is 877:7±0:7 sec [4]. In the standard nearly complete. Big Bang Nucleosynthesis, the neutron abundance when nucleosynthesis begins is given by [1] • Photon reheating: When T me the elec- tron/positron annihilation occurs, the entropy from t − t e± is fed into photons, leading to photon reheating. X (T ) = Xf exp − BBN f ≈ 0:13 ; (2) n BBN n τ 0 The already decoupled (frozen-out) neutrinos re- n main undisturbed. Therefore, after annihilation we arXiv:1805.06543v1 [nucl-th] 16 May 2018 f The normalizing neutron freeze-out yield Xn have two different temperatures in cosmic plasma: neutrino temperature Tν and the photon and pro- f f f f nn nn=np ton temperature T respectively. Xn ≡ f f = f f : (3) nn + np 1 + nn=np These two effect will be included in the following explo- f f ration of the neutron lifetime in the early universe as a where nn and np are freeze-out neutron and proton den- sities, respectively. The thermal equilibrium yield ratio function of T . We show how these effects alter the neu- is tron lifespan and obtain the modification of the neutron yield at the time of BBN. Yet another effect was consid- nf ered by Kuznetsova et al [2] which is due to time dilation n = exp (−Q=T ) ;Q = m − m ; (4) f f n p originating in particle thermal motion. In our case for np neutrons with T=m < 10−3 this effect is negligible. Be- assuming a instantanous freeze-out, depends on temper- low we will explicitly assume that the neutron decay is ature Tf at which neutrons decouple from the heat bath, studied in the neutron rest frame. 2 In section II we discuss the general form of the neutron In the following calculation, we assume the condition decay rate in a heat bath, including the Fermi suppres- Tk µν¯; mν , i.e. we consider the massless neutrino sion factors from the medium. In section III we derive in plasma. Substituting the distributions into the decay the relation between neutrino temperature and photon rate formula and using the neutron rest frame, the decay temperature in early universe. In section IV we com- rate can be written as pute the modified neutron abundance due to effects of 2 5 2 the medium and compare with the standard (vacuum) 1 GF Q Vud 2 0 = 3 (1 + 3gA) (1 + RC) (11) result. The implications of these results will be discussed τn 2π in section V. Z 1 ξ(1 − ξ)2 pξ2 − (m =Q)2 × dξ e ; me=Q exp (−Qξ=T ) + 1 exp (−Q(1 − ξ)=Tν ) + 1 II. DECAY RATE IN MEDIUM where Q was defined in Eq. (4) and we integrate using dimensionless variable ξ = Ee=Q. From Eq.(11), the The invariant matrix element for the neutron decay decay rate in vacuum can be written as Eq. (1) for non-relativistic neutron and proton is given 2 5 2 by 1 GF meVud 2 0 0 = 3 (1 + 3gA) (1 + RC) f ; (12) 2 2 2 2 τn 2π hjMj i ≈ 16 GF Vud mnmp(1 + 3gA)(1 + RC)Eν¯Ee; (5) where the phase space factor f 0 is given by where the Fermi constant is GF = 1:1663787 × 10−5 GeV−2, V = 0:97420 is an element of the ud Q 5 Z 1 Cabibbo-Kobayashi-Maskawa (CKM) matrix [5{7], and f 0 ≡ dξ ξ(1 − ξ)2pξ2 − (m =Q)2 m e gA = 1:2755 is the axial current constant for the nucle- e me=Q ons [5, 8]. We also consider the total effect of all radiative = 1:6360: (13) corrections relative to muon decay that have not been absorbed into Fermi constant GF . The most precise cal- The phase space factor is also modified by the Coulomb culation of this correction [6, 8] gives (1+RC) = 1:03886. correction between electron and proton, proton recoil, In the early universe the neutron decay rate in medium, nucleon size correction etc. This has been studied by at finite temperature can be written as [2] Wilkinson [10], and the phase space factor is given by [5, 7, 10] Z 3 3 3 1 1 d pν¯ d pp d pe 0 = 3 3 3 τn 2mn (2π) 2Eν¯ (2π) 2Ep (2π) 2Ee f = 1:6887: (14) 4 4 2 (2π) δ (pn − pp − pe − pν¯) hjMj i These effect amount to adding the factor F to our calcu- 1 − fp(pp) 1 − fe(pe) 1 − fν¯(pν¯) ; (6) lation where we consider this expression in the rest rest frame f F = 0 = 1:0322; (15) of neutron, i.e. pn = (mn; 0). The phase-space factors f (1 − fi) are Fermi suppression factors in the medium. The Fermi-Dirac distributions for electron and non- then the neutron lifespan can be written as relativistic proton are given by τ 0 τ Vacuum = n = 879:481 sec; (16) 1 n F fe = ; (7) eEe=T + 1 2 which compare well to the experiment result 877:7 ± −Ep=T −mp=T −p =2mpT fp = e = e e p : (8) 0:7 sec [4]. In the case of plasma medium, we do not ex- pect that these effect (Coulomb correction between elec- For neutrinos, after neutrino/antineutrino kinetic freeze- tron and proton, proton recoil, nucleon size correction out they become free streaming particles. If we assume etc) are modified in the cosmic plasma. Thus we adapt that kinetic freeze out occurs at some time t and temper- k the factor into our calculation and the neutron decay rate ature T , then for t > t the free streaming distribution k k in the cosmic plasma is given by function can be written as [9] 2 5 2 1 1 GF Q Vud 2 f = ; (9) = (1 + 3gA) (1 + RC) F (17) ν¯ q 2 2 2 Medium 3 E −mν mν µν¯ τn 2π exp T 2 + T 2 + T + 1 ν k k Z 1 ξ(1 − ξ)2 pξ2 − (m =Q)2 × dξ e : for antineutrinos and we define the effective neutrino me=Q exp (−Qξ=T ) + 1 exp (−Q(1 − ξ)=Tν ) + 1 temperature Tν as From Eq.(17) we see that the neuron decay rate in the a(tk) Tν ≡ Tk: (10) early universe depends on both the photon temperature a(t) T and the neutrino effective temperature Tν . 3 III. PHOTON REHEATING ± After neutrino free-out and when me T , the e be- comes non-relativistic and annihilate.