Dear Nishidate,<div id="yMail_cursorElementTracker_1640782920739">Thank you for your answer.</div><div id="yMail_cursorElementTracker_1640783669147"><br></div><div id="yMail_cursorElementTracker_1640783669329">If I understood correctly, you are telling me to plot the DOS using the EFermi to subtract the energy (first column of the DOS file), and set it as the zero. Is it?<br><div id="yMail_cursorElementTracker_1640782884321"><br></div><div id="yMail_cursorElementTracker_1640782884671">When I used the k-points 4x4x4 (SCF) and 8x8x8 (NSCF), the HOMO and EFermi were different, and using EFermi to subtract the values in the DOS file to set it equal to the zero changed the peaks, with the last peak of the Valence band ending close to -5 eV. On the other hand, when I used the HOMO value (1.04 eV) to "set the zero", this peak ended in the zero point. The same occur when I use the k-points 6x6x6 (SCF) and 12x12x12 (NSCF).</div><div id="yMail_cursorElementTracker_1640783404849"><br></div><div id="yMail_cursorElementTracker_1640783405007">However, when I use the k-points 4x4x4 or 6x6x6 (SCF) and 10x10x10 (NSCF), the EFermi is equal to HOMO (1.04 eV) and the plot of the DOS is "correct", with the last peak of Valence band ending in the zero point.</div><div id="yMail_cursorElementTracker_1640783865292"><br></div><div id="yMail_cursorElementTracker_1640783865432">I don't understand why only these two values of kpoints in SCF and 10x10x10 in NSCF gave me "correct" peaks, at the same time, why all the other combinations I've tried gave me a "strange" arrangement in the peaks. </div><div id="yMail_cursorElementTracker_1640784383934"><br></div><div id="yMail_cursorElementTracker_1640784384047">Best wishes.</div><div id="yMail_cursorElementTracker_1640784388705">And thank you again for answer my question</div><div id="yMail_cursorElementTracker_1640784409705"><br></div><div id="yMail_cursorElementTracker_1640784409836">José Xavier</div><div id="yMail_cursorElementTracker_1640783313599"><br><div id="ymail_android_signature"><a id="ymail_android_signature_link" href="https://go.onelink.me/107872968?pid=InProduct&c=Global_Internal_YGrowth_AndroidEmailSig__AndroidUsers&af_wl=ym&af_sub1=Internal&af_sub2=Global_YGrowth&af_sub3=EmailSignature">Enviado do Yahoo Mail no Android</a></div> <br> <blockquote style="margin: 0 0 20px 0;"> <div style="font-family:Roboto, sans-serif; color:#6D00F6;"> <div>Em qua., 29 29e dez. 29e 2021 às 2:13, Kazume NISHIDATE</div><div><nisidate@iwate-u.ac.jp> escreveu:</div> </div> <div style="padding: 10px 0 0 20px; margin: 10px 0 0 0; border-left: 1px solid #6D00F6;"> <div id="yiv6039729971"><title></title><div>
<div>
<div>Dear Xavier,<br clear="none">
<br clear="none">
> performed a calculation using 6x6x6 k-point in SCF (fixed<br clear="none">
> occupations) and 10x10x10 in NSCF (tetrahedra occupations). The E<br clear="none">
> fermi was equal to the HOMO value (1.04 eV). However, when I increase<br clear="none">
> (e.g. 8x8x8 and 12x12x12,<br clear="none">
<br clear="none">
It is normal that the Fermi level changes for the calculation to<br clear="none">
calculation with different set of k-points or cutoff energies.<br clear="none">
Because the energy standard (energy zero) can not be defined solely in<br clear="none">
the DFT calculation.<br clear="none">
<br clear="none">
> respectively) or decrease (e.g. 4x4x4 and 8x8x8) the number of<br clear="none">
> k-points in each calculation, the E fermi (6.03 eV) and HOMO (1.04<br clear="none">
> eV) are quite different.<br clear="none">
<br clear="none">
You should try to plot the DOS figures with shifting the energy<br clear="none">
scale. The easiest way to do this is to set the E fermi equal to 0.<br clear="none">
<br clear="none">
E(4x4x4) = E(4x4x4) - 1.04<br clear="none">
E(6x6x6) = E(6x6x6) - 6.03<br clear="none">
<br clear="none">
The resultant DOS figures of the two will be identical.</div>
</div>
<div><br clear="none">
<div class="yiv6039729971matchFont">
<div style="color:rgb(0, 0, 0);orphans:2;widows:2;" class="yiv6039729971">敬具 西館<br clear="none"></div>
<div style="color:rgb(0, 0, 0);orphans:2;widows:2;" class="yiv6039729971"><span style="">best regards</span><br clear="none"></div>
<div style="color:rgb(0, 0, 0);orphans:2;widows:2;" class="yiv6039729971"><br clear="none"></div>
<div style="color:rgb(0, 0, 0);orphans:2;widows:2;" class="yiv6039729971"><a rel="nofollow noopener noreferrer" shape="rect" ymailto="mailto:nisidate@iwate-u.ac.jp" target="_blank" href="mailto:nisidate@iwate-u.ac.jp" class="yiv6039729971">nisidate@iwate-u.ac.jp</a></div>
<div style="color:rgb(0, 0, 0);orphans:2;widows:2;" class="yiv6039729971"><span style="" class="yiv6039729971"><a rel="nofollow noopener noreferrer" shape="rect" ymailto="mailto:kazume.nishidate@gmail.com" target="_blank" href="mailto:kazume.nishidate@gmail.com" class="yiv6039729971">kazume.nishidate@gmail.com</a></span></div>
</div>
</div>
<div>2021年12月29日 5:17 +0900、José Xavier via users <users@lists.quantum-espresso.org> のメール:<div id="yiv6039729971yqtfd18897" class="yiv6039729971yqt5350300580"><br clear="none">
</div><blockquote type="cite" style="border-left-color:grey;border-left-width:thin;border-left-style:solid;margin:5px 5px;padding-left:10px;"><div id="yiv6039729971yqtfd34468" class="yiv6039729971yqt5350300580">Dear QE users,<br clear="none">
<br clear="none">
I performed a calculation using 6x6x6 k-point in SCF (fixed occupations) and 10x10x10 in NSCF (tetrahedra occupations). The E fermi was equal to the HOMO value (1.04 eV). However, when I increase (e.g. 8x8x8 and 12x12x12, respectively) or decrease (e.g. 4x4x4 and 8x8x8) the number of k-points in each calculation, the E fermi (6.03 eV) and HOMO (1.04 eV) are quite different.<br clear="none">
<br clear="none">
Could someone explain to me why is this happening?<br clear="none">
<br clear="none">
<br clear="none">
Sincerely,<br clear="none">
José Xavier<br clear="none">
Department of Biophysics and Pharmacology<br clear="none">
Federal University of Rio Grande do Norte<br clear="none">
Natal, Brazil</div><br clear="none">
_______________________________________________<br clear="none">
Quantum ESPRESSO is supported by MaX (www.max-centre.eu)<br clear="none">
users mailing list users@lists.quantum-espresso.org<br clear="none">
https://lists.quantum-espresso.org/mailman/listinfo/users</blockquote><div id="yiv6039729971yqtfd01272" class="yiv6039729971yqt5350300580">
</div></div><div id="yiv6039729971yqtfd21627" class="yiv6039729971yqt5350300580">
</div></div></div> </div> </blockquote></div></div>