Effects of cutoff functions of Tersoff potentials on molecular dynamics simulations of thermal transport

Past molecular dynamics studies of thermal transport have predominantly used Stillinger-Weber potentials. As materials continuously shrink, their properties increasingly depend on defect and surface effects. Unfortunately, Stillinger-Weber potentials are best used for diamond-cubic-like bulk crystals. They cannot represent the energies of many metastable phases, nor can they accurately predict the energetics of defective and surface regions. To study nanostructured materials, where these regions can dominate thermal transport, the accuracy of Tersoff potentials in representing these structures is more desirable. Based upon an analysis of thermal transport in a GaN system, we demonstrate that the cutoff function of the existing Tersoff potentials may lead to problems in determining the thermal conductivity. To remedy this issue, improved cutoff schemes are proposed and evaluated.