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.

Coulomb screening and weak interactions in a hot, magnetized plasma are investigated. Coulomb screening is evaluated in a relativistic thermal plasma in which electrons and positrons are in equilibrium. In addition to temperature effects, effects on weak screening from a strong external magnetic field are evaluated. In high fields, the electron transverse momentum components are quantized into Landau levels. The characteristic plasma screening length at high temperatures and at high magnetic fields is explored. In addition to changes to the screening length, changes in weak interaction rates are estimated. It is found that high fields can result in increased $\beta$-decay rates as the electron and positron spectra are dominated by Landau levels. Finally, the effects studied here are evaluated in a simple r-process model. It is found that relativistic Coulomb screening has a small effect on the final abundance distribution. While changes in weak interaction rates in strong magnetic fields can have an effect on the r-process evolution and abundance distribution, the field strength required to have a significant effect may be larger than what is currently thought to be typical of the r-process environment in collapsar jets or neutron star mergers. If r-process sites exist in fields $\gtrsim 10^{14}$ G effects from fields on weak decays could be significant.

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.

We explore screening effects arising from a relativistic magnetized plasma with applications to Big Bang nucleosynthesis (BBN). %Specifically, due to their small magnetic moments, energies of electrons and positrons can be easily quantized via Landau quantization. The screening potential which depends on the thermodynamics of charged particles in the plasma is altered by the magnetic field. We focus on the impact of screening on the electron capture interaction. Taking into account the correction in BBN arising from a homogeneous primordial magnetic field (PMF), we constrain the epoch at which the PMF was generated and its strength during BBN. Considering such screening corrections to the electron capture rates and using up-to-date observations of primordial elemental abundances, we also discuss the possibility of solving the problem of under-estimation of the deuterium abundance. We find for certain values of the PMF strength predicted D and $^4$He abundances are both consistent with the observational constraints.

We calculate the abundances of $^{7}$Li, $^{11}$B, $^{92}$Nb, $^{98}$Tc, $^{138}$La, and $^{180}$Ta produced by neutrino $(\nu)$ induced reactions in a core-collapse supernova explosion. We consider the modification by $\nu$ self-interaction ($\nu$-SI) near the neutrinosphere and the Mikheyev-Smirnov-Wolfenstein effect in outer layers for time-dependent neutrino energy spectra. Abundances of $^{7}$Li and heavy isotopes $^{92}$Nb, $^{98}$Tc and $^{138}$La are reduced by a factor of 1.5-2.0 by the $\nu$-SI. In contrast, $^{11}$B is relatively insensitive to the $\nu$-SI. We find that the abundance ratio of heavy to light nucleus, $^{138}$La/$^{11}$B, is sensitive to the neutrino mass hierarchy, and the normal mass hierarchy is more likely to be consistent with the solar abundances.

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 neutrino process ($\nu$-process) for the production of 7Li and 11B in core-collapse supernovae (SNe) is extensively investigated. Initial abundances of s-nuclei and other physical conditions are derived from an updated calculation of the SN 1987A progenitor. The nuclear reaction network including neutrino reactions is constructed with the variable order Bader-Deuflhard integration method. We find that yields of 7Li and 11B significantly depend on the stellar metallicity while they are independent of the weak s-process during the stellar evolution. When the metallicity is high, there are more neutron absorbers, i.e., 56Fe, 14N (from initial CNO nuclei), and 54Fe, and the neutron abundance is small during the $\nu$-process. Since 7Be is predominantly destroyed via 7Be(n,p)7Li, a change in the neutron abundance results in different 7Be yields. Then, the calculated yield ratio 7Li/11B=0.93 for the solar metallicity is larger than that for the SN 1987A 7Li/11B=0.80 by 16 % in the inverted mass hierarchy case. We analyze contributions of respective reactions as well as abundance evolution, and clarify the $\nu$-process of 7Li and 11B.

Adopting the 3+1 neutrino mixing parameters by the IceCube and shortbase line experiments, we investigate the sterile-active neutrino oscillation effects on the supernova neutrino process. For the sterile neutrino ($\nu_s$), we study two different luminosity models. First, we presume that the $\nu_s$ does not interact with other particles through the standard interactions apart from the oscillation with the active neutrinos. Second, we consider that $\nu_s$ can be directly produced by $\nu_e$ scattering with matter. In both cases, we find that the pattern of neutrino oscillations can be changed drastically by the $\nu_s$ in supernova environments. Especially multiple resonances occur, and consequently affect thermal neutrino-induced reaction rates. As a result, $^7$Li, $^7$Be, $^{11}$B, $^{11}$C, $^{92}$Nb, $^{98}$Tc and $^{138}$La yields in the $\nu$-process are changed. Among those nuclei, $^7$Li and $^{11}$B yields can be constrained by the analysis of observed SiC X grains. Based on the meteoritic data, we conclude that the second model can be allowed while first model is excluded. The viability of the second model depends on the sterile neutrino temperature and the neutrino mass hierarchy.

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.

