deepmodeling / abacus-develop

An electronic structure package based on either plane wave basis or numerical atomic orbitals.
http://abacus.ustc.edu.cn
GNU Lesser General Public License v3.0
164 stars 128 forks source link

Confused about Condon–Shortley phase related problem #4908

Closed Franklalalala closed 1 month ago

Franklalalala commented 1 month ago

Details

In the doc of abacus, it mentions that, the formula in wiki lacks the Condon–Shortley phase, therefore abacus team added the Condon–Shortley phase, and used this form in the lcao orbitals of ABACUS.

Here are two questions:

  1. It has been mentioned in another wiki, the real form of spherical harmonics can be infered from the complex form. In emphasis, quote, 'The Condon–Shortley phase convention is used here for consistency.' Is there a possibility that the real forms of spherical harmonics listed in the first wiki do contain the Condon–Shortley phase?

  2. Regardless of the correctness of the wiki form, does abacus team add another (-1)^m in their lcao orbitals? Since most developers follow the same form with wiki, such as pyscf, see doc. If abacus is different from wiki, we will add (-1)^m for communincation with pyscf. However, as far as I searched, there is no additional consideration for this, except a mention in test scipt. Here

Have you read FAQ on the online manual http://abacus.deepmodeling.com/en/latest/community/faq.html

Task list for Issue attackers (only for developers)

jinzx10 commented 1 month ago

Thanks for your question. Our documentation is indeed somewhat confusing and the relevant section will be revised to clarify this topic in a further PR. To answer your question 2, yes the real spherical harmonics adopted by abacus do differ from those listed in wikipedia by a factor of (-1)^m. To be more precise, real spherical harmonics in abacus takes the following form, in comparison with those in wiki:

image image

where $Y{lm}$ and $Y{l}^{m}$ represent real and conventional (complex) spherical harmonics respectively, and the complex spherical harmonics are defined to carry the Condon-Shortley phase. The order of m in abacus follows 0, +1, -1, +2, -2, ..., +l, -l, which is different from pyscf as well.

Question 1 seems to be a terminology issue. To my knowledge, Condon-Shortley phase is introduced to the complex spherical harmonics in order to simplify the treatment with ladder operators, while the sign convention of real spherical harmonics is a completely different matter motivated by computational consideration. The sign convention adopted in wiki does have an advantage that the resulting functions are signless in Cartesian coordinate (i.e., px/py/pz are simply x/r, y/r, z/r, while in abacus they are -x/r, -y/r, z/r), but I'm not sure if the extra (-1)^m are "Condon-Shortley phase" (in fact they act to cancel the original Condon-Shortley phase in the complex spherical harmonics). I think what wiki means by "The Condon–Shortley phase convention is used here for consistency" is probably that the complex spherical harmonics on the r.h.s. contains the phase, not the real spherical harmonics on the l.h.s. It's possible that I misunderstood wiki and missed some literatures on this terminology; I will appreciate if you could provide your feedback.

The convention adopted by wiki is indeed very popular, but the other convention can also be seen in a number of literatures, e.g.,

Homeier, H. H., & Steinborn, E. O. (1996). Some properties of the coupling coefficients of real spherical harmonics and their relation to Gaunt coefficients. Journal of Molecular Structure: THEOCHEM, 368, 31-37.

Hope this reply helps.

Franklalalala commented 1 month ago

Your reply has greatly clearified my confusion.

I have no question about the difference between abacus and wiki now.

For the terminology, here are some of my understandings to help with the discussion.

  1. The complex form of spherical harmonics takes the consideration of the Condon–Shortley phase factor. Once, not twice, with the associated Legendre polynomials takes no Condon–Shortley phase and an explicit form on the forehead. image

  2. The real form of spherical harmonics can be infered from the the complex form. image

  3. Point is, by obtaining the real form from the complex, the functionality of the (-1)^m used here. Is it the canceling or just follow the convention. Regardless of this point, there are two things we can confirm:

a. The complex form of spherical harmonics contains the Condon–Shortley phase as showed in wiki. b. The real form of spherical harmonics, after transformation, still contains a (-1)^m, we just don't know how to name it. see below: image

However, as many of the other online forums, some people take the (-1)^m (b) as Condon–Shortley phase for the real form, which leads to a simplified form as listed in another wiki. In this way, people conclude that

"the real form of spherical harmonics listed in the wiki follows the Condon–Shortley convention."

jinzx10 commented 1 month ago

Thanks for your feedback! I just have one more comment: while the complex spherical harmonics used by physicists almost always contain a Condon-Shortley phase, it still remains the question of whether to put this phase inside or outside the associated Legendre polynomial. In fact, I believe the $P_l^m$ in the equation you posted under Item 3.b already contain the Condon-Shortley phase (which is consistent with the wiki page of associated Legendre polynomial), so the hidden one will cancel the explicit one. This can be easily verified by taking $l=1, m=1$.

Franklalalala commented 1 month ago

Yes, me and another fellow has found this case 3 hours ago. The associated Legendre polynomial in 3.b indeed contains a phase, as can be verified by another website.

5215b1dd774ace134eb32229367913c

I think now we are very much clear about the whole progress, it has been a great pleasure to learn this new concept today. Thank you for your valuable discussion!