Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:cond-mat/9805255 · Physical Review Letters

Quantum Monte Carlo Method for Fermions, Free of Discretization Errors

S. Rombouts, K. Heyde, N. Jachowicz

Published May 24, 1999• 83 citations
arXivHTMLPDF

Abstract

In this Letter we present a novel quantum Monte Carlo method for fermions, based on an exact decomposition of the Boltzmann operator $\mathrm{exp}(\ensuremath{-}\ensuremath{\beta}\stackrel{^}{H})$. It can be seen as a synthesis of several related methods. It has the advantage that it is free of discretization errors, and applicable to general interactions, both for ground-state and finite-temperature calculations. The decomposition is based on low-rank matrices, which allows faster calculations. As an illustration, the method is applied to an analytically solvable model (pairing in a degenerate shell) and to the Hubbard model.

DOI
10.1103/physrevlett.82.4155

Related papers

Quantum Monte Carlo Method without Hubbard-Stratonovich Transformation for Many-body Systems

J. Wang, Weiwei Pan, D. Y. Sun

By precisely writing down the matrix element of the local Boltzmann operator ($e^(-{\tau}h$), where h is the paired operator formed by two interacting terms conjugate to each other), we have proposed a new path integral formulation for quantum field theory and developed a corresponding Monte Carlo algorithm. Our formula presents remarkable outcomes such that: 1) the Hubbard-Stratonovich transformation is not necessary, and the auxiliary field is no longer needed, which can improve the computational efficiency. 2) The new formula improves the accuracy of the Suzuki-Trotter decomposition. Our method also provides a potential theoretical framework for analyzing and solving the fermion sign problem. As an example, we studied the one-dimensional half-filled Hubbard model at a finite temperature. The obtained results were in excellent agreement with the known solutions. The new formula and Monte Carlo algorithm could be used in various studies.

•Oct 17, 2021•0 citations

Efficient continuous-time quantum Monte Carlo algorithm for fermionic lattice models

Mauro Iazzi, Matthias Troyer

Efficient continuous-time quantum Monte Carlo (CT-QMC) algorithms that do not suffer from time discretization errors have become the state of the art for most discrete quantum models. They have not been widely used yet for fermionic quantum lattice models, such as the Hubbard model, nor other fermionic lattice systems due to a suboptimal scaling of $O({\ensuremath{\beta}}^{3})$ with inverse temperature $\ensuremath{\beta}$, compared to the linear scaling of discrete-time algorithms. Here we present a CT-QMC algorithm for fermionic lattice systems that matches the scaling of discrete-time methods but is more efficient and free of time discretization errors. This provides an efficient simulation scheme that is free from the systematic errors opening an avenue to more precise studies of large systems at low and zero temperature.

arXiv:1411.0683•Physical Review B•Jun 30, 2015•32 citations

Entanglement of Interacting Fermions in Quantum Monte Carlo Calculations

Tarun Grover

Given a specific interacting quantum Hamiltonian in a general spatial dimension, can one access its entanglement properties, such as the entanglement entropy corresponding to the ground state wave function? Even though progress has been made in addressing this question for interacting bosons and quantum spins, as yet there exist no corresponding methods for interacting fermions. Here we show that the entanglement structure of interacting fermionic Hamiltonians has a particularly simple form-the interacting reduced density matrix can be written as a sum of operators that describe free fermions. This decomposition allows one to calculate the Renyi entropies for Hamiltonians which can be simulated via determinantal quantum Monte Carlo calculations, while employing the efficient techniques hitherto available only for free fermions. The method presented works for the ground state, as well as for the thermally averaged reduced density matrix.

arXiv:1307.1486•Physical Review Letters•Sep 24, 2013•119 citations

Gaussian Quantum Monte Carlo Methods for Fermions and Bosons

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.

arXiv:quant-ph/0404052•Physical Review Letters•Dec 20, 2004•104 citations

The fermion Monte Carlo revisited

Roland Assaraf, Michel Caffarel, Anatole Khélif

In this work we present a detailed study of the fermion Monte Carlo algorithm (FMC), a recently proposed stochastic method for calculating fermionic ground-state energies. A proof that the FMC method is an exact method is given. In this work the stability of the method is related to the difference between the lowest (bosonic-type) eigenvalue of the FMC diffusion operator and the exact Fermi energy. It is shown that within a FMC framework the lowest eigenvalue of the new diffusion operator is no longer the bosonic ground-state eigenvalue as in standard exact diffusion Monte Carlo (DMC) schemes but a modified value which is strictly greater. Accordingly, FMC can be viewed as an exact DMC method built from a correlated diffusion process having a reduced Bose–Fermi gap. As a consequence, the FMC method is more stable than any transient method (or nodal release-type approaches). It is shown that the most recent ingredient of the FMC approach (Kalos M H and Pederiva F 2000 Phys. Rev. Lett. 85 3547), namely the introduction of non-symmetric guiding functions, does not necessarily improve the stability of the algorithm. We argue that the stability observed with such guiding functions is in general a finite-size population effect disappearing for a very large population of walkers. The counterpart of this stability is a control population error which is different in nature from the standard diffusion Monte Carlo algorithm and which is at the origin of an uncontrolled approximation in FMC. We illustrate the various ideas presented in this work with calculations performed on a very simple model having only nine states but a full ‘sign problem’. Already for this toy model it is clearly seen that FMC calculations are inherently uncontrolled.

arXiv:cond-mat/0611264•Journal of Physics A Mathematical and Theoretical•Jan 23, 2007•19 citations

