Rotation of simulation box during NPT equilibrium

Dear LAMMPS users,

Hi,

During simulation of a ceramic Nanowire system, I observed that my simulation box start to rotate with respect to z axis during equilibrium. However, it stop as I use " rot yes " in " velocity all create 300 277387 ". But I think using this modification will affect my result. Any comment please ?

Thanks in advance,

Ali

Dear LAMMPS users,

Hi,

During simulation of a ceramic Nanowire system, I observed that my
simulation box start to rotate with respect to z axis during equilibrium.
However, it stop as I use " rot yes " in " velocity all create 300 277387 ".
But I think using this modification will affect my result. Any comment
please ?

it is *supposed* to have an effect.

however, the important question is:
which of the two setups is the one you desire? the one where the
system rotates, or the one where it doesn't?

and then the secondmost important question is:
what is the temperature, that you are interested in? the one with the
kinetic energy from the rotation included or the one without?

these are decisions, you have to make, and then what you have to do in
your input follows logically from that.

axel.