Previous studies of the reaction ${}^{3}\mathrm{He} +{}^{4}\mathrm{He} \rightarrow {}^{7}\mathrm{Be} + \gamma$ have focused on providing the best central value and error bar for the $S$ factor at solar energies. Measurements of this capture reaction, the ${}^{3}\mathrm{He}$-${}^{4}\mathrm{He}$ scattering phase shifts, as well as properties of ${}^{7}\mathrm{Be}$, have been used to constrain employed theoretical models. Here we show that much more information than was previously appreciated can be extracted from angle-integrated capture data alone. We use the next-to-leading-order (NLO) amplitude in an effective field theory (EFT) for the reaction to perform the extrapolation. At this order the EFT describes the reaction using an s-wave scattering length and effective range, the asymptotic properties of the final bound states, and short-distance contributions to the $E1$ capture amplitude. We extract the multi-dimensional posterior of all these parameters via a Bayesian analysis. We find that properties of the ${}^{7}\mathrm{Be}$ ground and excited states are well constrained. The total $S$ factor $S(0)= 0.578^{+0.015}_{-0.016}$ keV~b, while the branching ratio for excited- to ground-state capture at zero energy, $Br(0)=0.406^{+0.013}_{-0.011}$, both at 68\% degree of belief. This $S(0)$ is broadly consistent with other recent evaluations, including the previously recommended value $S(0)=0.56 \pm 0.03$ eV b, but has a smaller error bar. We also find significant constraints on the scattering parameters, and we obtain constraints on the $S(E)$'s angular dependence. The path forward seems to lie with better measurements of the scattering phase shift and $S(E)$'s angular dependence, together with better understanding of the asymptotic normalization coefficients of the ${}^7$Be bound states' wave functions. Data on these could further reduce $S(0)$'s uncertainty.

We employ an effective field theory (EFT) that exploits the separation of scales in the p-wave halo nucleus $^8\mathrm{B}$ to describe the process $^7\mathrm{Be}(p,\gamma)^8\mathrm{B}$ up to a center-of-mass energy of 500 keV. The calculation, for which we develop the lagrangian and power counting, is carried out up to next-to-leading order (NLO) in the EFT expansion. The power counting we adopt implies that Coulomb interactions must be included to all orders in $\alpha_{\rm em}$. We do this via EFT Feynman diagrams computed in time-ordered perturbation theory, and so recover existing quantum-mechanical technology such as the two-potential formalism for the treatment of the Coulomb-nuclear interference. Meanwhile the strong interactions and the E1 operator are dealt with via EFT expansions in powers of momenta, with a breakdown scale set by the size of the ${}^7$Be core, $\Lambda \approx 70$ MeV. Up to NLO the relevant physics in the different channels that enter the radiative capture reaction is encoded in ten different EFT couplings. The result is a model-independent parametrization for the reaction amplitude in the energy regime of interest. To show the connection to previous results we fix the EFT couplings using results from a number of potential model and microscopic calculations in the literature. Each of these models corresponds to a particular point in the space of EFTs. The EFT structure therefore provides a very general way to quantify the model uncertainty in calculations of $^7\mathrm{Be}(p,\gamma)^8\mathrm{B}$. We also demonstrate that the only N$^2$LO corrections in $^7\mathrm{Be}(p,\gamma)^8\mathrm{B}$ come from an inelasticity that is practically of N$^3$LO size in the energy range of interest, and so the truncation error in our calculation is effectively N$^3$LO. We also discuss the relation of our extrapolated $S(0)$ to the previous standard evaluation.

