The Neutron Mean Life and Big Bang Nucleosynthesis

We explore the effect of neutron lifetime and its uncertainty on standard big-bang nucleosynthesis (BBN). BBN describes the cosmic production of the light nuclides $^1{\rm H}$, ${\rm D}$, $^3{\rm H}$+$^3{\rm He}$, $^4{\rm He}$, and $^7{\rm Li}$+$^7{\rm Be}$ in the first minutes of cosmic time. The neutron mean life $\tau_n$ has two roles in modern BBN calculations: (1) it normalizes the matrix element for weak $n \leftrightarrow p$ interconversions, and (2) it sets the rate of free neutron decay after the weak interactions freeze out. We review the history of the interplay between $\tau_n$ measurements and BBN, and present a study of the sensitivity of the light element abundances to the modern neutron lifetime measurements. We find that $\tau_n$ uncertainties dominate the predicted $^4{\rm He}$ error budget, but these theory errors remain smaller than the uncertainties in $^4{\rm He}$ observations, even with the dispersion in recent neutron lifetime measurements. For the other light-element predictions, $\tau_n$ contributes negligibly to their error budget. Turning the problem around, we combine present BBN and cosmic microwave background (CMB) determinations of the cosmic baryon density to $\textit{predict}$ a"cosmologically preferred"mean life of $\tau_{n}({\rm BBN+CMB}) = 870 \pm 16 \ \rm sec$, which is consistent with experimental mean life determinations. We go on to show that if future astronomical and cosmological helium observations can reach an uncertainty of $\sigma_{\rm obs}(Y_p) = 0.001$ in the $^4{\rm He}$ mass fraction $Y_p$, this could begin to discriminate between the mean life determinations.

BBN calculations: (1) it normalizes the matrix element for weak n ↔ p interconversions, and (2) it sets the rate of free neutron decay after the weak interactions freeze out.We review the history of the interplay between τ n measurements and BBN, and present a study of the sensitivity of the light element abundances to the modern neutron lifetime measurements.We find that τ n uncertainties dominate the predicted 4 He error budget, but these theory errors remain smaller than the uncertainties in 4 He observations, even with the dispersion in recent neutron lifetime measurements.For the other light-element predictions, τ n contributes negligibly to their error budget.
Turning the problem around, we combine present BBN and cosmic microwave background (CMB) determinations of the cosmic baryon density to predict a "cosmologically preferred" mean life of τ n (BBN + CMB) = 870 ± 16 sec, which is consistent with experimental mean life determinations.
We go on to show that if future astronomical and cosmological helium observations can reach an uncertainty of σ obs (Y p ) = 0.001 in the 4 He mass fraction Y p , this could begin to discriminate between the mean life determinations.

