<html>
<head>
<style>
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
FONT-SIZE: 9pt;
FONT-FAMILY:Tahoma
}
</style>
</head>
<body class='hmmessage'>
<DIV class="ReadMsgBody BorderTop">
<DIV class=ExternalClass><FONT color=#339966>Dear Prof. Blaha;</FONT></DIV>
<DIV class=ExternalClass><FONT color=#339966>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;I am sorry that I did not express myself clearly&nbsp;the first time.</FONT><BR><FONT color=#339966>&nbsp;&nbsp;&nbsp;&nbsp; thanks for your kind reply last time, but I am so sorry to ask you somewhat entry-level </FONT><BR><FONT color=#339966>queations again.</FONT><BR><FONT color=#339966>&nbsp;&nbsp;&nbsp;&nbsp; now I guess the charge density in MT sphere is expanded in the basis of </FONT><BR><FONT color=#339966>real spherical harmonic.but I still have a question as the following:&nbsp;</FONT><BR><FONT color=#339966>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;I find some&nbsp;different treatments of the&nbsp;multipole moments&nbsp;between cubic and</FONT><BR><FONT color=#339966>non-cubic cases in lapw0.F&nbsp;when I read the source code file in SRC_lapw0.</FONT><BR><FONT color=#339966>Seemingly for the non-cubic case the&nbsp;coefficients of density&nbsp;is straightly transformed </FONT><BR><FONT color=#339966>f!
 rom real spherical haromonic representation to&nbsp;complex sphereical harmonic representation,</FONT><BR><FONT color=#339966>however, for cubic case, it go through some complex transformations distinctly from non-cubic&nbsp;</FONT><BR><FONT color=#339966>transformation as the code below indicats. </FONT><FONT color=#ff0000>I&nbsp;feel that for&nbsp;cubic case, the&nbsp;density is&nbsp;not<BR>expanded&nbsp;according to&nbsp;real spherical harmonics,but I don't know exactly how it is expanded for <BR>cubic cases,also,&nbsp;it is not expanded in the another sphere&nbsp;set of real&nbsp;spherical harmonics K(l,j),</FONT></DIV>
<DIV class=ExternalClass><FONT color=#ff0000>because based on this following code, there exists term&nbsp;like a1**2,a1*a2, a1*a3 and so on, not </FONT></DIV>
<DIV class=ExternalClass><FONT color=#ff0000>the&nbsp;case of&nbsp;directly transfroming K(l,j) to Y(L,M).K(l,j)&nbsp;is linear combination of y(l,m)&nbsp;<BR></FONT>&nbsp;<BR>&nbsp;&nbsp;&nbsp;<FONT color=#0000ff>## for non-cubic case</FONT><BR><FONT color=#0000ff>&nbsp;&nbsp;&nbsp;QQ(LLMM,JATOM)=QQ(LLMM,JATOM)*fc(llmm,jatom)</FONT><BR><FONT color=#0000ff>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</FONT><BR><FONT color=#0000ff>&nbsp;&nbsp;## for cubic case&nbsp;</FONT><BR><FONT color=#0000ff>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; c1=c_kub(ABS(lm(1,i,jatom)),lm(2,i,jatom))<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; c2=c_kub(ABS(lm(1,i,jatom)),lm(2,i,jatom)+4)<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; c3=c_kub(ABS(lm(1,i,jatom)),lm(2,i,jatom)+8)<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; qq(i,jatom)=a*c1*c1/sq1 + b*c1*c2/sq1 + &amp;<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n!
 bsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; c*c1*c3/sq1<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; qq(i+1,jatom)=a*c1*c2/sqrt2 + b*c2*c2/sqrt2&nbsp; &amp;<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; + c*c2*c3/sqrt2<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; qq(i+2,jatom)=a*c1*c3/sqrt2 + b*c2*c3/sqrt2&nbsp; &amp;<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; + c*c3*c3/sqrt2</FONT><BR>&nbsp;<BR>&nbsp;&nbsp;<FONT color=#339966> I have read the recommended papar in Userguide in lapw2 program about the lattice harmonic,</FONT><BR><FONT color=#339966>but the still cannot figure out this question,so, could you give me some your idea about this issue?</FONT><BR><FONT color=#339966>thanks very much !</FONT><BR><FONT color=#339966>&nbsp;&nbsp; </FONT><BR><FONT color=#339966>Regards!</FONT><BR><FONT color=#339966>Wenmei Ming</FONT><BR>&nb!
 sp;&nbsp;&nbsp; &nbsp;&nbsp;<BR>&nbsp;<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbs
p;&nbsp;&nbsp;<BR>&nbsp;<BR><BR>&gt; Date: Mon, 17 Dec 2007 21:21:10 +0100<BR>&gt; From: Peter Blaha &lt;pblaha@theochem.tuwien.ac.at&gt;<BR>&gt; Subject: Re: [Wien] charge density in CLM files<BR>&gt; To: A Mailing list for WIEN2k users &lt;wien@zeus.theochem.tuwien.ac.at&gt;<BR>&gt; Message-ID: &lt;4766DA36.9040605@theochem.tuwien.ac.at&gt;<BR>&gt; Content-Type: text/plain; charset=GB2312<BR>&gt; <BR>&gt; 1) I hope you agree: the electron density (psi^* psi) must be a "real"<BR>&gt; function.<BR>&gt; The usual definition of spherical harmonics gives them as a complex<BR>&gt; function. However, one can always form some linear combination and<BR>&gt; obtain "real"-spherical harmonics. This is what is used here. More<BR>&gt; specifically, we use "lattice" harmonics, i.e. a symmetry adapted<BR>&gt; selection of some spherical harmonics.<BR>&gt; <BR>&gt; 2) This "syntax" is connected with the selection of "lattice harmonics".<BR>&gt; Of course a "negative"-l does not mean that !
 l is negative, but indicates<BR>&gt; a particular form of linear combination. For more details see the UG or<BR>&gt; the original papers (cited in the UG).<BR>&gt; <BR>&gt; &gt; 1) why the density in MT sphere given in case.clmXX files are all <BR>&gt; &gt; real values<BR>&gt; &gt; rather than complex values for all radial points. Based on the <BR>&gt; &gt; userguide, these values<BR>&gt; &gt; represent the coefficients of spherical harmonic <BR>&gt; &gt; expansion,according to this,<BR>&gt; &gt; these coefficients should be complex values not real ones. so ,I <BR>&gt; &gt; want to know<BR>&gt; &gt; what the so-called "density" for MT sphere in clm file really <BR>&gt; &gt; mean !!<BR>&gt; &gt; <BR>&gt; &gt; 2) occasionaly,some of the lattice harmonic indexes for L listed in <BR>&gt; &gt; case.in2 are negative,<BR>&gt; &gt; and the possible LM listed in the table in the LAPW2 chapter in <BR>&gt; &gt; userguide <BR>&gt; &gt; also indicate that L can negative,furthermore, I r!
 ead the <BR>&gt; &gt; related source codes<BR>&gt; &gt; and find that 
some programs are firmly connected with this <BR>&gt; &gt; issue. thus, I am puzzled about this,because L should not be negative (&gt;=0)<BR>&gt; &gt; <BR>&gt; &gt; I am eager to figure out what is going on about the two issues, <BR>&gt; &gt; can you give me some ideas<BR>&gt; &gt; or let me know the details of both of them.<BR>&gt; &gt; <BR>&gt; &gt; thanks very much !<BR>&gt; &gt; <BR>&gt; &gt; Ming Wenmei<BR>&gt; &gt; <BR><BR></DIV></DIV><br /><hr />轻松把Hotmail下载到本地,试试 Windows Live Mail。 <a href='http://get.live.com/wl/all' target='_new'>立即尝试!</a></body>
</html>