Nucleosynthesis Predictions and High-Precision Deuterium Measurements

: Two new high-precision measurements of the deuterium abundance from absorbers along the line of sight to the quasar PKS1937–1009 were presented. The absorbers have lower neutral hydrogen column densities (N(HI) ≈ 18 cm − 2 ) than for previous high-precision measurements, boding well for further extensions of the sample due to the plenitude of low column density absorbers. The total high-precision sample now consists of 12 measurements with a weighted average deuterium abundance of D/H = 2.55 ± 0.02 × 10 − 5 . The sample does not favour a dipole similar to the one detected for the ﬁne structure constant. The increased precision also calls for improved nucleosynthesis predictions. For that purpose we have updated the public AlterBBN code including new reactions, updated nuclear reaction rates, and the possibility of adding new physics such as dark matter. The standard Big Bang Nucleosynthesis prediction of D/H = 2.456 ± 0.057 × 10 − 5 is consistent with the observed value within 1.7 standard deviations.


Introduction
Consisting only of a proton and a neutron, deuterium is the simplest element, and consequently it was the first one to form in the early universe. The nucleosynthesis of the light elements such as deuterium, helium and lithium took place during the first few minutes after Big Bang when the universe was still opaque, and today we can only probe the conditions indirectly through the resulting abundances. This requires high precision measurements combined with reliable predictions for the abundances. In Section 2 we provide a short overview of the nucleosynthesis, before we present new high precision measurements of the deuterium abundance in Section 3 and a new code for performing predictions in Section 4.

Background
In order to form bound nuclei it is necessary to have free neutrons and protons with temperatures below their binding energy. The decoupling of the weak interaction and consequently the appearance of free neutrons and protons happens already around one second after Big Bang, but the binding energy of deuterium is as low as 2.2 MeV. Since deuterium is the simplest element you can form, the consequence of the low binding energy is a delay of the nucleosynthesis until 2-3 min after Big Bang. Meanwhile the free neutrons decay, altering the final abundances of deuterium and helium.
During Big Bang Nucleosynthesis only light elements are formed because helium is very strongly bound (binding energy of 28.3 MeV) and there exist no stable nucleus with five nucleons. This means you cannot fuse helium and neutrons/protons into heavier elements, providing a severe bottleneck for the formation of heavier elements. Due to the decay of free neutrons, only the ones that become bound in nuclei survive. Most of the produced deuterium is converted to helium, leaving only of the order of 10 −5 deuterium nuclei per hydrogen nuclei.
The formation of the light elements is sensitive to the expansion rate (temperature) and timing of the events e.g., decoupling, and consequently the abundances of the light elements can probe the conditions during the nucleosynthesis.

Observational Deuterium Measurements
Neutral deuterium has absorption features similar to the Lyman series of hydrogen, but with a 82 km · s −1 off-set because of the heavier nucleus. Consequently we can measure the "fingerprint" of deuterium caused by absorption in gas clouds along the line of sight to distant bright sources such as quasars. However, due to the fairly small off-set and low abundance of deuterium relative to hydrogen (D/H = 10 −5 ), high resolution spectra (R ≈ 45000) with very good signal to noise are required. Given that there are no significant astrophysical sources of deuterium production [1,2] and the destruction rate in stars is low at the relevant redshifts and metallicities [3,4], medium-to-high redshift measurements will be very close to the primordial value.
Here we present two new high precision measurements along the line of sight to the quasar PKS1937-101. The absorber at z abs = 3.572 was previously analysed by Tytler et al. [5], Burles and Tytler [6], and the z abs = 3.256 absorber was analysed by Crighton et al. [7]. New spectra of the quasar have become available in the form of 63.4 ks exposures from the High Resolution Echelle Spectrometer (HIRES) at the Keck Telescope and 32.4 ks from the Ultraviolet and Visual Echelle Spectrograph (UVES) at the Very Large Telescope (VLT). A root-mean-square weighted stacking of the new spectra is shown in Figure 1. The spectra are reduced and normalised following standard procedures and the relevant atomic transitions are fitted with a model consisting of several Voigt-profiles that represent the individual sub-components of the absorbers. Each sub-component is characterised by a redshift, a column density, a temperature and a turbulent broadening. In addition, we allow for uncertainties in the continuum normalisation and calibration of the spectra. The details of the spectra and their analyses can be found in Riemer-Sørensen et al. [8,9], and here we only discuss the results.

