In recent years, the combination of precise quantum Monte Carlo (QMC) methods with realistic nuclear interactions and consistent electroweak currents, in particular those constructed within effective ...field theories (EFTs), has led to new insights in light and medium-mass nuclei, neutron matter, and electroweak reactions. For example, with the same chiral interactions, QMC calculations can reproduce binding energies and radii for light nuclei,
n
-α scattering phase shifts, and the neutron matter equation of state. This compelling new body of work has been made possible both by advances in QMC methods for nuclear physics, which push the bounds of applicability to heavier nuclei and to asymmetric nuclear matter, and by the development of local chiral EFT interactions up to next-to-next-to-leading order and minimally nonlocal interactions including Δ degrees of freedom. In this review, we discuss these recent developments and give an overview of the exciting results for nuclei, neutron matter and neutron stars, and electroweak reactions.
We carry out an ab initio calculation of the neutrino flux-folded inclusive cross sections measured on ^{12}C by the MiniBooNE and T2K Collaborations in the charged-current quasielastic regime. The ...calculation is based on realistic two- and three-nucleon interactions, and on a realistic nuclear electroweak current with one- and two-nucleon terms that are constructed consistently with these interactions and reproduce low-energy electroweak transitions. Numerically exact quantum Monte Carlo methods are utilized to compute the nuclear weak response functions, by fully retaining many-body correlations in the initial and final states and interference effects between one- and two-body current contributions. We employ a nucleon axial form factor of the dipole form with Λ_{A}=1.0 or 1.15 GeV, the latter more in line with a very recent lattice QCD determination. The calculated cross sections are found to be in good agreement with the neutrino data of MiniBooNE and T2K, and antineutrino MiniBooNE data, yielding a consistent picture of nuclei and their electroweak properties across a wide regime of energy and momenta.
In recent years local chiral interactions have been derived and implemented in quantum Monte Carlo methods in order to test to what extent the chiral effective field theory framework impacts our ...knowledge of few- and many-body systems. In this Letter, we present Green's function Monte Carlo calculations of light nuclei based on the family of local two-body interactions presented by our group in a previous paper in conjunction with chiral three-body interactions fitted to bound- and scattering-state observables in the three-nucleon sector. These interactions include Δ intermediate states in their two-pion-exchange components. We obtain predictions for the energy levels and level ordering of nuclei in the mass range A=4-12, accurate to ≤2% of the binding energy, in very satisfactory agreement with experimental data.
Nucleon momentum distributions calculated with a common one-body operator vary with the resolution scale (and scheme) of the Hamiltonian used. For high-resolution potentials such as Argonne v18 ...(AV18) there is a high-momentum tail, reflecting short-range correlations in the nuclear wave function, which is reduced or absent for softer, lower-resolution interactions. We explore if the similarity renormalization group (SRG) can be used to quantitatively reproduce the high-resolution distributions from variational Monte Carlo at all momenta using SRG-evolved operators and empirically fit single-particle orbitals rather than a full RG evolution of many-body wave functions. The goal of this approach is to enable calculations of high-resolution distributions for a wider range of nuclei as well as for other interactions, and provides connections to phenomenological analyses of experiments.
We extend the prediction range of Pionless Effective Field Theory with an analysis of the ground state of 16O in leading order. To renormalize the theory, we use as input both experimental data and ...lattice QCD predictions of nuclear observables, which probe the sensitivity of nuclei to increased quark masses. The nuclear many-body Schrödinger equation is solved with the Auxiliary Field Diffusion Monte Carlo method. For the first time in a nuclear quantum Monte Carlo calculation, a linear optimization procedure, which allows us to devise an accurate trial wave function with a large number of variational parameters, is adopted. The method yields a binding energy of 4He which is in good agreement with experiment at physical pion mass and with lattice calculations at larger pion masses. At leading order we do not find any evidence of a 16O state which is stable against breakup into four 4He, although higher-order terms could bind 16O.