Finite Temperature Auxiliary Field Quantum Monte Carlo in the Canonical Ensemble

TL;DR

This paper introduces a recursive AFQMC method in the canonical ensemble that eliminates the need for chemical potential tuning, enabling more efficient and direct simulations of many-body quantum systems at fixed particle number.

Contribution

A new recursive approach for AFQMC in the canonical ensemble that does not require chemical potentials, improving efficiency and benchmarking capabilities.

Findings

Converges more quickly to the ground state than grand canonical AFQMC.

Accurately benchmarks Bose and Fermi Hubbard models.

Facilitates direct comparison with other canonical ensemble methods.

Abstract

Finite temperature auxiliary field-based Quantum Monte Carlo methods, including Determinant Quantum Monte Carlo (DQMC) and Auxiliary Field Quantum Monte Carlo (AFQMC), have historically assumed pivotal roles in the investigation of the finite temperature phase diagrams of a wide variety of multidimensional lattice models and materials. Despite their utility, however, these techniques are typically formulated in the grand canonical ensemble, which makes them difficult to apply to condensates like superfluids and difficult to benchmark against alternative methods that are formulated in the canonical ensemble. Working in the grand canonical ensemble is furthermore accompanied by the increased overhead associated with having to determine the chemical potentials that produce desired fillings. Given this backdrop, in this work, we present a new recursive approach for performing AFQMC simulations in the canonical ensemble that does not require knowledge of chemical potentials. To derive this approach, we exploit the convenient fact that AFQMC solves the many-body problem by decoupling many-body propagators into integrals over one-body problems to which non-interacting theories can be applied. We benchmark the accuracy of our technique on illustrative Bose and Fermi Hubbard models and demonstrate that it can converge more quickly to the ground state than grand canonical AFQMC simulations. We believe that our novel use of HS-transformed operators to implement algorithms originally derived for non-interacting systems will motivate the development of a variety of other methods and anticipate that our technique will enable direct performance comparisons against other many-body approaches formulated in the canonical ensemble.