Hello everyone,
I am analyzing a few recent MD runs, and I've come across a few
strange results that I wanted to get some advice on from the pros.
I ran a 1ns production over the course of 5 200ps runs on a triplex
(no water, PME used, Zn2+ counterions, 1fs timestep, 300K, 5kcal DNA
and 25 kcal ion restraints throughout the run). When I plotted the
potential energy over the course of the calculation, it stayed steady
throughout each individual run, and dropped ~100 kcal/mol between
each restart. As a result, my graph looks like a set of stairs. Is
this common, or did I foul up somewhere?
Also, I ran a 1ns production over 5 200ps runs on a duplex in a
waterbox (5500 waters, PME used, Zn2+ counterions, 2fs timestep,
300K, no restraints). The PE was fine in this case, but the RMS-
deviation seemed to start at 0 and drift gradually up to 1-1.5A after
each restart. Is this supposed to happen?
Thanks!!
Best,
Mike Keller
UVA Chemistry Dept.
Received on Thu Aug 01 2002 - 08:01:07 PDT
