Regarding on heat output by using Nose-Hoover thermostat

[sjplimp]

The fix_modify command also adds the delta energy to the total
output potential eng (see the doc page). You can output the delta
as f_1 whether you do that or not.

Steve