&nbsp;&nbsp; <br>
<br>
&nbsp;&nbsp; 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;&nbsp;&nbsp;&nbsp;
I am&nbsp; doing optical calculation for&nbsp; fcc(001)&nbsp; Co&nbsp;
with&nbsp; [001]&nbsp;&nbsp; magnetisation. Changes in the files <br>
&nbsp; case.inso <br>
---------------------------------------------------------------------------------------<br>
&nbsp;&nbsp;&nbsp; WFFIL<br>
&nbsp;4&nbsp; 1&nbsp;
0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
llmax,ipr,kpot<br>
&nbsp;-10.0000&nbsp;&nbsp; 1.50000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; emin,emax (output energy window)<br>
&nbsp;&nbsp; 0.&nbsp; 0.&nbsp;
1.&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
direction of magnetization (lattice vectors)<br>
&nbsp;0&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;
number of atoms for which RLO is added<br>
&nbsp;0&nbsp; -4.97&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
0.005&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; atom number,e-lo,de
(case.in1), repeat NX times<br>
&nbsp;0 0 0 0
0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
number of atoms for which SO is switch off; atoms<br>
-----------------------------------------------------------------------------------------<br>
case.in2c<br>
----------------------------------------------------------------------------------------<br>
<br>
FERMI&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; (TOT,FOR,QTL,EFG,FERMI)<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -9.0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 15.0
0.50
0.05&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
EMIN, NE, ESEPERMIN, ES<br>
TETRA&nbsp;&nbsp;&nbsp;
101&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
(GAUSS,ROOT,TEMP,TETRA,ALL&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; eval)<br>
&nbsp; 0 0&nbsp; 2 0&nbsp; 4 0&nbsp; 4 4&nbsp; 6 0&nbsp; 6 4<br>
&nbsp;14.&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; GMAX<br>
FILE&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; FILE/NOFILE&nbsp; write recprlist<br>
~<br>
--------------------------------------------------------------------------------------------<br>
Changed&nbsp; structure files&nbsp; after running initso&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; case.struct<br>
<br>
-----------------------------------------------------------------------------------------------<br>
Co&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;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
s-o calc. M||&nbsp; 0.00&nbsp; 0.00&nbsp; 1.00<br>
F&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;
1 25_F<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; RELA<br>
&nbsp; 6.689633&nbsp; 6.689633&nbsp; 6.689633 90.000000 90.000000 90.000000<br>
ATOM&nbsp;&nbsp; 1: X=0.00000000 Y=0.00000000 Z=0.00000000<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; MULT= 1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ISPLIT= 2<br>
Co&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; NPT=&nbsp; 781&nbsp;
R0=.000100000 RMT=&nbsp;&nbsp; 2.35000&nbsp;&nbsp; Z:&nbsp; 27.00000<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>
&nbsp; 16&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; NUMBER OF SYMMETRY OPERATIONS<br>
-1 0 0 0.0000000<br>
&nbsp;0-1 0 0.0000000<br>
&nbsp;0 0-1 0.0000000<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1<br>
&nbsp;0 1 0 0.0000000<br>
-1 0 0 0.0000000<br>
&nbsp;0 0-1 0.0000000<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2<br>
&nbsp;0-1 0 0.0000000<br>
&nbsp;1 0 0 0.0000000<br>
&nbsp;0 0-1 0.0000000<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 3<br>
-----------------------------------------------------------------------------------------------<br>
<br>
&nbsp;After&nbsp; init_lapw&nbsp; i run<br>
&nbsp;&nbsp; runsp_lapw<br>
&nbsp;&nbsp;&nbsp; initso_lapw <br>
&nbsp;&nbsp;&nbsp; then i again run init_lapw&nbsp; due to change of
the structure files&nbsp; after this calculation i&nbsp; ran<br>
&nbsp;&nbsp;&nbsp; runsp_lapw -so <br>
<br>
&nbsp; for optic calculation <br>
&nbsp;&nbsp; <br>
<br>
&nbsp;&nbsp; change TOT to FERMI in case.in2c<br>
&nbsp;&nbsp; runsp_lapw -so -s lapw1 -e lcore<br>
&nbsp;&nbsp; x opticc -so -up<br>
&nbsp;&nbsp; x joint -up<br>
&nbsp;&nbsp; x kram -up<br>
<br>
here the three files<br>
&nbsp;case.inop<br>
-------------------------------------------------------------------------------------------<br>
244 1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; number of k-points, first k-point<br>
-5.0 2.2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Emin, Emax for matrix elements<br>
1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
number of choices (columns in *outmat) - 0: MME into case.mme<br>
1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Re xx<br>
OFF&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ON/OFF&nbsp;&nbsp; writes MME to unit 4<br>
<br>
Choices:<br>
1......Re &lt;x&gt;&lt;x&gt;<br>
2......Re &lt;y&gt;&lt;y&gt;<br>
3......Re &lt;z&gt;&lt;z&gt;<br>
4......Re &lt;x&gt;&lt;y&gt;<br>
5......Re &lt;x&gt;&lt;z&gt;<br>
6......Re &lt;y&gt;&lt;z&gt;<br>
7......Im &lt;x&gt;&lt;y&gt;<br>
8......Im &lt;x&gt;&lt;z&gt;<br>
9......Im &lt;y&gt;&lt;z&gt;<br>
----------------------------------------------------------------------------------------------------<br>
<br>
<br>
case.injoint&nbsp; 1&nbsp;&nbsp;
20&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
: LOWER AND UPPER BANDINDEX<br>
&nbsp;&nbsp; 0.0000&nbsp;&nbsp;&nbsp; 0.00100&nbsp;&nbsp; 1.0000 : EMIN DE EMAX FOR ENERGYGRID IN ryd<br>
eV&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;
: output units&nbsp; eV / ryd&nbsp; / cm-1<br>
&nbsp;&nbsp;&nbsp;&nbsp;
4&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
: SWITCH<br>
&nbsp;&nbsp;&nbsp;&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;
: NUMBER OF COLUMNS<br>
&nbsp;&nbsp; 0.1&nbsp; 0.1&nbsp;
0.3&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
: BROADENING (FOR DRUDE MODEL - switch 6,7 -<br>
ONLY)<br>
<br>
SWITCH:<br>
<br>
&nbsp;&nbsp; 0...JOINTDOS FOR EACH BAND COMBINATION<br>
&nbsp;&nbsp; 1...JOINTDOS AS SUM OVER ALL BAND COMBINATIONS<br>
&nbsp;&nbsp; 2...DOS FOR EACH BAND<br>
&nbsp;&nbsp; 3...DOS AS SUM OVER ALL BANDS<br>
&nbsp;&nbsp; 4...Im(EPSILON)<br>
&nbsp;&nbsp; 5...Im(EPSILON) for each band combination<br>
&nbsp;&nbsp; 6...INTRABAND contributions<br>
&nbsp;&nbsp; 7...INTRABAND contributions including band analysis<br>
<br>
-----------------------------------------------------------------------------------------------<br>
<br>
&nbsp;&nbsp; case.inkram<br>
------------------------------------------------------------------------------------------------<br>
----------------------------------------------------------------------------------------------------<br>
<br>
&nbsp;0.1&nbsp;&nbsp; Gamma: broadening of interband spectrum<br>
&nbsp; 1.0&nbsp;&nbsp;&nbsp; energy shift (scissors operator)<br>
&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; add intraband contributions? yes/no: 1/0<br>
&nbsp;12.60&nbsp;&nbsp; plasma frequencies&nbsp; (from joint, opt 6)<br>
&nbsp; 0.20&nbsp;&nbsp; Gammas for Drude terms<br>
~<br>
-----------------------------------------------------------------------------------------------------<br>
<br>
<br>
I think here are all the necessary changes are done as&nbsp;
manual.&nbsp; Then i plot case.sigmak&nbsp; for optical plot&nbsp; and
plot is not matches with published fig.&nbsp; I checked&nbsp; the
calculation very well but i can not what is the wrong here.&nbsp; <br>
&nbsp; <br>
<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; I have another question&nbsp; why we write atom no 0&nbsp; in the case.inso&nbsp; files .&nbsp; <br>
<br>
<br>
&nbsp;&nbsp;&nbsp;&nbsp; Any kind of suggestion regarding this problem highly appreciated.<br>
<br>
<br>
&nbsp;&nbsp;&nbsp; With best&nbsp; regards,<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>
<br>
<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
Swarup&nbsp;&nbsp; Saha<br>
<br>
<br>
<br>
&nbsp;&nbsp; <br>
<br>
&nbsp;&nbsp;&nbsp; <br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>