<div>Yes. For the first case, you need to increase the nstep.</div>
<div>&nbsp;</div>
<div><br><br>&nbsp;</div>
<div><span class="gmail_quote">2007/11/30, xu yuehua &lt;<a href="mailto:njuxuyuehua@gmail.com">njuxuyuehua@gmail.com</a>&gt;:</span>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">
<div>hi all :</div>
<div>i have a question about the out file of relax :</div>
<div>for example :sometimes :the end of file is like:</div>
<div>.........</div>
<p>&nbsp; convergence has been achieved</p>
<p>&nbsp;&nbsp;&nbsp;&nbsp; Forces acting on atoms (Ry/au):</p>
<p>&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp;&nbsp; 1 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp; -0.00016739&nbsp;&nbsp;&nbsp; 0.00006788&nbsp;&nbsp;&nbsp; 0.00006825<br>&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp;&nbsp; 2 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00000361&nbsp;&nbsp;&nbsp; 0.00001673&nbsp;&nbsp; -0.00004580<br>&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp;&nbsp; 3 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00007761&nbsp;&nbsp;&nbsp; 0.00005882&nbsp;&nbsp;&nbsp; 
0.00004672<br>&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp;&nbsp; 4 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00003052&nbsp;&nbsp;&nbsp; 0.00004707&nbsp;&nbsp;&nbsp; 0.00002272<br>&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp;&nbsp; 5 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00001339&nbsp;&nbsp;&nbsp; 0.00015088&nbsp;&nbsp; -0.00012165<br>&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp;&nbsp; 6 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp; -0.00010455&nbsp;&nbsp; -0.00002068
 &nbsp;&nbsp; -0.00010771<br>&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp;&nbsp; 7 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp; -0.00009487&nbsp;&nbsp; -0.00001394&nbsp;&nbsp;&nbsp; 0.00006473<br>&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp;&nbsp; 8 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp; -0.00000214&nbsp;&nbsp; -0.00007630&nbsp;&nbsp;&nbsp; 0.00001052<br>&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp;&nbsp; 9 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00002701
&nbsp;&nbsp;&nbsp; 0.00001038&nbsp;&nbsp;&nbsp; 0.00000479<br>&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 10 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00008230&nbsp;&nbsp;&nbsp; 0.00004823&nbsp;&nbsp;&nbsp; 0.00012935<br>&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 11 type&nbsp; 2&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00015243&nbsp;&nbsp; -0.00016928&nbsp;&nbsp; -0.00010298<br>&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 12 type&nbsp; 2&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp; -
0.00004574 &nbsp;&nbsp; -0.00002796&nbsp;&nbsp;&nbsp; 0.00000972<br>&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 13 type&nbsp; 2&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00003206&nbsp;&nbsp;&nbsp; 0.00001371&nbsp;&nbsp; -0.00009648<br>&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 14 type&nbsp; 2&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00000462&nbsp;&nbsp; -0.00000415&nbsp;&nbsp;&nbsp; 0.00003920<br>&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 15 type&nbsp; 2&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp; - 
0.00000887&nbsp;&nbsp; -0.00010141&nbsp;&nbsp;&nbsp; 0.00007861</p>
<p>&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp;&nbsp;&nbsp; 0.000506&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000024</p>
<p><br>&nbsp;&nbsp;&nbsp;&nbsp; entering subroutine stress ...</p>
<p>&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; -0.60<br>&nbsp; -0.00001045&nbsp; -0.00000009&nbsp;&nbsp; 0.00000000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -1.54&nbsp;&nbsp;&nbsp;&nbsp; -0.01&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00<br>&nbsp; -0.00000009&nbsp; -0.00001053&nbsp; -0.00000015&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -0.01&nbsp;&nbsp;&nbsp;&nbsp; - 
1.55&nbsp;&nbsp;&nbsp;&nbsp; -0.02<br>&nbsp;&nbsp; 0.00000000&nbsp; -0.00000015&nbsp;&nbsp; 0.00000867&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00&nbsp;&nbsp;&nbsp;&nbsp; -0.02&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1.27</p>
<p><br>&nbsp;&nbsp;&nbsp;&nbsp;<font color="#ff0000"> The maximum number of steps has been reached.</font></p>
<p><font color="#ff0000">&nbsp;&nbsp;&nbsp;&nbsp; End of BFGS Geometry Optimization</font></p>
<p><font color="#ff0000">and sometimes&nbsp; the same part of the outfile:<br></font></p>
<p>&nbsp;<font color="#000000">convergence has been achieved</font></p>
<p><font color="#000000">&nbsp;&nbsp;&nbsp;&nbsp; Forces acting on atoms (Ry/au):</font></p>
<p><font color="#000000">&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp;&nbsp; 1 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00000253&nbsp;&nbsp; -0.00000639&nbsp;&nbsp; -0.00002497<br>&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp;&nbsp; 2 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp; -0.00000977&nbsp;&nbsp;&nbsp; 0.00007654&nbsp;&nbsp;&nbsp; 0.00003837<br>&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp;&nbsp; 3 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00001686
 &nbsp;&nbsp; -0.00000393&nbsp;&nbsp; -0.00003744<br>&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp;&nbsp; 4 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp; -0.00002517&nbsp;&nbsp; -0.00008976&nbsp;&nbsp; -0.00005193<br>&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp;&nbsp; 5 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00002007&nbsp;&nbsp;&nbsp; 0.00000939&nbsp;&nbsp; -0.00004732<br>&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp;&nbsp; 6 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 
0.00007598&nbsp;&nbsp; -0.00004787&nbsp;&nbsp; -0.00003895<br>&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp;&nbsp; 7 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp; -0.00004187&nbsp;&nbsp; -0.00005210&nbsp;&nbsp; -0.00004987<br>&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp;&nbsp; 8 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00007160&nbsp;&nbsp;&nbsp; 0.00002922&nbsp;&nbsp; -0.00003445<br>&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp;&nbsp; 9 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp; - 
0.00002615&nbsp;&nbsp; -0.00003529&nbsp;&nbsp;&nbsp; 0.00005795<br>&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 10 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp; -0.00006129&nbsp;&nbsp;&nbsp; 0.00004927&nbsp;&nbsp;&nbsp; 0.00009998<br>&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 11 type&nbsp; 2&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp; -0.00005381&nbsp;&nbsp;&nbsp; 0.00007155&nbsp;&nbsp;&nbsp; 0.00007934<br>&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 12 type&nbsp; 2&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 
0.00000162&nbsp;&nbsp;&nbsp; 0.00002190&nbsp;&nbsp; -0.00004511<br>&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 13 type&nbsp; 2&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00000934&nbsp;&nbsp; -0.00002334&nbsp;&nbsp; -0.00004641<br>&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 14 type&nbsp; 2&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00003732&nbsp;&nbsp; -0.00001936&nbsp;&nbsp;&nbsp; 0.00003124<br>&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 15 type&nbsp; 2&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp; - 
0.00001727&nbsp;&nbsp;&nbsp; 0.00002015&nbsp;&nbsp;&nbsp; 0.00006958</font></p>
<p><font color="#000000">&nbsp;&nbsp;&nbsp;&nbsp; Total force =&nbsp;&nbsp;&nbsp;&nbsp; 0.000311&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000021</font></p>
<p><br><font color="#000000">&nbsp;&nbsp;&nbsp;&nbsp; entering subroutine stress ...</font></p>
<p><font color="#000000">&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.13<br>&nbsp; -0.00001084&nbsp;&nbsp; 0.00000003&nbsp;&nbsp; 0.00000019&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -1.59&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.01&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.03<br>&nbsp;&nbsp; 0.00000003&nbsp; -0.00001082&nbsp; -0.00000037
 &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.01&nbsp;&nbsp;&nbsp;&nbsp; -1.59&nbsp;&nbsp;&nbsp;&nbsp; -0.05<br>&nbsp;&nbsp; 0.00000019&nbsp; -0.00000037&nbsp; -0.00000140&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.03&nbsp;&nbsp;&nbsp;&nbsp; -0.05&nbsp;&nbsp;&nbsp;&nbsp; -0.21</font></p>
<p><br><font color="#000000">&nbsp;&nbsp;&nbsp;</font><font color="#ff0000">&nbsp; bfgs converged in&nbsp; 27 scf cycles and&nbsp; 24 bfgs steps</font></p>
<p><font color="#ff0000">&nbsp;&nbsp;&nbsp;&nbsp; End of BFGS Geometry Optimization</font></p>
<p><font color="#ff0000">&nbsp;&nbsp;&nbsp;&nbsp; Final energy =&nbsp;&nbsp;&nbsp; -172.2464671441 Ry</font></p>
<div>the red are different ,so i worried about the former is not &nbsp;convergence,because i set the step is 50.</div>
<div>is that so ?</div>
<div>i&nbsp;need your help .</div>
<div>thank you &nbsp;<br clear="all"><br>-- <br>Xu Yuehua<br>physics Department of Nanjing university<br>China </div><br>_______________________________________________<br>Pw_forum mailing list<br><a onclick="return top.js.OpenExtLink(window,event,this)" href="mailto:Pw_forum@pwscf.org">
Pw_forum@pwscf.org</a><br><a onclick="return top.js.OpenExtLink(window,event,this)" href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br><br></blockquote>
</div><br><br clear="all"><br>-- <br>------------------<br>Zhu feng<br>Department of Physics, <br>Anhui University, Hefei, China