球多极矩

源位置为 ${\displaystyle \mathbf {r} ^{\prime }}$  的点电荷 ${\displaystyle q}$  ，其电势 ${\displaystyle \Phi (\mathbf {r} )}$  在场位置 ${\displaystyle \mathbf {r} }$  为

${\displaystyle \Phi (\mathbf {r} )={\frac {q}{4\pi \varepsilon _{0}|\mathbf {r} -\mathbf {r^{\prime }} |}}={\frac {q}{4\pi \varepsilon _{0}}}{\frac {1}{\sqrt {r^{2}+r^{\prime 2}-2r^{\prime }r\cos \gamma }}}}$  ；

其中，${\displaystyle \varepsilon _{0}}$  是电常数，${\displaystyle \gamma }$  是 ${\displaystyle \mathbf {r} }$  与 ${\displaystyle \mathbf {r} ^{\prime }}$  之间的夹角。

假设 ${\displaystyle r'<r}$  ，场位置比源位置离原点更远，则此距离倒数函数 ${\displaystyle 1/|\mathbf {r} -\mathbf {r^{\prime }} |}$  以 ${\displaystyle r^{\prime }/r}$  的幂和勒壤得多项式展开为^[1] ：

${\displaystyle \Phi (\mathbf {r} )={\frac {q}{4\pi \varepsilon _{0}r}}\sum _{\ell =0}^{\infty }\left({\frac {r^{\prime }}{r}}\right)^{\ell }P_{\ell }(\cos \gamma )}$  。

应用球余弦定律（spherical law of cosine）， ${\displaystyle \cos \gamma }$  表示为

${\displaystyle \cos \gamma =\cos \theta \cos \theta ^{\prime }+\sin \theta \sin \theta ^{\prime }\cos(\phi -\phi ^{\prime })}$  。

这结果也可以直接用向量代数直接计算出来。

应用球谐函数加法定理，${\displaystyle P_{\ell }(\cos \gamma )}$  又表示为^[2]

${\displaystyle P_{\ell }(\cos \gamma )={\frac {4\pi }{2\ell +1}}\sum _{m=-\ell }^{\ell }Y_{\ell m}(\theta ,\phi )Y_{\ell m}^{*}(\theta ^{\prime },\phi ^{\prime })}$  ；

其中，${\displaystyle Y_{\ell m}}$  是球谐函数。

将这方程式代入电势的方程式，可以得到

${\displaystyle \Phi (\mathbf {r} )={\frac {q}{4\pi \varepsilon _{0}r}}\sum _{\ell =0}^{\infty }\left({\frac {r^{\prime }}{r}}\right)^{\ell }\left({\frac {4\pi }{2\ell +1}}\right)\sum _{m=-\ell }^{\ell }Y_{\ell m}(\theta ,\phi )Y_{\ell m}^{*}(\theta ^{\prime },\phi ^{\prime })}$  。

点电荷的“球多极矩” 定义为

${\displaystyle q_{\ell m}\ {\stackrel {\mathrm {def} }{=}}\ qr^{\prime \ell }Y_{\ell m}^{*}(\theta ^{\prime },\phi ^{\prime })}$  。

则电势的方程式又可写为

${\displaystyle \Phi (\mathbf {r} )={\frac {1}{\varepsilon _{0}}}\sum _{\ell =0}^{\infty }\sum _{m=-\ell }^{\ell }{\frac {q_{\ell m}Y_{\ell m}(\theta ,\phi )}{(2\ell +1)r^{\ell +1}}}}$  。

假设 ${\displaystyle r<r'}$  ，场位置比源位置离原点更近，则此距离倒数函数 ${\displaystyle 1/|\mathbf {r} -\mathbf {r^{\prime }} |}$  可以以 ${\displaystyle r/r^{\prime }}$  的幂和勒壤得多项式展开：

${\displaystyle \Phi (\mathbf {r} )={\frac {q}{4\pi \varepsilon _{0}r^{\prime }}}\sum _{\ell =0}^{\infty }\left({\frac {r}{r^{\prime }}}\right)^{\ell }\left({\frac {4\pi }{2\ell +1}}\right)\sum _{m=-\ell }^{\ell }Y_{\ell m}(\theta ,\phi )Y_{\ell m}^{*}(\theta ^{\prime },\phi ^{\prime })}$  。

点电荷的“内部球多极矩”（前述的球多极矩称为外部球多极矩）定义为

${\displaystyle I_{\ell m}\ {\stackrel {\mathrm {def} }{=}}\ {\frac {q}{\left(r^{\prime }\right)^{\ell +1}}}Y_{\ell m}^{*}(\theta ^{\prime },\phi ^{\prime })}$  。

则电势的方程式写为

${\displaystyle \Phi (\mathbf {r} )={\frac {1}{\varepsilon _{0}}}\sum _{\ell =0}^{\infty }\sum _{m=-\ell }^{\ell }{\frac {I_{\ell m}r^{\ell }Y_{\ell m}(\theta ,\phi )}{2\ell +1}}}$  。

前述多极展开方法可以推广至电荷密度分布。将点电荷 ${\displaystyle q}$  改换为微小电荷元素 ${\displaystyle \rho (\mathbf {r} ^{\prime })d\mathbf {r} ^{\prime }}$  ，然后积分，则可得到电势的方程式（假设 ${\displaystyle r'<r}$  ）：

${\displaystyle \Phi (\mathbf {r} )={\frac {1}{\varepsilon _{0}}}\sum _{\ell =0}^{\infty }\sum _{m=-\ell }^{\ell }{\frac {q_{\ell m}Y_{\ell m}(\theta ,\phi )}{(2\ell +1)r^{\ell +1}}}}$  ；

其中，电荷密度分布的球多极矩定义为 ${\displaystyle q_{\ell m}\ {\stackrel {\mathrm {def} }{=}}\ \int _{\mathbb {V} '}\rho (\mathbf {r} ^{\prime })\left(r^{\prime }\right)^{\ell }Y_{\ell m}^{*}(\theta ^{\prime },\phi ^{\prime })\ \mathrm {d} ^{3}\mathbf {r} ^{\prime }}$  ，${\displaystyle \mathbb {V} '}$  是积分体积。

特别注意，由于电势 ${\displaystyle \Phi (\mathbf {r} )}$  为实值，这展开式的复共轭也是同样正确的球多极展开式。然而，这样做会导致球多极矩的定义式含有 ${\displaystyle Y_{\ell m}}$  项目，而不是其复共轭数 ${\displaystyle Y_{\ell m}^{*}}$  。在某些领域，例如物理化学，这是一般常规。更详尽资料，请参阅条目分子多极矩（molecular multipole moment）。

内部球多极矩

类似地，假设 ${\displaystyle r<r'}$  ，场位置比源位置离原点更近，则电势的方程式为

${\displaystyle \Phi (\mathbf {r} )={\frac {1}{\varepsilon _{0}}}\sum _{\ell =0}^{\infty }\sum _{m=-\ell }^{\ell }{\frac {I_{\ell m}r^{\ell }Y_{\ell m}(\theta ,\phi )}{2\ell +1}}}$  ；

其中，电荷密度分布的内部球多极矩定义为 ${\displaystyle I_{\ell m}\ {\stackrel {\mathrm {def} }{=}}\ \int _{\mathbb {V} '}{\frac {\rho (\mathbf {r} ^{\prime })Y_{\ell m}^{*}(\theta ^{\prime },\phi ^{\prime })}{\left(r^{\prime }\right)^{\ell +1}}}\ \mathrm {d} ^{3}\mathbf {r} ^{\prime }}$  。

两个互不重叠，同心的电荷分布可以用简单公式来描述。设定第一个电荷分布 ${\displaystyle \rho _{1}}$  在第二个电荷分布 ${\displaystyle \rho _{2}}$  的内部，则由 ${\displaystyle \rho _{1}}$  所产生的电势 ${\displaystyle \Phi _{1}}$  ，因为作用于 ${\displaystyle \rho _{2}}$  而涉及的相互作用能 ${\displaystyle U}$  为

${\displaystyle U=\int _{\mathbb {V} }\rho _{2}(\mathbf {r} )\Phi _{1}(\mathbf {r} )\ \mathrm {d} ^{3}\mathbf {r} }$  。

电势 ${\displaystyle \Phi _{1}(\mathbf {r} )}$  可以以外部球多极矩展开为

${\displaystyle \Phi _{1}(\mathbf {r} )={\frac {1}{\varepsilon _{0}}}\sum _{\ell =0}^{\infty }\sum _{m=-\ell }^{\ell }{\frac {q_{1\ell m}Y_{\ell m}(\theta ,\phi )}{(2\ell +1)r^{\ell +1}}}}$  ；

其中，${\displaystyle q_{1\ell m}}$  是第一个电荷分布的 ${\displaystyle \ell m}$  外部球多极矩。

将这方程式代入相互作用能 ${\displaystyle U}$  的方程式，可以得到

${\displaystyle U={\frac {1}{\varepsilon _{0}}}\sum _{\ell =0}^{\infty }\sum _{m=-\ell }^{\ell }{\frac {q_{1\ell m}}{2\ell +1}}\int _{\mathbb {V} }{\frac {\rho _{2}(\mathbf {r} )Y_{\ell m}(\theta ,\phi )}{r^{\ell +1}}}\ \mathrm {d} ^{3}\mathbf {r} }$  。

注意到其积分项目等于 ${\displaystyle \rho _{2}(\mathbf {r} ^{\prime })}$  的内部球多极矩 ${\displaystyle I_{2\ell m}}$  的复共轭数，相互作用能 ${\displaystyle U}$  的方程式约化为简单形式

${\displaystyle U={\frac {1}{\varepsilon _{0}}}\sum _{\ell =0}^{\infty }\sum _{m=-\ell }^{\ell }{\frac {q_{1\ell m}I_{2\ell m}^{*}}{2\ell +1}}}$  。

这方程式可以用来计算，原子核产生的电势因为与其周围的原子轨域耦合而涉及的相互作用能。反过来，给定相互作用能与电子轨域的内部球多极矩，则可以计算原子核的外部球多极矩，从而得知其形状。

假设电荷密度为“轴对称”，即与方位角 ${\displaystyle \phi ^{\prime }}$  无关，则球多极展开式的形式很简单。在 ${\displaystyle q_{\ell m}}$  与 ${\displaystyle I_{\ell m}}$  的定义式内，对于 ${\displaystyle \phi ^{\prime }}$  积分，则可以发觉除了 ${\displaystyle m=0}$  球多极矩以外，其它球多极矩都等于零。应用数学恒等式^[2]

${\displaystyle P_{\ell }(\cos \theta )={\sqrt {\frac {4\pi }{2\ell +1}}}Y_{\ell 0}(\theta ,\phi )}$  ，

轴对称球多极矩定义为

${\displaystyle q_{\ell }\ {\stackrel {\mathrm {def} }{=}}\ \int _{\mathbb {V} '}{\sqrt {\frac {2\ell +1}{4\pi }}}\ \rho (\mathbf {r} ^{\prime })\left(r^{\prime }\right)^{\ell }P_{\ell }(\cos \theta ^{\prime })\ \mathrm {d} ^{3}\mathbf {r} ^{\prime }}$  ，

则外部球多极展开式为

${\displaystyle \Phi (\mathbf {r} )={\frac {1}{4\pi \varepsilon _{0}}}\sum _{\ell =0}^{\infty }{\sqrt {\frac {4\pi }{2\ell +1}}}\ {\frac {q_{\ell }P_{\ell }(\cos \theta )}{r^{\ell +1}}}}$  。

类似地，轴对称内部球多极矩定义为

${\displaystyle I_{\ell }\ {\stackrel {\mathrm {def} }{=}}\ \int _{\mathbb {V} '}{\sqrt {\frac {2\ell +1}{4\pi }}}\ {\frac {\rho (\mathbf {r} ^{\prime })}{\left(r^{\prime }\right)^{\ell +1}}}P_{\ell }(\cos \theta ^{\prime })\ \mathrm {d} ^{3}\mathbf {r} ^{\prime }}$  ，

内部球多极展开式为

${\displaystyle \Phi (\mathbf {r} )={\frac {1}{4\pi \varepsilon _{0}}}\sum _{\ell =0}^{\infty }{\sqrt {\frac {4\pi }{2\ell +1}}}I_{\ell }r^{\ell }P_{\ell }(\cos \theta )}$  。

注意到 ${\displaystyle q_{\ell (-m)}=(-1)^{m}q_{\ell m}^{*}}$  。以下列出几个最低阶的球多极矩的表达式，以及与笛卡儿多极矩之间的关系^[2]：

${\displaystyle {\begin{aligned}q_{00}&={\frac {1}{\sqrt {4\pi }}}\int _{\mathbb {V'} }\rho (\mathbf {r} ')\ \mathrm {d} ^{3}\mathbf {r} '&&={\frac {1}{\sqrt {4\pi }}}\ q\\q_{11}&=-{\sqrt {\frac {3}{8\pi }}}\int _{\mathbb {V'} }r'\sin {\theta '}\ e^{-i\phi '}\rho (\mathbf {r} ')\ \mathrm {d} ^{3}\mathbf {r} '&&=-{\sqrt {\frac {3}{8\pi }}}\ (p_{x}-ip_{y})\\q_{10}&={\sqrt {\frac {3}{4\pi }}}\int _{\mathbb {V'} }r'\cos {\theta '}\ \rho (\mathbf {r} ')\ \mathrm {d} ^{3}\mathbf {r} '&&=-{\sqrt {\frac {3}{4\pi }}}\ p_{z}\\q_{22}&={\sqrt {\frac {15}{32\pi }}}\int _{\mathbb {V'} }r^{\prime 2}\sin ^{2}{\theta '}\ e^{-2i\phi '}\rho (\mathbf {r} ')\ \mathrm {d} ^{3}\mathbf {r} '&&={\sqrt {\frac {15}{288\pi }}}\ (Q_{11}-2iQ_{12}-Q_{22})\\q_{21}&=-{\sqrt {\frac {15}{8\pi }}}\int _{\mathbb {V'} }r^{\prime 2}\sin {\theta '}\cos {\theta '}\ e^{-i\phi '}\rho (\mathbf {r} ')\ \mathrm {d} ^{3}\mathbf {r} '&&=-{\sqrt {\frac {15}{72\pi }}}\ (Q_{13}-iQ_{33})\\q_{20}&={\sqrt {\frac {5}{16\pi }}}\int _{\mathbb {V'} }r^{\prime 2}(\cos ^{2}{\theta '}-1)\rho (\mathbf {r} ')\ \mathrm {d} ^{3}\mathbf {r} '&&={\sqrt {\frac {5}{16\pi }}}\ Q_{33}\end{aligned}}}$  。

其中，${\displaystyle (p_{x},p_{y},p_{z})}$  是笛卡儿电偶极矩，${\displaystyle Q_{ij}}$  是笛卡儿电四极矩（electric quadruple moment）。

• 立体调和函数（solid harmonics）
• 圆柱多极矩