My understanding is that the v-site algorithm is used for the virtual sites and LINCS is used for bonds not involving v-sites (and also angles if you choose contraint=angles). <br><br><div class="gmail_quote">On Sat, Jun 18, 2011 at 4:40 PM, <a href="mailto:chris.neale@utoronto.ca">chris.neale@utoronto.ca</a> <span dir="ltr"><<a href="mailto:chris.neale@utoronto.ca">chris.neale@utoronto.ca</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">Thank you Roland.<br>
<br>
I did use:<br>
<br>
constraints = all-bonds<br>
lincs-iter = 1<br>
lincs-order = 6<br>
constraint_algorithm = lincs<br>
<br>
From looking at the manual, I figured that angle and bond constraints<br>
would all be done by LINCS if I had done (A):<br>
<br>
pdb2gmx -vsite none<br>
constraints = h-angles<br>
(a combination that I have never tried)<br>
<br>
But when I use (B):<br>
<br>
pdb2gmx -vsite hydrogen<br>
constraints = all-bonds<br>
<br>
It seems possible to me that LINCS is not used but instead the<br>
position of the atom is simply built from a mathematical function.<br>
Perhaps this all stems from my lack of thorough understanding of<br>
LINCS, but it seems to me that there need be no iteration to simply<br>
place an atoms based on virtual_sites3 (which are constructed by<br>
pdb2gms -hydrogen)<br>
<br>
For now, I'll simply add a line to state that I built virtual sites<br>
for hydrogen atoms to make it clear, but I'd still like to understand<br>
the difference between options A and B, above, if you have some time.<br>
<br>
Thank you again,<br>
Chris.<br>
<br>
<br>
<br>
<br>
<br>
On Sat, Jun 18, 2011 at 4:13 PM, chris.neale at <a href="http://utoronto.ca" target="_blank">utoronto.ca</a> <<br>
chris.neale at <a href="http://utoronto.ca" target="_blank">utoronto.ca</a>> wrote:<br>
<br>
> Dear Users:<br>
><br>
> If I create the topology of a peptide like this:<br>
><br>
> pdb2gmx -f protein.gro -vsite hydrogens<br>
><br>
> And then simulate it in vacuum, is lincs used at all? I believe that<br>
> it is, as if I use a timestep that is too large then I get LINCS<br>
> warnings about angles rotating more than 30 degrees, but that warning<br>
> message could possibly have been written with the assumption that I<br>
> used LINCS and not virtual hydrogens.<br>
><br>
Probably. To make sure check the constraint-algorithm selected in your mdp.<br>
BTW: If you want to use large timecheck you should normally use<br>
constraints=all-bonds and lincs-order=6.<br>
<br>
<br>
><br>
> Finally, is there a method that needs to be named or cited in relation<br>
> to the fact that the angles are now constrained? Is that also done<br>
> with P-LINCS?<br>
><br>
This is also done with P-LINCS. Not sure whether one should sign something<br>
regarding the construction/usage of v-sites.<br>
<br>
Roland<br>
<br>
<br>
><br>
> Thank you,<br>
> Chris.<br>
<font color="#888888"><br>
<br>
--<br>
gmx-users mailing list <a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>
<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/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/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">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/Support/Mailing_Lists</a><br>
<br>
<br>
</font></blockquote></div><br><br clear="all"><br>-- <br>ORNL/UT Center for Molecular Biophysics <a href="http://cmb.ornl.gov">cmb.ornl.gov</a><br>865-241-1537, ORNL PO BOX 2008 MS6309<br>