On Thu, Mar 16, 2006, Ilyas Yildirim wrote:
>
> 4. When SHAKE is turned on (see ./ntc_2_ntf_2_icfe_1_ntb_1), the system
> blows up.
Have you read note (5) on p. 156 of the Users' Manual. Amber right now
does
a check that the *number* of shaken bonds is the same in both states, but
it
does not really check that the list of shaken bonds is identical (with
identical equilibrium bond lengths.) Patches to further validate the
input
files would be welcome.
>
> The reason why it does not work when SHAKE is turned ON is because of
the
> DH bonds. I think SHAKE is just constraining the H-bonds, not DH-bonds.
> DH-bonds are also like H-bonds, but the atom type for these dummy atoms
> are DH, not H or H?.
It uses the element type in LEaP (I think). So you would have to tell
LEaP
that your dummy atoms are hydrogens. But using noshakemask would be a
better
solution. It would be great if you could try this (add noshakemask, with
ntc=2) since it has not been tested very much.
...thanks!...dave
--
==================================================================
David A. Case | e-mail: case.scripps.edu
Dept. of Molecular Biology, TPC15 | fax: +1-858-784-8896
The Scripps Research Institute | phone: +1-858-784-9768
10550 N. Torrey Pines Rd. | home page:
La Jolla CA 92037 USA | http://www.scripps.edu/case
==================================================================
Received on Wed Apr 05 2006 - 23:49:37 PDT
