<div>&nbsp;</div>
<div>about why i use etot_conv_thr<br>as well as &nbsp;force_conv_thr= 1.d-6 Ry</div>
<div>as usual,in vasp and castep software,about 1.d-4 ev&nbsp;is ok</div>
<div>why it is so different between these DFT software ?</div>
<div>i don not understand.</div>
<div>need your help</div>
<div><br><br>&nbsp;</div>
<div><span class="gmail_quote">2007/7/13, Stefano de Gironcoli &lt;<a onclick="return top.js.OpenExtLink(window,event,this)" href="mailto:degironc@sissa.it" target="_blank">degironc@sissa.it</a>&gt;:</span> 
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">a tipical bond stiffness, if I recall correctly, is in the range of 0.1-1.0<br>Ry/bohr^2 thus a threshold of 
1.d-6 Ry/bohr in the force means that you are<br>asking a numerical accuracy in your interatomic distances of the order of<br>1.d-5-1.d-6 bohr... are you sure you really need it ?<br>best regards,<br>stefano deGironcoli, SISSA &amp; DEMOCRITOS 
<br><br>Quoting xu yuehua &lt;<a onclick="return top.js.OpenExtLink(window,event,this)" href="mailto:njuxuyuehua@gmail.com" target="_blank">njuxuyuehua@gmail.com</a>&gt;:<br><br>&gt; hello all<br>&gt; i am doing relax.<br>
&gt; in my inputfile,<br>&gt; i use&nbsp;Bohr&gt; and the output is<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;atom&nbsp;&nbsp; 1 type&nbsp;&nbsp;1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00000290&nbsp;&nbsp;&nbsp;&nbsp;0.00000095&nbsp;&nbsp; -0.00000364<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;atom&nbsp;&nbsp; 2 type&nbsp;&nbsp;1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00000081&nbsp;&nbsp; -0.00000108&nbsp;&nbsp;&nbsp;&nbsp;0.00000033
<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;atom&nbsp;&nbsp; 3 type&nbsp;&nbsp;1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp;-0.00000435&nbsp;&nbsp; -0.00000425&nbsp;&nbsp; -0.00000160<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;atom&nbsp;&nbsp; 4 type&nbsp;&nbsp;1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp;- 0.00000312&nbsp;&nbsp; -0.00000321&nbsp;&nbsp; -0.00000080<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;atom&nbsp;&nbsp; 5 type&nbsp;&nbsp;1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp;-0.00000446
&nbsp;&nbsp; -0.00000036&nbsp;&nbsp; -0.00000273<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;atom&nbsp;&nbsp; 6 type&nbsp;&nbsp;1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp;-0.00000048&nbsp;&nbsp;&nbsp;&nbsp;0.00000230&nbsp;&nbsp;&nbsp;&nbsp;0.00000112<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;atom&nbsp;&nbsp; 7 type&nbsp;&nbsp;1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00000164&nbsp;&nbsp;&nbsp;&nbsp;0.00000118&nbsp;&nbsp;&nbsp;&nbsp;0.00000060<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;atom&nbsp;&nbsp; 8 type&nbsp;&nbsp;1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 
0.00000131&nbsp;&nbsp; -0.00000040&nbsp;&nbsp;&nbsp;&nbsp;0.00000145<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;atom&nbsp;&nbsp; 9 type&nbsp;&nbsp;2&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp;-0.00000362&nbsp;&nbsp;&nbsp;&nbsp;0.00000154&nbsp;&nbsp; -0.00000201<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;atom&nbsp;&nbsp;10 type&nbsp;&nbsp;2&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp;- 0.00000029&nbsp;&nbsp;&nbsp;&nbsp;0.00000191&nbsp;&nbsp;&nbsp;&nbsp;0.00000302<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;atom&nbsp;&nbsp;11 type&nbsp;&nbsp;2&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 
0.00000576&nbsp;&nbsp; -0.00000134&nbsp;&nbsp;&nbsp;&nbsp;0.00000194<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;atom&nbsp;&nbsp;12 type&nbsp;&nbsp;2&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00000390&nbsp;&nbsp;&nbsp;&nbsp;0.00000276&nbsp;&nbsp;&nbsp;&nbsp;0.00000233<br>&gt;<br>&gt; Total force =&nbsp;&nbsp;&nbsp;&nbsp; 0.000015&nbsp;&nbsp;&nbsp;&nbsp; Total SCF correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000054<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;SCF correction compared to forces is too large, reduce conv_thr
<br>&gt;<br>&gt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;The maximum number of steps has been reached.<br>&gt; my qustion is ,force at every atom is ok ,but why&nbsp;&nbsp;the program asks to <br>&gt; reduce conv_thr<br>&gt; my conv_thr is not enough small ?<br>&gt; or my total relax step (50) is not adequate?
<br>&gt; is there any principle about these parameter?<br>&gt; need your help<br>&gt; thanks a lot<br>&gt; --<br>&gt; Xu Yuehua<br>&gt; physics Department of Nanjing university<br>&gt; China<br>&gt;<br><br><br><br><br>----------------------------------------------------------------
<br>SISSA Webmail <a onclick="return top.js.OpenExtLink(window,event,this)" href="https://webmail.sissa.it/" target="_blank">https://webmail.sissa.it/</a><br>Powered by Horde <a onclick="return top.js.OpenExtLink(window,event,this)" href="http://www.horde.org/" target="_blank">
http://www.horde.org/</a><br><br>_______________________________________________<br>Pw_forum mailing list<br><a onclick="return top.js.OpenExtLink(window,event,this)" href="mailto:Pw_forum@pwscf.org" target="_blank">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></blockquote></div><br>
<br clear="all"><br>-- <br>Xu Yuehua<br>physics Department of Nanjing university <br>China