<html><head><style type="text/css"><!-- DIV {margin:0px;} --></style></head><body><div style="font-family:times new roman, new york, times, serif;font-size:12pt">Hi,<br>I have got a segmentation fault error when I run mdrun for 1ns. I would be glad if you let me know what is going on and how can I solve it.<br>Thank you very much.<br><br>Best,<br><br><div>&nbsp;</div>Jahanshah Ashkani,<br>PhD student of Biotechnology &amp; Genetics,<br>University of the Western Cape,<br>Biotechnology Department,<br>Private Bag X17,<br>7735 Bellville,<br>Cape Town,<br>South Africa<br>jashkani@mail.biotech.uwc.ac.za<div style="font-family: times new roman,new york,times,serif; font-size: 12pt;"><br><br><div style="font-family: times new roman,new york,times,serif; font-size: 12pt;">----- Original Message ----<br>From: Justin A. Lemkul &lt;jalemkul@vt.edu&gt;<br>To: JMandumpal &lt;jesbman@rediffmail.com&gt;<br>Cc: gmx-users@gromacs.org<br>Sent: Tuesday, December 4, 2007
 4:07:36 AM<br>Subject: Re: Re: Re: [gmx-users] Energy minimisation<br><br>
Quoting JMandumpal &lt;<a ymailto="mailto:jesbman@rediffmail.com" href="mailto:jesbman@rediffmail.com">jesbman@rediffmail.com</a>&gt;:<br><br>&gt;<br>&gt; Dear Justin,<br>&gt;<br>&gt; I made some changes in my .mdp file, but the result is same.<br>&gt;<br>&gt; My box size, as mentioned in the previous mail, is 20.1 A. The cut
 off I used<br>&gt; is 8 angstrom.<br>&gt;<br>&gt; I attach the input file here.<br>&gt; **************************************<br>&gt;<br>&gt; cpp&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = /lib/cpp<br>&gt; define&nbsp; &nbsp; &nbsp; &nbsp;  = -DPOSRES-DFLEX_SPC<br>&gt; constraints&nbsp; &nbsp; =&nbsp; none<br>&gt; morse&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; =&nbsp; no<br>&gt; integrator&nbsp; &nbsp;  = steep<br>&gt; nsteps&nbsp; &nbsp; &nbsp; &nbsp;  = 2000<br>&gt; ; Energy minimizing<br>&gt; emtol&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = 100<br>&gt; emstep&nbsp; &nbsp; &nbsp; &nbsp;  = 0.00001<br>&gt; nstcomm&nbsp; &nbsp; &nbsp; &nbsp; = 1<br>&gt; ns_type&nbsp; &nbsp; &nbsp; &nbsp; = grid<br>&gt; rlist&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = 0.8<br>&gt; coulumbtype&nbsp; &nbsp; = pme<br>&gt; rcoulomb&nbsp; &nbsp; &nbsp;  = 0.8<br>&gt; vdw-type&nbsp; &nbsp; &nbsp;  = cut-off<br>&gt; rvdw&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  = 0.8<br>&gt; nstenergy&nbsp; &nbsp; &nbsp;
 = 10<br>&gt; Tcoupl&nbsp; &nbsp; &nbsp; &nbsp;  = no<br>&gt; Pcoupl&nbsp; &nbsp; &nbsp; &nbsp;  = no<br>&gt; gen_vel&nbsp; &nbsp; &nbsp; &nbsp; = no<br>&gt; *************************************************<br>&gt;<br>&gt; And the output message is :<br>&gt;<br>&gt;
 ---------------------------------------------------------------------------<br>&gt; creating statusfile for 1 node...<br>&gt;<br>&gt; Back Off! I just backed up mdout.mdp to ./#mdout.mdp.5#<br>&gt; checking input for internal consistency...<br>&gt; calling /lib/cpp...<br>&gt; processing topology...<br>&gt; Generated 332520 of the 332520 non-bonded parameter combinations<br>&gt; Generating 1-4 interactions: fudge = 0.5<br>&gt; Generated 332520 of the 332520 1-4 parameter combinations<br>&gt; Excluding 2 bonded neighbours for SOL 258<br>&gt; NOTE:<br>&gt;&nbsp;  System has non-zero total charge: -2.580000e-01<br>&gt;<br>&gt; processing coordinates...<br>&gt; double-checking input for internal consistency...<br>&gt; ERROR: The cut-off length is longer than half the shortest box vector
 or<br>&gt; longer than the smallest box diagonal element. Increase the box size
 or<br>&gt; decrease rlist.<br>&gt; -------------------------------------------------------<br>&gt; Program grompp, VERSION 3.3.2<br>&gt; Source code file: grompp.c, line: 1111<br>&gt;<br>&gt; Fatal error:<br>&gt; There were 1 error(s) processing your input<br>&gt; -------------------------------------------------------<br>&gt;<br>&gt; "Stop Drinking My Beer !" (The Amps)<br>&gt;
&nbsp; ----------------------------------------------------------------------<br>&gt;<br>&gt;<br>&gt; HOW CAN I SETTLE THIS ISSUE?<br><br>Did you read Tsjerk's message from this morning?&nbsp; His suggestions
 regarding PBC<br>and verifying your box dimensions are good ones that I had not
 considered.<br><br>Also note that you probably have a broken topology, as the net charge
 on your<br>system is non-integral.&nbsp; This has nothing to do with the problem at
 hand, but<br>it will likely be problematic later.<br><br>-Justin<br><br>&gt;<br>&gt;<br>&gt; REGARDS<br>&gt; JESTIN<br>&gt;<br><br><br><br>========================================<br><br>Justin A. Lemkul<br>Graduate Research Assistant<br>Department of Biochemistry<br>Virginia Tech<br>Blacksburg, VA<br><a ymailto="mailto:jalemkul@vt.edu" href="mailto:jalemkul@vt.edu">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 ymailto="mailto:gmx-users@gromacs.org" href="mailto:gmx-users@gromacs.org">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't post (un)subscribe requests to the list. Use the<br>www interface or send it to <a ymailto="mailto:gmx-users-request@gromacs.org" 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></div><br></div></div><br>
      <hr size=1>Be a better sports nut! Let your teams follow you 
with Yahoo Mobile. <a href="http://us.rd.yahoo.com/evt=51731/*http://mobile.yahoo.com/sports;_ylt=At9_qDKvtAbMuh1G1SQtBI7ntAcJ
">Try it now.</a></body></html>