<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:Wingdings;
        panose-1:5 0 0 0 0 0 0 0 0 0;}
@font-face
        {font-family:"Nimbus Mono L";}
 /* Style Definitions */
 p.MsoNormal, li.MsoNormal, div.MsoNormal
        {margin:0in;
        margin-bottom:.0001pt;
        text-align:right;
        direction:rtl;
        unicode-bidi:embed;
        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;}
p.preformattedtext, li.preformattedtext, div.preformattedtext
        {margin:0in;
        margin-bottom:.0001pt;
        font-size:10.0pt;
        font-family:"Nimbus Mono L";}
span.EmailStyle18
        {mso-style-type:personal;
        font-family:Arial;
        color:windowtext;}
span.EmailStyle19
        {mso-style-type:personal-reply;
        font-family:Arial;
        color:navy;}
@page Section1
        {size:595.3pt 841.9pt;
        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 dir=RTL>

<p class=MsoNormal dir=LTR style='text-align:left;direction:ltr;unicode-bidi:
embed'><font size=3 face="Times New Roman"><span style='font-size:12.0pt'>Dear
wien2k users,<o:p></o:p></span></font></p>

<p class=MsoNormal dir=LTR style='text-align:left;direction:ltr;unicode-bidi:
embed'><font size=3 face="Times New Roman"><span style='font-size:12.0pt'><o:p>&nbsp;</o:p></span></font></p>

<p class=MsoNormal dir=LTR style='text-align:left;direction:ltr;unicode-bidi:
embed'><font size=3 face="Times New Roman"><span style='font-size:12.0pt'>I
wish to do calculations for UN compound with LDA+U potential. After studying
the UG and the relevant issues in the mailing list I did a calculations with SO
which converged without problems. Afterwards I tried to apply OP option to the
f- and d-orbitals of Uranium, but unsuccessfully (the result of this
&quot;experiment&quot; was error). Form the mailing list I understood that it
is impossible to apply OP option for both states. However I don't understand if
it is a technical or a physical reason for that. <o:p></o:p></span></font></p>

<p class=MsoNormal dir=LTR style='text-align:left;direction:ltr;unicode-bidi:
embed'><font size=3 face="Times New Roman"><span style='font-size:12.0pt'><o:p>&nbsp;</o:p></span></font></p>

<p class=preformattedtext dir=LTR style='margin-bottom:14.15pt'><font size=3
face="Times New Roman"><span style='font-size:12.0pt;font-family:"Times New Roman"'>When
I apply OP option only for d-orbitals, the SCF calculations work fine. But when
I change d-orbitals to f-orbitals in both case.inorb and case.indm files the
SCF calculations stop with the following error:<o:p></o:p></span></font></p>

<p class=preformattedtext dir=LTR><font size=3 face="Nimbus Mono L"><span
style='font-size:12.0pt'>start&nbsp;&nbsp; (Thu Sep&nbsp; 8 13:31:01 IDT 2005)
with lapw0 (30/40 to go)<o:p></o:p></span></font></p>

<p class=preformattedtext dir=LTR><font size=2 face="Nimbus Mono L"><span
style='font-size:10.0pt'><o:p>&nbsp;</o:p></span></font></p>

<p class=preformattedtext dir=LTR><font size=2 face="Nimbus Mono L"><span
style='font-size:10.0pt'>&nbsp;&nbsp;&nbsp; cycle 1&nbsp;&nbsp;&nbsp;&nbsp;
(Thu Sep&nbsp; 8 13:31:01 IDT 2005)&nbsp; (30/40 to go)<o:p></o:p></span></font></p>

<p class=preformattedtext dir=LTR><font size=2 face="Nimbus Mono L"><span
style='font-size:10.0pt'><o:p>&nbsp;</o:p></span></font></p>

<p class=preformattedtext dir=LTR><font size=2 face="Nimbus Mono L"><span
style='font-size:10.0pt'>&gt;&nbsp;&nbsp; lapw0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
(13:31:01) 3.152u 0.039s 0:03.25 97.8%&nbsp; 0+0k 0+0io 0pf+0w<o:p></o:p></span></font></p>

<p class=preformattedtext dir=LTR><font size=2 face="Nimbus Mono L"><span
style='font-size:10.0pt'>&gt;&nbsp;&nbsp; orb -up&nbsp;&nbsp;&nbsp;&nbsp;
(13:31:04) 0.005u 0.008s 0:00.01 0.0%&nbsp;&nbsp; 0+0k 0+0io 0pf+0w<o:p></o:p></span></font></p>

<p class=preformattedtext dir=LTR><font size=2 face="Nimbus Mono L"><span
style='font-size:10.0pt'><o:p>&nbsp;</o:p></span></font></p>

<p class=preformattedtext dir=LTR><font size=2 face="Nimbus Mono L"><span
style='font-size:10.0pt'>&gt;&nbsp;&nbsp; stop error<o:p></o:p></span></font></p>

<p class=preformattedtext dir=LTR><font size=2 face="Nimbus Mono L"><span
style='font-size:10.0pt'><o:p>&nbsp;</o:p></span></font></p>

<p class=preformattedtext dir=LTR><font size=2 face="Nimbus Mono L"><span
style='font-size:10.0pt'><o:p>&nbsp;</o:p></span></font></p>

<p class=MsoNormal dir=LTR style='text-align:left;direction:ltr;unicode-bidi:
embed'><font size=3 face="Times New Roman"><span style='font-size:12.0pt'>The
error is in the file uporb.error &quot;Error in Vorb&quot;.<o:p></o:p></span></font></p>

<p class=preformattedtext dir=LTR><font size=2 face="Nimbus Mono L"><span
style='font-size:10.0pt'><o:p>&nbsp;</o:p></span></font></p>

<p class=preformattedtext dir=LTR><font size=2 face="Nimbus Mono L"><span
style='font-size:10.0pt'><o:p>&nbsp;</o:p></span></font></p>

<p class=preformattedtext dir=LTR><font size=3 face="Times New Roman"><span
style='font-size:12.0pt;font-family:"Times New Roman"'>Here are case.inorb and
case.indm files in my case:<o:p></o:p></span></font></p>

<p class=preformattedtext dir=LTR><font size=2 face="Nimbus Mono L"><span
style='font-size:10.0pt'><o:p>&nbsp;</o:p></span></font></p>

<p class=preformattedtext dir=LTR style='margin-bottom:14.15pt'><font size=2
face="Nimbus Mono L"><span style='font-size:10.0pt'>UN_1.inorb<o:p></o:p></span></font></p>

<p class=MsoNormal dir=LTR style='text-align:left;direction:ltr;unicode-bidi:
embed'><font size=3 face="Times New Roman"><span style='font-size:12.0pt'>1&nbsp;
1&nbsp;
0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
nmod, natorb, ipr<o:p></o:p></span></font></p>

<p class=MsoNormal dir=LTR style='text-align:left;direction:ltr;unicode-bidi:
embed'><font size=3 face="Times New Roman"><span style='font-size:12.0pt'>PRATT&nbsp;
1.0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
BROYD/PRATT, mixing<o:p></o:p></span></font></p>

<p class=MsoNormal dir=LTR style='text-align:left;direction:ltr;unicode-bidi:
embed'><font size=3 face="Times New Roman"><span style='font-size:12.0pt'>&nbsp;
1 1 3&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;iatom
nlorb,
lorb&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;
</span></font><font face=Wingdings><span style='font-family:Wingdings'>&szlig;</span></font>for
2 instead of 3 the calculations work fine<o:p></o:p></p>

<p class=MsoNormal dir=LTR style='text-align:left;direction:ltr;unicode-bidi:
embed'><font size=3 face="Times New Roman"><span style='font-size:12.0pt'>&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;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
nsic 0..AFM, 1..SIC, 2..HFM &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;&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;<o:p></o:p></span></font></p>

<p class=MsoNormal dir=LTR style='text-align:left;direction:ltr;unicode-bidi:
embed'><font size=3 face="Times New Roman"><span style='font-size:12.0pt'>&nbsp;&nbsp;
0.8 &nbsp;0.00&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; U J (Ry)&nbsp;&nbsp;
Note: we recommend to use U_eff = U-J and J=0<o:p></o:p></span></font></p>

<p class=MsoNormal dir=LTR style='text-align:left;direction:ltr;unicode-bidi:
embed'><font size=3 face="Times New Roman"><span style='font-size:12.0pt'><o:p>&nbsp;</o:p></span></font></p>

<p class=MsoNormal dir=LTR style='text-align:left;direction:ltr;unicode-bidi:
embed'><font size=3 face="Times New Roman"><span style='font-size:12.0pt'><o:p>&nbsp;</o:p></span></font></p>

<p class=MsoNormal dir=LTR style='text-align:left;direction:ltr;unicode-bidi:
embed'><font size=3 face="Times New Roman"><span style='font-size:12.0pt'>&nbsp;UN_1.indm&nbsp;
<o:p></o:p></span></font></p>

<p class=MsoNormal dir=LTR style='text-align:left;direction:ltr;unicode-bidi:
embed'><font size=3 face="Times New Roman"><span style='font-size:12.0pt'>-9.&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
Emin cutoff energy<o:p></o:p></span></font></p>

<p class=MsoNormal dir=LTR style='text-align:left;direction:ltr;unicode-bidi:
embed'><font size=3 face="Times New Roman"><span style='font-size:12.0pt'>&nbsp;1&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 density matrix is calculated<o:p></o:p></span></font></p>

<p class=MsoNormal dir=LTR style='text-align:left;direction:ltr;unicode-bidi:
embed'><font size=3 face="Times New Roman"><span style='font-size:12.0pt'>&nbsp;1&nbsp;
1&nbsp; 3&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; index of 1st atom, number of L's,
L1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
&nbsp;&nbsp;&nbsp;&nbsp;</span></font><font face=Wingdings><span
style='font-family:Wingdings'>&szlig;</span></font>for 2 instead of 3 the
calculations work fine<o:p></o:p></p>

<p class=MsoNormal dir=LTR style='text-align:left;direction:ltr;unicode-bidi:
embed'><font size=3 face="Times New Roman"><span style='font-size:12.0pt'>&nbsp;0
0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; r-index,
(l,s)index&nbsp; <o:p></o:p></span></font></p>

<p class=MsoNormal dir=LTR style='text-align:left;direction:ltr;unicode-bidi:
embed'><font size=3 face="Times New Roman"><span style='font-size:12.0pt'><o:p>&nbsp;</o:p></span></font></p>

<p class=MsoNormal dir=LTR style='text-align:left;direction:ltr;unicode-bidi:
embed'><font size=3 face="Times New Roman"><span style='font-size:12.0pt'>I
removed .dmat files from the SCF and the output files before starting with the
calculations for f-states.<o:p></o:p></span></font></p>

<p class=MsoNormal dir=LTR style='text-align:left;direction:ltr;unicode-bidi:
embed'><font size=3 face="Times New Roman"><span style='font-size:12.0pt'><o:p>&nbsp;</o:p></span></font></p>

<p class=MsoNormal dir=LTR style='text-align:left;direction:ltr;unicode-bidi:
embed'><font size=3 face="Times New Roman"><span style='font-size:12.0pt'>Please
help me to solve this problem.<o:p></o:p></span></font></p>

<p class=MsoNormal dir=LTR style='text-align:left;direction:ltr;unicode-bidi:
embed'><font size=3 face="Times New Roman"><span style='font-size:12.0pt'><o:p>&nbsp;</o:p></span></font></p>

<p class=MsoNormal dir=LTR style='text-align:left;direction:ltr;unicode-bidi:
embed'><font size=3 face="Times New Roman"><span style='font-size:12.0pt'>Thank
you in advance,<o:p></o:p></span></font></p>

<p class=MsoNormal dir=LTR style='text-align:left;direction:ltr;unicode-bidi:
embed'><font size=3 face="Times New Roman"><span style='font-size:12.0pt'>Angelica.<o:p></o:p></span></font></p>

<p class=MsoNormal dir=LTR style='text-align:left;direction:ltr;unicode-bidi:
embed'><font size=3 face="Times New Roman"><span style='font-size:12.0pt'><o:p>&nbsp;</o:p></span></font></p>

<p class=preformattedtext dir=LTR><font size=2 face="Nimbus Mono L"><span
style='font-size:10.0pt'><o:p>&nbsp;</o:p></span></font></p>

<p class=preformattedtext dir=LTR><font size=2 face="Nimbus Mono L"><span
style='font-size:10.0pt'><o:p>&nbsp;</o:p></span></font></p>

<p class=preformattedtext dir=LTR style='margin-bottom:14.15pt'><font size=2
face="Nimbus Mono L"><span style='font-size:10.0pt'><o:p>&nbsp;</o:p></span></font></p>

<p class=MsoNormal dir=LTR style='text-align:left;direction:ltr;unicode-bidi:
embed'><font size=2 face=Arial><span style='font-size:10.0pt;font-family:Arial'>&nbsp;<o:p></o:p></span></font></p>

<p class=MsoNormal dir=LTR style='text-align:left;direction:ltr;unicode-bidi:
embed'><font size=2 face=Arial><span style='font-size:10.0pt;font-family:Arial'><o:p>&nbsp;</o:p></span></font></p>

<p class=MsoNormal dir=LTR style='text-align:left;direction:ltr;unicode-bidi:
embed'><font size=2 face=Arial><span style='font-size:10.0pt;font-family:Arial'><o:p>&nbsp;</o:p></span></font></p>

</div>

</body>

</html>