<html xmlns:v="urn:schemas-microsoft-com:vml" 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=us-ascii">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<!--[if !mso]><style>v\:* {behavior:url(#default#VML);}
o\:* {behavior:url(#default#VML);}
w\:* {behavior:url(#default#VML);}
.shape {behavior:url(#default#VML);}
</style><![endif]--><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:"Segoe UI Light";
        panose-1:2 11 5 2 4 2 4 2 2 3;}
@font-face
        {font-family:"Segoe UI";
        panose-1:2 11 5 2 4 2 4 2 2 3;}
/* 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
        {mso-style-priority:99;
        margin:0cm;
        margin-bottom:.0001pt;
        font-size:12.0pt;
        font-family:"Times New Roman",serif;}
p.msonormal0, li.msonormal0, div.msonormal0
        {mso-style-name:msonormal;
        margin:0cm;
        margin-bottom:.0001pt;
        font-size:12.0pt;
        font-family:"Times New Roman",serif;}
span.EmailStyle19
        {mso-style-type:personal-compose;
        font-family:"Calibri",sans-serif;
        color:black;}
.MsoChpDefault
        {mso-style-type:export-only;
        font-size:10.0pt;}
@page WordSection1
        {size:612.0pt 792.0pt;
        margin:72.0pt 72.0pt 72.0pt 72.0pt;}
div.WordSection1
        {page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="EN-US" link="blue" vlink="purple">
<div class="WordSection1">
<p><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black">Dear Prof. Blaha<o:p></o:p></span></p>
<p><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black">Thank you.<o:p></o:p></span></p>
<p><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black">In hexagonal structure, LuFeO3 has 6 f.u. per unit cell with 30 atoms.  Once I split Fe (6c Wyckoff) to Fe1 and Fe2, it changes space group to the one with lower symmetry and I am
 trying to even lower symmetry by splitting at least oxygen atoms. The final objective is core hole for EELS calculations, so I do not want multiplicity of atoms. <o:p></o:p></span></p>
<p><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black">The analysis of case.scf shows following results:<o:p></o:p></span></p>
<p><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:red">MMT: 0.00012
<o:p></o:p></span></p>
<p><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:red">MMI:<o:p></o:p></span></p>
<p><span style="font-size:10.0pt;font-family:"Calibri",sans-serif;color:red">case.scf::MMI001: MAGNETIC MOMENT IN SPHERE   1    =   -0.00350</span><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black"><o:p></o:p></span></p>
<div>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Calibri",sans-serif;color:red">case.scf::MMI002: MAGNETIC MOMENT IN SPHERE   2    =   -0.00613</span><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Calibri",sans-serif;color:red">case.scf::MMI003: MAGNETIC MOMENT IN SPHERE   3    =   -0.00768</span><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Calibri",sans-serif;color:red">case.scf::MMI004: MAGNETIC MOMENT IN SPHERE   4    =   -0.00385</span><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Calibri",sans-serif;color:red">case.scf::MMI005: MAGNETIC MOMENT IN SPHERE   5    =    4.02712
</span><span style="font-size:10.0pt;font-family:"Calibri",sans-serif;color:#5B9BD5">(Fe)</span><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Calibri",sans-serif;color:red">case.scf::MMI006: MAGNETIC MOMENT IN SPHERE   6    =   -0.07328</span><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Calibri",sans-serif;color:red">case.scf::MMI007: MAGNETIC MOMENT IN SPHERE   7    =    0.09103</span><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Calibri",sans-serif;color:red">case.scf::MMI008: MAGNETIC MOMENT IN SPHERE   8    =   -0.09647</span><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Calibri",sans-serif;color:red">case.scf::MMI009: MAGNETIC MOMENT IN SPHERE   9    =    0.10131</span><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Calibri",sans-serif;color:red">case.scf::MMI010: MAGNETIC MOMENT IN SPHERE  10    =   -0.27719</span><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Calibri",sans-serif;color:red">case.scf::MMI011: MAGNETIC MOMENT IN SPHERE  11    =    0.26978</span><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Calibri",sans-serif;color:red">case.scf::MMI012: MAGNETIC MOMENT IN SPHERE  12    =   -0.28291</span><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Calibri",sans-serif;color:red">case.scf::MMI013: MAGNETIC MOMENT IN SPHERE  13    =    0.29177</span><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Calibri",sans-serif;color:red">case.scf::MMI014: MAGNETIC MOMENT IN SPHERE  14    =   -4.01690
</span><span style="font-size:10.0pt;font-family:"Calibri",sans-serif;color:#5B9BD5">(Fe)</span><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black"><o:p> </o:p></span></p>
</div>
<p class="MsoNormal"><span style="font-family:"Calibri",sans-serif;color:black">Sincerely<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-family:"Calibri",sans-serif;color:black">Ali<o:p></o:p></span></p>
<p class="MsoNormal" style="margin-bottom:12.0pt"><span style="font-family:"Calibri",sans-serif;color:black"><o:p> </o:p></span></p>
<div>
<div>
<div class="MsoNormal" align="center" style="text-align:center"><span style="font-family:"Calibri",sans-serif;color:black">
<hr size="2" width="98%" align="center">
</span></div>
<div id="x_divRplyFwdMsg">
<p class="MsoNormal"><b><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black">From:</span></b><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black"> Wien <wien-bounces@zeus.theochem.tuwien.ac.at> on behalf of Peter
 Blaha <pblaha@theochem.tuwien.ac.at><br>
<b>Sent:</b> 04 March 2020 17:16<br>
<b>To:</b> wien@zeus.theochem.tuwien.ac.at<br>
<b>Subject:</b> Re: [Wien] LDA+U, spin up and down</span><span style="font-family:"Calibri",sans-serif;color:black">
<o:p></o:p></span></p>
<div>
<p class="MsoNormal"><span style="font-family:"Calibri",sans-serif;color:black"> <o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Calibri",sans-serif;color:black">LuFeO3 has 5 atoms per f.u. with one Fe.<br>
<br>
You have 2 Fe atoms (5 and 14), thus I wonder how you manage to have at <br>
least 14 different atoms per unit cell. Unless the Fe atoms have MULT=2 <br>
or higher, you can have at most 10 atoms/cell. So the question is: is <br>
your stoichiometry correct ?<br>
<br>
Also: What are the :MMI and :MMT ? Are the Lu non-magnetic ? Is :MMT~0<br>
<br>
If spin-up .ne. spin-dn it means that some atoms are not AFM ordered <br>
(even when the Fe are).<br>
<br>
<br>
<br>
Am 04.03.2020 um 15:50 schrieb Ali Baghizhadeh:<br>
> Dear WIEN2k users<br>
> <br>
> I am running spin polarized calculation on LuFeO3 structure, with Fe <br>
> ions having spin up and down in unit cell. I did volume <br>
> optimization/force minimization, and then I used LDA+U to open the <br>
> bandgap. It does work fine for spin up when I plot the density of <br>
> states, but for spin down, it is conductive. I cannot explain this <br>
> behavior, and I do not know how to look for the root/s of the problem. <br>
> Please any comment will be highly appreciated.<br>
> <br>
> Case.inorb file: (Fe1 is atom 5 and Fe2 is atom 14)<br>
> <br>
> 1  2  0                     nmod, natorb, ipr<br>
> <br>
> PRATT  1.0                    BROYD/PRATT, mixing<br>
> <br>
>    5 1 2                          iatom nlorb, lorb<br>
> <br>
>    14 1 2                          iatom nlorb, lorb<br>
> <br>
>    1                              nsic 0..AMF, 1..SIC, 2..HFM<br>
> <br>
>     0.367 0.073        U J (Ry)   Note: you can also use U_eff = U-J and J=0<br>
> <br>
>     0.367 0.073        U J<br>
> <br>
> Sincerely<br>
> <br>
> Ali<br>
> <br>
> <br>
> _______________________________________________<br>
> Wien mailing list<br>
> Wien@zeus.theochem.tuwien.ac.at<br>
> <a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" id="LPlnk740886">
http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
> SEARCH the MAILING-LIST at:  <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" id="LPlnk622964">
http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
> <br>
<br>
-- <br>
--------------------------------------------------------------------------<br>
Peter BLAHA, Inst.f. Materials Chemistry, TU Vienna, A-1060 Vienna<br>
Phone: +43-1-58801-165300             FAX: +43-1-58801-165982<br>
Email: blaha@theochem.tuwien.ac.at    WIEN2k: <a href="http://www.wien2k.at" id="LPlnk108844">
http://www.wien2k.at</a><br>
WWW: <br>
<a href="http://www.imc.tuwien.ac.at/tc_blaha-------------------------------------------------------------------------" id="LPlnk265101">http://www.imc.tuwien.ac.at/tc_blaha-------------------------------------------------------------------------</a>
<o:p></o:p></span></p>
<div style="margin-bottom:15.0pt;overflow:auto" id="LPBorder_GT_15833596784690.7628555230940461">
<table class="MsoNormalTable" border="1" cellspacing="0" cellpadding="0" width="90%" style="width:90.0%;background:white;border-top:dotted #C8C8C8 1.0pt;border-left:none;border-bottom:dotted #C8C8C8 1.0pt;border-right:none">
<tbody>
<tr>
<td valign="top" style="border:none;padding:0cm 0cm 0cm 0cm">
<div id="LPTitle_15833596784680.51465398305518">
<p class="MsoNormal" style="margin-top:15.0pt;mso-line-height-alt:15.75pt"><span style="font-size:16.0pt;font-family:"Segoe UI Light",sans-serif;color:#0078D7"><a href="http://www.imc.tuwien.ac.at/tc_blaha-------------------------------------------------------------------------" target="_blank"><span style="text-decoration:none">IMC
 :&nbsp;Prof. Dr. P. Blaha: Computational Materials Science - Home of WIEN2k</span></a><o:p></o:p></span></p>
</div>
<div style="margin-top:7.5pt;margin-bottom:12.0pt" id="LPMetadata_15833596784690.011727484701072433">
<p class="MsoNormal" style="margin-top:15.0pt;line-height:10.5pt"><span style="font-size:10.5pt;font-family:"Segoe UI",sans-serif;color:#666666">www.imc.tuwien.ac.at<o:p></o:p></span></p>
</div>
<div id="LPDescription_15833596784690.1949988564226628">
<p class="MsoNormal" style="margin-top:15.0pt;line-height:15.0pt"><span style="font-size:10.5pt;font-family:"Segoe UI",sans-serif;color:#666666">Homepage of Institute of Materials Chemistry<o:p></o:p></span></p>
</div>
</td>
</tr>
</tbody>
</table>
</div>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Calibri",sans-serif;color:black"><br>
<br>
<br>
_______________________________________________<br>
Wien mailing list<br>
Wien@zeus.theochem.tuwien.ac.at<br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" id="LPlnk536495">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
SEARCH the MAILING-LIST at:  <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" id="LPlnk650104">
http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><o:p></o:p></span></p>
</div>
</div>
</div>
</body>
</html>