<a name="vel">When we mention<br>gen_vel=no ;And provide pdb as input with no velocities<br></a><a name="vel"><br>As mentioned in the manual:</a><a name="vel"><br>"The velocities are set to zero
when there are no velocities in the input structure file"<br><br></a>Please elaborate what does this sentence mean.<br><br><div class="gmail_quote">On Fri, Apr 6, 2012 at 19:10, Justin A. Lemkul <span dir="ltr"><<a href="mailto:jalemkul@vt.edu">jalemkul@vt.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="im"><br>
<br>
bipin singh wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Also, if we give continuation=yes in mdp file and use input as pdb file as input instead of gro file, grompp never complains....I don't no how it reads velocities from pdb file (as no velocities are present in pdb files). Ideally it should complain that no velocities found in input file....<br>
<br>
</blockquote>
<br></div>
Again, the "continuation" keyword has nothing to do with velocities. If you have "gen_vel = no" and you provide a .pdb file with no velocities, then you do not preserve velocity information. Whatever the initial forces are govern the resulting motions. In any case, you do not preserve the previous simulation conditions.<br>
<br>
-Justin<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5">
On Fri, Apr 6, 2012 at 12:42, Peter C. Lai <<a href="mailto:pcl@uab.edu" target="_blank">pcl@uab.edu</a> <mailto:<a href="mailto:pcl@uab.edu" target="_blank">pcl@uab.edu</a>>> wrote:<br>
<br>
On 2012-04-06 11:05:51AM +0400, James Starlight wrote:<br>
> Dear Gromacs users!<br>
><br>
><br>
> I have small question about order of the runs and input data.<br>
><br>
> Ussually I do 2 equilibration phases and subsequent productive<br>
phase in the<br>
> conditions wich are equal to the last equilibration phase ( e.g<br>
often this<br>
> is npt ).<br>
><br>
> In the second equil.mdp and md.mdp there is option<br>
><br>
> continuation = yes<br>
><br>
> which means that there have been previous phases of the<br>
simulation from<br>
> wich coordinates and velocities should be taken.<br>
><br>
> As I understood the coordinates is taken from .gro file but from<br>
what file<br>
> the velocities must be providen ? Does it .cpt checkpoint file from<br>
> previous run? In some cases I've forgotten to define -t npt.cpt<br>
for my MD<br>
> run providing only coordinates in GRO file, topology and md.mdp<br>
but I have<br>
> not seen any errors in such simulation due to absence of that<br>
.cpt and<br>
> GROMPP never remind me of the absense of this file. What exactly<br>
is in that<br>
> .cpt file and from wich source the velocities from equilibration<br>
phase are<br>
> taken ?<br>
<br>
continuation = yes is telling LINCS that it is a continuation and<br>
it should not attempt to refit the constrained bonds on the first pass.<br>
<br>
The coords, velocities, state, and box information are taken either<br>
from the<br>
cpt file or you can specify the previous .trr and it will take the<br>
last frame<br>
from that and use it. If you used the output gro file without -t, then<br>
it will take coordinates and velocities from the .gro but the<br>
problem there<br>
is the limited precision (3 decimal points for each).<br>
<br>
--<br>
==============================<u></u>==============================<u></u>======<br>
Peter C. Lai | University of Alabama-Birmingham<br>
Programmer/Analyst | KAUL 752A<br>
Genetics, Div. of Research | 705 South 20th Street<br></div></div>
<a href="mailto:pcl@uab.edu" target="_blank">pcl@uab.edu</a> <mailto:<a href="mailto:pcl@uab.edu" target="_blank">pcl@uab.edu</a>> | Birmingham AL<div class="im"><br>
35294-4461<br>
(205) 690-0808 |<br>
==============================<u></u>==============================<u></u>======<br>
<br>
--<br>
gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br></div>
<mailto:<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>><div class="im"><br>
<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/<u></u>mailman/listinfo/gmx-users</a><br>
Please search the archive at<br>
<a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/<u></u>Support/Mailing_Lists/Search</a> before posting!<br>
Please don't post (un)subscribe requests to the list. Use the<br>
www interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a><br></div>
<mailto:<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@<u></u>gromacs.org</a>>.<div class="im"><br>
Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/<u></u>Support/Mailing_Lists</a><br>
<br>
<br>
<br>
<br>
-- <br>
-----------------------<br></div>
/Regards,/<br>
Bipin Singh<br>
<br><span class="HOEnZb"><font color="#888888">
</font></span></blockquote><span class="HOEnZb"><font color="#888888">
<br>
-- <br>
==============================<u></u>==========<br>
<br>
Justin A. Lemkul<br>
Ph.D. Candidate<br>
ICTAS Doctoral Scholar<br>
MILES-IGERT Trainee<br>
Department of Biochemistry<br>
Virginia Tech<br>
Blacksburg, VA<br>
jalemkul[at]<a href="http://vt.edu" target="_blank">vt.edu</a> | (540) 231-9080<br>
<a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin" target="_blank">http://www.bevanlab.biochem.<u></u>vt.edu/Pages/Personal/justin</a></font></span><div class="HOEnZb"><div class="h5"><br>
<br>
==============================<u></u>==========<br>
-- <br>
gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/<u></u>mailman/listinfo/gmx-users</a><br>
Please search the archive at <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/<u></u>Support/Mailing_Lists/Search</a> before posting!<br>
Please don't post (un)subscribe requests to the list. Use the www interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>.<br>
Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/<u></u>Support/Mailing_Lists</a><br>
</div></div></blockquote></div><br><br clear="all"><br>-- <br><span style="color:rgb(51,51,255)">-----------------------</span><br style="color:rgb(51,51,255)"><i><span style="color:rgb(51,51,255)">Regards,</span></i><br style="color:rgb(51,51,255)">
<span style="color:rgb(51,51,255);font-family:garamond,serif">Bipin Singh</span><br><br>