On Wed, Dec 01, 2010, Nilmeier, Jerome P. wrote:
>
> Thank you for the advice and assistance...I'll work through some of the
> equations some more. My current thought is that I can get it to work
> using either a static lambda or a gradient, as long as I'm careful about
> the bookkeeping. Perhaps I can put together a test run by patching
> some TI runs together prior to modifying the code. I'll dig into the
> technical side a bit more, and try to come up with a good plan.
If you really want to pursue this, check out the "iscale" variable in sander.
It allows "extra" dynamical variables beyond just the coordinates. Basically,
you could set iscale=1, and arrange to have your lambda value copied to
x(3*natom+1) and its force component to f(3*natom+1). The MD step will then
update these values along with the coordinates. You can see an example in how
the alignment tensor is handled for NMR residual dipolar coupling restraints.
...good luck...dac
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Thu Dec 02 2010 - 06:00:03 PST
