<div dir="ltr"><div dir="ltr"><div>Dear Giovanni and Paolo,</div><div><br></div><div>Thanks very much for your suggestions.</div><div><br></div><div>I reconsidered the definition of the local density of states (say, LD(x,y,z,E)). Since it is "local", the unit of it should be states/eV/bohr^3 or electrons/eV/bohr^3. Therefore, the integration of it within an energy window should lead to the charge density in this energy window: electrons/bohr^3 . Therefore, if we choose the energy window from the lowest energy to the Fermi energy, we should get exactly the total charge density.</div><div><br></div><div>So I did some tests following Giovanni's suggestion, using a simple case of momolayer MoS2, which has "number of electrons  =   26.00". The area of its xy plane is S. For those tests, the previous scf and nscf calculations are the same, with LDA USPP,  the occupation= 'fixed', and k-points of 40*40*1 for scf and 100*100*1 for nscf.</div><div><br></div><div>1. I  calculated the total charge density (rho(x,y,z)) using pp.x with plot_num = 0, and then calculate the planar average of it (rho_avg(z)). Then I integrated it by \int (S* rho_avg(z)) dz and I got 25.89. It is close to 26.00, maybe a more accurate value can be obtained by a calculation with denser k points.</div><div><br></div><div>2. I calculated the integrated local density of states  ( ILD(x,y,z)) from -65 eV (this energy is below the lower band) to Fermi energy using pp.x with plot_num =10, and and then calculate the planar average of it ( ILD_avg(z) ). When I integrated it by \int (S*

ILD_avg(z)) dz, I got 29.25, which is much larger than the total number of electrons of 26.</div><div><br></div><div> So it seems the the test 2 doesn't correctly reflect the  reality. I am not sure it is due to something happened with plot_num = 10 in pp.x, or just I understand this incorrectly. </div><div><br></div><div>Any suggestions? Thank you very much!</div><div><br></div><div>Best,</div><div><br></div><div>Ding-Fu</div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><blockquote style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex" class="gmail_quote">From: Giovanni Cantele <<a href="mailto:giovanni.cantele@spin.cnr.it" target="_blank">giovanni.cantele@spin.cnr.it</a>><br>To: Quantum Espresso users Forum <<a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a>><br>Cc: <br>Bcc: <br>Date: Fri, 26 Oct 2018 09:51:40 +0200<br>Subject: Re: [QE-users] Unit for the output of average.x<br>Dear Ding-Fu,<br>as far as I remember there is a surface factor that you need to adjust units.<br>For sure on the ascissa axis the coordinate is in bohr.<br>The planar average give you back a quantity with the same units as the averaged quantity<br>(e.g. if you star from charge density in electrons/bohr^3 you get an averaged electron density in electrons/bohr^3),<br>being defined as (let us suppose that you average in the plane defined by a1 and a2 vectors):<br>rho_avg(z) = ( 1 / S ) * integral( dx dy rho(x,y,z) )<br>That means that if you perform <br> integral( dz rho_avg(z) )<br>you get <br>number of electrons / S<br>If you need number of electrons than just multiply by S with<br>S = cross_product( a1, a2 )<br>(in bohr^2)<br>Just try, better if you do it with the total charge density, to check if the integral returns you<br>the number of electrons.<br>I’m sorry but I cannot check directly if I remember correctly at the moment, but<br>this should work.<br>Giovanni<br>-- <br>Giovanni Cantele, PhD<br>CNR-SPIN<br>c/o Dipartimento di Fisica<br>Universita' di Napoli "Federico II"<br>Complesso Universitario M. S. Angelo - Ed. 6<br>Via Cintia, I-80126, Napoli, Italy<br><a href="mailto:giovanni.cantele@spin.cnr.it" target="_blank">e-mail: giovanni.cantele@spin.cnr.it<br></a>            <a href="mailto:gcantele@gmail.com" target="_blank">gcantele@gmail.com<br></a>Phone: +39 081 676910<br>Skype contact: giocan74<br>Web page:<span class="gmail-m_4347110038969224771Apple-converted-space"> </span><a href="https://sites.google.com/view/giovanni-cantele" target="_blank">https://sites.google.com/view/giovanni-cantele</a><br>On 26 Oct 2018, at 03:31, Dingfu Shao <<a href="mailto:dingfu.shao@gmail.com" target="_blank">dingfu.shao@gmail.com</a>> wrote:<br>Dear QE developers and users:<br>I am wondering what should be the unit of the planar average data got from the average.x<br> I am calculating the planar average of charge density within a energy window. What I did is firstly using pp to get the integrated local density of states (ILDOS) of that energy window with plot_num=10, then using average.x to get the planar average. <br>In this case, what is the unit of the second column (say, rho(z)) of  the output file? I thought since the DOS has a unit of states/eV, the integration of DOS within a energy window should get some states or electrons. Then the unit of<span style="font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13.33px;line-height:19.99px"> rho(z)</span> should be electron/bohr. But seems it is not. In my case the energy window I concerned contains one electron, However, if I directly integrate <span style="font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13.33px;line-height:19.99px"> rho(z)</span>, I can only get a very small value. If I assume the unit is <span style="font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;font-size:13.33px;line-height:19.99px">electron/(bohr^3), the integretion of <span style="font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;line-height:19.99px"> rho(z)*A is also smaller than one (here A is the area of xy plane). </span><span style="font-variant-numeric:normal;font-variant-east-asian:normal;font-stretch:normal;line-height:19.99px"> </span></span><br>Can you help me about it? Thank you very much!<br>Best,<br>Ding-Fu<br><br><br clear="all"><span style="color:blue;font-size:12pt"><strong>Ding-Fu Shao, Ph. D.</strong></span><br><strong><font color="#0000ff" face="Verdana"><span style="font-size:16px">Department of Physics and Astronomy, University of Nebraska-Lincoln</span></font></strong><br><font color="#0000ff" style="font-family:Verdana"><span style="font-size:16px"><b>Lincoln, NE </b></span></font><font color="#0000ff" face="Verdana"><span style="font-size:16px"><b>68588-0299</b></span></font><br><b><span style="color:blue;font-size:12pt"><span style="color:rgb(192,192,192);font-size:small"><b><span style="color:blue;font-size:12pt">Email: </span></b><span style="color:blue;font-size:12pt"><font color="#0000ff"><a href="mailto:dingfu.shao@gmail.com" target="_blank"><b>dingfu.shao@gmail.com</b></a></font></span></span></span></b><br>_______________________________________________<br>users mailing list<br><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" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote><div style="overflow-wrap: break-word;"><blockquote type="cite"><br class="gmail-m_4347110038969224771Apple-interchange-newline"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div><b></b><i></i><u></u><strike></strike><b></b><i></i><u></u><strike></strike><b></b><i></i><u></u><strike></strike></div></div></div></div></div></div></div></div></div></div></div></blockquote></div></div></div>