The chemical separation of Li+ ions induced by a magnetic field during the hierarchical structure formation can reduce initial Li abundances in cosmic structures. It is shown that cosmological reionization of neutral Li atoms quickly completes as soon as the first star is formed. Since almost all Li is singly ionized during the main course of structure formation, it can efficiently separate from gravitationally collapsing neutral gas. The separation is more efficient in smaller structures which formed earlier. In the framework of the hierarchical structure formation, extremely metal-poor stars can have smaller Li abundances because of their earlier formations. It is found that the chemical separation by a magnetic field thus provides a reason that Li abundances in extremely metal-poor stars are lower than the Spite plateau and have a large dispersion as well as an explanation of the Spite plateau itself. In addition, the chemical separation scenario can explain Li abundances in NGC 6397 which are higher than the Spite plateau. Thus, Li abundances in metal-poor stars possibly keep information on the primordial magnetic field and the structure formation history.

We reinvestigate effects of neutrino oscillations on the production of 7Li and 11B in core-collapse supernovae (SNe). During the propagation of neutrinos from the proto-neutron star, their flavors change and the neutrino reaction rates for spallation of 12C and 4He are affected. In this work corrected neutrino spallation cross sections for 4He and 12C are adopted. Initial abundances involving heavy s-nuclei and other physical conditions are derived in a new calculation of the SN 1987A progenitor in which effects of the progenitor metallicity are included. A dependence of the SN nucleosynthesis and final yields of 7Li and 11B on the neutrino mass hierarchy are shown in several stellar locations. In the normal hierarchy case, the charged current reaction rates of electron neutrinos are enhanced, and yields of proton-rich nuclei, along with 7Be and 11C, are increased. In the inverted hierarchy case, the charged current reaction rates of electron antineutrinos are enhanced, and yields of neutron-rich nuclei, along with 7Li and 11B, are increased. We find that variation of the metallicity modifies the yields of 7Li, 7Be, 11B, and 11C. This effect is caused by changes in the neutron abundance during SN nucleosynthesis. Therefore, accurate calculations of Li and B production in SNe should take into account the metallicity of progenitor stars.

We investigate the effect on the Big Bang Nucleosynthesis (BBN) from the presence of a stochastic primordial magnetic field (PMF) whose strength is spatially inhomogeneous. We assume a uniform total energy density and a gaussian distribution of field strength. In this case, domains of different temperatures exist in the BBN epoch due to variations in the local PMF. We show that in such case, the effective distribution function of particle velocities averaged over domains of different temperatures deviates from the Maxwell-Boltzmann distribution. This deviation is related to the scale invariant strength of the PMF energy density $\rho_{\rm Bc}$ and the fluctuation parameter $\sigma_{\rm B}$. We perform BBN network calculations taking into account the PMF strength distribution, and deduce the element abundances as functions of the baryon-to-photon ratio $\eta$, $\rho_{\rm Bc}$, and $\sigma_{\rm B}$. We find that the fluctuations of the PMF reduces the $^7$Be production and enhances D production. We analyze the averaged thermonuclear reaction rates compared with those of a single temperature, and find that the averaged charged-particle reaction rates are very different. Finally, we constrain the parameters $\rho_{\rm Bc}$ and $\sigma_{\rm B}$ from observed abundances of $^4$He and D, and find that the $^7$Li abundance is significantly reduced. We also find that if the $\eta$ value during BBN was larger than the present-day value due to a dissipation of the PMF or a radiative decay of exotic particles after BBN or if the stellar depletion of $^7$Li occurred, abundances of all light elements can be consistent with observational constraints.

The distribution function of the relative velocity in a two-body reaction of nonrelativistic uncorrelated particles is derived for general cases of given distribution functions of single particle velocities. The distribution function is then used in calculations of thermonuclear reaction rates. As an example, we take the Tsallis non-Maxwellian distribution, and show that the distribution function of the relative velocity is different from the Tsallis distribution. We identify an inconsistency in previous studies of nuclear reaction rates within Tsallis statistics, and derive revised nuclear reaction rates. Utilizing the revised rates, accurate results of big bang nucleosynthesis are obtained for the Tsallis statistics. For this application it is more difficult to reduce the primordial 7Li abundance while keeping other nuclear abundances within the observational constraints. A small deviation from a Maxwell-Boltzmann distribution can increase the D abundance and slightly reduce 7Li abundance. Although it is impossible to realize a 7Li abundance at the level observed in metal-poor stars, a significant decrease is possible while maintaining a consistency with the observed D abundance.

Big bang nucleosynthesis (BBN) has been used as a probe of beyond-standard physics in the early Universe, which includes a time-dependent quark mass m_q. We investigate effects of a quark mass variation delta m_q on the cross sections of the 7Be(n, p)7Li reaction and primordial light element abundances taking into account roles of 8Be resonances in the reaction during BBN. This resonant reaction has not been investigated although behaviors of low-lying resonances are not trivial. It is found that a resonance at the resonance energy E_r=0.33 MeV enhances the reaction rate and lowers the 7Li abundance significantly when the quark mass variation is negative. Based upon up-to-date observational limits on primordial abundances of D, 4He and Li, the quark mass variation in the BBN epoch are derived. In a model in which the resonance energies of the reactions 3He(d, p)4He and 3H(d, n)4He are insensitive to the quark mass, we find that the Li abundance can be consistent with observations for delta m_q/m_q = (4--8) times 10^{-3}.

