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'.

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. We explore novel forms of the trial wave function and find one which leads to a ground state of the eight-particle system whose energy is almost equal to that of two four-particle systems. We investigate the clustering properties involved and also extrapolate to the zero-range limit. In addition to being experimentally testable, our results impact the prospects of developing nuclear physics as a perturbation around the unitary limit.

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. We use the single-particle spectrum to extract the effective mass in neutron matter. With a view to systematizing the error involved in this extraction, we carefully assess the impact of finite-size effects on the quasiparticle dispersion relation. We find an effective-mass ratio that drops from 1 as the density is increased. We conclude by connecting our results with the physics of ultracold gases as well as with energy-density functional theories of nuclei and neutron-star matter.

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. In order to remove systematic errors, we examine finite-size effects and the impact of the wave function ansatz. We also make contact with energy-density functional theories of nuclei and disentangle isovector gradient contributions from bulk properties. Finally, we calculate the static density-density linear response function of neutron matter and compare it with the response of other physical systems.

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). We also apply the new functional to few-particle harmonically trapped systems, comparing with previous calculations.

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. We also solve the Bogoliubov-de Gennes equations with variable constraints on the spin-singlet and spin-triplet pairing amplitudes. For nuclei that exhibit this new pairing behavior, the resulting energy surface can be rather soft, suggesting that there may be low-lying excitations associated with the spin mixing.

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. We calculate the ground-state energies for particles in a harmonic trap by using the nonlocal induced 3-body interaction and compare them to numerically accurate Green's function Monte Carlo calculations. Using no free parameters, we find that a nonlocality in the space domain provides a better description of the weak-coupling regime than the local-density approximation.

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. For a majority of heavy particles, we find an energy minimum for a normal state with a ratio of approximately 3ratio1 heavy to light particles. A slight increase in attraction to k_{F}a approximately 2.5 yields a ground state energy of nearly zero for this ratio. A cold unpolarized system in a harmonic trap at unitarity should phase separate into three regions, with a shell of unpolarized superfluid in the middle.