广义球谐函数

                     

贡献者: addis

预备知识 CG 系数,球谐函数

   有了 CG 系数的相位约定和球谐函数的相位约定,就可以由球谐函数定义广义球谐函数(Generalized Spherical Harmonics)

\begin{equation} \mathcal{Y}_{l_1,l_2}^{L,M}( \hat{\boldsymbol{\mathbf{r}}} _1, \hat{\boldsymbol{\mathbf{r}}} _2) = \sum_{m_1, m_2} \begin{bmatrix}l_1 & l_2 & L\\ m_1 & m_2 & M\end{bmatrix} Y_{l_1 m_1}( \hat{\boldsymbol{\mathbf{r}}} _1) Y_{l_2 m_2} ( \hat{\boldsymbol{\mathbf{r}}} _2)~. \end{equation}
由于只有 $m_1 + m_2 = M$ 时 CG 系数才不为零,公式中令 $m_2 = M - m_1$,求和就由双重求和变为了单个求和,而且只有有限项。$m_1$ 的上下限满足 $-l_1 \leqslant m_1 \leqslant l_1$,$-l_2 \leqslant M-m_1 \leqslant l_2$ 即
\begin{equation} \max\{-l_1, M-l_2\} \leqslant m_1 \leqslant \max\{l_1, M+l_1\}~. \end{equation}
式 1 可以看作一个酉矩阵乘以列向量。

   由于式 1 是一个正交变换,其逆变换为

\begin{equation} Y_{l_1 m_1}( \hat{\boldsymbol{\mathbf{r}}} _1) Y_{l_2 m_2} ( \hat{\boldsymbol{\mathbf{r}}} _2) = \sum_{L} \begin{bmatrix}l_1 & l_2 & L\\ m_1 & m_2 & M\end{bmatrix} \mathcal{Y}_{l_1,l_2}^{L,M}( \hat{\boldsymbol{\mathbf{r}}} _1, \hat{\boldsymbol{\mathbf{r}}} _2)~. \end{equation}
由于球谐函数一般是复函数,而 CG 系数是实数,所以 $\mathcal{Y}_{l_1,l_2}^{L,M}$ 一般也是复函数
未完成:什么时候是实数?

Mathematica 代码

   在 Mathematica 中可以把广义球谐函数定义为1

YY[l1_, l2_, L_, M_, θ1_, ϕ1_, θ2_, ϕ2_] := 
 Sum[ClebschGordan[{l1, m1}, {l2, M - m1}, {L, M}] SphericalHarmonicY[
    l1, m1, θ1, ϕ1] SphericalHarmonicY[l2, 
    M - m1, θ2, ϕ2], {m1, Max[-l1, -l2 + M], 
   Min[l2, l2 + M]}]

1. 微分方程

   广义球谐函数是量子力学中双粒子总角动量算符 $ \hat{L} ^2, \hat{L} _z$ 的本征函数。

\begin{equation} \hat{L} ^2 \mathcal Y_{l_1,l_2}^{L,M} = L(L+1) \hbar^2 \mathcal Y_{l_1,l_2}^{L,M}~, \qquad \hat{L} _z \mathcal Y_{l_1,l_2}^{L,M} = M\hbar \mathcal Y_{l_1,l_2}^{L,M}~. \end{equation}

2. 宇称

   宇称(parity)算符 $\Pi$ 作用在一个函数上,就是把函数所有自变量乘以 $-1$。其本征函数是所有中心对称或反对称的函数,本征值为分别为 $\pm 1$。

   球谐函数是宇称算符的本征矢,本征值为 $(-1)^l$(式 25 ),易得广义球谐函数也是宇称算符的本征矢,本征值为 $(-1)^{l_1+l_2}$

\begin{equation} \Pi \mathcal{Y}_{l_1,l_2}^{L,M}( \hat{\boldsymbol{\mathbf{r}}} _1, \hat{\boldsymbol{\mathbf{r}}} _2) = \mathcal{Y}_{l_1,l_2}^{L,M}(- \hat{\boldsymbol{\mathbf{r}}} _1, - \hat{\boldsymbol{\mathbf{r}}} _2) = (-1)^{l_1+l_2} \mathcal{Y}_{l_1,l_2}^{L,M}( \hat{\boldsymbol{\mathbf{r}}} _1, \hat{\boldsymbol{\mathbf{r}}} _2)~. \end{equation}