Type Ia supernovae (SNe Ia) are thought to be thermonuclear explosion of white dwarfs (WDs). Their progenitors are not well understood. One popular scenario is the double-degenerate (DD) scenario, which attributes SNe Ia to WD-WD binary mergers. The fates of the WD mergers depend on the rate of $^{12}$C+$^{12}$C reaction. Recently, the $^{12}$C+$^{12}$C cross sections have been measured and the analysis of the data using the Trojan Horse Method suggested that the astrophysical reaction rate is larger than conventional rates at astrophysical temperatures due to possible resonances. The resonance contribution results in a decrease of the carbon burning ignition temperature. Therefore accretion induced collapse occurs more easily and increases the birthrate of Galactic neutron stars with the contribution of the DD scenario to the SNe Ia rate becoming even smaller.

By assuming the existence of extra-dimensional sterile neutrinos in big bang nucleosynthesis (BBN) epoch, we investigate the sterile neutrino ($\nu_{\rm s}$) effects on the BBN and constrain some parameters associated with the $\nu_{\rm s}$ properties. First, for cosmic expansion rate, we take into account effects of a five-dimensional bulk and intrinsic tension of the brane embedded in the bulk, and constrain a key parameter of the extra dimension by using the observational element abundances. Second, effects of the $\nu_{\rm s}$ traveling on or off the brane are considered. In this model, the effective mixing angle between a $\nu_{\rm s}$ and an active neutrino depends on energy, which may give rise to a resonance effect on the mixing angle. Consequently, reaction rate of the $\nu_{\rm s}$ can be drastically changed during the cosmic evolution. We estimated abundances and temperature of the $\nu_{\rm s}$ by solving the rate equation as a function of temperature until the sterile neutrino decoupling. We then find that the relic abundance of the $\nu_{\rm s}$ is drastically enhanced by the extra-dimension and maximized for a characteristic resonance energy $E_{\rm res}\gtrsim 0.01$ GeV. Finally, some constraints related to the $\nu_{\rm s}$, mixing angle and mass difference, are discussed in detail with the comparison of our BBN calculations corrected by the extra-dimensional $\nu_{\rm s}$ to observational data on light element abundances.

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.

We review important reactions in the big bang nucleosynthesis (BBN) model involving a long-lived negatively charged massive particle, $X^-$, which is much heavier than nucleons. This model can explain the observed $^7$Li abundances of metal-poor stars, and predicts a primordial $^9$Be abundance that is larger than the standard BBN prediction. In the BBN epoch, nuclei recombine with the $X^-$ particle. Because of the heavy $X^-$ mass, the atomic size of bound states $A_X$ is as small as the nuclear size. The nonresonant recombination rates are then dominated by the $d$-wave $\rightarrow$ 2P transition for $^7$Li and $^{7,9}$Be. The $^7$Be destruction occurs via a recombination with the $X^-$ followed by a proton capture, and the primordial $^7$Li abundance is reduced. Also, the $^9$Be production occurs via the recombination of $^7$Li and $X^-$ followed by deuteron capture. The initial abundance and the lifetime of the $X^-$ particles are constrained from a BBN reaction network calculation. We estimate that the derived parameter region for the $^7$Li reduction is allowed in supersymmetric or Kaluza-Klein (KK) models. We find that either the selectron, smuon, KK electron or KK muon could be candidates for the $X^-$ with $m_X\sim {\mathcal O}(1)$ TeV, while the stau and KK tau cannot.

Primordial nucleosynthesis remains as one of the pillars of modern cosmology. It is the testing ground upon which many cosmological models must ultimately rest. It is our only probe of the universe during the important radiation-dominated epoch in the first few minutes of cosmic expansion. This chapter reviews the basic equations of space-time, cosmology, and big bang nucleosynthesis. We also summarize the current state of observational constraints on primordial abundances along with the key nuclear reactions and their uncertainties. We summarize which nuclear measurements are most crucial during the big bang. We also review various cosmological models and their constraints. In particular, we analyze the constraints that big bang nucleosynthesis places upon the possible time variation of fundamental constants, along with constraints on the nature and origin of dark matter and dark energy, long-lived supersymmetric particles, gravity waves, and the primordial magnetic field.

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.

The modified gravity is considered to be one of possible explanations of the accelerated expansions of the present and the early universe. We study effects of the modified gravity on big bang nucleosynthesis (BBN). If effects of the modified gravity are significant during the BBN epoch, they should be observed as changes of primordial light element abundances. We assume a $f(G)$ term with the Gauss-Bonnet term $G$, during the BBN epoch. A power-law relation of $df/dG \propto t^p$ where $t$ is the cosmic time was assumed for the function $f(G)$ as an example case. We solve time evolutions of physical variables during BBN in the $f(G)$ gravity model numerically, and analyzed calculated results. It is found that a proper solution for the cosmic expansion rate can be lost in some parameter region. In addition, we show that calculated results of primordial light element abundances can be significantly different from observational data. Especially, observational limits on primordial D abundance leads to the strongest constraint on the $f(G)$ gravity. We then derive constraints on parameters of the $f(G)$ gravity taking into account the existence of the solution of expansion rate and final light element abundances.

