I met a problem in simulating the movement of a group of atoms in LAMMPS.
As is shown in the input script “in.velocity” as attached, I created a bilayer plane of Au atoms and used the “velocity” command to set all the atoms a constant velocity in the z direction. However, the average velocities of these atoms in the z direction at different times were outputed in “velocity.txt” as attached and they are not constant at all. It seems that the “velocity” command is not valid in my simulation.
Does anyone know how to solve this problem?
Thank you in advance,
in.velocity (1002 Bytes)
velocity.txt (3.62 KB)