We report the discovery and analysis of the most metal-poor damped Lyman-alpha (DLA) system currently known, which also displays the Lyman series absorption lines of neutral deuterium. The average [O/H] abundance of this system is [O/H] = -2.804 +/- 0.015, which includes an absorption component with [O/H] = -3.07 +/- 0.03. Despite the unfortunate blending of many weak D I absorption lines, we report a precise measurement of the deuterium abundance of this system. Using the six highest quality and self-consistently analyzed measures of D/H in DLAs, we report tentative evidence for a subtle decrease of D/H with increasing metallicity. This trend must be confirmed with future high precision D/H measurements spanning a range of metallicity. A weighted mean of these six independent measures provides our best estimate of the primordial abundance of deuterium, 10^5 (D/H)_P = 2.547 +/- 0.033 (log_10 (D/H)_P = -4.5940 +/- 0.0056). We perform a series of detailed Monte Carlo calculations of Big Bang nucleosynthesis (BBN) that incorporate the latest determinations of several key nuclear cross sections, and propagate their associated uncertainty. Combining our measurement of (D/H)_P with these BBN calculations yields an estimate of the cosmic baryon density, 100 Omega_B,0 h^2(BBN) = 2.156 +/- 0.020, if we adopt the most recent theoretical determination of the d(p,gamma)3He reaction rate. This measure of Omega_B,0 h^2 differs by ~2.3 sigma from the Standard Model value estimated from the Planck observations of the cosmic microwave background. Using instead a d(p,gamma)3He reaction rate that is based on the best available experimental cross section data, we estimate 100 Omega_B,0 h^2(BBN) = 2.260 +/- 0.034, which is in somewhat better agreement with the Planck value. Forthcoming measurements of the crucial d(p,gamma)3He cross section may shed further light on this discrepancy.

We have studied the 7Be(p,photon)8B reaction in the Halo effective field theory (EFT) framework. The leading order (LO) results were published in Phys.Rev.C89,051602(2014) after the isospin mirror process, 7Li(n,photon)8Li, was addressed in Phys.Rev.C89,024613(2014). In both calculations, one key step was using the final shallow bound state asymptotic normalization coefficients (ANCs) computed by ab initio methods to fix the EFT couplings. Recently we have developed the next-to-LO (NLO) formalism (to appear soon), which could reproduce other model results by no worse than 1% when the 7Be-p energy was between 0 and 0.5 MeV. In our recent report (arXiv:1507.07239), a different approach from that in Phys.Rev.C89,051602(2014) was used. We applied Bayesian analysis to constrain all the NLO-EFT parameters based on measured S-factors, and found tight constraints on the S-factor at solar energies. Our S(E=0 MeV)= 21.3 + - 0.7 eV b. The uncertainty is half of that previously recommended. In this proceeding, we provide extra details of the Bayesian analysis, including the computed EFT parameters' probability distribution functions (PDFs) and how the choice of input data impacts final results.

We report an improved low-energy extrapolation of the cross section for the process Beryllium-7+proton -> Boron-8+photon, which determines the Boron-8 neutrino flux from the Sun. Our extrapolant is derived from Halo Effective Field Theory (EFT) at next-to-leading order. We apply Bayesian methods to determine the EFT parameters and the low-energy S-factor, using measured cross sections and scattering lengths as inputs. Asymptotic normalization coefficients of Boron-8 are tightly constrained by existing radiative capture data, and contributions to the cross section beyond external direct capture are detected in the data at E < 0.5 MeV. Most importantly, the S-factor at zero energy is constrained to be S(0)= 21.3 + - 0.7 eV b, which is an uncertainty smaller by a factor of two than previously recommended. That recommendation was based on the full range for S(0) obtained among a discrete set of models judged to be reasonable. In contrast, Halo EFT subsumes all models into a controlled low-energy approximant, where they are characterized by nine parameters at next-to-leading order. These are fit to data, and marginalized over via Monte Carlo integration to produce the improved prediction for S(E).

(abridged) In the presence of a light WIMP (mass m_chi < 30 MeV), there are degeneracies among the WIMP's nature, its couplings to standard model particles, its mass, and the number of equivalent (additional) neutrinos, Delta N_nu. These degeneracies cannot be broken by the CMB constraint on the effective number of neutrinos, N_eff. However, since big bang nucleosynthesis (BBN) is also affected by a light WIMP and equivalent neutrinos, complementary BBN and CMB constraints can break some of the degeneracy. In a previous paper BBN and CMB were combined to explore allowed ranges for m_chi, Delta N_nu, and N_eff for light WIMPs that annihilate electromagnetically (EM) to photons and/or electrons/positrons. In this paper BBN predictions with a light WIMP that only couples to neutrinos are calculated. Recent observed abundances of ^2H and ^4He are used to limit m_chi, Delta N_nu, N_eff, and the present-day baryon density. Allowing for a neutrino coupled light WIMP and nonzero Delta N_nu, combined BBN and CMB data give lower limits to m_chi, with a best fit m_chi > 35 MeV, equivalent to no light WIMP at all. All masses below 4--9 MeV (depending on spin) are excluded. Without any light WIMP, BBN alone prefers Delta N_nu = 0.50 +- 0.23, favoring neither Delta N_nu = 0, nor a fully thermalized sterile neutrino (Delta N_nu = 1). This result is consistent with the CMB constraint, N_eff = 3.30 +- 0.27, limiting "new physics" between BBN and recombination. Combining BBN and CMB data gives Delta N_nu = 0.35 +- 0.16 and N_eff = 3.40 +- 0.16; while BBN and the CMB combined require Delta N_nu > 0 at ~98% confidence, they disfavor Delta N_nu > 1 at > 99% confidence. Allowing a neutrino-coupled light WIMP extends the allowed range slightly downward for Delta N_nu and slightly upward for N_eff simultaneously, leaving best-fit values unchanged.

