no prob you wait for some time it will write topology file and the proceed further..<br><br><div class="gmail_quote">On Thu, Apr 24, 2008 at 5:04 PM,  &lt;<a href="mailto:gadies@fh.huji.ac.il">gadies@fh.huji.ac.il</a>&gt; wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Hi,<br>
<br>
The full mdrun output when the simulation blows up is:<br>
<br>
Step 0, time 0 (ps) &nbsp;LINCS WARNING<br>
relative constraint deviation after LINCS:<br>
max inf (between atoms 6 and 7) rms inf<br>
bonds that rotated more than 30 degrees:<br>
 &nbsp; &nbsp; &nbsp; &nbsp; .<br>
 &nbsp; &nbsp; &nbsp; &nbsp; .<br>
 &nbsp; &nbsp; &nbsp; &nbsp; .<br>
Warning: 1-4 interaction between 1 and 6 at distance 324490777485109120.000 which is larger than the 1-4 table size 1.000 nm<br>
These are ignored for the rest of the simulation<br>
This usually means your system is exploding,<br>
if not, you should increase table-extension in your mdp file<br>
<br>
Thanks Gadi<br>
<br>
Quoting &quot;Justin A. Lemkul&quot; &lt;<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@vt.edu</a>&gt;:<br>
<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<br>
Did your minimization converge to a nice, negative potential energy?  &nbsp;Also, what<br>
is the full mdrun output when the simulation blows up? &nbsp;You&#39;ve &nbsp;quoted the last<br>
bit, but it would be more informative to see the whole output, or at least a<br>
description of it (i.e., LINCS warnings).<br>
<br>
-Justin<br>
<br>
Quoting <a href="mailto:gadies@fh.huji.ac.il" target="_blank">gadies@fh.huji.ac.il</a>:<br>
<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<br>
Hi,<br>
<br>
I tried using no temperature coupling and no constraints, that did not help.<br>
If I run the simulation without the position restraints it works OK.<br>
Could there by something wrong with the posre file (I&#39;ll attach it at<br>
the end of the email).<br>
I want one of the atoms in the lipid to stay in one place and the rest<br>
to wiggle around it. Is there any other way to do that which is not<br>
position restraints.<br>
<br>
Thanks Gadi<br>
<br>
<br>
[ position_restraints ]<br>
; atom &nbsp;type &nbsp; &nbsp; &nbsp;fx &nbsp; &nbsp; &nbsp;fy &nbsp; &nbsp; &nbsp;fz<br>
 &nbsp; &nbsp; &nbsp;1 &nbsp; &nbsp; 1 &nbsp;0 &nbsp;0 &nbsp;0<br>
 &nbsp; &nbsp; &nbsp;2 &nbsp; &nbsp; 1 &nbsp;0 &nbsp;0 &nbsp;0<br>
 &nbsp; &nbsp; &nbsp;3 &nbsp; &nbsp; 1 &nbsp;0 &nbsp;0 &nbsp;0<br>
 &nbsp; &nbsp; &nbsp;4 &nbsp; &nbsp; 1 &nbsp;0 &nbsp;0 &nbsp;0<br>
 &nbsp; &nbsp; &nbsp;5 &nbsp; &nbsp; 1 &nbsp;0 &nbsp;0 &nbsp;0<br>
 &nbsp; &nbsp; &nbsp;6 &nbsp; &nbsp; 1 &nbsp;0 &nbsp;0 &nbsp;0<br>
 &nbsp; &nbsp; &nbsp;7 &nbsp; &nbsp; 1 &nbsp;0 &nbsp;0 &nbsp;0<br>
 &nbsp; &nbsp; &nbsp;8 &nbsp; &nbsp; 1 &nbsp;1000.0 &nbsp;1000.0 &nbsp;1000.0<br>
 &nbsp; &nbsp; &nbsp;9 &nbsp; &nbsp; 1 &nbsp;0 &nbsp;0 &nbsp;0<br>
 &nbsp; &nbsp; 10 &nbsp; &nbsp; 1 &nbsp;0 &nbsp;0 &nbsp;0<br>
 &nbsp; &nbsp; 11 &nbsp; &nbsp; 1 &nbsp;0 &nbsp;0 &nbsp;0<br>
 &nbsp; &nbsp; 12 &nbsp; &nbsp; 1 &nbsp;0 &nbsp;0 &nbsp;0<br>
 &nbsp; &nbsp; 13 &nbsp; &nbsp; 1 &nbsp;0 &nbsp;0 &nbsp;0<br>
 &nbsp; &nbsp; 14 &nbsp; &nbsp; 1 &nbsp;0 &nbsp;0 &nbsp;0<br>
 &nbsp; &nbsp; 15 &nbsp; &nbsp; 1 &nbsp;0 &nbsp;0 &nbsp;0<br>
 &nbsp; &nbsp; 16 &nbsp; &nbsp; 1 &nbsp;0 &nbsp;0 &nbsp;0<br>
 &nbsp; &nbsp; 17 &nbsp; &nbsp; 1 &nbsp;0 &nbsp;0 &nbsp;0<br>
 &nbsp; &nbsp; 18 &nbsp; &nbsp; 1 &nbsp;0 &nbsp;0 &nbsp;0<br>
 &nbsp; &nbsp; 19 &nbsp; &nbsp; 1 &nbsp;0 &nbsp;0 &nbsp;0<br>
 &nbsp; &nbsp; 20 &nbsp; &nbsp; 1 &nbsp;0 &nbsp;0 &nbsp;0<br>
 &nbsp; &nbsp; 21 &nbsp; &nbsp; 1 &nbsp;0 &nbsp;0 &nbsp;0<br>
 &nbsp; &nbsp; 22 &nbsp; &nbsp; 1 &nbsp;0 &nbsp;0 &nbsp;0<br>
 &nbsp; &nbsp; 23 &nbsp; &nbsp; 1 &nbsp;0 &nbsp;0 &nbsp;0<br>
 &nbsp; &nbsp; 24 &nbsp; &nbsp; 1 &nbsp;0 &nbsp;0 &nbsp;0<br>
 &nbsp; &nbsp; 25 &nbsp; &nbsp; 1 &nbsp;0 &nbsp;0 &nbsp;0<br>
 &nbsp; &nbsp; 26 &nbsp; &nbsp; 1 &nbsp;0 &nbsp;0 &nbsp;0<br>
 &nbsp; &nbsp; 27 &nbsp; &nbsp; 1 &nbsp;0 &nbsp;0 &nbsp;0<br>
 &nbsp; &nbsp; 28 &nbsp; &nbsp; 1 &nbsp;0 &nbsp;0 &nbsp;0<br>
 &nbsp; &nbsp; 29 &nbsp; &nbsp; 1 &nbsp;0 &nbsp;0 &nbsp;0<br>
 &nbsp; &nbsp; 30 &nbsp; &nbsp; 1 &nbsp;0 &nbsp;0 &nbsp;0<br>
 &nbsp; &nbsp; 31 &nbsp; &nbsp; 1 &nbsp;0 &nbsp;0 &nbsp;0<br>
 &nbsp; &nbsp; 32 &nbsp; &nbsp; 1 &nbsp;0 &nbsp;0 &nbsp;0<br>
 &nbsp; &nbsp; 33 &nbsp; &nbsp; 1 &nbsp;0 &nbsp;0 &nbsp;0<br>
 &nbsp; &nbsp; 34 &nbsp; &nbsp; 1 &nbsp;0 &nbsp;0 &nbsp;0<br>
 &nbsp; &nbsp; 35 &nbsp; &nbsp; 1 &nbsp;0 &nbsp;0 &nbsp;0<br>
 &nbsp; &nbsp; 36 &nbsp; &nbsp; 1 &nbsp;0 &nbsp;0 &nbsp;0<br>
 &nbsp; &nbsp; 37 &nbsp; &nbsp; 1 &nbsp;0 &nbsp;0 &nbsp;0<br>
 &nbsp; &nbsp; 38 &nbsp; &nbsp; 1 &nbsp;0 &nbsp;0 &nbsp;0<br>
 &nbsp; &nbsp; 39 &nbsp; &nbsp; 1 &nbsp;0 &nbsp;0 &nbsp;0<br>
 &nbsp; &nbsp; 40 &nbsp; &nbsp; 1 &nbsp;0 &nbsp;0 &nbsp;0<br>
 &nbsp; &nbsp; 41 &nbsp; &nbsp; 1 &nbsp;0 &nbsp;0 &nbsp;0<br>
 &nbsp; &nbsp; 42 &nbsp; &nbsp; 1 &nbsp;0 &nbsp;0 &nbsp;0<br>
 &nbsp; &nbsp; 43 &nbsp; &nbsp; 1 &nbsp;0 &nbsp;0 &nbsp;0<br>
 &nbsp; &nbsp; 44 &nbsp; &nbsp; 1 &nbsp;0 &nbsp;0 &nbsp;0<br>
 &nbsp; &nbsp; 45 &nbsp; &nbsp; 1 &nbsp;0 &nbsp;0 &nbsp;0<br>
 &nbsp; &nbsp; 46 &nbsp; &nbsp; 1 &nbsp;0 &nbsp;0 &nbsp;0<br>
 &nbsp; &nbsp; 47 &nbsp; &nbsp; 1 &nbsp;0 &nbsp;0 &nbsp;0<br>
 &nbsp; &nbsp; 48 &nbsp; &nbsp; 1 &nbsp;0 &nbsp;0 &nbsp;0<br>
 &nbsp; &nbsp; 49 &nbsp; &nbsp; 1 &nbsp;0 &nbsp;0 &nbsp;0<br>
 &nbsp; &nbsp; 50 &nbsp; &nbsp; 1 &nbsp;0 &nbsp;0 &nbsp;0<br>
