Dear friends from CAS and Italy,<br><br>Thank you for your advice, I will try to adjust the parameters that you just mentioned above.<br><br>At last, I found QE has been upgraded to 4.1.1, but the QE I used is 4.0.3. I wonder whether the new vision of QE could do better than the older one.<br>
<br>Thank you again for your time.<br><br>best,<br><br>vega <br><br><div class="gmail_quote">On Fri, Nov 6, 2009 at 9:00 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 06 novembre 2009 alle ore 04:35:39, vega lew<br>
<<a href="mailto:quantumdft@gmail.com">quantumdft@gmail.com</a>> ha scritto:<br>
<div class="im">> During my NEB calculations, I found there are two possible situation.<br>
> First, the calculation could be converged after about 100 iterations.<br>
> Second, the calculation never converged in more than 300 iterations.<br>
<br>
</div>Dear Vega,<br>
NEB calculation are a bit tricky in general and require extreme care to be<br>
setup correctly. NEB also takes easily hunders of iteration to converge,<br>
of course depending on the number of atoms and of images. Here is some<br>
free advice I can give, although only practice make prefection.<br>
<br>
1. don't use climbing image at the beginning, CI makes convergence slower;<br>
especially if the special image changes during the convergence process (it<br>
may happen if CI_scheme='auto') this does ususally mess up everything.<br>
Converge your calculation than restart from the last configuration with<br>
the climbing image (note that this will *increase* the barrier!)<br>
<br>
2. carefully choose the initial path, remember that QE assumens continuity<br>
between the first and last image at the initial condition. In other words,<br>
periodic images are NOT used, you may have to manually translate an atom<br>
by one or more unit cell base vectors in order to have a meaningful<br>
initial path.<br>
<br>
3. you can try to start the NEB process with most of the atoms fixed in<br>
position, in order to have the more "problematic" one to converge, before<br>
freeing all of them.<br>
<br>
4. especially for larger systems, you can start the NEB with lower<br>
accuracy (less k-points, lower cutoff) and then increase it when it has<br>
converged to refine your calculation.<br>
<br>
<br>
5. use the Broyden algorithm instead of the default one, it is a bit more<br>
fragile, but it removes the oscillation problem. If this happens, and you<br>
cannot afford more images, focus to a smaller problem, decompose it to<br>
pieces.<br>
<br>
6. you can grossly estimate the number of iterations required as<br>
number_of_images * number_of_atoms * 3 (the atoms that don't move should<br>
not be counted). It may take half that many iterations, or twice as many,<br>
but more or less that's the order of magnitude, unless you start from a<br>
very accurate or very bad initial guess.<br>
<br>
good luck<br>
<br>
--<br>
Lorenzo Paulatto<br>
SISSA  &  DEMOCRITOS (Trieste)<br>
phone: +39 040 3787 511<br>
skype: paulatz<br>
www:   <a href="http://people.sissa.it/%7Epaulatto/" target="_blank">http://people.sissa.it/~paulatto/</a><br>
<br>
     *** save italian brains ***<br>
  <a href="http://saveitalianbrains.wordpress.com/" target="_blank">http://saveitalianbrains.wordpress.com/</a><br>
_______________________________________________<br>
<div><div></div><div class="h5">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>
</div></div></blockquote></div><br><br clear="all"><br>-- <br>==================================================================================<br>Vega Lew ( weijia liu)<br>Graduate student<br>State Key Laboratory of Materials-oriented Chemical Engineering<br>
College of Chemistry and Chemical Engineering<br>Nanjing University of Technology, 210009, Nanjing, Jiangsu, China<br>******************************************************************************************************************<br>
Email: <a href="mailto:vegalew@gmail.com">vegalew@gmail.com</a><br>Office: Room A705, Technical Innovation Building, Xinmofan Road 5#, Nanjing, Jiangsu, China  <br>****************************************************************************************************************** <br>