Hi, Fellows,

I am now trying to use lammps to calculate thermal conductivity (TC) of a metal using NEMD.

I have a question. According to the manual, TC should be calculated this way:" the total kinetic energy transferred by these swaps is computed by the fix and can be output. Dividing this quantity by time and the cross-sectional area of the simulation box yields a heat flux. The ratio of heat flux to the slope of the temperature profile is the thermal conductivity of the fluid, in appopriate units. "

I have calculated the heat flux from the output and I have calculated the temperature slope. I am wondering **whether I should divided by 2** when using the ratio of the heat flux and the slope to get TC?

Can some one who knows it tell us?

Thanks,

Shiyu