Deuterium in the z abs = 3.256 Absorber
We fit the absorption with Voigt profiles using VPFIT [10] to obtain the gas cloud characteristics including the deuterium and hydrogen abundances.
We find a deuterium abundance of D/H = 2.45 ± 0.28 × 10 −5 to be compared with the previous measurement of D/H = 1.6 ± 0.30 × 10 −5 [7]. If naively interpreted as a statistical fluctuation this corresponds to a 3σ difference. However, it origins in a systematical model difference, which is hard to quantify. The measurement from the new spectra is not more precise despite the increased signal-to-noise ratio, but we believe it to be more robust, as the better data allow us to resolve the internal structure of the absorber better. This is illustrated in Figure 2. The best fit model is consistent between the UVES and HIRES spectra. Figure 2. The normalised stacked flux spectrum (thick black) of the z abs = 3.256 absorber, and the best fit model (green), as well as the normalised fit residuals (above each fit). The stacked spectrum is only used for visualisation, as the model fitting is always performed on the individual spectra simultaneously. The fitted velocity components are marked by vertical dot-dashed lines denoted by A to F (HI in orange, DI in blue, metals in cyan), and interloping HI lines are shown as vertical dotted lines (light grey). It is seen that the structure of the absorber is more complicated than the old model from Crighton et al. [7] (orange), which has a significant impact on the final result. Figure from

Deuterium in the z abs = 3.572 Absorber
For the high redshift absorber we find a ratio of D/H = 2.62 ± 0.03 × 10 −5 to be compared with the previous measurement of D/H = 3.3 ± 0.3 × 10 −5 [6]. In contrast to Burles and Tytler [6] we include the heavy element transitions in the model fitting, resulting in a more complicated substructure for the absorber as shown in Figure 3. The interesting point here is the very high precision we obtain, even though the hydrogen column density of log(N(H I)) = 17.9 is lower than for previous similar high-precision measurements [11]. This is important for future measurements because the neutral hydrogen column density distribution in quasar absorption systems is a steep power law, with lower column density systems being more common. A statistically large sample of measurements is therefore feasible and at the same time strictly necessary in order to reveal a plateau of primordial values as a function of, e.g., metallicity. Despite the complexity of the absorber, the very high precision on the D/H ratio is obtained due to a combination of (1) very high signal to noise ratio in the spectra; (2) fitting 9 Lyman transitions for which the oscillator strengths span more than two orders of magnitude; and (3) the presence of several heavy element species that allow us to constrain the velocity structure of the system. Figure 3. The normalised stacked flux spectrum (thick black) of the z abs = 3.572 absorber, and the best fit model (green), as well as the normalised fit residuals (above each fit). The stacked spectrum is only used for visualisation, as the model fitting is always performed on the individual spectra simultaneously. The fitted velocity components are marked by vertical dot-dashed lines denoted by A to D (HI in red, DI i in blue, metals in cyan), and interloping HI lines are shown as vertical dotted lines (light grey). The model from Tytler et al. [5] (without blends) is over-plotted in orange, illustrating the difference between the derived model structures, particularly for the metals.

Deuterium Measurement Sample
For standard Big Bang Nucleosynthesis the inferred value of the primordial deuterium abundance can be used to infer the baryon density of the universe, Ω b . Table 1 provide a list of current measurements, for which the weighted average is 100Ω b h 2 = 2.17 ± 0.024 which deviates by 1.9σ from the value of 2.225 ± 0.016 inferred from the measurements of the cosmic microwave background [12]. We consider the two results to be consistent.
We find no correlations between deuterium abundance and redshift, metallicity or hydrogen column density (as tentatively claimed by Cooke et al. [13]). Neither do we find any evidence for a dipole in the deuterium measurements [9]. Such a dipole would be expected if there were variations in fundamental constants like the fine structure constant, the hadronic masses or binding energies [14]. As an example we consider a dipole with the same direction as the potential dipole in the fine structure constant [15,16], for which the preferred slopes are close to zero with uncertainties larger than the preferred value and consequently consistent both with a small dipole and with no dipole.

