Re: UNABLE TO FIND ANGLE PARAMETERS FOR HT HT OT

From: Peter Freddolino (petefred_at_ks.uiuc.edu)
Date: Thu Dec 22 2005 - 14:35:37 CST

Hi Shirley,
just for the list's sake, I'll point out to everyone that the problem
with the CN7 atoms that Shirley originally posted is just because the
angle consisting of 3 CN7 angles is not in the standard parameter file;
it doesn't matter if you've used CN7 atoms before because each angle is
parameterized differently.

The issue with the HT HT OT angle that you're experiencing means that
the waters that you have contain an extra bond between the two
hydrogens. This is present in some topology files and is a legacy of
SHAKE and charmm, which needed this bond to accurately maintain the
planarity of water molecules. This bond is not needed anymore for NAMD,
and I generally advise people to remove it. How did you generate your
psf file? The easiest way to do this is usually to modify your topology
file and then regenerate your psf. This is done by changing the RESI
TIP3 entry in your topology file to read

RESI TIP3 0.000 ! tip3p water model, generate using noangle nodihedral
GROUP
ATOM OH2 OT -0.834
ATOM H1 HT 0.417
ATOM H2 HT 0.417
!BOND OH2 H1 OH2 H2 H1 H2 ! the last bond is needed for shake
BOND OH2 H1 OH2 H2
ANGLE H1 OH2 H2 ! required
ACCEPTOR OH2
PATCHING FIRS NONE LAST NONE

(This is the form now found in the namd tutorial).

Best,
Peter

This archive was generated by hypermail 2.1.6 : Wed Feb 29 2012 - 15:41:27 CST