<p style="margin: 0px; padding: 0px;" mce_style="margin: 0px; padding: 0px;">Dear All,</p><p style="margin: 0px; padding: 0px;" mce_style="margin: 0px; padding: 0px;"> I have a doubt I am not able to fix by myself. I am doing an "scf" spin polarized calculation for both the ZB and the WZ form of ZnO, at HSE level, with the development version of PWscf. For what regard the pseudo of Zn, I am using one with 2e only in the valence (PRB 81, 125207 (2010)). What I notice is that the DeltaE between the two systems is huge (4.4 Ry), with a macroscopic error on the bandgap of ZB-ZnO. Can this only be ascribed to the small number of electrons in the valence? I here report the two input files for both systems. I moreover think that the downsampling of the Fock exchange operator (q) can not be so much influent. But maybe I am wrong</p><p style="margin: 0px; padding: 0px;" mce_style="margin: 0px; padding: 0px;"><br></p><p style="margin: 0px; padding: 0px;" mce_style="margin: 0px; padding: 0px;">Thanks in advance for any hint you may give me.<br></p><p style="margin: 0px; padding: 0px;" mce_style="margin: 0px; padding: 0px;"><br></p><p style="margin: 0px; padding: 0px;" mce_style="margin: 0px; padding: 0px;">Best, <br></p><p style="margin: 0px; padding: 0px;" mce_style="margin: 0px; padding: 0px;">Giacomo</p><p style="margin: 0px; padding: 0px;" mce_style="margin: 0px; padding: 0px;">---------------------------------<br></p><p style="margin: 0px; padding: 0px;" mce_style="margin: 0px; padding: 0px;"><br></p><p style="margin: 0px; padding: 0px;" mce_style="margin: 0px; padding: 0px;"><br></p><p style="margin: 0px; padding: 0px;" mce_style="margin: 0px; padding: 0px;"><br></p><p style="margin: 0px; padding: 0px;" mce_style="margin: 0px; padding: 0px;"><br></p><p style="margin: 0px; padding: 0px;" mce_style="margin: 0px; padding: 0px;"><br></p><p style="margin: 0px; padding: 0px;" mce_style="margin: 0px; padding: 0px;">Giacomo Giorgi<br></p><p style="margin: 0px; padding: 0px;" mce_style="margin: 0px; padding: 0px;"><br></p><p style="margin: 0px; padding: 0px;" mce_style="margin: 0px; padding: 0px;"><tt>Department of Chemical System Engineering,<br>
School of Engineering, The University of Tokyo<br>
7-3-1 Hongo, Bunkyo-ku, Tokyo 113-8656, Japan<br>
E-mail: <a target="_self">giacomo@tcl.t.u-tokyo.ac.jp</a></tt></p><p style="margin: 0px; padding: 0px;" mce_style="margin: 0px; padding: 0px;"><br></p><p style="margin: 0px; padding: 0px;" mce_style="margin: 0px; padding: 0px;"><br></p><p style="margin: 0px; padding: 0px;" mce_style="margin: 0px; padding: 0px;"><br></p><p style="margin: 0px; padding: 0px;" mce_style="margin: 0px; padding: 0px;"><br></p><p style="margin: 0px; padding: 0px;" mce_style="margin: 0px; padding: 0px;"><br></p><p style="margin: 0px; padding: 0px;" mce_style="margin: 0px; padding: 0px;"><br></p><p style="margin: 0px; padding: 0px;" mce_style="margin: 0px; padding: 0px;"><br></p><p style="margin: 0px; padding: 0px;" mce_style="margin: 0px; padding: 0px;"><br></p><p style="margin: 0px; padding: 0px;" mce_style="margin: 0px; padding: 0px;"><br></p><p style="margin: 0px; padding: 0px;" mce_style="margin: 0px; padding: 0px;">1) WZ-ZnO</p><p style="margin: 0px; padding: 0px;" mce_style="margin: 0px; padding: 0px;"><br></p><p style="margin: 0px; padding: 0px;" mce_style="margin: 0px; padding: 0px;">&SYSTEM<br> ibrav= 4, celldm(1) =6.140475627, celldm(3)=1.6014648858250753985351141749246, nat= 4, ntyp= 2,<br> ecutwfc =80.0, nbnd = 20,force_symmorphic=.true., nspin=2, tot_magnetization=0,<br> input_dft='hse', nqx1 = 2, nqx2 = 2, nqx3 = 1,<br> x_gamma_extrapolation = .TRUE.<br>/<br>&ELECTRONS<br> diagonalization='david'<br> mixing_mode = 'plain'<br> mixing_beta = 0.7<br> conv_thr = 1.0d-8<br>/<br>ATOMIC_SPECIES<br> Zn 65.409 Zn_pz-2_nlcc.cpi.UPF<br> O 15.999 O.pz-mt.UPF<br>ATOMIC_POSITIONS (crystal)<br> Zn 0.33333 0.66667 0.00000<br> Zn 0.66667 0.33333 0.50000<br> O 0.33333 0.66667 0.38210<br> O 0.66667 0.33333 0.88210<br>K_POINTS (automatic)<br>8 8 6 0 0 0<br></p><p style="margin: 0px; padding: 0px;" mce_style="margin: 0px; padding: 0px;"><br></p><p style="margin: 0px; padding: 0px;" mce_style="margin: 0px; padding: 0px;">for this one I get (and it seems a reasonable result, at least from the gap)<br></p><p style="margin: 0px; padding: 0px;" mce_style="margin: 0px; padding: 0px;"><br></p><p style="margin: 0px; padding: 0px;" mce_style="margin: 0px; padding: 0px;"> highest occupied, lowest unoccupied level (ev): 1.4320 5.2338</p><p style="margin: 0px; padding: 0px;" mce_style="margin: 0px; padding: 0px;"> total energy = -110.48522966 Ry<br></p><p style="margin: 0px; padding: 0px;" mce_style="margin: 0px; padding: 0px;"><br></p><p style="margin: 0px; padding: 0px;" mce_style="margin: 0px; padding: 0px;"><br></p><p style="margin: 0px; padding: 0px;" mce_style="margin: 0px; padding: 0px;"><br></p><p style="margin: 0px; padding: 0px;" mce_style="margin: 0px; padding: 0px;"><br></p><p style="margin: 0px; padding: 0px;" mce_style="margin: 0px; padding: 0px;"><br></p><p style="margin: 0px; padding: 0px;" mce_style="margin: 0px; padding: 0px;"><br></p><p style="margin: 0px; padding: 0px;" mce_style="margin: 0px; padding: 0px;"><br></p><p style="margin: 0px; padding: 0px;" mce_style="margin: 0px; padding: 0px;"><br></p><p style="margin: 0px; padding: 0px;" mce_style="margin: 0px; padding: 0px;">2) ZB-ZnO</p><p style="margin: 0px; padding: 0px;" mce_style="margin: 0px; padding: 0px;"><br></p><p style="margin: 0px; padding: 0px;" mce_style="margin: 0px; padding: 0px;">&system<br> ibrav= 2, celldm(1) =8.447075169<br> nat= 8, ntyp= 2, force_symmorphic=.true., nspin=2, tot_magnetization=0,<br> ecutwfc =80., nbnd = 50,<br> input_dft='hse', nqx1=2, nqx2=2, nqx3=2,<br> x_gamma_extrapolation = .TRUE.<br>/<br>&electrons<br> diagonalization='david'<br> mixing_mode = 'plain'<br> mixing_beta = 0.7<br> conv_thr = 1.0d-8<br>/<br>ATOMIC_SPECIES<br> Zn 65.38 Zn_pz-2_nlcc.cpi.UPF<br> O 15.9994 O.pz-mt.UPF<br>ATOMIC_POSITIONS (crystal)<br> Zn 0.00000 0.00000 0.00000<br> Zn 0.00000 0.50000 0.50000<br> Zn 0.50000 0.00000 0.50000<br> Zn 0.50000 0.50000 0.00000<br> O 0.25000 0.25000 0.25000<br> O 0.25000 0.75000 0.75000<br> O 0.75000 0.75000 0.25000<br> O 0.75000 0.25000 0.75000<br>K_POINTS automatic<br>8 8 8 0 0 0<br></p><p style="margin: 0px; padding: 0px;" mce_style="margin: 0px; padding: 0px;"><br></p><p style="margin: 0px; padding: 0px;" mce_style="margin: 0px; padding: 0px;">For this one I get (!!)<br></p><p style="margin: 0px; padding: 0px;" mce_style="margin: 0px; padding: 0px;"><br></p><p style="margin: 0px; padding: 0px;" mce_style="margin: 0px; padding: 0px;">highest occupied, lowest unoccupied level (ev): 36.4937 31.9513 (!!!) <br></p><p style="margin: 0px; padding: 0px;" mce_style="margin: 0px; padding: 0px;">total energy = -200.32746186 Ry<br></p><div style="position: absolute; left: -10000px; top: 0px; width: 1px; height: 1px; overflow: hidden;" id="_mcePaste"> </div>