<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" 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 12 (filtered medium)">
<style>
<!--
/* Font Definitions */
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Times New Roman","serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Times New Roman","serif";
color:windowtext;
font-weight:normal;
font-style:normal;}
.MsoChpDefault
{mso-style-type:export-only;}
@page Section1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.Section1
{page:Section1;}
-->
</style>
<!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang=EN-US link=blue vlink=purple>
<div class=Section1>
<p class=MsoNormal>Dear Wien2k users, <o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>I am trying to get DOS for the following structure (already
optimized wrt volume, c/a, b/a and atomic positions) <o:p></o:p></p>
<p class=MsoNormal>UZr2<o:p></o:p></p>
<p class=MsoNormal>P
3 51 Pmma<o:p></o:p></p>
<p class=MsoNormal>
RELA<o:p></o:p></p>
<p class=MsoNormal> 9.411209 5.762730 16.626706 90.000000 90.000000
90.000000<o:p></o:p></p>
<p class=MsoNormal>ATOM -1: X=0.75000000 Y=0.00000000 Z=0.25252548<o:p></o:p></p>
<p class=MsoNormal> MULT=
2 ISPLIT= 8<o:p></o:p></p>
<p class=MsoNormal> -1: X=0.25000000 Y=0.00000000
Z=0.74747452<o:p></o:p></p>
<p class=MsoNormal>Zr1 NPT=
781 R0=0.00001000 RMT= 2.50000 Z: 40.00000<o:p></o:p></p>
<p class=MsoNormal>LOCAL ROT MATRIX: 1.0000000 0.0000000
0.0000000<o:p></o:p></p>
<p class=MsoNormal>
0.0000000 1.0000000 0.0000000<o:p></o:p></p>
<p class=MsoNormal> 0.0000000
0.0000000 1.0000000<o:p></o:p></p>
<p class=MsoNormal>ATOM -2: X=0.25000000 Y=0.50000000 Z=0.07157794<o:p></o:p></p>
<p class=MsoNormal> MULT=
2 ISPLIT= 8<o:p></o:p></p>
<p class=MsoNormal> -2: X=0.75000000 Y=0.50000000
Z=0.92842206<o:p></o:p></p>
<p class=MsoNormal>U 1 NPT=
781 R0=0.00000500 RMT= 2.50000 Z: 92.00000<o:p></o:p></p>
<p class=MsoNormal>LOCAL ROT MATRIX: 1.0000000 0.0000000
0.0000000<o:p></o:p></p>
<p class=MsoNormal>
0.0000000 1.0000000 0.0000000<o:p></o:p></p>
<p class=MsoNormal>
0.0000000 0.0000000 1.0000000<o:p></o:p></p>
<p class=MsoNormal>ATOM -3: X=0.75000002 Y=0.50000000 Z=0.59371330<o:p></o:p></p>
<p class=MsoNormal> MULT=
2 ISPLIT= 8<o:p></o:p></p>
<p class=MsoNormal> -3: X=0.24999998 Y=0.50000000
Z=0.40628670<o:p></o:p></p>
<p class=MsoNormal>Zr2 NPT=
781 R0=0.00000500 RMT= 2.50000 Z: 40.00000<o:p></o:p></p>
<p class=MsoNormal>LOCAL ROT MATRIX: 1.0000000 0.0000000
0.0000000<o:p></o:p></p>
<p class=MsoNormal>
0.0000000 1.0000000 0.0000000<o:p></o:p></p>
<p class=MsoNormal>
0.0000000 0.0000000 1.0000000<o:p></o:p></p>
<p class=MsoNormal> 8 NUMBER OF
SYMMETRY OPERATIONS<o:p></o:p></p>
<p class=MsoNormal> 1 0 0 0.00000000<o:p></o:p></p>
<p class=MsoNormal> 0 1 0 0.00000000<o:p></o:p></p>
<p class=MsoNormal> 0 0 1 0.00000000<o:p></o:p></p>
<p class=MsoNormal> 1<o:p></o:p></p>
<p class=MsoNormal>-1 0 0 0.50000000<o:p></o:p></p>
<p class=MsoNormal> 0-1 0 0.00000000<o:p></o:p></p>
<p class=MsoNormal> 0 0 1 0.00000000<o:p></o:p></p>
<p class=MsoNormal> 2<o:p></o:p></p>
<p class=MsoNormal>-1 0 0 0.00000000<o:p></o:p></p>
<p class=MsoNormal> 0 1 0 0.00000000<o:p></o:p></p>
<p class=MsoNormal> 0 0-1 0.00000000<o:p></o:p></p>
<p class=MsoNormal> 3<o:p></o:p></p>
<p class=MsoNormal> 1 0 0 0.50000000<o:p></o:p></p>
<p class=MsoNormal> 0-1 0 0.00000000<o:p></o:p></p>
<p class=MsoNormal> 0 0-1 0.00000000<o:p></o:p></p>
<p class=MsoNormal> 4<o:p></o:p></p>
<p class=MsoNormal>-1 0 0 0.00000000<o:p></o:p></p>
<p class=MsoNormal> 0-1 0 0.00000000<o:p></o:p></p>
<p class=MsoNormal> 0 0-1 0.00000000<o:p></o:p></p>
<p class=MsoNormal> 5<o:p></o:p></p>
<p class=MsoNormal> 1 0 0 0.50000000<o:p></o:p></p>
<p class=MsoNormal> 0 1 0 0.00000000<o:p></o:p></p>
<p class=MsoNormal> 0 0-1 0.00000000<o:p></o:p></p>
<p class=MsoNormal> 6<o:p></o:p></p>
<p class=MsoNormal> 1 0 0 0.00000000<o:p></o:p></p>
<p class=MsoNormal> 0-1 0 0.00000000<o:p></o:p></p>
<p class=MsoNormal> 0 0 1 0.00000000<o:p></o:p></p>
<p class=MsoNormal> 7<o:p></o:p></p>
<p class=MsoNormal>-1 0 0 0.50000000<o:p></o:p></p>
<p class=MsoNormal> 0 1 0 0.00000000<o:p></o:p></p>
<p class=MsoNormal> 0 0 1 0.00000000<o:p></o:p></p>
<p class=MsoNormal> 8<o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>The system has been run with both non-spin-spin-orbit and
spin-polarized spin-orbit cases. <o:p></o:p></p>
<p class=MsoNormal>It has been run with RK_max=8.00 k-points=20000 with –cc
0.00001 –ec 0.00000001 convergence criteria<o:p></o:p></p>
<p class=MsoNormal>The DOS commands <o:p></o:p></p>
<p class=MsoNormal>Non-Spin <o:p></o:p></p>
<p class=MsoNormal>x lapw1 <o:p></o:p></p>
<p class=MsoNormal>x qtl <o:p></o:p></p>
<p class=MsoNormal>configure_int_lapw <o:p></o:p></p>
<p class=MsoNormal>x tetra <o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>works perfect…I get the DOS <o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>Spin-polarized case <o:p></o:p></p>
<p class=MsoNormal>x lapw1 –up <o:p></o:p></p>
<p class=MsoNormal>x lapw1 –dn <o:p></o:p></p>
<p class=MsoNormal>x qtl –up <o:p></o:p></p>
<p class=MsoNormal>x qtl –dn <o:p></o:p></p>
<p class=MsoNormal>configure_int_lapw <o:p></o:p></p>
<p class=MsoNormal>x tetra –up <o:p></o:p></p>
<p class=MsoNormal>x tetra –dn <o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>works perfect….I get the DOS <o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>However, when I want to get the DOS for the projected
f-states, I have my case.inso file as <o:p></o:p></p>
<p class=MsoNormal>WFFIL<o:p></o:p></p>
<p class=MsoNormal>4 0
0
llmax,ipr,kpot<o:p></o:p></p>
<p class=MsoNormal>-10
5.0
Emin, Emax<o:p></o:p></p>
<p class=MsoNormal> 0 0
1
h,k,l (direction of magnetization)<o:p></o:p></p>
<p class=MsoNormal> 3
number of atoms with RLO<o:p></o:p></p>
<p class=MsoNormal>1 -2.07 0.002
CONT atom-number,
E-param for RLO<o:p></o:p></p>
<p class=MsoNormal>2 -1.27 0.002
CONT
atom-number, E-param for RLO<o:p></o:p></p>
<p class=MsoNormal>3 -2.07 0.002
CONT
atom-number, E-param for RLO<o:p></o:p></p>
<p class=MsoNormal>0 0 number of atoms without
SO, atomnumbers<o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>The case.in1 file is <o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>WFFIL EF=.8507115155 (WFFIL, WFPRI, ENFIL,
SUPWF)<o:p></o:p></p>
<p class=MsoNormal> 8.00
10 4 (R-MT*K-MAX; MAX L IN WF, V-NMT<o:p></o:p></p>
<p class=MsoNormal> .64996 5
0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES,
global APW/LAPW)<o:p></o:p></p>
<p class=MsoNormal> 0 .64996 0.000
CONT 1<o:p></o:p></p>
<p class=MsoNormal> 0 -3.72
0.001 STOP 1<o:p></o:p></p>
<p class=MsoNormal> 1 -2.07 0.002
CONT 1<o:p></o:p></p>
<p class=MsoNormal> 1 .64996 0.000
CONT 1<o:p></o:p></p>
<p class=MsoNormal> 2 .64996 0.005
CONT 1<o:p></o:p></p>
<p class=MsoNormal> .64996 6
0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES,
global APW/LAPW)<o:p></o:p></p>
<p class=MsoNormal> 0 .64996 0.000
CONT 1<o:p></o:p></p>
<p class=MsoNormal> 0 -3.23
0.001 STOP 1<o:p></o:p></p>
<p class=MsoNormal> 1 -1.27
0.002 CONT 1<o:p></o:p></p>
<p class=MsoNormal> 1 .64996 0.000
CONT 1<o:p></o:p></p>
<p class=MsoNormal> 3 .64996 0.005
CONT 1<o:p></o:p></p>
<p class=MsoNormal> 2 .64996 0.005
CONT 1<o:p></o:p></p>
<p class=MsoNormal> .64996 5
0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES,
global APW/LAPW)<o:p></o:p></p>
<p class=MsoNormal> 0 .64996 0.000
CONT 1<o:p></o:p></p>
<p class=MsoNormal> 0 -3.72
0.001 STOP 1<o:p></o:p></p>
<p class=MsoNormal> 1 -2.07
0.002 CONT 1<o:p></o:p></p>
<p class=MsoNormal> 1 .64996 0.000
CONT 1<o:p></o:p></p>
<p class=MsoNormal> 2 .64996 0.005
CONT 1<o:p></o:p></p>
<p class=MsoNormal>K-VECTORS FROM UNIT:4
-9.0 5.0 157 emin /
de (emax=Ef+de) / nband<o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>The case.inq file is <o:p></o:p></p>
<p class=MsoNormal>-9.0
3.0 Emin Emax<o:p></o:p></p>
<p class=MsoNormal>
3
number of atoms<o:p></o:p></p>
<p class=MsoNormal> 1 -2 0
0 iatom,qsplit,symmetrize,locrot<o:p></o:p></p>
<p class=MsoNormal>3 0 1
2 nL, l-values<o:p></o:p></p>
<p class=MsoNormal> 2 -2 0
0 iatom,qsplit,symmetrize,locrot<o:p></o:p></p>
<p class=MsoNormal>4 0 1 2
3 nL, l-values<o:p></o:p></p>
<p class=MsoNormal> 3 -2 0
0 iatom,qsplit,symmetrize,locrot<o:p></o:p></p>
<p class=MsoNormal>3 0 1
2 nL, l-values<o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>Now, if QSPLIT is changed to 0 or -1, and assuming I copy
case.cf_drel from templates to case.cf1 and case.cf3 and case.cf_f_rel to
case.cf2, please suggest on the last line of the case.inq file for “new
axis z” Should it be same as the 3<sup>rd</sup> line of case.inso ? <o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>My second question is<o:p></o:p></p>
<p class=MsoNormal>As a test case, if I wish to get the projected DOS of 6d-U
and 5f-U both, how should the case.cf* file be changed for U. <o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>Thanks <o:p></o:p></p>
<p class=MsoNormal>Suddhasattwa <o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
</div>
</body>
<a></a><p class=""avgcert"" align="left" color="#000000">No virus found in this message.<br>
Checked by AVG - <a href='http://www.avg.com'>www.avg.com</a><br>
Version: 2014.0.4744 / Virus Database: 4007/8025 - Release Date: 08/12/14<br>
Internal Virus Database is out of date.</p></html>