<div dir="ltr"><div>The two figures are exactly the same, only the origin of the x axis is shifted.<br></div><div><br></div><div>Paolo<br></div><div dir="ltr"><br></div><div dir="ltr">On Tue, Jan 4, 2022 at 10:29 AM José Xavier via users <<a href="mailto:users@lists.quantum-espresso.org">users@lists.quantum-espresso.org</a>> wrote:<br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr" id="gmail-m_5301897612338113415yMail_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="gmail-m_5301897612338113415yMail_cursorElementTracker_1641287642135"><br></div><div dir="ltr" id="gmail-m_5301897612338113415yMail_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="gmail-m_5301897612338113415yMail_cursorElementTracker_1641288215288"><br></div><div dir="ltr" id="gmail-m_5301897612338113415yMail_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="gmail-m_5301897612338113415yMail_cursorElementTracker_1641288144811">I wish you a happy new year.</div><div dir="ltr" id="gmail-m_5301897612338113415yMail_cursorElementTracker_1641288191468"><br></div><div dir="ltr" id="gmail-m_5301897612338113415yMail_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" href="mailto:nisidate@iwate-u.ac.jp" target="_blank">nisidate@iwate-u.ac.jp</a>> escreveu:<div id="gmail-m_5301897612338113415yqtfd22236"><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" href="mailto:nisidate@iwate-u.ac.jp" target="_blank">nisidate@iwate-u.ac.jp</a><br clear="none"><br clear="none"><a shape="rect" href="mailto:kazume.nishidate@gmail.com" target="_blank">kazume.nishidate@gmail.com</a><br clear="none"></div></div>_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" rel="noreferrer" target="_blank">www.max-centre.eu</a>)<br>
users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote></div><br clear="all"><br>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>Univ. Udine, via delle Scienze 206, 33100 Udine, Italy<br>Phone +39-0432-558216, fax +39-0432-558222<br><br></div></div></div></div></div></div>