<div dir="ltr"><div><div><div>Dear Pietro.<br></div>I have tried your input with QE 5.0, it works fine for me, try it.<br></div>The problem could also arrise from diagonalization method though it is strange... <br></div>Try to set diagonalization = 'cg'. <br>
<div><div>good luck,<br></div><div>Alexander<br><br></div></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">2014-04-29 17:18 GMT+02:00 Pietro Bonfa' <span dir="ltr"><<a href="mailto:pietro.bonfa@fis.unipr.it" target="_blank">pietro.bonfa@fis.unipr.it</a>></span>:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear Dr. Alexander Smogunov,<br>
<br>
the short digression on the non-linear core corrections was intended to<br>
explain why I need DFT+U to correctly reproduce the insulating and<br>
magnetic ground state of La2CuO4.<br>
<br>
I attached a modified version of example 11 (that I called example12).<br>
I set all the parameters to zero just to check if everything is working<br>
properly. I tried also with realistic values for U and J but I obtained<br>
the same errors.<br>
<div class=""><br>
Thanks,<br>
kind regards,<br>
Pietro Bonfa'<br>
<br>
<br>
<br>
</div><div class="">On 04/29/2014 04:13 PM, Alex Smogunov wrote:<br>
> Dear Pietro.<br>
> It should work, in principle, at least for US PPs ... even though<br>
> I do not understand what is the relation between  non-linear core<br>
> corrections and<br>
> the need for DFT+U together with spin-orbit, moreover you set up all the<br>
> parameters<br>
> essentially to zero. Anyway the code should run even in this case. Could<br>
> you provide the complete<br>
> input file, please?<br>
> regards,<br>
> Alexander<br>
><br>
><br>
><br>
> 2014-04-28 10:52 GMT+02:00 Pietro Bonfa' <<a href="mailto:pietro.bonfa@fis.unipr.it">pietro.bonfa@fis.unipr.it</a><br>
</div>> <mailto:<a href="mailto:pietro.bonfa@fis.unipr.it">pietro.bonfa@fis.unipr.it</a>>>:<br>
<div><div class="h5">><br>
>     Dear All,<br>
><br>
>     I would like to compare the total energies of various non-collinear<br>
>     magnetic configurations in La2CuO4 (LCO).<br>
><br>
>     I started working with pseudopotentials from the PSLibrary, the GBRV<br>
>     library and some home brewed norm conserving Cu pseudos. After a short<br>
>     analysis I realized that the experimentally observed insulating and<br>
>     antiferromagnetic ground state is obtained only if I remove non-linear<br>
>     core corrections. This holds true for both my NC pseudos and the Cu<br>
>     pseudo from PSL 1.0.0 (I didn't try GBRV).<br>
><br>
>     Nonetheless, since studying the magnetic ground state of LCO without<br>
>     using non-linear core corrections does not make much sense to me, I<br>
>     tried DFT+U with noncolin and lspinorb flags set to true.<br>
><br>
>     More precisely, I added the following lines to example 11 in the<br>
>     QE-5.1rc1 release<br>
><br>
>         lda_plus_u=.true.<br>
>         lda_plus_u_kind=1<br>
>         Hubbard_U(1)=1.d-10<br>
>         Hubbard_J(1,1)=1.d-10<br>
>         Hubbard_J(2,1)=0.0<br>
><br>
>     Depending on the pseudopotential type (PAW, US or NC) I get errors like:<br>
><br>
>          from cdiaghg : error #        19<br>
>          eigenvectors failed to converge<br>
>     or<br>
>          from cdiaghg : error #       109<br>
>          S matrix not positive definite<br>
><br>
>     So my question is: is lda+u and non-collinear magnetism still not<br>
>     implemented, is this a bug or me doing something utterly wrong?<br>
><br>
>     Thanks,<br>
>     kind regards,<br>
>     Pietro Bonfa'<br>
><br>
>     --<br>
>     Pietro Bonfa' - PhD student<br>
>     Dipartimento di Fisica e Scienze della Terra "Macedonio Melloni"<br>
>     Viale delle Scienze 7A<br>
>     43124 Parma - Italy<br>
>     _______________________________________________<br>
>     Pw_forum mailing list<br>
</div></div>>     <a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a> <mailto:<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>
<div class="HOEnZb"><div class="h5">><br>
><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>
<br>
<br>
--<br>
Pietro Bonfa' - PhD student<br>
Dipartimento di Fisica e Scienze della Terra "Macedonio Melloni"<br>
Viale delle Scienze 7A<br>
43124 Parma - Italy<br>
</div></div><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></blockquote></div><br></div>