Quantum Monte Carlo methods have proved 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. The nuclear interactions and currents are reviewed along with a description of 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. A variety of results are presented, including the low-lying spectra of light nuclei, nuclear form factors, and transition matrix elements. 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 are also described. A coherent picture of nuclear structure and dynamics emerges based upon rather simple but realistic interactions and currents.

The symmetry energy contribution to the nuclear equation of state impacts various phenomena in nuclear astrophysics, nuclear structure, and nuclear reactions. Its determination is a key objective of contemporary nuclear physics, with consequences for the understanding of dense matter within neutron stars. We examine the results of laboratory experiments that have provided initial constraints on the nuclear symmetry energy and on its density dependence at and somewhat below normal nuclear matter density. Even though some of these constraints have been derived from properties of nuclei while others have been derived from the nuclear response to electroweak and hadronic probes, within experimental uncertainties-they are consistent with each other. We also examine the most frequently used theoretical models that predict the symmetry energy and its slope parameter. By comparing existing constraints on the symmetry pressure to theories, we demonstrate how contributions of three-body forces, which are essential ingredients in neutron matter models, can be determined.

The dense matter equation of state (EOS) determines neutron star (NS) structure but can be calculated reliably only up to one to two times the nuclear saturation density, using accurate many-body methods that employ nuclear interactions from chiral effective field theory constrained by scattering data. In this work, we use physically motivated ansatzes for the speed of sound c _S at high density to extend microscopic calculations of neutron-rich matter to the highest densities encountered in stable NS cores. We show how existing and expected astrophysical constraints on NS masses and radii from X-ray observations can constrain the speed of sound in the NS core. We confirm earlier expectations that c _S is likely to violate the conformal limit of , possibly reaching values closer to the speed of light c at a few times the nuclear saturation density, independent of the nuclear Hamiltonian. If QCD obeys the conformal limit, we conclude that the rapid increase of c _S required to accommodate a 2 M _⊙ NS suggests a form of strongly interacting matter where a description in terms of nucleons will be unwieldy, even between one and two times the nuclear saturation density. For typical NSs with masses in the range of 1.2–1.4 M _⊙, we find radii between 10 and 14 km, and the smallest possible radius of a 1.4 M _⊙ NS consistent with constraints from nuclear physics and observations is 8.4 km. We also discuss how future observations could constrain the EOS and guide theoretical developments in nuclear physics.

We calculate the equation of state of neutron matter with realistic two- and three-nucleon interactions using quantum Monte Carlo techniques and illustrate that the short-range three-neutron interaction determines the correlation between neutron matter energy at nuclear saturation density and higher densities relevant to neutron stars. Our model also makes an experimentally testable prediction for the correlation between the nuclear symmetry energy and its density dependence—determined solely by the strength of the short-range terms in the three-neutron force. The same force provides a significant constraint on the maximum mass and radius of neutron stars.

We present quantum Monte Carlo calculations of light nuclei, neutron-$\alpha$ scattering, and neutron matter using local two- and three-nucleon ($3N$) interactions derived from chiral effective field theory up to next-to-next-to-leading order (${\mathrm{N}}^2\mathrm{LO}$). The two undetermined $3N$ low-energy couplings are fit to the ${}^4\mathrm{He}$ binding energy and, for the first time, to the spin-orbit splitting in the neutron-$\alpha$ $P$-wave phase shifts. Furthermore, we investigate different choices of local $3N$-operator structures and find that chiral interactions at ${\mathrm{N}}^2\mathrm{LO}$ are able to simultaneously reproduce the properties of $A=3,4,5$ systems and of neutron matter, in contrast to commonly used phenomenological $3N$ interactions.

We present the first quantum Monte Carlo (QMC) calculations with chiral effective field theory (EFT) interactions. To achieve this, we remove all sources of nonlocality, which hamper the inclusion in QMC calculations, in nuclear forces to next-to-next-to-leading order. We perform auxiliary-field diffusion Monte Carlo (AFDMC) calculations for the neutron matter energy up to saturation density based on local leading-order, next-to-leading order, and next-to-next-to-leading order nucleon-nucleon interactions. Our results exhibit a systematic order-by-order convergence in chiral EFT and provide nonperturbative benchmarks with theoretical uncertainties. For the softer interactions, perturbative calculations are in excellent agreement with the AFDMC results. This work paves the way for QMC calculations with systematic chiral EFT interactions for nuclei and nuclear matter, for testing the perturbativeness of different orders, and allows for matching to lattice QCD results by varying the pion mass.