Hey everyone,
Does anyone know how to fix the error "Per-atom viral was not tallied
on needed timestep"? I'm using free boundary conditions in the X and Y
directions and periodic boundary conditions in the Z direction for my
system. During simulation, the actual volume changes (dimensions in
the X and Y directions both decrease nonuniformly). However, by
default, LAMMPS uses for pressure calculation the total volume of the
simulation box, which is always larger than the actual volume in my
case. So I need to dump the local stress tensor for each atom to
calculate the actual stress. I've already used the command compute
stress/atom for per-atom virial. But why does the error occur?
Thank you.
Best,
Ting