Quantum Monte Carlo methods have proved very valuable to study the structure and reactions of light nuclei and nucleonic matter starting from realistic nuclear interactions and currents. These ab-initio calculations reproduce many low-lying states, moments and transitions in light nuclei, and simultaneously predict many properties of light nuclei and neutron matter over a rather wide range of energy and momenta. We review the nuclear interactions and currents, and describe the continuum Quantum Monte Carlo methods used in nuclear physics. These methods are similar to those used in condensed matter and electronic structure but naturally include spin-isospin, tensor, spin-orbit, and three-body interactions. We present a variety of results including the low-lying spectra of light nuclei, nuclear form factors, and transition matrix elements. We also describe low-energy scattering techniques, studies of the electroweak response of nuclei relevant in electron and neutrino scattering, and the properties of dense nucleonic matter as found in neutron stars. A coherent picture of nuclear structure and dynamics emerges based upon rather simple but realistic interactions and currents.

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.

Spurred by the recent complete determination of the weak currents in two-nucleon systems up to ${\cal O}(Q^3)$ in heavy-baryon chiral perturbation theory, we carry out a parameter-free calculation of the threshold $S$-factors for the solar $pp$ (proton-fusion) and $hep$ processes in an effective field theory that {\it combines} the merits of the standard nuclear physics method and systematic chiral expansion. The power of the EFT adopted here is that one can correlate in a unified formalism the weak-current matrix elements of two-, three- and four-nucleon systems. Using the tritium $\beta$-decay rate as an input to fix the only unknown parameter in the theory, we can evaluate the threshold $S$ factors with drastically improved precision; the results are $S_{pp}(0) = 3.94\times(1 \pm 0.004) \times 10^{-25} {MeV-b}$ and $S_{hep}(0) = (8.6\pm 1.3)\times 10^{-20} {keV-b}$. The dependence of the calculated $S$-factors on the momentum cutoff parameter $\Lambda$ has been examined for a physically reasonable range of $\Lambda$. This dependence is found to be extremely small for the $pp$ process, and to be within acceptable levels for the $hep$ process, substantiating the consistency of our calculational scheme.

Using effective field theory, we calculate the S-factor for the hep process in a totally parameter-free formulation. The transition operators are organized according to chiral counting, and their matrix elements are evaluated using the realistic nuclear wave functions obtained in the correlated-hyperspherical-harmonics method. Terms of up to next-to-next-to-next-to-leading order in heavy-baryon chiral perturbation theory are considered. Fixing the only parameter in the theory by fitting the tritium \beta-decay rate, we predict the hep S-factor with accuracy better than \sim 20 %.

Spurred by the recent complete determination of the weak currents in two-nucleon systems up to ${\cal O}(Q^3)$ in heavy-baryon chiral perturbation theory, we carry out a parameter-free calculation of the solar proton fusion rate in an effective field theory that combines the merits of the standard nuclear physics method and systematic chiral expansion. Using the tritium beta-decay rate as an input to fix the only unknown parameter in the effective Lagrangian, we can evaluate with drastically improved precision the ratio of the two-body contribution to the well established one-body contribution; the ratio is determined to be (0.86\pm 0.05) %. This result is essentially independent of the cutoff parameter for a wide range of its variation (500 MeV \le \Lambda \le 800 MeV), a feature that substantiates the consistency of the calculation.

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.

The astrophysical S-factor for the proton weak capture on 3He is calculated with correlated-hyperspherical-harmonics bound and continuum wave functions corresponding to realistic Hamiltonians consisting of the Argonne v14 or Argonne v18 two-nucleon and Urbana-VIII or Urbana-IX three-nucleon interactions. The nuclear weak charge and current operators have vector and axial-vector components, that include one- and many-body terms. All possible multipole transitions connecting any of the p 3He S- and P-wave channels to the 4He bound state are considered. The S-factor at a p 3He center-of-mass energy of 10 keV, close to the Gamow-peak energy, is predicted to be 10.1 10^{-20} keV b with the AV18/UIX Hamiltonian, a factor of about 4.5 larger than the value adopted in the standard solar model. The P-wave transitions are found to be important, contributing about 40 % of the calculated S-factor. The energy dependence is rather weak: the AV18/UIX zero-energy S-factor is 9.64 10^{-20} keV b, only 5 % smaller than the 10 keV result quoted above. The model dependence is also found to be weak: the zero-energy S-factor is calculated to be 10.2 10^{-20} keV b with the older AV14/UVIII model, only 6 % larger than the AV18/UIX result. Our best estimate for the S-factor at 10 keV is therefore (10.1 \pm 0.6) 10^{-20} keV b, when the theoretical uncertainty due to the model dependence is included. This value for the calculated S-factor is not as large as determined in fits to the Super-Kamiokande data in which the hep flux normalization is free. However, the precise calculation of the S-factor and the consequent absolute prediction for the hep neutrino flux will allow much greater discrimination among proposed solar neutrino oscillation solutions.

The astrophysical factor for the proton weak capture on 3He is calculated with correlated-hyperspherical-harmonics bound and continuum wave functions corresponding to a realistic Hamiltonian consisting of the Argonne v18 two-nucleon and Urbana-IX three-nucleon interactions. The nuclear weak charge and current operators have vector and axial-vector components, that include one- and many-body terms. All possible multipole transitions connecting any of the p-3He S- and P-wave channels to the 4He bound state are considered. The S-factor at a p-3He center-of-mass energy of 10 keV, close to the Gamow-peak energy, is predicted to be 10.1 10^{-20} keV b, a factor of five larger than the standard-solar-model value. The P-wave transitions are found to be important, contributing about 40 % of the calculated S-factor.