<div dir="ltr">dV/dR is by definition the limit of Delta V/Delta R for Delta R=>0. If this does not converge to a well defined function for dR->0, there is something wrong in your calculation <br></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Jun 15, 2015 at 4:56 PM, Jiqiang Li <span dir="ltr"><<a href="mailto:jqli14@fudan.edu.cn" target="_blank">jqli14@fudan.edu.cn</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div>
<div><span></span>Dear developers,</div><div>        I want to ask a trivial question about total local potential Vtot (V_Hartree + Vxc + V_ion), i.e., the variable <b>vrs </b>in electron.f90 (correct me if I am wrong). When we move a specific atom along a specific direction, does the derivative quantity dVtot/dR remain unchanged ?</div><div>       I try 0.01~0.1 Angstrom as the amplitude of the displacements, and find that dVtot/dR is not independent of dR. does the differece arise from the numerical error or theoretically <span style="font-size:10.5pt;line-height:1.5;background-color:window">should</span><span style="font-size:10.5pt;line-height:1.5;background-color:window"> </span><span style="font-size:10.5pt;line-height:1.5;background-color:window">dVtot/dR depend on dR? If the former is right, then how long is the optimal value of dR for this calculation generally?</span></div><div><br></div><div>Best regards!</div><span class="HOEnZb"><font color="#888888"><div><br></div><div>Jiqiang Li</div><blockquote style="margin-top:0px;margin-bottom:0px;margin-left:0.5em"><div><div><div class="gmail_extra"><div></div>
</div>
</div></div></blockquote>
</font></span></div></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature"><div dir="ltr"><span><span><font color="#888888">Paolo Giannozzi, Dept. Chemistry&Physics&Environment,<br>
Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>
hone <a href="tel:%2B39-0432-558216" value="+390432558216" target="_blank">+39-0432-558216</a>, fax <a href="tel:%2B39-0432-558222" value="+390432558222" target="_blank">+39-0432-558222</a></font></span></span></div></div>
</div>