<DIV>
<P>Dear Dr. Cottenier,</P><FONT face="Arial Unicode MS" size=2>
<P>Hello. Thank you for considering my question about LDA+U. By the sentence " I did whatever was said in the userguide! " </FONT><FONT lang=JA face="Arial Unicode MS" size=2>I did not mean that the userguide is not complete but I cannot understand it. Please excuse me for that. Here let me explain you ( for the special case of gadolinium) what I did.</FONT><FONT face="Arial Unicode MS" size=2>First I ran the scf and after that i</FONT><FONT face=Arial color=#333333 size=2>n order to run LDA+U, I generated the two files <U>Gd.indm</U> and <U>Gd.inorb</U> in the corresponding directory, and then changed the data as follows:</P></FONT><FONT face="Courier New" size=2>
<P>&nbsp;</P>
<P>Gd.indm</P>
<P>-9.&nbsp;&nbsp;&nbsp; Emin cutoff energy</P>
<P>1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; number of atoms for which density matrix is calculated</P>
<P>1 1 3&nbsp; index of 1st atom, number of L's, L1</P>
<P>0 0&nbsp;&nbsp;&nbsp; r-index, (l,s)index </P>
<P></P>
<P>&nbsp;</P>
<P>Gd.inorb</P>
<P>1 1 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nmod, natorb, ipr</P>
<P>BROYD 0.2&nbsp;&nbsp;&nbsp; BROYD/PRATT, mixing</P>
<P>1 1 3&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; iatom nlorb, lorb</P>
<P>1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;nsic 0..AFM, 1..SIC, 2..HFM</P>
<P>0.7 0.05&nbsp;&nbsp;&nbsp;&nbsp; U J (Ry)</P>
<P>The values written for U and J are not the exact ones.Next according to the instructions given in the userguide I tried to run the program by (for inctance)the following command:</P>
<P></FONT><FONT face="Times New Roman">&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; runsp_lapw -cc 0.0001 -orb </P>
<P align=justify>and the program just stopped in first cycle of <U>scf</U> giving empty lapwdm. I also ran x lapwup/dn. I cannot find where I'm wrong.</P>
<P align=justify>Two more questions:</P>
<P align=justify>1. Is it O.K. we say that since in Gd L=0 then because of "orbital angular momentum quenching" the spin-orbit coupling does not play a major role in calculations?</P>
<P align=justify>2. When including SO it asks for the case.struture to be modified. How should this be done?</P>
<P align=justify>Thank you very much again.</P>
<P align=justify>Ali Afshar</P>
<P align=justify><A href="mailto:afshar_130@yahoo.com">afshar_130@yahoo.com</A></P></FONT></DIV><p>
        
                <hr size=1>Do you Yahoo!?<br> 
Check out the new Yahoo! Front Page. <a href="http://www.yahoo.com">www.yahoo.com</a>