Re: AMBER: Implicit MD problem

From: Carlos Simmerling <carlos.ilion.bio.sunysb.edu>
Date: Wed, 14 Jul 2004 09:42:18 -0400

It is very hard to help without seeing the actual files that you used.
If your batches were 1ns, why does this input file say nstlim=500000?
How did you get it to stop at exactly 1ns?
Also, please send the scripts so we can see how you did the restarting
and what was used for the refc for restraints, etc.
Please include them in the mail as plain text, since I do not open Word
documents in email.

  ----- Original Message -----
  From: yen li
  To: amber.scripps.edu
  Sent: Wednesday, July 14, 2004 7:01 AM
  Subject: AMBER: Implicit MD problem


  Hello amber users,
  I am doing the md implicit solvent md simulations on a small protein by putting restraints on the helices, using amber7 suite of programs. The input file for the production phase is given below

   &cntrl
          imin=0,
          ntx=5, irest=1,
          ntwx=5000, ntpr=200,
          ntf=2, ntb=0,dielc=4,cut=300.0,
          ntr=1, nscm=200,
          nstlim=500000,dt=0.002,temp0=300,tempi=300,
          ntc=2,ntt=1,tautp=2.0
   &end
   &ewald
          eedmeth=5,
   &end
  restraint 1
   25.0
  RES 9 17
  RES 22 27
  RES 32 43
  END
  END

  I am doing the simulation in baches of 1 ns each.The problem is that the temp v/s runlength, total energy v/s runlength and pot energy v/s run length plots shows kinks or drops at the points where i restart the next batch of the production phase. The plots for the 16 ns are attached.
  Is there any problem in the restarting of the simulation.


------------------------------------------------------------------------------
  Do you Yahoo!?
  Take Yahoo! Mail with you! Get it on your mobile phone.
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Wed Jul 14 2004 - 14:53:01 PDT
Custom Search