Big bang nucleosynthesis in a modified gravity model of $f(R)\propto R^n$ is investigated. The only free parameter of the model is a power-law index $n$. We find cosmological solutions in a parameter region of $1< n \leq (4+\sqrt{6})/5$. We calculate abundances of $^4$He, D, $^3$He, $^7$Li, and $^6$Li during big bang nucleosynthesis. We compare the results with the latest observational data. It is then found that the power-law index is constrained to be $(n-1)=(-0.86\pm 1.19)\times 10^{-4}$ (95 % C.L.) mainly from observations of deuterium abundance as well as $^4$He abundance.

During the structure formation, charged and neutral chemical species may have separated from each other at the gravitational contraction in primordial magnetic field (PMF). A gradient in the PMF in a direction perpendicular to the field direction leads to the Lorentz force on the charged species. Resultantly, an ambipolar diffusion occurs, and charged species can move differently from neutral species, which collapses gravitationally during the structure formation. We assume a gravitational contraction of neutral matter in a spherically symmetric structure, and calculate fluid motions of charged and neutral species. It is shown that the charged fluid, i.e., proton, electron and $^7$Li$^+$, can significantly decouple from the neutral fluid depending on the field amplitude. The charged species can, therefore, escape from the gravitational collapse. We take the structure mass, the epoch of the gravitational collapse, and the comoving Lorenz force as parameters. We then identify a parameter region for an effective chemical separation. This type of chemical separation can reduce the abundance ratio of Li/H in early structures because of inefficient contraction of $^7$Li$^+$ ion. Therefore, it may explain Li abundances of Galactic metal-poor stars which are smaller than the prediction in standard big bang nucleosynthesis model. Amplitudes of the PMFs are controlled by a magneto-hydrodynamic turbulence. The upper limit on the field amplitude derived from the turbulence effect is close to the value required for the chemical separation.

The primordial lithium abundance inferred from observations of metal-poor stars is ~3 times smaller than the theoretical value in standard big bang nucleosynthesis (BBN) model. We assume a simple model including a sterile neutrino nu_H with mass of O(10) MeV which decays long after BBN. We then investigate cosmological effects of a sterile neutrino decay. We formulate the injection spectrum of nonthermal photons induced by electrons and positrons generated at the nu_H decay, as a function of the nu_H mass and the photon temperature. We then consistently solve (1) the cosmic thermal history, (2) nonthermal nucleosynthesis induced by the nonthermal photons, (3) the baryon-to-photon ratio eta, and (4) the effective neutrino number N_eff. Amounts of energy injection at the nu_H decay are constrained from limits on primordial D and 7Li abundances, the N_eff value, and the cosmic microwave background energy spectrum. We find that 7Be is photodisintegrated and the Li problem is partially solved for the lifetime 10^4-10^5 s and the mass >~ 14 MeV. 7Be destruction by more than a factor of 3 is not possible because of an associated D over-destruction. In the parameter region, the eta value is decreased slightly, while the N_eff value is increased by a factor of <~ 1. In this study, errors in photodisintegration cross sections of 7Be(g, a)3He and 7Li(g, a)3H that have propagated through literatures are corrected. It is then found that the new photodisintegration rates are 2.3 to 2.5 times smaller than the old rates, so that efficiencies of 7Be and 7Li photodisintegration are significantly smaller.

The $^7$Li abundance calculated in BBN with the baryon-to-photon ratio fixed from fits to the CMB power spectrum is inconsistent with the observed lithium abundances on the surface of metal-poor halo stars. Previous cosmological solutions proposed to resolve this $^7$Li problem include photon cooling (possibly via the Bose-Einstein condensation of a scalar particle) or the decay of a long-lived $X-$particle (possibly the next-to-lightest supersymmetric particle). In this paper we reanalyze these solutions, both separately and in concert. We also introduce the possibility of a primordial magnetic field (PMF) into these models. We constrain the $X-$particles and the PMF parameters by the observed light element abundances using a likelihood analysis to show that the inclusion of all three possibilities leads to an optimum solution to the lithium problem. We deduce allowed ranges for the $X-$particle parameters and energy density in the PMF that can solve $^7$Li problem.

We extensively reanalyze effects of a long-lived negatively charged massive particle, X-, on big bang nucleosynthesis (BBN). The BBN model with an X- particle was originally motivated by the discrepancy between 6,7Li abundances predicted in standard BBN model and those inferred from observations of metal-poor stars. In this model 7Be is destroyed via the recombination with an X- particle followed by radiative proton capture. We calculate precise rates for the radiative recombinations of 7Be, 7Li, 9Be, and 4He with X-. In nonresonant rates we take into account respective partial waves of scattering states and respective bound states. The finite sizes of nuclear charge distributions cause deviations in wave functions from those of point-charge nuclei. For a heavy X- mass, m_X \gtrsim 100 GeV, the d-wave --> 2P transition is most important for 7Li and 7,9Be, unlike recombination with electrons. Our new nonresonant rate of the 7Be recombination for m_X=1000 GeV is more than 6 times larger than the existing rate. Moreover, we suggest a new important reaction for 9Be production: the recombination of 7Li and X- followed by deuteron capture. We derive binding energies of X-nuclei along with reaction rates and Q-values. We then calculate BBN and find that the amount of 7Be destruction depends significantly on the charge distribution of 7Be. Finally, updated constraints on the initial abundance and the lifetime of the X- are derived in the context of revised upper limits to the primordial 6Li abundance. Parameter regions for the solution to the 7Li problem are revised, and the primordial 9Be abundances is revised.

