Hi<br /><br />I am trying to make a rigid toluene molecule (CH3C6H5) based on the Gromos <br />53a5 force field (all hydrogen except for the CH3). It's not a problem to <br />contrain all bonds, but I am having some problems with the angles. If I try <br />using:<br />&nbsp; constraints = all-angles<br />&nbsp; constraint-algorithm = Shake<br />Gromacs (ver. 3.3) halts with the line:<br />&nbsp; turning all bonds and angles into constraints...<br /><br />I have not been able the constrain angles with Shake/Lincs by removing angles. <br />At the moment, I am representing hydrogen as virtual sites, but I am still <br />having &quot;soft&quot; angles/dihedrals between the carbon atoms.<br /><br />My question is, what is the best way to make a rigid toluene with Gromacs?<br /><br />&nbsp;&nbsp; Ulf R. Pedersen<br /><BR>