<div dir="ltr"><div>The nonlocal contribution involves derivatives of j_l(kr). I vaguely remember that the limit k=0 is tricky, but I was convinced that it had been properly implemented many years ago</div><div><br></div><div>Paolo<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, May 16, 2020 at 8:12 AM L Y <<a href="mailto:xifntx@gmail.com">xifntx@gmail.com</a>> wrote:<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 dir="ltr">Dear Quantum ESPRESSO developers,<div><br></div><div>My question was posted initially to the user mail list, and I already obtained help from Dr. Iurii. It is abit technical, so I formulate it here again to see if anyone has ideas or suggestions. I think it could be a bug.</div><div><br></div><div>I calculated the the momentum matrix element between the valence and conduction bands in MgO, along a k-path going through the Gamma point. At the Gamma point, there is a clear discontinuity, see attachment "p_with_nonlocal.png" (Gamma point at ik=31), while the band energies are all continuous. I tried many different types of pseudopotentials (norm conserving, ultrasoft, PAW) and functionals (LDA, PBE), all showing this discontinuity, (although the SG15 ONCV pseudopotentials seem to be better). </div><div><br></div><div>When I commented out the part in "compute_ppsi.f90" involving the non-local pseudopotential, the discontinuity is gone, see attached "p_without_nonlocal.png". This leads me to believe that there is potentially a bug in the calculation of the commutator i [r, V_nl] for the Gamma point. I was wondering whether anyone has ideas or suggestions on why this occurs and where in the code I should look at to try to find the potential "bug".</div><div><br></div><div>Sorry for the long mail.</div><div><div>Best regards,</div><div>Lun Yue</div><div>Louisiana State University</div></div></div>
_______________________________________________<br>
developers mailing list<br>
<a href="mailto:developers@lists.quantum-espresso.org" target="_blank">developers@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/developers" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/developers</a><br>
</blockquote></div><br clear="all"><br>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>Phone +39-0432-558216, fax +39-0432-558222<br><br></div></div></div></div></div>