Thank you, Axel.<br>

<br>

My calculation is modified from example31, and the script is attached below.<br>

The field strength i set for calculation is 1.0e-3 a.u.,  the same as

all examples implemented in Q.E.   I want to do a 'relax' calculation

after a 'scf' calculation under finite-field, but  the calculation

always cannot reach a successful convergence after one or more 'bfgs'

steps.<br>

<br>

 Is there any trick to accelerate or reach  scf convergence ?  I have

tried several times with different schemes  , such as increasing bands

number ,   adding smearing , or decreasing /increasing  mixing

parameters.  Unfortranately ,  there seems unsuccessful. <br>

<br>

<br>

Regards,<br>

<br>

H.P<br><br><div><span class="gmail_quote">On 7/5/07, <b class="gmail_sendername">Axel Kohlmeyer</b> <<a href="mailto:akohlmey@cmm.chem.upenn.edu">akohlmey@cmm.chem.upenn.edu</a>> wrote:</span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">

On Thu, 5 Jul 2007, lan haiping wrote: HPL> Dear All, HPL> HPL>   I am studying electric-enthalpy method implemented in pw.x with HPL> example31.  There are some problems i cannot solve by myself.

HPL> Firstly,  when i set "ecutrho" explictly, calculations under finite-field HPL> are failed without any complain information, such as CRASH report. My server please provide an input that can reproduce this.

<br><br>HPL> is of 8GB memory, and no such crash for other scf calculations . I had tried<br>HPL> 'cg' and 'david' methods, and cannot solve this  problems.<br>HPL> Secondly,  when i set occupations to 'mv' smearing scheme with degauss=

0.03 HPL> a.u. , the scf calculation of finite-field situation HPL> is very difficult to reach convergence.  I have tried to increase or calculations under external field _are_ more difficult to converge.

<br>how does your field strength compare to the value in example31?<br><br>HPL> decrease mixing_ndim,  but failed.  Besides, i have tried to increase<br>HPL> nbnd , which seems no effect on scf convergence.  Would you please give me

<br>HPL> some suggestion on parameters setting for scf calculation under finite-field<br>HPL> ?<br><br>it would be much easier to detect problems or make suggestions,<br>if you'd post the relevant part of your input file (or another

<br>input that has the exact same problem).<br><br>cheers,<br>   axel.<br><br>HPL><br>HPL> Thanks for your time and attention . Any suggestion or comment is<br>HPL> appreciated .<br>HPL><br>HPL> Best Wishes,

<br>HPL> H.P<br>HPL><br>HPL><br>HPL><br>HPL><br><br>--<br>=======================================================================<br>Axel Kohlmeyer   <a href="mailto:akohlmey@cmm.chem.upenn.edu">akohlmey@cmm.chem.upenn.edu

</a>   <a href="http://www.cmm.upenn.edu">http://www.cmm.upenn.edu</a><br>   Center for Molecular Modeling   --   University of Pennsylvania<br>Department of Chemistry, 231 S.34th Street, Philadelphia, PA 19104-6323<br>tel: 1-215-898-1582,  fax: 1-215-573-6233,  office-tel: 1-215-898-5425

<br>=======================================================================<br>If you make something idiot-proof, the universe creates a better idiot.<br><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">http://www.democritos.it/mailman/listinfo/pw_forum</a><br></blockquote></div><br><br clear="all">

<br>-- <br>Hai-Ping Lan <br>Department of Electronics ,<br>Peking University , Bejing, 100871<br><a href="mailto:lanhaiping@gmail.com">lanhaiping@gmail.com</a>, <a href="mailto:hplan@pku.edu.cn">hplan@pku.edu.cn</a>