<html xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<META HTTP-EQUIV="Content-Type" CONTENT="text/html; charset=us-ascii">
<meta name=Generator content="Microsoft Word 11 (filtered medium)">
<style>
<!--
/* Font Definitions */
@font-face
{font-family:Batang;
panose-1:2 3 6 0 0 1 1 1 1 1;}
@font-face
{font-family:"\@Batang";
panose-1:2 3 6 0 0 1 1 1 1 1;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman";}
a:link, span.MsoHyperlink
{color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{color:purple;
text-decoration:underline;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:Arial;
color:windowtext;}
@page Section1
{size:8.5in 11.0in;
margin:1.0in 1.25in 1.0in 1.25in;}
div.Section1
{page:Section1;}
-->
</style>
</head>
<body lang=EN-US link=blue vlink=purple>
<div class=Section1>
<p class=MsoNormal><font size=2 face=Arial><span style='font-size:10.0pt;
font-family:Arial'>Dear users,<o:p></o:p></span></font></p>
<p class=MsoNormal><font size=2 face=Arial><span style='font-size:10.0pt;
font-family:Arial'><o:p> </o:p></span></font></p>
<p class=MsoNormal><font size=2 face=Arial><span style='font-size:10.0pt;
font-family:Arial'>I’m trying to calculate the optical property of an
semiconductor material, and I’m following the user’s guide for
non-spin polarized calculation with spin-orbit (SO) interaction (page
127). When I didn’t include the SO switch, the resultant optical spectra
( the imaginary part of dielectric matrix) seemed to be reasonable, but when I
did include the SO, the spectra seemed to be very strange. I guess it resulted
from the erroneous estimation of Fermi level in executing lapw2 with
–fermi and –so switch to generate weight files. When SO is
included, the number of band becomes doubled because in this case the spin-up
and spin-down states merge. Then, the number of occupied band should be also
doubled, but I found that the lapw2 execution with –fermi and –so
switch also just gave me the same number of occupied band with the case without
–so switch, and hence giving erroneously low value of Fermi level.
I tried to fix this problem by goring through the source codes, but I couldn’t.
Do you think this is definitely a bug of the program, or I did something wrong?
I attach the first few lines (for one k-point) of ‘case.weight’
files with (case.weight1) and without (case.weight2) –SO flag in
executing lapw2 with –fermi flag. You’ll see the evident error in
determining the Fermi level in the latter case comparing to the former one.
I’m using wien2k version 3 released on 18/01/2003.<o:p></o:p></span></font></p>
<p class=MsoNormal><font size=2 face=Arial><span style='font-size:10.0pt;
font-family:Arial'><o:p> </o:p></span></font></p>
<p class=MsoNormal><font size=2 face=Arial><span style='font-size:10.0pt;
font-family:Arial'>Thanks,<o:p></o:p></span></font></p>
<p class=MsoNormal><font size=2 face=Arial><span style='font-size:10.0pt;
font-family:Arial'><o:p> </o:p></span></font></p>
<p class=MsoNormal><font size=2 face=Arial><span style='font-size:10.0pt;
font-family:Arial'>Chang-Youn Moon<o:p></o:p></span></font></p>
<p class=MsoNormal><font size=2 face=Arial><span style='font-size:10.0pt;
font-family:Arial'>National Renewable Energy Laboratory<o:p></o:p></span></font></p>
<p class=MsoNormal><font size=2 face=Arial><span style='font-size:10.0pt;
font-family:Arial'><o:p> </o:p></span></font></p>
</div>
</body>
</html>