Highly efficient path-integral molecular dynamics simulations with GPUMD using neuroevolution potentials: Case studies on thermal properties of materials

Abstract

Path-integral molecular dynamics (PIMD) simulations are crucial for accurately capturing nuclear quantum effects in materials. However, their computational intensity and reliance on multiple software packages often limit their applicability at large scales. Here, we present an integration of PIMD methods, including thermostatted ring-polymer molecular dynamics (TRPMD), into the open-source GPUMD package, combined with highly accurate and efficient machine-learned neuroevolution potential (NEP) models. This approach achieves almost the accuracy of first-principles calculations with the computational efficiency of empirical potentials, enabling large-scale atomistic simulations that incorporate nuclear quantum effects. We demonstrate the efficacy of the combined NEP-PIMD approach by examining various thermal properties of diverse materials, including lithium hydride (LiH), three porous metal-organic frameworks (MOFs), liquid water, and elemental aluminum. For LiH, our NEP-PIMD simulations successfully capture the isotope effect, reproducing the experimentally observed dependence of the lattice parameter on the reduced mass. For MOFs, our results reveal that achieving good agreement with experimental data requires consideration of both nuclear quantum effects and dispersive interactions. For water, our PIMD simulations capture the significant impact of nuclear quantum effects on its microscopic structure. For aluminum, the TRPMD method effectively captures thermal expansion and phonon properties, aligning well with quantum mechanical predictions. This efficient NEP-PIMD approach opens new avenues for exploring complex material properties influenced by nuclear quantum effects, with potential applications across a broad range of materials.

16 pages, 9 figures in the main text; 1 table and 8 figures in the SI

References in corpus (44)

• Van der Waals Density Functional for General Geometries
• Implementation strategies in phonopy and phono3py
• An exchange functional that tests the robustness of the plasmon description of the van der Waals density functional
• Nuclear quantum effects enter the mainstream
• Nuclear quantum effects in water
• The hiphive package for the extraction of high-order force constants by machine learning
• GPUMD: A package for constructing accurate machine-learned potentials and performing highly efficient atomistic simulations
• Efficient stochastic thermostatting of path integral molecular dynamics
• i-PI 2.0: A Universal Force Engine for Advanced Molecular Simulations
• Ab initio thermodynamics of liquid and solid water
• The Phase Diagram of a Deep Potential Water Model
• Efficient molecular dynamics simulations with many-body potentials on graphics processing units
• How to remove the spurious resonances from ring polymer molecular dynamics
• The first-principles phase diagram of monolayer nanoconfined water
• Deep Potential: a general representation of a many-body potential energy surface
• i-PI: A Python interface for ab initio path integral molecular dynamics simulations
• A foundation model for atomistic materials chemistry
• General-purpose machine-learned potential for 16 elemental metals and their alloys
• SCAN-based hybrid and double-hybrid density functionals from models without fitted parameters
• Improving the accuracy of the neuroevolution machine learning potential for multi-component systems
• Molecular dynamics simulations of heat transport using machine-learned potentials: A mini review and tutorial on GPUMD with neuroevolution potentials
• Quantum-corrected thickness-dependent thermal conductivity in amorphous silicon predicted by machine-learning molecular dynamics simulations
• DYNASOR -- A tool for extracting dynamical structure factors and current correlation functions from molecular dynamics simulations
• Isotope Effects in Liquid Water via Deep Potential Molecular Dynamics
• Efficient Anharmonic Lattice Dynamics Calculations of Thermal Transport in Crystalline and Disordered Solids
• Quantum-mechanical exploration of the phase diagram of water
• Phase transitions in inorganic halide perovskites from machine learning potentials
• i-PI 3.0: a flexible and efficient framework for advanced atomistic simulations
• 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
• Cayley modification for strongly stable path-integral and ring-polymer molecular dynamics
• 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
• Tuning the lattice thermal conductivity in van-der-Waals structures through rotational (dis)ordering
• Revealing the free energy landscape of halide perovskites: Metastability and transition characters in CsPbBr${}_3$ and MAPbI${}_3$
• Tensorial properties via the neuroevolution potential framework: Fast simulation of infrared and Raman spectra
• Combining linear-scaling quantum transport and machine-learning molecular dynamics to study thermal and electronic transports in complex materials
• Combining the D3 dispersion correction with the neuroevolution machine-learned potential
• Understanding correlations in BaZrO3:Structure and dynamics on the nano-scale
• Impact of Surface Adsorbates and Dimensionality on Templating of Halide Perovskites
• Solute segregation in polycrystalline aluminum from hybrid Monte Carlo and molecular dynamics simulations with a unified neuroevolution potential
• Isotope engineering of carrier mobility via Fröhlich electron-phonon interaction
• Size dependent solid-solid crystallization of halide perovskites
• Lattice distortion leads to glassy thermal transport in crystalline Cs${}_3$Bi${}_2$I${}_6$Cl${}_3$

Cited by in corpus (9)

• Advances in modeling complex materials: The rise of neuroevolution potentials
• PYSED: A tool for extracting kinetic-energy-weighted phonon dispersion and lifetime from molecular dynamics simulations
• NEP-MB-pol: A unified machine-learned framework for fast and accurate prediction of water's thermodynamic and transport properties
• Ultrahigh interfacial thermal conductance for cooling gallium oxide electronics using cubic boron arsenide
• GPU acceleration of ab initio simulations of large-scale identical particles based on path integral molecular dynamics
• Predicting neutron experiments from first principles: A workflow powered by machine learning
• Consistency between the Green-Kubo formula and Lorentz model for predicting the infrared dielectric function of polar materials
• Thermal conductivities of monolayer graphene oxide from machine learning molecular dynamics simulations
• Benchmarking Chemically Scalable Machine-Learning Interatomic Potentials for Large-Scale Simulations of Multicomponent Alloys