please reply to the mailing list not directly to me.
I think Axel has given a better suggestion on this.
Nonetheless, I want to highlight that your group command is weird, where you listed same atom type on different groups. I suggest you to review the documentation. http://lammps.sandia.gov/doc/group.html