Dear Paulatto,<br><br>I compared the source codes under PW directory between 4.0.4 and 4.0.5. I only found a minor relevant change in compute_dip.f90. I don't think this change will affect the force.<br><br>Best regards<br>
<br>zah<br><br><div class="gmail_quote">On Tue, Apr 21, 2009 at 10:21 PM, Lorenzo Paulatto <span dir="ltr"><<a href="mailto:paulatto@sissa.it">paulatto@sissa.it</a>></span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
In data 21 aprile 2009 alle ore 15:55:29, Aihua Zhang <<a href="mailto:zah7903@gmail.com">zah7903@gmail.com</a>><br>
ha scritto:<br>
<div class="im">> I had expected the force exerted on nucleus by the field  would<br>
> cancel the force on ion due to change in electron charge density, so the<br>
> total force should vanish since the system is neutral. I also varied the<br>
> cutoff energy, but failed to see a convergence behavior as shown below.<br>
> Could anyone give me a clue? Thank you!<br>
<br>
</div>dear Aihua Zhang,<br>
that part of the code has been "fixed" more than once recently, hence the<br>
answer depends on the exact version of QE you are using.<br>
<br>
regards<br>
<font color="#888888"><br>
<br>
--<br>
Creato con il rivoluzionario client e-mail di Opera:<br>
<a href="http://www.opera.com/mail/" target="_blank">http://www.opera.com/mail/</a><br>
_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
</font></blockquote></div><br><br clear="all"><br>-- <br>------------------------------------------------------------------<br>Research Fellow,<br>Dept. of Physics, National Univ. of Singapore<br>10 Kent Ridge Crescent, Singapore, 119260<br>
Tel: +65 6516 7844<br>------------------------------------------------------------------<br><br>