Re: AMBER: Temperature Drift in MD (NAB)

From: David A. Case <>
Date: Sat, 19 May 2007 21:03:24 -0700

On Fri, May 18, 2007, wrote:

> Only atoms in the bases 19-26 were allowed to move by fitting
> restraints in the mme_init command.
> The typical input file goes as
> mm_options("cut=25.0, nsnb=100, gamma_ln=50, temp0=300, tempi=50,
> gb=1, kappa=0.19");
> mme_init(molMD, ":19-26:", "::ZZZ", dummy, NULL);
> md(3*molMD.natoms, 100000, xyz, f, v, mme);
> setmol_from_xyz(molMD, NULL, xyz);
> putpdb("test.pdb",molMD);
> >>Hi, in my MD simulation (freezing 1141 atoms, moving 255) using NAB
> >>and options shown below, the temperature drifts from the set 300K to
> >>an average of 290.5K.

As I said, I have not tested the Langevin implementation for frozen atoms;
there could be something wrong there.


The AMBER Mail Reflector
To post, send mail to
To unsubscribe, send "unsubscribe amber" to
Received on Wed May 23 2007 - 06:07:02 PDT
Custom Search