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.

arXiv comment23 pages. 7 figures

References in corpus (31)

• Many-Body Physics with Ultracold Gases
• Quantum ESPRESSO toward the exascale
• Computational complexity and fundamental limitations to fermionic quantum Monte Carlo simulations
• Solving Lattice QCD systems of equations using mixed precision solvers on GPUs
• Worm Algorithm and Diagrammatic Monte Carlo: A New Approach to Continuous-Space Path Integral Monte Carlo Simulations
• Worm Algorithm for Continuous-space Path Integral Monte Carlo Simulations
• The Uniform Electron Gas at Warm Dense Matter Conditions
• The Fermion Sign Problem in Path Integral Monte Carlo Simulations: Quantum Dots, Ultracold Atoms, and Warm Dense Matter
• Establishing the Quantum Supremacy Frontier with a 281 Pflop/s Simulation
• Multi-GPU Accelerated Multi-Spin Monte Carlo Simulations of the 2D Ising Model
• First principles simulations of dense hydrogen
• Path Integral Molecular Dynamics for Bosons
• Real-space density functional theory on graphical processing units: computational approach and comparison to Gaussian basis set methods
• Quantum Simulation of an Extended Fermi-Hubbard Model Using a 2D Lattice of Dopant-based Quantum Dots
• Phaseless Auxiliary-Field Quantum Monte Carlo on Graphical Processing Units
• Path Integral Molecular Dynamics for Fermions: Alleviating the Sign Problem with the Bogoliubov Inequality
• On the thermodynamic properties of fictitious identical particles and the application to fermion sign problem
• Prediction of a supersolid phase in high-pressure deuterium
• Graphics processing units accelerated semiclassical initial value representation molecular dynamics
• 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
• Data Assimilation using a GPU Accelerated Path Integral Monte Carlo Approach
• Fast quantum Monte Carlo on a GPU
• 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
• Accelerated free energy estimation in ab initio path integral Monte Carlo simulations
• Path integral molecular dynamics for thermodynamics and Green's function of ultracold spinor bosons

Cited by in corpus (1)

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