© 2026 Aqora Quantum S.A.S.

Terms Privacy Legal Notice

arXiv:2401.00274 · Physical review. E

Quadratic scaling path integral molecular dynamics for fictitious identical particles and its application to fermion systems

Yunuo Xiong, Shujuan Liu, Hongwei Xiong

Published Dec 17, 2024• 13 citations

Abstract

Recently, fictitious identical particles have provided a promising way to overcome the fermion sign problem and have been used in path integral Monte Carlo to accurately simulate warm dense matter with up to 1000 electrons [T. Dornheim et al., J. Phys. Chem. Lett. 15, 1305 (2024)1948-718510.1021/acs.jpclett.3c03193]. The inclusion of fictitious identical particles in path integral molecular dynamics (PIMD) can provide another way to simulate fermion systems. In a recent paper [Y. M. Y. Feldman et al., J. Chem. Phys. 159, 154107 (2023)0021-960610.1063/5.0173749], Feldman and Hirshberg improved the recursive formula for PIMD of identical bosons, significantly reducing the computational complexity. In this paper, we extend this latest recursive formula for bosons to PIMD of fictitious identical particles to improve the efficiency of simulating fermion systems. We also provide the virial estimator for calculating energy of fictitious identical particles by using the new recursive technique, which can help suppress fluctuations in energy calculation. As an example, we use the quadratic scaling PIMD for fictitious identical particles to study the simulation of hundreds and even thousands of fermions in a two-dimensional periodic potential, in the hope of providing a simulation tool for the large-scale Fermi-Hubbard model and other strongly correlated fermion systems, such as the simulation of ultracold fermionic gases in optical lattices and the extended Fermi-Hubbard model using quantum dots.

DOI: 10.1103/physreve.110.065303

Related papers

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

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

Reweighting estimator for ab initio path integral Monte Carlo simulations of fictitious identical particles

Tobias Dornheim, Pontus Svensson, Paul Hamann +4 more

The fermion sign problem constitutes one of the most fundamental obstacles in quantum many-body theory. Recently, it has been suggested to circumvent the sign problem by carrying out path integral simulations with a fictitious quantum statistics variable ξ, which allows for a smooth interpolation between the bosonic and fermionic limits [Xiong and Xiong, J. Chem. Phys.157, 094112 (2022)]. This ξ-extrapolation method has subsequently been applied to a variety of systems and has facilitated the analysis of an x-ray scattering measurement taken at the National Ignition Facility with unprecedented accuracy [Dornheim et al., Nat. Commun. 16, 5103 (2025)]. Yet, it comes at the cost of performing an additional 10-20 simulations, which, in combination with the required small error bars, can pose a serious practical limitation. Here, we remove this bottleneck by presenting a new reweighting estimator, which allows the study of the full ξ-dependence from a single path integral Monte Carlo (PIMC) simulation. This is demonstrated for various observables of the uniform electron gas and also warm dense beryllium. We expect our study to be useful for future PIMC simulations of Fermi systems, including ultracold atoms, electrons in quantum dots, and warm dense quantum plasmas.

arXiv:2508.12323•The Journal of Chemical Physics•Oct 15, 2025•4 citations

Path integral molecular dynamics simulations for Green’s function in a system of identical bosons

Yunuo Xiong, Hongwei Xiong

Path integral molecular dynamics (PIMD) has been successfully applied to perform simulations of large bosonic systems in a recent study [Hirshberg et al., Proc. Natl. Acad. Sci. U. S. A. 116, 21445 (2019)]. In this work, we extend PIMD techniques to study Green's function for bosonic systems. We demonstrate that the development of the original PIMD method enables us to calculate Green's function and extract momentum distribution from our simulations. We also apply our method to systems of identical interacting bosons to study Berezinskii-Kosterlitz-Thouless transition around its critical temperature.

arXiv:2203.09919•The Journal of Chemical Physics•Apr 5, 2022•12 citations

Numerical calculation of Green’s function and momentum distribution for spin-polarized fermions by path integral molecular dynamics

Yunuo Xiong, Hongwei Xiong

Most recently, path integral molecular dynamics (PIMD) has been successfully applied to perform simulations of identical bosons and fermions by Hirshberg et al. In this work, we demonstrate that PIMD can be developed to calculate Green's function and extract momentum distributions for spin-polarized fermions. In particular, we show that the momentum distribution calculated by PIMD has potential applications to numerous quantum systems, e.g., ultracold fermionic atoms in optical lattices.

