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:2412.06796

Introduction of Additive Particle Theory for Path Integral Approaches

Ken‐ichi Amano

Published Nov 22, 2024• 0 citations
arXivHTMLPDF

Abstract

Path integral approaches have been used for boson and fermion systems. The path integral approach has been successful in the many-boson system. However, in the many-fermion system, the path integral approach is not feasible due to the sign problem. In this letter, I introduce additive particle (AP) theory in order to generate an approximation method that avoids the sign problem. The AP theory considers one electron as a string polymer, and virtual particles are added into the system. The AP theory is an approximation, but it is constructed to be able to generate the pair distribution function between free electrons and the density of states of the free electrons at an arbitrary temperature. Hence, when the electrostatic interactions are decreased, the AP theory converges to the free electron system. On the other hand, it deviates from the actual system when the electrostatic interactions are increased. Star polymer approximation and extended star polymer approximation are also introduced.

DOI
10.48550/arxiv.2412.06796

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 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

Interacting electrons in one dimension: a path integral Monte Carlo study

Alexander P. Lyubartsev

Path integral Monte Carlo simulations have been carried out for a system of interacting fermionic particles in one dimension. Due to the fact that the sign problem can be completely eliminated for one-dimensional systems, such simulations allow one to obtain accurate energies and particle densities in arbitrary external potentials. Two cases were considered: electrons in the field of a uniform neutralizing background and in an effective field of atoms on a lattice. In the latter case, a clear asymmetry of the charge density distribution of conducting electrons and holes has been observed.

Journal of Physics A Mathematical and Theoretical•Jun 12, 2007•6 citations

Fermionic Free Energies from Ab Initio Path Integral Monte Carlo Simulations of Fictitious Identical Particles

Tobias Dornheim, Zhandos A. Moldabekov, Sebastian Schwalbe +2 more

High Resolution Image Download MS PowerPoint Slide We combine the recent η-ensemble path integral Monte Carlo approach to the free energy [Dornheim et al. Phys. Rev. B 2025 111, L041114] with a recent fictitious partition function technique based on inserting a continuous variable that interpolates between the bosonic and Fermionic limits [Xiong and Xiong J. Chem. Phys. 2022 157, 094112] to deal with the Fermion sign problem. As a practical example, we apply our setup to the warm, dense, uniform electron gas over a broad range of densities and temperatures. We obtain accurate results for the exchange–correlation free energy down to half the Fermi temperature and find excellent agreement with the state-of-the-art parametrization by Groth et al. [ Phys. Rev. Lett. 2017 119, 135001]. Our work opens up new avenues for the future study of a host of interacting Fermi systems, including warm dense matter, ultracold atoms, and electrons in quantum dots, and for Fermionic free energy calculations with unprecedented system size.

Journal of Chemical Theory and Computation•Jul 24, 2025•12 citations

Solution of the 'sign problem' in pair approximation

Alexander S. Larkin, Vladimir S. Filinov, Vladimir E. Fortov

The main difficulty for path integral Monte Carlo studies of Fermi systems results from the requirement of antisymmetrization of the density matrix and is known in literature as the 'sign problem'. To overcome this issue the new numerical version of the %{\large {\tiny {\Huge }}} Wigner approach to quantum mechanics for treatment thermodynamic properties of degenerate systems of fermions has been developed. The new path integral representation of quantum Wigner function in the phase space has been obtained for canonical ensemble. Explicit analytical expression of the Wigner function accounting for Fermi statistical effects by effective pair pseudopotential has been proposed. Derived pseudopotential depends on coordinates, momenta and degeneracy parameter of fermions and takes into account Pauli blocking of fermions in phase space. The new quantum Monte-Carlo method for calculations of average values of arbitrary quantum operators has been proposed. To test the developed approach calculations of the momentum distribution function of the degenerate ideal system of Fermi particles has been carried out in a good agreement with analytical Fermi distributions. Generalization of this approach for studies influence of interparticle interaction on momentum distribution functions of strongly coupled Fermi system is in progress.

arXiv:1702.04091•Feb 14, 2017•0 citations
Path integral approaches have been used for boson and fermion systems. The path integral approach has been successful in the many-boson system. However, in the many-fermion system, the path integral approach is not feasible due to the sign problem. In this letter, I introduce additive particle (AP) theory in order to generate an approximation method that avoids the sign problem. The AP theory considers one electron as a string polymer, and virtual particles are added into the system. The AP theory is an approximation, but it is constructed to be able to generate the pair distribution function between free electrons and the density of states of the free electrons at an arbitrary temperature. Hence, when the electrostatic interactions are decreased, the AP theory converges to the free electron system. On the other hand, it deviates from the actual system when the electrostatic interactions are increased. Star polymer approximation and extended star polymer approximation are also introduced.
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.
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.
Path integral Monte Carlo simulations have been carried out for a system of interacting fermionic particles in one dimension. Due to the fact that the sign problem can be completely eliminated for one-dimensional systems, such simulations allow one to obtain accurate energies and particle densities in arbitrary external potentials. Two cases were considered: electrons in the field of a uniform neutralizing background and in an effective field of atoms on a lattice. In the latter case, a clear asymmetry of the charge density distribution of conducting electrons and holes has been observed.
High Resolution Image Download MS PowerPoint Slide We combine the recent η-ensemble path integral Monte Carlo approach to the free energy [Dornheim et al. Phys. Rev. B 2025 111, L041114] with a recent fictitious partition function technique based on inserting a continuous variable that interpolates between the bosonic and Fermionic limits [Xiong and Xiong J. Chem. Phys. 2022 157, 094112] to deal with the Fermion sign problem. As a practical example, we apply our setup to the warm, dense, uniform electron gas over a broad range of densities and temperatures. We obtain accurate results for the exchange–correlation free energy down to half the Fermi temperature and find excellent agreement with the state-of-the-art parametrization by Groth et al. [ Phys. Rev. Lett. 2017 119, 135001]. Our work opens up new avenues for the future study of a host of interacting Fermi systems, including warm dense matter, ultracold atoms, and electrons in quantum dots, and for Fermionic free energy calculations with unprecedented system size.
The main difficulty for path integral Monte Carlo studies of Fermi systems results from the requirement of antisymmetrization of the density matrix and is known in literature as the 'sign problem'. To overcome this issue the new numerical version of the %{\large {\tiny {\Huge }}} Wigner approach to quantum mechanics for treatment thermodynamic properties of degenerate systems of fermions has been developed. The new path integral representation of quantum Wigner function in the phase space has been obtained for canonical ensemble. Explicit analytical expression of the Wigner function accounting for Fermi statistical effects by effective pair pseudopotential has been proposed. Derived pseudopotential depends on coordinates, momenta and degeneracy parameter of fermions and takes into account Pauli blocking of fermions in phase space. The new quantum Monte-Carlo method for calculations of average values of arbitrary quantum operators has been proposed. To test the developed approach calculations of the momentum distribution function of the degenerate ideal system of Fermi particles has been carried out in a good agreement with analytical Fermi distributions. Generalization of this approach for studies influence of interparticle interaction on momentum distribution functions of strongly coupled Fermi system is in progress.