广义球谐函数

                     

贡献者: 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$ 有可能是半整数,这时就不能变为加号了。

                     

© 小时科技 保留一切权利