<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>&nbsp;&amp;system<br>&nbsp;&nbsp;&nbsp; ibrav=0, celldm(1)=0 , nat=22, ntyp=4,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ecutwfc =37.0, occupations='smearing', smearing='gauss', degauss=0.01, nelec=156<br>&nbsp;&nbsp;&nbsp; nspin= 1,<br>&nbsp;/<br>&nbsp;&amp;electrons<br>&nbsp; electron_maxstep= 300<br>&nbsp;&nbsp;&nbsp; mixing_mode = 'plain'<br>&nbsp;&nbsp;&nbsp; mixing_beta = 0.2<br>&nbsp;&nbsp;&nbsp; conv_thr =&nbsp; 1.0d-6<br>&nbsp;&nbsp;&nbsp; diagonalization= 'david'<br>&nbsp;/<br>--------------<br>My pseudopotentials are norm conversing, which are generated by&nbsp; ld1.x and their transferability is OK.<br>here is the output<br>----------<br>the Fermi energy is&nbsp;&nbsp;&nbsp;&nbsp; 3.4391 ev<br><br>!&nbsp;&nbsp;&nbsp; total
 energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -1213.21835062 Ry<br>&nbsp;&nbsp;&nbsp;&nbsp; Harris-Foulkes estimate&nbsp;&nbsp; =&nbsp;&nbsp; -1213.21835085 Ry<br>&nbsp;&nbsp;&nbsp;&nbsp; estimated scf accuracy&nbsp;&nbsp;&nbsp; &lt;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00000034 Ry<br><br>&nbsp;&nbsp;&nbsp;&nbsp; The total energy is the sum of the following terms:<br><br>&nbsp;&nbsp;&nbsp;&nbsp; one-electron contribution =&nbsp;&nbsp;&nbsp; -711.80653324 Ry<br>&nbsp;&nbsp;&nbsp;&nbsp; hartree contribution&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp;&nbsp; 352.57845405 Ry<br>&nbsp;&nbsp;&nbsp;&nbsp; xc contribution&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp; -222.78424774 Ry<br>&nbsp;&nbsp;&nbsp;&nbsp; ewald contribution&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp; -631.18916691 Ry<br>&nbsp;&nbsp;&nbsp;&nbsp; smearing contrib. (-TS)&nbsp;&nbsp;
 =&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -0.01685678 Ry<br><br>&nbsp;&nbsp;&nbsp;&nbsp; convergence has been achieved in&nbsp; 32 iterations<br>----------------------<br>and Dos<br>3.000&nbsp; 0.2551E+02&nbsp; 0.1178E+03<br>&nbsp; 3.100&nbsp; 0.2931E+02&nbsp; 0.1207E+03<br>&nbsp; 3.200&nbsp; 0.3391E+02&nbsp; 0.1241E+03<br>&nbsp; 3.300&nbsp; 0.3375E+02&nbsp; 0.1275E+03<br>&nbsp; 3.400&nbsp; 0.2724E+02&nbsp; 0.1302E+03<br>&nbsp; 3.500&nbsp; 0.2299E+02&nbsp; 0.1325E+03<br>&nbsp; 3.600&nbsp; 0.2335E+02&nbsp; 0.1348E+03<br>&nbsp; 3.700&nbsp; 0.2043E+02&nbsp; 0.1369E+03<br>&nbsp; 3.800&nbsp; 0.1297E+02&nbsp; 0.1381E+03<br>&nbsp; 3.900&nbsp; 0.9369E+01&nbsp; 0.1391E+03<br>&nbsp; 4.000&nbsp; 0.8969E+01&nbsp; 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&gt;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>