Hi Justin,<br>I&#39;ve tried with no restraint and the minimization stopped at the same time. I actually suspect that the restraint is not working at all. My top file is as below:-<br>#include &quot;ffG53a6_lipid.itp&quot;<br>
#include &quot;popc.itp&quot;<br>#ifdef POSRES<br>#include &quot;posres_popc.itp&quot;<br>#endif<br>#ifdef FLEX_SPC<br>#include &quot;flexspc.itp&quot;<br>#else<br>#include &quot;spc.itp&quot;<br>#endif<br><br>[ system ]<br>
; name<br>Pure POPC bilayer with 128 lipids 2460 water<br><br>[ molecules ]<br>; name  number<br>POPC    128<br>SOL    2460<br><br>my posres_popc.itp is as below:-<br>[ position_restraints ]<br>; atom  type      fx      fy      fz<br>
     1     1  1000  1000  1000<br>.....<br>    51     1  1000  1000  1000<br>    52     1  1000  1000  1000<br><br>my mdp file as below :-<br>define = -DPOSRES<br>integrator = steep<br>emtol = 1000<br>emstep = 0.01<br>nsteps = 100<br>
rlist = 1.2<br>coulombtype = PME<br>rcoulomb = 1.2<br>rvdw = 1.2<br>constraints = none<br><br>May i know that is the top, posres and also the mdp file correctly done?<br>And by looking at the log file in the previous mail, is it confirm that it manage to read the restraint file? <br>
The starting structure is taken from Peter Tielemen website and what i did is editconf to make box and i didn&#39;t solvate because the starting structure come together with water. I&#39;m not sure this is the correct way of doing it. Please advice.<br>
<br><br>Thanks in advance.<br>Bing<br><br><br><div class="gmail_quote">On Tue, Jun 23, 2009 at 7:16 PM, Justin A. Lemkul <span dir="ltr">&lt;<a href="mailto:jalemkul@vt.edu">jalemkul@vt.edu</a>&gt;</span> 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><div></div><div class="h5"><br>
<br>
Bing Bing 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! I&#39;m a postgraduate student from Malaysia and currently facing energy minimization problem for my membrane system.<br>
I&#39;m trying to perform energy minimization for an empty POPC membrane with position restraint on the lipid. Minimization stop due to inf Fmax. And I noticed that in the log file as below:-<br>
  Initiating Steepest Descents<br>
Max number of connections per atom is 25<br>
Total number of connections is 84976<br>
Max number of graph edges per atom is 4<br>
Total number of graph edges is 22896<br>
Going to use C-settle (2460 waters)<br>
wo = 0.333333, wh =0.333333, wohh = 3, rc = 0.08165, ra = 0.0384897<br>
rb = 0.0192448, rc2 = 0.1633, rone = 1, dHH = 0.1633, dOH = 0.1<br>
Started Steepest Descents on node 0 Tue Jun 23 17:03:43 2009<br>
Steepest Descents:<br>
   Tolerance (Fmax)   =  1.00000e+03<br>
   Number of steps    =          100<br>
Grid: 11 x 11 x 11 cells<br>
           Step           Time         Lambda<br>
              0        0.00000        0.00000<br>
   Energies (kJ/mol)<br>
           Bond          Angle    Proper Dih. Ryckaert-Bell.  Improper Dih.<br>
    1.64082e+05    1.71875e+04    4.06560e+03    6.45104e+03    8.36392e+03<br>
          LJ-14     Coulomb-14        LJ (SR)   Coulomb (SR)   Coul. recip.<br>
   2.12067e-314   2.12067e-314   -1.58963e+04   -1.05392e+05   -1.15292e+05<br>
 Position Rest.      Potential Pressure (bar)<br>
    0.00000e+00   -3.64297e+04            nan<br>
<br>
</blockquote>
<br></div></div>
>From the looks of some of these energies (i.e., LJ-14, Coulomb-14), it looks like something is severely wrong with your starting structure, and applying position restraints probably isn&#39;t helping.<br>
<br>
Try the minimization without restraints and see if it works.<br>
<br>
Also, have you set the appropriate &quot;define&quot; line in the .mdp file, corresponding to what is prescribed in the .top?<br>
<br>
-Justin<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="im">
-------------------------------------------------------<br>
Program mdrun_d, VERSION 4.0.5<br>
Source code file: nsgrid.c, line: 357<br>
Range checking error:<br>
Explanation: During neighborsearching, we assign each particle to a grid<br>
based on its coordinates. If your system contains collisions or parameter<br>
errors that give particles very high velocities you might end up with some<br>
coordinates being +-Infinity or NaN (not-a-number). Obviously, we cannot<br>
put these on a grid, so this is usually where we detect those errors.<br>
Make sure your system is properly energy-minimized and that the potential<br>
energy seems reasonable before trying again.<br>
Variable ci has value -2147483648. It should have been within [ 0 .. 1331 ]<br>
  the position restraint energy is 0 as highlighted. Here&#39;s the position restraint file :-<br>
 [ position_restraints ]<br>
; atom  type      fx      fy      fz<br>
     1     1  1000  1000  5000<br>
     5     1  1000  1000  1000<br>
     6     1  1000  1000  1000<br>
     7     1  1000  1000  1000<br>
     8     1  1000  1000  1000<br>
     9     1  1000  1000  1000<br>
    10     1  1000  1000  1000<br>
......<br>
<br>
whereby the restraint energy is 1000kcal/mol in x, y and z plane. In the log file the position restraint energy is stated as 0.0000e+00 (as highlighted above).<br>
I think that the position restraint is not working fine here. Is there anyways to solve this? did I miss out something here?<br>
I would be grateful if anyone can advice me on this.<br>
   <br>
 <br>
<br></div>
------------------------------------------------------------------------<br>
<br>
_______________________________________________<br>
gmx-users mailing list    <a href="mailto:gmx-users@gromacs.org" target="_blank">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&#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>
</blockquote>
<br>
-- <br>
========================================<br>
<br>
Justin A. Lemkul<br>
Ph.D. Candidate<br>
ICTAS Doctoral Scholar<br>
Department of Biochemistry<br>
Virginia Tech<br>
Blacksburg, VA<br>
jalemkul[at]<a href="http://vt.edu" target="_blank">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    <a href="mailto:gmx-users@gromacs.org" target="_blank">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&#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>
</blockquote></div><br>