<br>
Quoting Mark Abraham &lt;<a href="mailto:Mark.Abraham@anu.edu.au" target="_blank">Mark.Abraham@anu.edu.au</a>&gt;:<br>
<br>
&gt; <a href="mailto:gadies@fh.huji.ac.il" target="_blank">gadies@fh.huji.ac.il</a> wrote:<br>
&gt;<br>
&gt; Please generate replies to the mailing list sensibly. It&#39;s very hard to<br>
&gt; work out who you are quoting where in this email. It&#39;d be easy just to<br>
&gt; ignore it, and that&#39;s the last thing you should want.<br>
&gt;<br>
&gt;&gt; Hello,<br>
&gt;&gt;<br>
&gt;&gt; I&#39;m trying to run a simulation with one lipid molecule in a fixed<br>
&gt;&gt; simulation box without pbc.<br>
&gt;&gt; I set position restraints on one of the atoms and I seem to keep on<br>
&gt;&gt; getting the following error:<br>
&gt;<br>
&gt; Why do you want MD of a single lipid and a position restraint on one<br>
&gt; atom? If you just want to see a lipid wiggle around, then don&#39;t create<br>
&gt; more numerical complexity.<br>
&gt;<br>
&gt;&gt; The lipid is DPPC.The .itp file I got from a coworker that used it<br>
&gt;&gt; before for membrane MD simulations.<br>
&gt;&gt; To minimize the system I used steepest descent with the following<br>
&gt;&gt; parameters:<br>
&gt;<br>
&gt; Well if your grompp and mdrun completed successfully and without<br>
&gt; warnings then you&#39;ve probably got an OK topology and structure.<br>
&gt;<br>
&gt; The combination of all bond constraints, a single lipid in vacuo,<br>
&gt; temperature coupling and a single position restraint sounds like a<br>
&gt; recipe for a numerical disaster, i.e. blowing up. Try with fewer of<br>
&gt; those things applied, or get a more physically-reasonable system, like<br>
&gt; a whole membrane.<br>
&gt;<br>
&gt; Mark<br>
&gt; _______________________________________________<br>
&gt; gmx-users mailing list &nbsp; &nbsp;<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
&gt; <a href="http://www.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://www.gromacs.org/mailman/listinfo/gmx-users</a><br>
&gt; Please search the archive at <a href="http://www.gromacs.org/search" target="_blank">http://www.gromacs.org/search</a> before posting!<br>
&gt; Please don&#39;t post (un)subscribe requests to the list. Use the www<br>
&gt; interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>.<br>
&gt; Can&#39;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>
----------------------------------------------------------------<br>
This message was sent using IMP, the Internet Messaging Program.<br>
<br>
_______________________________________________<br>
gmx-users mailing list &nbsp; &nbsp;<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<a href="http://www.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://www.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&#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" target="_blank">gmx-users-request@gromacs.org</a>.<br>
Can&#39;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>
</blockquote>
<br>
<br>
<br>
========================================<br>
<br>
Justin A. Lemkul<br>
Graduate Research Assistant<br>
Department of Biochemistry<br>
Virginia Tech<br>
Blacksburg, VA<br>
<a href="mailto:jalemkul@vt.edu" target="_blank">jalemkul@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>
gmx-users mailing list &nbsp; &nbsp;<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<a href="http://www.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://www.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&#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" target="_blank">gmx-users-request@gromacs.org</a>.<br>
Can&#39;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>
</blockquote>
<br>
<br>
<br>
----------------------------------------------------------------<br>
This message was sent using IMP, the Internet Messaging Program.<br>
<br>
_______________________________________________<br>
gmx-users mailing list &nbsp; &nbsp;<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<a href="http://www.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://www.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&#39;t post (un)subscribe requests to the list. Use thewww interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>.<br>
Can&#39;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>