Question about thermo output frequency

Hello all

I’m trying to calculate the thermal conductivity using M-P method with lammps reax/c. I followed the KAPPA example, like

variable kb equal 1.3806504e-23
variable avogadro equal 6.02214e23
variable kCal2J equal 4186.0/${avogadro}

pair_style reax/c NULL
pair_coeff * * ffield O Si

You don't have tdiff calculated at step 0.


29.10.2014, 00:58, "JeJoon Yeon" <[email protected]>:

At step 10000, to be more correct.

29.10.2014, 18:04, "Oleg Sergeev" <[email protected]>:

Thank you. Could you give me brief explanation regarding this? How can I escape from this situation?

The brief explanation is that your f2 fix accumulates values for 1000 steps before calculating average (that you then use to find tdiff). So, don’t print tdiff before step 11000.


JeJoon Yeon <[email protected]…24…> 31 октября 2014 г. 0:00:19 написал: