fix bond/create

Well this is my last diagnosis. The problem is not with fix bond/create.

I am guessing this is with fix set force.

After the fix bond/create command set the velocities of the atoms in the group 4 to zero
and then use

fix 6 GP setforce 0.0 0.0 0.0


Sorry, but that also doesn’t work.

If they 1 correctly changes to type 4, since I don’t integrate for type 4, type 4 should stopped.
But the fact that type 4 still move means that Lammps doesn’t recognize the change when integrating, despite Lammps recognize the change when dumping trajectories.