<html xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:Verdana;
panose-1:2 11 6 4 3 5 4 4 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman",serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
p.msonormal0, li.msonormal0, div.msonormal0
{mso-style-name:msonormal;
mso-margin-top-alt:auto;
margin-right:0cm;
mso-margin-bottom-alt:auto;
margin-left:0cm;
font-size:12.0pt;
font-family:"Times New Roman",serif;}
span.EmailStyle18
{mso-style-type:personal;
font-family:"Calibri",sans-serif;
color:#1F497D;}
span.EmailStyle19
{mso-style-type:personal-compose;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri",sans-serif;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:72.0pt 72.0pt 72.0pt 72.0pt;}
div.WordSection1
{page:WordSection1;}
--></style>
</head>
<body lang="EN-US" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D">Dear Prof. Marc<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D">Thank you for your comments. For comment 3, actually I did not pay attention, and for comment 4, I will try other AFM ordering as I have 4 types for.
<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D">Thank you again.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D">Ali<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><b><span style="font-size:11.0pt;font-family:"Calibri",sans-serif">From:</span></b><span style="font-size:11.0pt;font-family:"Calibri",sans-serif"> Wien [mailto:wien-bounces@zeus.theochem.tuwien.ac.at]
<b>On Behalf Of </b>Laurence Marks<br>
<b>Sent:</b> Wednesday, 4 March 2020 17:10<br>
<b>To:</b> A Mailing list for WIEN2k users <wien@zeus.theochem.tuwien.ac.at><br>
<b>Subject:</b> Re: [Wien] LDA+U, spin up and down<o:p></o:p></span></p>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<div>
<p class="MsoNormal"><span style="font-family:"Verdana",sans-serif;color:black">Let me add to Fabien's statement:<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Verdana",sans-serif;color:#222222">"There may be two reasons:<br>
1) You started the calculation with an electron density that is not AFM.<br>
2) The Fe atoms (5 and 14) are inequivalent also crystallographically.</span><span style="font-family:"Verdana",sans-serif;color:black"> "<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Verdana",sans-serif;color:black"><o:p> </o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Verdana",sans-serif;color:black">3) The Lu 4f electrons are not well below the Fermi energy.<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-family:"Verdana",sans-serif;color:black">4) You have chosen the wrong AFM ordering. Assuming that you started with the Pbnm perovskite there are many ways to have the AFM ordering. If you pick the wrong one then the results
will be wrong -- you need to search over the subgroups.<o:p></o:p></span></p>
</div>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<div>
<p class="MsoNormal">On Wed, Mar 4, 2020 at 9:49 AM Tran, Fabien <<a href="mailto:fabien.tran@tuwien.ac.at">fabien.tran@tuwien.ac.at</a>> wrote:<o:p></o:p></p>
</div>
<blockquote style="border:none;border-left:solid #CCCCCC 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-right:0cm">
<p class="MsoNormal">Ok but what are the values of :GAP in case.scf? More precisely, what are the last two values (i.e., spin up and spin down for the last SCF iteration) of<br>
:GAP (this spin):<br>
in case.scf.<br>
<br>
________________________________________<br>
From: Wien <<a href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at" target="_blank">wien-bounces@zeus.theochem.tuwien.ac.at</a>> on behalf of Ali Baghizhadeh <<a href="mailto:ali.baghizhadeh@ua.pt" target="_blank">ali.baghizhadeh@ua.pt</a>><br>
Sent: Wednesday, March 4, 2020 4:41 PM<br>
To: A Mailing list for WIEN2k users<br>
Subject: Re: [Wien] LDA+U, spin up and down<br>
<br>
Dear Fabian<br>
The magnetic moment for Fe5 is 4.02712 and for Fe14 is -4.01690 and gap is 1.523 eV, very close to experimental value.<br>
Regards<br>
Ali<br>
<br>
-----Original Message-----<br>
From: Wien [mailto:<a href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at" target="_blank">wien-bounces@zeus.theochem.tuwien.ac.at</a>] On Behalf Of Tran, Fabien<br>
Sent: Wednesday, 4 March 2020 16:35<br>
To: A Mailing list for WIEN2k users <<a href="mailto:wien@zeus.theochem.tuwien.ac.at" target="_blank">wien@zeus.theochem.tuwien.ac.at</a>><br>
Subject: Re: [Wien] LDA+U, spin up and down<br>
<br>
What are the magnetic moments on the atoms (grep for :MMI in case.scf)? Do atoms 5 and 14 have moments with opposite sign?<br>
What is the value of the band gap of spin up (grep for :GAP in case.scf)?<br>
<br>
________________________________________<br>
From: Wien <<a href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at" target="_blank">wien-bounces@zeus.theochem.tuwien.ac.at</a>> on behalf of Ali Baghizhadeh <<a href="mailto:ali.baghizhadeh@ua.pt" target="_blank">ali.baghizhadeh@ua.pt</a>><br>
Sent: Wednesday, March 4, 2020 4:30 PM<br>
To: A Mailing list for WIEN2k users<br>
Subject: Re: [Wien] LDA+U, spin up and down<br>
<br>
Dear Fabian<br>
Thank you Fabian. I did not use AFM option, instead I used spin polarized option. For this, I did assig the spin up and down in instgen_lapw and following steps as you mentioned. Before conducting force minimization, the Fe1 and Fe2 are crystallographically
equivalent. But after force minimization, there is slight change in some coordinates of both Fe, no longer are equivalent. I did not consider it as negative point, but apparently it is affecting my calculation.<br>
Thank you again.<br>
Ali<br>
<br>
-----Original Message-----<br>
From: Wien [mailto:<a href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at" target="_blank">wien-bounces@zeus.theochem.tuwien.ac.at</a>] On Behalf Of Tran, Fabien<br>
Sent: Wednesday, 4 March 2020 16:06<br>
To: A Mailing list for WIEN2k users <<a href="mailto:wien@zeus.theochem.tuwien.ac.at" target="_blank">wien@zeus.theochem.tuwien.ac.at</a>><br>
Subject: Re: [Wien] LDA+U, spin up and down<br>
<br>
Hi,<br>
From your previous emails, I understood that you want to do a AFM calculations on LuFeO3.<br>
If the gaps of both spins are not the same, then it means that the resulting magnetic phase is not AFM, but something else (ferromagnetic or ferrimagnetic). There may be two reasons:<br>
1) You started the calculation with an electron density that is not AFM.<br>
2) The Fe atoms (5 and 14) are inequivalent also crystallographically.<br>
<br>
To produce a starting AFM electron density execute:<br>
instgen_lapw (and select u,d or n appropriately) x lstart x dstart x dstart -up x dstart -dn<br>
<br>
FT<br>
<br>
From: Wien <<a href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at" target="_blank">wien-bounces@zeus.theochem.tuwien.ac.at</a>> on behalf of Ali Baghizhadeh <<a href="mailto:ali.baghizhadeh@ua.pt" target="_blank">ali.baghizhadeh@ua.pt</a>><br>
Sent: Wednesday, March 4, 2020 3:50 PM<br>
To: A Mailing list for WIEN2k users<br>
Subject: [Wien] LDA+U, spin up and down<br>
<br>
<br>
Dear WIEN2k users<br>
I am running spin polarized calculation on LuFeO3 structure, with Fe ions having spin up and down in unit cell. I did volume optimization/force minimization, and then I used LDA+U to open the bandgap. It does work fine for spin up when I plot the density of
states, but for spin down, it is conductive. I cannot explain this behavior, and I do not know how to look for the root/s of the problem. Please any comment will be highly appreciated.<br>
<br>
Case.inorb file: (Fe1 is atom 5 and Fe2 is atom 14)<br>
1 2 0 nmod, natorb, ipr PRATT 1.0 BROYD/PRATT, mixing<br>
5 1 2 iatom nlorb, lorb<br>
14 1 2 iatom nlorb, lorb<br>
1 nsic 0..AMF, 1..SIC, 2..HFM<br>
0.367 0.073 U J (Ry) Note: you can also use U_eff = U-J and J=0<br>
0.367 0.073 U J<br>
<br>
Sincerely<br>
Ali<br>
<br>
<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.com/v3/__http:/zeus.theochem.tuwien.ac.at/mailman/listinfo/wien__;!!Dq0X2DkFhyF93HkjWTBQKhk!EGe9noI1LB2I7rFh231p1z10mr9orJDxf0byHoAksgzoyaF_G51px7k3KV9QpaZ49dLZ0w$" target="_blank">https://urldefense.com/v3/__http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien__;!!Dq0X2DkFhyF93HkjWTBQKhk!EGe9noI1LB2I7rFh231p1z10mr9orJDxf0byHoAksgzoyaF_G51px7k3KV9QpaZ49dLZ0w$</a>
<br>
SEARCH the MAILING-LIST at: <a href="https://urldefense.com/v3/__http:/www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html__;!!Dq0X2DkFhyF93HkjWTBQKhk!EGe9noI1LB2I7rFh231p1z10mr9orJDxf0byHoAksgzoyaF_G51px7k3KV9Qpabnw4N9Ew$" target="_blank">
https://urldefense.com/v3/__http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html__;!!Dq0X2DkFhyF93HkjWTBQKhk!EGe9noI1LB2I7rFh231p1z10mr9orJDxf0byHoAksgzoyaF_G51px7k3KV9Qpabnw4N9Ew$</a>
<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.com/v3/__http:/zeus.theochem.tuwien.ac.at/mailman/listinfo/wien__;!!Dq0X2DkFhyF93HkjWTBQKhk!EGe9noI1LB2I7rFh231p1z10mr9orJDxf0byHoAksgzoyaF_G51px7k3KV9QpaZ49dLZ0w$" target="_blank">https://urldefense.com/v3/__http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien__;!!Dq0X2DkFhyF93HkjWTBQKhk!EGe9noI1LB2I7rFh231p1z10mr9orJDxf0byHoAksgzoyaF_G51px7k3KV9QpaZ49dLZ0w$</a>
<br>
SEARCH the MAILING-LIST at: <a href="https://urldefense.com/v3/__http:/www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html__;!!Dq0X2DkFhyF93HkjWTBQKhk!EGe9noI1LB2I7rFh231p1z10mr9orJDxf0byHoAksgzoyaF_G51px7k3KV9Qpabnw4N9Ew$" target="_blank">
https://urldefense.com/v3/__http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html__;!!Dq0X2DkFhyF93HkjWTBQKhk!EGe9noI1LB2I7rFh231p1z10mr9orJDxf0byHoAksgzoyaF_G51px7k3KV9Qpabnw4N9Ew$</a>
<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.com/v3/__http:/zeus.theochem.tuwien.ac.at/mailman/listinfo/wien__;!!Dq0X2DkFhyF93HkjWTBQKhk!EGe9noI1LB2I7rFh231p1z10mr9orJDxf0byHoAksgzoyaF_G51px7k3KV9QpaZ49dLZ0w$" target="_blank">https://urldefense.com/v3/__http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien__;!!Dq0X2DkFhyF93HkjWTBQKhk!EGe9noI1LB2I7rFh231p1z10mr9orJDxf0byHoAksgzoyaF_G51px7k3KV9QpaZ49dLZ0w$</a>
<br>
SEARCH the MAILING-LIST at: <a href="https://urldefense.com/v3/__http:/www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html__;!!Dq0X2DkFhyF93HkjWTBQKhk!EGe9noI1LB2I7rFh231p1z10mr9orJDxf0byHoAksgzoyaF_G51px7k3KV9Qpabnw4N9Ew$" target="_blank">
https://urldefense.com/v3/__http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html__;!!Dq0X2DkFhyF93HkjWTBQKhk!EGe9noI1LB2I7rFh231p1z10mr9orJDxf0byHoAksgzoyaF_G51px7k3KV9Qpabnw4N9Ew$</a>
<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.com/v3/__http:/zeus.theochem.tuwien.ac.at/mailman/listinfo/wien__;!!Dq0X2DkFhyF93HkjWTBQKhk!EGe9noI1LB2I7rFh231p1z10mr9orJDxf0byHoAksgzoyaF_G51px7k3KV9QpaZ49dLZ0w$" target="_blank">https://urldefense.com/v3/__http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien__;!!Dq0X2DkFhyF93HkjWTBQKhk!EGe9noI1LB2I7rFh231p1z10mr9orJDxf0byHoAksgzoyaF_G51px7k3KV9QpaZ49dLZ0w$</a>
<br>
SEARCH the MAILING-LIST at: <a href="https://urldefense.com/v3/__http:/www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html__;!!Dq0X2DkFhyF93HkjWTBQKhk!EGe9noI1LB2I7rFh231p1z10mr9orJDxf0byHoAksgzoyaF_G51px7k3KV9Qpabnw4N9Ew$" target="_blank">
https://urldefense.com/v3/__http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html__;!!Dq0X2DkFhyF93HkjWTBQKhk!EGe9noI1LB2I7rFh231p1z10mr9orJDxf0byHoAksgzoyaF_G51px7k3KV9Qpabnw4N9Ew$</a>
<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.com/v3/__http:/zeus.theochem.tuwien.ac.at/mailman/listinfo/wien__;!!Dq0X2DkFhyF93HkjWTBQKhk!EGe9noI1LB2I7rFh231p1z10mr9orJDxf0byHoAksgzoyaF_G51px7k3KV9QpaZ49dLZ0w$" target="_blank">https://urldefense.com/v3/__http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien__;!!Dq0X2DkFhyF93HkjWTBQKhk!EGe9noI1LB2I7rFh231p1z10mr9orJDxf0byHoAksgzoyaF_G51px7k3KV9QpaZ49dLZ0w$</a>
<br>
SEARCH the MAILING-LIST at: <a href="https://urldefense.com/v3/__http:/www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html__;!!Dq0X2DkFhyF93HkjWTBQKhk!EGe9noI1LB2I7rFh231p1z10mr9orJDxf0byHoAksgzoyaF_G51px7k3KV9Qpabnw4N9Ew$" target="_blank">
https://urldefense.com/v3/__http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html__;!!Dq0X2DkFhyF93HkjWTBQKhk!EGe9noI1LB2I7rFh231p1z10mr9orJDxf0byHoAksgzoyaF_G51px7k3KV9Qpabnw4N9Ew$</a>
<o:p></o:p></p>
</blockquote>
</div>
<p class="MsoNormal"><br clear="all">
<o:p></o:p></p>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<p class="MsoNormal">-- <o:p></o:p></p>
<div>
<div>
<p class="MsoNormal">Professor Laurence Marks<br>
Department of Materials Science and Engineering<br>
Northwestern University<br>
<a href="http://www.numis.northwestern.edu/" target="_blank">www.numis.northwestern.edu</a><o:p></o:p></p>
<div>
<p class="MsoNormal">Corrosion in 4D: <a href="http://www.numis.northwestern.edu/MURI" target="_blank">
www.numis.northwestern.edu/MURI</a><br>
Co-Editor, Acta Cryst A<br>
"Research is to see what everybody else has seen, and to think what nobody else has thought"<br>
Albert Szent-Gyorgi<o:p></o:p></p>
</div>
</div>
</div>
</div>
</body>
</html>