Efficient molecular dynamics simulations with many-body potentials on graphics processing units

Abstract

Graphics processing units have been extensively used to accelerate classical molecular dynamics simulations. However, there is much less progress on the acceleration of force evaluations for many-body potentials compared to pairwise ones. In the conventional force evaluation algorithm for many-body potentials, the force, virial stress, and heat current for a given atom are accumulated within different loops, which could result in write conflict between different threads in a CUDA kernel. In this work, we provide a new force evaluation algorithm, which is based on an explicit pairwise force expression for many-body potentials derived recently [Phys. Rev. B 92 (2015) 094301]. In our algorithm, the force, virial stress, and heat current for a given atom can be accumulated within a single thread and is free of write conflicts. We discuss the formulations and algorithms and evaluate their performance. A new open-source code, GPUMD, is developed based on the proposed formulations. For the Tersoff many-body potential, the double precision performance of GPUMD using a Tesla K40 card is equivalent to that of the LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator) molecular dynamics code running with about 100 CPU cores (Intel Xeon CPU X5670 @ 2.93 GHz).

8 pages, 4 figures

References in corpus (5)

• Strong scaling of general-purpose molecular dynamics simulations on GPUs
• Thermal conductivity decomposition in two-dimensional materials: Application to graphene
• Amorphized graphene: A stiff material with low thermal conductivity
• Multiscale modeling of polycrystalline graphene: A comparison of structure and defect energies of realistic samples from phase field crystal models
• Accelerated molecular dynamics force evaluation on graphics processing units for thermal conductivity calculations

Cited by in corpus (66)

