<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none"><!-- p { margin-top: 0px; margin-bottom: 0px; } @font-face {
font-family: Cambria Math;
}
@font-face {
font-family: Calibri;
}
@font-face {
font-family: Verdana;
}
p.MsoNormal, li.MsoNormal, div.MsoNormal { margin: 0cm 0cm 0pt; font-family: "Times New Roman",serif; font-size: 12pt; } a:link, span.MsoHyperlink { color: blue; text-decoration: underline; } a:visited, span.MsoHyperlinkFollowed { color: purple; text-decoration: underline; } p.msonormal0, li.msonormal0, div.msonormal0 { font-family: "Times New Roman",serif; font-size: 12pt; margin-right: 0cm; margin-left: 0cm; } span.EmailStyle18 { color: rgb(31, 73, 125); font-family: "Calibri",sans-serif; } .MsoChpDefault { font-family: "Calibri",sans-serif; }--></style>
</head>
<body dir="ltr" style="font-size:12pt;color:#000000;background-color:#FFFFFF;font-family:Calibri,Arial,Helvetica,sans-serif;">
<p>No, it's not related to option AFM specifically. If I understood correctly, your system is AFM (atoms 5 and 14 are Fe with opposite signs of the magnetic moment). If yes, then the U and J in case.inorb for both atoms should be the same:</p>
<p class="MsoNormal"><span style="color: rgb(31, 73, 125); font-family: "Calibri",sans-serif; font-size: 11pt;">1 2 0 nmod, natorb, ipr</span></p>
<p class="MsoNormal"><span style="color: rgb(31, 73, 125); font-family: "Calibri",sans-serif; font-size: 11pt;">PRATT 1.0 BROYD/PRATT, mixing</span></p>
<p class="MsoNormal"><span style="color: rgb(31, 73, 125); font-family: "Calibri",sans-serif; font-size: 11pt;">5 1 2 iatom nlorb, lorb</span></p>
<p class="MsoNormal"><span style="color: rgb(31, 73, 125); font-family: "Calibri",sans-serif; font-size: 11pt;">14 1 2 iatom nlorb, lorb</span></p>
<p class="MsoNormal"><span style="color: rgb(31, 73, 125); font-family: "Calibri",sans-serif; font-size: 11pt;">1 </span></p>
<p><font color="#1f497d" face="Calibri">0.367 0.073 U J (Ry) Note: you can also use U_eff = U-J and J=0 (for atom 5)</font></p>
<p><span><font color="#1f497d" face="Calibri">0.367 0.073 U J (Ry) Note: you can also use U_eff = U-J and J=0<font color="#1f497d"> (for atom 14)</font></font></span></p>
<p> </p>
<div style="color: rgb(0, 0, 0); font-family: Calibri,Arial,Helvetica,sans-serif; font-size: 12pt; background-color: rgb(255, 255, 255);" dir="ltr">
<hr tabindex="-1" style="width: 98%; display: inline-block;">
<div id="divRplyFwdMsg" dir="ltr"><font color="#000000" face="Calibri, sans-serif" style="font-size: 11pt;"><b>From:</b> Wien <wien-bounces@zeus.theochem.tuwien.ac.at> on behalf of Ali Baghizhadeh <ali.baghizhadeh@ua.pt><br>
<b>Sent:</b> Thursday, February 20, 2020 5:16 PM<br>
<b>To:</b> A Mailing list for WIEN2k users<br>
<b>Subject:</b> Re: [Wien] Large band gap in LDA+U</font>
<div> </div>
</div>
<div>
<div style="margin: 0px; padding: 0px; color: black; font-family: sans-serif; font-size: 11pt; direction: ltr;" dir="auto">
Dear Fabian<br>
</div>
<div style="margin: 0px; padding: 0px; color: black; font-family: sans-serif; font-size: 11pt; direction: ltr;" dir="auto">
I did not use AFM option. I meant that my magnetic system is AFM, by simply considering one Fe spin up and another down. I guess what you commented is for AFM option. I will cross check it again in manual. Thank you very much.<br>
</div>
<div style="margin: 0px; padding: 0px; color: black; font-family: sans-serif; font-size: 11pt; direction: ltr;" dir="auto">
Ali<br>
<br>
</div>
<div style="margin: 0px; padding: 0px; color: black; font-family: sans-serif; font-size: 11pt; direction: ltr;" dir="auto">
<span id="OutlookSignature">
<div style="margin: 0px; padding: 0px; color: black; font-family: sans-serif; font-size: 11pt; direction: ltr;" dir="auto">
Get <a href="https://aka.ms/ghei36">Outlook for Android</a></div>
</span><br>
</div>
<hr tabindex="-1" style="width: 98%; display: inline-block;">
<div id="divRplyFwdMsg" dir="ltr"><font color="#000000" face="Calibri, sans-serif" style="font-size: 11pt;"><b>From:</b> Wien <wien-bounces@zeus.theochem.tuwien.ac.at> on behalf of Tran, Fabien <fabien.tran@tuwien.ac.at><br>
<b>Sent:</b> Thursday, February 20, 2020 5:08:33 PM<br>
<b>To:</b> A Mailing list for WIEN2k users <wien@zeus.theochem.tuwien.ac.at><br>
<b>Subject:</b> Re: [Wien] Large band gap in LDA+U</font>
<div> </div>
</div>
<div>
<p>You mentioned that you are doing AFM calculations, but case.inorb and case.indmc do not correspond to that. In these files, there should be two lines (one for Fe (spin-up) and the other for Fe(spin-down)), as well as the corresponding two lines for the values
of U and J. Look at the example of AFM NiO in this file:</p>
<p>$WIENROOT/SRC_templates/case.inorb<br>
</p>
<p><br>
</p>
<div style="color: rgb(33, 33, 33);">
<hr tabindex="-1" style="width: 98%; display: inline-block;">
<div id="divRplyFwdMsg" dir="ltr"><font color="#000000" face="Calibri, sans-serif" style="font-size: 11pt;"><b>From:</b> Wien <wien-bounces@zeus.theochem.tuwien.ac.at> on behalf of Ali Baghizhadeh <ali.baghizhadeh@ua.pt><br>
<b>Sent:</b> Thursday, February 20, 2020 4:37 PM<br>
<b>To:</b> A Mailing list for WIEN2k users<br>
<b>Subject:</b> Re: [Wien] Large band gap in LDA+U</font>
<div> </div>
</div>
<div>
<div class="WordSection1">
<p class="MsoNormal"><span style="color: rgb(31, 73, 125); font-family: "Calibri",sans-serif; font-size: 11pt;">Dear Prof. Laurence</span></p>
<p class="MsoNormal"><span style="color: rgb(31, 73, 125); font-family: "Calibri",sans-serif; font-size: 11pt;">Thank you.</span></p>
<p class="MsoNormal"><span style="color: rgb(31, 73, 125); font-family: "Calibri",sans-serif; font-size: 11pt;">My lattice is (Lu0.5Sc0.5)FeO3, for LuFeO3 everything works fine (LDA+U).
</span></p>
<p class="MsoNormal"><span style="color: rgb(31, 73, 125); font-family: "Calibri",sans-serif; font-size: 11pt;">My case.inorb is as following, with atom 5 and 14 as Fe:</span></p>
<p class="MsoNormal"><span style="color: rgb(31, 73, 125); font-family: "Calibri",sans-serif; font-size: 11pt;">1 2 0 nmod, natorb, ipr</span></p>
<p class="MsoNormal"><span style="color: rgb(31, 73, 125); font-family: "Calibri",sans-serif; font-size: 11pt;">PRATT 1.0 BROYD/PRATT, mixing</span></p>
<p class="MsoNormal"><span style="color: rgb(31, 73, 125); font-family: "Calibri",sans-serif; font-size: 11pt;"> 5 1 2 iatom nlorb, lorb</span></p>
<p class="MsoNormal"><span style="color: rgb(31, 73, 125); font-family: "Calibri",sans-serif; font-size: 11pt;"> 14 1 2 iatom nlorb, lorb</span></p>
<p class="MsoNormal"><span style="color: rgb(31, 73, 125); font-family: "Calibri",sans-serif; font-size: 11pt;"> 1 nsic 0..AMF, 1..SIC, 2..HFM</span></p>
<p class="MsoNormal"><span style="color: rgb(31, 73, 125); font-family: "Calibri",sans-serif; font-size: 11pt;"> 0.367 0.073 U J (Ry) Note: you can also use U_eff = U-J and J=0</span></p>
<p class="MsoNormal"><span style="color: rgb(31, 73, 125); font-family: "Calibri",sans-serif; font-size: 11pt;"> 5.00 1.00 U J</span></p>
<p class="MsoNormal"><span style="color: rgb(31, 73, 125); font-family: "Calibri",sans-serif; font-size: 11pt;"> </span></p>
<p class="MsoNormal"><span style="color: rgb(31, 73, 125); font-family: "Calibri",sans-serif; font-size: 11pt;">And case.indmc</span></p>
<p class="MsoNormal"><span style="color: rgb(31, 73, 125); font-family: "Calibri",sans-serif; font-size: 11pt;">-12. Emin cutoff energy</span></p>
<p class="MsoNormal"><span style="color: rgb(31, 73, 125); font-family: "Calibri",sans-serif; font-size: 11pt;">2 number of atoms for which density matrix is calculated</span></p>
<p class="MsoNormal"><span style="color: rgb(31, 73, 125); font-family: "Calibri",sans-serif; font-size: 11pt;">5 1 2 index of 1st atom, number of L's, L1</span></p>
<p class="MsoNormal"><span style="color: rgb(31, 73, 125); font-family: "Calibri",sans-serif; font-size: 11pt;">14 1 2 ditto for 2nd atom, repeat NATOM times</span></p>
<p class="MsoNormal"><span style="color: rgb(31, 73, 125); font-family: "Calibri",sans-serif; font-size: 11pt;">0 0 r-index, (l,s)index
</span></p>
<p class="MsoNormal"><span style="color: rgb(31, 73, 125); font-family: "Calibri",sans-serif; font-size: 11pt;"> </span></p>
<p class="MsoNormal"><span style="color: rgb(31, 73, 125); font-family: "Calibri",sans-serif; font-size: 11pt;">Best regards</span></p>
<p class="MsoNormal"><span style="color: rgb(31, 73, 125); font-family: "Calibri",sans-serif; font-size: 11pt;">Ali</span></p>
<p class="MsoNormal"><span style="color: rgb(31, 73, 125); font-family: "Calibri",sans-serif; font-size: 11pt;"> </span></p>
<p class="MsoNormal"><b><span style="font-family: "Calibri",sans-serif; font-size: 11pt;">From:</span></b><span style="font-family: "Calibri",sans-serif; font-size: 11pt;"> Wien [mailto:wien-bounces@zeus.theochem.tuwien.ac.at]
<b>On Behalf Of </b>Laurence Marks<br>
<b>Sent:</b> Thursday, 20 February 2020 16:31<br>
<b>To:</b> A Mailing list for WIEN2k users <wien@zeus.theochem.tuwien.ac.at><br>
<b>Subject:</b> Re: [Wien] Large band gap in LDA+U</span></p>
<p class="MsoNormal"> </p>
<div>
<div>
<p class="MsoNormal"><span style="color: black; font-family: "Verdana",sans-serif;">For LuFeO3 you need to do something for the Lu (ScFeO3 is simpler).</span></p>
</div>
<div>
<p class="MsoNormal"><span style="color: black; font-family: "Verdana",sans-serif;"> </span></p>
</div>
<div>
<p class="MsoNormal"><span style="color: black; font-family: "Verdana",sans-serif;">Beyond this, I suspect you setup the LDA+U incorrectly. However, this is a guess.</span></p>
</div>
</div>
<p class="MsoNormal"> </p>
<div>
<div>
<p class="MsoNormal">On Thu, Feb 20, 2020 at 9:25 AM Ali Baghizhadeh <<a href="mailto:ali.baghizhadeh@ua.pt">ali.baghizhadeh@ua.pt</a>> wrote:</p>
</div>
<blockquote style="border-width: medium medium medium 1pt; border-style: none none none solid; border-color: currentColor currentColor currentColor rgb(204, 204, 204); padding: 0cm 0cm 0cm 6pt; margin-right: 0cm; margin-left: 4.8pt;">
<div>
<div>
<p class="MsoNormal">Dear Wien2k users</p>
<p class="MsoNormal">I am trying to calculate band gap of h-(Lu,Sc)FeO3 oxide using WieN2k 19, with PBE potential. To do so, I used spin polarized AFM calculation by breaking symmetry of lattice to have opposite spin directions for Fe. Then after convergence
of scf, I did optimize atomic positions. In PDOS plot, it shows metallic behavior, while experimentally, I expect something like 1.6 eV band gap. Then I used LDA+U, from 1 ev to 6 ev, the PDOS shows band gap of 3.2 eV no matter which U chooses for spin up
Fe, and again metallic phase for spin down Fe. I do not know how to resolve my problem, and what parameters or steps might be wrong and I have to check? Any comment will be really appreciated.</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">Sincerely</p>
<p class="MsoNormal">Ali Baghi zadeh</p>
<p class="MsoNormal"> </p>
</div>
</div>
<p class="MsoNormal">_______________________________________________<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=hq5MSDWdxJEPwQ5co4rYk2RDNIxVXkxIsYp-77rIXZA&s=8IwEoi8Urh-vGGjLg2nIG8ftikyK1H_UJWqM0g9oW-Q&e=" 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=hq5MSDWdxJEPwQ5co4rYk2RDNIxVXkxIsYp-77rIXZA&s=8IwEoi8Urh-vGGjLg2nIG8ftikyK1H_UJWqM0g9oW-Q&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=hq5MSDWdxJEPwQ5co4rYk2RDNIxVXkxIsYp-77rIXZA&s=8Y-51p-BNaqEEKL6yAOGdggmPb7YgHtc-41fmRKFLeQ&e=" 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=hq5MSDWdxJEPwQ5co4rYk2RDNIxVXkxIsYp-77rIXZA&s=8Y-51p-BNaqEEKL6yAOGdggmPb7YgHtc-41fmRKFLeQ&e=</a>
</p>
</blockquote>
</div>
<p class="MsoNormal"><br clear="all">
</p>
<div>
<p class="MsoNormal"> </p>
</div>
<p class="MsoNormal">-- </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></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</p>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</body>
</html>