<div dir="ltr">Hello,<br><br>The simulation does not run even for one step and therefore, it is not possible for me to check any of the components such as energy, temperature etc. There is no error message except for &quot;segmentation fault&quot;. It does the same thing if I try to energy minimize the system instead of running an MD simulation. The corresponding mdp file is pasted below:<br>
cpp&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; /lib/cpp<br>constraints&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; none<br>integrator&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; steep<br>nstcgsteep&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp; =&nbsp; 500<br>emtol&nbsp;&nbsp;&nbsp; &nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp; =&nbsp; 500.0<br>emstep&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 0.001<br>vdwtype&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; user<br>
energygrps &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; POL<br>energygrp_table&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; POL POL<br><br>The commands I use are as follows:<br>grompp -f em.mdp -p 25mer_wca_vac.top -c 3.5ns.gro<br>mdrun -s topol.tpr -table table.xvg <br><br>Regards<br>Sapna<br>
<br><br><div class="gmail_quote">On Wed, Jul 16, 2008 at 10:08 AM, David van der Spoel &lt;<a href="mailto:spoel@xray.bmc.uu.se">spoel@xray.bmc.uu.se</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;">
<div class="Ih2E3d">sapna sarupria wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Hello,<br>
<br>
Thanks David for your response. Actually I have used the same tables for simulations of the polymer in water and have had no problem with them. Those simulations run for 4 ns without a problem. So the tables are correct and I am sure of that. I am not using the CVS version and so I give the energy and second derivative in the table. &nbsp;I have used user-defined tables in the past and so the setup is correct in terms of changing the mdp file and the top file. The problem seems to be more system specific (meaning vacuum) than user error specific. Can you suggest any other thing that may be the problem.<br>

<br>
</blockquote></div>
Have you checked energy components, temperature etc.?<br>
<br>
<br>
<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><div class="Ih2E3d">
Thanks a lot for your help.<br>
<br>
Regards<br>
Sapna<br>
<br>
<br></div><div class="Ih2E3d">
On Wed, Jul 16, 2008 at 9:23 AM, David van der Spoel &lt;<a href="mailto:spoel@xray.bmc.uu.se" target="_blank">spoel@xray.bmc.uu.se</a> &lt;mailto:<a href="mailto:spoel@xray.bmc.uu.se" target="_blank">spoel@xray.bmc.uu.se</a>&gt;&gt; wrote:<br>

<br>
 &nbsp; &nbsp;sapna sarupria wrote:<br>
<br>
<br>
<br>
 &nbsp; &nbsp; &nbsp; &nbsp;---------- Forwarded message ----------<br>
 &nbsp; &nbsp; &nbsp; &nbsp;From: *sapna sarupria* &lt;<a href="mailto:sapna.sarupria@gmail.com" target="_blank">sapna.sarupria@gmail.com</a><br>
 &nbsp; &nbsp; &nbsp; &nbsp;&lt;mailto:<a href="mailto:sapna.sarupria@gmail.com" target="_blank">sapna.sarupria@gmail.com</a>&gt;<br></div><div class="Ih2E3d">
 &nbsp; &nbsp; &nbsp; &nbsp;&lt;mailto:<a href="mailto:sapna.sarupria@gmail.com" target="_blank">sapna.sarupria@gmail.com</a> &lt;mailto:<a href="mailto:sapna.sarupria@gmail.com" target="_blank">sapna.sarupria@gmail.com</a>&gt;&gt;&gt;<br>
 &nbsp; &nbsp; &nbsp; &nbsp;Date: Tue, Jul 15, 2008 at 9:16 AM<br>
 &nbsp; &nbsp; &nbsp; &nbsp;Subject: Using user-tables for simulations in vacuum<br>
 &nbsp; &nbsp; &nbsp; &nbsp;To: Discussion list for GROMACS users &lt;<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br></div><div><div></div><div class="Wj3C7c">
 &nbsp; &nbsp; &nbsp; &nbsp;&lt;mailto:<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>&gt; &lt;mailto:<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
 &nbsp; &nbsp; &nbsp; &nbsp;&lt;mailto:<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>&gt;&gt;&gt;<br>
<br>
<br>
 &nbsp; &nbsp; &nbsp; &nbsp;Hello users,<br>
