[gmx-users] the -extend function of the tpbconv command

Justin A. Lemkul jalemkul at vt.edu
Wed Apr 4 02:36:31 CEST 2012

Acoot Brett wrote:
> Dear All,
> I have completed a 1ns production MD. Then I intended to extend it for 
> another 1 ns.
> I use command
> "tpbconv -s previous.tpr -extend timetoextendby -o next.tpr
> mdrun -s next.tpr -cpi previous.cpt"
> I think the results will be automatically attended to previous 
> md_0_1.tpr and md_0_1_noPBC.xtc.

Only if you then run mdrun -append and the files are appended.  tpbconv doesn't 
do any of this for you, it just creates a new run input file capable of running 
a longer simulation.  The fact that you have a file called "md_0_1_noPBC.xtc" 
seems like you've manipulated it in some way, so you shouldn't (and likely 
can't) append to it.  You should append to your original files.

> Then I use "g_rms -s md_0_1.tpr -f md_0_1_noPBC.xtc -o rmsd.xvg -tu ns" 
> to get the rmsd.xvg file. Then I use "xmgrace -nxy potential.xvg" to 
> display it. However I find for the  curve the x-axis is from 0 ns to 1 
> ns, not from 0 ns to 2 ns.

If you haven't run a new simulation, you'll only have the data that already exists.

> The definition of "-extend" of the command "tpbconv" is "Extend runtime 
> by this amount (ps)". Thus the commands I used are correct.

You haven't shown us your actual tpbconv command, so there's no way to verify 
that.  In any case, tpbconv doesn't actually extend the simulation, mdrun does.

> Will you please tell me how to process for the extension MD so that the 
> curve of the rmsd.xvg (and all the other xvg files) will have x-axis 
> from the 0 to the 2 ns (2 ns is the total MD time, including the 1 ns 
> extended part)?

There are two steps here:


You need to do both.



Justin A. Lemkul
Ph.D. Candidate
ICTAS Doctoral Scholar
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080


More information about the gromacs.org_gmx-users mailing list