(Abridged) In the presence of a light WIMP (< 30 MeV), there are degeneracies among the nature of the WIMP, its couplings to standard model particles, its mass m_chi, and the number of equivalent neutrinos beyond the standard model, Delta N_nu. These degeneracies cannot be broken by the CMB constraint on the effective number of neutrinos, N_eff. However, big bang nucleosynthesis (BBN) is affected by a light WIMP and equivalent neutrinos, so the combination of BBN and CMB can break some of the degeneracies. Here, BBN predictions for abundances in the presence of a light WIMP and equivalent neutrinos are explored, and estimates of their observationally determined relic abundances are used to limit m_chi, Delta N_nu, and the present Universe baryon density (Omega_B h^2). These constraints are explored for Majorana and Dirac fermion WIMPs, as well as for real and complex scalar WIMPs that couple to electrons, positrons, and photons. In a separate paper this analysis is repeated for WIMPs that couple only to the standard model neutrinos, and constraints for the two cases are contrasted. Without a light WIMP but allowing for equivalent neutrinos, the combined BBN and CMB constraints favor N_eff = 3.46 +/- 0.17, Omega_B h^2 = 0.0224 +/- 0.0003, and Delta N_nu = 0.40 +/- 0.17. In this case, standard BBN (Delta N_nu = 0) is disfavored at ~98% confidence, and the presence of one sterile neutrino (Delta N_nu = 1) is disfavored at > 99% confidence. Allowing a light WIMP and Delta N_nu equivalent neutrinos together, the combined BBN and CMB data provide lower limits to the WIMP masses (> 0.5 - 5 MeV) that depend on the nature of the WIMP, favor m_chi ~8 MeV slightly over no WIMP, and loosen constraints on equivalent neutrinos, Delta N_nu = 0.65+0.46-0.35. While Delta N_nu = 0 is still disfavored at ~95% confidence with a light WIMP, Delta N_nu = 1 is now allowed.

Very light WIMPs (chi), thermal relics that annihilate late in the early Universe, change the energy and entropy densities at BBN and at recombination. BBN, in combination with the CMB, can remove some of the degeneracies among light WIMPs and equivalent neutrinos, constraining the existence and properties of each. Depending on the nature of the light WIMP (Majorana or Dirac fermion, real or complex scalar) the joint BBN + CMB analyses set lower bounds to m_chi in the range 0.5 - 5 MeV (m_chi/m_e > 1 - 10), and they identify best fit values for m_chi in the range 5 - 10 MeV. The joint BBN + CMB analysis finds a best fit value for the number of equivalent neutrinos, Delta N_nu ~ 0.65, nearly independent of the nature of the WIMP. In the absence of a light WIMP (m_chi > 20 MeV), N_eff = 3.05(1 + Delta N_nu /3). In this case, there is excellent agreement between BBN and the CMB, but the joint fit reveals Delta N_nu = 0.40+-0.17, disfavoring standard big bang nucleosynthesis (SBBN) (Delta N_nu = 0) at ~ 2.4 sigma, as well as a sterile neutrino (Delta N_nu = 1) at ~ 3.5 sigma. The best BBN + CMB joint fit disfavors the absence of dark radiation (Delta N_nu = 0 at ~ 95% confidence), while allowing for the presence of a sterile neutrino (Delta N_nu = 1 at ~ 1 sigma). For all cases considered here, the lithium problem persists. These results, presented at the TAUP 2013 Conference, are based on Nollett & Steigman (2013).

