<div dir="ltr">Dear all,<div><br></div><div>I have been attempting to calculate the band structure (and consequently the band gap) of FeO, using HSE XC and ONCV pseudopotentials, in QE v. 6.2.1 and 6.4.1 (plus Wannier90 to interpolate the band structure). However, the underlying SCF calculations have proven quite volatile, frequently but inconsistently running into "dexx < 0" errors for particular values of `nk` and `nq` – which makes a proper convergence of the results more difficult.</div><div><br></div><div>I've tried a fair number of variations to circumvent the issue, with minimal success (among them: changing the energy cutoff, convergence threshold, smearing type, or EXX adaptive threshold – both varying initial threshold and removing it altogether.) Occasionally, a different `adaptive_thr` will avoid dexx < 0, but more often it makes little difference (and sometimes does worse).</div><div><br></div><div>In addition, I find that changing from smearing to tetrahedral occupations leads to noticeable qualitative differences in energy levels and resulting band structure. While the results with tetrahedral occupation seem more reasonable, the SCF calculations then also tend to converge slower and/or be more prone to dexx < 0.</div><div><br></div><div>Is this at all expected behavior? Any suggestions to improve these results would be greatly appreciated.</div><div><br></div><div>Sample Input (fails with dexx<0 for either occupation method):</div><div> &control<br></div><div>    calculation = 'scf'<br>    restart_mode='from_scratch',<br>    prefix='FeO',<br>    pseudo_dir = '[...]',<br>    outdir = '[...]',<br>    verbosity = 'high'<br> /<br> &system<br>    ibrav=  0,<br>! lattice constant taken from GGA optimization<br>    A = 4.3100,<br>    nat=  4, ntyp= 3, nbnd=30,<br>    nspin = 2,<br>    starting_magnetization(1)= 0.0,<br>    starting_magnetization(2)= 0.5,<br>    starting_magnetization(3)=-0.5,<br>!    occupations='tetrahedra',   ! alt.</div><div>    occupations='smearing', smearing='mv', degauss=0.01,<br>    ecutwfc = 90.0,<br>    ecutrho = 360.0,<br>    input_dft='hse',<br>    exxdiv_treatment='gygi-baldereschi',<br>    x_gamma_extrapolation = .true.,<br>    nqx1=2, nqx2=2, nqx3=2<br> /<br> &electrons<br>    mixing_mode = 'plain'<br>    mixing_beta = 0.3<br>    conv_thr =  1.0d-6<br>    mixing_fixed_ns = 0,<br>    adaptive_thr = .true.<br> /<br></div><div>! Rhombohedral 4-atom AFM cell - cf. PW example08<br>CELL_PARAMETERS {alat}<br>0.50 0.50 1.00<br>0.50 1.00 0.50<br>1.00 0.50 0.50<br>! SG15 ONCV pseudopotentials (scalar-rel)<br>ATOMIC_SPECIES<br> O1   1.  O_ONCV_PBE-1.0.upf<br> Fe1  1.  Fe_ONCV_PBE-1.0.upf<br> Fe2  1.  Fe_ONCV_PBE-1.0.upf<br>ATOMIC_POSITIONS {crystal}<br> O1  0.25 0.25 0.25<br> O1  0.75 0.75 0.75<br> Fe1 0.0  0.0  0.0<br> Fe2 0.5  0.5  0.5<br>K_POINTS {automatic}<br>6 6 6 0 0 0<br></div><div><br></div><div>[End of sample input]<br></div><div><br></div><div>Peter Stoeckl</div><div>Department of Physics</div><div>University of Minnesota</div></div>