[Wien] Problems with FM+U+SOC with natorb=1 nlorb=2

[Laurence Marks]

Be aware that +U for 4f in my experience shoves them out of the way,
whereas -eece does not so is considerably better. Depending on your ARPES
energies this can be a significant issues. (The tabulated cross sections
for 4f are really bad, but that is another problem.)

N.B., also, be aware that with large RMTs there are issues with the tails
of adjacent atom states (e.g. O 2p) mimicking d/f states inside a cation
RMT.

---
Prof Laurence Marks
"Research is to see what everyone else has seen, and to think what nobody
else has thought", Albert Szent-Gyorgi
www.numis.northwestern.edu

On Fri, Jul 24, 2020, 11:51 pluto <pluto at physics.ucdavis.edu> wrote:

> Dear Prof. Marks, Prof. Blaha, dear All,
>
> Thank you for the rapid answer! I appreciate!
>
> I think this resolves the issue. For practical reasons I will just use
> +U only for 4f levels.
>
> I agree with Prof. Marks who has reservations regarding the GGA+U for
> 4f. However, for practical reasons of checking the trends of band
> dispersions in ARPES it seems the only quick solution. It's great that
> it is included despite all the deficiencies!
>
> Best,
> Lukasz
>
>
>
>
>
> On 7/24/2020 6:34 PM, Peter Blaha wrote:
> > I do NOT officially support 2 l values / atom for GGA+U. It may/may not
> > work in various programs.
> >
> > The reason is simply, that in WIEN2k the U is applied only inside the
> > sphere and the V_orb potential is calculated from the occupations of an
> > orbital, varying from 0 to 1.
> > I.e. if n=1 (an orbital is fully occupied), it gets a potential which
> > shifts this state by U/2 down, on the other hand, if a particular
> > orbital is empty, it will be shifted up by U/2.
> >
> > The problems come with relatively "delocalized" orbitals like the 5d
> > states of a 4f element.
> >
> > Check out the "fractional charge" within the sphere as printed in
> > case.outputst.
> >   4F*     -0.590624     -0.590624  2.50  2.50    0.9844  F
> >   5D*     -0.208649     -0.208649  0.50  0.50    0.4864  F
> >
> > You can find that Sm-4f is nearly completely confined in the sphere
> > (0.9844) and applicationsof U is "well defined", but one 5d electron
> > has less than 0.5 e inside the sphere !!!
> > Depending on the rare earth and on RMT it means that even a fully
> > occupied 5d state does NOT get shifted down due to GGA+U, but basically
> > all d-orbitals will be shifted up.
> >
> >
> >
> > Am 24.07.2020 um 18:03 schrieb pluto:
> >> Dear All,
> >>
> >> I am trying to calculate a simple 4f compound with GGA+U,
> >> ferromagnetic exchange (FM) and SOC. FM calculation without SOC using
> >> the literature U settings looks reasonable, I am getting a gapped band
> >> structure similar to anyone else.
> >>
> >> But running SCF with GGA+U and SOC removes the gap, it also seems to
> >> change the exchange splitting of 4f, which basically tells me that U
> >> on 4f is not in effect anymore. Actually the problem happens with this
> >> inorb file, when natorb=1 and nlorb=2:
> >>
> >>    1  1  0                     nmod, natorb, ipr
> >> PRATT  1.0                    BROYD/PRATT, mixing
> >>    1 2 2 3                          iatom nlorb, lorb
> >>    1                              nsic 0..AMF, 1..SIC, 2..HFM
> >>     0.25 0.00        U J           U on d levels
> >>     0.42 0.00        U J           U on f levels
> >>
> >> When nlorb is reduced to 1, and U is used only for 4f levels things
> >> start to look "normal" also with SOC. So, it seems SOC cannot handle
> >> nlorb>1 (perhaps just a bug).
> >>
> >> I am using couple of years old version of WIEN2k, perhaps this has
> >> been fixed in the meantime. Otherwise, any comment would be
> >> appreciated!
> >>
> >> Best,
> >> Lukasz
> >>
> _______________________________________________
> Wien mailing list
> Wien at zeus.theochem.tuwien.ac.at
>
> https://urldefense.com/v3/__http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien__;!!Dq0X2DkFhyF93HkjWTBQKhk!F91fPf2jikDWIIB4vyHwju5DwqJSDcqrBeFomZcKf68DtGoOD8p1NkeKAypznfuZiXcK-g$
> SEARCH the MAILING-LIST at:
> https://urldefense.com/v3/__http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html__;!!Dq0X2DkFhyF93HkjWTBQKhk!F91fPf2jikDWIIB4vyHwju5DwqJSDcqrBeFomZcKf68DtGoOD8p1NkeKAypznfvP6tte3g$
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://zeus.theochem.tuwien.ac.at/pipermail/wien/attachments/20200724/1b0b85d9/attachment.html>