Recent updates to the observational determinations of the primordial abundances of helium and deuterium are compared to the predictions of BBN to infer the universal ratio of baryons to photons (or, the present Universe baryon mass density parameter Omega_B h^2), as well as to constrain the effective number of neutrinos (N_eff) and the number of equivalent neutrinos (Delta N_nu). These BBN results are compared to those derived independently from the Planck CMB data. In the absence of a light WIMP (chi), N_eff = 3.05(1 + Delta N_nu/3). In this case, there is excellent agreement between BBN and the CMB, but the joint fit finds that Delta N_nu = 0.40 +/- 0.17, disfavoring standard big bang nucleosynthesis (SBBN: Delta N_nu = 0) at 2.4 sigma, as well as a sterile neutrino (Delta N_nu = 1) at 3.5 sigma. In the presence of a light WIMP, the relation between N_eff and Delta N_nu depends on the WIMP mass, leading to degeneracies among N_eff, Delta N_nu, and m_chi. The complementary and independent BBN and CMB data can break some of these degeneracies. Depending on the nature of the light WIMP (Majorana or Dirac fermion, real or complex scalar) the joint BBN + CMB analyses set a lower bound to m_chi in the range from 0.5 to 5 MeV, and they identify best fit values for m_chi in the range from 5 to 10 MeV. The joint BBN + CMB analyses find a best fit value for the number of equivalent neutrinos, Delta N_nu = 0.65, nearly independent of the nature of the WIMP. The best fit still disfavors the absence of dark radiation (Delta N_nu = 0 at 95% confidence), while allowing for the presence of a sterile neutrino (Delta N_nu = 1 at less than 1 sigma). For all cases considered here, the lithium problem persists. These results, presented at the 2013 Rencontres de l'Observatoire de Paris - ESO Workshop, are based on Nollett & Steigman 2013 (arXiv:1312.5725 [astro-ph.CO]).

We report a leading-order (LO) calculation of $^7\mathrm{Be}(p,\gamma)^8\mathrm{B}$ in a low-energy effective field theory. $^8\mathrm{B}$ is treated as a shallow proton$+^7\mathrm{Be}$ core and proton$+^7\mathrm{Be}^{*}$ (core excitation) $p$-wave bound state. The couplings are fixed using measured binding energies and proton-$^7\mathrm{Be}$ $s$-wave scattering lengths, together with $^8\mathrm{B}$ asymptotic normalization coefficients from ab initio calculations. We obtain a zero-energy $S$-factor of $18.2 \pm 1.2~({\rm ANC~only})$ eV b. Given that this is a LO result it is consistent with the recommended value $S(0)=20.8\pm1.6$ eV b. Our computed $S(E)$ compares favorably with experimental data on $^7\mathrm{Be}(p,\gamma)^8\mathrm{B}$ for $E <0.4$ MeV. We emphasize the important role of proton-$^7\mathrm{Be}$ scattering parameters in determining the energy dependence of $S(E)$, and demonstrate that their present uncertainties significantly limit attempts to extrapolate these data to stellar energies.

We report a leading-order calculation of radiative ${}^{7}\mathrm{Li}$ neutron captures to both the ground and first excited state of ${}^{8}\mathrm{Li}$ in the framework of a low-energy effective field theory (Halo-EFT). Each of the possible final states is treated as a shallow bound state composed of both $n+{}^{7}\mathrm{Li}$ and $n+{}^{7}\mathrm{Li}^{*}$ (core excitation) configurations. The {\it ab initio} variational Monte Carlo method is used to compute the asymptotic normalization coefficients of these bound states, which are then used to fix couplings in our EFT. We calculate the total and partial cross sections in the radiative capture process using this calibrated EFT. Fair agreement with measured total cross sections is achieved and excellent agreement with the measured branching ratio between the two final states is found. In contrast,a previous Halo-EFT calculation [G. Rupak and R. Higa, Phys. Rev. Lett {\bf 106}, 222501 (2011)] assumes that the $n$-${}^{7}\mathrm{Li}$ couplings in different spin channels are equal, fits the $P$-wave "effective-range" parameter to the threshold cross section for ${}^7{\rm Li} + n \rightarrow {}^8{\rm Li} + \gamma$, and assumes the core excitation is at high enough energy scale that it can be integrated out.

