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