Wow, the  convergence is good when the "dt" flag is further
reduced to 30 considering my special system in which it only contains
some light elements. Thanks Dr. Axel and Cesar again for your
suggestion -:)<br>
<br>
==========<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; -7.01<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; -6.21<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; -5.28<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; -4.38<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; -3.39<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; -2.40<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.49<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; -0.61<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;&nbsp; 0.24<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;&nbsp; 0.59<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; -0.00<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;&nbsp; 0.14<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;&nbsp; 0.15<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;&nbsp; 0.21<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;&nbsp; 0.24<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;&nbsp; 0.19<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;&nbsp; 0.15<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;&nbsp; 0.13<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;&nbsp; 0.10<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;&nbsp; 0.07<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;&nbsp; 0.07<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;&nbsp; 0.07<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;&nbsp; 0.07<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;&nbsp; 0.03<br>
...<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Total force
=&nbsp;&nbsp;&nbsp;&nbsp; 0.001176&nbsp;&nbsp;&nbsp;&nbsp; Total SCF
correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000004<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Total force
=&nbsp;&nbsp;&nbsp;&nbsp; 0.001374&nbsp;&nbsp;&nbsp;&nbsp; Total SCF
correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000004<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Total force
=&nbsp;&nbsp;&nbsp;&nbsp; 0.001858&nbsp;&nbsp;&nbsp;&nbsp; Total SCF
correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000006<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Total force
=&nbsp;&nbsp;&nbsp;&nbsp; 0.001973&nbsp;&nbsp;&nbsp;&nbsp; Total SCF
correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000003<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Total force
=&nbsp;&nbsp;&nbsp;&nbsp; 0.001762&nbsp;&nbsp;&nbsp;&nbsp; Total SCF
correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000004<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Total force
=&nbsp;&nbsp;&nbsp;&nbsp; 0.001475&nbsp;&nbsp;&nbsp;&nbsp; Total SCF
correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000004<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Total force
=&nbsp;&nbsp;&nbsp;&nbsp; 0.001355&nbsp;&nbsp;&nbsp;&nbsp; Total SCF
correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000004<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Total force
=&nbsp;&nbsp;&nbsp;&nbsp; 0.001877&nbsp;&nbsp;&nbsp;&nbsp; Total SCF
correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000005<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Total force
=&nbsp;&nbsp;&nbsp;&nbsp; 0.001823&nbsp;&nbsp;&nbsp;&nbsp; Total SCF
correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000003<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Total force
=&nbsp;&nbsp;&nbsp;&nbsp; 0.001150&nbsp;&nbsp;&nbsp;&nbsp; Total SCF
correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000005<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Total force
=&nbsp;&nbsp;&nbsp;&nbsp; 0.000929&nbsp;&nbsp;&nbsp;&nbsp; Total SCF
correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000003<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Total force
=&nbsp;&nbsp;&nbsp;&nbsp; 0.000917&nbsp;&nbsp;&nbsp;&nbsp; Total SCF
correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000002<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Total force
=&nbsp;&nbsp;&nbsp;&nbsp; 0.000931&nbsp;&nbsp;&nbsp;&nbsp; Total SCF
correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000003<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Total force
=&nbsp;&nbsp;&nbsp;&nbsp; 0.001003&nbsp;&nbsp;&nbsp;&nbsp; Total SCF
correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000003<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Total force
=&nbsp;&nbsp;&nbsp;&nbsp; 0.000736&nbsp;&nbsp;&nbsp;&nbsp; Total SCF
correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000002<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Total force
=&nbsp;&nbsp;&nbsp;&nbsp; 0.000728&nbsp;&nbsp;&nbsp;&nbsp; Total SCF
correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000002<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Total force
=&nbsp;&nbsp;&nbsp;&nbsp; 0.000749&nbsp;&nbsp;&nbsp;&nbsp; Total SCF
correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000001<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Total force
=&nbsp;&nbsp;&nbsp;&nbsp; 0.000747&nbsp;&nbsp;&nbsp;&nbsp; Total SCF
correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000002<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Total force
=&nbsp;&nbsp;&nbsp;&nbsp; 0.000692&nbsp;&nbsp;&nbsp;&nbsp; Total SCF
correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000003<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Total force
=&nbsp;&nbsp;&nbsp;&nbsp; 0.000751&nbsp;&nbsp;&nbsp;&nbsp; Total SCF
correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000002<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Total force
=&nbsp;&nbsp;&nbsp;&nbsp; 0.000872&nbsp;&nbsp;&nbsp;&nbsp; Total SCF
correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000002<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Total force
=&nbsp;&nbsp;&nbsp;&nbsp; 0.000671&nbsp;&nbsp;&nbsp;&nbsp; Total SCF
correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000003<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Total force
=&nbsp;&nbsp;&nbsp;&nbsp; 0.000723&nbsp;&nbsp;&nbsp;&nbsp; Total SCF
correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000002<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Total force
=&nbsp;&nbsp;&nbsp;&nbsp; 0.000715&nbsp;&nbsp;&nbsp;&nbsp; Total SCF
correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000003<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Total force
=&nbsp;&nbsp;&nbsp;&nbsp; 0.000620&nbsp;&nbsp;&nbsp;&nbsp; Total SCF
correction =&nbsp;&nbsp;&nbsp;&nbsp; 0.000002<br>
...<br>
<br>
Best regards,<br>
C.H. <br clear="all"><br>-- <br>============================<br>C.H. Hu<br>Postdoctoral fellow<br>Chimie et Physico-Chimie appliquees<br>Institut Francais du Petrole (IFP)<br>Rueil-Malmaison, France<br>E_mail: <a href="mailto:chaohao.mse@gmail.com">
chaohao.mse@gmail.com</a><br>============================