<div>The system is Li2O, and the input file and the result is :</div>
<div>&nbsp;</div>
<div>&nbsp;&nbsp;&nbsp; ---!&nbsp;&nbsp; convergence scf caculation for Li2O&nbsp;&nbsp; !---</div>
<div>&nbsp;ecut=80 :<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -61.89999202 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; 0.59<br>&nbsp;ecut=100 :<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -61.90651146 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; -5.39<br>&nbsp;ecut=140 :<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -61.90331161 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; 
0.97<br>&nbsp;ecut=180 :<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -61.90482904 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; 10.64<br>&nbsp;ecut=200 :<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -61.90542701
 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; -4.55<br>&nbsp;ecut=240 :<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -61.90875043 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.92<br>&nbsp;ecut=280 :<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -61.90854533 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; 1.68<br>&nbsp;ecut=300 :<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -61.90926981
 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; -1.18<br>&nbsp;ecut=320 :<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -61.90990283 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; 
2.17<br>&nbsp;ecut=350 :<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -61.91013980 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.51<br>&nbsp;ecut=380 :<br>!&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -61.91023915
 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.47<br>&nbsp;</div>
<div>---------------------------------------------------------</div>
<div>#!/bin/sh&nbsp;<br>#run from directory where this script is<br>cd `echo $0 | sed &#39;s/\(.*\)\/.*/\1/&#39;` # extract pathname<br>&nbsp;<br># check whether echo has the -e option<br>if test &quot;`echo -e`&quot; = &quot;-e&quot; ; then ECHO=echo ; else ECHO=&quot;echo -e&quot; ; fi
<br>&nbsp;<br>touch totResult.out<br>&nbsp;<br>cat &gt; totResult.out &lt;&lt;EOF<br>&nbsp;&nbsp;&nbsp; ---!&nbsp;&nbsp; convergence scf caculation for Li2O&nbsp;&nbsp; !---<br>EOF<br>$ECHO<br>$ECHO &quot;-----$PWD : starting&quot;<br>date<br>$ECHO <br>$ECHO &quot;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ---!&nbsp;&nbsp; convergence scf caculation for Li2O&nbsp;&nbsp; !---&quot;
<br>&nbsp;<br>for i in 80 100 120 .......<br>do<br>&nbsp;&nbsp; cat &gt; li2O.ecut$i.scf.in &lt;&lt;EOF<br>&nbsp;&amp;control<br>&nbsp;&nbsp;&nbsp; calculation=&#39;scf&#39;<br>&nbsp;&nbsp;&nbsp; restart_mode=&#39;from_scratch&#39;,<br>&nbsp;&nbsp;&nbsp; tstress = .true.<br>&nbsp;&nbsp;&nbsp; tprnfor = .true.
<br>&nbsp;&nbsp;&nbsp; prefix=&#39;li2O&#39;,<br>&nbsp;&nbsp;&nbsp; pseudo_dir = &#39;/home/puma02/sunviva/qgr/program/pseudo_all&#39;,<br>&nbsp;&nbsp;&nbsp; outdir=&#39;/home/puma02/sunviva/qgr/work/pwscf_li2O/pw91-van_ak/cvg_ecut/temp&#39;<br>&nbsp;/<br>&nbsp;&amp;system&nbsp;&nbsp;&nbsp; 
<br>&nbsp;&nbsp;&nbsp; ibrav=&nbsp; 2,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; A =&nbsp; 4.61, <br>&nbsp;&nbsp;&nbsp;&nbsp; nat =&nbsp; 3,<br>&nbsp;&nbsp;&nbsp; ntyp =&nbsp; 2,<br>&nbsp;ecutwfc = $i<br>&nbsp;/<br>&nbsp;&amp;electrons<br>&nbsp;&nbsp;&nbsp; conv_thr =&nbsp; 1.0d-10<br>&nbsp;&nbsp;&nbsp; mixing_beta = 0.7<br>&nbsp;/<br>ATOMIC_SPECIES<br>&nbsp;Li&nbsp;&nbsp;&nbsp; 6.941&nbsp;&nbsp;&nbsp;&nbsp; Li.pw91-s-van_ak.UPF
<br>&nbsp; O&nbsp;&nbsp; 15.9994&nbsp;&nbsp;&nbsp;&nbsp; O.pw91-van_ak.UPF<br>ATOMIC_POSITIONS<br>&nbsp; O&nbsp; 0.00 0.00 0.00 <br>&nbsp;Li&nbsp; 0.25 0.25 0.25 <br>&nbsp;Li&nbsp; 0.75 0.75 0.75<br>K_POINTS {automatic}<br>&nbsp;12&nbsp; 12&nbsp; 12&nbsp;&nbsp; 0 0 0<br>&nbsp;<br>EOF<br>&nbsp;<br>$ECHO <br>#-----------------$ECHO &quot; SCF calculation for Li2O&nbsp;&nbsp; kpt = $i x $i x $i&quot; 
<br>$ECHO &quot; SCF calculation for Li2O&nbsp;&nbsp;&nbsp; ecut = $i &quot;<br>$ECHO<br>&nbsp;<br>/home/puma02/sunviva/qgr/program/espresso-3.2/bin/pw.x &lt; li2O.ecut$i.scf.in &gt; li2O.ecut$i.scf.out <br>&nbsp;<br>cat &gt;&gt; totResult.out&lt;&lt;EOF
<br>&nbsp;ecut=$i :<br>EOF<br>&nbsp;<br>&nbsp;grep -e &#39;!&#39;&nbsp; &lt; li2O.ecut$i.scf.out &gt;&gt; totResult.out<br>&nbsp;grep -e &#39;P=&#39; &lt; li2O.ecut$i.scf.out &gt;&gt; totResult.out<br>done<br>&nbsp;<br>$ECHO &quot;------END TIME: &quot; 
<br>date<br>$ECHO<br>$ECHO&nbsp;<br>&nbsp;<br>&nbsp;<br>&nbsp;<br>&nbsp;</div>