<font face="arial,helvetica,sans-serif">Dear Guido,<br>       Thank you very much for your advice.<br>        Yes, it is a slab. The problem is it is difficult to converge if the dipole correction is added because in our system the charge density is changing too abruptly. Therefore, I was wondering whether adding a dipole correction which is keep constant during the scf convergence can improve the convergence problem. Of course, this will depends very much on the initial value of the dipole field.<br>
        Fenggong <br></font><br><div class="gmail_quote">On Tue, Oct 9, 2012 at 3:41 AM, Guido Fratesi <span dir="ltr"><<a href="mailto:fratesi@mater.unimib.it" target="_blank">fratesi@mater.unimib.it</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0pt 0pt 0pt 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Dear Fenggong,<br>
<br>
I assume you are computing properties for a 2D system such as a slab, or<br>
0D such as molecule-in-a-box.<br>
<br>
You can use variables in input (Doc/INPUT_PW):<br>
  tefield (necessary in any case below)<br>
  dipfield for a scf correction (which is determined by density)<br>
  eamp /= 0 to add a constant value added.<br>
<br>
I you want to dig into the code, PW/src/add_efield.f90 is the routine of<br>
your interest.<br>
<br>
Hope this helps,<br>
Guido<br>
<div><div class="h5"><br>
On 10/08/2012 06:34 PM, Fenggong Wang wrote:<br>
> Dear Everyone,<br>
>         I am relatively new to Quantum Espresso.<br>
>         Could we keep the dipole field at a constant value during the<br>
> scf interations or constaint it to be in a value range when adding the<br>
> dipole correction in pwscf? Which part of the code should we refer to<br>
> for this? Any guidance or comments on this?<br>
>        Thank you very much.<br>
> Fenggong<br>
><br>
> --<br>
> Dr. Fenggong Wang<br>
> Postdoctoral Researcher<br>
><br>
> Department of Chemistry<br>
> University of Pennsylvania<br>
> 231 South 34th Street, BOX 314<br>
> Philadelphia, PA 19104-6323<br>
</div></div>> e-mail: <a href="mailto:wangfenggong@gmail.com">wangfenggong@gmail.com</a> <mailto:<a href="mailto:wangfenggong@gmail.com">wangfenggong@gmail.com</a>>;<br>
> <a href="mailto:fenggong@sas.upenn.edu">fenggong@sas.upenn.edu</a> <mailto:<a href="mailto:fenggong@sas.upenn.edu">fenggong@sas.upenn.edu</a>><br>
><br>
><br>
> _______________________________________________<br>
> Pw_forum mailing list<br>
> <a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
> <a href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br>
><br>
<span class="HOEnZb"><font color="#888888"><br>
--<br>
Guido Fratesi<br>
<br>
Dipartimento di Scienza dei Materiali<br>
Universita` degli Studi di Milano-Bicocca<br>
via Cozzi 53, 20125 Milano, Italy<br>
<br>
_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br>
</font></span></blockquote></div><br><br clear="all"><br>-- <br><span>Dr. Fenggong Wang</span><br><span>Postdoctoral Researcher</span><br><br><span>Department of Chemistry</span><br><span>University of Pennsylvania</span><br>
<span>231 South 34th Street, BOX 314</span><br><span>Philadelphia, PA 19104-6323</span><br><font color="#222222" face="arial, sans-serif">e-mail: <a href="mailto:wangfenggong@gmail.com" target="_blank">wangfenggong@gmail.com</a>; <a href="mailto:fenggong@sas.upenn.edu" target="_blank">fenggong@sas.upenn.edu</a></font><br>