<html>
  <head>
    <meta content="text/html; charset=ISO-8859-1"
      http-equiv="Content-Type">
  </head>
  <body bgcolor="#CCCCFF" text="#330000">
    Hey Yuehua,<br>
    <br>
    I think that means the program was not able to optimized the
    structure with the parameter you provided.<br>
    Try to reduce etot_conv_thr or to increase the maximum number of
    steps by setting nstep.<br>
    Hope it helps,<br>
    <br>
    &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &Eacute;ric.<br>
    <br>
    On 10/11/2011 10:53 PM, xu yuehua wrote:
    <blockquote
cite="mid:CABZKPKD9mLX_abWv-ySW+i61FnD28kzSL0Ro0OoyR5X0aC14Kg@mail.gmail.com"
      type="cite">
      <div>Hi, all</div>
      <div>&nbsp;</div>
      <div>I&nbsp;relaxed a system, my relaxation criterion is <br
          clear="all">
        &nbsp;etot_conv_thr=1.D-5<br>
        &nbsp;forc_conv_thr=1.D-4, </div>
      <div>and conv_thr&nbsp;&nbsp;&nbsp; = 1.D-8 , </div>
      <div>&nbsp;</div>
      <div>and in my output file,&nbsp; i see the force_conv_thr is met, and
        the conv_thr seems also met.</div>
      <div>for in the last relaxation step, the printed information is :</div>
      <div>&nbsp;</div>
      <div><br>
        <em>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp; -551.13047812 Ry<br>
          &nbsp;&nbsp;&nbsp;&nbsp; Harris-Foulkes estimate&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp; -551.13047807 Ry<br>
          &nbsp;&nbsp;&nbsp;&nbsp; estimated scf accuracy&nbsp;&nbsp;&nbsp; &lt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2.4E-11 Ry<br>
          &nbsp;<br>
          &nbsp;&nbsp;&nbsp;&nbsp; The total energy is the sum of the following terms:<br>
          &nbsp;<br>
          &nbsp;&nbsp;&nbsp;&nbsp; one-electron contribution =&nbsp;&nbsp; -2144.14675345 Ry<br>
          &nbsp;&nbsp;&nbsp;&nbsp; hartree contribution&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp; 1097.77738769 Ry<br>
          &nbsp;&nbsp;&nbsp;&nbsp; xc contribution&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp; -134.85772630 Ry<br>
          &nbsp;&nbsp;&nbsp;&nbsp; ewald contribution&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp; 630.09661393 Ry<br>
          &nbsp;<br>
          &nbsp;&nbsp;&nbsp;&nbsp; convergence has been achieved in&nbsp;&nbsp; 9 iterations<br>
          &nbsp;<br>
          &nbsp;&nbsp;&nbsp;&nbsp; Forces acting on atoms (Ry/au):<br>
          &nbsp;<br>
          &nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp;&nbsp; 1 type&nbsp; 2&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00006652&nbsp;&nbsp;&nbsp;
          0.00001074&nbsp;&nbsp;&nbsp; 0.00001717<br>
          &nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp;&nbsp; 2 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp; -0.00001149&nbsp;&nbsp;
          -0.00003650&nbsp;&nbsp;&nbsp; 0.00002296<br>
          &nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp;&nbsp; 3 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00002092&nbsp;&nbsp;&nbsp;
          0.00003527&nbsp;&nbsp;&nbsp; 0.00001095<br>
          &nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp;&nbsp; 4 type&nbsp; 2&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp; -0.00006106&nbsp;&nbsp;
          -0.00005036&nbsp;&nbsp;&nbsp; 0.00002893<br>
          &nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp;&nbsp; 5 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00012199&nbsp;&nbsp; -0.00005785&nbsp;&nbsp;
          -0.00011993<br>
          &nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp;&nbsp; 6 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp; -0.00005472&nbsp;&nbsp;&nbsp;
          0.00002358&nbsp;&nbsp;&nbsp; 0.00009707<br>
          &nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp;&nbsp; 7 type&nbsp; 2&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00003025&nbsp;&nbsp;&nbsp;
          0.00003532&nbsp;&nbsp;&nbsp; 0.00004727<br>
          &nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp;&nbsp; 8 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00001390&nbsp;&nbsp;&nbsp;
          0.00009541&nbsp;&nbsp;&nbsp; 0.00000283<br>
          &nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp;&nbsp; 9 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp; -0.00007213&nbsp;&nbsp; -0.00012498&nbsp;&nbsp;
          -0.00005854<br>
          &nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 10 type&nbsp; 2&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp; -0.00000989&nbsp;&nbsp; -0.00001883&nbsp;&nbsp;
          -0.00002476<br>
          &nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 11 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00002862&nbsp;&nbsp;
          -0.00010630&nbsp;&nbsp;&nbsp; 0.00001950<br>
          &nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 12 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00000813&nbsp;&nbsp;&nbsp; 0.00006367&nbsp;&nbsp;
          -0.00002737<br>
          &nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 13 type&nbsp; 2&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp; -0.00010953&nbsp;&nbsp;&nbsp;
          0.00005997&nbsp;&nbsp;&nbsp; 0.00012796<br>
          &nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 14 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00012468&nbsp;&nbsp;&nbsp;
          0.00004819&nbsp;&nbsp;&nbsp; 0.00006199<br>
          &nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 15 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp; -0.00003227&nbsp;&nbsp; -0.00002568&nbsp;&nbsp;
          -0.00012529<br>
          &nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 16 type&nbsp; 2&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00000131&nbsp;&nbsp;&nbsp;
          0.00002155&nbsp;&nbsp;&nbsp; 0.00006799<br>
          &nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 17 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00002056&nbsp;&nbsp;&nbsp;
          0.00002091&nbsp;&nbsp;&nbsp; 0.00006192<br>
          &nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 18 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp; -0.00004895&nbsp;&nbsp; -0.00000474&nbsp;&nbsp;
          -0.00012101<br>
          &nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 19 type&nbsp; 2&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp; -0.00000208&nbsp;&nbsp;
          -0.00000446&nbsp;&nbsp;&nbsp; 0.00003292<br>
          &nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 20 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp; -0.00003014&nbsp;&nbsp;&nbsp;
          0.00001046&nbsp;&nbsp;&nbsp; 0.00005040<br>
          &nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 21 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp; -0.00000665&nbsp;&nbsp;&nbsp; 0.00000012&nbsp;&nbsp;
          -0.00002565<br>
          &nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 22 type&nbsp; 2&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00004366&nbsp;&nbsp;&nbsp; 0.00001375&nbsp;&nbsp;
          -0.00003139<br>
          &nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 23 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00001038&nbsp;&nbsp;&nbsp; 0.00000199&nbsp;&nbsp;
          -0.00004611<br>
          &nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 24 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00000001&nbsp;&nbsp;
          -0.00002044&nbsp;&nbsp;&nbsp; 0.00003004<br>
          &nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 25 type&nbsp; 2&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp; -0.00000668&nbsp;&nbsp;
          -0.00002989&nbsp;&nbsp;&nbsp; 0.00007813<br>
          &nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 26 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00003097&nbsp;&nbsp;
          -0.00004437&nbsp;&nbsp;&nbsp; 0.00003757<br>
          &nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 27 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00001552&nbsp;&nbsp;&nbsp; 0.00000595&nbsp;&nbsp;
          -0.00005710<br>
          &nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 28 type&nbsp; 2&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00001236&nbsp;&nbsp;&nbsp;
          0.00004340&nbsp;&nbsp;&nbsp; 0.00012950<br>
          &nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 29 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp; -0.00002372&nbsp;&nbsp;&nbsp; 0.00002505&nbsp;&nbsp;
          -0.00007569<br>
          &nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 30 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00003122&nbsp;&nbsp;&nbsp;
          0.00000291&nbsp;&nbsp;&nbsp; 0.00000276<br>
          &nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 31 type&nbsp; 2&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp; -0.00001390&nbsp;&nbsp; -0.00010264&nbsp;&nbsp;
          -0.00009180<br>
          &nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 32 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00004590&nbsp;&nbsp;&nbsp;
          0.00003207&nbsp;&nbsp;&nbsp; 0.00003011<br>
          &nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 33 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp; -0.00002671&nbsp;&nbsp;&nbsp; 0.00005738&nbsp;&nbsp;
          -0.00001122<br>
          &nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 34 type&nbsp; 2&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp; -0.00001380&nbsp;&nbsp; -0.00001243&nbsp;&nbsp;
          -0.00005321<br>
          &nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 35 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp; -0.00002033&nbsp;&nbsp;&nbsp;
          0.00001199&nbsp;&nbsp;&nbsp; 0.00006508<br>
          &nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 36 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00001161&nbsp;&nbsp; -0.00000439&nbsp;&nbsp;
          -0.00003033<br>
          &nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 37 type&nbsp; 2&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp; -0.00008127&nbsp;&nbsp;&nbsp;
          0.00002129&nbsp;&nbsp;&nbsp; 0.00005130<br>
          &nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 38 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp; -0.00004667&nbsp;&nbsp; -0.00008260&nbsp;&nbsp;
          -0.00004117<br>
          &nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 39 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00000793&nbsp;&nbsp;&nbsp; 0.00013767&nbsp;&nbsp;
          -0.00000675<br>
          &nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 40 type&nbsp; 2&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00001123&nbsp;&nbsp;
          -0.00004604&nbsp;&nbsp;&nbsp; 0.00006549<br>
          &nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 41 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00003825&nbsp;&nbsp;&nbsp; 0.00001125&nbsp;&nbsp;
          -0.00000180<br>
          &nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 42 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp; -0.00005610&nbsp;&nbsp;&nbsp; 0.00001922&nbsp;&nbsp;
          -0.00008202<br>
          &nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 43 type&nbsp; 2&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp; -0.00001184&nbsp;&nbsp;&nbsp; 0.00000509&nbsp;&nbsp;
          -0.00002456<br>
          &nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 44 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp; -0.00005112&nbsp;&nbsp;
          -0.00001419&nbsp;&nbsp;&nbsp; 0.00000597<br>
          &nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 45 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00009312&nbsp;&nbsp; -0.00000567&nbsp;&nbsp;
          -0.00002162<br>
          &nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 46 type&nbsp; 2&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp; -0.00002442&nbsp;&nbsp;
          -0.00003492&nbsp;&nbsp;&nbsp; 0.00009560<br>
          &nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 47 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp; -0.00005855&nbsp;&nbsp;&nbsp; 0.00002530&nbsp;&nbsp;
          -0.00010962<br>
          &nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 48 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00008502&nbsp;&nbsp; -0.00001219&nbsp;&nbsp;
          -0.00005446<br>
          &nbsp;<br>
          &nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp;&nbsp;&nbsp; 0.000643&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp;
          0.000037</em></div>
      <div><em></em>&nbsp;</div>
      <div><em>entering subroutine stress ...<br>
          &nbsp;<br>
          &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; total&nbsp;&nbsp; stress&nbsp; (Ry/bohr**3)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
          (kbar)&nbsp;&nbsp;&nbsp;&nbsp; P=&nbsp;&nbsp; -1.92<br>
          &nbsp; -0.00001798&nbsp;&nbsp; 0.00000011&nbsp;&nbsp; 0.00000000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -2.64&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
          0.02&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00<br>
          &nbsp;&nbsp; 0.00000011&nbsp; -0.00001708&nbsp;&nbsp; 0.00000001&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.02&nbsp;&nbsp;&nbsp;&nbsp;
          -2.51&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00<br>
          &nbsp;&nbsp; 0.00000000&nbsp;&nbsp; 0.00000001&nbsp; -0.00000407&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
          0.00&nbsp;&nbsp;&nbsp;&nbsp; -0.60<br>
          &nbsp;<br>
          &nbsp;<br>
          &nbsp;&nbsp;&nbsp;&nbsp; The maximum number of steps has been reached.<br>
          &nbsp;<br>
          &nbsp;&nbsp;&nbsp;&nbsp; End of BFGS Geometry Optimization</em></div>
      <div>&nbsp;</div>
      <div>&nbsp;Ie. &nbsp;the code saying <em>The maximum number of steps has
          been reached.</em> the relaxation is not&nbsp;convergence.&nbsp;I am
        puzzled&nbsp;the&nbsp;&nbsp;etot_conv_thr does not meet ?? The&nbsp;etot_conv_thr&nbsp;
        is not high, and the nstep is already 200.&nbsp;&nbsp;try a lower
        etot_conv_thr?</div>
      <div>&nbsp;</div>
      <div>i wish someone could have a look at the output. Thanks a lot
        !</div>
      <div>&nbsp;</div>
      <div>&nbsp;</div>
      <div>&nbsp;</div>
      <div>&nbsp;</div>
      <div><br>
        -- <br>
        Yuehua Xu<br>
        Group of Computational Condensed Matter Physics, <br>
        National Laboratory of Solid State Microstructures and
        Department of Physics, <br>
        Nanjing University,<br>
        &nbsp;Nanjing 210093,<br>
        &nbsp;P. R. China<br>
        <br>
      </div>
      <br>
      <fieldset class="mimeAttachmentHeader"></fieldset>
      <br>
      <pre wrap="">_______________________________________________