Predictions from Nucleosynthesis
The improved observational precision calls for nucleosynthesis predictions with similar accuracy. We have updated and expanded the existing AlterBBN code with nuclear rates and non-standard physics options. AlterBBN is a publicly available C code for evaluating abundances [27], which is based on the Wagoner code [28] and similar to NUC123 ( [29], also known as the Kawano code) and PArthENoPE [30] 1 . AlterBBN can be run as an independent code or included in multi-parameter analyses using e.g., CosmoMC 2 or Montepython 3 . Until the improvements described here are implemented in the official version of AlterBBN, the code and further details can be obtained from github.com/espensem/AlterBBN. AlterBBN evaluates the light element abundances as a function of time during the Big Bang Nucleosynthesis (BBN). It starts from a thermodynamic equilibrium, whereafter the weak interaction and the neutrinos decouple, leading to the decay of free neutrons. Eventually the electrons and positrons will decouple and annihilate, which will reheat the photons (and any other particles with an electromagnetic coupling). Until the photon temperature has dropped below the binding energy of deuterium, the protons and neutrons cannot form any deuterium.
The resulting abundances depend on the nuclear reaction rates and the timing/expansion in the early universe. Without any non-standard physics, the baryon to photon ratio is the only free parameter in BBN. Since it can be determined very precisely from the cosmic microwave background [12], the main uncertainty in the predictions comes from the nuclear reaction rates.

Updating Nuclear Reaction Rates
We have extended the nuclear network in AlterBBN from 88 to 100 reactions, and updated six important rates with new theoretical/experimental determinations [26]. The resulting helium 1 PArthENoPE is also publicly available, but it requires expensive fortran libraries to compile 2 http://cosmologist.info/cosmomc 3 http://baudren.github.io/montepython.html and hydrogen-3 abundances are relatively unaffected by these changes, but D/H change by 4.5%, and the uncertainty double. Lithium-7 and beryllium-7 increase by 18% thereby worsening the lithium problem.
An example of the updated rates is the formation of helium-3 via d(p,γ) 3 He. For this rate there is only very scarce experimental data at the relevant energy range [31]. Instead the rates have been determined theoretically or from fitting various polynomials or theoretical models to an extended energy range and extrapolating to the range of interest. The newest theoretical calculation by Marcucci et al. [32] is inconsistent with the experimental rate measured around 0.1 MeV. To reconcile the discrepancy Coc et al. [26] rescaled the theoretical model from Viviani et al. [33], Marcucci et al. [34] to the experimental rates for the full range of 0.002-2 Mev. Using AlterBBN with the theoretically determined rate from [32] leads to D/H = 2.49 ± 0.03 ± 0.03 × 10 −5 (the two errors are due to nuclear rate uncertainties and uncertainties in the baryon density, respectively) while the experimental rescaling from [26] leads to D/H = 2.45 ± 0.057 × 10 −5 . The two results deviate by approximately 1σ.
Another option is to use the experimental data directly (as in e.g., [35,36]) 4 , which will further increase the deviation between the various predictions from AlterBBN. At this point we do not want to advocate one choice over another, but simply want to emphasise that the choice of rates matters, and unless the specific choice is discussed in the context of each result (which is usually not the case (e.g., [8,11,25])), this deviation between rates should be included as a systematic uncertainty. Better experimental determination of the rates at BBN energies will solve the issue.

