<div dir="ltr"><div>It seems gmail did not send my last email, so I am trying again:</div><div><br></div>Dear All,<div><br></div><div>A bit more information on my lapw7 issue.</div><div><br></div><div>First of all, my degenerate states are not orthogonal. Perhaps their linear combinations are.</div><div><br></div><div>I calculate charge density for a particular k and a particular state number. I work with a symmetric slab of the 1T TMDC lattice. For a surface state, this means that what is spin-up on side A of the slab should be spin-dn on side B of the slab, |Psi_k_A(r)| = |Psi_k_B(-r)| (sides A and B are symmetric, r=0 is e.g. in the middle of the slab).</div><div><br></div><div>Now, I have 2 degenerate eigenvalues at some energy, in my case e.g. eigenvalue #1461 and #1462. In principle one of them should host spin-up state on side A of the slab, and the other should host spin-dn state at side B of the slab. However, this is not the case, and the charge of the surface states is shared by the eigenvalues (this also makes the states non-orthogonal).</div><div><br></div><div>lapw7 can calculate the full wave function with the phase. Now the following DOES NOT hold:</div><div><br></div><div>NOT TRUE: |Psi_k_#1461_A_up(r) + Psi_k_#1462_A_up(r)| = |Psi_k_#1461_B_up(-r) <a class="gmail_plusreply" id="m_631210883791230767plusReplyChip-1">+</a> Psi_k_#1462_B_up(-r)|</div><div><br></div><div>But the following holds:</div><div><br></div><div><div>TRUE: |Psi_k_#1461_A_up(r) + Psi_k_#1462_A_up(r)| = |Psi_k_#1461_B_up(-r) - Psi_k_#1462_B_up(-r)|</div><div><br></div></div><div>This means that the spin-dn wave functions need to be added out-of-phase to obtain what is expected. Actually, from the 2D plot of the phase of Psi one can see that, for the two degenerate eigenvalues, spin-up wave functions are in-phase, and spin-dn wave functions are out-of-phase.</div><div><br></div><div>Can you comment and advise how to deal with that?</div><div><br></div><div>Best,</div><div>Lukasz</div></div>