We present constraints on the number of relativistic species from a joint analysis of cosmic microwave background (CMB) fluctuations and light element abundances (helium and deuterium) compared to big bang nucleosynthesis (BBN) predictions. Our BBN calculations include updates of nuclear rates in light of recent experimental and theoretical information, with the most significant change occuring for the d(p,gamma)^3He cross section. We calculate a likelihood function for BBN theory and observations that accounts for both observational errors and nuclear rate uncertainties and can be easily embedded in cosmological parameter fitting. We then demonstrate that CMB and BBN are in good agreement, suggesting that the number of relativistic species did not change between the time of BBN and the time of recombination. The level of agreement between BBN and CMB, as well as the agreement with the standard model of particle physics, depends somewhat on systematic differences among determinations of the primordial helium abundance. We demonstrate that interesting constraints can be derived combining only CMB and D/H observations with BBN theory, suggesting that an improved D/H constraint would be an extremely valuable probe of cosmology.

We summarize and critically evaluate the available data on nuclear fusion cross sections important to energy generation in the Sun and other hydrogen-burning stars and to solar neutrino production. Recommended values and uncertainties are provided for key cross sections, and a recommended spectrum is given for 8B solar neutrinos. We also discuss opportunities for further increasing the precision of key rates, including new facilities, new experimental techniques, and improvements in theory. This review, which summarizes the conclusions of a workshop held at the Institute for Nuclear Theory, Seattle, in January 2009, is intended as a 10-year update and supplement to Reviews of Modern Physics 70 (1998) 1265.

The fraction of matter that is in the form of baryons or dark matter could have spatial fluctuations in the form of baryon-dark matter isocurvature fluctuations. We use big bang nucleosynthesis calculations compared with observed light element abundances as well as galaxy cluster gas fractions to constrain cosmological variations in the baryon fraction. Light element abundances constrain spatial variations to be less than 26-27%, while a sample of "relaxed" galaxy clusters shows spatial variations in gas fractions less than 8%. Larger spatial variations could cause differential screening of the primary cosmic microwave background anisotropies, leading to asymmetries in the fluctuations and ease some tension with the halo-star 7Li abundance. Fluctuations within our allowed bounds can lead to "B-mode" CMB polarization anisotropies at a non-negligible level.

It is known that there must be some weak form of transport (called cool bottom processing, or CBP) acting in low mass RGB and AGB stars, adding nuclei, newly produced near the hydrogen-burning shell, to the convective envelope. We assume that this extra-mixing originates in a stellar dynamo operated by the differential rotation below the envelope, maintaining toroidal magnetic fields near the hydrogen-burning shell. We use a phenomenological approach to the buoyancy of magnetic flux tubes, assuming that they induce matter circulation as needed by CBP models. This establishes requirements on the fields necessary to transport material from zones where some nuclear burning takes place, through the radiative layer, and into the convective envelope. Magnetic field strengths are determined by the transport rates needed by CBP for the model stellar structure of a star of initially 1.5 solar mass, in both the AGB and RGB phases. The field required for the AGB star in the processing zone is B_0 ~ 5x10^6 G; at the base of the convective envelope this yields an intensity B_E < 10^4 G (approximately). For the RGB case, B_0 ~ 5x10^4 to 4x10^5 G, and the corresponding B_E are ~ 450 to 3500 G. These results are consistent with existing observations on AGB stars. They also hint at the basis for high field sources in some planetary nebulae and the very large fields found in some white dwarfs. It is concluded that transport by magnetic buoyancy should be considered as a possible mechanism for extra mixing through the radiative zone, as is required by both stellar observations and the extensive isotopic data on circumstellar condensates found in meteorites.

(Abridged) We review abundances of short-lived nuclides in the early solar system (ESS) and the methods used to determine them. We compare them to the inventory for a uniform galactic production model. Within a factor of two, observed abundances of several isotopes are compatible with this model. I-129 is an exception, with an ESS inventory much lower than expected. The isotopes Pd-107, Fe-60, Ca-41, Cl-36, Al-26, and Be-10 require late addition to the solar nebula. Be-10 is the product of particle irradiation of the solar system as probably is Cl-36. Late injection by a supernova (SN) cannot be responsible for most short-lived nuclei without excessively producing Mn-53; it can be the source of Mn-53 and maybe Fe-60. If a late SN is responsible for these two nuclei, it still cannot make Pd-107 and other isotopes. We emphasize an AGB star as a source of nuclei, including Fe-60 and explore this possibility with new stellar models. A dilution factor of about 4e-3 gives reasonable amounts of many nuclei. We discuss the role of irradiation for Al-26, Cl-36 and Ca-41. Conflict between scenarios is emphasized as well as the absence of a global interpretation for the existing data. Abundances of actinides indicate a quiescent interval of about 1e8 years for actinide group production in order to explain the data on Pu-244 and new bounds on Cm-247. This interval is not compatible with Hf-182 data, so a separate type of r-process is needed for at least the actinides, distinct from the two types previously identified. The apparent coincidence of the I-129 and trans-actinide time scales suggests that the last actinide contribution was from an r-process that produced actinides without fission recycling so that the yields at Ba and below were governed by fission.

