<html><head><meta http-equiv="Content-Type" content="text/html; charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><span style="caret-color: rgb(255, 255, 255); font-family: UICTFontTextStyleBody; font-size: 17px;" class="">Dear Salman, to give you an exhaustive response regarding the behaviour of your DOS, outputs and/or plots would be necessary. However, the only thing I would say you can take for sure is that the values you get from the code are in eV. </span><div style="caret-color: rgb(255, 255, 255); font-family: UICTFontTextStyleBody; font-size: 17px;" class=""><br class=""></div><div style="caret-color: rgb(255, 255, 255); font-family: UICTFontTextStyleBody; font-size: 17px;" class="">The differences of the values you observe just by changing the cell size are due to a change of the physics and chemistry of your system with the cell side. You have made a virtual experiment in which you expand the side of the cell in all three direction by as much as 100%. It would have been strange if you had get very similar values.</div><div style="caret-color: rgb(255, 255, 255); font-family: UICTFontTextStyleBody; font-size: 17px;" class=""><br class=""></div><div style="caret-color: rgb(255, 255, 255); font-family: UICTFontTextStyleBody; font-size: 17px;" class="">Just as a check, verify that the three integrated DOS are the same.</div><div style="caret-color: rgb(255, 255, 255); font-family: UICTFontTextStyleBody; font-size: 17px;" class=""><br class=""></div><div style="caret-color: rgb(255, 255, 255); font-family: UICTFontTextStyleBody; font-size: 17px;" class="">Giovanni</div><div><br class=""><blockquote type="cite" class=""><div class="">On 14 Nov 2021, at 16:37, Salman Zarrini <<a href="mailto:salman.zarrini@gmail.com" class="">salman.zarrini@gmail.com</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div dir="ltr" class=""><div dir="ltr" class=""><div dir="ltr" class=""><div dir="ltr" class=""><div dir="ltr" class=""><div dir="ltr" class=""><div dir="ltr" class=""><div dir="ltr" class=""><div dir="ltr" class=""><div dir="ltr" class=""><div dir="ltr" class=""><div dir="ltr" class=""><div dir="ltr" class=""><div dir="ltr" class=""><div dir="ltr" class=""><div dir="ltr" class=""><div dir="ltr" class=""><div dir="ltr" class=""><div dir="ltr" class=""><div dir="ltr" class=""><div dir="ltr" class=""><div dir="ltr" class=""><div dir="ltr" class=""><div dir="ltr" class=""><div dir="ltr" class=""><div dir="ltr" class="">Thanks all for your responses.</div><div dir="ltr" class=""><br class=""><div class="">@Giovanni:</div><div class="">Actually, as an example I had considered the FCC Ni in three hypothetical systems with primitive cells with celldm(1) of A0 = 3.511 angstrom, A1 = 1.5*A0 and A2 = 2*A0, respectively. So, only the celldm increases from one system to another one. Please note that they are not super-cell.</div><div class="">And ran Quantum-Espresso to calculate the density of states for each of the three systems using similar convergence criterias. To avoid any complication regarding the spin, I hypothetically just used the spin--nonpolarized mode in my calculations.</div><div class=""><br class=""></div><div class="">The results for density of states at the Fermi level (<i class="">N</i>(<i class="">E</i>_F)) for each system is as follows:</div><div class="">A0 :> <i class="">N</i>(<i class="">E</i>_F) = 5 </div><div class="">A1 = 1.5*A0 :> <i class="">N</i>(<i class="">E</i>_F) = 33.78<br class=""></div><div class="">A2 = 2*A0:> <i class="">N</i>(<i class="">E</i>_F) = 159</div><div class=""><br class=""></div><div class="">And then the ratio of <i class="">N</i>(<i class="">E</i>_F) / V where V is the volume of primitive cell of each system is as follows:</div><div class="">A0 :> <i class="">N</i>(<i class="">E</i>_F) / V0 = 0.462<br class=""></div><div class="">A1 = 1.5A0 :> <i class="">N</i>(<i class="">E</i>_F) / V1 = 0.924</div><div class="">A2 = 2*A0:> <i class="">N</i>(<i class="">E</i>_F) / V2 = 1.836</div><div class=""><br class=""></div><div class="">So, comparing the density of states at the Fermi level, which indeed is a critical quantity in some concepts, shows that the <i class="">N</i>(<i class="">E</i>_F) / V change from one to another around two times from system A0 to system A1 and foure times from system A0 to system A2. Although the general <i class="">N</i>(<i class="">E</i>) / V are not as similar as I plotted them. </div><div class=""><br class=""></div><div class="">I understand the chemistry point of view as by enlarging the cell dimension, the overlap of <i class="">d</i> wave functions decreases so consequently the band width decreases and that increases the <i class="">N</i>(<i class="">E</i>_F) ultimately. However, I still have the problem in understanding the unit(s) used for density of states.</div><div class=""><br class=""></div><div class=""><div class=""><div dir="ltr" class=""><div dir="ltr" class=""><div class=""><div dir="ltr" class="">Regards,<br class="">
Salman Zarrini</div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div><br class=""><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sun, Nov 14, 2021 at 6:23 AM Giovanni Cantele <<a href="mailto:giovanni.cantele@spin.cnr.it" target="_blank" class="">giovanni.cantele@spin.cnr.it</a>> wrote:<br class=""></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div class="">Dear Salman,<div class=""><br class=""></div><div class="">if a quantity is extensive, so the larger the volume/mass/size the larger that quantity, it should NOT have Volume^-1 in its units.</div><div class=""><br class=""></div><div class="">Indeed, the DOS as calculated by Quantum-ESPRESSO is in eV^-1. Let us suppose that you calculate the ground state of a given crystal within its primitive cell and obtain a certain DOS. If you now compute the ground state of the SAME crystal but with twice the unit cell, the number of electrons doubles and as does the DOS. On the other hand, if you divide both DOSs by the respective unit cell volumes, you’ll get a quantity with Volume^-1 in its units that will be the exactly same for both calculations (provided the convergence of both calculations with respect to the used parameters is the same).</div><div class=""><br class=""></div><div class="">Giovanni<br class=""><div class=""><br class=""><blockquote type="cite" class=""><div class="">On 13 Nov 2021, at 21:06, Salman Zarrini <<a href="mailto:salman.zarrini@gmail.com" target="_blank" class="">salman.zarrini@gmail.com</a>> wrote:</div><br class=""><div class=""><div dir="ltr" class=""><div class=""><div dir="ltr" class=""><div dir="ltr" class=""><div class=""><div class="">So, that means Quantum-Espresso gives an extensive density of states, right? if so, then it should have a Volume^-1 in its unit.</div><div dir="ltr" class=""><br class=""></div><div dir="ltr" class="">Regards,<br class="">
Salman<br class=""></div></div></div></div></div><br class=""></div><br class=""><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, Nov 13, 2021 at 2:46 PM Stefano Baroni <<a href="mailto:baroni@sissa.it" target="_blank" class="">baroni@sissa.it</a>> wrote:<br class=""></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div dir="auto" class="">
it depends on the volume of the unit cell. once you divide by it, you get an intensive (volume-independent) quantity. sb<br class="">
<br class="">
<div dir="ltr" class="">
<div class=""><span style="background-color:rgba(255,255,255,0)" class="">___</span></div>
<div class=""><span style="background-color:rgba(255,255,255,0)" class="">Stefano Baroni, Trieste -- <a href="http://stefano.baroni.me/" target="_blank" class="">http://stefano.baroni.me</a></span></div>
</div>
<div dir="ltr" class=""><br class="">
<blockquote type="cite" class="">On 13 Nov 2021, at 20:29, Salman Zarrini <<a href="mailto:salman.zarrini@gmail.com" target="_blank" class="">salman.zarrini@gmail.com</a>> wrote:<br class="">
<br class="">
</blockquote>
</div>
<blockquote type="cite" class="">
<div dir="ltr" class="">
<div dir="ltr" class="">
<div dir="ltr" class="">
<div class=""><br class="">
</div>
<div class="">Dear Giovanni,</div>
<div class=""><br class="">
</div>
<div class="">Thanks for your response.</div>
<div class=""><br class="">
</div>
<div class="">Then, considering the density of states in an electronic system and what Quantum-Espresso calculates as the density of states, should we expect to have a volume-independent quantity? if I understood you correctly!</div>
<div class=""><br class="">
</div>
<div class="">
<div dir="ltr" class="">
<div dir="ltr" class="">
<div class="">
<div dir="ltr" class="">Regards,<br class="">
Salman<br class="">
</div>
</div>
</div>
</div>
</div>
<br class="">
</div>
</div>
<br class="">
<div class="gmail_quote">
<div dir="ltr" class="gmail_attr">On Sat, Nov 13, 2021 at 1:30 PM Giovanni Cantele <<a href="mailto:giovanni.cantele@spin.cnr.it" target="_blank" class="">giovanni.cantele@spin.cnr.it</a>> wrote:<br class="">
</div>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
Dear Salman,<br class="">
<br class="">
Actually, the two definitions are not mutually exclusive. The first you speak about, is the density of states per unit volume and, as you correctly mention, has units Energy^-1 Volume^-1. However, the definition of density of states a system of electrons and
has units Energy^-1:<br class="">
<br class="">
DOS(E) = sum_i Dirac_delta(E-E_i)<br class="">
<br class="">
Integral( dE DOS(E) ) = number of electrons<br class="">
<br class="">
What Quantum-Espresso calculates, is the density of states of the electron system in the unit cell of a given Bravais lattice (due to periodicity, you refer to the primitive cell). If you plot it as is, you should give it units eV^-1. However, you could need
the density of states per unit volume. In that case, you can easily obtain the unit cell volume of your system, divide the computed density of states by it, and then the resulting density-of-states-per-unit-volume has units eV^-1 au^-3 (if you express the
volume in au^3).<br class="">
<br class="">
In this case, if you integrate over the energy, you obtain number of electrons per unit volume, that is, electron density.<br class="">
<br class="">
Giovanni<br class="">
<br class="">
> On 13 Nov 2021, at 19:14, Salman Zarrini <<a href="mailto:salman.zarrini@gmail.com" target="_blank" class="">salman.zarrini@gmail.com</a>> wrote:<br class="">
> <br class="">
> Dearl all,<br class="">
> <br class="">
> As the density of states's definition implies, the electronic density of states has a unit of "Number of electronic states per Energy per Volume" or simply Volume^-1 Energy^-1. However, the "Volume^-1" is apparently missing in the unit of density of states
in literatures as well as here in manual/tutorials of Quantum-Espresso. So that the Energy^-1 is used as the unit for total density of states, atomic site projected density of states and orbital projected density of states.
<br class="">
> <br class="">
> I guess it is just a misunderstanding from my side, so, I would be thankful if one could elaborate further on that.<br class="">
> <br class="">
> Regards,<br class="">
> Salman<br class="">
> _______________________________________________<br class="">
> Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu/" rel="noreferrer" target="_blank" class="">www.max-centre.eu</a>)<br class="">
> users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank" class="">
users@lists.quantum-espresso.org</a><br class="">
> <a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank" class="">
https://lists.quantum-espresso.org/mailman/listinfo/users</a><br class="">
<br class="">
_______________________________________________<br class="">
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu/" rel="noreferrer" target="_blank" class="">www.max-centre.eu</a>)<br class="">
users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank" class="">
users@lists.quantum-espresso.org</a><br class="">
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank" class="">https://lists.quantum-espresso.org/mailman/listinfo/users</a><br class="">
</blockquote>
</div>
<span class="">_______________________________________________</span><br class="">
<span class="">Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu/" target="_blank" class="">www.max-centre.eu</a>)</span><br class="">
<span class="">users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank" class="">users@lists.quantum-espresso.org</a></span><br class="">
<span class=""><a href="https://lists.quantum-espresso.org/mailman/listinfo/users" target="_blank" class="">https://lists.quantum-espresso.org/mailman/listinfo/users</a></span></div>
</blockquote>
</div>
_______________________________________________<br class="">
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu/" rel="noreferrer" target="_blank" class="">www.max-centre.eu</a>)<br class="">
users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank" class="">users@lists.quantum-espresso.org</a><br class="">
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank" class="">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote></div>
_______________________________________________<br class="">Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu/" target="_blank" class="">www.max-centre.eu</a>)<br class="">users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank" class="">users@lists.quantum-espresso.org</a><br class=""><a href="https://lists.quantum-espresso.org/mailman/listinfo/users" target="_blank" class="">https://lists.quantum-espresso.org/mailman/listinfo/users</a></div></blockquote></div><br class=""></div></div>_______________________________________________<br class="">
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu/" rel="noreferrer" target="_blank" class="">www.max-centre.eu</a>)<br class="">
users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank" class="">users@lists.quantum-espresso.org</a><br class="">
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank" class="">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote></div>
_______________________________________________<br class="">Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" class="">www.max-centre.eu</a>)<br class="">users mailing list <a href="mailto:users@lists.quantum-espresso.org" class="">users@lists.quantum-espresso.org</a><br class=""><a href="https://lists.quantum-espresso.org/mailman/listinfo/users" class="">https://lists.quantum-espresso.org/mailman/listinfo/users</a></div></blockquote></div><br class=""></body></html>