<div dir="ltr"><div><div>Dear all,<br></div>           I was calculating band gap with hybrid functionals for a semiconductor. If I don't use smearing (i.e. occupations = 'fixed') the dos  seems ok but whenever I use smearing and plot dos, at some energies I get negative density of states ! I can't understand why this is happening. Can/should one use smearing for hybrid functionals? Following is the important part of the input:<br>
<br> &system<br>    ibrav = 4,<br>    celldm(1) = 6.474424, celldm(3) = 5.837554613 ,<br>    nat = 10, ntyp = 3,<br>    ecutwfc = 40, <br>    london = .true.,<br>    nbnd = 25,<br>    occupations = 'smearing', smearing = 'mv', degauss = 0.001,<br>
    input_dft = 'hse',<br>    exxdiv_treatment = 'gygi-baldereschi',<br>    x_gamma_extrapolation = .true., <br>    nqx1 = 4, nqx2 = 4, nqx3 = 1<br><br>...<br><br></div>    Another problem is that whenever I use smearing the output doesn't give HOMO-LUMO as it considers it as metal and prints only the fermi energy. So if I'm studying semiconductor to metal transition there's no way one can check the band gap change. <br>
<div>                 <br></div><div>    I'd be greatful if anyone can shine some light on above problems.<br><br></div><div>Thnaks,<br></div><div>Rajdeep Banerjee<br></div><div>Ph. D. student<br></div><div>JNCASR<br>
</div><div>Bangalore, India <br></div></div>