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
Arun
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.
Makoto