Thanks Charles, I wish to use a few angle constraints to define the orientation between different monomers. I guess the bondAngle term is the better option. Do I need to write this potential myself or is it already coded in an example?
best regards, Jakob Hello Jakob-- > > I wish to use angles constaints constraining three atoms, which are > not related by bonds, to a certain angle. It seems it is only > possible to constrain dihedral angles, and guess adding to the > parameter file is not a good option. > If you group three atoms together in a rigid body (using the group() method of the IVM) the angle will be preserved. Some care must be taken if you are also using torsion angle dynamics and the atoms are involved with bonds to other atoms. It is also possible to write a simple Python energy term using bondAngle.BondAngle if you'd like a few angle-only restraints If things aren't 100% clear please provide further details. best regards-- Charles
