<div dir="ltr"><div>Thank you for reporting this. I heard about this problem and discussed it with Stefano dG who knows well that part of code. He thinks that the current implementation is correct and that the discrepancy with respect to analytical formulas is only apparent: in the limit |grad rho| =>0 it is convenient to have a single diverging

quantity 

(dq0_dgradrho)

 instead of two (dq0_dgradrho and grad rho/mod_grad), or something like that</div><div><br></div><div>Paolo<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Apr 24, 2019 at 8:55 PM <<a href="mailto:tran@theochem.tuwien.ac.at">tran@theochem.tuwien.ac.at</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">Hello,<br>
<br>
In xc_rVV10.f90, the formula for dq0_dgradrho (line 621) contains an error.<br>
Instead of<br>
dq0_dgradrho(i_grid) = dq0_dq  / ( 2.0_dp*k*w0 ) * 4.0_dp*wg2 / (mod_grad**2)<br>
this should be<br>
dq0_dgradrho(i_grid) = dq0_dq  / ( 2.0_dp*k*w0 ) * 4.0_dp*wg2 / mod_grad<br>
<br>
The authors of PBEsol+rVV10s<br>
<a href="https://journals.aps.org/prb/abstract/10.1103/PhysRevB.98.214108" rel="noreferrer" target="_blank">https://journals.aps.org/prb/abstract/10.1103/PhysRevB.98.214108</a><br>
told me that they arrived at the same conclusion.<br>
<br>
Best regards,<br>
<br>
F. Tran (TU Wien)<br>
_______________________________________________<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>