Re: [AMBER] question about scmask

From: <steinbrt.rci.rutgers.edu>
Date: Fri, 6 May 2011 04:23:28 -0400 (EDT)

Hi,

> I agree with your comment. I think what puzzles me is the scmask
> keyword. For example if I use scmask='x,y' for state 0 and scmask='x'
> for state 1, or use scmask='y' for state 0 and scmask=' ' for state 1,
> assuming everything else is the same, am I going to get the same dv/dl
> for each transformation step?

correct, with sufficient sampling that should be true. But in all normal
cases, including the one we discuss here, everything else will *not* be
the same, since the transformation including x in the scmasks will start
from a molecule with different atomic partial charges than the
transformation in which x is a common atom (namely differing in the charge
on x), leading to a different conformational ensemble.

Kind Regards,

Dr. Thomas Steinbrecher
formerly at the
BioMaps Institute
Rutgers University
610 Taylor Rd.
Piscataway, NJ 08854

_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Fri May 06 2011 - 01:30:04 PDT
Custom Search