• Neuroevolution machine learning potentials: Combining high accuracy and low cost in atomistic simulations and application to heat transport
• GPUMD: A package for constructing accurate machine-learned potentials and performing highly efficient atomistic simulations
• Homogeneous nonequilibrium molecular dynamics method for heat transport and spectral decomposition with many-body potentials
• Modeling heat transport in crystals and glasses from a unified lattice-dynamical approach
• Influence of Boundaries and Thermostatting on Nonequilibrium Molecular Dynamics Simulations of Heat Conduction in Solids
• General-purpose machine-learned potential for 16 elemental metals and their alloys
• Improving the accuracy of the neuroevolution machine learning potential for multi-component systems
• N-Graphdiyne two-dimensional nanomaterials: Semiconductors with low thermal conductivity and high stretchability
• Heat conduction theory including phonon coherence
• A revisit to phonon-phonon scattering in single-layer graphene
• Molecular dynamics simulations of heat transport using machine-learned potentials: A mini review and tutorial on GPUMD with neuroevolution potentials
• DeePMD-kit v3: A Multiple-Backend Framework for Machine Learning Potentials
• Equivalence of the equilibrium and the nonequilibrium molecular dynamics methods for thermal conductivity calculations: From bulk to nanowire silicon
• Spectral Decomposition of Thermal Conductivity: Comparing Velocity Decomposition Methods in Homogeneous Molecular Dynamics Simulations
• Quantum-corrected thickness-dependent thermal conductivity in amorphous silicon predicted by machine-learning molecular dynamics simulations
• Thermal Transport in MoS${}_2$ from Molecular Dynamics using Different Empirical Potentials
• Anisotropic and high thermal conductivity in monolayer quasi-hexagonal fullerene: A comparative study against bulk phase fullerene
• Kapitza thermal resistance across individual grain boundaries in graphene
• Magic angle in thermal conductivity of twisted bilayer graphene
• Sub-micrometer phonon mean free paths in metal-organic frameworks revealed by machine-learning molecular dynamics simulations
• Accurate prediction of heat conductivity of water by a neuroevolution potential
• Mechanisms of temperature-dependent thermal transport in amorphous silica from machine-learning molecular dynamics
• Limits of the phonon quasi-particle picture at the cubic-to-tetragonal phase transition in halide perovskites
• A generalized decay law for particle- and wave-like thermal phonons
• Reduced Thermal Conductivity of Supported and Encased Monolayer and Bilayer MoS${}_2$
• Tuning the lattice thermal conductivity in van-der-Waals structures through rotational (dis)ordering
• Bimodal grain-size scaling of thermal transport in polycrystalline graphene from large-scale molecular dynamics simulations
• Heat transport in pristine and polycrystalline single-layer hexagonal boron nitride
• Correcting force error-induced underestimation of lattice thermal conductivity in machine learning molecular dynamics
• Ultrahigh convergent thermal conductivity of carbon nanotubes from comprehensive atomistic modeling
• Thermal transport properties of single-layer black phosphorous from extensive molecular dynamics simulations
• Thermal and electronic transport characteristics of highly stretchable graphene kirigami
• Thermal conductivity reduction in carbon nanotube by fullerene encapsulation: A molecular dynamics study
• Advances in modeling complex materials: The rise of neuroevolution potentials
• Revealing the free energy landscape of halide perovskites: Metastability and transition characters in CsPbBr${}_3$ and MAPbI${}_3$
• Highly efficient path-integral molecular dynamics simulations with GPUMD using neuroevolution potentials: Case studies on thermal properties of materials
• Phonon Confinement and Transport in Ultrathin Films
• Insight into the effect of force error on the thermal conductivity from machine-learned potentials
• Combining linear-scaling quantum transport and machine-learning molecular dynamics to study thermal and electronic transports in complex materials
• NEP-MB-pol: A unified machine-learned framework for fast and accurate prediction of water's thermodynamic and transport properties
• Effect of the accuracy of interatomic force constants on the prediction of lattice thermal conductivity
• Variable thermal transport in black, blue, and violet phosphorene from extensive atomistic simulations with a neuroevolution potential
• Mode Localization and Suppressed Heat Transport in Amorphous Alloys
• Inter-layer and Intra-layer Heat Transfer in Bilayer/Monolayer Graphene van der Waals Heterostructure: Is There a Kapitza Resistance Analogous?
• Combining the D3 dispersion correction with the neuroevolution machine-learned potential
• Dissimilar thermal transport properties in $\kappa$-Ga${}_2$O${}_3$ and $\beta$-Ga${}_2$O${}_3$ revealed by machine-learning homogeneous nonequilibrium molecular dynamics simulations
• Understanding correlations in BaZrO3:Structure and dynamics on the nano-scale
• A minimal Tersoff potential for diamond silicon with improved descriptions of elastic and phonon transport properties
• Heat transport across graphene/hexagonal-BN tilted grain boundaries from phase-field crystal model and molecular dynamics simulations
• Interpretation of apparent thermal conductivity in finite systems from equilibrium molecular dynamics simulations
• Impact of Surface Adsorbates and Dimensionality on Templating of Halide Perovskites
• Anharmonicity of the antiferrodistortive soft mode in barium zirconate BaZrO${}_3$
• GPUQT: An efficient linear-scaling quantum transport code fully implemented on graphics processing units
• Materials Informatics: Emergence To Autonomous Discovery In The Age Of AI
• Exactly equivalent thermal conductivity in finite systems from equilibrium and nonequilibrium molecular dynamics simulations
• Untangling the Raman spectrum of cubic and tetragonal BaZrO${}_3$
• Ultrahigh interfacial thermal conductance for cooling gallium oxide electronics using cubic boron arsenide
• Grain extraction and microstructural analysis method for two-dimensional poly and quasicrystalline solids
• Accelerating the Discovery of Materials with Expected Thermal Conductivity via a Synergistic Strategy of DFT and Interpretable Deep Learning
• A neural network-based four-body potential energy surface for parahydrogen
• Modular hybrid machine learning and physics-based potentials for scalable modeling of van der Waals heterostructures
• Predicting the Thermal Conductivity Collapse in SWCNT Bundles: The Interplay of Symmetry Breaking and Scattering Revealed by Machine-Learning-Driven Quantum Transport
• Substrate-Dependence of Monolayer MoS${}_2$ Thermal Conductivity and Thermal Boundary Conductance
• Optimizing AIREBO: Navigating the Journey from Complex Legacy Code to High Performance
• Interface phonon modes governing the ideal limit of thermal transport across diamond/cubic boron nitride interfaces
• A deep learning approach to the structural analysis of proteins