We report benchmark calculations of the energy per particle of pure neutron matter as a function of the baryon density using three independent many-body methods: Brueckner-Bethe-Goldstone, Fermi hypernetted chain/single-operator chain, and auxiliary-field diffusion Monte Carlo. Significant technical improvements are implemented in the latter two methods. The calculations are made for two distinct families of realistic coordinate-space nucleon-nucleon potentials fit to scattering data, including the standard Argonne $v_{18}$ interaction and two of its simplified versions, and four of the new Norfolk $\Delta$-full chiral effective field theory potentials. The results up to twice nuclear matter saturation density show some divergence among the methods, but improved agreement compared to earlier work. We find that the potentials fit to higher-energy nucleon-nucleon scattering data exhibit a much smaller spread of energies.

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 review the calculation of the equation of state of pure neutron matter using quantum Monte Carlo (QMC) methods. QMC algorithms permit the study of many-body nuclear systems using realistic two- and three-body forces in a nonperturbative framework. We present the results for the equation of state of neutron matter, and focus on the role of three-neutron forces at supranuclear density. We discuss the correlation between the symmetry energy, the neutron star radius and the symmetry energy. We also combine QMC and theoretical models of the three-nucleon interactions, and recent neutron star observations to constrain the value of the symmetry energy and its density dependence.

We examine how the binding of light ($A\leq 8$) nuclei depends on possible variations of hadronic masses, including meson, nucleon, and nucleon-resonance masses. Small variations in hadronic masses may have occurred over time; the present results can help evaluate the consequences for big bang nucleosynthesis. Larger variations may be relevant to current attempts to extrapolate properties of nucleon-nucleon interactions from lattice QCD calculations. Results are presented as derivatives of the energy with respect to the different masses so they can be combined with different predictions of the hadronic mass-dependence on the underlying current-quark mass $m_q$. As an example, we employ a particular set of relations obtained from a study of hadron masses and sigma terms based on Dyson-Schwinger equations and a Poincar\'{e}-covariant Faddeev equation for confined quarks and diquarks. We find that nuclear binding decreases moderately rapidly as the quark mass increases, with the deuteron becoming unbound when the pion mass is increased by $\sim$60% (corresponding to an increase in $X_q=m_q/\Lambda_{QCD}$ of 2.5). In the other direction, the dineutron becomes bound if the pion mass is decreased by $\sim$15% (corresponding to a reduction of $X_q$ by $\sim$30%). If we interpret the disagreement between big bang nucleosynthesis calculations and measurements to be the result of variation in $X_q$, we obtain an estimate $\delta X_q/X_q=K \cdot (0.013 \pm 0.002)$ where $K \sim 1$ (the expected accuracy in $K$ is about a factor of 2). The result is dominated by $^7$Li data.

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.