
<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">

<HTML><HEAD>

<META http-equiv=Content-Type content="text/html; charset=us-ascii">

<META content="MSHTML 6.00.6000.16608" name=GENERATOR></HEAD>

<BODY>

<DIV dir=ltr align=left><FONT face=Arial size=2><SPAN 

class=748280117-12022008>Hi,</SPAN></FONT></DIV>

<DIV><FONT face=Arial size=2><SPAN class=748280117-12022008>I am using 

Wien2k_08.1. I'd like to calculate the optical properties of GaAs using the 

OPTIC&nbsp;function in Wien. I have run the calculations without including the 

interband transitions and it is fine but when I include the interband 

contributions, the programme crashes at<SPAN class=443402020-18022008><FONT 

color=#0000ff><FONT color=#000000>&nbsp;the</FONT> </FONT></SPAN>Kram stage. I 

think the reason for&nbsp;it is that when the interband contributions are 

included the eps1 and eps2 outputs of joint are no longer valid and hence Kram 

crashes. This is the procedure that&nbsp;I have followed:</SPAN></FONT></DIV>

<DIV><FONT face=Arial size=2><SPAN 

class=748280117-12022008></SPAN></FONT>&nbsp;</DIV>

<DIV><FONT face=Arial><SPAN class=748280117-12022008>&nbsp;<FONT 

face="Courier New" size=2>1) run_lapw <BR>&nbsp;&nbsp; <BR>&nbsp; 2) x kgen 

(with 20000 k points)<BR>&nbsp;&nbsp; <BR>&nbsp; 3) edit case.in2c, changing TOT 

to FERMI and writing 101.0 instead of 0.000 in TETRA<BR>&nbsp;&nbsp; <BR>&nbsp; 

4) x lapw1<BR>&nbsp;&nbsp; <BR>&nbsp; 5) x lapw2<BR>&nbsp;&nbsp; <BR>&nbsp; 6) 

Modify case.inop:<BR>&nbsp;&nbsp;&nbsp;<BR>&nbsp;&nbsp;99999 

1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; number of k-points, first 

k-point<BR>&nbsp;&nbsp;&nbsp;<BR>&nbsp; 0.0 2.2&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; 

Emin, Emax for matrix elements<BR>&nbsp;&nbsp; <BR>&nbsp; 

2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; number 

of choices (columns in *outmat) - 0: MME into case.mme<BR>&nbsp;&nbsp; 

<BR>&nbsp; 

1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Re 

xx<BR>&nbsp;&nbsp; <BR>&nbsp; 

3&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Re 

zz<BR>&nbsp;&nbsp; <BR>&nbsp; 

OFF&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

ON/OFF&nbsp;&nbsp; writes MME to unit 4<BR>&nbsp;&nbsp; <BR>&nbsp; 7)&nbsp; x 

opticc<BR>&nbsp;&nbsp; <BR>&nbsp; 8)&nbsp; Modify case.injoint<BR>&nbsp;&nbsp; 

<BR>&nbsp; 1&nbsp;&nbsp;9999 

9999&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; : LOWER,UPPER 

and (optional) UPPER-VAL BANDINDEX<BR>&nbsp;&nbsp; <BR>&nbsp;&nbsp; 

0.0000&nbsp;&nbsp;&nbsp; 0.00100&nbsp;&nbsp; 2.0000 : EMIN DE EMAX FOR 

ENERGYGRID IN ryd<BR>&nbsp;&nbsp; <BR>&nbsp; 

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;<BR>&nbsp;&nbsp;6&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; <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;&nbsp; 

: NUMBER OF COLUMNS<BR>&nbsp;&nbsp; <BR>&nbsp; 0.1&nbsp; 0.1&nbsp; 

0.3&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; : 

BROADENING (FOR DRUDE MODEL - switch 6,7 -ONLY<BR>&nbsp;&nbsp; <BR>&nbsp; 9) x 

joint<BR>&nbsp;&nbsp; <BR>&nbsp; 10) Modify case.inkram<BR>&nbsp;&nbsp; 

<BR>&nbsp; 11) It crashes immediately giving the error 

message:</FONT></SPAN></FONT></DIV>

<DIV><FONT face="Courier New" size=2><SPAN 

class=748280117-12022008></SPAN></FONT>&nbsp;</DIV>

<DIV><FONT face=Arial><SPAN class=748280117-12022008><FONT size=2><FONT 

face="Courier New">&nbsp;&nbsp;&nbsp; forrtl: severe (64): input conversion 

error, unit 10,&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; file 

/home/sk701/Wien2k/GaAs_12-02-08/GaAs_12-02-08.joint<BR>Image&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

PC&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

Routine&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

Line&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

Source&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

<BR>kram&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

000000000044BE97&nbsp; 

Unknown&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

Unknown&nbsp; 

Unknown<BR>kram&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

000000000044A25E&nbsp; 

Unknown&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

Unknown&nbsp; 

Unknown<BR>kram&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

000000000042E274&nbsp; 

Unknown&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

Unknown&nbsp; 

Unknown<BR>kram&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

00000000004094E3&nbsp; 

Unknown&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

Unknown&nbsp; 

Unknown<BR>kram&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

000000000040914C&nbsp; 

Unknown&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

Unknown&nbsp; 

Unknown<BR>kram&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

0000000000418218&nbsp; 

Unknown&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

Unknown&nbsp; 

Unknown<BR>kram&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

0000000000402CD0&nbsp; 

Unknown&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

Unknown&nbsp; 

Unknown<BR>kram&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

000000000040272A&nbsp; 

Unknown&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

Unknown&nbsp; 

Unknown<BR>libc.so.6&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

00002AAAAAE4B5AA&nbsp; 

Unknown&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

Unknown&nbsp; 

Unknown<BR>kram&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

000000000040266A&nbsp; 

Unknown&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

Unknown&nbsp; Unknown<BR>0.000u 0.000s 0:00.00 0.0%&nbsp;0+0k 0+0io 

0pf+0w<BR>error: command&nbsp;&nbsp; /usr/local/wien2k_new/kram 

kram.def&nbsp;&nbsp; failed</FONT></FONT></SPAN></FONT></DIV>

<DIV><FONT face="Courier New" size=2><SPAN 

class=748280117-12022008></SPAN></FONT>&nbsp;</DIV>

<DIV><FONT face=Arial size=2><SPAN class=748280117-12022008>I have&nbsp;searched 

through the mailing list and other people seemed to have had the same problem 

but I couldn't find an answer to anyone's problem posted on the list. I have run 

the calculations for TiC as well but the same problem occurs. I would appreciate 

your suggestions on how to solve this problem.</SPAN></FONT></DIV>

<DIV><FONT face=Arial size=2><SPAN 

class=748280117-12022008>Regards,</SPAN></FONT></DIV>

<DIV><FONT face=Arial size=2><SPAN 

class=748280117-12022008>Shima</SPAN></FONT></DIV>

<DIV><FONT face=Arial><SPAN class=748280117-12022008><FONT size=2><FONT 

face="Courier New"></FONT></FONT></SPAN></FONT>&nbsp;</DIV>

<DIV><FONT face=Arial><SPAN class=748280117-12022008><FONT size=2><FONT 

face="Courier New"><FONT face=Arial></FONT>&nbsp;</DIV>

<DIV><BR></DIV></FONT></FONT></SPAN></FONT></BODY></HTML>