PYSED: A tool for extracting kinetic-energy-weighted phonon dispersion and lifetime from molecular dynamics simulations
arXiv:2505.00353 · doi:10.1063/5.0278798
Abstract
Machine learning potential-driven molecular dynamics (MD) simulations have significantly enhanced the predictive accuracy of thermal transport properties across diverse materials. However, extracting phonon-mode-resolved insights from these simulations remains a critical challenge. Here, we introduce PYSED, a Python-based package built on the spectral energy density (SED) method, designed to efficiently compute kinetic-energy-weighted phonon dispersion and extract phonon lifetime from large-scale MD simulation trajectories. By integrating high-accuracy machine-learned neuroevolution potential (NEP) models, we validate and showcase the effectiveness of the implemented SED method across systems of varying dimensionalities. Specifically, the NEP-driven MD-SED accurately reveals how phonon modes are affected by strain in carbon nanotubes, as well as by interlayer coupling strengths and the twist angles in two-dimensional molybdenum disulfide. For three-dimensional systems, the SED method effectively establishes the thermal transport regime diagram for metal-organic frameworks, distinguishing between particlelike and wavelike propagation regions. Moreover, using bulk silicon as an example, we show that phonon SED can efficiently capture quantum dynamics based on path-integral trajectories. The PYSED package bridges MD simulations with detailed phonon-mode insights, delivering a robust tool for investigating thermal transport properties with detailed mechanisms across various materials.
16 pages in main text; 9 figures in main text, 5 figures in SI
References in corpus (12)
- A Higher-Accuracy van der Waals Density Functional
- How to remove the spurious resonances from ring polymer molecular dynamics
- An Accurate and Transferable Machine Learning Potential for Carbon
- Molecular dynamics simulations of heat transport using machine-learned potentials: A mini review and tutorial on GPUMD with neuroevolution potentials
- Mechanical and Tribological Properties of Layered Materials Under High Pressure: Assessing the Importance of Many-Body Dispersion Effects
- 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
- Advances in modeling complex materials: The rise of neuroevolution potentials
- Moiré-Driven Interfacial Thermal Transport in Twisted Transition Metal Dichalcogenides
- Highly efficient path-integral molecular dynamics simulations with GPUMD using neuroevolution potentials: Case studies on thermal properties of materials
- Dynasor 2: From Simulation to Experiment Through Correlation Functions
- Accurate Modeling of Interfacial Thermal Transport in van der Waals Heterostructures via Hybrid Machine Learning and Registry-Dependent Potentials