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

The Journal of Chemical Physics

Improved methods for path integral Monte Carlo integration in fermionic systems

William H. Newman, Atsuo Kuki

Published Jan 15, 1992• 16 citations
Publisher

Abstract

We generalize the discretized Feynman path integral expansion by replacing the path through pure states with a path through idempotent density matrices. The transformed expression converges to the ordinary path integral, but is computationally more flexible than the ordinary form. By introducing a particular choice of these idempotent density matrices, based on rotational averaging around the two-particle center of mass, we greatly reduce the sign problem for systems of two fermions in three dimensions. In the ordinary path integral Monte Carlo approach, low temperature simulation of fermions is inefficient as the integral decays exponentially faster than the integrand (and its variance) as the temperature decreases. The new rotationally averaged algorithm dramatically retards this relative decay. Fermionic simulations of the model system of Kestner and Sinanoglu, and of the harmonic oscillator, demonstrate this improvement, as integrals calculated at temperatures much lower than the singlet–triplet splitting display negligible decay of the desired averages relative to the partition function of the modified path ensemble.

DOI
10.1063/1.462176

Related papers

Local Coarse-Grained Approximation to Path Integral Monte Carlo Integration for Fermion Systems

Deyan Sun

An approximate treatment of exchange in finite-temperature path integral Monte Carlo simulations for fermions has been proposed. In this method, some of the fine details of density matrix due to permutations have been smoothed over or averaged out by using the coarse-grained approximation. The practical usefulness of the method is tested for interacting fermions in a three dimensional harmonic well. The results show that, the present method not only reduces the sign fluctuation of the density matrix, but also avoid the fermion system collapsing into boson system at low temperatures. The method is substantiated to be exact when applied to free particles.

arXiv:0808.2329•Journal of the Physical Society of Japan•Feb 25, 2010•1 citation

Calculation of Canonical Properties of a Quantum System by Path Integral Numerical Methods

Mikhail Voznesenskiy, P. N. Vorontsov‐Velyaminov

Abstract We propose a new approach in obtaining the partition function of a system of several interacting particles (fermions) in external field within path integral Monte Carlo method based on direct averaging of the exchange contributions over the positive weight determined by non‐closed trajectories. This procedure is performed at fixed temperatures, ranging to rather low values. The complete partition function for each temperature is then obtained within an expanded ensemble procedure. We found a good agreement of data obtained by the new method with results of previously proposed path‐integral‐expanded‐ensemble Monte Carlo calculations [1, 2]. The new approach allows to reach observably lower temperatures compared to method [1, 2], thus providing a significant reduce of the sign problem. It yields an independent way to treat thermal properties of quantum systems, so the good agreement with previous data [1, 2] allows us to test and state the validity of both approaches. Simulations for systems of 2 up to 5 non‐interacting particles and particles interacting with Coulomb repulsion in 3D harmonic field were performed. The presented approach allows one to carry out calculations for low temperatures that makes it possible to extract data for the ground‐state energy and low‐temperature thermodynamics. (© 2013 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)

Contributions to Plasma Physics•May 1, 2013•0 citations

Path Integral Monte Carlo Methods for Fermions

Ethan Brown

In general, Quantum Monte Carlo methods suffer from a sign problem when simulating fermionic systems. This causes the efficiency of a simulation to decrease exponentially with the number of particles and inverse temperature. To circumvent this issue, a nodal constraint is often implemented, restricting the Monte Carlo procedure from sampling paths that cause the many-body density matrix to change sign. Unfortunately, this high-dimensional nodal surface is not a priori known unless the system is exactly solvable, resulting in uncontrolled errors. We will discuss two possible routes to extend the applicability of finite-temperatue path integral Monte Carlo. First we extend the regime where signful simulations are possible through a novel permutation sampling scheme. Afterwards, we discuss a method to variationally improve the nodal surface by minimizing a free energy during simulation. Applications of these methods will include both free and interacting electron gases, concluding with discussion concerning extension to inhomogeneous systems.

•Jan 1, 2014•0 citations

The Feynman path integral

Mark E. Tuckerman

Abstract The focus of Chapter 12 is the Feynman path integral formulations of quantum mechanics and quantum statistical mechanics. The connection between real-time quantum propagator and the density matrix is established via the Wick rotation, and a Trotter expansion is used to derive discrete path integral expressions for both. The continuous limit of the discrete path integral is derived to yield the functional integral formulation of the path integral. The construction of estimators and their averages for computing quantum observables is discussed. Next, the path integral is extended to incorporate Boltzmann, Fermi-Dirac, and Bose-Einstein statistics, and a useful recursive formula for the partition function is derived that accelerates convergence by building up N-body partition functions from those of smaller systems. Quantum free energies are discussed, and algorithms for the evaluation of path integrals using molecular dynamics or Monte Carlo are presented. The chapter concludes with a discussion of higher-order factorizations and ring-polymer contraction methods to reduce the computational overhead of path-integral calculations

•Aug 1, 2023•0 citations

Path integral molecular dynamics for fermions: Alleviating the sign problem with the Bogoliubov inequality

Barak Hirshberg, Michele Invernizzi, Michele Parrinello

We present a method for performing path integral molecular dynamics (PIMD) simulations for fermions and address its sign problem. PIMD simulations are widely used for studying many-body quantum systems at thermal equilibrium. However, they assume that the particles are distinguishable and neglect bosonic and fermionic exchange effects. Interacting fermions play a key role in many chemical and physical systems, such as electrons in quantum dots and ultracold trapped atoms. A direct sampling of the fermionic partition function is impossible using PIMD since its integrand is not positive definite. We show that PIMD simulations for fermions are feasible by employing our recently developed method for bosonic PIMD and reweighting the results to obtain fermionic expectation values. The approach is tested against path integral Monte Carlo (PIMC) simulations for up to seven electrons in a two-dimensional quantum dot for a range of interaction strengths. However, like PIMC, the method suffers from the sign problem at low temperatures. We propose a simple approach for alleviating it by simulating an auxiliary system with a larger average sign and obtaining an upper bound to the energy of the original system using the Bogoliubov inequality. This allows fermions to be studied at temperatures lower than would otherwise have been feasible using PIMD, as demonstrated in the case of a three-electron quantum dot. Our results extend the boundaries of PIMD simulations of fermions and will hopefully stimulate the development of new approaches for tackling the sign problem.

arXiv:2003.10317•The Journal of Chemical Physics•May 5, 2020•54 citations

Path-integral–expanded-ensemble Monte Carlo method in treatment of the sign problem for fermions

Mikhail Voznesenskiy, P. N. Vorontsov‐Velyaminov, Alexander P. Lyubartsev

Expanded-ensemble Monte Carlo method with Wang-Landau algorithm was used for calculations of the ratio of partition functions for classes of permutations in the problem of several interacting quantum particles (fermions) in an external field. Simulations for systems consisting of 2 up to 7 interacting particles in harmonic or Coulombic field were performed. The presented approach allows one to carry out calculations for low enough temperatures that makes it possible to extract data for the ground-state energy and low-temperature thermodynamics.

Physical Review E•Dec 3, 2009•13 citations
We generalize the discretized Feynman path integral expansion by replacing the path through pure states with a path through idempotent density matrices. The transformed expression converges to the ordinary path integral, but is computationally more flexible than the ordinary form. By introducing a particular choice of these idempotent density matrices, based on rotational averaging around the two-particle center of mass, we greatly reduce the sign problem for systems of two fermions in three dimensions. In the ordinary path integral Monte Carlo approach, low temperature simulation of fermions is inefficient as the integral decays exponentially faster than the integrand (and its variance) as the temperature decreases. The new rotationally averaged algorithm dramatically retards this relative decay. Fermionic simulations of the model system of Kestner and Sinanoglu, and of the harmonic oscillator, demonstrate this improvement, as integrals calculated at temperatures much lower than the singlet–triplet splitting display negligible decay of the desired averages relative to the partition function of the modified path ensemble.
An approximate treatment of exchange in finite-temperature path integral Monte Carlo simulations for fermions has been proposed. In this method, some of the fine details of density matrix due to permutations have been smoothed over or averaged out by using the coarse-grained approximation. The practical usefulness of the method is tested for interacting fermions in a three dimensional harmonic well. The results show that, the present method not only reduces the sign fluctuation of the density matrix, but also avoid the fermion system collapsing into boson system at low temperatures. The method is substantiated to be exact when applied to free particles.
Abstract We propose a new approach in obtaining the partition function of a system of several interacting particles (fermions) in external field within path integral Monte Carlo method based on direct averaging of the exchange contributions over the positive weight determined by non‐closed trajectories. This procedure is performed at fixed temperatures, ranging to rather low values. The complete partition function for each temperature is then obtained within an expanded ensemble procedure. We found a good agreement of data obtained by the new method with results of previously proposed path‐integral‐expanded‐ensemble Monte Carlo calculations [1, 2]. The new approach allows to reach observably lower temperatures compared to method [1, 2], thus providing a significant reduce of the sign problem. It yields an independent way to treat thermal properties of quantum systems, so the good agreement with previous data [1, 2] allows us to test and state the validity of both approaches. Simulations for systems of 2 up to 5 non‐interacting particles and particles interacting with Coulomb repulsion in 3D harmonic field were performed. The presented approach allows one to carry out calculations for low temperatures that makes it possible to extract data for the ground‐state energy and low‐temperature thermodynamics. (© 2013 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)
In general, Quantum Monte Carlo methods suffer from a sign problem when simulating fermionic systems. This causes the efficiency of a simulation to decrease exponentially with the number of particles and inverse temperature. To circumvent this issue, a nodal constraint is often implemented, restricting the Monte Carlo procedure from sampling paths that cause the many-body density matrix to change sign. Unfortunately, this high-dimensional nodal surface is not a priori known unless the system is exactly solvable, resulting in uncontrolled errors. We will discuss two possible routes to extend the applicability of finite-temperatue path integral Monte Carlo. First we extend the regime where signful simulations are possible through a novel permutation sampling scheme. Afterwards, we discuss a method to variationally improve the nodal surface by minimizing a free energy during simulation. Applications of these methods will include both free and interacting electron gases, concluding with discussion concerning extension to inhomogeneous systems.
Abstract The focus of Chapter 12 is the Feynman path integral formulations of quantum mechanics and quantum statistical mechanics. The connection between real-time quantum propagator and the density matrix is established via the Wick rotation, and a Trotter expansion is used to derive discrete path integral expressions for both. The continuous limit of the discrete path integral is derived to yield the functional integral formulation of the path integral. The construction of estimators and their averages for computing quantum observables is discussed. Next, the path integral is extended to incorporate Boltzmann, Fermi-Dirac, and Bose-Einstein statistics, and a useful recursive formula for the partition function is derived that accelerates convergence by building up N-body partition functions from those of smaller systems. Quantum free energies are discussed, and algorithms for the evaluation of path integrals using molecular dynamics or Monte Carlo are presented. The chapter concludes with a discussion of higher-order factorizations and ring-polymer contraction methods to reduce the computational overhead of path-integral calculations
We present a method for performing path integral molecular dynamics (PIMD) simulations for fermions and address its sign problem. PIMD simulations are widely used for studying many-body quantum systems at thermal equilibrium. However, they assume that the particles are distinguishable and neglect bosonic and fermionic exchange effects. Interacting fermions play a key role in many chemical and physical systems, such as electrons in quantum dots and ultracold trapped atoms. A direct sampling of the fermionic partition function is impossible using PIMD since its integrand is not positive definite. We show that PIMD simulations for fermions are feasible by employing our recently developed method for bosonic PIMD and reweighting the results to obtain fermionic expectation values. The approach is tested against path integral Monte Carlo (PIMC) simulations for up to seven electrons in a two-dimensional quantum dot for a range of interaction strengths. However, like PIMC, the method suffers from the sign problem at low temperatures. We propose a simple approach for alleviating it by simulating an auxiliary system with a larger average sign and obtaining an upper bound to the energy of the original system using the Bogoliubov inequality. This allows fermions to be studied at temperatures lower than would otherwise have been feasible using PIMD, as demonstrated in the case of a three-electron quantum dot. Our results extend the boundaries of PIMD simulations of fermions and will hopefully stimulate the development of new approaches for tackling the sign problem.
Expanded-ensemble Monte Carlo method with Wang-Landau algorithm was used for calculations of the ratio of partition functions for classes of permutations in the problem of several interacting quantum particles (fermions) in an external field. Simulations for systems consisting of 2 up to 7 interacting particles in harmonic or Coulombic field were performed. The presented approach allows one to carry out calculations for low enough temperatures that makes it possible to extract data for the ground-state energy and low-temperature thermodynamics.