I’ve encountered a problem.
I have a rigid layer of atoms in my system, on the right of the layer there are some molecules interacting with them, while on the left of which there is a constant force to every atom in the layer by using the command:
fix 1 slab rigid group 2 slab_1 slab_2 force * off off on
fix 2 slab_1 setforce 0 0 NULL
fix 3 slab_1 aveforce 0 0 0.05
Everything goes well in the calculation, by which I mean the layer only move in a rigid body in the z direction, and every atom in the layer have the same z coordination, then I write a restart file for later use.
But when I restart the file, (and of course using the fix rigid ,setforce and aveforce), in the first frame all the atoms have the same z coordination while, as the simulation proceeds their z coordination become different, does that mean either they are not moving as a rigid body or there are force not in the z direction.
I tried to use restart2data tool to make a data file and avoid using read_restart command, but when I read the data file generated and start, nothing happened, only on screen: