I am learning to calculate thermal conductivity of silicon using Green-Kubo method, but my result doesn’t agree with literature values. When I track back the source of the problem, a question about how LAMMPS calculates heat flux comes. In the manual about the command compute stress/atom, the stress tensor equation only includes pair-wise energy contribution as a second term. Does this term include three-body interaction like in Stillinger-Webber or Tersoff potential? Thanks.