On Fri, May 18, 2007, adelene.stanford.edu 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.
...dac
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Wed May 23 2007 - 06:07:02 PDT
