Hi all.
I have prepared a MD simulation of a protein with solvent and counterions in an octahedron box using the same parameters set given in the amber manual for a "plain MD simulation". Here is my input file:
&cntrl
nstlim=5000000, dt=0.002, ntx=5, irest=1, ntpr=500, ntwr=5000,
ntwx=5000,
ntt=3, gamma_ln = 5.0, tempi = 300.0, temp0=300.0,
ntb=2, ntp=1, taup=2.0
ntc=2, ntf=2,
nrespa=1,
/
By looking the output file I can see that the pressure have great fluctuation, while T fluctuations , nevertheless the reported average pressure at the end of the simulation is 1.0.
Shouldn't pressure fluctuations be smaller??
I have no much experience in MD simulations, so, sorry if this is a RTFM question.
Here is part of the output file (only the line regarding pressure information
NSTEP = 500 TIME(PS) = 121.000 TEMP(K) = 300.88 PRESS = -129.1
NSTEP = 1000 TIME(PS) = 122.000 TEMP(K) = 300.23 PRESS = -130.1
NSTEP = 1500 TIME(PS) = 123.000 TEMP(K) = 299.49 PRESS = 16.1
NSTEP = 2000 TIME(PS) = 124.000 TEMP(K) = 299.41 PRESS = -124.0
NSTEP = 2500 TIME(PS) = 125.000 TEMP(K) = 300.10 PRESS = -24.4
NSTEP = 3000 TIME(PS) = 126.000 TEMP(K) = 300.66 PRESS = 40.4
NSTEP = 3500 TIME(PS) = 127.000 TEMP(K) = 299.67 PRESS = -30.6
...............................................
Thanks in advance for any answer.
yunierkis
Servicio de Correos de la Universidad Central "Marta Abreu" de Las Villas.
http://www.uclv.edu.cu
-VI Conferencia Internacional de Ciencias Empresariales, del 16 al 18 de octubre de 2008, Cayo Santa Mar?a, Cuba.
http://economia-publica.uab.es/VIcubaCICE.doc
-V Conferencia Cient?fica Internacional de Ingenier?a Mec?nica, COMEC 2008, del 4 al 6 de noviembre de 2008, UCLV, Cuba.
http://eventos.fim.uclv.edu.cu/comec
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" (in the *body* of the email)
to majordomo.scripps.edu
Received on Wed Oct 15 2008 - 05:07:31 PDT
