Quadratic Scaling Bosonic Path Integral Molecular Dynamics

Abstract

Bosonic exchange symmetry leads to fascinating quantum phenomena, from exciton condensation in quantum materials to the superfluidity of liquid Helium-4. Unfortunately, path integral molecular dynamics (PIMD) simulations of bosons are computationally prohibitive beyond $\sim 100$ particles, due to a cubic scaling with the system size. We present an algorithm that reduces the complexity from cubic to quadratic, allowing the first simulations of thousands of bosons using PIMD. Our method is orders of magnitude faster, with a speedup that scales linearly with the number of particles and the number of imaginary time slices (beads). Simulations that would have otherwise taken decades can now be done in days. In practice, the new algorithm eliminates most of the added computational cost of including bosonic exchange effects, making them almost as accessible as PIMD simulations of distinguishable particles.

References in corpus (24)

• Competing quantum effects in the dynamics of a flexible water model
• Nuclear quantum effects enter the mainstream
• 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
• i-PI 2.0: A Universal Force Engine for Advanced Molecular Simulations
• Helium-4 Luttinger liquids in nanopores
• Path Integral Molecular Dynamics for Bosons
• Predicting the phase behaviors of superionic water at planetary conditions
• Inverse Temperature Dependence of Nuclear Quantum Effects in DNA Base Pairs
• Competing quantum effects in the free energy profiles and diffusion rates of hydrogen and deuterium molecules through clathrate hydrates
• Stability of Complex Biomolecular Structures: Vander Waals, Hydrogen Bond Cooperativity, and Nuclear Quantum Effects
• Interacting bosons in one dimension and Luttinger liquid theory
• Path Integral Molecular Dynamics for Fermions: Alleviating the Sign Problem with the Bogoliubov Inequality
• Attenuating the fermion sign problem in path integral Monte Carlo simulations using the Bogoliubov inequality and thermodynamic integration
• On the thermodynamic properties of fictitious identical particles and the application to fermion sign problem
• Prediction of a supersolid phase in high-pressure deuterium
• Progress and Challenges in Ab Initio Simulations of Quantum Nuclei in Weakly Bonded Systems
• Finite Temperature Auxiliary Field Quantum Monte Carlo in the Canonical Ensemble
• Theory of Non-Interacting Fermions and Bosons in the Canonical Ensemble
• Routine Molecular Dynamics Simulations Including Nuclear Quantum Effects: from Force Fields to Machine Learning Potentials
• Neural network representation of electronic structure from $ab$ $initio$ molecular dynamics
• Quantum symmetry from enhanced sampling methods
• A Stable, Recursive Auxiliary Field Quantum Monte Carlo Algorithm in the Canonical Ensemble: Applications to Thermometry and the Hubbard Model
• Path integral molecular dynamics for thermodynamics and Green's function of ultracold spinor bosons

Cited by in corpus (6)

• First principles simulations of dense hydrogen
• i-PI 3.0: a flexible and efficient framework for advanced atomistic simulations
• Quadratic scaling path integral molecular dynamics for fictitious identical particles and its application to fermion systems
• 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