Quantitatively Analyzing Phonon Spectral Contribution of Thermal Conductivity Based on Non-Equilibrium Molecular Dynamics Simulation I: From Space Fourier Transform

TL;DR

This paper introduces a new method called TDDDM for analyzing phonon spectral contributions to thermal conductivity directly from NEMD simulations, enabling size effect studies that were previously difficult.

Contribution

The paper presents TDDDM, a novel approach based on NEMD and lattice dynamics, to predict phonon mode specific thermal conductivity, validated against existing methods.

Findings

TDDDM accurately predicts phonon contributions in LJ Argon and Si.

TDDDM reveals size effects by showing phonons with mean free path larger than system size contribute less.

The method confirms the physical origin of size effects in NEMD thermal conductivity predictions.

Abstract

Probing detailed spectral dependence of phonon transport properties in bulk materials is critical to improve the function and performance of structures and devices in a diverse spectrum of technologies. Currently, such information can only be provided by the phonon spectral energy density (SED) or equivalently time domain normal mode analysis (TDNMA) methods in the framework of equilibrium molecular dynamics simulation (EMD), but has not been realized in non-equilibrium molecular dynamics simulations (NEMD) so far. In this paper we generate a new scheme directly based on NEMD and lattice dynamics theory, called time domain direct decomposition method (TDDDM), to predict the phonon mode specific thermal conductivity. Two benchmark cases of Lennard-Jones (LJ) Argon and Stillinger-Weber (SW) Si are studied by TDDDM to characterize contributions of individual phonon modes to overall thermal conductivity and the results are compared with that predicted using SED and TDNMA. Excellent agreements are found for both cases, which confirm the validity of our TDDDM approach. The biggest advantage of TDDDM is that it can be used to investigate the size effect of individual phonon modes in NEMD simulations, which cannot be tackled by SED and TDNMA in EMD simulations currently. We found that the phonon modes with mean free path larger than the system size are truncated in NEMD and contribute little to the overall thermal conductivity. The TDDDM provides direct physical origin for the well-known strong size effects in thermal conductivity prediction by NEMD.