<div dir="ltr"><div>Dear Paolo,</div><div><br></div><div>Thank you for the quick response and pointing me in a direction! </div><div>I see that in "compute_ppsi.f90", the subroutines "gen_us_dj" and "gen_us_dy" are called, inside of which the derivatives of the spherical Bessel and spherical harmonics are calculated. I think the Gamma point is special in the sense that it is the only point that can satisfy k+G=0, so I thought perhaps the problem could be due to "dvkb" and "dvkb1" being treated inaccurately at k+G=0. I checked that in "compute_ppsi.f90", for k+G=0, the array "work(ig,:)=dvkb1(ig,:)+dvkb(ig,:)*gk(ipol,ig)" is zero. The second term is obviously zero due to |k+G|=0, but I cannot see why "dvkb1(ig,:)" is zero. </div><div>I think I am abit too inexperienced with the theory and code to find the "bug", I only pieced together part of the puzzle by looking at <span style="white-space:pre-wrap">Iurii's thesis Sec 6.13. </span>I dunno if anyone has time to look at this problem?<span style="white-space:pre-wrap"> I think in my emerging physical </span>discipline<span style="white-space:pre-wrap"> (attosecond strong-field physics in solids), the correct calculation of the momentum matrix elements will be needed alot in the future (to be used in the TDSE or time-dependent density-matrix equations). But maybe this is not so urgent at the moment, since only the Gamma point is an issue, so I can just circumvent this by not choosing that point. </span></div><div><br></div><div>Best regards,</div><div>Lun Yue</div><div>Louisiana State University <br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Den lør. 16. maj 2020 kl. 01.44 skrev Paolo Giannozzi <<a href="mailto:p.giannozzi@gmail.com">p.giannozzi@gmail.com</a>>:<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"><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" target="_blank">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"><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>
</blockquote></div>