球谐函数与球谐函数导数乘积的积分求解咨询

嘿，这个问题在量子力学、计算电磁学或者角动量耦合的计算里经常碰到，我来一步步给你拆解清楚。

先铺垫关键公式：球谐函数对θ的导数展开 #

首先，球谐函数对极角$\theta$的偏导数可以表示为相邻阶数球谐函数的线性组合，这个标准公式是：
$$
\partial_\theta Y_l^m(\theta,\phi) = \sqrt{\frac{(l-m)(l+m+1)}{4(l+1)}} Y_{l+1}^m(\theta,\phi) - \sqrt{\frac{(l+m)(l-m+1)}{4l}} Y_{l-1}^m(\theta,\phi)
$$
（注：当$l=0$时，第二项的系数为0，因为分母有$l$，所以$\partial_\theta Y_0^0=0$，这也符合直观）

第一个积分：$\int_\Omega d\Omega, Y_{l}^m Y_{l^\prime}{m^\prime} \partial_\theta Y_{l^{{\prime\prime}}}{m^{\prime\prime}}$ #

把上面的导数展开式代入积分，我们可以把它拆成两个积分的和：
$$
\begin{align*}
&\int_\Omega d\Omega, Y_{l}^m Y_{l^\prime}{m^\prime} \partial_\theta Y_{l^{{\prime\prime}}}{m^{\prime\prime}} \
=& \sqrt{\frac{(l''-m'')(l''+m''+1)}{4(l''+1)}} \int_\Omega d\Omega, Y_l^m Y_{l'}^m' Y_{l''+1}^{m''} \
&- \sqrt{\frac{(l''+m'')(l''-m''+1)}{4l''}} \int_\Omega d\Omega, Y_l^m Y_{l'}^m' Y_{l''-1}^{m''}
\end{align*}
$$

接下来，我们利用**三个球谐函数的积分（Gaunt系数）**的性质：

• 这个积分只有当满足以下两个条件时才不为零：
  1. 阶数满足三角条件：$|l - l'| \leq l''' \leq l + l'$（这里$l'''$是第三个球谐函数的阶数，也就是$l''+1$或$l''-1$）
  2. 磁量子数满足：$m + m' + m''' = 0$（这里$m'''=m''$，所以要求$m + m' + m'' = 0$）
• 非零的情况下，积分值可以用Gaunt系数$G(l,l',l''';m,m',m''')$表示：
  $$
  \int_\Omega d\Omega, Y_l^m Y_{l'}^m' Y_{l'''}^m''' = G(l,l',l''';m,m',m''') \sqrt{\frac{(2l+1)(2l'+1)(2l'''+1)}{(4\pi)^2}}
  $$

所以总结一下，第一个积分的结果：

• 当$m + m' + m'' \neq 0$，或者阶数不满足三角条件（对$l''+1$或$l''-1$）时，积分值为0；
• 否则，就是上面两个系数分别乘以对应的Gaunt系数相关项的和。

第二个积分：$\int_\Omega d\Omega, Y_{l}^m \partial_\theta Y_{l^\prime}{m^\prime} \partial_\theta Y_{l^{{\prime\prime}}}{m^{\prime\prime}}$ #

同样的思路，先把两个导数都用球谐函数展开：
$$
\partial_\theta Y_{l'}^m' = \sqrt{\frac{(l'-m')(l'+m'+1)}{4(l'+1)}} Y_{l'+1}^{m'} - \sqrt{\frac{(l'+m')(l'-m'+1)}{4l'}} Y_{l'-1}^{m'}
$$
$$
\partial_\theta Y_{l''}^{m''} = \sqrt{\frac{(l''-m'')(l''+m''+1)}{4(l''+1)}} Y_{l''+1}^{m''} - \sqrt{\frac{(l''+m'')(l''-m''+1)}{4l''}} Y_{l''-1}^{m''}
$$

把这两个展开式相乘，会得到四个项的和，然后每个项和$Y_l^m$做曲面积分，就变成四个三个球谐函数的积分：
$$
\begin{align*}
&\int_\Omega d\Omega, Y_l^m \partial_\theta Y_{l'}^m' \partial_\theta Y_{l''}^{m''} \
=& C_1 \int_\Omega d\Omega, Y_l^m Y_{l'+1}^{m'} Y_{l''+1}^{m''} + C_2 \int_\Omega d\Omega, Y_l^m Y_{l'+1}^{m'} Y_{l''-1}^{m''} \
&+ C_3 \int_\Omega d\Omega, Y_l^m Y_{l'-1}^{m'} Y_{l''+1}^{m''} + C_4 \int_\Omega d\Omega, Y_l^m Y_{l'-1}^{m'} Y_{l''-1}^{m''}
\end{align*}
$$
其中$C_1,C_2,C_3,C_4$是四个乘积项的系数（比如$C_1$就是两个展开式中第一项系数的乘积，以此类推）。

每个积分项的非零条件和第一个积分一样：

• 磁量子数满足$m + m' + m'' = 0$（因为每个积分里的磁量子数是$m, m', m''$）
• 阶数满足三角条件：$|l - (l'\pm1)| \leq (l''\pm1) \leq l + (l'\pm1)$

非零的时候，每个积分项同样用Gaunt系数来计算，最后把四个非零的项加起来就是结果。

内容的提问来源于stack exchange，提问作者Fluid