Molecular-Dynamic Simulations of Many-Particle Systems: New Faces on Old Problems

Christian F. A. Negre, Anders M. N. Niklasson, Antonio Redondo

Bradley P. Feuston

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.

K. Esselink

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.

Søren Toxværd