<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;">Hi there, <br>I have a problem related to DOS at fermi level. Let me start with my scf file first<br>-------------<br> &system<br>    ibrav=0, celldm(1)=0 , nat=22, ntyp=4,<br>      ecutwfc =37.0, occupations='smearing', smearing='gauss', degauss=0.01, nelec=156<br>    nspin= 1,<br> /<br> &electrons<br>  electron_maxstep= 300<br>    mixing_mode = 'plain'<br>    mixing_beta = 0.2<br>    conv_thr =  1.0d-6<br>    diagonalization= 'david'<br> /<br>--------------<br>My pseudopotentials are norm conversing, which are generated by  ld1.x and their transferability is OK.<br>here is the output<br>----------<br>the Fermi energy is     3.4391 ev<br><br>!    total
 energy              =   -1213.21835062 Ry<br>     Harris-Foulkes estimate   =   -1213.21835085 Ry<br>     estimated scf accuracy    <       0.00000034 Ry<br><br>     The total energy is the sum of the following terms:<br><br>     one-electron contribution =    -711.80653324 Ry<br>     hartree contribution      =     352.57845405 Ry<br>     xc contribution           =    -222.78424774 Ry<br>     ewald contribution        =    -631.18916691 Ry<br>     smearing contrib. (-TS)  
 =      -0.01685678 Ry<br><br>     convergence has been achieved in  32 iterations<br>----------------------<br>and Dos<br>3.000  0.2551E+02  0.1178E+03<br>  3.100  0.2931E+02  0.1207E+03<br>  3.200  0.3391E+02  0.1241E+03<br>  3.300  0.3375E+02  0.1275E+03<br>  3.400  0.2724E+02  0.1302E+03<br>  3.500  0.2299E+02  0.1325E+03<br>  3.600  0.2335E+02  0.1348E+03<br>  3.700  0.2043E+02  0.1369E+03<br>  3.800  0.1297E+02  0.1381E+03<br>  3.900  0.9369E+01  0.1391E+03<br>  4.000  0.8969E+01  0.1400E+03<br><br>--------------<br>My question are concerned to:<br>+ DOS is still not zero at fermi level E=3.4391<br>+ Integrated DOS is not equal number of electron (156) , even I try DeltaE=0.001<br>-------------<br>I have thought about the question<br>+First,
 My compound is insulator, an oxide with optical bandgap of 4 eV. But in scf input, I use occupations='smearing'. Shoud I change it to occupations='fixed' and increase nbnd>156/2 and delete smearing variable.<br>------------<br>Is there anyone had this problem before and please help me to fix that.<br><br>Sau Nguyen<br>U of Houston<br><br></td></tr></table><br>