Dear Steve and users,
I simply tested the example given in lammps(solid Ar conductivity calculation, http://lammps.sandia.gov/doc/compute_heat_flux.html) with two different variations on computer myStress line.
I removed virial to just ONLY NULL as an option. The final conductivity came out different, 0.4528 compared to with viral term, i.e., NULL virial as options. Basically, my question is how does Lammps calculate the stress term without having any terms for stress? Is this zero or something else?
I used the system temperature (solid Ar temperature) as temp ID to compare with the results with NULL. The thermal results are the same. Now, I thought that NULL could make some difference as the part of the stress term is related to kinetic energy, but it turns out that NULL does not make any difference. Do you have any ideas why?