Pw_forum mailing list
<a class="moz-txt-link-abbreviated" href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a>
<a class="moz-txt-link-freetext" href="http://www.democritos.it/mailman/listinfo/pw_forum">http://www.democritos.it/mailman/listinfo/pw_forum</a>
</pre>
    </blockquote>
    <br>
    <div class="moz-signature">-- <br>
      <center>
        <font color="#946e8c">
          <em>
            Be the change you wish to see in the world<br>
          </em>
          <small>
            <font color="#946e8c">
              &#8212; Mahatma Gandhi &#8212;<br>
            </font></small>
          <br>
          <a
            href="http://www.linkedin.com/pub/%C3%A9ric-germaneau/30/931/986">Dr.
            &Eacute;ric Germaneau</a><br>
          <br>
          Graduate University of Chinese Academy of Sciences<br>
          College of Physical Sciences<br>
          Yuquan Road 19A<br>
          Beijing 100049<br>
          China<br>
          <br>
          <small>
            <em>
              <font style="color: #498848" face="Times">
                Please, if possible, don't send me MS Word or PowerPoint
                attachments<br>
                Why? See: <a
                  href="http://www.gnu.org/philosophy/no-word-attachments.html"
                  style="text-decoration:none"><font color="#498848">http://www.gnu.org/philosophy/no-word-attachments.html</font></a>
              </font></em>
          </small>
        </font></center>
    </div>
  </body>
</html>