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:2203.09919 · The Journal of Chemical Physics

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

Yunuo Xiong, Hongwei Xiong

Published Apr 5, 2022• 12 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1063/5.0086118

Related papers

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

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

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

Yunuo Xiong, Shujuan Liu, Hongwei Xiong

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.

arXiv:2401.00274•Physical review. E•Dec 17, 2024•13 citations

Path integral molecular dynamics for thermodynamics and Green’s function of ultracold spinor bosons

Y. C. Yu, Shujuan Liu, Hongwei Xiong +1 more

Most recently, the path integral molecular dynamics has been successfully used to consider the thermodynamics of single-component identical bosons and fermions. In this work, the path integral molecular dynamics is developed to simulate thermodynamics, Green's function, and momentum distribution of two-component bosons in three dimensions. As an example of our general method, we consider the thermodynamics of up to 16 bosons in a three-dimensional harmonic trap. For noninteracting spinor bosons, our simulation shows a bump in the heat capacity. As the repulsive interaction strength increases, however, we find the gradual disappearance of the bump in the heat capacity. We believe that this simulation result can be tested by ultracold spinor bosons with optical lattices and magnetic-field Feshbach resonance to tune the inter-particle interaction. We also calculate Green's function and momentum distribution of spinor bosons. Our work facilitates the exact numerical simulation of spinor bosons, whose property is one of the major problems in ultracold Bose gases.

arXiv:2207.07653•The Journal of Chemical Physics•Jul 14, 2022•8 citations

Open quantum systems beyond equilibrium: Lindblad equation and path integral molecular dynamics

Benedikt M. Reible, Somayeh Ahmadkhani, Luigi Delle Site

The Lindblad equation determines the time evolution of the density operator of open quantum systems. While valid for any system size, its use is, in practice, restricted to prototype/surrogate models with the aim of tackling specific aspects of the overall quantum complexity of a multi-atomic system. Path integral molecular dynamics (PIMD) instead provides static and dynamical quantum statistical averages of physical observables for systems in equilibrium composed of up to thousands of atoms over timescales up to nanoseconds, under the condition that short-time quantum coherence is not relevant for the properties of interest. PIMD relies on the well-established technique of molecular dynamics (MD) with its associated classical trajectories. However, it cannot describe a direct time evolution of a system and its convergence to a stationary state in situations out of equilibrium. In this work, we analyze the link between the Lindblad equation and PIMD; specifically, we will discuss how PIMD can actually be used to calculate the time evolution of ensemble-averaged physical observables and their convergence to a stationary state for situations out of equilibrium, bypassing the need of explicitly solving the Lindblad equation. Yet, at the same time, the Lindblad equation and PIMD are linked to one another through a formal relation of equivalence, which provides an argument for the consistency of PIMD results, namely the positivity of the density operator at any time. A numerical study of a prototype system, which is of interest in chemical physics, will be used to showcase the method.

quant-pharXiv:2603.10839•Physical Review A•Mar 12, 2026•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 (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.
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.
Most recently, the path integral molecular dynamics has been successfully used to consider the thermodynamics of single-component identical bosons and fermions. In this work, the path integral molecular dynamics is developed to simulate thermodynamics, Green's function, and momentum distribution of two-component bosons in three dimensions. As an example of our general method, we consider the thermodynamics of up to 16 bosons in a three-dimensional harmonic trap. For noninteracting spinor bosons, our simulation shows a bump in the heat capacity. As the repulsive interaction strength increases, however, we find the gradual disappearance of the bump in the heat capacity. We believe that this simulation result can be tested by ultracold spinor bosons with optical lattices and magnetic-field Feshbach resonance to tune the inter-particle interaction. We also calculate Green's function and momentum distribution of spinor bosons. Our work facilitates the exact numerical simulation of spinor bosons, whose property is one of the major problems in ultracold Bose gases.
The Lindblad equation determines the time evolution of the density operator of open quantum systems. While valid for any system size, its use is, in practice, restricted to prototype/surrogate models with the aim of tackling specific aspects of the overall quantum complexity of a multi-atomic system. Path integral molecular dynamics (PIMD) instead provides static and dynamical quantum statistical averages of physical observables for systems in equilibrium composed of up to thousands of atoms over timescales up to nanoseconds, under the condition that short-time quantum coherence is not relevant for the properties of interest. PIMD relies on the well-established technique of molecular dynamics (MD) with its associated classical trajectories. However, it cannot describe a direct time evolution of a system and its convergence to a stationary state in situations out of equilibrium. In this work, we analyze the link between the Lindblad equation and PIMD; specifically, we will discuss how PIMD can actually be used to calculate the time evolution of ensemble-averaged physical observables and their convergence to a stationary state for situations out of equilibrium, bypassing the need of explicitly solving the Lindblad equation. Yet, at the same time, the Lindblad equation and PIMD are linked to one another through a formal relation of equivalence, which provides an argument for the consistency of PIMD results, namely the positivity of the density operator at any time. A numerical study of a prototype system, which is of interest in chemical physics, will be used to showcase the method.
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.