<div dir="ltr"><div>Hi,</div><div><br></div><div>generally speaking, it is expected that, as one moves to higher energies (in the unoccupied manifold), the sum of the projected DOS over all atoms/species becomes smaller than the total DOS.<br>This is because projwfc.x projects Kohn–Sham states onto atomic wavefunctions defined in the pseudopotential files. These atomic functions are constructed to describe well the valence states and, </div><div>to some extent, low-lying unoccupied states. However, higher-energy states may contain components that are not well represented within this atomic-like basis, and therefore are only partially captured</div><div> by the projections. You can directly observe this behavior in the projwfc.x output: the quantity reported after each eigenvalue (the sum of projections, i.e. |psi|^2 onto the atomic basis) is typically very close to </div><div>1 for occupied states, while it starts to decrease for unoccupied states and continues to do so with increasing energy.</div><div><br>That said, it is important to understand whether your specific observation has a physical origin or is instead a numerical artifact. A few points you may want to check:<br>- while a 36×36×1 grid for NSCF should already give a reasonable DOS, I would still test denser meshes. The presence of many peaks in your DOS </div><div>could indicate insufficient Brillouin zone sampling rather than a physical effect;</div><div><br>- identify the k-points and energies corresponding to the states at the conduction band edge in the strained system, and visualize the corresponding |psi|^2 (e.g. with pp.x). </div><div>If these states are more delocalized in the interstitial or vacuum region, rather than localized on B/N atoms or bonds, then it is expected that atomic projections do not fully capture them;</div><div><br>- is this bulk or monolayer hBN? In the monolayer case, vacuum regions can play an important role, and low-lying conduction states may have significant weight outside the atomic regions.</div><div>Did you check, if it applies, convergence with respect to the vacuum?</div><div><br>- make sure the strain has been applied consistently and that atomic positions have been properly relaxed. Geometrical inconsistencies can lead to spurious states near the band edges;</div><div><br>- sharing input files and band structure plots would help to better assess whether the observed features are physical.</div><div><br></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: <a href="mailto:giovanni.cantele@unina.it" target="_blank">giovanni.cantele@unina.it</a><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 gio 26 mar 2026 alle ore 09:35 Ms. Chandrika K. via users <<a href="mailto:users@lists.quantum-espresso.org">users@lists.quantum-espresso.org</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 class="msg4116693296035537118">




<div dir="ltr">
<div style="margin-top:1em;margin-bottom:1em;font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Hello everyone,</div>
<div style="margin-top:1em;margin-bottom:1em;font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
I am performing DFT calculations on a pristine and uniaxially strained hexagonal boron nitride (hBN) unit cell using Quantum ESPRESSO with the following computational parameters:</div>
<ul>
<li style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<b>Pseudopotentials:</b> B.pbe-n-rrkjus_psl.1.0.0.UPF and N.pbe-n-rrkjus_psl.1.0.0.UPF (USPP, PBE, scalar relativistic)</li><li style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<b>ecutwfc:</b> 60 Ry | <b>ecutrho:</b> 480 Ry</li><li style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<div role="presentation"><b>K-mesh:</b> 12×12×1 (SCF), 36×36×1 (NSCF)</div>
</li><li style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<b>Broadening:</b> ngauss = 0, degauss = 0.004 Ry, DeltaE = 0.004 eV</li><li style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<div role="presentation"><b>DOS tool:</b> projwfc.x only (both total DOS and PDOS from same run)</div>
</li><li style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<div role="presentation"><b>Basis:</b> Only s and p orbitals present for both B and N</div>
</li></ul>
<ul style="margin-top:0px;margin-bottom:0px;list-style-type:disc">
<li style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<div role="presentation" style="margin-top:1em;margin-bottom:1em">
In the pristine unit cell DOS (Figure 1), the total DOS and the sum of B and N projected DOS match upto certain point in the conduction band.</div>
</li></ul>
<ul>
<li style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<div role="presentation" style="margin-top:1em;margin-bottom:1em">
However, in the strained unit cell (Figure 2), additional states at the edge of conduction band are visible  and also the total DOS is
<b>not accounted </b> by either the B or N projected DOS individually or collectively.</div>
<div role="presentation">Any insights on the origin of the extra conduction band states and the enhanced PDOS mismatch in the strained system would be greatly appreciated.</div>
<p role="presentation" style="text-align:left;text-indent:0px;background-color:rgb(255,255,255);margin:0px">
Chandrika Yadav K</p>
<p role="presentation" style="text-align:left;text-indent:0px;background-color:rgb(255,255,255);margin:0px">
<br>
</p>
<div role="presentation" style="text-align:left;text-indent:0px;background-color:rgb(255,255,255);margin:0px">
Research Scholar</div>
<p role="presentation" style="text-align:left;text-indent:0px;background-color:rgb(255,255,255);margin:0px">
NanoElectronics Laboratory </p>
<p role="presentation" style="text-align:left;text-indent:0px;background-color:rgb(255,255,255);margin:0px">
Amrita Viswa Vidyapeetham </p>
<p role="presentation" style="text-align:left;text-indent:0px;background-color:rgb(255,255,255);margin:0px">
Bangalore</p>
<div role="presentation"><br>
</div>
</li></ul>
<div style="font-family:Aptos,Aptos_EmbeddedFont,Aptos_MSFontService,Calibri,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
</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></div></blockquote></div>