<div dir="ltr"><div style="font-family: arial, sans-serif; font-size: 12.8px;">I have monolayers of various materials which have been simulated multiple times and with different software.</div><div style="font-family: arial, sans-serif; font-size: 12.8px;">There is a vacuum in z direction 15 Ang. PBE, GGA is used.</div><div style="font-family: arial, sans-serif; font-size: 12.8px;">Let's take the below scenario where the different levels are given in eV .</div><div style="font-family: arial, sans-serif; font-size: 12.8px;"><br></div><span style="font-family: arial, sans-serif; font-size: 12.8px;">VBM     = -3.88        Valence band</span><div style="font-family: arial, sans-serif; font-size: 12.8px;">CBM     = -2.255     conduction band</div><div style="font-family: arial, sans-serif; font-size: 12.8px;">Ef         = -1.44       Fermi energy from cp2k output</div><div style="font-family: arial, sans-serif; font-size: 12.8px;">Ef         = -2.66       (VBM + (CBM-VBM)/2) mid-gap</div><div style="font-family: arial, sans-serif; font-size: 12.8px;"><br></div><div style="font-family: arial, sans-serif; font-size: 12.8px;">The fermi level is in the bandgap.</div><div style="font-family: arial, sans-serif; font-size: 12.8px;"><br></div><div style="font-family: arial, sans-serif; font-size: 12.8px;">The error in the two fermi levels above :</div><div style="font-family: arial, sans-serif; font-size: 12.8px;">delta Ef = 0.40</div><div style="font-family: arial, sans-serif; font-size: 12.8px;"><br></div><div style="font-family: arial, sans-serif; font-size: 12.8px;">All my systems are monolayers but I have simulated these monolayers earlier using supercells. I know just using Gamma points doesn't give</div><div style="font-family: arial, sans-serif; font-size: 12.8px;">a good resolution but still there shouldn't be 0.4 eV difference. Would you suggest just taking an average of VBM and CBM in that case.</div><div style="font-family: arial, sans-serif; font-size: 12.8px;">VBM and CBM are identified by filling total electrons with 2 per level. so basically nelectrons/2 th energy level is my valence band and the next</div><div style="font-family: arial, sans-serif; font-size: 12.8px;">one with 0 occupation my conduction band.</div><div style="font-family: arial, sans-serif; font-size: 12.8px;"><br></div><div style="font-family: arial, sans-serif; font-size: 12.8px;"><br></div><div style="font-family: arial, sans-serif; font-size: 12.8px;">Thanks for any suggestions.</div><br>On Wednesday, 2 August 2017 10:33:46 UTC+2, Matt W  wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Hi,<div><br></div><div>can you define 'slight offset'? The Fermi level is determined numerically, so it is bound to not be precisely at the centre of the gap.</div><div><br></div><div>(Curiosity - should the Fermi level be precisely mid gap for a real material, or is the result just for simple 1D bands and no degeneracy?)</div><div><br></div><div>Matt</div><div><br></div><div><br><br>On Tuesday, August 1, 2017 at 5:48:45 PM UTC+1, ashish dabral wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi,<div><br></div><div>I have calculated PDOS for semiconducting systems. While plotting dos I realized that for the semiconductor, </div><div>the fermi energy is not lying in the midgap. No spin polarization is used. I just take the fermi energy and </div><div>check if the mid of CB min. and VB maximum give me the fermi energy(mid-gap). There is a slight offset towards</div><div>conduction band. The system is definitely semiconducting. I can calculate the fermi energy by taking mid value</div><div>from CB and VB but why is fermi energy from cp2k  calculation is not matching? Only Gamma point is used as</div><div>the systems are large (100-200 atoms)</div><div><br></div><div>Any pointers?</div><div><br></div><div>Thanks</div></div></blockquote></div></div></blockquote></div>