Hi Steve,
Thanks for your reply, we've tried the same way. Inside every loop the atom velocities are reset into zero. Then the energy vibrates instead of showing convergence.
Regards,
Xiaoyin Ji
Graduate Research Assistant
Materials Science and Engineering
North Carolina State University