The injections of energetic hadrons could have occurred in the early universe by decays of hypothetical long-lived exotic particles. The injections induce the showers of nonthermal hadrons via nuclear scattering. Neutrons generated at these events can react with $^7$Be nuclei and reduce $^7$Be abundance solving a problem of the primordial $^7$Li abundance. We suggest that thermal neutron injection is a way to derive a model independent conservative limit on the relation between abundances of D and $^7$Li in a hadronic energy injection model. We emphasize that an uncertainty in cross sections of inelastic $n+p$ scattering affects the total number of induced neutrons, which determines final abundances of D and $^7$Li. In addition, the annihilations of antinucleons with $^4$He result in higher D abundance and trigger nonthermal $^6$Li production. It is concluded that a reduction of $^7$Li abundance from a value in the standard big bang nucleosynthesis (BBN) model down to an observational two $\sigma$ upper limit is necessarily accompanied by an undesirable increase of D abundance up to at least an observational 12 $\sigma$ upper limit from observations of quasi-stellar object absorption line systems. The effects of antinucleons and secondary particles produced in the hadronic showers always lead to a severer constraint. The BBN models involving any injections of extra neutrons are thus unlikely to reproduce a small $^7$Li abundance consistent with observations.

The existence of an exotic long-lived negatively charged massive particle, i.e., X-, during big bang nucleosynthesis can affect primordial light element abundances. Especially, the final abundance of 7Li, mainly originating from the electron capture of 7Be, has been suggested to reduce by the 7Be destruction via the radiative X- capture of 7Be followed by the radiative proton capture of the bound state of 7Be and X- (7BeX). We suggest a new route of 7BeX formation, that is the 7Be charge exchange at the reaction of 7Be3+ ion and X-. The formation rate depends on the number fraction of 7Be3+ ion, the charge exchange cross section of 7Be3+, and the probability that produced excited states 7BeX* are converted to the ground state. We estimate respective quantities affecting the 7BeX formation rate, and find that this reaction pathway can be more important than ordinary radiative recombination of 7Be and X-. The effect of the charge exchange reaction is then shown in a latest nuclear reaction network calculation. Quantum physical model calculations for related reactions are needed to precisely estimate the efficiency of this pathway in future.

Using the primordial helium abundance, an upper limit to the magnetic moments for Dirac neutrinos had been provided by imposing restrictions on the number of the additional helicity states. Considering non-thermal photons produced in the decay of the heavy sterile mass eigenstates due to the neutrino magnetic moment, we explore the constraints imposed by the observed abundances of all the light elements produced during the Big Bang nucleosynthesis.

Spectroscopic observations of metal-poor stars have indicated possible 6Li abundances that are much larger than the primordial abundance predicted in the standard big bang nucleosynthesis model. Possible mechanisms of 6Li production in metal-poor stars include pregalactic and cosmological cosmic ray (CR) nucleosynthesis and nucleosynthesis by flare-accelerated nuclides. We study 9Be production via two-step alpha-fusion reactions of CR or flare-accelerated 3,4He through 6He and 6,7Li, in pregalactic structure, intergalactic medium, and stellar surfaces. We solve transfer equations of CR or flare particles and calculate nuclear yields of 6He, 6,7Li, and 9Be taking account of probabilities of processing 6He and 6,7Li into 9Be via fusions with alpha particles. Yield ratios, i.e., 9Be/6Li, are then calculated for the CR and flare nucleosynthesis models. We suggest that the future observations of 9Be in metal-poor stars may find enhanced abundances originating from metal-poor CR or flare activities.

Big bang nucleosynthesis (BBN) is affected by the energy density of a primordial magnetic field (PMF). For an easy derivation of constraints on models for PMF generations, we assume a PMF with a power law (PL) distribution in wave number defined with a field strength, a PL index, and maximum and minimum scales at a generation epoch. We then show a relation between PL-PMF parameters and the scale invariant (SI) strength of PMF for the first time. We perform a BBN calculation including PMF effects, and show abundances as a function of baryon to photon ratio $\eta$. The SI strength of the PMF is constrained from observational constraints on abundances of $^4$He and D. The minimum abundance of $^7$Li/H as a function of $\eta$ slightly moves to a higher $^7$Li/H value at a larger $\eta$ value when a PMF exists during BBN. We then discuss degeneracies between the PL-PMF parameters in the PMF effect. In addition, we assume a general case in which both the existence and the dissipation of PMF are possible. It is then found that an upper limit on the SI strength of the PMF can be derived from a constraint on $^4$He abundance, and that a lower limit on the allowed $^7$Li abundance is significantly higher than those observed in metal-poor stars.

