Hi Carla,<br><br>I&#39;ve seen this behavior too but I didn&#39;t look at the code. I figured with the hydrogen included it checks the angle too or if the hydrogen is between the two heavy atoms. So the first pass just looks at distances between donor and acceptor atoms and doesn&#39;t consider that the hydrogen may be completely on the wrong side of the donor atom.<br>

<br>Did you look at any of the bonds that were left out in the 2nd pass?<br><br>Oliver<br><br><div class="gmail_quote">On 9 November 2010 14:08, Carla Jamous <span dir="ltr">&lt;<a href="mailto:carlajamous@gmail.com">carlajamous@gmail.com</a>&gt;</span> wrote:<br>

<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">Hi everyone,<br><br>I ran g_hbond on a trajectory. When g_hbond asks for two groups in the index file, I give:<br>

group1: &quot;Protein&quot;<br>group2: &quot;ligand&quot;<br>The output .ndx file contains 44 Hbonds.<br><br>
In order to verify my result on each and every Hbond, I ran g_hbond on the same trajectory but this time in my index file, I gave atom triplets, so when g_hbond asks for two groups, I give for example:<br>group1: &quot;a4810_a4811_a1857&quot;<br>


group2: &quot;a4810_a4811_a1857&quot;<br>Surprisingly, 10 of the 44 Hbonds are not found using this method.<br><br>Am I doing something wrong or is there a problem with g_hbond?<br><br>Thank you,<br><font color="#888888"><br>

Carla<br>
</font><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&#39;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&#39;t post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br></blockquote></div><br>