Re: AMBER: GB with Langevin dynamics problem 2

From: Andreas Svrcek-Seiler <svrci.tbi.univie.ac.at>
Date: Fri, 8 Apr 2005 19:55:48 +0200 (CEST)

Hi,
> The result shows constant and large translational movement of the
> molecule.
....What is large and what is constant? With Langevin-dynamics, COM
translation gets what it deserves according to the equipartition theorem
(on the average). I guess if it is a small ("light-weight") molecule, it would
travel around quite a bit due to diffusion.
COM rotation behaves like COM translation. So, periodically removing COM
translation and rotation one would expect downwards jumps in KE Trans and
KE Rot. In my opinion, one should not remove COM motion when using Langevin
dynamics. For keeping the molecule "in place", one could now and then
just re-set the COM to the origin (that's somewhat like following the
molecule's overall motion with a camera).
Alternatively, one could remove COM motion every timestep. It depends
on the way Langevin dynamics is implemented whether that's easily possible
in a physically sound way.

> Lastly, has anyone carried out long and reasonable
> Langevin simulation with Amber 8?
> ... not with amber8, but with NAB (5.0 is just out!). Three 15ns
simulations for testing purposes ran without any problem. I just didn't
remove COM motion. To be sure I just took the power spectrum of E_Kin
(from part of one run).
Result: More or less white noise with a moderate low-frequency peak as
expected. Amber8 and nab5 should give identical results (besides different
random seeds) I hope.

Of course I can't guarantee there are no bugs, but I don't think so.


I hope that helps,
Andreas

-- 
             )))))
             (((((
            ( O O )
-------oOOO--(_)--OOOo-----------------------------------------------------
               o        Wolfgang Andreas Svrcek-Seiler
               o        (godzilla)
                        svrci.tbi.univie.ac.at
       .oooO            Tel.:01-4277-52733
       (   )   Oooo. 
-------\ (----(   )--------------------------------------------------------
         \_)    ) /
               (_/
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Fri Apr 08 2005 - 19:53:00 PDT
Custom Search