A HOMOGENEOUS NON-EQUILIBRIUM MOLECULAR DYNAMICS METHOD FOR CALCULATING THERMAL CONDUCTIVITY WITH A THREE-BODY POTENTIAL

K.K. Mandadapu, R.E. Jones and P. Papadopoulos
J. Chem. Phys., 130, 204106-1-9, (2009)



Abstract

In this work, a homogeneous non-equilibrium 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, the 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.


(If your institution subscribes to the electronic version of the journal, click here for a copy of this article.)