A new upper limit on the amplitude of primordial magnetic field (PMF) is derived by a comparison between a calculation of elemental abundances in big bang nucleosynthesis (BBN) model and the latest observational constraints on the abundances. Updated nuclear reaction rates are adopted in the calculation. Effects of PMF on the abundances are consistently taken into account in the numerical calculation with the precise formulation of changes in physical variables. We find that abundances of 3He and 6Li increase while that of 7Li decreases when the PMF amplitude increases, in the case of the baryon-to-photon ratio determined from the measurement of cosmic microwave background radiation. We derive a constraint on the present amplitude of PMF, i.e., B(0)<1.5 micro G [corresponding to the amplitude less than 2.0x10^{11} G at BBN temperature of T=10^9 K] based on the rigorous calculation.

Following a recent suggestion of axion cooling of photons between the nucleosynthesis and recombination epochs in the Early Universe, we investigate a hybrid model with both axions and relic supersymmetric particles. In this model we demonstrate that the 7Li abundance can be consistent with observations without destroying the important concordance of deuterium abundance.

Hypothetical long-lived massive colored particles (MCPs or Ys) would be confined in colorless exotic strongly interacting massive particles (SIMPs) at color confinement temperature of T_C /sim 200 MeV. Two long-lived MCPs form a bound state (Y /bar{Y}) at collisions of two SIMPs. We study sensitivities of MCP annihilation to decay properties of resonances (Y /bar{Y}), and binding energies or energy levels of exotic SIMPs. The (Y /bar{Y}) formation is assumed to dominantly proceeds through resonances of (Y/bar{Y}) in this paper. We make a toy model of the effective cross section for Y /bar{Y} annihilation. Abundances of SIMPs are then calculated for different sets of parameters specifying properties of (Y /bar{Y}) resonances, binding energies of SIMPs, the initial abundance and the mass of MCP. Calculated relic abundances for respective SIMP species are 2 x 10^{-8}--3 x 10^{-4} times that of baryon. They can be much higher but cannot be much smaller than the previous estimate. The abundances can be consistent depending on parameters with the possible scenario that SIMPs bind to nuclei and subsequent exotic nuclear reactions reduce the primordial abundance of 7Li or enhance those of 9Be and/or B in the early universe. A unique information on the quark-hadron phase transition in the early universe may become available in future by elaborated studies on the annihilation process with light element abundances as observables.

We study the excitation of fine structure levels of C I, C II and O I by ultraviolet (UV) photons around strong UV sources which are also ionizing sources of the cosmological reionization at redshift of $\sim 10$. The evolutions of ionized regions around a point source are calculated by solving rate equations for non-equilibrium chemistry. Signals of UV photons through the fine structure lines are considered to be stronger at locations of more abundant chemical species of C I, C II and O I. Such environments would be realized where strong fluxes of non-ionizing UV line photons available for the pumping up of fine structure levels exist, and simultaneously ionizing UV photons are effectively shielded by dense H I regions. Signals from H I regions of moderately large densities induced by redshifted UV photons emitted at the point sources are found to be dominantly large over those of others. We discuss the detectability of the signals, and show that signals from idealized environments will be possibly detected by radio observations with next-generation arrays to come after the Atacama Large Millimeter/submillimeter Array (ALMA).

We reinvestigate the constraints from primordial nucleosynthesis on a possible time-dependent quark mass. The limits on such quark-mass variations are particularly sensitive to the adopted observational abundance constraints. Hence, in the present study we have considered updated light-element abundances and uncertainties deduced from observations. We also consider new nuclear reaction rates and an independent analysis of the influence of such quark-mass variations on the resonance properties of the important 3He(d,p)4He reaction. We find that the updated abundance and resonance constraints imply a narrower range on the possible quark-mass variations in the early universe. We also find that, contrary to previous investigations, the optimum concordance region reduces to a (95% C.L.) value of -0.005 < delta m_q/m_q < 0.007 consistent with no variation in the averaged quark mass.

We identify reactions which destroy 7Be and 7Li during big bang nucleosynthesis (BBN) in the scenario of BBN catalyzed by a long-lived sub-strongly interacting massive particle (sub-SIMP or X particle). The destruction associated with non radiative X captures of the nuclei can be realized only if the interaction strength between an X particle and a nucleon is properly weaker than that between two nucleons to a degree depending on the mass of X. Binding energies of nuclei to an X particle are estimated taking the mass and the interaction strength to nuclei of the X as input parameters. Nuclear reaction rates associated with the X are estimated naively, and adopted in calculating evolutions of nuclear abundances. We suggest that the 7Li problem, which might be associated with as-yet-unrecognized particle processes operating during BBN, can be solved if the X particle interacts with nuclei strongly enough to drive 7Be destruction but not strongly enough to form a bound state with 4He of relative angular momentum L=1. Justifications of this scenario by rigorous calculations of reaction rates using quantum mechanical many-body models are highly desirable since this result involves many significant uncertainties.

Tabulated rates for astrophysical photodisintegration reactions make use of Boltzmann statistics for the photons involved as well as the interacting nuclei. Here we derive analytic corrections for the Planck-spectrum quantum statistics of the photon energy distribution. These corrections can be deduced directly from the detailed-balance condition without the assumption of equilibrium as long as the photons are represented by a Planck spectrum. Moreover we show that these corrections affect not only the photodisintegration rates but also modify the conditions of nuclear statistical equilibrium as represented in the Saha equation. We deduce new analytic corrections to the classical Maxwell-Boltzmann statistics which can easily be added to the reverse reaction rates of existing reaction network tabulations. We show that the effects of quantum statistics, though generally quite small, always tend to speed up photodisintegration rates and are largest for nuclei and environments for which Q/kT ~ 1. As an illustration, we examine possible effects of these corrections on the r-process, the rp-process, explosive silicon burning, the $\gamma$-process and big bang nucleosynthesis. We find that in most cases one is quite justified in neglecting these corrections. The correction is largest for reactions near the drip line for an r-process with very high neutron density, or an rp-process at high-temperature.

