<html>
  <head>
    <meta content="text/html; charset=ISO-8859-1"
      http-equiv="Content-Type">
  </head>
  <body bgcolor="#FFFFFF" text="#000000">
    On 28/10/2011 2:30 AM, James Starlight wrote:
    <blockquote
cite="mid:CAALQopwUKWABUOhe9tvELq4p7Ow37vzD-_jWkkty3WmMUaCzkg@mail.gmail.com"
      type="cite">That's energy ouptut from minimization with that
      parametries ( there is also 1 step of steep minimization before
      that )<br>
      <br>
      integrator&nbsp;&nbsp;&nbsp; = l-bfgs<br>
      emtol&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; = 0.001 &nbsp;&nbsp;&nbsp; <br>
      emstep&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0.001&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; Energy step size<br>
      nsteps&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; = 5000000000&nbsp;&nbsp;&nbsp; &nbsp; &nbsp;&nbsp;&nbsp; ; Maximum number of
      (minimization) steps to perform<br>
      <br>
      <br>
      ; Parameters describing how to find the neighbors of each atom and
      how to calculate the interactions<br>
      nstlist&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; = 1&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; ; Frequency to update the neighbor list
      and long range forces<br>
      ns_type&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; = grid&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; ; Method to determine neighbor list
      (simple, grid)<br>
      rlist&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; = 1.2&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; ; Cut-off for making neighbor list
      (short range forces)<br>
      coulombtype&nbsp;&nbsp;&nbsp; = Shift&nbsp;&nbsp;&nbsp; <br>
      rcoulomb&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1.0<br>
      rcoulomb_switch&nbsp;&nbsp;&nbsp; = 0.7<br>
      vdwtype&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = Shift<br>
      rvdw&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; = 1.0<br>
      rvdw_switch&nbsp;&nbsp;&nbsp;&nbsp; = 0.7<br>
      <br>
      <br>
      Output:<br>
      <br>
      Low-Memory BFGS Minimizer converged to Fmax &lt; 0.001 in 3646
      steps<br>
      Potential Energy&nbsp; =&nbsp; 2.28946300988746e+03<br>
      Maximum force&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 8.47298141911052e-04 on atom 100<br>
      Norm of force&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 3.90346504086500e-04<br>
      <br>
      I'm not sure about the succses of that minimization due to the big
      Epot. &nbsp;&nbsp; <br>
      <br>
      <fieldset class="mimeAttachmentHeader"></fieldset>
      <br>
    </blockquote>
    <br>
    &nbsp;You look like you are shuffling deck chairs on the Titanic. I
    suggested you actually go and visualize your input and output
    trajectories. Sadly that doesn't seem to have happened.<br>
    <br>
    Mark<br>
  </body>
</html>