Quantum Monte Carlo calculations of light nuclei using chiral potentials

Auxiliary field Quantum Monte Carlo for dilute neutrons on the lattice

We employ constrained path Auxiliary Field Quantum Monte Carlo (AFQMC) in the pursuit of studying physical nuclear systems using a lattice formalism. Since AFQMC has been widely used in the study of condensed-matter systems such as the Hubbard model, we benchmark our method against published results for both one- and two-dimensional Hubbard model calculations. We then turn our attention to cold atomic and nuclear systems. We use an onsite contact interaction that can be tuned in order to reproduce the known scattering length and effective range of a given interaction. Developing this machinery allows us to extend our calculations to study nuclear systems within a lattice formalism. We perform initial calculations for a range of nuclear systems from two- to few-body neutron systems. This article is part of the theme issue 'The liminal position of Nuclear Physics: from hadrons to neutron stars'.

Trends of Neutron Skins and Radii of Mirror Nuclei from First Principles

The neutron skin of atomic nuclei impacts the structure of neutron-rich nuclei, the equation of state of nucleonic matter, and the size of neutron stars. Here we predict the neutron skin of selected light- and medium-mass nuclei using coupled-cluster theory and the auxiliary field diffusion Monte Carlo method with two- and three-nucleon forces from chiral effective field theory. We find a linear correlation between the neutron skin and the isospin asymmetry in agreement with the liquid-drop model and compare with data. We also extract the linear relationship that describes the difference between neutron and proton radii of mirror nuclei and quantify the effect of charge symmetry breaking terms in the nuclear Hamiltonian. Our results for the mirror-difference charge radii and binding energies per nucleon agree with existing data.

Short-Range Correlations and the EMC Effect in Effective Field Theory

We show that the empirical linear relation between the magnitude of the EMC effect in deep inelastic scattering on nuclei and the short-range correlation scaling factor a_{2} extracted from high-energy quasielastic scattering at x≥1 is a natural consequence of scale separation and derive the relationship using effective field theory. While the scaling factor a_{2} is a ratio of nuclear matrix elements that individually depend on the calculational scheme, we show that the ratio is independent of this choice. We perform Green's function Monte Carlo calculations with both chiral and Argonne-Urbana potentials to verify this and determine the scaling factors for light nuclei. The resulting values for ^{3}He and ^{4}He are in good agreement with experimental values. We also present results for ^{9}Be and ^{12}C extracted from variational Monte Carlo calculations.

Chiral Three-Nucleon Interactions in Light Nuclei, Neutron-α Scattering, and Neutron Matter

We present quantum Monte Carlo calculations of light nuclei, neutron-α 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 (N(2)LO). The two undetermined 3N low-energy couplings are fit to the (4)He binding energy and, for the first time, to the spin-orbit splitting in the neutron-α P-wave phase shifts. Furthermore, we investigate different choices of local 3N-operator structures and find that chiral interactions at N(2)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.