Skip to main content
Open Access Publications from the University of California

UC Davis

UC Davis Previously Published Works bannerUC Davis

Force and heat current formulas for many-body potentials in molecular dynamics simulations with applications to thermal conductivity calculations


We derive expressions of interatomic force and heat current for many-body potentials such as the Tersoff, the Brenner, and the Stillinger-Weber potential used extensively in molecular dynamics simulations of covalently bonded materials. Although these potentials have a many-body nature, a pairwise force expression that follows Newton's third law can be found without referring to any partition of the potential. Based on this force formula, a stress applicable for periodic systems can be unambiguously defined. The force formula can then be used to derive the heat current formulas using a natural potential partitioning. Our heat current formulation is found to be equivalent to most of the seemingly different heat current formulas used in the literature, but to deviate from the stress-based formula derived from two-body potential. We validate our formulation numerically on various systems described by the Tersoff potential, namely three-dimensional silicon and diamond, two-dimensional graphene, and quasi-one-dimensional carbon nanotube. The effects of cell size and production time used in the simulation are examined.

Many UC-authored scholarly publications are freely available on this site because of the UC's open access policies. Let us know how this access is important for you.

Main Content
For improved accessibility of PDF content, download the file to your device.
Current View