hi,<br>&nbsp; I performed other USPPs tests (blyp-van_ak, pw91-van_ak, pz-van_ak) and found that all these calculations with USPPs were not convergence both in press and Etot.For example, with blyp-van_ak( result is at the bottom of the mail).All these USPPs are downloaded from the homepage of PWSCF/pseudopotential. 
<br>&nbsp; However, this unconvergence situation disappeared with NCPP (PBE-TM-GGA, converted from fhi format PP using fhi2upf code). The result is at the bottom as well.<br>&nbsp; As Axel Kohlmeyer mentioned:<br>&nbsp; &quot;please try a different Li pseudopotential (
i.e. a non-ak one).<br>the one you are using seems to have some not yet understood<br>problems for higher-order properties: lattice constant from<br>EOS fit, MD and geometry optimizations were fine, phonons<br>were bad.&quot;
<br>&nbsp; And can anyone tell me what does the &#39;ak&#39;(in blyp-van_ak) mean? What kind of systems are these &#39;ak&#39; USPPs suitable for calculation? <br><br>WITH blyp-van_ak USPP:------------------------------------------------------
<br>&nbsp;&nbsp;&nbsp; ---!&nbsp;&nbsp; convergence scf caculation for Li2O&nbsp;&nbsp; !---<br>&nbsp;ecutwfc = 100 Ry,&nbsp; ecutrho = 320 Ry :<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -61.82021122 Ry<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; 
19.03<br>&nbsp;ecutwfc = 100 Ry,&nbsp; ecutrho = 360 Ry :<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -61.82317722 Ry<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.75<br>&nbsp;ecutwfc = 100 Ry,&nbsp; ecutrho = 400 Ry :
<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -61.82694596 Ry<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; 4.67<br>&nbsp;ecutwfc = 100 Ry,&nbsp; ecutrho = 600 Ry :<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -61.82089409
 Ry<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; 4.29<br>&nbsp;ecutwfc = 100 Ry,&nbsp; ecutrho = 800 Ry :<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -61.81983678 Ry<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; 
19.40<br>&nbsp;ecutwfc = 100 Ry,&nbsp; ecutrho = 1000 Ry :<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -61.82031614 Ry<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; 22.33<br>&nbsp;ecutwfc = 200 Ry,&nbsp; ecutrho = 800 Ry :
<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -61.82432402 Ry<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; 9.13<br>&nbsp;ecutwfc = 200 Ry,&nbsp; ecutrho = 1000 Ry :<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -61.82480674
 Ry<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; 12.09<br>&nbsp;ecutwfc = 200 Ry,&nbsp; ecutrho = 1600 Ry :<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -61.82437317 Ry<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; 
6.83<br>&nbsp;ecutwfc = 200 Ry,&nbsp; ecutrho = 1800 Ry :<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -61.82423710 Ry<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; 5.72<br>-------------------------------------------------------------------------------
<br>&nbsp; &nbsp;<br>WITH PBE-MT-GGA NCPP:-----------------------------------------------------------<br>&nbsp;&nbsp;&nbsp; ---!&nbsp;&nbsp; convergence scf caculation for Li2O&nbsp;&nbsp; !---<br>&nbsp;ecutwfc = ecut = 50 Ry :<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -33.53719310
 Ry<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= -329.31<br>&nbsp;ecutwfc = ecut = 80 Ry :<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -33.60362083 Ry<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; -
46.81<br>&nbsp;ecutwfc = ecut = 100 Ry :<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -33.60420031 Ry<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; -48.44<br>&nbsp;ecutwfc = ecut = 140 Ry :<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -
33.60518304 Ry<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; -46.36<br>&nbsp;ecutwfc = ecut = 180 Ry :<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -33.60525632 Ry<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; -
45.97<br>&nbsp;ecutwfc = ecut = 200 Ry :<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -33.60552604 Ry<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; -46.30<br>&nbsp;ecutwfc = ecut = 240 Ry :<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -
33.60553853 Ry<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; -46.24<br>&nbsp;ecutwfc = ecut = 280 Ry :<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -33.60553917 Ry<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; -
46.23<br>&nbsp;ecutwfc = ecut = 300 Ry :<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -33.60553970 Ry<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; -46.24<br>&nbsp;ecutwfc = ecut = 340 Ry :<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -
33.60553987 Ry<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; -46.23<br>&nbsp;ecutwfc = ecut = 400 Ry :<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -33.60554013 Ry<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; -
46.23<br>---------------------------------------------------------------------------------