The Journal of Chemical Physics•May 11, 2022•10 citations

Numerical calculation of Green's function and momentum distribution for spin-polarized fermions by path integral molecular dynamics

Xiong Yunu, Hongwei Xiong

Most recently, path integral molecular dynamics (PIMD) has been successfully applied to perform simulations of identical bosons and fermions by B. Hirshberg et al.. In this work, we demonstrate that PIMD can be developed to calculate Green's function and extract momentum distribution for spin-polarized fermions. In particular, we show that the momentum distribution calculated by PIMD has potential application to numerous quantum systems, such as cold atom simulation of Mott insulator in Fermi-Hubbard model.

arXiv:2204.02401•Apr 5, 2022•0 citations

Recently, fictitious identical particles have provided a promising way to overcome the fermion sign problem and have been used in path integral Monte Carlo to accurately simulate warm dense matter with up to 1000 electrons [T. Dornheim et al., J. Phys. Chem. Lett. 15, 1305 (2024)1948-718510.1021/acs.jpclett.3c03193]. The inclusion of fictitious identical particles in path integral molecular dynamics (PIMD) can provide another way to simulate fermion systems. In a recent paper [Y. M. Y. Feldman et al., J. Chem. Phys. 159, 154107 (2023)0021-960610.1063/5.0173749], Feldman and Hirshberg improved the recursive formula for PIMD of identical bosons, significantly reducing the computational complexity. In this paper, we extend this latest recursive formula for bosons to PIMD of fictitious identical particles to improve the efficiency of simulating fermion systems. We also provide the virial estimator for calculating energy of fictitious identical particles by using the new recursive technique, which can help suppress fluctuations in energy calculation. As an example, we use the quadratic scaling PIMD for fictitious identical particles to study the simulation of hundreds and even thousands of fermions in a two-dimensional periodic potential, in the hope of providing a simulation tool for the large-scale Fermi-Hubbard model and other strongly correlated fermion systems, such as the simulation of ultracold fermionic gases in optical lattices and the extended Fermi-Hubbard model using quantum dots.

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.

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.

The fermion sign problem constitutes one of the most fundamental obstacles in quantum many-body theory. Recently, it has been suggested to circumvent the sign problem by carrying out path integral simulations with a fictitious quantum statistics variable ξ, which allows for a smooth interpolation between the bosonic and fermionic limits [Xiong and Xiong, J. Chem. Phys.157, 094112 (2022)]. This ξ-extrapolation method has subsequently been applied to a variety of systems and has facilitated the analysis of an x-ray scattering measurement taken at the National Ignition Facility with unprecedented accuracy [Dornheim et al., Nat. Commun. 16, 5103 (2025)]. Yet, it comes at the cost of performing an additional 10-20 simulations, which, in combination with the required small error bars, can pose a serious practical limitation. Here, we remove this bottleneck by presenting a new reweighting estimator, which allows the study of the full ξ-dependence from a single path integral Monte Carlo (PIMC) simulation. This is demonstrated for various observables of the uniform electron gas and also warm dense beryllium. We expect our study to be useful for future PIMC simulations of Fermi systems, including ultracold atoms, electrons in quantum dots, and warm dense quantum plasmas.

Path integral molecular dynamics (PIMD) has been successfully applied to perform simulations of large bosonic systems in a recent study [Hirshberg et al., Proc. Natl. Acad. Sci. U. S. A. 116, 21445 (2019)]. In this work, we extend PIMD techniques to study Green's function for bosonic systems. We demonstrate that the development of the original PIMD method enables us to calculate Green's function and extract momentum distribution from our simulations. We also apply our method to systems of identical interacting bosons to study Berezinskii-Kosterlitz-Thouless transition around its critical temperature.

Most recently, path integral molecular dynamics (PIMD) has been successfully applied to perform simulations of identical bosons and fermions by B. Hirshberg et al.. In this work, we demonstrate that PIMD can be developed to calculate Green's function and extract momentum distribution for spin-polarized fermions. In particular, we show that the momentum distribution calculated by PIMD has potential application to numerous quantum systems, such as cold atom simulation of Mott insulator in Fermi-Hubbard model.