(Abridged) We examine the effects of cool bottom processing (CBP) on several isotopic ratios in the convective envelope during the TP-AGB phase of evolution in a 1.5 M_sun initial-mass star of solar initial composition. We use a parametric model which treats extra mixing by introducing mass flow between the convective envelope and the underlying radiative zone. The parameters of this model are the mass circulation rate (Mdot) and the maximum temperature (T_P) experienced by the circulating material. The effects of nuclear reactions in the flowing matter were calculated using a set of structures of the radiative zone selected from a complete stellar evolution calculation. The compositions of the flowing material were obtained and the resulting changes in the envelope determined. Abundant ^26Al was produced by CBP for log T_P > 7.65. While ^26Al/^27Al depends on T_P, the isotopic ratios in CNO elements depend dominantly on the circulation rate. The correspondence is shown between models of CBP as parameterized by a diffusion formalism within the stellar evolution model and those using the mass-flow formalism employed here. The isotopic ratios are compared with the data on circumstellar dust grains. It is found that the ratios ^{18}O/^{16}O, ^{17}O/^{16}O, and ^26Al/^27Al observed for oxide grains formed at C/O < 1 are reasonably well-understood. However, the ^15N/^14N, ^12C/^13C, and ^26Al/^27Al in carbide grains (C/O > 1) require many stellar sources with ^14N/^15N at least a factor of 4 below solar. The rare grains with ^12C/^13C < 10 cannot be produced by any red-giant or AGB source.

We compute primordial light-element abundances for cases with fine structure constant alpha different from the present value, including many sources of alpha dependence neglected in previous calculations. Specifically, we consider contributions arising from Coulomb barrier penetration, photon coupling to nuclear currents, and the electromagnetic components of nuclear masses. We find the primordial abundances to depend more weakly on alpha than previously estimated, by up to a factor of 2 in the case of ^7Li. We discuss the constraints on variations in alpha from the individual abundance measurements and the uncertainties affecting these constraints. While the present best measurements of primordial D/H, ^4He/H, and ^7Li/H may be reconciled pairwise by adjusting alpha and the universal baryon density, no value of alpha allows all three to be accommodated simultaneously without consideration of systematic error. The combination of measured abundances with observations of acoustic peaks in the cosmic microwave background favors no change in alpha within the uncertainties.

Recent studies of heavy r-process elements in low [Fe/H] halo stars have suggested that an initial population of metal-free very massive stars (VMSs) may be required to provide early Fe enrichment without coproducing heavy r nuclei. We find similar abundance trends in alpha-elements (which should be copiously produced by VMSs), but not in other elements such as carbon (which should not), in agreement with this hypothesis. We then combine the corresponding level of prompt initial enrichment with models of VMS nucleosynthetic yields and spectra to estimate the corresponding ionizing fluxes. The result suggests that there may have been enough VMS activity to reionize the universe. The unusually hard spectrum of VMSs would imply a different reionization history from canonical models. HeII could have been reionized at high redshift, only to recombine as a subsequent generation of stars formed with a ``normal'' initial mass function.

The determination of the primeval deuterium abundance has opened a precision era in big-bang nucleosynthesis (BBN), making accurate predictions more important than ever before. We present in analytic form new, more precise predictions for the light-element abundances and their error matrix. Using our predictions and the primeval deuterium abundance we infer a baryon density of Omega_B h^2 = 0.020 pm 0.002 (95% cl) and find no evidence for stellar production (or destruction) of 3He beyond burning D to 3He. Conclusions about 4He and 7Li currently hinge upon possible systematic error in their measurements.

We report the first calculations of cross sections for the radiative capture reactions t(alpha,gamma)^7Li and ^3He(alpha,gamma)^7Be below 2 MeV which use wave functions derived from realistic nucleon-nucleon interactions by the variational Monte Carlo technique. After examining several small corrections to the dominant E1 operator, we find energy dependences for the low-energy S-factors which agree reasonably with experimental measurements. There is no contradiction with the previous theoretical understanding of these processes, but the zero-energy derivative of the t(alpha,gamma)^7Li S-factor is smaller than in most models. While this method can in principle predict cross section normalizations, the normalizations of our results are mostly too low.

