<div dir="ltr"><div><br></div><div>I'm pretty sure forces and stresses were tested in all +U flavors, at least for the version of the code where they were originally implemented.</div><div>It is true that +U typically leads to the creation of several local minima and maybe this makes the convergence harder. However, I have never seen </div><div>a total force staying so constant (and big) as you show. Are the ions moving at all? How big is your trust radius? </div><div>Does the relaxation converge if you set Hubbard_J0 = 0? </div><div>Also, is "exxdiv_treatment = 'gygi-baldereschi' " needed?</div><div><br></div>Matteo<div><br><div class="gmail_extra"><br><div class="gmail_quote">On Fri, May 5, 2017 at 10:37 AM, Lorenzo Paulatto <span dir="ltr"><<a href="mailto:lorenzo.paulatto@impmc.upmc.fr" target="_blank">lorenzo.paulatto@impmc.upmc.fr</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><span class="gmail-">On Friday, 5 May 2017 09:26:53 CEST Q.J.Wang wrote:<br>
>              lda_plus_u_kind = 0 ,<br>
>                 Hubbard_U(3) = 4.5,<br>
>                Hubbard_J0(3) = 0.5,<br>
</span>>...<br>
<span class="gmail-">> The total force seems nearly invariable. I also use QE5.4.0 to run the input<br>
> file, the total force still can't  converge after 300 scf cycles and  2 bfgs<br>
> steps.<br>
<br>
</span>My guess is that there is an inconsistency between LDA+U energy and forces,<br>
which prevents the BFGS calculation from converging.<br>
<br>
kind regards<br>
<span class="gmail-HOEnZb"><font color="#888888"><br>
--<br>
Dr. Lorenzo Paulatto<br>
IdR @ IMPMC -- CNRS & Université Paris 6<br>
phone: <a href="tel:%2B33%20%280%291%20442%2079822" value="+33144279822">+33 (0)1 442 79822</a> / skype: paulatz<br>
www:   <a href="http://www-int.impmc.upmc.fr/~paulatto/" rel="noreferrer" target="_blank">http://www-int.impmc.upmc.fr/~<wbr>paulatto/</a><br>
mail:  23-24/423 Boîte courrier 115, 4 place Jussieu 75252 Paris Cédex 05<br>
</font></span><div class="gmail-HOEnZb"><div class="gmail-h5"><br>
______________________________<wbr>_________________<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" rel="noreferrer" target="_blank">http://pwscf.org/mailman/<wbr>listinfo/pw_forum</a><br>
</div></div></blockquote></div><br></div></div></div>