We calculate nucleosynthesis of proton-rich isotopes in the carbon-deflagration model for Type Ia supernovae (SNe Ia). The seed abundances are obtained by calculating the s-process nucleosynthesis that is expected to occur in the repeating helium shell flashes on the carbon-oxygen (CO) white dwarf (WD) during mass accretion from a binary companion. When the deflagration wave passes through the outer layer of the CO WD, p-nuclei are produced by photodisintegration reactions on s-nuclei in a region with the peak temperature ranging from 1.9 to 3.6 x10^9 K. We confirm the sensitivity of the p-process on the initial distribution of s-nuclei. We show that the initial C/O ratio in the WD does not much affect the yield of p-nuclei. On the other hand, the abundance of 22Ne left after the s-processing has a large influence on the p-process via 22Ne(alpha,n) reaction. We find that about 50% of p-nuclides are co-produced when normalized to their solar abundances in all adopted cases of seed distribution. Mo and Ru, which are largely underproduced in Type II supernovae (SNe II), are produced more than in SNe II although they are underproduced with respect to the yield levels of other p-nuclides. The ratios between p-nuclei and iron in the ejecta are larger than the solar ratios by a factor of 1.2. We also compare the yields of oxygen, iron, and p-nuclides in SNe Ia and SNe II and suggest that SNe Ia could make a larger contribution than SNe II to the solar system content of p-nuclei.

We study effects of relic long-lived strongly interacting massive particles (X particles) on big bang nucleosynthesis (BBN). The X particle is assumed to have existed during the BBN epoch, but decayed long before detected. The interaction strength between an X and a nucleon is assumed to be similar to that between nucleons. Rates of nuclear reactions and beta decay of X-nuclei are calculated, and the BBN in the presence of neutral charged X^0 particles is calculated taking account of captures of X^0 by nuclei. As a result, the X^0 particles form bound states with normal nuclei during a relatively early epoch of BBN leading to the production of heavy elements. Constraints on the abundance of X^0 are derived from observations of primordial light element abundances. Particle models which predict long-lived colored particles with lifetimes longer than about 200 s are rejected. This scenario prefers the production of 9Be and 10B. There might, therefore, remain a signature of the X particle on primordial abundances of those elements. Possible signatures left on light element abundances expected in four different models are summarized.

It has been proposed that the apparent discrepancies between the inferred primordial abundances of 6Li and 7Li and the predictions of big bang nucleosynthesis (BBN) can be resolved by the existence of a negatively-charged massive unstable supersymmetric particle (X-) during the BBN epoch. Here, we present new BBN calculations with an X- particle utilizing an improved nuclear reaction network including captures of nuclei by the particle, nuclear reactions and beta-decays of normal nuclei and nuclei bound to the X- particles (X-nuclei), and new reaction rates derived from recent rigorous quantum many-body dynamical calculations. We find that this is still a viable model to explain the observed 6Li and 7Li abundances. However, contrary to previous results, neutral X-nuclei cannot significantly affect the BBN light-element abundances. We also show that with the new rates the production of heavier nuclei is suppressed and there is no signature on abundances of nuclei heavier than Be in the X--particle catalyzed BBN model as has been previously proposed. We also consider the version of this model whereby the X- particle decays into the present cold dark matter. We analyze the this paradigm in light of the recent constraints on the dark-matter mass deduced from the possible detected events in the CDMS-II experiment. We conclude that based upon the inferred range for the dark-matter mass, only X- decay via the weak interaction can achieve the desired 7Li destruction while also reproducing the observed 6Li abundance.

A recent measurement of $^4$He photodisintegration reactions, $^4$He($\gamma$,$p$)$^3$H and $^4$He($\gamma$,$n$)$^3$He with laser-Compton photons shows smaller cross sections than those estimated by other previous experiments at $E_\gamma \lesssim 30$ MeV. We study big-bang nucleosynthesis with the radiative particle decay using the new photodisintegration cross sections of $^4$He as well as previous data. The sensitivity of the yields of all light elements D, T, $^3$He, $^4$He, $^6$Li, $^7$Li and $^7$Be to the cross sections is investigated. The change of the cross sections has an influence on the non-thermal yields of D, $^3$He and $^4$He. On the other hand, the non-thermal $^6$Li production is not sensitive to the change of the cross sections at this low energy, since the non-thermal secondary synthesis of $^6$Li needs energetic photons of $E_\gamma \gtrsim 50$ MeV. The non-thermal nucleosynthesis triggered by the radiative particle decay is one of candidates of the production mechanism of $^6$Li observed in metal-poor halo stars (MPHSs). In the parameter region of the radiative particle lifetime and the emitted photon energy which satisfies the $^6$Li production above the abundance level observed in MPHSs, the change of the photodisintegration cross sections at $E_\gamma \lesssim 30$ MeV as measured in the recent experiment leads to $\sim 10$% reduction of resulting $^3$He abundance, whereas the $^6$Li abundance does not change for this change of the cross sections of $^4$He($\gamma$,$p$)$^3$H and $^4$He($\gamma$,$n$)$^3$He. The $^6$Li abundance, however, could show a sizable change and therefore the future precise measurement of the cross sections at high energy $E_\gamma \gtrsim$ 50 MeV is highly required.

