<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD>
<META http-equiv=Content-Type content="text/html; charset=iso-8859-1">
<META content="MSHTML 6.00.6000.16674" name=GENERATOR>
<STYLE></STYLE>
</HEAD>
<BODY bgColor=#ffffff>
<DIV><FONT face=Arial size=2>Dear All:</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>After I sent out the a question about <FONT
size=1>Optical calculation(see <A
href="">http://zeus.theochem.tuwien.ac.at/pipermail/wien/2008-July/011200.html</A> )</FONT>,
I </FONT><FONT face=Arial size=2>have read the following link: </FONT><FONT
face=Arial size=2><A
href="">http://zeus.theochem.tuwien.ac.at/pipermail/wien/2008-February/010503.html</A></FONT></DIV>
<DIV><FONT face=Arial size=2>I follow the procedure as follows:</FONT></DIV>
<DIV><FONT face=Arial><EM> 1) run_lapw<BR></EM><I> 2) x kgen (with
100000 k points)</I></FONT></DIV>
<DIV><FONT face=Arial><EM> 3) edit case.in2c, changing TOT to FERMI and
writing 101.0 instead of <BR></EM><EM>0.000 in TETRA<BR></EM><EM> 4)
run_lapw -p</EM></FONT></DIV>
<DIV><FONT face=Arial><I> 5) x lapw2 -qtl -p<BR></I><I>
6) edit case.inop: </I></FONT></DIV>
<DIV><FONT
face=Arial><I>_______________________________________</I></FONT></DIV>
<DIV><FONT face=Arial><I>2456 1 number of
k-points, first k-point <BR>-9.0 4.0 Emin, Emax
for matrix
elements<BR>2
number of choices (columns in *outmat) - 0: MME into
case.mme<BR>1
Re
xx<BR>3
Re zz<BR>ON
ON/OFF writes MME to unit 4</I></FONT></DIV>
<DIV><EM><FONT
face=Arial>_______________________________________</FONT></EM></DIV>
<DIV><EM><FONT face=Arial> 7) X OPTIC -P</FONT></EM></DIV>
<DIV><EM><FONT face=Arial> </FONT></EM></DIV>
<DIV><EM><FONT face=Arial>8) edit case.injoint</FONT></EM></DIV>
<DIV><EM><FONT
face=Arial>_______________________________________</FONT></EM></DIV>
<DIV><FONT face=Arial size=2> 1 2456
2456
: LOWER,UPPER and (optional) UPPER-VAL BANDINDEX<BR>
-9.0000 0.0100 4.0000 : EMIN DE EMAX FOR
ENERGYGRID IN ryd
<BR>eV
: output units eV / ryd / cm-1<BR>
4
: SWITCH <BR>
2
: NUMBER OF COLUMNS<BR> 0.1 0.1
0.1
: BROADENING (FOR DRUDE MODEL - switch 6,7 - <BR>ONLY)<BR><EM><FONT
size=3>_______________________________________</FONT></EM></FONT></DIV>
<DIV><EM><FONT face=Arial></FONT></EM> </DIV>
<DIV><EM><FONT face=Arial>9) x joint</FONT></EM></DIV>
<DIV><EM><FONT face=Arial>10)edit case.inkram</FONT></EM></DIV>
<DIV><EM><FONT
face=Arial>_______________________________________</FONT></EM></DIV>
<DIV><EM><FONT face=Arial> 0.05 Gamma: broadening of
interband spectrum<BR> 0.4 energy shift (scissors
operator)<BR> 1 add intraband contributions?
yes/no: 1/0<BR> 11.02 plasma frequencies (from joint, opt
6)<BR> 0.20 Gammas for Drude terms </FONT></EM></DIV>
<DIV><EM><FONT
face=Arial>_______________________________________</FONT></EM></DIV>
<DIV><EM><FONT face=Arial>11) x kram</FONT></EM></DIV>
<DIV><EM><FONT face=Arial>then I got following error message:</FONT></EM></DIV>
<DIV><EM><FONT
face=Arial>_______________________________________</FONT></EM></DIV>
<DIV>xx<BR>zz<BR> Energy units: [eV]<BR> Lorentzian broadening
with gamma: 5.000000000000000E-002 [eV]<BR>forrtl: severe
(64): input conversion error, unit 10, file
/home/dongsu/calculation/MgO1/MgO1.joint<BR>Image
PC
Routine
Line
Source
<BR>kram
000000000044AB7A
Unknown
Unknown
Unknown<BR>kram
000000000044905E
Unknown
Unknown
Unknown<BR>kram
000000000042DCAC
Unknown
Unknown
Unknown<BR>kram
00000000004090BA
Unknown
Unknown
Unknown<BR>kram
0000000000408D1B
Unknown
Unknown
Unknown<BR>kram
00000000004179D7
Unknown
Unknown
Unknown<BR>kram
0000000000403364
Unknown
Unknown
Unknown<BR>kram
000000000040266A
Unknown
Unknown
Unknown<BR>libc.so.6
0000003E8401C3FB
Unknown
Unknown
Unknown<BR>kram
00000000004025AA
Unknown
Unknown Unknown<BR>0.001u 0.002s 0:01.94 0.0% 0+0k 0+0io 1pf+0w<BR>error:
command /home/dongsu/wien2k.08/kram kram.def
failed<BR><BR><FONT face=Arial
size=2>_______________________________________________________</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>Could any one tell me what is the mistake I
had?</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>And I still looking forward to hear
the answer my last e-mail: o<FONT size=1>ptical calculation on the case
without spin-orbit(<A
href="">http://zeus.theochem.tuwien.ac.at/pipermail/wien/2008-July/011200.html</A>)</FONT></FONT></DIV>
<DIV><FONT face=Arial size=1></FONT> </DIV>
<DIV><FONT face=Arial size=1>Thanks a lot!!</FONT></DIV>
<DIV><FONT face=Arial size=1></FONT> </DIV>
<DIV><FONT face=Arial size=1>Dong</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>Dong Su<BR>Center for Functional
Nanomaterials<BR>Brookhaven National Laboratory<BR>Building 735<BR>Upton, New
York 11973<BR><A href="">dsu@bnl.gov</A> </FONT></DIV></BODY></HTML>