<div>Thanks Andreas.</div>
<div> </div>
<div>Is there any difference force-wise between this and our little workaround?</div>
<div> </div>
<div>Alex</div>
<div><br> </div>
<div class="gmail_quote">On Mon, May 3, 2010 at 2:17 PM, Kukol, Andreas <span dir="ltr"><<a href="mailto:a.kukol@herts.ac.uk">a.kukol@herts.ac.uk</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">Hi,<br><br>The POPC_53a6 topology works only with Gromacs versions lower than 4. Please find attached the topology compatible with Gromacs4.<br>
<br>I shall update files on the Gromacs website. Sorry for causing confusion.<br><br>Best wishes<br>Andreas<br><br>________________________________________<br>From: <a href="mailto:gmx-users-bounces@gromacs.org">gmx-users-bounces@gromacs.org</a> [<a href="mailto:gmx-users-bounces@gromacs.org">gmx-users-bounces@gromacs.org</a>] On Behalf Of Alex Smolyanitsky [<a href="mailto:asmolya1@asu.edu">asmolya1@asu.edu</a>]<br>
Sent: 03 May 2010 18:22<br>To: Discussion list for GROMACS users<br>Subject: Re: [gmx-users] POPC 53a6 topology (dihedral multiplicity)<br>
<div class="im"><br>Yup, those are indeed constants associated with phi2 (Fig. 1(B)), thanks.<br><br>Just in case anyone else is faced with the same issue, as far as forces are concerned, these dihedrals can be commented out. When considering the energy associated with dihedrals, just shift it appropriately.<br>
</div>
<div class="im">On Mon, May 3, 2010 at 10:57 AM, XAvier Periole <<a href="mailto:x.periole@rug.nl">x.periole@rug.nl</a><mailto:<a href="mailto:x.periole@rug.nl">x.periole@rug.nl</a>>> wrote:<br><br>HAve look at the paper describing the topology ...<br>
<br><br>On May 3, 2010, at 6:52 PM, Alex Smolyanitsky wrote:<br><br>Hello everyone,<br><br>I am trying to include a POPC 53a6 topology from LipidsForGro96_53a6.zip ( <a href="http://www.gromacs.org/index.php?title=Download_%26_Installation/User_contributions/Molecule_topologies" target="_blank">http://www.gromacs.org/index.php?title=Download_%26_Installation/User_contributions/Molecule_topologies</a> ) for a simple lipid in water sim (gromacs version 4.0.5).<br>
<br>I am using the lipid coordinates from the same archive and here's my topology:<br>***<br>#include "ffG53a6.itp"<br>#include "popc_53a6.itp"<br>#include "spc.itp"<br><br>[ system ]<br>
;Name<br>test<br><br>[ molecules ]<br>POPC 128<br>SOL 3585<br>***<br><br>grompp throws a fatal error: " Value of 'multiplicity' in Proper Dih. is 0, which is smaller than the minimum of 1" without any prior warnings.<br>
The error has to do with two dihedral angles with zero multiplicity:<br><br>...<br>22 23 24 25 1 90 -5.685 0<br>...<br>24 25 26 27 1 90 -5.685 0<br>...<br><br>So I looked at an ancient POPC topology for GROMOS87 and these dihedrals have zero equilibrium angle and a multiplicity of 3:<br>
<br>...<br>22 23 24 25 1 0.0 5.858 3<br>...<br>24 25 26 27 1 0.0 5.858 3<br>...<br><br>Did the author of the 53a6 topology want to simply remove these dihedrals by setting a multiplicity of zero? Or, should these two be set according to the older topology? I know similar questions have been asked here in the past, but I really want to make sure my topology is sane. Thanks a lot.<br>
<br>Alex<br>--<br></div>gmx-users mailing list <a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><mailto:<a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a>><br>
<div class="im"><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/search" target="_blank">http://www.gromacs.org/search</a> before posting!<br>
Please don't post (un)subscribe requests to the list. Use the<br></div>www interface or send it to <a href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a><mailto:<a href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>>.<br>
<div class="im">Can't post? Read <a href="http://www.gromacs.org/mailing_lists/users.php" target="_blank">http://www.gromacs.org/mailing_lists/users.php</a><br><br><br>--<br></div>gmx-users mailing list <a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><mailto:<a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a>><br>
<div class="im"><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/search" target="_blank">http://www.gromacs.org/search</a> before posting!<br>
Please don't post (un)subscribe requests to the list. Use the<br></div>www interface or send it to <a href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a><mailto:<a href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>>.<br>
<div>
<div></div>
<div class="h5">Can't post? Read <a href="http://www.gromacs.org/mailing_lists/users.php" target="_blank">http://www.gromacs.org/mailing_lists/users.php</a><br><br></div></div><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/search" target="_blank">http://www.gromacs.org/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/mailing_lists/users.php" target="_blank">http://www.gromacs.org/mailing_lists/users.php</a><br>
</blockquote></div><br>