The $^6$Li abundance observed in metal poor halo stars appears to exhibit a plateau as a function of metallicity similar to that for $^7$Li, suggesting a big bang origin. However, the inferred primordial abundance of $^6$Li is $\sim$1000 times larger than that predicted by standard big bang nucleosynthesis for the baryon-to-photon ratio inferred from the WMAP data. Also, the inferred $^7$Li primordial abundance is 3 times smaller than the big bang prediction. We here describe in detail a possible simultaneous solution to both the problems of underproduction of $^6$Li and overproduction of $^7$Li in big bang nucleosynthesis. This solution involves a hypothetical massive, negatively-charged leptonic particle that would bind to the light nuclei produced in big bang nucleosynthesis, but would decay long before it could be detected. We consider only the $X$-nuclear reactions and assume that the effect of decay products is negligible, as would be the case if lifetime were large or the mass difference between the charged particle and its daughter were small. An interesting feature of this paradigm is that, because the particle remains bound to the existing nuclei after the cessation of the usual big bang nuclear reactions, a second longer epoch of nucleosynthesis can occur among $X$-nuclei. We confirm that reactions in which the hypothetical particle is transferred can occur that greatly enhance the production of $^6$Li while depleting $^7$Li. We also identify a new reaction that destroys large amounts of $^7$Be, and hence reduces the ultimate $^7$Li abundance. Thus, big-bang nucleosynthesis in the presence of these hypothetical particles, together with or without an event of stellar processing, can simultaneously solve the two Li abundance problems.

I calculate the evolution of Be and B abundances produced by cosmic rays generated by massive stars in the pregalactic phase of the universe. The inputs for calculation, i.e. the star formation rate and the nuclear abundances of cosmic rays, which I assume to be the same as those of the ISM, are taken from the results of a detailed cosmic chemical evolution model with its parameters best fitted from several items of observational information including an early reionization of the IGM by $z\sim 15$. I found that when the $^6$Li plateau abundance observed in metal-poor halo stars originated in the pregalactic cosmological cosmic ray nucleosynthesis, Be and B simultaneously produced with $^6$Li amount to the lowest levels ever detected in metal-poor halo stars. It is desirable to observe Be and B abundances in metal-poor halo stars with [Fe/H]$\leq -3$ in order to elucidate the possibility of early $^6$LiBeB production by pregalactic supernova cosmic ray nucleosynthesis.

The $^6$Li abundance observed in metal poor halo stars exhibits a plateau similar to that for $^7$Li suggesting a primordial origin. However, the observed abundance of $^6$Li is a factor of $10^3$ larger and that of $^7$Li is a factor of 3 lower than the abundances predicted in the standard big bang when the baryon-to-photon ratio is fixed by WMAP. Here we show that both of these abundance anomalies can be explained by the existence of a long-lived massive, negatively-charged leptonic particle during nucleosynthesis. Such particles would capture onto the synthesized nuclei thereby reducing the reaction Coulomb barriers and opening new transfer reaction possibilities, and catalyzing a second round of big bang nucleosynthesis. This novel solution to both of the Li problems can be achieved with or without the additional effects of stellar destruction.

Recent spectroscopic observations of metal poor stars have indicated that both $^7$Li and $^6$Li have abundance plateaus with respect to the metallicity. Abundances of $^7$Li are about a factor three lower than the primordial abundance predicted by standard big-bang nucleosynthesis (SBBN), and $^6$Li abundances are $\sim 1/20$ of $^7$Li, whereas SBBN predicts negligible amounts of $^6$Li compared to the detected level. These discrepancies suggest that $^6$Li has another cosmological or Galactic origin than the SBBN. Furthermore, it would appear that $^7$Li (and also $^6$Li) has been depleted from its primordial abundance by some post-BBN processes. In this paper we study the possibility that the radiative decay of long-lived particles has affected the cosmological lithium abundances. We calculate the non-thermal nucleosynthesis associated with the radiative decay and explore the allowed region of the parameters specifying the properties of long-lived particles. We also impose constraints from observations of the CMB energy spectrum. We find that non-thermal nucleosynthesis produces $^6$Li at the level detected in metal poor halo stars (MPHSs), when the lifetime of the unstable particles is of the order $\sim 10^8-10^{12}$~s and their initial abundance with respect to that of the photons is $\sim (10^{-13}-10^{-12}~{\rm GeV})/E_{\gamma 0}$, where $E_{\gamma 0}$ is the emitted photon energy in the radiative decay. We conclude that the most likely nucleosynthetic scenario involves two different processes. First, a non-thermal cosmological nucleosynthesis associated with the radiative decay of unstable particles; and second, the stellar depletion of both of the primordial lithium isotopic abundances.