A Stable, Recursive Auxiliary Field Quantum Monte Carlo Algorithm in the Canonical Ensemble: Applications to Thermometry and the Hubbard Model

Abstract

Many experimentally-accessible, finite-sized interacting quantum systems are most appropriately described by the canonical ensemble of statistical mechanics. Conventional numerical simulation methods either approximate them as being coupled to a particle bath, or use projective algorithms which may suffer from non-optimal scaling with system size or large algorithmic prefactors. In this paper, we introduce a highly stable, recursive Auxiliary Field Quantum Monte Carlo approach that can directly simulate systems in the canonical ensemble. We apply the method to the fermion Hubbard model in one and two spatial dimensions in a regime known to exhibit a significant "sign" problem and find improved performance over existing approaches including rapid convergence to ground state expectation values. The effects of excitations above the ground state are quantified using an estimator-agnostic approach including studying the temperature dependence of the purity and overlap fidelity of the canonical and grand canonical density matrices. As an important application, we show that thermometry approaches often exploited in ultra-cold atoms that employ an analysis of the velocity distribution in the grand canonical ensemble may be subject to errors leading to an under-estimation of extracted temperatures with respect to the Fermi temperature.

14 pages, 7 figures, and Supplemental Information with key derivations. For associated data and code repository see: https://github.com/DelMaestroGroup/papers-code-CanEnsAFQMC

References in corpus (14)

• Measuring entanglement growth in quench dynamics of bosons in an optical lattice
• Path Integral Monte Carlo Simulation of the Warm-Dense Homogeneous Electron Gas
• Homogeneous Atomic Fermi Gases
• Quantum Gases in Optical Boxes
• Coupled cluster theory in materials science
• Entanglement constrained by superselection rules
• Cooling and thermometry of atomic Fermi gases
• An alternative quantum fidelity for mixed states of qudits
• Heavy deformed nuclei in the shell model Monte Carlo method
• Stable Quantum Monte Carlo Simulations for Entanglement Spectra of Interacting Fermions
• Renyi Entropies of Interacting Fermions from Determinantal Quantum Monte Carlo Simulations
• Finite-Temperature Auxiliary-Field Quantum Monte Carlo for Bose-Fermi Mixtures
• Finite-size effects in canonical and grand-canonical quantum Monte Carlo simulations for fermions
• Dynamical tuning of the chemical potential to achieve a target particle number in grand canonical Monte Carlo simulations

Cited by in corpus (6)

• Quadratic Scaling Bosonic Path Integral Molecular Dynamics
• Quantum sensing with ultracold simulators in lattice and ensemble systems: a review
• Disentangling the Physics of the Attractive Hubbard Model via the Accessible and Symmetry-Resolved Entanglement Entropies
• Exact kinetic propagators for coherent state complex Langevin simulations
• A projected complex Langevin sampling method for bosons in the canonical and microcanonical ensembles
• Fermionic-propagator and alternating-basis quantum Monte Carlo methods for correlated electrons on a lattice