Improved BBN constraints on the variation of the gravitational constant

Big Bang Nucleosynthesis (BBN) is very sensitive to the cosmological expansion rate. If the gravitational constant G took a different value during the nucleosynthesis epoch than today, the primordial abundances of light elements would be affected. In this work, we improve the bounds on this variation using recent determinations of the primordial element abundances, updated nuclear and weak reaction rates and observations of the Cosmic Microwave Background (CMB). When combining the measured abundances and the baryon density from CMB observations by Planck, we find GBBN/G0=0.99-0.05+0.06\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$G_\mathrm {BBN}/G_0 = 0.99^{+0.06}_{-0.05}$$\end{document} at 2σ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$2\sigma $$\end{document} confidence level. If the variation of G is linear in time, we find G˙/G0=0.7-4.3+3.8×10-12year-1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\dot{G}/G_0 = 0.7^{+3.8}_{-4.3}\times 10^{-12} \, \mathrm {year}^{-1}$$\end{document}, again at 2σ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$2\sigma $$\end{document}. These bounds are significantly stronger than those from previous primordial nucleosynthesis studies, and are comparable and complementary to CMB, stellar, solar system, lunar laser ranging, pulsar timing and gravitational wave constraints.


Introduction
Ever since the Large Number Hypothesis of Dirac, physicists have wondered whether the constants of nature may evolve over time [1]. Since then the mathematical framework for such a variation has been developed. This was done first with ideas such as Kaluza-Klein theory where the coupling between gravity and matter sectors was set by the size of a compact dimension [2,3], and then by advancements such as Brans-Dicke theory where the coupling between matter and gravity is endowed with dynamics [4]. A key problem in the ongoing attempt to unify the four forces of nature using string theory is figuring out how to stabilise the higher a e-mail: james.alvey@kcl.ac.uk b e-mail: nashwan.sabti@kcl.ac.uk (corresponding author) c e-mail: miguel.escudero@kcl.ac.uk d e-mail: malcolm.fairbairn@kcl.ac.uk dimensions, and consequently the value of the low energy couplings, including the gravitational constant G [5].
A theory where G is allowed to vary almost always involves (by definition) promoting G to be related to the expectation value of some dynamical scalar field. For consistency, the scalar field will then have a kinetic term and a potential. In order to observe any dynamical changes of such a field over cosmological time, the curvature of that potential, and therefore the mass of this field, has to be very small. Constraints from solar system observations, like radar ranging of Mars, place very tight limits on the contribution to gravity due to such light scalars in the Universe today [6]. 1 On the other hand, we are in possession of much less information about the detailed behaviour of gravity in the early Universe, in particular during the period of Big Bang Nucleosynthesis at which point the Universe is expected to be extremely uniform. Most scenarios where the value of G is set by a scalar field obtaining an expectation value result in an effective mass for that field that is much larger than the energy scales relevant even at this earlier epoch. As such, we do not have any good reason to expect that G is likely to be different. Nevertheless, if we have the technology to place new limits on the variation of a fundamental constant like this one, we should take the opportunity to obtain those limits to the best of our abilities. The subject may also become more important in the coming years given upcoming atomic interferometry experiments which will test extensions of gravity [8].
Strictly speaking, since G is dimensionful, we should be very careful when we discuss its time variation. This is because any physically meaningful change in couplings should be written as a change in the dimensionless ratio between two dimensionful quantities. We are assuming a scenario where G = M −2 Pl changes while none of the particle physics energy scales vary. For example, consider the case where all the dimensionful parameters of the Standard Model (like the Brout-Englert-Higgs expectation value v and QCD ) are determined by dimensional transmutation of couplings set at a single energy scale M GUT . If gravitational corrections to the running of those couplings are irrelevant, then discussing variations in M Pl /M GUT becomes better defined [9]. This is the kind of situation we are considering. In what follows, we will refer to G as the value of Newton's gravitational constant over time and G 0 as its value as measured today.
Big Bang Nucleosynthesis (BBN) is sensitive to modifications of the expansion history in the early Universe. Indeed, predictions for the primordial element abundances are strongly dependent on the Hubble rate H and subsequently on the value of the gravitational constant G during the relevant epoch. Recent improvements in measurements of these early Universe abundances [10][11][12][13][14][15][16] allow us to improve constraints on the gravitational physics during the BBN era. With some mild assumptions detailed in the Sect. 3, this can be translated into bounds on the variation in the value of G between nucleosynthesis and today.
It is the nucleosynthesis bounds which we look to update in this work. We do this by (i) using up-to-date measurements of the primordial element abundances reported by the PDG [38], (ii) including a weak determination of the baryon density [17] to alleviate the degeneracy between b h 2 and G in the deuterium abundance, (iii) accurately accounting for incomplete neutrino decoupling following [39,40], and (iv) making use of the state-of-the-art BBN code PRIMAT [41] which has updated nuclear reaction rates and accounts for many corrections to the weak reaction rates. The methodology applied here follows that in [42].

Cosmological implications
We can understand the effect on the primordial abundances of a different value of G during nucleosynthesis, G BBN = G 0 , in terms of the Hubble expansion rate H ∝ √ G. The dominant effect of an increased (decreased) expansion rate is to alter the time at which various weak and nuclear processes freeze-out. In particular, the proton-to-neutron conversion processes and the p + n ↔ D + γ reaction will freezeout earlier (later). This leads to an over(under)-production of both helium and deuterium compared to the case where G BBN = G 0 . This is seen clearly in Fig. 1. From just these figures and the errors in the measurements of the abundances (indicated by the grey bands), we expect that our analysis will be able to constrain variations in G below the 10% level.

