<div dir="ltr"><div dir="ltr"><div class="gmail_default" style="font-family:tahoma,sans-serif">Hi All,</div><div class="gmail_default" style="font-family:tahoma,sans-serif"><br></div><div class="gmail_default" style="font-family:tahoma,sans-serif">I am studying monolayer black phosphorus using version 6.0. The primitive unit cell of the crystal contains four atoms each with four valence orbitals (s, px, py, pz). So, in total there are 16 orbitals in the unit cell. <br></div><div class="gmail_default" style="font-family:tahoma,sans-serif">I run scf and nscf calculations by setting nbnd=16. Then using projwfc module I get the projection of the wavefunctions over atomic orbitals at k=0 by setting kresolved=.true.. <br></div><div class="gmail_default" style="font-family:tahoma,sans-serif">Also, I am using a norm-conserving scalar relativistic pseudopotential. <br></div><div class="gmail_default" style="font-family:tahoma,sans-serif">Here are the results for the 16 energy bands at the Gamma point (k=0). <br></div><div class="gmail_default" style="font-family:tahoma,sans-serif"><br></div><div class="gmail_default" style="font-family:tahoma,sans-serif"><span style="font-family:monospace,monospace">==== e(   1) =   -19.25235 eV ==== |psi|^2 = 0.997<br>==== e(   2) =   -16.76376 eV ==== |psi|^2 = 0.994<br>==== e(   3) =   -13.46989 eV ==== |psi|^2 = 0.996<br>==== e(   4) =    -9.99107 eV ==== |psi|^2 = 0.977<br>==== e(   5) =    -8.95201 eV ==== |psi|^2 = 0.996<br>==== e(   6) =    -8.74968 eV ==== |psi|^2 = 0.983<br>==== e(   7) =    -7.87978 eV ==== |psi|^2 = 0.987<br>==== e(   8) =    -6.19935 eV ==== |psi|^2 = 0.977<br>==== e(   9) =    -5.92184 eV ==== |psi|^2 = 0.981<br>==== e(  10) =    -4.28819 eV ==== |psi|^2 = 0.979<br>==== e(  11) =    -3.24853 eV ==== |psi|^2 = 0.945<br>==== e(  12) =    -2.54804 eV ==== |psi|^2 = 0.735<br>==== e(  13) =    -1.77330 eV ==== |psi|^2 = 0.905<br>==== e(  14) =    -1.10489 eV ==== |psi|^2 = 0.978<br>==== e(  15) =     0.18973 eV ==== |psi|^2 = 0.008<br>==== e(  16) =     0.25455 eV ==== |psi|^2 = 0.009</span></div></div><div dir="ltr"><div style="font-family:tahoma,sans-serif" class="gmail_default"></div><div style="font-family:tahoma,sans-serif" class="gmail_default">As you can see not all the bands are normalized, especially the last two ones (<span style="font-family:monospace,monospace">|psi|^2 << 1</span>). <br></div><div style="font-family:tahoma,sans-serif" class="gmail_default">I should mention that scf is converging and there are no diagonalization issues. <br></div><div style="font-family:tahoma,sans-serif" class="gmail_default"><br>My question is that why are these wavefunctions not normalized? <br></div><div style="font-family:tahoma,sans-serif" class="gmail_default">Do we expect them to be normalized at all? Why not? <br></div><div style="font-family:tahoma,sans-serif" class="gmail_default"><br></div><div style="font-family:tahoma,sans-serif" class="gmail_default">Thanks in advance,</div><div style="font-family:tahoma,sans-serif" class="gmail_default">Seyed <br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div><div><div><font size="4"><span style="font-family:garamond,serif"><i>S. M. Farzaneh</i></span></font><br></div><font size="1"><span style="font-family:georgia,serif">Graduate Student<br></span></font></div><div><font size="1"><span style="font-family:georgia,serif">ECE Department<br></span></font></div><div><font size="1" face="georgia, serif">New York University</font></div><font size="1"><span style="font-family:georgia,serif"><br></span></font></div><div><font size="1"></font></div><span style="font-family:georgia,serif"></span><br></div></div></div></div></div></div></div></div></div>