 This thread is locked.Only browsing is available.
 Top Page > Browsing LCAO coefficients in real spherical harmonics, trasform into spherical harmonics Date: 2017/04/17 11:16 Name: P Hello Prof. T. Ozaki,I am looking at the LCAO coefficients. It is not clear to me which convention is chosen for the real spherical harmonics to be written as spherical harmonics. I would like to know the convention for l=1,2,3. Page: Re: LCAO coefficients in real spherical harmonics, trasform into spherical harmonics ( No.1 ) Date: 2017/04/18 00:01 Name: Artem Pulkin Spherical Harmonics come in the following order:"s","px","py","pz","dz2","d(x2-y2)","dxy","dxz","dyz","fz3","fxz2","fyz2","fz(x2-y2)","fxyz","fx(x2-3y2)","fy(3x2-y2)",The corresponding expressions are here:https://en.wikipedia.org/wiki/Table_of_spherical_harmonics#Real_spherical_harmonics  Re: LCAO coefficients in real spherical harmonics, trasform into spherical harmonics ( No.2 ) Date: 2017/04/19 07:08 Name: P Hello Artem Pulkin,Thank you for your answer. I understand how real spherical harmonics are simply a unitary transformation of regular spherical harmonics. However the relative sings are not uniquely defined. I have seen many cases where, depending on author, real spherical harmonics will differ by a minus sing in there transformation into regular spherical harmonics. I was interested in the exact definition used by the code.Except if I misunderstood your answer, and you mean the code uses the definition in the wiki page.  Re: LCAO coefficients in real spherical harmonics, trasform into spherical harmonics ( No.3 ) Date: 2017/04/19 18:04 Name: Artem Pulkin Well, the phase of a wavefunction is not uniquely defined so your request is a bit misleading. Let me formulate it this way: I was able to reproduce the overlap matrix exactly by using notation from the wikipedia page with two exceptions:"dxy" comes with the minus sign;"fxyz" comes with the minus sign.This information is sufficient to perform rotations you are looking for. If this is still not enough for you here is the corresponding code from OpenMX: /* Angular */ siQ = sin(Q); coQ = cos(Q); siP = sin(P); coP = cos(P); for (L0=0; L0<=Spe_MaxL_Basis[wan]; L0++){ if (L0==0){ AF = 0.282094791773878; } else if (L0==1){ dum = 0.48860251190292*siQ; AF = dum*coP; AF = dum*siP; AF = 0.48860251190292*coQ; } else if (L0==2){ dum1 = siQ*siQ; dum2 = 1.09254843059208*siQ*coQ; AF = 0.94617469575756*coQ*coQ - 0.31539156525252; AF = 0.54627421529604*dum1*(1.0 - 2.0*siP*siP); AF = 1.09254843059208*dum1*siP*coP; AF = dum2*coP; AF = dum2*siP; } else if (L0==3){ AF = 0.373176332590116*(5.0*coQ*coQ*coQ - 3.0*coQ); AF = 0.457045799464466*coP*siQ*(5.0*coQ*coQ - 1.0); AF = 0.457045799464466*siP*siQ*(5.0*coQ*coQ - 1.0); AF = 1.44530572132028*siQ*siQ*coQ*(coP*coP-siP*siP); AF = 2.89061144264055*siQ*siQ*coQ*siP*coP; AF = 0.590043589926644*siQ*siQ*siQ*(4.0*coP*coP*coP - 3.0*coP); AF = 0.590043589926644*siQ*siQ*siQ*(3.0*siP - 4.0*siP*siP*siP); } else if (4<=L0){  Re: LCAO coefficients in real spherical harmonics, trasform into spherical harmonics ( No.4 ) Date: 2017/04/20 17:07 Name: T. Ozaki Hi, The transformation matrices from the complex to the real can be found in Set_Comp2Real() in SetPara_DFT.c. Regards, TO Page: