GROMACS version: 4.5.5
GROMACS modification: No
Hi everybody,
I would like to continue a trajectory of gromacs from the last frame, therefore I made the state.cpt file of it. Furthermore, I made the same trajectory, this time with the full length I desired in order to understand the correct way to lengthen the first mentioned: in that case, indeed, the prosecution will be the same of the one I see in the full one. Some numbers: 1st traj stops at frame 14000, whereas the full one goes to 40000.
For several reasons, I have to extend the first trajectory via gromacs 4.5.5 (both trajectories have been generated with it, of course) and, looking at the gromacs website and also at various forums, I am doing the steps that will be listed in the following. I please ask you what is wrong, since the first frame of the new trajectory generated (i.e. the prosectuon of the 1st, to be compared to the full length one) is not the same of the 14001 of the full trajectory I run, neither the last one of the 1st trajectory (which could have been a possibility, indeed). Here are the instructions I follow:
-
grompp -c template_system_structure.gro -f mdpfile.mdp -t state.cpt -p system.top -n system.ndx -o system.tpr
-
mdrun -s system_new.tpr -t precedent_traj.trr -append
Indeed, I’ve seen that if I use the .cpt file to generate the system_new.tpr, I don’t need to state the .cpt file again in the second instruction. Moreover, I stress that I’m doing these two commands because others do not work in gromacs 4.5.5, so please don’t tell “try doing this other way ….” because I really, really have tried the whole of them.
Of course, I changed the .mdp file taking care in the prosecution, taking care of the subsequent parameters:
init_step = 10400000 (and also tried 10401000. I have dt = 0.001 ps, so that’s why the 3 additional zeros)
gen_vel = no
continuation = yes
And, just to be sure, I tried tinit = 0, 10400 and 10401 (which are correct values and do not have to be of the same order of init_step, since they are referred to the frame, not the step. Therefore 10400000 would be an enormous error). The rest of the mdp is identic to the first I used for the run of the 1st trajectory.
I tried all the combinations of init_step and tinit. I am sure the state.cpt is the last frame of the 1st traj (I checked via check -f state.cpt). Anyone knowing how to exactly set this continuation of the run? Thanks in advance!
Here is the link of one of the .mdp I tried for the prosecution: mdout.mdp - Google Drive
Best Regards,
Jacopo