<div dir="ltr" id="yMail_cursorElementTracker_1641287215579">Dear Nishidate,<br clear="none"><br clear="none">Thank you for your answer.<br clear="none"><br clear="none">I don't know if I understood... Could you explain me better what should I do, please?<br clear="none"><br clear="none">I'm sending a figure of my DOS plotted using the HOMO and the EFermi as the zero, respectively, to make my doubt clearer. <br clear="none"><br clear="none">I learned that I should use a "large" number of k-points to converge the system. However, I only found the first figure when I used the 10x10x10 k-points in the NSCF calculation, and 2x2x2, 4x4x4 or 6x6x6 in the SCF. When I increased or decreased the k-points in the NSCF, it showed me the second figure (that is different from the reference one).</div><div dir="ltr" id="yMail_cursorElementTracker_1641287642135"><br></div><div dir="ltr" id="yMail_cursorElementTracker_1641287642280">I tried to use the smearing occupations. In this test, the EFermi in SCF and NSCF calculations were the same (2.74 eV), except when I used 10x10x10 k-points in the NSCF. In this case, the EFermi was 2.74 in SCF and 1.04 eV in the NSCF (it is the same value that I found previously as the HOMO in SCF calculation). </div><div dir="ltr" id="yMail_cursorElementTracker_1641288215288"><br></div><div dir="ltr" id="yMail_cursorElementTracker_1641288215451">My main question is: Why is this 10x10x10 a "magic number"? <br clear="none"><br clear="none">Best wishes.</div><div dir="ltr" id="yMail_cursorElementTracker_1641288144811">I wish you a happy new year.</div><div dir="ltr" id="yMail_cursorElementTracker_1641288191468"><br></div><div dir="ltr" id="yMail_cursorElementTracker_1641287780250">José Xavier<br clear="none"><br clear="none"><br clear="none">>  <br clear="none">>  <br clear="none">> Em qua., 29 29e dez. 29e 2021 às 18:22, Kazume NISHIDATE<br clear="none">> <<a shape="rect" ymailto="mailto:nisidate@iwate-u.ac.jp" href="mailto:nisidate@iwate-u.ac.jp">nisidate@iwate-u.ac.jp</a>> escreveu:<div class="yqt0949207044" id="yqtfd22236"><br clear="none"><br clear="none">> <br clear="none">> <br clear="none">>  <br clear="none">> <br clear="none"><br clear="none"><br clear="none"><br clear="none">Dear Xavier<br clear="none"><br clear="none">> When I used the k-points 4x4x4 (SCF) and 8x8x8 (NSCF), the HOMO and<br clear="none">> EFermi were different, and using EFermi to subtract the values in the<br clear="none"><br clear="none">If the k-point mesh you selected covers the typical energy peaks of<br clear="none">the valence band maximum (VBM or the HOMO in your description) depends<br clear="none">on the shape of Brillouin zone (BZ) of your system.<br clear="none"><br clear="none">You should try to plot the band structure and assess the VBM and find<br clear="none">out its k-point location. If the k-point mesh you selected does not<br clear="none">cover the k-point, then the VBM peak will not contribute to your DOS<br clear="none"><br clear="none">plot.<br clear="none"><br clear="none"><br clear="none"><br clear="none"><br clear="none"><br clear="none">敬具　西館<br clear="none"><br clear="none">best regards<br clear="none"><br clear="none"><br clear="none"><br clear="none"><a shape="rect" ymailto="mailto:nisidate@iwate-u.ac.jp" href="mailto:nisidate@iwate-u.ac.jp">nisidate@iwate-u.ac.jp</a><br clear="none"><br clear="none"><a shape="rect" ymailto="mailto:kazume.nishidate@gmail.com" href="mailto:kazume.nishidate@gmail.com">kazume.nishidate@gmail.com</a><br clear="none"></div></div>