Big-bang nucleosynthesis (BBN) is valuable as a means to constrain the physics of the early universe and it is the only probe of the radiation-dominated epoch. A fundamental assumption in BBN is that the nuclear velocity distributions obey Maxwell-Boltzmann (MB) statistics as they do in stars. Specifically, the BBN epoch is characterized by a dilute baryon plasma for which the velocity distribution of nuclei is mainly determined by the dominant Coulomb elastic scattering with mildly relativistic electrons. One must therefore deduce the momentum distribution for reacting nuclei from the multi-component relativistic Boltzmann equation. However, the full multi-component relativistic Boltzmann equation has only recently been analyzed and its solution has only been worked out in special cases. Moreover, a variety of schemes have been proposed that introduce non-thermal components into the BBN environment which can alter the thermal distribution of reacting nuclei. Here, we construct the relativistic Boltzmann equation in the context of BBN. We also derive a Langevin model and perform relativistic Monte-Carlo simulations which clarify the baryon distribution during BBN and can be used to analyze any relaxation from a non-thermal injection. We show by these analyses that the thermalization process leads to a nuclear distribution function that remains very close to MB statistics even during the most relativistic environment relevant to BBN. Hence, the predictions of standard BBN remain unchanged.

We report on the development and application of a multi dimensional relativistic Monte-Carlo code to explore the thermalization process in a relativistic multi-component environment. As an illustration we simulate the fully relativistic three dimensional Brownian-motion-like solution to the thermalization of a high mass particle in a bath of relativistic low-mass particles. We follow the thermalization and ultimate equilibrium distribution of the Brownian-like particle as can happen in the cosmic plasma during Big bang nucleosynthesis.

This paper is superseded by Arxiv:1911.07334. Big-bang nucleosynthesis (BBN) is a valuable tool to constrain the physics of the early universe and is the only probe of the radiation-dominated epoch. A fundamental assumption in BBN is that the nuclear velocity distributions obey Maxwell-Boltzmann statistics as they do in stars. In this letter, however, we point out that there is a fundamental difference between stellar reaction rates and BBN reaction rates. Specifically, the BBN epoch is characterized by a dilute baryon plasma for which the velocity distribution of nuclei is mainly determined by the dominant Coulomb scattering with mildly relativistic electrons. This modifies the nuclear velocity distributions and significantly alters the thermonuclear reaction rates, and hence, the light-element abundances. We show that this novel result alters all previous calculations of light-element abundances from BBN, and indeed exacerbates the discrepancies between BBN and inferred primordial light-element abundances possibly suggesting the need for new physics in the early universe.

The abundance of primordial lithium is derived from the observed spectroscopy of metal-poor stars in the galactic halo. However, the observationally inferred abundance remains at about a factor of three below the abundance predicted by standard big bang nucleosynthesis (BBN). The resolution of this dilemma can be either astrophysical (stars destroy lithium after BBN), nuclear (reactions destroy lithium during BBN), or cosmological, i.e. new physics beyond the standard BBN is responsible for destroying lithium. Here, we overview a variety of possible cosmological solutions, and their shortcomings. On the one hand, we examine the possibility of physical processes that modify the velocity distribution of particles from the usually assumed Maxwell-Boltzmann statistics. A physical justification for this is an inhomogeneous spatial distribution of domains of primordial magnetic field strength as a means to reduce the primordial lithium abundance. Another possibility is that scattering with the mildly relativistic electrons in the background plasma alters the baryon distribution to one resembling a Fermi-Dirac distribution. We show that neither of these possibilities can adequately resolve the lithium problem. A number of alternate hybrid models are discussed including a mix of neutrino degeneracy, unified dark matter, axion cooling, and the presence of decaying and/or charged supersymmetric particles.