<br>
 &nbsp; &nbsp; &nbsp; &nbsp;I am trying to do a simulation of a polymer chain (which is<br>
 &nbsp; &nbsp; &nbsp; &nbsp;simply a bead of unified methane molecules) in vacuum using<br>
 &nbsp; &nbsp; &nbsp; &nbsp;user-defined tables. The interactions are<br>
 &nbsp; &nbsp; &nbsp; &nbsp;Weeks-Chandler-Andersen instead of Lennard-Jones. However, &nbsp;when<br>
 &nbsp; &nbsp; &nbsp; &nbsp;I run the simulations, mdrun gives me a segmentation fault. When<br>
 &nbsp; &nbsp; &nbsp; &nbsp;I run simulations for the exact same configurations without the<br>
 &nbsp; &nbsp; &nbsp; &nbsp;tables (and therefore using VdW) the simulations run fine. In<br>
 &nbsp; &nbsp; &nbsp; &nbsp;addition, the starting configurations were obtained after 3ns<br>
 &nbsp; &nbsp; &nbsp; &nbsp;simulations of the polymer in water.<br>
<br>
 &nbsp; &nbsp; &nbsp; &nbsp;So I was wondering if there is any issue with using user-defined<br>
 &nbsp; &nbsp; &nbsp; &nbsp;tables with simulations in vacuum. If not do you have any idea<br>
 &nbsp; &nbsp; &nbsp; &nbsp;what could be going wrong.<br>
 &nbsp; &nbsp; &nbsp; &nbsp;Some details of the mdp file:<br>
 &nbsp; &nbsp; &nbsp; &nbsp;1. pbc is turned off.<br>
 &nbsp; &nbsp; &nbsp; &nbsp;2. there is no pressure coupling.<br>
 &nbsp; &nbsp; &nbsp; &nbsp;3. center of mass removal is set to angular.<br>
 &nbsp; &nbsp; &nbsp; &nbsp;4. temperature is 298 K and berendsen thermostat is used.<br>
 &nbsp; &nbsp; &nbsp; &nbsp;5. energy groups and table are defined (correctly).<br>
 &nbsp; &nbsp; &nbsp; &nbsp;6. the cut-offs are set to 1.0nm (box size is larger than 4 nm).<br>
 &nbsp; &nbsp; &nbsp; &nbsp;7. no constraints are being used.<br>
<br>
<br>
 &nbsp; &nbsp;Most obvious error source is the table itself. Do you use the CVS<br>
 &nbsp; &nbsp;version? In that case you should provide energy and force, otherwise<br>
 &nbsp; &nbsp;energy and second derivative. Distance units are in nm. Compare to<br>
 &nbsp; &nbsp;existing tables in share/gromacs/top<br>
<br>
<br>
 &nbsp; &nbsp; &nbsp; &nbsp;Thank you<br>
<br>
 &nbsp; &nbsp; &nbsp; &nbsp;Regards<br>
 &nbsp; &nbsp; &nbsp; &nbsp;Sapna<br>
<br>
<br>
<br>
 &nbsp; &nbsp; &nbsp; &nbsp;--  &nbsp; &nbsp; &nbsp; &nbsp;Sapna Sarupria<br>
 &nbsp; &nbsp; &nbsp; &nbsp;Ph.D. Student - Chemical Engineering<br>
 &nbsp; &nbsp; &nbsp; &nbsp;Rensselaer Polytechnic Institute<br>
 &nbsp; &nbsp; &nbsp; &nbsp;Troy, New York 12180<br>
 &nbsp; &nbsp; &nbsp; &nbsp;U.S.A.<br>
 &nbsp; &nbsp; &nbsp; &nbsp;Ph#: (518)276-3031<br>
 &nbsp; &nbsp; &nbsp; &nbsp;Life isn&#39;t about finding yourself. Life is about creating yourself.<br>
 &nbsp; &nbsp; &nbsp; &nbsp;George Bernard Shaw.<br>
 &nbsp; &nbsp; &nbsp; &nbsp;Dare to Dream<br>
<br>
<br>
 &nbsp; &nbsp; &nbsp; &nbsp;------------------------------------------------------------------------<br>
