A Homogeneous Nonequilibrium Molecular Dynamics Method for Calculating Thermal Conductivity with a Three-body Potential

In this work, Evans' homogeneous nonequilibrium molecular dynamics method for estimating thermal conductivity is extended to systems employing three-body potentials. This extension is put on a firm theoretical basis and applied to a silicon lattice modeled by the Stillinger-Weber potential. Two new methods are suggested for estimating the thermal conductivity based on a range of values of the fictitious force. Also, kinetic theory is used to estimate the linear range of the fictitious force necessary to bias the heat flow, thereby potentially reducing the number of simulations needed to estimate thermal conductivity.