Efficient perturbation-tracking method for directly probing the spectral phonon properties from molecular dynamics...

Area of Science:

• Computational Materials Science
• Condensed Matter Physics
• Statistical Mechanics

Background:

• Traditional methods for spectral phonon property extraction from molecular dynamics (MD) simulations, such as normal mode analysis (NMA) and spectral energy density analysis, demand extensive simulation times for convergence.
• The computational expense of these existing methods limits their application to small systems and empirical potentials, hindering detailed investigations.
• Accurate phonon anharmonic frequencies and lifetimes are crucial for understanding thermal transport properties.

Purpose of the Study:

• To develop a novel, efficient, and accurate method for directly probing mode-wise phonon anharmonic frequencies and lifetimes from MD simulations.
• To overcome the computational limitations of existing techniques, enabling applications to larger systems and more demanding simulations like ab initio MD.
• To provide a straightforward and easily implementable approach for analyzing phonon dynamics.

Main Methods:

• Introduction of the perturbation-tracking (PT) method for direct phonon property analysis.
• Validation of the PT method using silicon (Si) as a model system, comparing results with conventional normal mode analysis (NMA).
• Assessment of computational efficiency and accuracy improvements over existing NMA approaches.

Main Results:

• The PT method yields results in excellent agreement with conventional NMA.
• PT demonstrates significant efficiency gains, requiring two to three orders of magnitude less simulation time for comparable accuracy in phonon lifetime determination.
• The method preserves phonon mode dynamics, allowing for the study of transient thermal transport in non-equilibrium systems.

Conclusions:

• The perturbation-tracking method offers a highly efficient and accurate alternative for extracting spectral phonon properties from MD simulations.
• Its computational advantages make it suitable for large-scale systems and ab initio MD simulations, broadening the scope of phonon transport studies.
• The PT method's simplicity and flexibility in selecting specific phonon modes enhance its utility for diverse research applications.