Re: AMBER: My model loses its motion in MD

From: Andreas Svrcek-Seiler <svrci.tbi.univie.ac.at>
Date: Thu, 24 Nov 2005 17:07:42 +0100 (CET)

Hi,
> I have been using sander for years, mostly for classical MD in PBC boxes
> using the Ewald method and never seen something that happened in my last run.
>
> A nucleic acid (RNA) model in aqueous solution, using Ewald etc. was
> subjected to MD at 300K and nothing strange happened during the first 6
> nanoseconds or so. Then, the mobility of the model started to slowly decrease
> which didn't draw my attention first, but after 11 ns, it was already very
> striking. At the end of the 20 ns run, the trajectory looked like during a
> heavily restrained MD - although no restraints were used after the initial
> equilibration 20 ns earlier. The output file showed energy and temperature
> parameters that didn't look strange at all - when you look at the entire
> trajectory it really looks like after a normal dynamics the model slowly
> stopped moving...
....What follows is pure speculation: You might have used a berendsen
thermostat with tight thermal coupling (=small tautp) *and* periodically
removed center of mass (COM) motion. If this is the case, the phenomenon
is called 'flying icecube'. If not, I honestly have no idea.

I hope that helps
best regards
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 Thu Nov 24 2005 - 16:53:00 PST
Custom Search