> You can also change the bond force constants, but TIP3P is defined
> strangely (O-H, O-H, and H-H bonds to maintain geometry), with no
> VDW parameters on the H atoms. They are designed to be treated with
> SHAKE, or you'll probably get H-atoms collapsing into a strongly
> negative atom.
Alternatively the vdw-less H might extend beyond the vdw of the O into
another negatively-charged atom, similarly leading to fusion and
divide-by-0 or other mathematical artifact. TIP3 waters are designed
to use SHAKE to avoid such problems.
Bill
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Mon Nov 22 2010 - 10:00:04 PST
