Appendix J Phase convention of real spherical harmonics

In quantum chemistry applications, an Condon-Shortley phase can be added to the real spherical harmonics (RSH) based on the developers’ choice.

If, the commonly used Condon-Shortley phase is added to the RSH, the definition of RSH becomes

Ylm={(1)m2Klmcos(mϕ)Plm(cosθ),m>0(1)m2Klmsin(|m|ϕ)Pl|m|(cosθ),m<0KlmPlm(cosθ),m=0 (J.4)

where Plm is the associated Legendre polynomials without the Condon-Shortley phase

Plm(cosθ)=(sinθ)mdmd(cosθ)m(Pl(cosθ)), (J.5)

and the coefficient

Klm=(2l+1)(l|m|)!4π(l+|m|)! (J.6)

The RSH in the FHI-aims code adopt a partly Condon-Shortley phase. Thus, the definition of RSH becomes

Ylm={(1)m2Klmcos(mϕ)Plm(cosθ),m>02Klmsin(|m|ϕ)Pl|m|(cosθ),m<0KlmPlm(cosθ),m=0 (J.10)

The code related to these definition can be found in file basis_sets/increment_ylm_original.f90 or basis_sets/SHEval.f90. In some subroutines, the real spherical harmonics are generated using external/ylm_real.f90. The definition of RSH therein seems a bit different

Ylm={2ReY~lm,m>02ImY~lm,m<0Y~lm,m=0 (J.14)

where Y~lm is the complex spherical harmonics (CSH)

Y~lm(θ,ϕ)=(1)m2l+14π(lm)!(l+m)!Plm(cosθ)eimϕ (J.15)

It can be shown that the two conventions are equivalent. it is trivial for m=0. For m>0,

(1)m2Klmcos(mϕ)Plm(cosθ)= (1)m2Re[eimϕ]KlmPlm(cosθ)
= 2Re[(1)meimϕKlmPlm(cosθ)]
= 2ReY~lm

For m<0

2Klmsin(|m|ϕ)Pl|m|(cosθ)
= 2Im[ei|m|ϕ]KlmPl|m|(cosθ)
= 2Im[eimϕ](2l+1)(l+m)!4π(lm)!Plm(cosθ)
= 2Im[eimϕ(2l+1)(lm)!4π(l+m)!(1)mPlm(cosθ)]
= 2ImY~lm

where in the third equality we have used the following property

Plm=(1)m(lm)!(l+m)!Plm. (J.16)

With the property of CSH

Y~lm=(1)mY~lm, (J.17)

we can derive another equivalent expression from Eq. (J.14)

Ylm={12[Y~lm+(1)mY~lm],m>0i2[Y~lm(1)mY~lm],m<0Y~lm,m=0 (J.21)

leading to the representation matrix of FHI-aims RSH in CSH. This can be useful to connect between RSHs from FHI-aims and other codes.