|
|||||||||||||||||||||||||||||||||
AMBER Archive (2000)Subject: Decreasing PE during MD
From: Peter Gannett (pgannett_at_hsc.wvu.edu)
Dear Amber Users:
We have been doing long dynamics runs (2-3 ns) under constant pressure . During these runs everything looks pretty good for the first 1ns (temp, pressure, energy all settle down). After the first 1ns, the potential energy begins to fall at an increasing rate so that after 2ns it has dropped off about 10%. And, if you look closely at this term over 0-1ns you can see a slight downward slope. Nothing else appears to be changing (temp, pressure, volume, or kinetic energy. Can anyone tell me what is going on? In case it helps, the mdin file is given below.
Thanks,
cgcgcgcgcg, no_mod, warm it up
timlim = 999999., nmropt = 0,
ntf = 2, ntb = 2, idiel = 1, dielc = 1.0,
isftrp = 0, rwell = 0.0,
ipol = 0,
ibelly = 0, ntr = 0,
imin = 0,
nrun = 1,
temp0 = 300.0, tempi = 300.0,
tauv0 = 0.0, tauv = 0.1, vzero = 0.0,
ntp = 1, pres0 = 1.0, comp = 44.6,
ntc = 2, tol = 0.0005,
imgslt = 0, iftres = 1,
ivcap = 0, matcap = 0, fcap = 1.5,
iewald = 1,
| |||||||||||||||||||||||||||||||||
|