Heat flux for many-body interactions: Corrections to LAMMPS.

The virial stress tensor-based instantaneous heat flux, which is used by LAMMPS, is only valid for the small subset of simulations that contain only pairwise interactions. For systems that contain many-body interactions using 3- or 4-body potentials, a more complete derivation is required. We have created a software patch to LAMMPS that implements the correct heat flux calculation approach for 3- and 4-body potentials, based on the derivation by Torii et al. Using two example systems, the error in the uncorrected code for many-body potential heat flux is shown to be significant and reaches nearly 100% of the many-body potential heat flux for the systems we studied; hence, the error of the total heat flux calculation is proportional to the fraction of the total heat flux transferred through the many-body potentials. This error in may have consequences for calculating thermal conductivities calculated using the Green-Kubo method or any NEMD method that uses the instantaneous heat flux. We recommend that all researchers using LAMMPS for heat flux calculations where significant heat is transferred via the many-body potentials adopt the corrected code.

[1]  C. Wilmer,et al.  Thermal Transport in Interpenetrated Metal–Organic Frameworks , 2018 .

[2]  K. Pipe,et al.  Thermal Conductance in Cross-linked Polymers: Effects of Non-Bonding Interactions. , 2017, The journal of physical chemistry. B.

[3]  A. Henry,et al.  Phonon transport in amorphous carbon using Green-Kubo modal analysis , 2016 .

[4]  B. Hong,et al.  Length-dependent thermal conductivity in suspended single-layer graphene , 2014, Nature Communications.

[5]  Marcus D. Hanwell,et al.  Avogadro: an advanced semantic chemical editor, visualization, and analysis platform , 2012, Journal of Cheminformatics.

[6]  T. Ohara,et al.  Heat conduction in chain polymer liquids: molecular dynamics study on the contributions of inter- and intramolecular energy transfer. , 2011, The Journal of chemical physics.

[7]  Aidan P Thompson,et al.  General formulation of pressure and stress tensor for arbitrary many-body interaction potentials under periodic boundary conditions. , 2009, The Journal of chemical physics.

[8]  José Mario Martínez,et al.  PACKMOL: A package for building initial configurations for molecular dynamics simulations , 2009, J. Comput. Chem..

[9]  T. Ohara,et al.  Contribution of inter- and intramolecular energy transfers to heat conduction in liquids. , 2008, The Journal of chemical physics.

[10]  Juan J. de Pablo,et al.  ON THE SIMULATION OF VAPOR-LIQUID EQUILIBRIA FOR ALKANES , 1998 .

[11]  F. Müller-Plathe.  A simple nonequilibrium molecular dynamics method for calculating the thermal conductivity , 1997 .

[12]  Steve Plimpton,et al.  Fast parallel algorithms for short-range molecular dynamics , 1993 .

[13]  R. Hardy,et al.  Formulas for determining local properties in molecular‐dynamics simulations: Shock waves , 1982 .

[14]  Melville S. Green,et al.  Markoff Random Processes and the Statistical Mechanics of Time‐Dependent Phenomena. II. Irreversible Processes in Fluids , 1954 .

[15]  J. Kirkwood,et al.  The Statistical Mechanical Theory of Transport Processes. IV. The Equations of Hydrodynamics , 1950 .