<br>
 &nbsp; &nbsp; &nbsp; &nbsp;_______________________________________________<br>
 &nbsp; &nbsp; &nbsp; &nbsp;gmx-users mailing list &nbsp; &nbsp;<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
 &nbsp; &nbsp; &nbsp; &nbsp;&lt;mailto:<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>&gt;<br></div></div><div class="Ih2E3d">
 &nbsp; &nbsp; &nbsp; &nbsp;<a href="http://www.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://www.gromacs.org/mailman/listinfo/gmx-users</a><br>
 &nbsp; &nbsp; &nbsp; &nbsp;Please search the archive at <a href="http://www.gromacs.org/search" target="_blank">http://www.gromacs.org/search</a><br>
 &nbsp; &nbsp; &nbsp; &nbsp;before posting!<br>
 &nbsp; &nbsp; &nbsp; &nbsp;Please don&#39;t post (un)subscribe requests to the list. Use the<br>
 &nbsp; &nbsp; &nbsp; &nbsp;www interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a><br></div>
 &nbsp; &nbsp; &nbsp; &nbsp;&lt;mailto:<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>&gt;.<div class="Ih2E3d"><br>
 &nbsp; &nbsp; &nbsp; &nbsp;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>
 &nbsp; &nbsp;--  &nbsp; &nbsp;David van der Spoel, Ph.D., Professor of Biology<br>
 &nbsp; &nbsp;Molec. Biophys. group, Dept. of Cell &amp; Molec. Biol., Uppsala University.<br>
 &nbsp; &nbsp;Box 596, 75124 Uppsala, Sweden. Phone: &nbsp;+46184714205. Fax: +4618511755.<br></div>
 &nbsp; &nbsp;<a href="mailto:spoel@xray.bmc.uu.se" target="_blank">spoel@xray.bmc.uu.se</a> &lt;mailto:<a href="mailto:spoel@xray.bmc.uu.se" target="_blank">spoel@xray.bmc.uu.se</a>&gt; &nbsp; &nbsp; &nbsp; <a href="mailto:spoel@gromacs.org" target="_blank">spoel@gromacs.org</a> &lt;mailto:<a href="mailto:spoel@gromacs.org" target="_blank">spoel@gromacs.org</a>&gt; &nbsp; <a href="http://folding.bmc.uu.se" target="_blank">http://folding.bmc.uu.se</a><div class="Ih2E3d">
<br>
 &nbsp; &nbsp;_______________________________________________<br>
 &nbsp; &nbsp;gmx-users mailing list &nbsp; &nbsp;<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br></div><div class="Ih2E3d">
 &nbsp; &nbsp;&lt;mailto:<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>&gt;<br></div><div class="Ih2E3d">
 &nbsp; &nbsp;<a href="http://www.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://www.gromacs.org/mailman/listinfo/gmx-users</a><br>
 &nbsp; &nbsp;Please search the archive at <a href="http://www.gromacs.org/search" target="_blank">http://www.gromacs.org/search</a> before<br>
 &nbsp; &nbsp;posting!<br>
 &nbsp; &nbsp;Please don&#39;t post (un)subscribe requests to the list. Use the www<br>
 &nbsp; &nbsp;interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a><br></div>
 &nbsp; &nbsp;&lt;mailto:<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>&gt;.<div class="Ih2E3d"><br>
 &nbsp; &nbsp;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>
-- <br>
Sapna Sarupria<br>
Ph.D. Student - Chemical Engineering<br>
Rensselaer Polytechnic Institute<br>
Troy, New York 12180<br>
U.S.A.<br>
Ph#: (518)276-3031<br>
Life isn&#39;t about finding yourself. Life is about creating yourself.<br>
George Bernard Shaw.<br>
Dare to Dream<br>
<br>
<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 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>
</div></blockquote><div><div></div><div class="Wj3C7c">
<br>
<br>
-- <br>
David van der Spoel, Ph.D., Professor of Biology<br>
Molec. Biophys. group, Dept. of Cell &amp; Molec. Biol., Uppsala University.<br>
Box 596, 75124 Uppsala, Sweden. Phone: &nbsp;+46184714205. Fax: +4618511755.<br>
<a href="mailto:spoel@xray.bmc.uu.se" target="_blank">spoel@xray.bmc.uu.se</a> &nbsp; &nbsp;<a href="mailto:spoel@gromacs.org" target="_blank">spoel@gromacs.org</a> &nbsp; <a href="http://folding.bmc.uu.se" target="_blank">http://folding.bmc.uu.se</a><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 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>
</div></div></blockquote></div><br><br clear="all"><br>-- <br>Sapna Sarupria<br>Ph.D. Student - Chemical Engineering<br>Rensselaer Polytechnic Institute<br>Troy, New York 12180<br>U.S.A.<br>Ph#: (518)276-3031<br>Life isn&#39;t about finding yourself. Life is about creating yourself. <br>
George Bernard Shaw. <br>Dare to Dream
</div>