Re: AMBER: Question about dummy atoms

From: David A. Case <case.scripps.edu>
Date: Tue, 29 Jun 2004 10:24:56 -0700

On Tue, Jun 29, 2004, Sébastien Osborne wrote:

> Hi all, I have a little question about dummy atoms used in FEP runs.
> When a hydrogen atom disappears to become a dummy atom (type "DH"), what
> should the parameters of this dummy atom be? From my understanding,
> you give it zero vdw parameters and equilibrium values equal to regular
> hydrogen atoms (for bonds, angles, torsions) with the exception that I
> gave all those parameters very small force constants. I guess my
> question is should I use the small force constants or the regular force
> constants?

I recommend using "regular" force constants. These contributions will
cancel on two sides of a thermodynamic cycle.

....good luck...dac

-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Tue Jun 29 2004 - 18:53:00 PDT
Custom Search