Methodology
To derive the bounds on the variation of the gravitational constant, we follow the approach presented in [42]. We make use of the publicly available code NUDEC_BSM [39,40] to compute the background cosmology, including the effects of non-instantaneous neutrino decoupling. The relevant cosmological parameters are subsequently forwarded to the state-of-the-art BBN code PRIMAT [41], which takes care of the nuclear reaction network and time evolution of primordial abundances. In all simulations we take the neutron lifetime to be the default value in PRIMAT: τ n = 879.5 s. The abundances scale with G/G 0 in a way that closely matches the semi-analytic relations given in [43].
To quantify the effect of a time variation in G we choose to parametrise the evolution as a slowly evolving linear function of time t [35], Here t 0 is the current age of the Universe such that (t 0 − t BBN ) 13.8 Gyr [44]. This slowly varying function ensures that the assumption of a constant G during the cosmologically very short period of primordial nucleosynthesis is an excellent approximation. In the Sect. 5, we will quote bounds on both G BBN /G 0 anḋ G/G 0 , where the latter can be derived from Eq. (1).

Data analysis
We make use of the observed primordial abundances of helium and deuterium as reported by the PDG [38]. At 68% CL, these are: D/H| P = (2.569 ± 0.027) × 10 −5 .
We also include theoretical errors in the predictions of Y P and D/H| P due to uncertainties in the various nuclear reaction rates and the neutron lifetime [43]; σ (D/H| P ) Theo = 0.13 × 10 −5 .
To quantify deviations from the measured primordial abundances due to changes in the gravitational constant, we construct a χ 2 for BBN as follows: We are also interested in including a conservative determination of the baryon density to lift the degeneracy between b h 2 and G. Using directly the posterior values from the baseline CDM Planck 2018 analysis will not be satisfactory, because G is kept constant there. Instead, we use the results of [17], who carry out a Planck likelihood analysis including variations in G CMB . In Table I of [17], they find that the mean baryon density exactly matches that of the base Planck 2018 TTTEEE + lowE analysis within CDM [44], albeit with twice as large error bars: This allows us to define an extended χ 2 for BBN+ b h 2 , For both the pure BBN analysis and the extended scenario, we compute the relevant χ 2 on a grid of ( b h 2 , G BBN /G 0 ). We then marginalize over the baryon density to find a 1-D χ 2 (G BBN /G 0 ). To rule out values of G BBN /G 0 , we compare the computed statistics to critical values of the 1-D χ 2 distribution. In particular, at 2σ , we rule out a scenario if

Results
We show the cosmological implications of a variation in the gravitational constant on the primordial helium and deuterium abundances in Fig. 1. It is evident that measurements of both primordial abundances are able to strongly constrain a deviation of G from today's value. Moreover, we see that the impact of a higher value of G BBN /G 0 can be compensated by a higher value of the baryon density b h 2 , which introduces a degeneracy in the ( b h 2 , G BBN /G 0 ) plane. This is because changes in the primordial deuterium abundance are linearly proportional to changes in the baryon density ( D/H| P ∝ − b h 2 ), while the primordial helium abundance is only logarithmically dependent on b h 2 [41].
In Fig. 2 we show the 1σ and 2σ confidence intervals in the ( b h 2 , G BBN /G 0 ) plane and χ 2 as a function of the variation in G for the two benchmark analyses considered here. At 95.4% CL, we obtain: Alternatively, assuming a linear time evolution of the gravitational constant as described in Eq. (1), these bounds can be translated into a constraint onĠ/G 0 , which at 95.4% CL reads: Moreover, the BBN + b h 2 analysis highly disfavours a fractional deviation of G larger than ∼10% (at more than 5σ ), while a similar conclusion can be drawn in the BBNonly analysis for deviations larger than ∼20%. This difference arises because a restriction on b h 2 lifts the degeneracy with G regarding the primordial deuterium abundance. Note that our constraints improve on the previous 1σ primordial nucleosynthesis bounds [33,34] by approximately a factor of 10.
As a final comment, other references, such as [41], quote an error on the deuterium value that is approximately 4 times smaller than the one used here. Running our analysis with this lower error, we find no change in the BBN only bounds while the BBN + b h 2 constraints strengthen by ∼30%.

Conclusions
Big Bang Nucleosynthesis is sensitive to changes in the value of the gravitational constant in the early Universe. By using current measurements of the primordial abundances of helium and deuterium we have shown that at 95.4% CL, G BBN /G 0 = 0.98 +0.06 −0.06 and G BBN /G 0 = 0.99 +0.06 −0.05 if b h 2 measurements from Planck are also accounted for. Assuming a very slow linear time evolution of G, these constraints map into a bound on the time variation at 95.4% CL ofĠ/G 0 = 1.4 +4.4 −4.7 × 10 −12 yr −1 andĠ/G 0 = 0.7 +3.8 −4.3 × 10 −12 yr −1 respectively. These constraints are competitive and complementary to those from CMB, stellar, solar system, pulsar timing, gravitational waves and lunar laser ranging measurements.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: This article is based on data presented in Ref. [38].] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .