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:0808.2329 · Journal of the Physical Society of Japan

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

Deyan Sun

Published Feb 25, 2010• 1 citation
arXivHTMLPDF

Abstract

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.

DOI
10.1143/jpsj.79.034002

Related papers

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

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

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 molecular dynamics approximations of quantum canonical observables

Xin Huang, Petr Plecháč, Mattias Sandberg +1 more

Mean-field molecular dynamics based on path integrals is used to approximate canonical quantum observables for particle systems consisting of nuclei and electrons. A computational bottleneck is the Monte Carlo sampling from the Gibbs density of the electron operator, which due to the fermion sign problem has a computational complexity that scales exponentially with the number of electrons. In this work, we construct an algorithm that approximates the mean-field Hamiltonian by path integrals for fermions. The algorithm is based on the determinant of a matrix with components built on Brownian bridges connecting permuted electron coordinates. The computational work for n electrons is O ( n 3 ) , which reduces the computational complexity associated with the fermion sign problem. We analyze a bias resulting from this approximation and provide a rough computational error indicator. It remains to rigorously explain the surprisingly high accuracy for high temperatures. The method becomes infeasible at low temperatures due to a large sample variance.

arXiv:2311.17333•Journal of Computational Physics•Nov 28, 2024•1 citation

Overcoming the Fermion Sign Problem in Homogeneous Systems

Jonathan L. DuBois, Ethan Brown, Berni J. Alder

Explicit treatment of many-body Fermi statistics in path integral Monte Carlo (PIMC) results in exponentially scaling computational cost due to the near cancellation of contributions to observables from even and odd permutations. Through direct analysis of exchange statistics we find that individual exchange probabilities in homogeneous systems are, except for finite size effects, independent of the configuration of other permutations present. For two representative systems, 3-He and the homogeneous electron gas, we show that this allows the entire antisymmetrized density matrix to be generated from a simple model depending on only a few parameters obtainable directly from a standard PIMC simulation. The result is a polynomial scaling algorithm and up to a 10 order of magnitude increase in efficiency in measuring fermionic observables for the systems considered.

arXiv:1409.3262•Feb 23, 2017•20 citations

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
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.
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.
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.
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.
Mean-field molecular dynamics based on path integrals is used to approximate canonical quantum observables for particle systems consisting of nuclei and electrons. A computational bottleneck is the Monte Carlo sampling from the Gibbs density of the electron operator, which due to the fermion sign problem has a computational complexity that scales exponentially with the number of electrons. In this work, we construct an algorithm that approximates the mean-field Hamiltonian by path integrals for fermions. The algorithm is based on the determinant of a matrix with components built on Brownian bridges connecting permuted electron coordinates. The computational work for n electrons is O ( n 3 ) , which reduces the computational complexity associated with the fermion sign problem. We analyze a bias resulting from this approximation and provide a rough computational error indicator. It remains to rigorously explain the surprisingly high accuracy for high temperatures. The method becomes infeasible at low temperatures due to a large sample variance.
Explicit treatment of many-body Fermi statistics in path integral Monte Carlo (PIMC) results in exponentially scaling computational cost due to the near cancellation of contributions to observables from even and odd permutations. Through direct analysis of exchange statistics we find that individual exchange probabilities in homogeneous systems are, except for finite size effects, independent of the configuration of other permutations present. For two representative systems, 3-He and the homogeneous electron gas, we show that this allows the entire antisymmetrized density matrix to be generated from a simple model depending on only a few parameters obtainable directly from a standard PIMC simulation. The result is a polynomial scaling algorithm and up to a 10 order of magnitude increase in efficiency in measuring fermionic observables for the systems considered.
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)