Dear wien2k user,<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
I want to reproduce DOS&nbsp; of FeCr2S4 which is
published&nbsp;&nbsp;&nbsp;&nbsp; <span style="font-weight: bold;">Phys. Rev. B 59, 10018 - 10024 (1999)</span><br style="font-weight: bold;">
<span style="font-weight: bold;">[Issue 15 – April 1999]&nbsp;&nbsp; &nbsp;&nbsp; </span>&nbsp; <a href="http://prola.aps.org/abstract/PRB/v59/i15/p10018_1">http://prola.aps.org/abstract/PRB/v59/i15/p10018_1</a> <br>
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; The LSDA calculation
yeilds the ferrimagnetic half metallic states for FeCr2S4 which I get
both for LMTO and LAPW calculation.<br>
Using LSDA+U I get insulator for LMTO method&nbsp; but LAPW method it
remain half metal&nbsp; .I also used gga+u . I tried with different U
and J parameter&nbsp; for&nbsp; LAPW&nbsp; method.<br>
&nbsp; <br>
&nbsp;&nbsp;&nbsp; The necessary files which i used for LAPW calculation are <br>
<br>
&nbsp;FeCr2S4.struct<br>
F&nbsp;&nbsp; LATTICE,NONEQUIV.ATOMS:&nbsp; 3227_Fd-3m<br>
MODE OF CALC=RELA unit=ang<br>
&nbsp;18.876732 18.876732 18.876732 90.000000 90.000000 90.000000<br>
ATOM&nbsp;&nbsp; 1: X=0.37500000 Y=0.37500000 Z=0.37500000<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; MULT= 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ISPLIT= 2<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1: X=0.62500000 Y=0.62500000 Z=0.62500000<br>
Fe&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; NPT=&nbsp; 781&nbsp;
R0=0.00010000 RMT=&nbsp;&nbsp;&nbsp; 2.3000&nbsp;&nbsp; Z: 26.0<br>
LOCAL ROT MATRIX:&nbsp;&nbsp;&nbsp; 1.0000000 0.0000000 0.0000000<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
0.0000000 1.0000000 0.0000000<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
0.0000000 0.0000000 1.0000000<br>
ATOM&nbsp; -2: X=0.00000000 Y=0.00000000 Z=0.00000000<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; MULT= 4&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ISPLIT= 4<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -2: X=0.00000000 Y=0.25000000 Z=0.25000000<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -2: X=0.25000000 Y=0.25000000 Z=0.00000000<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -2: X=0.25000000 Y=0.00000000 Z=0.25000000<br>
Cr&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; NPT=&nbsp; 781&nbsp;
R0=0.00010000 RMT=&nbsp;&nbsp;&nbsp; 2.3000&nbsp;&nbsp; Z: 24.0<br>
LOCAL ROT MATRIX:&nbsp;&nbsp;&nbsp; 0.4082483-0.7071068 0.5773503<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
0.4082483 0.7071068 0.5773503<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
-0.8164966 0.0000000 0.5773503<br>
ATOM&nbsp; -3: X=0.75900000 Y=0.75900000 Z=0.75900000<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; MULT= 8&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ISPLIT= 4<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -3: X=0.24100000 Y=0.24100000 Z=0.24100000<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -3: X=0.24100000 Y=0.00900000 Z=0.00900000<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -3: X=0.75900000 Y=0.49100000 Z=0.49100000<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -3: X=0.49100000 Y=0.49100000 Z=0.75900000<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -3: X=0.00900000 Y=0.00900000 Z=0.24100000<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -3: X=0.00900000 Y=0.24100000 Z=0.00900000<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -3: X=0.49100000 Y=0.75900000 Z=0.49100000<br>
S&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; NPT=&nbsp;
781&nbsp; R0=0.00010000 RMT=&nbsp;&nbsp;&nbsp; 1.8000&nbsp;&nbsp; Z:
16.0<br>
LOCAL ROT MATRIX:&nbsp;&nbsp;&nbsp; 0.4082483-0.7071068 0.5773503<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
0.4082483 0.7071068 0.5773503<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
-0.8164966 0.0000000 0.5773503<br>
&nbsp; 48&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; NUMBER OF SYMMETRY OPERATIONS<br>
-1 0 0 0.00000000<br>
&nbsp;0-1 0 0.00000000<br>
&nbsp;0 0-1 0.00000000<br>
<br>
FeCr2S4.inorb<br>
&nbsp; 1&nbsp; 2&nbsp;
0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
nmod, natorb, ipr<br>
PRATT&nbsp;
0.3&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
BROYD/PRATT, mixing<br>
&nbsp; 1 1
2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
iatom nlorb, lorb<br>
&nbsp; 2 1
2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
iatom nlorb, lorb<br>
&nbsp;
1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
nsic 0..AFM, 1..SIC, 2..HFM<br>
&nbsp;&nbsp; 0.2205&nbsp; 0.0588&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; U J
(Ry)&nbsp;&nbsp; Note: we recommend to use U_eff = U-J and J=0<br>
&nbsp;&nbsp; 0.1470&nbsp; 0.0588&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; U J<br>
<br>
&nbsp;FeCr2S4.indm<br>
-9.&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
Emin cutoff energy<br>
&nbsp;2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
number of atoms for which density matrix is calculated<br>
&nbsp;1&nbsp; 1&nbsp; 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; index of 1st atom, number of L&#39;s, L1<br>
&nbsp;2&nbsp; 1&nbsp; 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; dtto for 2nd atom, repeat NATOM times<br>
&nbsp;0 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; r-index, (l,s)index<br>
<br>
LAPW magnetic moment which I get which is quite ok but To make it as a insulator I am not able to find out . <br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; MAGNETIC MOMENTS OF MIXED CHARGE DENSITY<br>
:MMINT: MAGNETIC MOMENT IN INTERSTITIAL =&nbsp;&nbsp;&nbsp; 0.25471<br>
:MMI001: MAGNETIC MOMENT IN SPHERE&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -3.29972<br>
:MMI002: MAGNETIC MOMENT IN SPHERE&nbsp;&nbsp; 2&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp; 2.72983<br>
:MMI003: MAGNETIC MOMENT IN SPHERE&nbsp;&nbsp; 3&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp; -0.07183<br>
:MMTOT: TOTAL MAGNETIC MOMENT IN CELL&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp; 3.99994<br>
&nbsp;<br>
I used U and J parameter for this calculation<br>
&nbsp;<br>
&nbsp;Cr&nbsp;&nbsp; U&nbsp;&nbsp;&nbsp;&nbsp; 3&nbsp; 5&nbsp;&nbsp;
7&nbsp;&nbsp; 10&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; J&nbsp; 1&nbsp;
1&nbsp;&nbsp; 1&nbsp;&nbsp; 1<br>
Fe&nbsp;&nbsp;&nbsp; U&nbsp;&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp; 4&nbsp;
6&nbsp;&nbsp;&nbsp;&nbsp; 9&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; J&nbsp;
1&nbsp;&nbsp; 1&nbsp;&nbsp; 1&nbsp;&nbsp; 1<br>
&nbsp;&nbsp; Also I used&nbsp; U_effective <br>
&nbsp;Is it possible to get insulator using LAPW calculation. Is there any silly mistake .<br>
&nbsp;&nbsp; Any kind of help is highly appreciated. If any farther information is requierd you can mail me .&nbsp;  <br>
<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
Swarup.<br>
&nbsp;<br>
&nbsp;<br>
 &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>