3. 交换对称性

   广义球谐函数不具有交换对称性2(除非 $l_1 = l_2$)

\begin{equation} P_{12} \mathcal{Y}_{l_1,l_2}^{L,M}( \hat{\boldsymbol{\mathbf{r}}} _1, \hat{\boldsymbol{\mathbf{r}}} _2) = \mathcal{Y}_{l_1,l_2}^{L,M}( \hat{\boldsymbol{\mathbf{r}}} _2, \hat{\boldsymbol{\mathbf{r}}} _1) = (-1)^{l_1+l_2-L} \mathcal{Y}_{l_2,l_1}^{L,M}( \hat{\boldsymbol{\mathbf{r}}} _1, \hat{\boldsymbol{\mathbf{r}}} _2)~, \end{equation}
也可以记为
\begin{equation} \mathcal{Y}_{l_1,l_2}^{L,M}( \hat{\boldsymbol{\mathbf{r}}} _1, \hat{\boldsymbol{\mathbf{r}}} _2) = (-1)^{l_1+l_2-L} \mathcal{Y}_{l_2,l_1}^{L,M}( \hat{\boldsymbol{\mathbf{r}}} _2, \hat{\boldsymbol{\mathbf{r}}} _1)~. \end{equation}
但是我们可以轻易构造对称(+)或反对称(-)的函数
\begin{equation} \mathcal{Y}_{l_1,l_2}^{L,M}( \hat{\boldsymbol{\mathbf{r}}} _1, \hat{\boldsymbol{\mathbf{r}}} _2) \pm (-1)^{l_1+l_2-L}\mathcal{Y}_{l_2,l_1}^{L,M}( \hat{\boldsymbol{\mathbf{r}}} _1, \hat{\boldsymbol{\mathbf{r}}} _2)~, \end{equation}
\begin{equation} \Psi( \boldsymbol{\mathbf{r}} _1, \boldsymbol{\mathbf{r}} _2) = R(r_1, r_2)\mathcal{Y}_{l_1,l_2}^{L,M}( \hat{\boldsymbol{\mathbf{r}}} _1, \hat{\boldsymbol{\mathbf{r}}} _2) - R(r_1, r_2)\mathcal{Y}_{l_2,l_1}^{L,M}( \hat{\boldsymbol{\mathbf{r}}} _1, \hat{\boldsymbol{\mathbf{r}}} _2)~, \qquad (l_1+l_2-L = odd). \end{equation}

   证明可以用 CG 系数的对称性(交换左边两列,式 7

\begin{equation} \begin{aligned} \mathcal{Y}_{l_1,l_2}^{L,M}( \hat{\boldsymbol{\mathbf{r}}} _2, \hat{\boldsymbol{\mathbf{r}}} _1) &= \sum_{m_1, m_2} \begin{bmatrix}l_1 & l_2 & L\\ m_1 & m_2 & M\end{bmatrix} Y_{l_1, m_1}( \hat{\boldsymbol{\mathbf{r}}} _2)Y_{l_2, m_2}( \hat{\boldsymbol{\mathbf{r}}} _1) \\ &= (-1)^{l_1+l_2-L} \sum_{m_1, m_2} \begin{bmatrix}l_2 & l_1 & L\\ m_2 & m_1 & M\end{bmatrix} Y_{l_2, m_2}( \hat{\boldsymbol{\mathbf{r}}} _1) Y_{l_1, m_1}( \hat{\boldsymbol{\mathbf{r}}} _2)\\ & = (-1)^{l_1+l_2-L} \mathcal{Y}_{l_2,l_1}^{L,M}( \hat{\boldsymbol{\mathbf{r}}} _1, \hat{\boldsymbol{\mathbf{r}}} _2)~, \end{aligned} \end{equation}
式 6 也可以得到对易关系
\begin{equation} \left[L^2, P_{12}\right] = \left[L_z, P_{12}\right] = 0~, \end{equation}
即交换两个粒子不改变总角动量。

   操作类算符(宇称,平移,交换)如果与某物理量算符对易,就说明波函数经过该操作,改物理量守恒。如果哈密顿中的某一项算符(例如 $z$ 方向的电场的 dipole)如果与某物理量算符(例如 $L_z$)对易,就说明波函数经过该传播子传播,该物理量守恒。

4. 正交性

   将式 10 用这里的符号表示,就是

\begin{equation} \int \mathcal{Y}_{l_1', l_2'}^{L', M'}( \hat{\boldsymbol{\mathbf{r}}} _1, \hat{\boldsymbol{\mathbf{r}}} _2) \mathcal{Y}_{l_1, l_2}^{L, M}( \hat{\boldsymbol{\mathbf{r}}} _1, \hat{\boldsymbol{\mathbf{r}}} _2) \,\mathrm{d}{\Omega_1} \,\mathrm{d}{\Omega_2} = \delta_{l_1, l_1'} \delta_{l_2, l_2'} \delta_{L, L'} \delta_{M, M'}~. \end{equation}

5. 其他性质

\begin{equation} \mathcal{Y}_{l_1,l_2}^{L,-M}( \hat{\boldsymbol{\mathbf{r}}} _1, \hat{\boldsymbol{\mathbf{r}}} _2) = (-1)^{l_1 + l_2 + L + M} \mathcal{Y}_{l_1,l_2}^{L,M}( \hat{\boldsymbol{\mathbf{r}}} _1, \hat{\boldsymbol{\mathbf{r}}} _2)^*~, \end{equation}
其中 $*$ 表示复共轭。推导如下,
\begin{equation} \begin{aligned} \mathcal{Y}_{l_1,l_2}^{L,-M}( \hat{\boldsymbol{\mathbf{r}}} _1, \hat{\boldsymbol{\mathbf{r}}} _2) &= \sum_{m_1+m_2 = M} \begin{bmatrix}l_1 & l_2 & L \\ -m_1 & -m_2 & -M\end{bmatrix} Y_{l_1,-m_1}( \hat{\boldsymbol{\mathbf{r}}} _1) Y_{l_2, -m_2} ( \hat{\boldsymbol{\mathbf{r}}} _2)\\ &= (-1)^{l_1+l_2+L} \sum_{m_1+m_2 = M} \begin{bmatrix}l_1 & l_2 & L \\ m_1 & m_2 & M\end{bmatrix} Y_{l_1,-m_1}( \hat{\boldsymbol{\mathbf{r}}} _1) Y_{l_2, -m_2} ( \hat{\boldsymbol{\mathbf{r}}} _2)\\ &= (-1)^{l_1+l_2+L+M} \sum_{m_1+m_2 = M} \begin{bmatrix}l_1 & l_2 & L \\ m_1 & m_2 & M\end{bmatrix} Y_{l_1,m_1}^*( \hat{\boldsymbol{\mathbf{r}}} _1) Y_{l_2, m_2}^* ( \hat{\boldsymbol{\mathbf{r}}} _2)\\ &= (-1)^{l_1+l_2+L+M} \mathcal{Y}_{l_1,l_2}^{L,M}( \hat{\boldsymbol{\mathbf{r}}} _1, \hat{\boldsymbol{\mathbf{r}}} _2)^*~. \end{aligned} \end{equation}


1. ^ 使用时出现的 “not triangular” 警告可忽略。
2. ^ 式中 $l_1+l_2-L$ 前面的负号可以变为加号,但习惯上仍然写成负号,因为虽然球谐函数要求 $L$ 为整数,但在涉及到电子自旋时 $L$ 有可能是半整数,这时就不能变为加号了。


致读者: 小时百科一直以来坚持所有内容免费无广告,这导致我们处于严重的亏损状态。 长此以往很可能会最终导致我们不得不选择大量广告以及内容付费等。 因此,我们请求广大读者热心打赏 ,使网站得以健康发展。 如果看到这条信息的每位读者能慷慨打赏 20 元,我们一周就能脱离亏损, 并在接下来的一年里向所有读者继续免费提供优质内容。 但遗憾的是只有不到 1% 的读者愿意捐款, 他们的付出帮助了 99% 的读者免费获取知识, 我们在此表示感谢。

                     

友情链接: 超理论坛 | ©小时科技 保留一切权利