Continuous-time quantum Monte Carlo method for fermions

Alexey N. Rubtsov, V. V. Savkin, Alexander I. Lichtenstein

We present a numerically exact continuous-time quantum Monte Carlo algorithm for fermions with a general interaction nonlocal in space-time. The new determinantal grand-canonical scheme is based on a stochastic series expansion for the partition function in the interaction representation. The method is particularly applicable for multiband, time-dependent correlations since it does not invoke the Hubbard-Stratonovich transformation. The test calculations for exactly solvable models, as well results for the Green function and for the time-dependent susceptibility of the multiband supersymmetric model with a spin-flip interaction are discussed.

arXiv:cond-mat/0411344•Physical Review B•Jul 18, 2005•663 citations
In this Letter we present a novel quantum Monte Carlo method for fermions, based on an exact decomposition of the Boltzmann operator \mathrm{exp}(\ensuremath{-}\ensuremath{\beta}\stackrel{^}{H}). It can be seen as a synthesis of several related methods. It has the advantage that it is free of discretization errors, and applicable to general interactions, both for ground-state and finite-temperature calculations. The decomposition is based on low-rank matrices, which allows faster calculations. As an illustration, the method is applied to an analytically solvable model (pairing in a degenerate shell) and to the Hubbard model.
By precisely writing down the matrix element of the local Boltzmann operator (e(−τhe^(-{\tau}he(−τh), where h is the paired operator formed by two interacting terms conjugate to each other), we have proposed a new path integral formulation for quantum field theory and developed a corresponding Monte Carlo algorithm. Our formula presents remarkable outcomes such that: 1) the Hubbard-Stratonovich transformation is not necessary, and the auxiliary field is no longer needed, which can improve the computational efficiency. 2) The new formula improves the accuracy of the Suzuki-Trotter decomposition. Our method also provides a potential theoretical framework for analyzing and solving the fermion sign problem. As an example, we studied the one-dimensional half-filled Hubbard model at a finite temperature. The obtained results were in excellent agreement with the known solutions. The new formula and Monte Carlo algorithm could be used in various studies.
Efficient continuous-time quantum Monte Carlo (CT-QMC) algorithms that do not suffer from time discretization errors have become the state of the art for most discrete quantum models. They have not been widely used yet for fermionic quantum lattice models, such as the Hubbard model, nor other fermionic lattice systems due to a suboptimal scaling of $O({\ensuremath{\beta}}^{3})$ with inverse temperature $\ensuremath{\beta}$, compared to the linear scaling of discrete-time algorithms. Here we present a CT-QMC algorithm for fermionic lattice systems that matches the scaling of discrete-time methods but is more efficient and free of time discretization errors. This provides an efficient simulation scheme that is free from the systematic errors opening an avenue to more precise studies of large systems at low and zero temperature.
Given a specific interacting quantum Hamiltonian in a general spatial dimension, can one access its entanglement properties, such as the entanglement entropy corresponding to the ground state wave function? Even though progress has been made in addressing this question for interacting bosons and quantum spins, as yet there exist no corresponding methods for interacting fermions. Here we show that the entanglement structure of interacting fermionic Hamiltonians has a particularly simple form-the interacting reduced density matrix can be written as a sum of operators that describe free fermions. This decomposition allows one to calculate the Renyi entropies for Hamiltonians which can be simulated via determinantal quantum Monte Carlo calculations, while employing the efficient techniques hitherto available only for free fermions. The method presented works for the ground state, as well as for the thermally averaged reduced density matrix.
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.
In this work we present a detailed study of the fermion Monte Carlo algorithm (FMC), a recently proposed stochastic method for calculating fermionic ground-state energies. A proof that the FMC method is an exact method is given. In this work the stability of the method is related to the difference between the lowest (bosonic-type) eigenvalue of the FMC diffusion operator and the exact Fermi energy. It is shown that within a FMC framework the lowest eigenvalue of the new diffusion operator is no longer the bosonic ground-state eigenvalue as in standard exact diffusion Monte Carlo (DMC) schemes but a modified value which is strictly greater. Accordingly, FMC can be viewed as an exact DMC method built from a correlated diffusion process having a reduced Bose–Fermi gap. As a consequence, the FMC method is more stable than any transient method (or nodal release-type approaches). It is shown that the most recent ingredient of the FMC approach (Kalos M H and Pederiva F 2000 Phys. Rev. Lett. 85 3547), namely the introduction of non-symmetric guiding functions, does not necessarily improve the stability of the algorithm. We argue that the stability observed with such guiding functions is in general a finite-size population effect disappearing for a very large population of walkers. The counterpart of this stability is a control population error which is different in nature from the standard diffusion Monte Carlo algorithm and which is at the origin of an uncontrolled approximation in FMC. We illustrate the various ideas presented in this work with calculations performed on a very simple model having only nine states but a full ‘sign problem’. Already for this toy model it is clearly seen that FMC calculations are inherently uncontrolled.
We present a numerically exact continuous-time quantum Monte Carlo algorithm for fermions with a general interaction nonlocal in space-time. The new determinantal grand-canonical scheme is based on a stochastic series expansion for the partition function in the interaction representation. The method is particularly applicable for multiband, time-dependent correlations since it does not invoke the Hubbard-Stratonovich transformation. The test calculations for exactly solvable models, as well results for the Green function and for the time-dependent susceptibility of the multiband supersymmetric model with a spin-flip interaction are discussed.