Gaussian quantum Monte Carlo methods for fermions

Abstract

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, we calculate finite-temperature properties of the two dimensional Hubbard model.

4 pages, 3 figures, Revised version has expanded discussion, simplified mathematical presentation, and application to 2D Hubbard model