Hi<br><br><div class="gmail_quote">2012/5/17 Lara Bunte <span dir="ltr"><<a href="mailto:lara.bunte@yahoo.de" target="_blank">lara.bunte@yahoo.de</a>></span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi<br>
<br>
Just for correct understanding: In this "empty" [ dihedrals ] bock are 11 lines without equilibrium angle and force constant, that looks like:<br>
<br>
<br>
[ dihedrals ]<br>
; ai aj ak al funct c0 c1 c2 c3 c4 c5<br>
2 1 19 8 2<br>
6 8 9 10 2<br>
<br>
<br>
Did this mean, that this are really physical dihedrals, that were given by pdb2gmx and the only reason that they are "empty" is, that they are not parametrized? <br>
Is this correct?<br>
<div class="im"><br></div></blockquote><div><br></div><div>No, it means that defaukt values (i.e. listed in [ dihedraltypes ] section ) will be applied</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div class="im">
Greetings<br>
Lara<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
----- Ursprüngliche Message -----<br>
</div><div class="im">Von: Justin A. Lemkul <<a href="mailto:jalemkul@vt.edu">jalemkul@vt.edu</a>><br>
An: Lara Bunte <<a href="mailto:lara.bunte@yahoo.de">lara.bunte@yahoo.de</a>>; Discussion list for GROMACS users <<a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a>><br>
CC:<br>
</div>Gesendet: 17:02 Donnerstag, 17.Mai 2012<br>
<div class="im HOEnZb">Betreff: Re: [gmx-users] Two [ dihedrals ] sections in topology<br>
<br>
<br>
<br>
</div><div class="HOEnZb"><div class="h5">On 5/17/12 10:48 AM, Lara Bunte wrote:<br>
> Hi<br>
><br>
> grompp -f em.mdp -p topol.top -c solvated.gro -o em.tpr<br>
><br>
> gives me errors like<br>
><br>
> ERROR 1 [file topol.top, line 233]:<br>
> No default Improper Dih. types<br>
><br>
> and so on. This line 233 is the starting of my empty [ dihedrals ] block. I erase it by hand out of my topology, so that only the correct [ dihedrals ] out of the .rtp and ffbonded.itp file are in the topology. After this grompp works without error and even without any warning.<br>
><br>
><br>
> So the question how to force pdb2gmx to use only my [ dihedrals ] block is important for me, because I don't like to erase it always by hand out of the topol.top file.<br>
><br>
<br>
This goes back to something that has been stated several times now, and was discussed the other day. If pdb2gmx finds 4 atoms connected via bonds, it assigns a dihedral to it. That's just physically real. If your force field then doesn't have parameters for that particular interaction, you get the error shown above. That doesn't mean you should go deleting things until grompp "works" because what will result will likely be senseless.<br>
<br>
If your "empty" (which is not really empty, as we've been saying) dihedrals are causing the problems, then you need to add parameters for them in ffbonded.itp or in the topology itself, not hack the topology apart.<br>
<br>
-Justin<br>
<br>
-- ========================================<br>
<br>
Justin A. Lemkul, Ph.D.<br>
Research Scientist<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.vt.edu/Pages/Personal/justin</a><br>
<br>
========================================<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>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br>Cordiali saluti, Dr.Oteri Francesco<br>