I. INTRODUCTION
The origin of the light elements ( 1 H, D, 3 He, 4 He, and 7 Li) is well explained by big bang nucleosynthesis (BBN) taking place in the early Universe when the temperature scale was roughly 1 MeV and below, i.e., the energies associated with nuclear reactions [1][2][3][4][5].Because the physics at these energies is well known, testable predictions of BBN are possible given reliable nuclear cross sections and the subject of this contribution, a well-measured neutron lifetime [6].
Standard BBN is the theoretical framework that implements these cross sections in the context of both the Standard Models of particle physics and cosmology, and is the model we will adopt for this paper.For cosmology, we will thus assume ΛCDM.For BBN, this means that the Universe was dominated by photon and neutrino radiation during BBN, with a baryon density consistent with that determined from measurements of the cosmic microwave background (CMB) anisotropies [7]; it is useful to parameterize the baryon density in terms of baryon-to-photon ratio η ≡ n b /n γ , with the density of blackbody photons fixed by the temperature: n γ ∝ T 3 .From the Standard model of electroweak interactions, we will assume that the number of light neutrino flavors is N ν = 3 and that there are no other relativistic contributions to the energy density.
Historically, progress in BBN has been tightly linked to progress in the measurement of the neutron lifetime.Until the 1980's, the uncertainties in the prediction of the light element abundances by BBN was dominated by the uncertainty in these three parameters (η, N ν , τ n ); this was prior to the CMB determination of the baryon density, the measurement of the number of neutrino flavors at LEP/SLAC [8], and the vast improvement in the measurement of the neutron lifetime [9,10].For example, the primordial helium mass fraction is sensitive to the variation in each of these quantities as can be seen in Fig. 1 taken from the 1981 analysis in Ref. [9].In this figure, the helium mass fraction, Y p = ρ( 4 He)/ρ baryon is plotted as a function of the baryon-to-photon ratio η for three values of N ν and three values for the neutron lifetime.At the time, BBN provided the best estimate for η and N ν (this was pre-LEP/SLAC and direct limits from accelerators on N ν were quite poor [11]).The values of τ n used here corresponded to what was the "accepted" value for the neutron half-life, τ 1/2 = 10.61 ± 0.16 min corresponding to a mean life τ n = 918 ± 14 s that came from a 1972 paper by Christensen et al. [12]; and more recent though very different determinations of τ 1/2 = 10.13 ± 0.09 min or τ n = 877 ± 8 s by Bondarenko et al. [13]; and τ 1/2 = 10.82 ± 0.20 min or τ n = 937 ± 17 s [14] by Byrne et al.
From Table 6 we see that a primordial helium abundance as low as Y P =0.22 may still be consistent with a nucleon-dominated universe and the standard model.For a nonbaryon-dominated universe (e.g., neutrino dominated), the constraints on rj are less stringent.Using the mass in galaxies inferred from mass-to-light ratios (tj>0.29X 10 _1 °), the standard model can produce Y P implying that nucleons alone cannot close the universe.With the assumption that most of the entropy in the universe today is in the photon and neutrino backgrounds (with N p = 3), the baryon-to-photon ratio translates into a ratio of the baryon number to the specific entropy of kn B /s=l0~l 0S±l .Since this ratio remains constant when the baryon number is effectively conserved, it is the input parameter for models of cosmological baryon generation.
If one assumes that the mass which binds BSG is primarily in the form of baryons, then only three  [17] which gave τ n = 875 ± 95 s.The Bondarenko et al. measurement was excluded as it was "in significant disagreement with the other two precise direct mean life measurements and the inferred value given by Wilkinson 80" [18].
In the 1984 edition of the Review of Particle Properties [19] the Bondarenko measurement was included replacing Kosvintsev 1980 giving τ n = 898 ± 16 s, noting that the origin of the discrepancies between the measurements was not known.Because of the discrepancies, the uncertainty in the mean was inflated by a scale factor of 2.4.
It wasn't until 1990, when Mampe et al. [20] presented a significantly more accurate result τ n = 887.6 ± 3.0 s using ultra cold neutrons in a fluid-walled bottle.This had a strong effect on the world average which was now 888.6 ± 3.5 based on seven measurements.The discrepancy was largely gone and the scale factor for the uncertainty was reduced to 1.3, in part due to the withdrawal of the 1980 measurement by Byrne et al.
The evolution of the world average is plotted in Fig. 2 showing the average mean lifetime since 1960.One clearly sees the marked drop in the value and uncertainty following the Mampe et al. measurement.The impact of this measurement on BBN was immediate [21] since, as described in more detail below, the uncertainty in the predicted helium abundance is very sensitive to the uncertainty in the neutron mean lifetime [22], meaning that a drop in the uncertainty from 16 to 3.5 reduced the uncertainty in Y P from 0.0032 to 0.0007 assuming Y p = 0.247.
By the time of the 2002 edition of the Review of Particle Physics [23] there were enough high quality direct measurements allowing the Particle Data Group (PDG) to drop measurements with uncertainties larger than 10s.This led to another impressive drop in the uncertainty of the mean, with τ n = 885.7 ± 0.8 s.This drop was largely driven by the measurement of Arzumanov et al. [24] and further dropped the uncertainty in Y P to 0.0002 making the dominant source of uncertainty (at the time) other nuclear rates [25].
At about this time, WMAP provided the first all-sky measurements of anisotropies in the CMB temperature [26].The these fluctuations encode a wealth of cosmological parameters, including the baryon density or baryon-to-photon ratio η.Using this input from the CMB along with the improvements in the measurements of τ n effectively made BBN a parameterfree theory [27].
However it is certainly reasonable to question whether the quoted uncertainties are truly reflective of systematic errors included.In 2005, Serebrov et al. [28] published a result which was in fact systematically lower the previous world average, τ n = 878.5 ± 0.7 ± 0.3.
This value was used in a BBN analysis in Ref. [29].As a sole and severely discrepant value, it was not included in the PDG average until 2012 [30] when additional low values were reported.The 2012 average was now τ n = 880.1 ± 1.1 where the uncertainty has been scaled by a factor of 1.8.This seismic shift in τ n is clearly seen in Fig. 2.An ideogram of the seven measurements leading to this average is shown in Fig. 3. Here, we see the onset of a discrepancy in τ n .The vertical blue rectangle corresponds to the 2012 mean and 1σ spread.
A review including a discussion of the methods used to obtain τ n can be found in [31].
Subsequent changes in the world average have been relatively minor-certainly from the perspective of BBN calculations [15], where the current value of based on eight measurements [6,28,[32][33][34][35][36][37][38] is used.As the uncertainties in individual measurements continue to drop, the discrepancy appears more pronounced and the uncertainty already includes a scale factor of 1.8 [6].This can be seen in the current ideogram shown in Fig. 4.
Despite the tightness of the world average in Eq. ( 2), the neutron lifetime remains an outstanding puzzle.The results used in Eq. ( 2) are based on confined ultracold neutrons, with earlier experiments using material bottles, and later experiments using magnetic and gravitational traps.The measurements in traps give consistent results, with recent very precise UCNτ determination [38] dominating the weighted average among them and leading to the sharp peak in Fig. 4. The most recent bottle measurements give lifetimes systematically longer than the trap measurements, leading to the shoulder to the right of the peak in Fig. 4.
The resulting dispersion between the trap and bottle measurements is of order ∆τ n ∼ 5 sec, and could stem from systematic differences between the methods.There is in addition an in-beam measurement by Yue et al. [39] with a relatively high value of τ n = 887.7 ± 1.2 ± 1.9 s.This differs from the best trap measurement by ∆τ n ∼ 10 sec.We will not here enter further into the current discrepancy, noting only that for a recent discussion of the differences between the ultracold neutron results (the ones used in the PDG average) and the in-beam results, see [40] as well as other articles in this Special Issue.
In what follows, we will first briefly review the role of the neutron lifetime in standard BBN calculations in section II.We will also review the current results of BBN based on the latest input nuclear cross-sections (including τ n ) and abundance data.In section III, we will discuss the sensitivity of the results to τ n .We will also test the potential effect of the τ n discrepancy by instead of Monte Carlo sampling of the Gaussian corresponding to the mean value of τ n (and its uncertainty), but as a test, we will use the ideogram in Fig. 4 to sample values of τ n -keeping in mind that this procedure is not fully rigorous.In section IV, we will take a novel approach to τ n and treat it again as a parameter for which we generate a likelihood function and make a BBN prediction for τ n .Finally in section V we summarize the current state of BBN with respect to the neutron lifetime.

II. STANDARD BBN
Standard BBN (SBBN) is built upon Standard Model of nuclear and particle physics, in the background of the Friedmann-Robertson-Walker cosmological model based on Einstein gravity.We assume only the standard set of nuclear and particle interactions and nuclear and particle content2 , in particular, with N ν = 3.Furthermore, in SBBN, we assume a radiation dominated Universe during the epoch of nucleosynthesis.The radiation density can be expressed as taking into account the contributions of photons, electrons and positrons, and neutrino flavors appropriate for temperatures T > 1 MeV.The expansion rate of the Universe, is determined by the Hubble parameter which can be expressed as where G N is Newton's constant and scales as H ∝ G 1/2 N T 2 in a radiation dominated universe.The cosmic evolution of the light nuclides is plotted in Fig. 5; many critical features can be understood analytically, as we now summarize.At temperatures T > ∼ 1 MeV, weak interactions between neutrons and protons maintain equilibrium.These are: As one might expect, the weak interaction rates scale as Γ wk ∝ G 2 F T 5 , where G F is the Fermi constant.These reactions freeze-out when their interaction rates become slower than the expansion rate of the Universe determined by the Hubble parameter, or in other words, when the mean time between interactions is longer than the age of the Universe, determined by H −1 .Thus, the freeze-out condition is set by Weak freeze-out occurs at T f 0.84 MeV, and corresponds to a flattening of the neutron and proton curves in is the neutron-proton mass difference.The resulting mass fraction of neutrons at freezeout is thus After freeze-out, free neutron decays reduce the ratio slightly, which leads to the shallow slope in the neutron and proton curves in Fig. 5. Accounting for these decays, X n,BBN = X n,f e −t BBN /τn [n/(n + p)] f e −t BBN /τn before nucleosynthesis begins at time t BBN .For t BBN = 322 s determined by the time when the photo-destruction rates of deuterium freeze-out, and τ n = 878.4s, we have X n,BBN = 1/8.17or (n/p) BBN = 1/7.17.A useful and more elaborate semi-analytic description of freeze-out can be found in [44,45].
At t BBN , nucleosynthesis activity reaches a crescendo.Deuterons are formed via np → dγ, and then undergo a series of strong reactions culminating in 4 He production.These rapid reactions lead to a dramatic drop in the free n abundance, as seen in Fig. 5. Finally, the strong reactions also freeze out, largely due to the inability to overcome the Coulomb barrier in the expanding and cooling plasma.Then the stable nuclides plateau, while the radioactive species ultimately decay-as most evidently seen by the free neutron decay in Fig. 5 at late times.
The neutron lifetime plays two roles in BBN calculations.First, it is used to normalize the zero-temperature matrix element for weak n − p interconversions (5).The rates for these reactions scale as Γ n↔p ∝ 1/τ n .Thus τ n affects the determination of the freeze-out temperature of these weak interactions.Freeze-out is determined from the competition between the weak interaction rates and the Hubble expansion as in Eq. 6.Thus an increase in the neutron lifetime leads to an increase in the freeze-out temperature leading to more neutrons at freeze-out (recall the n/p ratio scales as e −∆m/T prior to freeze-out) and hence more 4 He as seen in Figs. 1 and 6.
Secondly the neutron lifetime controls the rate of free neutron decay, which occurs between weak freezeout at t f ∼ 1 sec and the end of the D bottleneck at t d ≈ 320 sec.These decays lower the n/p ratio.So an increase in the neutron lifetime leads again to more neutrons left over and hence a higher 4 He abundance.
Similar arguments can be made for the dependence of Y p on N ν and η.As one can see from Fig. 6, Y p increases with increases in all three inputs: η, N ν , and τ n .It is interesting to compare Fig. 6 with Fig. 1.Qualitatively, they are similar.However there are several jarring differences.First, as remarked earlier, the domain in η, previously did not even cover the current best fit of η 10 = 6.115.Second, as has been stressed in the Introduction, the relevant values of τ n are very different.Whereas in Fig. 1 values of τ n between 877 s and 937 s, were deemed reasonable.In Fig. 6, the 3σ spread in curves shown for N ν = 3, cover 876.8 s to 879.8 s.Finally, whilst the uncertainty in Y p from N ν (for ∆N ν = ±1) and τ n were comparable in Fig. 1, as one can plainly see, for the same spread in N ν , the uncertainty due a 3σ spread in τ n is almost negligible.

III. ABUNDANCE SENSITIVITIES TO τ n
In this section, we review the dependence of the light element abundances to τ n .As we have discussed earlier, we expect the abundance of 4 He to increase with increasing τ n .It is instructive and relatively straightforward to estimate this dependence.As is well known,  the dominant isotope emerging from BBN is 4 He with a mass fraction of roughly 25%.In contrast, the abundance (by number) of D and 3 He is only of order 10 −5 .Therefore to a good approximation, we can assume that after the deuterium bottleneck, all neutrons will eventually go to 4 He and we can write for the 4 He mass fraction, where the second and third expressions are evaluated at the time t BBN at the end of the deuterium bottleneck.Recalling that the neutron-to-proton ratio at freeze-out is fixed by the Boltzmann factor e −∆m/T f , when BBN begins, the n/p ratio is lowered by neutron decays so that which explicitly shows the direct τ n dependence as well as its indirect influence through the freezeout temperature T f .The latter is set by the equating the weak interaction rates Γ np ∝ τ −1 n T 5 and the Hubble rate H ∼ G 1/2 N T 2 , which gives We then can trace the effect of τ n shifts by looking at the linearized response to changes in the neutron lifetime: Then using which can be compared with the full numerical dependence (at the best fit value of η) of [5] given in Eq. ( 1).Recall that the contributions are comparable and work in the same direction: a longer τ n leads to more 4 He due to (1) a earlier weak freezeout and (2) more free neutrons surviving decay.
The numerical scaling of the light element abundances to η, N ν and τ n , is where other inputs such as nuclear cross section have been set at their mean values.These fits are normalized at the Standard Model value of N ν = 3, the best fit value of η 10 = 6.115 [15] and Eq. ( 2) for τ n .
The sensitivity of the 4 He abundance to the neutron mean life is shown in Fig. 7. Typically, because the uncertainty in the neutron mean life is small, BBN predictions for the light element abundances are obtained by combining the BBN likelihood function, L BBN (η; X i ), where the abundances X i cover D, 3 He, 4 He and 7 Li, with a CMB likelihood function taken from Planck data [7], L CMB (η, Y p ), where we include the dependence on Y p and do not assume any a priori relation between Y p and η [4,5].Throughout, we are assuming N ν = 3.
We can, however, include the explicit dependence on τ n to form the following likelihood function where L τn (τ n ) can either be derived from the Gaussian with mean and uncertainty given in Eq. ( 2) or the ideogram in Fig. 4. The projection of this likelihood function (for X i = Y p ) onto the (τ n , Y p ) plane is shown in Fig. 7. Here, we show the 1, 2, and 3σ contours from a Monte Carlo scan over the neutron mean life assuming Eq. ( 2) with a Gaussian distributed uncertainty.The star indicates the peak value of the likelihood function.The tightness of the ellipses is a consequence of the small uncertainty in τ n .This plot is an update of that in [4].The expected correlation between the neutron mean lifetime and 4 He abundance prediction is clear.It is not perfectly linear because other reaction rate uncertainties significantly contribute to the total uncertainty in 4 He as sampled in our Monte Carlo.Note that marginalizing over τ n determines the theoretical 4 He likelihood function discussed below.
If instead of taking a Gaussian distribution with mean and uncertainty given by Eq. ( 2), we use the distribution indicated by the ideogram in Fig. 4, we obtain the elongated prediction for Y p shown in Fig. 8.While the ideogram distribution should not be taken as statistically rigorous, we remind the reader that there is somewhat significant dispersion in the experimental results for τ n (a scale factor of 1.8 was already applied to obtain σ τn = 0.5 s) and the Gaussian may mask the true uncertainty. 3The tail end of the ideogram tends up to lifetimes of about 885 s, leading to significantly more 4 He. 3 Including the in-beam measurement would further increase the dispersion requiring a scale factor of 2.2.
PDG22 Weighted Average FIG. 7. The sensitivity of the 4 He abundance to the neutron mean lifetime, assuming a Gaussian distribution for τ n with mean and uncertainty given by Eq. ( 2).addition to the peak of the likelihood, denoted by a star, we show the 1, 2 and 3σ contours.
It is also common to compare the predictions of BBN to the observations by examining the theoretical and observational likelihood functions [4,5,22].The theoretical likelihood function can be expressed as a convolution of the BBN theory, dependent on η (we again fix ) and the CMB likelihood functions.The combined likelihood is defined by The second proportionality indicates that the combined BBN and CMB likelihood function is given by marginalizing Eq. ( 17) over τ n .We normalize each of the likelihood functions so that at their peaks L = 1. Figure 9 shows the comparison of these likelihood functions for (a) Y p (upper left), (b) D/H (upper right), (c) 3 He/H (lower left), and (d) 7 Li/H (lower right).In the case of 4 He, we show all three likelihood functions.The combined CMB-BBN likelihood from Eq. ( 18), L CMB−BBN (Y ), is shaded purple.The observational likelihood determined from [46,47] Y p,obs = 0.2448 ± 0.0033 .
is shaded yellow.The CMB-only likelihood, given by is shaded cyan.Given the observational and CMB uncertainties, the agreement is quite good.
The agreement between the CMB and BBN likelihoods as seen in panel b) of Fig. 9 is a major success for early Universe cosmology.We have not shown the observational likelihood function for 3 He as there are no unambiguous observations which can be associated with the primordial abundance.Similarly, it has been argued that the long-standing problem associated with 7 Li [56], is no longer well-founded [57].The main arguments for associating the observational abundance with the primordial abundance relied on the unambiguous observation of 6 Li in halo stars, as well as the lack of dispersion in 7 Li abundances at low metallicity.Both of these arguments are now suspect.We show nevertheless, the 7 Li abundance from low metallicity by the unshaded likelihood curve at 7 Li/H = (1.6 ± 0. The most striking effect between the two distributions, is the increase by a factor of 2 in the predicted uncertainty in Y p . Finally, we consider the case in which the τ n errors are dominated by systematics.In this case, it would be inappropriate to use the world average procedure we and the PDG have adopted.Using Eq. ( 1) with ∆τ n = 5 s gives ∆Y p = 0.0010 ∼ 0.3σ obs (Y p ), where σ obs (Y p ) is the observed 4 He error in eq. ( 19).We see that while this shift is large compared to the Y p theoretical error in eq. ( 22) and ( 26), it is small compared to the observed uncertainty.
Even if we consider a ∆t n = 10 s discrepancy between the in-beam measurement and UCNτ trap measurements, the 4 He change is ∆Y p = 0.0021 = 0.6σ obs (Y p ), and again the shift in the prediction is modest compare to the observed 4 He errors.The lesson is that the neutron lifetime does not impact the basic BBN concordance, but also that future improvements in both Y p and τ n will make this comparison more interesting, as we now see.

IV. BBN 'PREDICTIONS' OF τ n
It is also possible to combine the likelihood function in Eq. ( 17) with an observational likelihood function L OBS (Y p ) from Eq. ( 19) to obtain a likelihood function for τ n A similar exercise was performed in [59].This likelihood function (again normalized so that the peak is at L = 1) is shown in Fig. 11 by the blue curve.It is characterized by an approximate Gaussian given by τ n (BBN + CMB) = 870.2± 15.8s for σ Yp = 0.0033 .
In Fig. 11, we compare this BBN 'prediction' with the Gaussian given by Eq. 2 (red dotdashed curve), and that given by the ideogram (thin cyan curve).
Although we did not discuss it here, the combined BBN and CMB analysis is capable of setting strong constraints on the number of relativistic degrees of freedom at the time of BBN characterized by N ν .Currently the joint analysis gives [15] N ν = 2.989 ± 0.141 or a 2σ upper limit N ν < 3.180.The S4 CMB-only sensitivity [60] is expected to reach σ(N ν ) = 0.07.That means that if BBN is to remain competitive in the determination of N ν , improvements in 4 He observations are needed so that the uncertainty in Y p is lowered to 0.001.In that case, we obtain a better determination of τ n shown by the green dashed curve in Fig. 11.In this case, we have τ n (BBN + CMB) = 869.8± 4.8s for σ Yp = 0.001 .
This can be compared to likelihood from Eq. ( 2) shown by the dot-dashed red curve.It is no surprise that it will be very difficult for BBN + CMB to compete with direct experimental measurements of the neutron mean life.[61].We see that the two independent results are consistent within errors, and that interestingly the errorbars are almost identical.We note that unlike the lunar result, the cosmologically-preferred central value lies somewhat below the laboratory results, though the offset is not statistically significant.
Upcoming experiments and observations can improve the current situation and sharpen the link between the neutron lifetime and cosmology.This is important not only for finding the resolution of the τ n puzzle, but also for BBN.As cosmological measurements become even more precise, it remains important to improve the determination of τ n and other nuclear inputs to the BBN predictions.We look forward to progress on several fronts.
• New neutron lifetime measurements are planned.These include (a) both an upgrade magneto-gravitational trap experiment UCNτ +, and (b) an upgraded pulsed beam experiment, Beam Lifetime 3 (BL3) [62].These can shed new light on and perhaps resolve the τ n puzzle.
• The next generation CMB measurements from CMB-S4 will significantly improve both the determination of Y p and N ν from the CMB [60].Improved τ n measurements will be important for BBN to fully exploit these results, particularly N ν .
• The ongoing effort to improve astronomical Y p determinations continues.As we have discussed here and elsewhere [15], reaching the ambitious goal of σ obs (Y p ) = 0.001 would open a new window on new physics generally and τ n in particular.
As these new measurements come in, they will help deepen the links among nuclear physics, particle physics, and early Universe cosmology.

FIG. 2 .
FIG. 2. The average mean lifetime of the neutron as compiled by the Review of Particle Properties/Physics.Note that the uncertainty in the mean life for more recent measurements is smaller than the symbol showing the mean.

FIG. 3 .
FIG. 3. Ideogram for the seven measurements contributing to the PDG average neutron mean life in 2012.

Neutron
FIG. 4. Ideogram for the eight measurements contributing to the PDG average neutron mean life in 2022.

Fig. 5 .FIG. 5 .
FIG. 5. Time evolution of the light element abundances during BBN.Note that time (upper axis) increases to the right, and so the temperature is shown decreasing to the right.

10 − 2 baryon density Ω b h 2 FIG. 6 .
FIG.6.The helium mass fraction as a function of the baryon-to-photon ratio for three choices of N ν = 2, 3, and 4 and for a spread in values of τ 1/2 up to ±3σ about the mean with N ν = 3.

FIG. 8 .
FIG.8.As in Fig.7we have assumed the distribution of τ n taken from the ideogram in Fig.4.

10 10 × 7
FIG. 9. Light element abundance likelihood functions.Shown are likelihoods for each of the light nuclides.The solid-lined, dark-shaded (purple) curves are the BBN+CMB predictions, based on Planck inputs as discussed in the text.The dashed-lined, light-shaded (yellow) curves show astronomical measurements of the primordial abundances, for all but 3 He where reliable primordial abundance measures do not exist.For 4 He, the dotted-lined, medium-shaded (cyan) curve shows the independent CMB determination of 4 He.

− 3
The continued improvements in τ n measurements have played a key factor in the present robustness and precision of BBN.In the 1980s the neutron lifetime was a limiting factor in BBN predictions, and it remains the case that the τ n uncertainty dominates the uncertainty in the4 He predictions.But as we have shown, modern τ n data lead to tight predictions of the light element abundances-even given the systematic uncertainties in the neutron lifetime measurements.We have shown that combining neutron lifetime measurements as in the PDG ideogram somewhat broadens the errors in the4 He prediction, but this change remains far below the observed 4 He errors.Even a ∆τ n = 10 sec systematic uncertainty would change the Y p prediction by less than the 1σ uncertainty in the astronomical observations, and the change in the other light elements is negligible.This robustness of the BBN predictions represents a triumph of nuclear experiment, and allows BBN to play a key role in cosmology and particle physics.Indeed, the precision of BBN nuclear inputs-and that of CMB observations-allows us to turn the problem around and make a joint BBN+CMB prediction of the neutron lifetime.Our result, τ n (BBN + CMB) = 870.2± 15.8 sec (eq.31) is consistent within errors with all experimental determinations.This concordance represents a success of the big bang cosmology, but the agreement arises because the errors in the BBN+CMB predictions are too large (at present) to discriminate among the experiments.We can compare this cosmological measurement of the neutron lifetime with another recent astrophysical determination based on the NASA Lunar Prospector measurements of the neutron flux around the Moon: τ Moon n = 887 ± 14 stat +7 syst sec