Finite-size effects in canonical and grand-canonical quantum Monte Carlo simulations for fermions

Abstract

We introduce a quantum Monte Carlo method at finite temperature for interacting fermionic models in the canonical ensemble, where the conservation of the particle number is enforced. Although general thermodynamic arguments ensure the equivalence of the canonical and the grand-canonical ensembles in the thermodynamic limit, their approach to the infinite-volume limit is distinctively different. Observables computed in the canonical ensemble generically display a finite-size correction proportional to the inverse volume, whereas in the grand-canonical ensemble the approach is exponential in the ratio of the linear size over the correlation length. We verify these predictions by quantum Monte Carlo simulations of the Hubbard model in one and two dimensions in the grand-canonical and the canonical ensemble. We prove an exact formula for the finite-size part of the free energy density, energy density and other observables in the canonical ensemble and relate this correction to a susceptibility computed in the corresponding grand-canonical ensemble. This result is confirmed by an exact computation of the one-dimensional classical Ising model in the canonical ensemble, which for classical models corresponds to the so-called fixed-magnetization ensemble. Our method is useful for simulating finite systems which are not coupled to a particle bath, such as in nuclear or cold atom physics.

arXiv comment10 pages, 7 figs; v2: 11 pages, 9 figs, new results on the doped Hubbard model; v3, 12 pages, 9 figs, generalized formula for the finite-size correction in the canonical ensemble, published version

References in corpus (12)

• Inequivalence of ensembles in a system with long range interactions
• Fermionic quantum criticality in honeycomb and $\pi$-flux Hubbard models: Finite-size scaling of renormalization-group-invariant observables from quantum Monte Carlo
• Demonstration of the auxiliary-field Monte Carlo approach for sd-shell nuclei
• The ALF (Algorithms for Lattice Fermions) project release 1.0. Documentation for the auxiliary field quantum Monte Carlo code
• Ensemble inequivalence in systems with long-range interactions
• Corrections to Finite-Size Scaling in the Lattice N-Vector Model for Infinite N
• Critical adsorption and critical Casimir forces in the canonical ensemble
• Optimization of finite-size errors in finite-temperature calculations of unordered phases
• Statistical field theory with constraints: application to critical Casimir forces in the canonical ensemble
• Three-dimensional Ising model in the fixed-magnetization ensemble: a Monte Carlo study
• Continuum Limits and Exact Finite-Size-Scaling Functions for One-Dimensional $O(N)$-Invariant Spin Models
• Stabilizing Canonical-Ensemble Calculations in the Auxiliary-Field Monte Carlo Method

Cited by in corpus (17)

• Entanglement Hamiltonian of Interacting Fermionic Models
• The Semimetal-Mott Insulator Quantum Phase Transition of the Hubbard Model on the Honeycomb Lattice
• Finite Temperature Auxiliary Field Quantum Monte Carlo in the Canonical Ensemble
• The pseudogap regime in the unitary Fermi gas
• An efficient method for grand-canonical twist averaging in quantum Monte Carlo calculations
• Projection of Infinite-$U$ Hubbard Model and Algebraic Sign Structure
• Exact expressions for the partition function of the one-dimensional Ising model in the fixed-$M$ ensemble
• Mutual information in heavy-fermion systems
• Study of the superconducting order parameter in the negative-$U$ 2D-Hubbard model by grand-canonical twist-averaged boundary conditions
• A Stable, Recursive Auxiliary Field Quantum Monte Carlo Algorithm in the Canonical Ensemble: Applications to Thermometry and the Hubbard Model
• Classical and quantum liquids induced by quantum fluctuations
• The ALF (Algorithms for Lattice Fermions) project release 2.4. Documentation for the auxiliary-field quantum Monte Carlo code
• Entanglement studies of interacting fermionic models
• Finite-Size Scaling at fixed Renormalization-Group invariant
• Phase diagram of the SU($N$) antiferromagnet of spin $S$ on a square lattice
• Precision thermodynamics of the strongly interacting Fermi gas in two dimensions
• Constraint correlation functions of the one-dimensional Ising model in the scaling limit