<div dir="ltr"><div>Hi,</div><div><br></div><div>from the header of PP/src/average.f90</div><div>  !      Note that if Q is a charge density whose integral is Z_v:<br>  !         Z_v = \int p(z) dV = \sum_k p(k) \Omega/N_3<br>  !      where \Omega is the size of the unit cell (or supercell)</div><div><br></div><div>So, if I remember, the units should be au^-3, while the horizontal axis is au.</div><div><br></div><div>Just to cross check, take the planar average of a given system (e.g. your eterostructure) and integrate the 1D function as is.</div><div>Supposing that you chose idir = 3 and that a3 is orthogonal to a1,a2, you should get the number of electrons divider by the surface of</div><div>the unit cell defined by a1 and a2.</div><div><br></div><div>Let me know.</div><div>Giovanni</div><div><br></div><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><span style="color:rgb(0,0,0)">-- <br><br>Prof. Giovanni Cantele<br>Dipartimento di Fisica "Ettore Pancini"<br>Universita' degli Studi di Napoli "Federico II"<br>Complesso Universitario M. S. Angelo - Ed. 6<br>Via Cintia, I-80126, Napoli, Italy<br>e-mail: giovanni.cantele@unina<br>Phone: +39 081 676910<br><br>Web page: </span><a href="https://sites.google.com/view/giovanni-cantele/home" style="color:rgb(17,85,204)" target="_blank">https://sites.google.com/view/giovanni-cantele/home</a><br></div></div></div><br></div><br><div class="gmail_quote gmail_quote_container"><div dir="ltr" class="gmail_attr">Il giorno mer 8 ott 2025 alle ore 11:26 jkmodi22 <<a href="mailto:jkmodi22@iitk.ac.in">jkmodi22@iitk.ac.in</a>> ha scritto:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div style="font-size:10pt;font-family:Verdana,Geneva,sans-serif">
<p>Dear QE users,<br><br>I was trying to calculate the planar-averaged charge density difference for heterostructure CdCl2-SiH. I was trying to reproduce Figure 5b) of the reference paper "Versatility of the SiH−CdCl2Heterostructure: Piezoelectricity, Photocatalysis, and Transistor Applications".<br><br>I calculated the planar averaged charge density for the hetero-structure, then for monolayer SiH, then for CdCl2.<br><br>Then delta_rho(z) was calculated = planar averaged charge density for the hetero-structure - planar averaged charge density for SiH - planar averaged charge density for CdCl2.<br><br>My question is, <strong>what is the unit of the third column of average.out files?</strong> <strong>is it e/Bohr?</strong> I know the first column is in Bohr. <br><br>pp.x page says this : </p>
<pre style="color:rgb(34,34,34);font-size:14px;background-color:rgb(255,255,255)">All output quantities are in ATOMIC (RYDBERG) UNITS unless
otherwise explicitly specified.
All charge densities integrate to the NUMBER of electrons
not to the total charge.
All potentials have the dimension of an energy (e*V, not V).</pre>
<p><br>Therefore, I guessed the third column to be e/Bohr by equation no. 5 of the reference paper, which calculates delta_Q(z) as an integral of delta_rho(z)</p>
<p><br><br>Thanking you.<br><br>Sincerely,<br><br>Jainandan Modi<br>Ph.D. scholar<br>EE Dept.<br>IIT Kanpur<br>India<br><br><br><br><br><br></p>

</div>
_______________________________________________________________________________<br>
The Quantum ESPRESSO Foundation stands in solidarity with all civilians worldwide who are victims of terrorism, military aggression, and indiscriminate warfare.<br>
--------------------------------------------------------------------------------<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>