Implementing New Physics
Non-standard physics such as extra relativistic species, modified gravity or lepton asymmetry may change the expansion rate in the early universe or the order of the events governing the nucleosynthesis [37][38][39]. These effects could potentially origin in production/destruction of dark matter or neutrino-like species. At present and, likely for the foreseeable future, BBN provides the only window to a universal lepton asymmetry. The most obvious way to generate lepton asymmetry is via additional neutrino species We have updated AlterBBN with the possibility of adding the presence of equivalent neutrinos (non-interacting relativistic species) and Weakly Interacting Massive Particles (WIMPs) during the nucleosynthesis. The WIMPs are generic and their effect is purely determined by their mass, the number of internal degrees of freedom, their nature (fermion or boson) as well as any coupling to the standard model particles. The WIMPs can be either electromagnetically coupled or coupled to neutrinos (standard model neutrinos and any equivalent species) (following [38,39]). They are not required to provide 100% of the dark matter but could be a subcomponent. We consider four types of WIMPs: (1) Real scalars that are self-conjugate and have one degree of freedom; (2) Complex scalars that have two degrees of freedom; (3) Majorana fermions which are self-conjugate and have two degrees of freedom; and (4) Dirac fermions with four degrees of freedom.
The presence of such non-standard model particles will alter the expansion rate because they contribute to the total energy density. They might also increase the temperature of photons or neutrinos due to annihilation, and by mimicking neutrinos, they might change the inferred number of relativistic species often parametrised as a change in the number of neutrino species N eff = 3.046 + ∆N ν where ∆N ν is the equivalent number of neutrinos that would provide the same effect (e.g., [40,41]). The relativistic particles contribute to the radiation energy density so that ρ rad → ρ rad = ρ rad + ∆N ν ρ ν . Combining with the Friedman equation [42] we can relate the expansion factor directly to N eff as 4 To be implemented in the next version of AlterBBN where H /H is the relative expansion rate with and without the extra species. Similarly, lepton asymmetry can be parametrised in terms of N eff (see e.g., [37]), and an increase in entropy (e.g., from decaying dark matter) will affect the number of degrees of freedom, which can also be expressed as a change in N eff . In Figure 4 we show the probability contours for simultaneously varying the baryon density (η = n b /n γ ) and the additional number of neutrinos for otherwise standard BBN and comparing to observational constraints on D I/H I and 4 He. As seen the constraints from BBN alone are consistent with those from the Planck Surveyor and of comparable precision. The combined D I/H I and 4 He constraint on ∆N ν is ∆N ν = 0.06 +0. 27 −0.36 .  Including WIMPs will change the radiation and/or matter energy density and pressure directly, as well as possible changes in the entropy and lepton asymmetry depending on the specific properties. The expansion rate governs the particle number densities and consequently the interaction rates for nucleosynthesis. Figure 5 illustrate the changes of the D/H ratio when including WIMPs in the mass range from 0.01 to 100 MeV (heavier or lighter WIMPs will provide linear extensions of the abundance in Figure 5 as they decouple completely before/after the nucleosynthesis). The results are compared to those of [38,39]. The deviations are due to different choices of nuclear rates (e.g., the [38,39] results are based on the theoretical calculation of d(p,γ) 3 He from [33] while we have used the experimentally rescaled rate from [26]). Comparing to the observational constraints, we see that for electromagnetically coupled WIMPs we can exclude masses below 10 MeV and for neutrino coupled WIMPs only masses around 10 MeV are allowed. The left plot is for electromagnetically coupled WIMPs (coupling to electrons or photons) and the right plot is for neutrino coupling. The grey lines are for comparison to [38,39]. The deviations are due to different choices of nuclear rates. The yellow horisontal bands are the observed ratio of D/H [8]. Figure from [44].
Neither non-standard physics nor the new rates alleviate the Lithium problem.

Summary
We have presented new high precision measurement of the deuterium to hydrogen ratio from low column density absorbers (N(H I) = 17.9 and 18.1 cm −2 ) along the line of sight to the quasar PKS1937-101. This is important for future measurements because the neutral hydrogen column density distribution in quasar absorption systems is a steep power law, with lower column density systems being more common. A statistically large sample of measurements is therefore feasible.
With the acquired precision on the observational measurements, it is timely to similarly improve the theoretical predictions from Big Bang nucleosynthesis. We have updated the publicly available code, AlterBBN, with new reaction rates and non-standard physics. It is clear that the choice between theoretical and measured rates at relevant energies leads to different results, providing a source of systematic uncertainty which is often neglected.
Supplementary Materials: The updates of AlterBBN is available here github.com/espensem/AlterBBN and will be included in the next official and distributed release of AlterBBN.