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.

Clustering of Four-Component Unitary Fermions.

Ab initio nuclear physics tackles the problem of strongly interacting four-component fermions. The same setting could foreseeably be probed experimentally in ultracold atomic systems, where two- and three-component experiments have led to major breakthroughs in recent years. Both due to the problem's inherent interest and as a pathway to nuclear physics, in this Letter we study four-component fermions at unitarity via the use of quantum Monte Carlo methods.

Nonperturbative Extraction of the Effective Mass in Neutron Matter.

We carry out nonperturbative calculations of the single-particle excitation spectrum in strongly interacting neutron matter. These are microscopic quantum Monte Carlo computations of many-neutron energies at different densities as well as several distinct excited states. As input, we employ both phenomenological and chiral two- and three-nucleon interactions.

Static Response of Neutron Matter.

We generalize the problem of strongly interacting neutron matter by adding a periodic external modulation. This allows us to study from first principles a neutron system that is extended and inhomogeneous, with connections to the physics of both neutron-star crusts and neutron-rich nuclei. We carry out fully nonperturbative microscopic quantum Monte Carlo calculations of the energy of neutron matter at different densities, as well as different strengths and periodicities of the external potential.

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.

Quantum Monte Carlo calculations of light nuclei using chiral potentials.

We present the first Green's function Monte Carlo calculations of light nuclei with nuclear interactions derived from chiral effective field theory up to next-to-next-to-leading order. Up to this order, the interactions can be constructed in a local form and are therefore amenable to quantum Monte Carlo calculations. We demonstrate a systematic improvement with each order for the binding energies of A=3 and A=4 systems.

Quantum Monte Carlo calculations with chiral effective field theory 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.

Mixed-spin pairing condensates in heavy nuclei.

We show that the Bogoliubov-de Gennes equations for nuclear ground-state wave functions support solutions in which the condensate has a mixture of spin-singlet and spin-triplet pairing. We find that such mixed-spin condensates do not occur when there are equal numbers of neutrons and protons, but only when there is an isospin imbalance. Using a phenomenological Hamiltonian, we predict that such nuclei may occur in the physical region within the proton dripline.

Resonantly interacting fermions in a box.

We use two fundamental theoretical frameworks to study the finite-size (shell) properties of the unitary gas in a periodic box: (1) an ab initio quantum Monte Carlo (QMC) calculation for boxes containing 4 to 130 particles provides a precise and complete characterization of the finite-size behavior, and (2) a new density functional theory (DFT) fully encapsulates these effects. The DFT predicts vanishing shell structure for systems comprising more than 50 particles, and allows us to extrapolate the QMC results to the thermodynamic limit, providing the tightest bound to date on the ground-state energy of the unitary gas: ξ(S)≤0.383(1).

Effective 3-body interaction for mean-field and density-functional theory.

Density functionals for nuclei usually include an effective 3-body interaction that depends on a fractional power of the density. Using insights from the many-body theory of the low-density two-component Fermi gas, we consider a new, nonlocal, form for the energy functional that is consistent with the Fock-space representation of interaction operators. In particular, there is a unique spatially nonlocal generalization of the contact form of the interaction that preserves the ρ(7/3) density dependence required by the many-body theory.

Heavy-light fermion mixtures at unitarity.

We investigate fermion pairing in the unitary regime for a mass ratio corresponding to a ;{6}Li-;{40}K mixture using quantum Monte Carlo methods. The ground-state energy and the average light- and heavy-particle excitation spectrum for the unpolarized superfluid state are nearly independent of the mass ratio. In the majority light system, the polarized superfluid is close to the energy of a phase separated mixture of nearly fully polarized normal and unpolarized superfluid.