<div dir="ltr"><div><div><div><div><div><div><div>Dear Pietro.<br></div>It should work, in principle, at least for US PPs ... even though<br></div>I do not understand what is the relation between  non-linear core corrections and<br>
</div>the need for DFT+U together with spin-orbit, moreover you set up all the parameters<br></div>essentially to zero. Anyway the code should run even in this case. Could you provide the complete<br></div>input file, please?<br>
</div>regards,<br></div>Alexander<br><div><div><br></div></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">2014-04-28 10:52 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 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>
<span class="HOEnZb"><font color="#888888"><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>
<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></div>