球多極矩

給予在源位置 ${\displaystyle \mathbf {r} '=(r',\theta ',\phi ')}$ 的電荷分佈，計算在場位置 ${\displaystyle \mathbf {r} =(r,\theta ,\phi )}$ 產生的電勢。

源位置為 ${\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 )}$ 。

應用球餘弦定律（）， ${\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}^{*}}$ 。在某些領域，例如物理化學，這是一般常規。更詳盡資料，請參閱條目分子多極矩（）。

兩個互不重疊，同心的電荷分佈可以用簡單公式來描述。設定第一個電荷分佈 ${\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}}$ 是笛卡兒電四極矩（）。

• 立體調和函數（）
• 圓柱多極矩