<div dir="ltr"><div class="gmail_default" style="font-family:verdana,sans-serif;color:#000000">Quite a few things are going on "Under the hood" in these systems, so let me explain a little.</div><div class="gmail_default" style="font-family:verdana,sans-serif;color:#000000"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;color:#000000">For heavy atoms with large RMTs the initial pseudocharge (PW density inside the RMTs) from dstart is quite bad. In addition, unless you are careful to chose a "good" initial spin state in case.inst, the spins can be quite far from the fixed point solution. While it might appear that the density and energy (:DIS and :ENE) are oscillating, in reality the pseudocharge and the spin/orbital moments are rotating slowly in a N-Dimensional space. Rotations are non-linear, and the mixer (all versions in all DFT codes to my knowledge) use linear expansions. It can therefore take a long time for these rotations tocomplete. What you will often see if you look at the pseudocharge (grep :CPC in 18.1, the pseudocharge was not monited in earlier versions as its role was not known) you will often see that it is slowly "walking" -- similarly the spin moments etc.</div><div class="gmail_default" style="font-family:verdana,sans-serif;color:#000000"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;color:#000000">When you add SOC and/or +U (or -eece) there will be large rotations of the spins and/or spin moments.</div><div class="gmail_default" style="font-family:verdana,sans-serif;color:#000000"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;color:#000000">In an ideal world some transformation would be applied to the Wien2k variables to separate the pseudocharge and convert the rotations to linear variables. I do not know how to do this, and simple ideas such as using the sum and difference of spins don't work. (As expected, since this is a linear transformation.) I am open to (mathematical) suggestions/inspirations. It has to be math not fudge -- the mixer only works when the math is right!</div></div><br><div class="gmail_quote"><div dir="ltr">On Thu, Nov 1, 2018 at 9:34 AM Xavier Rocquefelte <<a href="mailto:xavier.rocquefelte@univ-rennes1.fr">xavier.rocquefelte@univ-rennes1.fr</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">I have exactly the same experience. It is really a matter of cooking to <br>
my point of view and strongly depends on the system and what is <br>
happening at the Fermi level.<br>
<br>
In my case I am using DFT+U+SOC for such systems. The hubbard term <br>
indeed allows to separate the 4f states. However, the SOC is also <br>
essential and leads to large effects in these systems.<br>
In the case of HoPtBi, all the elements requires SOC. In addition it is <br>
an intermetallic compound, meaning that you have both localized and <br>
delocalized electrons!<br>
<br>
Did you try the following procedure:<br>
- DFT+SOC in a self-consistent manner<br>
- DFT+SOC+U on top of the first calculation<br>
<br>
Sometimes the lift of degeneracy due to SOC is sufficient to converge <br>
more efficiently DFT+SOC and then after DFT+SOC+U. And sometimes it is <br>
the other way around;)<br>
<br>
Best wishes<br>
Xavier<br>
<br>
<br>
-- <br>
------------------------<br>
Institut des Sciences Chimiques de Rennes (ISCR)<br>
Univ Rennes - CNRS - UMR6226, France<br>
<a href="https://urldefense.proofpoint.com/v2/url?u=https-3A__iscr.univ-2Drennes1.fr_cti_people_permanent-2Dstaff_rocquefelte-2Dxavier&d=DwICAg&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=V1_Z6_elFWMevjRRzb6Er-1CwhwNLWJ5nHzKWh26Fo0&s=adjkiyhXqoBS6b7SN1a8EitYiDIMwrzO8SLBzLqxJ0M&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.com/v2/url?u=https-3A__iscr.univ-2Drennes1.fr_cti_people_permanent-2Dstaff_rocquefelte-2Dxavier&d=DwICAg&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=V1_Z6_elFWMevjRRzb6Er-1CwhwNLWJ5nHzKWh26Fo0&s=adjkiyhXqoBS6b7SN1a8EitYiDIMwrzO8SLBzLqxJ0M&e=</a><br>
------------------------<br>
ICAMM2019 : VASP Workshop and International Materials Modelling Conf<br>
June 26-July 3 2019, Rennes France<br>
<a href="http://icamm.2019.sciencesconf.org" rel="noreferrer" target="_blank">icamm.2019.sciencesconf.org</a>  : Registration opening soon, see you there!<br>
<br>
_______________________________________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a><br>
<a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__zeus.theochem.tuwien.ac.at_mailman_listinfo_wien&d=DwICAg&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=V1_Z6_elFWMevjRRzb6Er-1CwhwNLWJ5nHzKWh26Fo0&s=bHKBKCqE5r9Ol5EZoCphMsvUqWp9VPpa4qmZjlZ76MY&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.com/v2/url?u=http-3A__zeus.theochem.tuwien.ac.at_mailman_listinfo_wien&d=DwICAg&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=V1_Z6_elFWMevjRRzb6Er-1CwhwNLWJ5nHzKWh26Fo0&s=bHKBKCqE5r9Ol5EZoCphMsvUqWp9VPpa4qmZjlZ76MY&e=</a><br>
SEARCH the MAILING-LIST at:  <a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__www.mail-2Darchive.com_wien-40zeus.theochem.tuwien.ac.at_index.html&d=DwICAg&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=V1_Z6_elFWMevjRRzb6Er-1CwhwNLWJ5nHzKWh26Fo0&s=QFcXH13Qr3d3q8iYQbINXbHncLwDP5UPJWJ9QJbhqIA&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.com/v2/url?u=http-3A__www.mail-2Darchive.com_wien-40zeus.theochem.tuwien.ac.at_index.html&d=DwICAg&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=V1_Z6_elFWMevjRRzb6Er-1CwhwNLWJ5nHzKWh26Fo0&s=QFcXH13Qr3d3q8iYQbINXbHncLwDP5UPJWJ9QJbhqIA&e=</a><br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><span style="font-size:12.8px">Professor Laurence Marks</span><br></div><div dir="ltr"><span style="font-size:12.8px">"Research is to see what everybody else has seen, and to think what nobody else has thought", </span><span style="font-size:12.8px">Albert Szent-Gyorgi</span><br><a href="http://www.numis.northwestern.edu" target="_blank">www.numis.northwestern.edu</a> ; <span style="font-size:12.8px">Corrosion in 4D: </span><a href="http://MURI4D.numis.northwestern.edu" style="font-size:12.8px" target="_blank">MURI4D.numis.northwestern.edu</a><div><span style="font-size:12.8px">Partner of the CFW 100% program for gender equity, </span><a href="http://www.cfw.org/100-percent" style="font-size:12.8px" target="_blank">www.cfw.org/100-percent</a></div><div>Co-Editor, Acta Cryst A</div></div></div></div></div></div></div></div></div></div></div>