Related Experiment Video
Updated: Jan 21, 2026

Measurements of CO2 Fluxes at Non-Ideal Eddy Covariance Sites
Published on: June 24, 2019
Heat Flux for Many-Body Interactions: Corrections to LAMMPS
Paul Boone1, Hasan Babaei1,2, Christopher E Wilmer1
1Department of Chemical and Petroleum Engineering , University of Pittsburgh , Pittsburgh , Pennsylvania 15261 , United States.
Abstract:
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. ( J. Chem. Phys. 2008 , 128 , 044504 ) 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 may have consequences for 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.
Related Concept Videos
Distance Corrections
Power Factor Correction
Electric Flux
Quantifying Heat
Aqueous Solutions and Heats of Hydration
When ionic compounds dissolve in water, the ions in the solid separate and disperse uniformly throughout the solution because water molecules surround and solvate the ions, reducing the strong electrostatic forces between them. This process...
Specific Heat
For example, increasing the temperature of one gram of water by 1°C requires one calorie of heat energy and can be written as 1 cal/g-°C, or...