We have computed the cross section for the process d + alpha -> 6Li + gamma at the low energies relevant for primordial nucleosynthesis and comparison with laboratory data. The final state is a six-body wave function generated by the variational Monte Carlo method from the Argonne v_18 and Urbana IX potentials, including improved treatment of large-particle-separation behavior. The initial state is built up from the alpha-particle and deuteron ground-state solutions for these potentials, with phenomenological descriptions of scattering and cluster distortions. The dominant E2 cross section is in reasonable agreement with the laboratory data. Including center-of-energy and other small corrections, we obtain an E1 contribution which is larger than the measured contribution at 2 MeV by a factor of 7. We calculate explicitly the impulse-approximation M1 contribution, which is expected to be very small, and obtain a result consistent with zero. We find little reason to suspect that the cross section is large enough to produce significant 6Li in the big bang.

Together, the standard theory of big-bang nucleosynthesis (BBN) and the primeval deuterium abundance now very precisely peg the baryon density. Based upon our analysis of the deuterium data and the theoretical uncertainties associated with the BBN predictions, we determine Omega_B h^2 = 0.020\pm 0.002 (95% C.L.), with the uncertainty from the measured deuterium abundance about twice that from the predicted abundance. We discuss critically the reliability of the BBN baryon density, and in light of possible systematic uncertainties also derive a very conservative range. We conclude that within the standard cosmology and standard theory of BBN a baryon density Omega_B h^2=0.032 (the central value implied by recent CMB anisotropy measurements) simply cannot be accommodated.

We present a Monte Carlo method for direct incorporation of nuclear inputs in primordial nucleosynthesis calculations. This method is intended to remedy shortcomings of current error estimation, by eliminating intermediate data evaluations and working directly with experimental data, allowing error estimation based solely on published experimental uncertainties. This technique also allows simple incorporation of new data and reduction of errors with the introduction of more precise data. Application of our method indicates that previous error estimates on the calculated abundances were too large by as much as a factor of three. Since uncertainties in the BBN calculation currently dominate inferences drawn from light-element abundances, the re-estimated errors have important consequences for cosmic baryon density, neutrino physics, and lithium depletion in halo stars. Our direct method allows detailed discussion of the status of the nuclear inputs, by identifying clearly the places where improved cross section measurements would be most useful.

We summarize the history, theory, observational status, and implications of big-bang nucleosynthesis.

Motivated by the recent measurement of the primeval abundance of deuterium, we re-examine the nuclear inputs to big-bang nucleosynthesis (BBN). Using Monte-Carlo realization of the nuclear cross-section data to directly estimate the theoretical uncertainties for the yields of D, 3-He and 7-Li, we show that previous estimates were a factor of 2 too large. We sharpen the BBN determination of the baryon density based upon deuterium, rho_B = (3.6 +/- 0.4) * 10^{-31} g/cm^3 (Omega_B h^2 = 0.019 +/- 0.0024), which leads to a predicted 4-He abundance, Y_P = 0.246 +/- 0.0014 and a stringent limit to the equivalent number of light neutrino species: N_nu < 3.20 (all at 95% cl). The predicted 7-Li abundance, 7-Li/H = (3.5 + 1.1 - 0.9) * 10^{-10}, is higher than that observed in pop II stars, (1.7 +/- 0.3) * 10^{-10} (both, 95% cl). We identify key reactions and the energies where further work is needed.

We examine the possibility that Big Bang Nucleosynthesis (BBN) may produce non-trivial amounts of $^6$Li. If a primordial component of this isotope could be observed, it would provide a new fundamental test of Big-Bang cosmology, as well as new constraints on the baryon density of the universe. At present, however, theoretical predictions of the primordial $^6$Li abundance are extremely uncertain due to difficulties in both theoretical estimates and experimental determinations of the deuterium-plus-alpha radiative capture reaction cross-section. We also argue that present observational capabilities do not yet allow the detection of primeval $^6$Li in very metal-poor stars of the galactic halo. However, if the critical cross section is towards the upper end of its plausible range, then improvements in $^6$Li detection capabilities may allow the establishment of $^6$Li as another product of BBN. It is also noted that a primordial $^6$Li detection could help resolve current concerns about the extragalactic D/H determination.