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

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 (PIMC) to accurately simulate warm dense matter with up to 1000 electrons (T. Dornheim et al., arXiv:2311.08098 (2023)). The inclusion of fictitious identical particles in path integral molecular dynamics (PIMD) can provide another way to simulate fermion systems. In a recent paper (J. Chem. Phys. 159, 154107 (2023)), Feldman and Hirshberg improved the recursive formula for PIMD of N identical bosons, significantly reducing the computational complexity from $O(P{N}^3)$ to $O(N^2+PN)$. 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 by using the recursive technique. As an example, we use the quadratic scaling PIMD for fictitious identical particles to study the simulation of hundreds of fermions in a two-dimensional periodic potential, in the hope of providing a simulation tool for two-dimensional Fermi-Hubbard model and other strongly correlated fermion systems, such as the simulation of ultracold fermionic gases in optical lattices.

30 pages, 8 figures

References in corpus (12)

• Many-Body Physics with Ultracold Gases
• Computational complexity and fundamental limitations to fermionic quantum Monte Carlo simulations
• Nearsightedness of Electronic Matter
• Worm Algorithm and Diagrammatic Monte Carlo: A New Approach to Continuous-Space Path Integral Monte Carlo Simulations
• Bold Diagrammatic Monte Carlo: When Sign Problem is Welcome
• Critical Temperature Curve in the BEC-BCS Crossover
• A Phaseless Auxiliary-Field Quantum Monte Carlo Perspective on the Uniform Electron Gas at Finite Temperatures: Issues, Observations, and Benchmark Study
• Towards ab initio thermodynamics of the electron gas at strong degeneracy
• From Density Response to Energy Functionals and Back: An ab initio perspective on Matter Under Extreme Conditions
• Path-integral Monte Carlo worm algorithm for Bose systems with periodic boundary conditions
• Path integral molecular dynamics simulations for Green's function in a system of identical bosons
• Path integral molecular dynamics for thermodynamics and Green's function of ultracold spinor bosons

Cited by in corpus (5)

• Normal liquid ${}^3$He studied by Path Integral Monte Carlo with a parametrized partition function
• GPU acceleration of ab initio simulations of large-scale identical particles based on path integral molecular dynamics
• Revisiting the Fermion Sign Problem from the Structure of Lee-Yang Zeros. I. The Form of Partition Function for Indistinguishable Particles and Its Zeros at 0~K
• Periodic Boundary Conditions for Bosonic Path Integral Molecular Dynamics
• A Pseudo-Fermion Propagator Approach to the Fermion Sign Problem