<DIV>Dear Prof. Blaha,</DIV>
<DIV>&nbsp;</DIV>
<DIV>Thanks a lot to your patient explanations. However, I still have problems with it. According to Table7.39 in the UG, for example, there are four possible choices for l=1 with point group "M", (+-l,l-2m), they are (1,1),(-1,1) for m=0 and (1,-1),(-1,-1) for m=1. It is obvious that only (1,1) and (-1,1) are listed in the case.in2 file. I don't understand why (1,-1)(-1,-1) are not included? I think I have to understand this first before I try to&nbsp; translate the lm definitions up to higher combinations given&nbsp;in Table7.39.</DIV>
<DIV>&nbsp;</DIV>
<DIV>Best regards,</DIV>
<DIV>&nbsp;</DIV>
<DIV>Stargmoon</DIV>
<DIV>&nbsp;</DIV>
<DIV>P.S. Could you please specify what kind of values "l,m" (symbols in Table7.39)&nbsp;can take? (l=0,1,2,...;m=0,1,2,...l or l=0,1,2...;m=-l,-l+1,...0,...l or something else?)<BR><BR><B><I>Peter Blaha &lt;pblaha@zeus.theochem.tuwien.ac.at&gt;</I></B> wrote:</DIV>
<BLOCKQUOTE class=replbq style="PADDING-LEFT: 5px; MARGIN-LEFT: 5px; BORDER-LEFT: #1010ff 2px solid">As I mentioned before, the negative l is just a method to specify <BR>one of the four possible linear combinations.<BR>To expand the LM list in case.in2, just type in manually the higher <BR>combinations. In the cubic cases (positive atomic numbers in the <BR>case.struct file) the order of LM pairs matters and higher LM values are <BR>given explicitely in the UG. For other pointgroups either take the <BR>"intelligence test" and continue the series in analogy to the l=0-6 <BR>values, or use the Table in the UG and translate the lm definitions given <BR>there. <BR>&gt; Yes, for the spherical harmonics, "l" never takes negative values. <BR>But I don't understand why there are "plus and minus" symbols shown before "l" in <BR>Table 7.39 in the manual, and the real LM list generated used in WIEN2k has negative "l" <BR>(see Table 7.38)? Is it just an indicator of the possible combin!
 ations?
 <BR>How can I add the extra LM combinations beyond maximum value 6 by hand? <BR>Does the order of the LM combinations matter in the calculation?<BR>&gt; <BR>&gt; Best regards,<BR>&gt; <BR>&gt; Stargmoon<BR>&gt; <BR>&gt; Peter Blaha <PBLAHA@ZEUS.THEOCHEM.TUWIEN.AC.AT>wrote:<BR>&gt; Te charge density is expanded in "real spherical harmonics". These <BR>&gt; functions can be constructed from the (more common) "complex spherical <BR>&gt; harmonics" by various linear combinations. As you all know, a complex <BR>&gt; spherical harmonics is Y_lm = P_lm(cos theta) . exp(i m phi)<BR>&gt; <BR>&gt; The LM list defines linear combinations, which produce the required "real <BR>&gt; harmonics".<BR>&gt; <BR>&gt; Obviously, when m=0, the imaginary part (exp(im.phi) is always zero, thus <BR>&gt; there is only ONE (l 0) pair, but never a negative l.<BR>&gt; <BR>&gt; Of course also in the general case, negative l's do not exists, but the <BR>&gt; notation of the 4 possible combinations to pro!
 duce real
 harmonics are <BR>&gt; indicated by (l m),(-l m),(l -m) and (-l -m). The possible linear <BR>&gt; combinations are<BR>&gt; <BR>&gt; Y_lm + Y_l-m<BR>&gt; Y_lm - Y_l-m<BR>&gt; (Y_lm - Y_l-m)*i<BR>&gt; (Y-lm - Y_lm-m)*i<BR>&gt; <BR>&gt; <BR>&gt; &gt; I think the lm list looks okay, for (+-l,l-2m) the sets should be (1,1),(-1,1) <BR>&gt; &gt; for m=0 and for m=1 gives (1,-1) and (-1,-1)... I'm not an expert but I've <BR>&gt; &gt; negative m values are not valid and (-1,0), (-2,0),...(-n,0) never seem to be <BR>&gt; &gt; included in the lm-list, though I have no idea why. Maybe someone can <BR>&gt; &gt; explain?<BR>&gt; &gt; <BR>&gt; &gt; Javier<BR>&gt; &gt; <BR>&gt; &gt; &gt;===== Original Message From "Jorissen Kevin" =====<BR>&gt; &gt; &gt;Hello,<BR>&gt; &gt; &gt;the lm-list is usually read and used explicitly, which sounds like the order <BR>&gt; &gt; doesn't matter, but I'm not sure (perhaps the last value is used as 'maximum' <BR>&gt; &gt; somehow)? Better play safe and pu!
 t them in
 the right order, unless someone <BR>&gt; &gt; else tells us it's okay not to.<BR>&gt; &gt; &gt;How I got them : 1/ I remember finding something about it in the old archives <BR>&gt; &gt; which gave some examples; 2/ for pointgroup "1" it's not very difficult, <BR>&gt; &gt; obviously ;-).<BR>&gt; &gt; &gt;As for the second atom, it would surprise me to find a mistake in the l=1 <BR>&gt; &gt; group, since I definitely took that from wien and only added l=7 tot 10 <BR>&gt; &gt; myself. But who knows, it's a file I haven't used in a while.<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt;Kevin Jorissen<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt;EMAT - Electron Microscopy for Materials Science <BR>&gt; &gt; (http://webhost.ua.ac.be/emat/)<BR>&gt; &gt; &gt;Dept. of Physics<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt;UA - Universiteit Antwerpen<BR>&gt; &gt; &gt;Groenenborgerlaan 171<BR>&gt; &gt; &gt;B-2020 Antwerpen<BR>&gt; &gt; &gt;Belgium<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt;!
 tel +32 3
 2653249<BR>&gt; &gt; &gt;fax + 32 3 2653257<BR>&gt; &gt; &gt;e-mail kevin.jorissen@ua.ac.be<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt;________________________________<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt;Van: wien-admin@zeus.theochem.tuwien.ac.at namens stargmoon<BR>&gt; &gt; &gt;Verzonden: do 20-1-2005 3:17<BR>&gt; &gt; &gt;Aan: wien@zeus.theochem.tuwien.ac.at<BR>&gt; &gt; &gt;Onderwerp: RE: [Wien] LMMAX again<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt;Dear Kevin,<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt;First of all, thanks a lot for your reply. But I still have questions.<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt;In your example, the second sort has point group "m", so its LM combination <BR>&gt; &gt; should take (+-l,l-2m). As "l" takes value 1, there are supposed to have four <BR>&gt; &gt; sets of LM, that is, (1,0),(1,1),(-1,0) and (-1,1), I think. However I can't <BR>&gt; &gt; find (-1,0) in the corresponding line (second line) in your .in2 file. Am I <BR>&gt; !
 &gt;
 right?<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt;By the way, how do you get the LM combinations up to L=10? Do you just change <BR>&gt; &gt; the default value 6 to 10 in the source code directly? If I add it by hand, <BR>&gt; &gt; does the order of the LM list matter?<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt;Looking forward to your reply!<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt;Best,<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt;Stargmoon<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt;Jorissen Kevin wrote:<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt; Did you enter them in 'fixed format'?<BR>&gt; &gt; &gt; Here's an example for TiO:<BR>&gt; &gt; &gt; 0 0 4 0 4 4 6 0 6 4 8 0 8 4 8 8<BR>&gt; &gt; &gt; 0 0 4 0 4 4 6 0 6 4 8 0 8 4 8 8<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt; And this is a (working) example for KC60 :<BR>&gt; &gt; &gt; tsm2.cmi.ua.ac.be&gt; more vijf.in2<BR>&gt; &gt; &gt; TOT (TOT,FOR,QTL,EFG,FERMI)<BR>&gt; &gt; &gt; -9.0 249.0 0.50 0.05 EMIN, NE,
 ESEPERMIN, ESEPER0<BR>&gt; &gt; &gt; TETRA 0.000 (GAUSS,ROOT,TEMP,TETRA,ALL eval)<BR>&gt; &gt; &gt; 0 0 4 0 4 4 6 0 6 4 6 2 6 6 8 0 8 4 8 8 10 0 10 2 10 4 10 6 10 8 1010<BR>&gt; &gt; &gt; 0 0 1 1 -1 1 2 0 2 2 -2 2 3 1 -3 1 3 3 -3 3 4 0 4 2 -4 2 4 4 -4 4 5 1<BR>&gt; &gt; &gt; -5 1 5 3 -5 3 5 5 -5 5 6 0 6 2 -6 2 6 4 -6 4 6 6 -6 6 -7 7 -7 5 -7 3 -7 1<BR>&gt; &gt; &gt; 7 1 7 3 7 5 7 7 -8 8 -8 6 -8 4 -8 2 8 0 8 2 8 4 8 6 8 8 -9 9 -9 7 -9 5<BR>&gt; &gt; &gt; -9 3 -9 1 9 1 9 3 9 5 9 7 9 9-1010-10 8-10 6-10 4-10 2 10 0 10 2 10 4 10 6<BR>&gt; &gt; &gt; 10 8 1010<BR>&gt; &gt; &gt; 0 0 1 0 1 1 -1 1 2 0 2 1 -2 1 2 2 -2 2 3 0 3 1 -3 1 3 2 -3 2 3 3 -3 3<BR>&gt; &gt; &gt; 4 0 4 1 -4 1 4 2 -4 2 4 3 -4 3 4 4 -4 4 5 0 5 1 -5 1 5 2 -5 2 5 3! -5 3<BR>&gt; &gt; &gt; 5 4 -5 4 5 5 -5 5 6 0 6 1 -6 1 6 2 -6 2 6 3 -6 3 6 4 -6 4 6 5 -6 5 6 6<BR>&gt; &gt; &gt; -6 6 -7 7 -7 6 -7 5 -7 4 -7 3 -7 2 -7 1 7 0 7 1 7 2 7 3 7 4 7 5 7 6 7 7<BR>&gt; &gt; &gt; -8 8 -8 7 -8 6 -8 5 -8 4 -8 3 -8 2 -8 1 8 0 8 1 8 2 8!
  3 8 4 8
 5 8 6 8 7<BR>&gt; &gt; &gt; 8 8 -9 9 -9 8 -9 7 -9 6 -9 5 -9 4 -9 3 -9 2 -9 1 9 0 9 1 9 2 9 3 9 4 9 5<BR>&gt; &gt; &gt; 9 6 9 7 9 8 9 9-1010-10 9-10 8-10 7-10 6-10 5-10 4-10 3-10 2-10 1 10 0 10 1<BR>&gt; &gt; &gt; 10 2 10 3 10 4 10 5 10 6 10 7 10 8 10 9 1010<BR>&gt; &gt; &gt; 0 0 1 0 1 1 -1 1 2 0 2 1 -2 1 2 2 -2 2 3 0 3 1 -3 1 3 2 -3 2 3 3 -3 3<BR>&gt; &gt; &gt; 4 0 4 1 -4 1 4 2 -4 2 4 3 -4 3 4 4 -4 4 5 0 5 1 -5 1 5 2 -5 2 5 3 -5 3<BR>&gt; &gt; &gt; 5 4 -5 4 5 5 -5 5 6 0 6 1 -6 1 6 2 -6 2 6 3 -6 3 6 4 -6 4 6 5 -6 5 6 6<BR>&gt; &gt; &gt; -6 6 -7 7 -7 6 -7 5 -7 4 -7 3 -7 2 -7 1 7 0 7 1 7 2 7 3 7 4 7 5 7 6 7 7<BR>&gt; &gt; &gt; -8 8 -8 7 -8 6 -8 5 -8 4 -8 3 -8 2 -8 1 8 0 8 1 8 2 8 3 8 4 8 5 8 6 8 7<BR>&gt; &gt; &gt; 8 8 -9 9 -9 8 -9 7 -9 6 -9 5 -9 4 -9 3 -9 2 -9 1 9 0 9 1 9 2 9 3 9 4 9 5<BR>&gt; &gt; &gt; 9 6 9 7 9 8 9 9-1010-10 9-10 8-10 7-10 6-10 5-10 4-10 3-10 2-10 1 10 0 10 1<BR>&gt; &gt; &gt; 10 2 10 3 10 4 10 5 10 6 10 7 ! 10 8 10 9 1010<BR>&gt; &gt; &gt; 20. GMAX<BR>&gt;!
  &gt;
 &gt; FILE FILE/NOFILE write recprlist<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt; where the point symmetry of the four positions is the following (from <BR>&gt; &gt; case.outputsgroup) :<BR>&gt; &gt; &gt; Sort number: 1<BR>&gt; &gt; &gt; Names of point group: m-3 2/m-3 Th<BR>&gt; &gt; &gt; Sort number: 2<BR>&gt; &gt; &gt; Names of point group: m m Cs<BR>&gt; &gt; &gt; Sort number: 3<BR>&gt; &gt; &gt; Names of point group: 1 1 C1<BR>&gt; &gt; &gt; Sort number: 3<BR>&gt; &gt; &gt; Names of point group: 1 1 C1<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt; The exact number of spaces is important! In my example for KC60, lines end <BR>&gt; &gt; only with 1010 (other 'new lines' due to copy/paste, sorry).<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt; Hope this will do,<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt; Kevin Jorissen<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt; EMAT - Electron Microscopy for Materials Science <BR>&gt; &gt; (http://webhost.ua.ac.be/emat/)<BR>&gt; &gt; &gt; Dept. of Physics<BR>&gt; &gt; &gt;<BR>&gt;!
  &gt;
 &gt; UA - Universiteit Antwerpen<BR>&gt; &gt; &gt; Groenenborgerlaan 171<BR>&gt; &gt; &gt; B-2020 Antwerpen<BR>&gt; &gt; &gt; Belgium<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt; tel +32 3 2653249<BR>&gt; &gt; &gt; fax + 32 3 2653257<BR>&gt; &gt; &gt; e-mail kevin.jorissen@ua.ac.be<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt; ________________________________<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt; Van: wien-admin@zeus.theochem.tuwien.ac.at namens stargmoon<BR>&gt; &gt; &gt; Verzonden: wo 19-1-2005 18:27<BR>&gt; &gt; &gt; Aan: Wien@zeus.theochem.tuwien.ac.at<BR>&gt; &gt; &gt; Onderwerp: [Wien] LMMAX again<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt; Dear WIEN community,<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt; A couple of days ago, I asked a question about LMMAX here and Kevin's answer <BR>&gt; &gt; helped me a lot. However, when I try to add LM list manually, I meet problems. <BR>&gt; &gt; For example, with a site symmetry "2/M", as listed in Table7.39 in usersguide, <BR>&g!
 t; &gt;
 the LM combination should be (+-2l,2m). So I try to set l=0,1,2,3, (m=0,1,...l <BR>&gt; &gt; for each l) in order to reproduce the default lmmax up to 6, and I get 19 LM <BR>&gt; &gt; combinations manually. However, I only got 16 LM combinations after I run <BR>&gt; &gt; "SYMMETRY". Those three (-2,0) (-4,0) and (-6,0) do not show up in the LM list <BR>&gt; &gt; obtained from SYMMETRY. What's the problem here? Or I did something wrong? <BR>&gt; &gt; Your reply will be highly appreciated.<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt; Best,<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt; Stargmoon<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt; ________________________________<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt; Do you Yahoo!?<BR>&gt; &gt; &gt; Yahoo! Mail - You care! about security. So do we.<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt; &gt; ATTACHMENT part 2 application/ms-tnef name=winmail.dat<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt;<BR>&gt; &gt; &gt;________________________________<BR>&gt; &gt; &gt;<B!
 R>&gt;
 &gt; &gt;Do you Yahoo!?<BR>&gt; &gt; &gt;The all-new My Yahoo! - What will yours do?<BR>&gt; &gt; <BR>&gt; &gt; <BR>&gt; &gt; _______________________________________________<BR>&gt; &gt; Wien mailing list<BR>&gt; &gt; Wien@zeus.theochem.tuwien.ac.at<BR>&gt; &gt; http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien<BR>&gt; &gt; <BR>&gt; <BR>&gt; <BR>&gt; P.Blaha<BR>&gt; --------------------------------------------------------------------------<BR>&gt; Peter BLAHA, Inst.f. Materials Chemistry, TU Vienna, A-1060 Vienna<BR>&gt; Phone: +43-1-58801-15671 FAX: +43-1-58801-15698<BR>&gt; Email: blaha@theochem.tuwien.ac.at WWW: http://info.tuwien.ac.at/theochem/<BR>&gt; --------------------------------------------------------------------------<BR>&gt; <BR>&gt; _______________________________________________<BR>&gt; Wien mailing list<BR>&gt; Wien@zeus.theochem.tuwien.ac.at<BR>&gt; http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien<BR>&gt; <BR>&gt; <BR>&gt;
 ---------------------------------<BR>&gt; Do you Yahoo!?<BR>&gt; Yahoo! Search presents - Jib Jab's 'Second Term'<BR><BR><BR>P.Blaha<BR>--------------------------------------------------------------------------<BR>Peter BLAHA, Inst.f. Materials Chemistry, TU Vienna, A-1060 Vienna<BR>Phone: +43-1-58801-15671 FAX: +43-1-58801-15698<BR>Email: blaha@theochem.tuwien.ac.at WWW: http://info.tuwien.ac.at/theochem/<BR>--------------------------------------------------------------------------<BR><BR>_______________________________________________<BR>Wien mailing list<BR>Wien@zeus.theochem.tuwien.ac.at<BR>http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien<BR></BLOCKQUOTE><p>
                <hr size=1>Do you Yahoo!?<br> 
Yahoo! Search presents - <a href="http://us.rd.yahoo.com/evt=30648/*http://movies.yahoo.com/movies/feature/jibjabinaugural.html">Jib Jab's 'Second Term'</a>