<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
  <meta content="text/html;charset=ISO-8859-7" http-equiv="Content-Type">
  <title></title>
</head>
<body bgcolor="#ffffff" text="#000066">
Hi<br>
<br>
Jussi Lehtola wrote:
<blockquote
 cite="mid:1231518641.3927.8.camel@politzer.theorphys.helsinki.fi"
 type="cite">
  <pre wrap="">Hi,


is it possible to simulate rigid molecules with GROMACS?
  </pre>
</blockquote>
Yes, rigid up to shake accuracy, but only using constraints in topol.<br>
<blockquote
 cite="mid:1231518641.3927.8.camel@politzer.theorphys.helsinki.fi"
 type="cite">
  <pre wrap="">It seems that the run option
        constraints = all-angles
does this, but using it requires SHAKE, which cannot be used in energy
minimization.
  </pre>
</blockquote>
This, except dihedrals that are not constrained by
constraints=all-angles,<br>
could lead to wrong results in case it seems to work - degrees of
freedom issue.<br>
<blockquote
 cite="mid:1231518641.3927.8.camel@politzer.theorphys.helsinki.fi"
 type="cite">
  <pre wrap="">
How can one generate suitable starting configurations of rigid
molecules? At least using a configuration minimized without any
constraints seems to crash in the first step since SHAKE cannot satisfy
the constraints in the allowed number of steps..
  </pre>
</blockquote>
Write your own utility that creates a box of rigid molecules or<br>
use a configuration from another program that handles rigid molecules
(after converting format)<br>
<blockquote
 cite="mid:1231518641.3927.8.camel@politzer.theorphys.helsinki.fi"
 type="cite">
  <pre wrap="">
Also: what is the reference used for the constraints? Is it the force
field parameters, or the starting configuration..?
  </pre>
</blockquote>
<br>
</body>
</html>