<div>Because p -1/2 ''relativistic' LOs is not supported on OPTICS program, atom no 0&nbsp; in the case.inso&nbsp; file is used.</div>  <div>&nbsp;</div>  <div>Haiquan Hu</div>  <div>Liaocheng University</div>  <div>Shandong Provence</div>  <div>China</div>  <div>&nbsp;&nbsp;</div>  <div><BR><BR><B><I>swarup saha &lt;saha18swarup@gmail.com&gt;</I></B> wrote:</div>  <BLOCKQUOTE class=replbq style="PADDING-LEFT: 5px; MARGIN-LEFT: 5px; BORDER-LEFT: #1010ff 2px solid">&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>_______________________________________________<BR>Wien mailing list<BR>Wien@zeus.theochem.tuwien.ac.at<BR>http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien<BR></BLOCKQUOTE><BR><p>&#32;
        



<hr size=1>Everyone is raving about <a href="http://us.rd.yahoo.com/evt=42297/*http://advision.webevents.yahoo.com/mailbeta">the all-new Yahoo! Mail beta.</a>