GPU acceleration of ab initio simulations of large-scale identical particles based on path integral molecular dynamics

Abstract

Path integral Monte Carlo (PIMC) and path integral molecular dynamics (PIMD) provide the golden standard for the ab initio simulations of identical particles. In this work, we achieved significant GPU acceleration based on PIMD, which is equivalent to PIMC in the ab initio simulations, and developed an open-source PIMD code repository that does not rely on any other third party library. Numerical experiments show that for a system of 1600 interacting identical bosons in a harmonic trap, using a single GPU and a single CPU, it only takes two hours to achieve satisfactory simulation accuracy. With the increase of the number of identical particles, the advantage of GPU acceleration over CPU becomes more obvious, making it possible to simulate tens of thousands of identical particles from first principles using a single GPU. For example, for a system of 10000 non-interacting bosons, numerical experiments show that it takes 23 hours to obtain a simulation that is highly consistent with the exact results. Our study shows that GPU acceleration can lay a solid foundation for the wide application of PIMD simulations for extremely large-scale identical particle quantum systems with more than 10,000 particles. Numerical experiments show that a 24GB GPU can simulate up to 40000 identical particles from first principles, and the GPU acceleration leads to a roughly linear relationship between the computation time and the number of identical particles. In addition, we have also successfully implemented simulations for fictitious identical particle thermodynamics using GPU to overcome the Fermion sign problem, which makes it promising to efficiently and accurately simulate tens of thousands of fermions based on GPU.

23 pages. 7 figures

References in corpus (17)

• Many-Body Physics with Ultracold Gases
• Quantum ESPRESSO toward the exascale
• Computational complexity and fundamental limitations to fermionic quantum Monte Carlo simulations
• Worm Algorithm and Diagrammatic Monte Carlo: A New Approach to Continuous-Space Path Integral Monte Carlo Simulations
• Multi-GPU Accelerated Multi-Spin Monte Carlo Simulations of the 2D Ising Model
• First principles simulations of dense hydrogen
• On the thermodynamic properties of fictitious identical particles and the application to fermion sign problem
• Highly efficient path-integral molecular dynamics simulations with GPUMD using neuroevolution potentials: Case studies on thermal properties of materials
• Ab initio path integral Monte Carlo simulations of hydrogen snapshots at warm dense matter conditions
• On the thermodynamics of fermions at any temperature based on parametrized partition function
• From Density Response to Energy Functionals and Back: An ab initio perspective on Matter Under Extreme Conditions
• Normal liquid ${}^3$He studied by Path Integral Monte Carlo with a parametrized partition function
• Quadratic Scaling Bosonic Path Integral Molecular Dynamics
• Quadratic scaling path integral molecular dynamics for fictitious identical particles and its application to fermion systems
• 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
• Accelerated free energy estimation in ab initio path integral Monte Carlo simulations

Cited by in corpus (1)

• A Pseudo-Fermion Propagator Approach to the Fermion Sign Problem