<div dir="ltr"><div>Hi,</div><div><br></div><div>there is not a problem. </div><div><br></div><div>CP2K reports the fermi level at the same energy as the HOMO. Your smearing-pdos script broadens the energy of the highest level so it looks like there is DOS slightly above the reported fermi level.  But this is just an artifact of the broadening applied for plotting - e'd guess this might be set by the 'ww = 0.08' in your script - try reducing this and see.</div><div><br></div><div>Matt</div><div><br></div><div><br></div><br><br>On Friday, August 4, 2017 at 3:26:53 PM UTC+1, zhj...@gmail.com wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr"><blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><font size="4">Hi, Matt,</font></div></div></blockquote></div></div></blockquote><div><font size="4"><br></font></div><div><font size="4">Thank you for your reply.</font></div><div><font size="4">Now the Fermi level is not in the band gap, can you tell me how can I solve this problem? </font></div><div><font size="4"><br></font></div></div></blockquote></div>