<div>Hi Dear User</div>
<div>I think you should use LSDA approximation for running LDA+U and after if you read the user guide you could see </div>
<div>That you shoud find the existing file and then constract some </div>
<div>information ,depend on your compound.</div>
<div>and after that you should run the SCF.(+ lO).</div>
<div>Bye<br><br>&nbsp;</div>
<div><span class="gmail_quote">On 6/1/06, <b class="gmail_sendername"><a href="mailto:zhylimin@sohu.com">zhylimin@sohu.com</a></b> &lt;<a href="mailto:zhylimin@sohu.com">zhylimin@sohu.com</a>&gt; wrote:</span>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">
<div>
<div>Dear Wien users &amp; developers, 
<p>I have recently done some calcuations on AlN doped with Er. LDA+U method 
<p>
<p>is used to treate Er 4f electrons. Here are my operations step by step: 
<p>
<p>1. A directory named AlN is created; 
<p>
<p>2. input-files such as AlN.struct, AlN.inst, AlN.indm and AlN.inorb is created. 
<p>
<p>3. init_lapw (GGA, no add-inversion and spin-polarized are selected in this step); 
<p>
<p>4. [hust@icesword AlN]$ runsp_lapw -orb 
<p>LAPW0 END 
<p>LAPW1 END 
<p>LAPW1 END 
<p>LAPW2 END 
<p>LAPW2 END 
<p>
<p>after LAPW2, SCF terminated without LAPWDM and LCORE...etc steps. I do know if my operations 
<p>
<p>
<p>is correct. Could you tell me the detail steps on how to us LDA+U method? 
<p>
<p>
<p>Thank you for your helpl 
<p>
<p>
<p>Regards, 
<p>
<p>
<p>Zhang 
<p>
<p>Appendix:some input-files: 
<p>AlN.indm: 
<p>------------------------------------------ 
<p>-9. Emin cutoff energy 
<p>1 number of atoms for which density matrix is calcul! ated 
<p>2 1 3 index of 1st atom, number of L's, L1 
<p>0 0 r-index, (l,s)index 
<p>----------------------------------------- 
<p>
<p>AlN.inorb: 
<p>----------------------------------------- 
<p>1 1 0 nmod, natorb, ipr 
<p>BROYD 0.3 BROYD/PRATT, mixing 
<p>2 1 3 iatom nlorb, lorb 
<p>1 nsic 0..AFM, 1..SIC, 2..HFM 
<p>2.59 0.7 U J (Ry) 
<p>--------------------------------------- 
<p>
<p>
<p>一步登天就是这么简单！ 
<p>
<p>
<p>
<p>
<p>
<hr size="1">

<table border="0">
<tbody>
<tr>
<td height="18"><a onclick="return top.js.OpenExtLink(window,event,this)" href="http://goto.sohu.com/goto.php3?code=mailadt-ta1" target="_blank"><font color="red">免费1G网络U盘，简简单单U来U去！</font></a> </td></tr>
<tr>
<td><br>
<tr>
<td height="18"><a onclick="return top.js.OpenExtLink(window,event,this)" href="http://goto.sohu.com/goto.php3?code=mailadt-ta" target="_blank">一步登天就是这么简单！</a> </td></tr></td></tr></tbody></table>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p>
<p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p></p>
</p></p></div></div><br>_______________________________________________<br>Wien mailing list<br><a onclick="return top.js.OpenExtLink(window,event,this)" href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at
</a><br><a onclick="return top.js.OpenExtLink(window,event,this)" href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br><br><br></blockquote>
</div><br>