<DIV>&nbsp;</DIV><DIV>&nbsp;</DIV><SPAN>On 02/11/11, <B class=name>James Starlight </B>&lt;jmsstarlight@gmail.com&gt; wrote:</SPAN>
<BLOCKQUOTE style="BORDER-LEFT: #00f 1px solid; PADDING-LEFT: 13px; MARGIN-LEFT: 0px" class=iwcQuote cite=mid:CAALQopxsVniKq45UsAkbMA_AEi-3334FkenSn9NpVrJc2=C7kA@mail.gmail.com type="cite">Justin, Mark,<br /><br />Ok<br />not I've included posre only to the lipid.itp wich consist of parameters for just 1 lipid<br /><br />#ifdef POSRES_LIPID<br />#include "posre_lipid.itp"<br />#endif<br /><br />to the minim.mdp I've included<br />define&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; = -DPOSRES_LIPID&nbsp;&nbsp;&nbsp; ; position restrain the lipids<br /><br />As the result I have not obtained any error in grompp but after minimization my lipid bilayer was perturbed. I've used very strong posres. Why it have not worked ?</BLOCKQUOTE>
<DIV>&nbsp;</DIV><DIV>We&nbsp;can't know,&nbsp;because we don't know what's in any files, nor why you think the bilayer was perturbed. If you would like&nbsp;useful help, you need to be exhaustive in your description of what you have done, what you have observed and what you think about those. It is easier for people who might volunteer free help to do something else with their time than to ask follow-up questions when few details were given originally.</DIV><DIV>&nbsp;</DIV><DIV>Mark&nbsp;</DIV><DIV>&nbsp;</DIV><BLOCKQUOTE style="BORDER-LEFT: #00f 1px solid; PADDING-LEFT: 13px; MARGIN-LEFT: 0px" class=iwcQuote cite=mid:CAALQopxsVniKq45UsAkbMA_AEi-3334FkenSn9NpVrJc2=C7kA@mail.gmail.com type="cite">
<DIV class="mimepart text html"><br /><br />
<DIV class=gmail_quote>2011/10/31 Justin A. Lemkul <SPAN dir=ltr>&lt;<a href="mailto:jalemkul@vt.edu" >jalemkul@vt.edu</A>&gt;</SPAN><br />
<BLOCKQUOTE style="BORDER-LEFT: #ccc 1px solid; MARGIN: 0px 0px 0px 0.8ex; PADDING-LEFT: 1ex" class=gmail_quote>
<DIV>
<DIV></DIV><DIV class=h5><br /><br />James Starlight wrote:<br />
<BLOCKQUOTE style="BORDER-LEFT: #ccc 1px solid; MARGIN: 0px 0px 0px 0.8ex; PADDING-LEFT: 1ex" class=gmail_quote>I've done all of that<br />i have "posre_lipid.itp" - for 1 lipid ( posres for each of 50 atoms )<br /><br /><br />I've included this in the topology of the bilayer<br /><br />; Include chain topologies<br />#include "gromos53a6_lipid.ff/<U></U>forcefield.itp"<br />#include "dppc.itp"<br /><br />; Include water topology<br />#include "gromos53a6_lipid.ff/spc.itp"<br /><br />; Include ion topologies<br />#include "gromos53a6_lipid.ff/ions.itp"<br /><br />#ifdef POSRES_LIPID<br />#include "posre_lipid.itp"<br />#endif<br /><br />; System specifications<br />[ system ]<br />128-Lipid DMPC Bilayer in water<br /><br />[ molecules ]<br />; molecule name nr.<br />DPPC 64<br />SOL &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;1193<br />SOL &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 690<br /><br /><br />Also I've tried to make posre for whole system ( large posres ) but it also was finished with same error<br /><br />Fatal error:<br />[ file posre_lipid.itp, line 6 ]:<br />Atom index (2) in position_restraints out of bounds (1-1).<br />This probably means that you have inserted topology section "position_restraints"<br />in a part belonging to a different molecule than you intended to.<br />In that case move the "position_restraints" section to the right molecule.<br />For more information and tips for troubleshooting, please check the GROMACS<br />website at <a href="http://www.gromacs.org/Documentation/Errors" target=1 >http://www.gromacs.org/<U></U>Documentation/Errors</A><br /><br />What's most true for such multy system. Generate posre only for 1 part or for whole system?<br /><br /></BLOCKQUOTE><br /></DIV></DIV>Position restraints can only be applied on a [moleculetype] basis. &nbsp;Thus, order matters a lot, and the whole system cannot be restrained in one file. &nbsp;See the example here:<br /><br /><a href="http://www.gromacs.org/Documentation/Errors#Atom_index_n_in_position_restraints_out_of_bounds" target=1 >http://www.gromacs.org/<U></U>Documentation/Errors#Atom_<U></U>index_n_in_position_<U></U>restraints_out_of_bounds</A><br /><br />As it stands now, the topology you've shown tries to apply lipid position restraints after the ions have been #included, which makes no sense at all.
<DIV>
<DIV></DIV><DIV class=h5><br /><br />-Justin<br /><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=1 >vt.edu</A> | (540) 231-9080<br /><a href="http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin" target=1 >http://www.bevanlab.biochem.<U></U>vt.edu/Pages/Personal/justin</A><br /><br />==============================<U></U>==========<br />-- <br />gmx-users mailing list &nbsp; &nbsp;<a href="mailto:gmx-users@gromacs.org" target=1 >gmx-users@gromacs.org</A><br /><a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target=1 >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=1 >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=1 >gmx-users-request@gromacs.org</A>.<br />Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target=1 >http://www.gromacs.org/<U></U>Support/Mailing_Lists</A><br /></DIV></DIV></BLOCKQUOTE></DIV><br /></DIV></BLOCKQUOTE>
<DIV>&nbsp;</DIV>