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 carry out multifractal analyses of multiple tracers namely the main galaxy sample, the LRG sample and the quasar sample from the SDSS to test the assumption of cosmic homogeneity and identify the scale of transition to homogeneity, if any. We consider the behaviour of the scaled number counts and the scaling relations of different moments of the galaxy number counts in spheres of varying radius $R$ to calculate the spectrum of the Minkowski-Bouligand general dimension $D_{q} (R)$ for $-4 \leq q \leq 4$. The present analysis provides us the opportunity to study the spectrum of the generalized dimension $D_{q}(R)$ for multiple tracers of the cosmic density field over a wide range of length scales and allows us to confidently test the validity of the assumption of cosmic homogeneity. Our analysis indicates that the SDSS main galaxy sample is homogeneous on a length scales of $80\, h^{-1}\, {\rm Mpc} $ and beyond whereas the SDSS quasar sample and the SDSS LRG sample show transition to homogeneity on an even larger length scales at $\sim 150\, h^{-1}\, {\rm Mpc}$ and $\sim 230\, h^{-1}\, {\rm Mpc}$ respectively. These differences in the scale of homogeneity arise due to the effective mass and redshift scales probed by the different tracers in a Universe where structures form hierarchically. Our results reaffirm the validity of cosmic homogeneity on large scales irrespective of the tracers used and strengthens the foundations of the Standard Model of Cosmology.