
<html><head><style type="text/css"><!-- DIV {margin:0px;} --></style></head><body><div style="font-family:times new roman,new york,times,serif;font-size:12pt"><div>&nbsp;Dear all<br>for noncubic cell ,in EXX calculation How we can choose the nq points?<br>For TiO2 I used this input file but the code stop in midway without any message? and it didn't print total energy?<br>&amp;control<br>&nbsp;&nbsp;&nbsp; calculation='scf'<br>&nbsp;&nbsp;&nbsp; restart_mode='from_scratch',<br>&nbsp;&nbsp;&nbsp; pseudo_dir = './',<br>&nbsp;&nbsp;&nbsp; prefix='tio2'<br>&nbsp;&nbsp;&nbsp; outdir='/p5/batch/kazempou/'<br>/<br>&nbsp;&amp;system<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ibrav = 6,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; celldm(1) =

 8.79,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; celldm(3) = 0.639,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nat = 6,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ntyp = 2,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ecutwfc = 25 ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ecutrho = 400 , input_dft='pbe0', nqx1 = 2, nqx2 = 2, nqx3 = 4,<br>&nbsp;/<br>&nbsp;&amp;electrons<br>&nbsp;&nbsp;&nbsp; diagonalization='david',<br>&nbsp;&nbsp;&nbsp; conv_thr =&nbsp; 1.0d-4<br>&nbsp;&nbsp;&nbsp; mixing_beta =

 0.5,<br>&nbsp;/<br>ATOMIC_SPECIES<br>&nbsp;&nbsp; Ti&nbsp;&nbsp; 47.86700&nbsp; ti.optgga2.fhi.UPF<br>&nbsp;&nbsp;&nbsp; O&nbsp;&nbsp; 15.99940&nbsp; o.optgga1.fhi.UPF<br>ATOMIC_POSITIONS crystal<br>Ti&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.<span class="yshortcuts" id="lw_1261481063_0">000000000</span>&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp;&nbsp; 0.000000000<br>Ti&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.500000000&nbsp;&nbsp;&nbsp; 0.500000000&nbsp;&nbsp;&nbsp; 0.500000000<br>O&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.305131381&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.305131381&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.000000000<br>O&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -0.305131381&nbsp;&nbsp; -0.305131381&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.000000000<br>O&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.805131381&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.194868619&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;

 0.5<br>O&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;

 0.194868619&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.805131381&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.5<br>K_POINTS automatic<br>2 2 4&nbsp;&nbsp; 1 1 1<br>------------------------------<br>output file<br><br><br><br><br><br></div><br>&nbsp;&nbsp;&nbsp;&nbsp; total cpu time spent up to now is&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1.60 secs<br><br>&nbsp;&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp; -355.97704102 <span class="yshortcuts" id="lw_1261481063_1">Ry</span><br>&nbsp;&nbsp;&nbsp;&nbsp; Harris-Foulkes estimate&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp; -355.97835356 Ry<br>&nbsp;&nbsp;&nbsp;&nbsp; estimated scf accuracy&nbsp;&nbsp;&nbsp; &lt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00523443 Ry<br><br>&nbsp;&nbsp;&nbsp;&nbsp; iteration #&nbsp; 5&nbsp;&nbsp;&nbsp;&nbsp; ecut=&nbsp;&nbsp;&nbsp; 25.00 Ry&nbsp;&nbsp;&nbsp;&nbsp; beta=0.50<br>&nbsp;&nbsp;&nbsp;&nbsp;

 Davidson diagonalization with

 overlap<br>&nbsp;&nbsp;&nbsp;&nbsp; ethr =&nbsp; 1.09E-05,&nbsp; avg # of iterations =&nbsp; 4.0<br><br>&nbsp;&nbsp;&nbsp;&nbsp; total cpu time spent up to now is&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1.84 secs<br><br>&nbsp;&nbsp;&nbsp;&nbsp; total energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp; -355.97795182 Ry<br>&nbsp;&nbsp;&nbsp;&nbsp; Harris-Foulkes estimate&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp; -355.97903987 Ry<br>&nbsp;&nbsp;&nbsp;&nbsp; estimated scf accuracy&nbsp;&nbsp;&nbsp; &lt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00242696 Ry<br><br>&nbsp;&nbsp;&nbsp;&nbsp; iteration #&nbsp; 6&nbsp;&nbsp;&nbsp;&nbsp; ecut=&nbsp;&nbsp;&nbsp; 25.00 Ry&nbsp;&nbsp;&nbsp;&nbsp; beta=0.50<br>&nbsp;&nbsp;&nbsp;&nbsp; Davidson diagonalization with overlap<br>&nbsp;&nbsp;&nbsp;&nbsp; ethr =&nbsp; 5.06E-06,&nbsp; avg # of iterations =&nbsp; 3.0<br><br>&nbsp;&nbsp;&nbsp;&nbsp; total cpu time spent up to now

 is&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2.02 secs<br><br>&nbsp;&nbsp;&nbsp;&nbsp; End of self-consistent calculation<br><br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; k = 0.2500 0.2500 0.1956 (&nbsp;&nbsp; 920 PWs)&nbsp;&nbsp; bands (ev):<br><br>&nbsp;&nbsp; -50.3131 -50.2750 -26.6056 -26.5256 -26.5170 -26.5064 -26.4805 -26.4632<br>&nbsp;&nbsp;&nbsp; -9.2828&nbsp; -8.6924&nbsp; -8.2771&nbsp; -8.1594&nbsp;&nbsp; 4.0836&nbsp;&nbsp; 4.7890&nbsp;&nbsp; 5.3145&nbsp;&nbsp; 5.8431<br>&nbsp;&nbsp;&nbsp;&nbsp; 6.1571&nbsp;&nbsp; 6.8145&nbsp;&nbsp; 7.1870&nbsp;&nbsp; 7.4771&nbsp;&nbsp; 7.7110&nbsp;&nbsp; 8.6006&nbsp;&nbsp; 8.7972&nbsp;&nbsp; 9.2575<br><br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; k = 0.2500 0.2500 0.5869 (&nbsp;&nbsp; 908 PWs)&nbsp;&nbsp; bands (ev):<br><br>&nbsp;&nbsp; -50.2617 -50.2509 -26.7111 -26.6144 -26.5061 -26.4959 -26.4896 -26.4409<br>&nbsp;&nbsp;&nbsp; -8.7021&nbsp; -8.4646&nbsp; -8.4175&nbsp; -8.2262&nbsp;&nbsp;

 4.7053&nbsp;&nbsp; 5.0674&nbsp;&nbsp; 5.0763&nbsp;&nbsp; 5.4017<br>&nbsp;&nbsp;&nbsp;&nbsp; 6.0629&nbsp;&nbsp; 6.5874&nbsp;&nbsp; 7.1698&nbsp;&nbsp; 7.2399&nbsp;&nbsp; 7.5577&nbsp;&nbsp; 8.4242&nbsp;&nbsp; 8.4760&nbsp;&nbsp; 9.0196<br>&nbsp;0.892062058076385500 0.892062058076385500<br>EXX divergence (&nbsp;&nbsp; 2)=&nbsp;&nbsp;&nbsp; -108.3050&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.4000<br>&nbsp;&nbsp;&nbsp;&nbsp; exx_div&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; :&nbsp;&nbsp;&nbsp;&nbsp; 0.01s CPU<br>&nbsp; ! EXXALFA SET TO&nbsp; 0.250000000000000000<br><br><br><br><br><br>thanks a lot<br><br><br><br><br>Ali Kazempour<br><br><div><br>Fritz-Haber-Institut              fax   : ++49-30-8413 4701<br>der Max-Planck-Gesellschaft<br>Faradayweg 4-6                    e-mail: <span style="border-bottom: 1px dashed rgb(0, 102, 204); background: transparent none repeat scroll 0% 0%; cursor: pointer; -moz-background-clip: -moz-initial; -moz-background-origin: -moz-initial;

 -moz-background-inline-policy: -moz-initial;" class="yshortcuts" id="lw_1261481063_2">kazempou@fhi-berlin.mpg.de</span><br>D-14 195 Berlin-Dahlem / German<div><br></div></div><div>&nbsp;</div><div><br><div><br></div></div>

<!-- cg34.c4.mail.gq1.yahoo.com compressed/chunked Sat Dec 19 03:36:30 PST 2009 -->

</div><br>


      </body></html>