Quantum Monte Carlo simulations of the Fermi-polaron problem and bosons with Gaussian interactions

C. Lavalle, Marcos Rigol, M. Feldbacher +3 more

Lode Pollet

This is a review of recent developments in Monte Carlo methods in the field of ultracold gases. For bosonic atoms in an optical lattice we discuss path-integral Monte Carlo simulations with worm updates and show the excellent agreement with cold atom experiments. We also review recent progress in simulating bosonic systems with long-range interactions, disordered bosons, mixtures of bosons and spinful bosonic systems. For repulsive fermionic systems, determinantal methods at half filling are sign free, but in general no sign-free method exists. We review the developments in diagrammatic Monte Carlo for the Fermi polaron problem and the Hubbard model, and show the connection with dynamical mean-field theory. We end the review with diffusion Monte Carlo for the Stoner problem in cold gases.

Michael Kolodrubetz, Bryan K. Clark

Finding the ground state of a fermionic Hamiltonian using quantum Monte Carlo (QMC) is a very difficult problem, due to the Fermi sign problem. While still scaling exponentially, full configuration-interaction Monte Carlo (FCIQMC) mitigates some of the exponential variance by allowing annihilation of noise---whenever two walkers arrive at the same configuration with opposite signs, they are removed from the simulation. While FCIQMC has been quite successful for quantum chemistry problems, its application to problems in condensed systems has been limited. In this paper, we apply the FCIQMC algorithm to the Fermi polaron problem, which provides an ideal test bed for improving the algorithm, although since we restrict the number of many-body excitations, our algorithm is more precisely a QMC implementation of the configuration-interaction method (CIQMC). In its simplest form, CIQMC is unstable for even a fairly small system sizes. However, with a series of algorithmic improvements, we are able to significantly increase its effectiveness. We modify fixed node QMC to work in these systems, introduce a well chosen importance sampled trial wave function, a partial node approximation, and a variant of release node. Finally, we develop a way to perform CIQMC directly in the thermodynamic limit.

P. E. Kornilovitch

I. Morgenstern

Joel Frederick Corney, P. D. Drummond

We introduce a new class of quantum Monte Carlo methods, based on a Gaussian quantum operator representation of fermionic states. The methods enable first-principles dynamical or equilibrium calculations in many-body Fermi systems, and, combined with the existing Gaussian representation for bosons, provide a unified method of simulating Bose-Fermi systems. As an application relevant to the Fermi sign problem, we calculate finite-temperature properties of the two dimensional Hubbard model and the dynamics in a simple model of coherent molecular dissociation.