I was trying to use CG method to relax a structure. My starting configuration is the last snapshot from a MD run.
I know CG method shouldn't do anything with the initial velocity. But looks like LAMMPS read in velocities and used them.
What I observed is the velocities remains constant during CG relaxation, and total energy simply fluctuates instead of droping.
Does this mean a bug in LAMMPS?