I would like to constrain all the bond lengths in a polymer while still allowing angles and dihedral rotation. I’ve looked into the shake fix but it seems it only allows small molecule constrains and the rigid fix seems it locks the entire molecule. Am I missing something? What would be the proper way to lock just bonds?
Thanks for the help,