We summarize some applications of big bang nucleosythesis (BBN) and the cosmic microwave background (CMB) to constrain the first moments of the creation of matter in the universe. We review the basic elements of BBN and how it constraints physics of the radiation-dominated epoch. In particular, how the existence of higher dimensions impacts the cosmic expansion through the projection of curvature from the higher dimension in the "dark radiation" term. We summarize current constraints from BBN and the CMB on this brane-world dark radiation term. At the same time, the existence of extra dimensions during the earlier inflation impacts the tensor to scalar ratio and the running spectral index as measured in the CMB. We summarize how the constraints on inflation shift when embedded in higher dimensions. Finally, one expects that the universe was born out of a complicated multiverse landscape near the Planck time. In these moments the energy scale of superstrings was obtainable during the early moments of chaotic inflation. We summarize the quest for cosmological evidence of the birth of space-time out of the string theory landscape. We will explore the possibility that a superstring excitations may have made itself known via a coupling to the field of inflation. This may have left an imprint of "dips" in the power spectrum of temperature fluctuations in the cosmic microwave background. The identification of this particle as a superstring is possible because there may be evidence for different oscillator states of the same superstring that appear on different scales on the sky. It will be shown that from this imprint one can deduce the mass, number of oscillations, and coupling constant for the superstring. Although the evidence is marginal, this may constitute the first observation of a superstring in Nature.

Big Bang Nucleosynthesis (BBN) explores the first few minutes of nuclei formation after the Big Bang. We present updates that result in new constraints at the 2{\sigma} level for the abundances of the four primary light nuclides - D,3He,4He, and 7Li - in BBN. A modified standard BBN code was used in a Monte Carlo analysis of the nucleosynthesis uncertainty as a function of baryon-to-photon ratio. Reaction rates were updated to those of NACRE and REACLIB, and R-Matrix calculations. The results are then used to derive a new constraint on the effective number of neutrinos.

The term dark radiation is used both to describe a noninteracting neutrino species and as a correction to the Friedmann Equation in the simplest five-dimensional RS-II brane-world cosmology. In this paper we consider the constraints on both meanings of dark radiation based upon the newest results for light-element nuclear reaction rates, observed light-element abundances and the power spectrum of the Cosmic Microwave Background (CMB). Adding dark radiation during big bang nucleosynthesis (BBN) alters the Friedmann expansion rate causing the nuclear reactions to freeze out at a different temperature. This changes the final light element abundances at the end of BBN. Its influence on the CMB is to change the effective expansion rate at the surface of last scattering. We find that the BBN constraint reduces the allowed range for both types of dark radiation at 10 Mev to between $-12.1\%$ and $+6.2\%$ of the {\bf total} background energy density at 10 Mev. Combining this result with fits to the CMB power spectrum, produces different results for particle vs. brane-world dark radiation. In the brane-world, the range decreases to $-6.0\%$ to $+6.2\%$. Thus, we find, that the ratio of dark radiation to the background total relativistic mass energy density $\rho_{\rm DR}/\rho$ is consistent with zero although there remains a very slight preference for a positive (rather than negative) contribution.

A dark radiation term arises as a correction to the energy momentum tensor in the simplest five-dimensional RS-II brane-world cosmology. In this paper we revisit the constraints on dark radiation based upon the newest results for light-element nuclear reaction rates, observed light-element abundances and the power spectrum of the Cosmic Microwave Background (CMB). Adding dark radiation during big bang nucleosynthesis alters the Friedmann expansion rate causing the nuclear reactions to freeze out at a different temperature. This changes the final light element abundances at the end of BBN. Its influence on the CMB is to change the effective expansion rate at the surface of last scattering. We find that our adopted BBN constraints reduce the allowed range for dark radiation to between $-12.1\%$ and $+6.2\%$ of the ambient background energy density. Combining this result with fits to the CMB power spectrum, the range decreases to $-6.0\%$ to $+6.2\%$. Thus, we find, that the ratio of dark radiation to the background total relativistic mass energy density $\rho_{\rm DR}/\rho$ is consistent with zero although in the BBN analysis there could be a slight preference for a negative contribution. However, the BBN constraint depends strongly upon the adopted primordial helium abundance.