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
