Decomposition model for phonon thermal conductivity of a monatomic lattice

An analytical treatment of decomposition of the phonon thermal conductivity of a crystal with a monatomic unit cell is developed on the basis of a two-stage decay of the heat current autocorrelation function observed in molecular dynamics simulations. It is demonstrated that the contributions from the acoustic short- and long-range phonon modes to the total phonon thermal conductivity can be presented in the form of simple kinetic formulas, consisting of products of the heat capacity and the average relaxation time of the considered phonon modes as well as the square of the average phonon velocity. On the basis of molecular dynamics calculations of the heat current autocorrelation function, this treatment allows for a self-consistent numerical evaluation of the aforementioned variables. In addition, the presented analysis allows, within the Debye approximation, for the identification of the temperature range where classical molecular dynamics simulations can be employed for the prediction of phonon thermal transport properties. As a case example, Cu is considered.