Dears,
I have reruning simulations using gmx3.2 and found out some deferences
between the original and rerun energies. I found the posts made by David
Mobley and Mickael Shirts about this issue. Accordingly I downloaded
gmx3.3 and compiled it with the fixed md.c for the rerun, found in the
bugzila page (thanks to D. Mobley for the info).
The present email is about the fact that I see differences between the
original
energies and the rerun ones but in a different manner as previously
reported.
I'd like to understand the reason of this.
From this I have been reruning a piece of trajectory (first 200 ps)
generated
with gmx3.2. I rerun it on both trr and xtc files (generated with
trjconv and the
original xtc and trr files), using gmx3.3 (with the md.c fixed) and gmx3.2,
and using tpr files build with gmx3.2 with the initial coordinate file
and the final
coordinate file, and the same mdp file (including: nstenergy = 1).
The corresponding energies are in the attached file for the original
energy file
and the 8 reruns performed.
The energies indicate that:
0) I see a systematic difference with the original energies (not
significant), but using
the trr file helps getting better.
1) I do not see ANY difference between the rerun using the two tpr file
(Unexpected ?!?!)
2) I do not see significant differences between gmx3.2 and gmx3.3 (New
?!?!) in all cases.
3) I see significant differences in energies when generated with xtc and
trr files (New ?!?!).
trr files help.
All this seems in agreement with the past remarks of David and Mickael
until I compare
the energies obtained with the twp tpr files ?! Am I missing something ?
Note that my system is quite big (in cubic (!) box about 12.0, 12.0,
8.0, in nm).
Does this make sense to any of you ?
Best,
XAvier
--
----------------------------------
Xavier Periole - Ph.D.
Dept. of Biophysical Chemistry / MD Group
Univ. of Groningen
Nijenborgh 4
9747 AG Groningen
The Netherlands
Tel: +31-503634329
Fax: +31-503634398
email: [EMAIL PROTECTED]
web-page: http://md.chem.rug.nl/~periole
----------------------------------
_______________________________________________
gmx-users mailing list gmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php
