<div dir="ltr">Dear Prof blaha and wien users<br><br> I have a problem in optic with spin orbit coupling. I am using the latest version of wien2k (wien2k10.1). <br><br>I run the following script to calculate the optical properties <br>
<br><br>runsp -so -orb -s lapw1 -e lcore (with large k points)<br>x lapw2 -fermi -so -up<br>x opticc -so -up<br><br>case.inso file ...........................<br><br>WFFIL<br>4 0 0 llmax,ipr,kpot<br>-10 1.5 Emin, Emax<br>
0 0 1 h,k,l (direction of magnetization)<br> 1 number of atoms with RLO<br>1 -3.96 0.001 STOP atom-number, E-param for RLO<br>0 0 number of atoms without SO, atomnumbers<br>
<br><br>case.mat_diagup file looks like ..............<br><br><br> 3 Re <x><x> Re <z><z> Im <x><y> ALL<br><br> KP: 1 NEMIN NEMAX : 1 18 dE:-5.00 2.20 K: 1<br>
<br> 1 1 NaN NaN NaN<br> 2 2 NaN NaN NaN<br> 3 3 NaN NaN NaN<br> 4 4 NaN NaN NaN<br> 5 5 NaN NaN NaN<br>
6 6 NaN NaN NaN<br> 7 7 NaN NaN NaN<br> 8 8 NaN NaN NaN<br> 9 9 NaN NaN NaN<br> 10 10 NaN NaN NaN<br>
11 11 NaN NaN NaN<br> 12 12 NaN NaN NaN<br> 13 13 NaN NaN NaN<br> 14 14 NaN NaN NaN<br> 15 15 NaN NaN NaN<br>
16 16 NaN NaN NaN<br> 17 17 NaN NaN NaN<br> 18 18 NaN NaN NaN<br><br> KP: 2 NEMIN NEMAX : 1 18 dE:-5.00 2.20 K: 2<br><br> 1 1 NaN NaN NaN<br>
2 2 NaN NaN NaN<br> 3 3 NaN NaN NaN<br> 4 4 NaN NaN NaN<br> 5 5 NaN NaN NaN<br> 6 6 NaN NaN NaN<br>
7 7 NaN NaN NaN<br> 8 8 NaN NaN NaN<br> 9 9 NaN NaN NaN<br> 10 10 NaN NaN NaN<br> 11 11 NaN NaN NaN<br>
12 12 NaN NaN NaN<br> 13 13 NaN NaN NaN<br> 14 14 NaN NaN NaN<br> 15 15 NaN NaN NaN<br> 16 16 NaN NaN NaN<br>
17 17 NaN NaN NaN<br> 18 18 NaN NaN NaN<br><br> KP: 3 NEMIN NEMAX : 1 18 dE:-5.00 2.20 K: 3<br><br> 1 1 NaN NaN NaN<br> 2 2 NaN NaN NaN<br>
3 3 NaN NaN NaN<br> 4 4 NaN NaN NaN<br> 5 5 NaN NaN NaN<br> 6 6 NaN NaN NaN<br> 7 7 NaN NaN NaN<br>
"feso.mat_diagup" 2273L, 107420C <br>and continue..................<br><br><br>looking forward to hear from you <br>Thanking you<br><br>with regards<br>nsingh<br></div>