球面調和関数

球面調和関数（きゅうめんちょうわかんすう、英: spherical harmonics^[1]）あるいは球関数（きゅうかんすう、英: spherical functions^[2]）は以下のいずれかを意味する関数である：

1. n 次元ラプラス方程式の解となる斉次多項式を単位球面に制限する事で得られる関数。
2. 次元 n が 3 の場合の 1 の意味での球面調和関数で、球面座標 (r, θ, φ) で書いたラプラス方程式の変数分離解を記述するのに用いる事ができる関数 Y n
   k (θ, φ).

本項では 1 及び 2 双方の意味の球面調和関数について述べるが、特に断りがない限り、「球面調和関数」という言葉を 1 の意味で用いる。

R を実数全体の集合とし、C を複素数全体の集合とし、n 個の実数からなる組の集合を Rⁿ とし、Rⁿ の元を (x₁, …, x_n) ∈ Rⁿ と書き表すことにする。

Rⁿ 上の複素数値関数

φ: Rⁿ → C

が2階微分可能なとき、Δφ を

${\displaystyle \Delta \phi =\sum _{i=1}^{n}{\frac {\partial ^{2}}{\partial x_{i}{}^{2}}}\phi }$

と定義し、Δ をラプラス作用素という。さらに Rⁿ 上の（複素数値の）多項式 p(x₁, …, x_n) で

Δp = 0

を満たすものを調和多項式（英語版）という^[3]。なおラプラス作用素は回転行列 R に対し、

Δp(R(x)) = R(Δp(x))

を満たすので^[4]、調和多項式の定義は座標系のとり方に依存しない。

調和多項式 p が k 次の斉次多項式であるとき、p を単位球面

${\displaystyle S^{n-1}=\{(x_{1},\ldots ,x_{n})\in \mathbf {R} ^{n}\mid \sum _{i}x_{i}{}^{2}=1\}}$

(P1)

に制限した制限写像

${\displaystyle Y=p|_{S^{n-1}}\colon S^{n-1}\to \mathbf {C} }$

を k 次の球面調和関数という^[5]。

（n 次元空間 Rⁿ における）k 次の球面調和関数全体の集合を H_k とすると、H_k は C 上のベクトル空間であり、

${\displaystyle \operatorname {dim} _{\mathbf {C} }{\mathcal {H}}_{k}={\frac {(n+2k-2)(n+k-3)!}{(n-2)!k!}}}$

(P2)

である^[6]。

e_n を Rⁿ 上のベクトル

e_n = (0, ..., 0, 1) ∈ Rⁿ

とする。

次の事実が成立する^[7]。

帯球関数を具体的に書き表す為、記号を導入する。自然数 n と非負の実数 x に対しポッホハマー記号 (x)_n を

${\displaystyle (x)_{n}:={\frac {\Gamma (x+n)}{\Gamma (x)}}}$

により定義する。ここで Γ(x) はガンマ関数である。さらにガウスの超幾何関数を

${\displaystyle {}_{2}F_{1}(-k,b;c;z)=\sum _{i=0}^{k}(-1)^{i}{\binom {k}{i}}{\frac {(b)_{i}}{(c)_{i}}}z^{i}}$

により定義し^{[注 1]}、さらに超球多項式を

${\displaystyle P_{k}^{(\alpha )}(z):=\,_{2}F_{1}\!\left(-k,2\alpha +k;\alpha +{\frac {1}{2}};{\frac {1-z}{2}}\right)}$

により定義する^{[注 2]}。このとき、次が成立する。

${\displaystyle Z_{k}^{n}(x_{1},\ldots ,x_{n}):=P_{k}^{((n-2)/2)}(x_{n})}$ は k 次の帯球関数である^[8]。

すでに述べたように、k 次の帯球関数は定数倍を除いて一意なので、全ての k 次帯球関数は上述したものの定数倍として表記可能である。

3次元空間 R³ の場合、R³ を球面座標 (r, θ, φ) で表すと、下記の Y m
k (θ, φ) が球面調和関数になる事が知られている。

${\displaystyle Y_{k}^{m}(\theta ,\phi )=(-1)^{(m+|m|)/2}{\sqrt {{\frac {2k+1}{4\pi }}{\frac {(k-|m|)!}{(k+|m|)!}}\,}}\,P_{k}^{|m|}(\cos \theta )\,e^{im\phi }.}$

(B1)

ここで

であり、P m
k (t) はルジャンドルの陪多項式^[9]

${\displaystyle P_{k}{}^{m}(t)={\frac {1}{2^{k}}}(1-t^{2})^{\tfrac {m}{2}}\sum _{j=0}^{\lfloor (k-m)/2\rfloor }{(-1)^{j}(2k-2j)! \over j!(k-j)!(k-2j-m)!}t^{k-2j-m}}$

(B3)

である。すなわち P m
k (t) はルジャンドルの陪微分方程式

${\displaystyle (1-t^{2})y''(t)-2t\,y'(t)+\left[k(k+1)-{\frac {m^{2}}{1-t^{2}}}\right]y(t)=0}$

の解である。なお、ルジャンドルの陪微分方程式は条件 (B2) を満たすとき、およびそのときだけ解を持つことが知られている。また、Y m
k (θ, φ) の定義における係数は、後述するノルムが1になるよう選んだものである。

Y m
k (θ, φ) が球面調和関数の定義を満たすことは自明ではないが、p を p(r, θ, φ) = r^k Y m
k (θ, φ) と定義した上で直交座標に変換すると p が斉次多項式になっている事を確認できる。

なお、本項では、「球面調和関数」という言葉をラプラス方程式の解となる斉次多項式（の球面への制限）一般を指す用語として用いるが、物理の教科書などでは上述した Y m
k (θ, φ) のみを球面調和関数と呼んでいるものも多い。

Y m
k (θ, φ) は斉次多項式に関する3次元空間のラプラス方程式を変数分離で解く事で自然に得られる。k 次の斉次多項式 p に対し、変数分離形

p(r, θ, φ) = R(r) Θ(θ) Φ(φ)

でラプラス方程式 Δp = 0 を解くと、変数分離形の解は必ず

${\displaystyle p(r,\theta ,\phi )=r^{k}(AY_{k}{}^{m}(\theta ,\phi )+BY_{k}{}^{-m}(\theta ,\phi )),}$ m は整数で k ≤ |m|

と書ける事を証明できる。

証明

Y m
k (θ, φ) は斉次多項式に関する3次元空間のラプラス方程式を変数分離で解く事で自然に得られるものである。このことを見るために3次元空間 R³ を球面座標 (r, θ, φ) でラプラス作用素を表記すると、

${\displaystyle \Delta ={\frac {1}{r^{2}}}(\Delta _{r}+\Delta _{s})}$

となる。ここで、

${\displaystyle \Delta _{r}={\frac {\partial }{\partial r}}\left(r^{2}{\frac {\partial }{\partial r}}\right),}$

${\displaystyle \Delta _{s}={\frac {1}{\sin \theta }}{\frac {\partial }{\partial \theta }}\sin \theta {\frac {\partial }{\partial \theta }}+{\frac {1}{\sin ^{2}\theta }}{\frac {\partial ^{2}}{\partial \varphi ^{2}}}}$

である。定義より次数 k の球面調和関数は、k 次の斉次多項式 p を単位球面上に制限したものとして表現可能である。p が k 次の斉次多項式である事から、p の極座標表示は

p(r, θ, φ) = r^kY(θ, φ)

(A1)

の形に書ける。ラプラス方程式 Δp = 0 の以下、変数分離解

Y(θ, φ) = Θ(θ) Φ(φ)

を求める。R(r) = r^k とすれば、

${\displaystyle {\frac {1}{R}}\Delta _{r}R={\frac {1}{R}}{\frac {\mathrm {d} }{\mathrm {d} r}}\left(r^{2}{\frac {\mathrm {d} R}{\mathrm {d} r}}\right)=k(k+1)}$

なので、変数分離解をラプラス方程式の極座標表示に代入することで、

${\displaystyle {\frac {1}{Y}}\Delta _{s}Y=-k(k+1)}$

が成立する。上式に対してさらに変数分離を適用する事で、複素数 m を適切に選べば

${\displaystyle {\frac {1}{\Phi }}{\frac {\mathrm {d} ^{2}\Phi }{\mathrm {d} \varphi ^{2}}}=-m^{2}}$

(A2)

${\displaystyle k(k+1)\sin ^{2}\theta +{\frac {\sin \theta }{\Theta }}{\frac {\mathrm {d} }{\mathrm {d} \theta }}\left(\sin \theta {\frac {\mathrm {d} \Theta }{\mathrm {d} \theta }}\right)=m^{2}}$

(A3)

が成立する事がわかる。以下、m が定数である場合の解を求める。

(A2) は初等的に解くことができ、一般解

${\displaystyle \Phi (\phi )=A\mathrm {e} ^{im\phi }+B\mathrm {e} ^{-im\phi }}$

(A4)

を得られる。ここで i は虚数単位である。それに対しスツルム＝リウヴィル型の微分方程式 (A3) は t = cos θ と変数変換すると、y(t) = Θ(arccos t) はルジャンドルの陪微分方程式

${\displaystyle (1-t^{2})y''(t)-2t\,y'(t)+\left(k(k+1)-{\frac {m^{2}}{1-t^{2}}}\right)y(t)=0}$

を満たす。よってルジャンドルの陪多項式 P m
k (t) を (B3) のように定義すれば、結論として

${\displaystyle \Theta (\theta )=P_{k}^{|m|}(\cos {\theta })}$

(A5)

がわかる。ここで k は (B2) の条件を満たす整数である。そこで Y m
k (θ, φ) を

${\displaystyle Y_{k}^{m}(\theta ,\phi )=(-1)^{(m+|m|)/2}{\sqrt {{\frac {2k+1}{4\pi }}{\frac {(k-|m|)!}{(k+|m|)!}}\,}}\,P_{k}^{|m|}(\cos \theta )\,e^{im\phi }}$

と定義すれば、(A1), (A4), (A5), (B2) より、変数分離形の k 次の調和多項式 p は必ず

${\displaystyle p(r,\theta ,\phi )=r^{k}(AY_{k}{}^{m}(\theta ,\phi )+BY_{k}{}^{-m}(\theta ,\phi )),}$ m は整数で k ≤ |m|,

と書ける事になる。なお、p を直交座標に変換すると p が斉次多項式になっている事を確認できる。

また、3次元空間の場合、k 次球面調和関数全体のなすベクトル空間 H_k の次元は、(P2) より

${\displaystyle \operatorname {dim} _{\mathbf {C} }{\mathcal {H}}_{k}=2k+1}$

なので、(B2) より、以下の結論が得られる：

本節では、球面調和関数の空間に内積を定義し、球面調和関数がこの内積に関して完全直交性を満たすことを示す。

n 次元空間 Rⁿ の単位球面 S^{n − 1} を (P1) のように定義し、dS を Sⁿ⁻¹ 上の面素とし、S^{n − 1} 上定義された2つの球面調和関数 f, g の内積を

${\displaystyle \langle f\mid g\rangle _{S^{n-1}}:=\int _{S^{n-1}}f(\mathbf {x} )g(\mathbf {x} )\,\mathrm {d} S}$

(C1)

により定義する。なお、面素 dS は球面座標 (r, θ₁, …, θ_{n − 1}) を

${\displaystyle x_{j}={\begin{cases}r\sin \theta _{1}\cdots \sin \theta _{j-1}\cos \theta _{j}&{\text{if }}1\leq j\leq n-1\\r\sin \theta _{1}\cdots \sin \theta _{n}&{\text{if }}j=n\end{cases}}}$

を用いて

${\displaystyle \mathrm {d} S=\prod _{j=1}^{n-1}\sin ^{j-1}\theta _{j}\,\mathrm {d} \theta _{1}\cdots \mathrm {d} \theta _{n-1}}$

と書ける^[10]。特に 3 次元空間の場合は球面座標 (r, θ, φ) に対し、

${\displaystyle \mathrm {d} S=\sin \theta \,\mathrm {d} \theta \,\mathrm {d} \varphi }$

である。

k 次球面調和関数全体のなすベクトル空間を H_k とすると、以上のように定義された内積に対し、以下の事実が成立する事が知られている。

特に 3 次元空間では次が成立する。

H_k が更に強い性質を満たすことも証明可能である。S^{n − 1} 上の自乗可積分函数全体の空間

L²(S^{n − 1}) = { f: S^{n − 1} → C | f は可測かつ ⟨ f | f ⟩_{S^{n − 1}} が有限値 }

は H_k を使って直交分解可能である^[11]：

これを言い換えると、以下の系が従う：

特に 3 次元の場合は、上述の事実と定理1から以下が成立する：

3 次元空間 R³ の球面座標 (r, θ, φ) に対し、

${\displaystyle \mathrm {d} x\,\mathrm {d} y\,\mathrm {d} z=r^{2}\sin \theta \,\mathrm {d} r\,\mathrm {d} \theta \,\mathrm {d} \varphi }$

が成立する。そこで、R 上の関数 χ, ξ に対し、χ, ξ の内積を

${\displaystyle \langle \chi \mid \xi \rangle _{R}:=\int _{-\infty }^{\infty }\chi (r)\xi (r)r^{2}\mathrm {d} r}$

(D1)

により定義し、さらに R³ の関数 f₁, f₂ の内積を

${\displaystyle \langle f_{1}\mid f_{2}\rangle :=\int _{\mathbf {R} ^{3}}f_{1}(\mathbf {x} )f_{2}(\mathbf {x} )\,\mathrm {d} x\,\mathrm {d} y\,\mathrm {d} z}$

(D2)

とする。f₁, f₂ が

f₁(x) = χ₁(r) Y₁(θ, φ), f₂(x) = χ₂(r) Y₂(θ, φ)

と変数分離形で書けていた場合には、(C1), (D1), (D2) で定義した内積は以下の性質を満たす。

${\displaystyle \langle f_{1}|f_{2}\rangle =\langle \chi _{1}|\chi _{2}\rangle _{R}\langle Y_{1}|Y_{2}\rangle _{S}}$

(C1), (D1), (D2) の内積を用いて自乗可積分な関数全体の集合をそれぞれ L²(S², sin θ dθ dφ), L²(R, r² dr), L²(R³, dx dy dz) と書くと、ヒルベルト空間の一般論から、次が成立する^{[注 3]}。

上述した定理と定理1から、以下の結論が従う。

いくつかの球面調和関数の具体的な表式を示す。

${\displaystyle {\begin{aligned}Y_{0}^{0}&={\frac {1}{\sqrt {4\pi }}}\\Y_{1}^{0}&={\sqrt {\frac {3}{4\pi }}}\,\cos \theta \\Y_{1}^{\pm 1}&=\mp {\sqrt {\frac {3}{8\pi }}}\,\sin \theta \,e^{\pm i\varphi }\\Y_{2}^{0}&={\sqrt {\frac {5}{16\pi }}}\,(3\cos ^{2}\theta -1)\\Y_{2}^{\pm 1}&=\mp {\sqrt {\frac {15}{8\pi }}}\,\sin \theta \cos \theta \,e^{\pm i\varphi }\\Y_{2}^{\pm 2}&={\sqrt {\frac {15}{32\pi }}}\,\sin ^{2}\theta \,e^{\pm 2i\varphi }\\Y_{3}^{0}&={\sqrt {\frac {7}{16\pi }}}\,(5\cos ^{3}\theta -3\cos \theta )\\Y_{3}^{\pm 1}&=\mp {\sqrt {\frac {21}{64\pi }}}\,\sin \theta \,(5\cos ^{2}\theta -1)\,e^{\pm i\varphi }\\Y_{3}^{\pm 2}&={\sqrt {\frac {105}{32\pi }}}\,\sin ^{2}\theta \cos \theta \,e^{\pm 2i\varphi }\\Y_{3}^{\pm 3}&=\mp {\sqrt {\frac {35}{64\pi }}}\,\sin ^{3}\theta \,e^{\pm 3i\varphi }\end{aligned}}}$

球面調和関数には「加法定理」と呼ばれる性質がある。これは三角関数における加法定理

${\displaystyle \cos(\theta '-\theta )=\cos \theta '\cos \theta +\sin \theta \sin \theta '}$

を一般化したものと捉えることができる。上式の右辺は球面調和関数に、左辺はルジャンドル多項式に置き換えられる。

二つの単位ベクトル x および y を考え、それらの球面座標をそれぞれ (θ, φ) および (θ′, φ′) とする。このとき、加法定理は以下のように表すことができる^[12]：

${\displaystyle P_{\ell }({\boldsymbol {x}}\cdot {\boldsymbol {y}})={\frac {4\pi }{2\ell +1}}\sum _{m=-\ell }^{\ell }Y_{\ell m}(\theta ',\varphi ')\,Y_{\ell m}^{*}(\theta ,\varphi ).}$

(1)

ここで P_ℓ は ℓ 次のルジャンドル多項式である。この表式は実数調和関数・虚数調和関数の双方について成り立つ^{[注 4]}。この結果は単位球面上のポアソン核の性質を用いて、あるいはベクトル y を z 軸に沿うように幾何的に回転させたのちに右辺を直接計算することにより解析的に証明することができる^[13]。

特に、x = y の場合はウンゼルトの定理^[14]

${\displaystyle \sum _{m=-\ell }^{\ell }Y_{\ell m}^{*}(\theta ,\varphi )\,Y_{\ell m}(\theta ,\varphi )={\frac {2\ell +1}{4\pi }}}$

に帰着する。この式は一次元の三角関数における恒等式 cos² θ + sin²θ = 1 を二次元に拡張したものとみなすことができる。

式 (1) の左辺 P_ℓ(x⋅y) は ℓ 次の帯球調和関数（英語版）の定数倍である。この観点から、より高次元の場合にも次のように一般化することができる。Y_j を n 次元超球面上の ℓ 次の球面調和関数の張る空間 H_ℓ の任意の正規直交基底とする。このとき、単位ベクトル x に対応する ℓ 次の帯球調和関数 Z (ℓ)
x  は以下のように書き下せる^[15]。

${\displaystyle Z_{\boldsymbol {x}}^{(\ell )}({\boldsymbol {y}})=\sum _{j=1}^{\dim(\mathbf {H} _{\ell })}{\overline {Y_{j}({\boldsymbol {x}})}}\,Y_{j}({\boldsymbol {y}}).}$

(2)

さらに、帯球調和関数 Z (ℓ)
x (y) は適切なゲーゲンバウアー多項式の定数倍として表すことができる：

${\displaystyle Z_{\boldsymbol {x}}^{(\ell )}({\boldsymbol {y}})=C_{\ell }^{((n-1)/2)}({\boldsymbol {x}}\cdot {\boldsymbol {y}}).}$

(3)

x および y が球面座標で表される場合、(2) および (3) を組み合わせると (1) が得られる。最後に、x = y の場合を評価すると次の恒等式が得られる：

${\displaystyle {\frac {\dim \mathbf {H} _{\ell }}{\omega _{n-1}}}=\sum _{j=1}^{\dim(\mathbf {H} _{\ell })}|Y_{j}({\boldsymbol {x}})|^{2}.}$

ここで ω_{n − 1} は (n − 1) 次元超球の体積である。

クレブシュ–ゴルダン係数とは、二つの球面調和関数の積を球面調和関数の線形結合で展開する際の展開係数である。ウィグナーの3-j記号やラカー係数、スレーター積分など様々な計算方法があるが、本質は同じである。抽象的には、クレブシュ–ゴルダン係数は二つの回転群の既約表現のテンソル積を既約表現の和で表わすときの係数と見ることができる。よって、適切に正規化すれば多重度と一致する。

原点に対する点対称操作で符号が替わらない（偶関数）かあるいは符号が逆になる（奇関数）かに依って、球面調和関数に対する「パリティ」が定義される。原点を不動点とする点対称操作は PΨ(r→) = Ψ(−r→) と表わせる。立体角で表わせば、{θ, φ} を {π − θ, π + φ} に置き換える操作になる。ルジャンドル陪多項式（Associated Legendre polynomials）はパリティとして (−1)^{ℓ + m} を、指数関数は (−1)^m を与えるので、両者を併せると球面調和関数のパリティは（mには依らずに） (−1)^ℓ となる。

${\displaystyle Y_{\ell }^{m}(\theta ,\phi )\;\rightarrow \;Y_{\ell }^{m}(\pi -\theta ,\pi +\phi )=(-1)^{\ell }Y_{\ell }^{m}(\theta ,\phi )}$

このことは、高次元に一般化した場合にも成り立つ。ℓ 次の球面調和関数に点対称操作を施した場合、符号の変化は (−1)^ℓ となる。 （これは調和多項式が次数の偶・奇に併せて空間反転で偶関数・奇関数であること、球面調和函数が調和多項式の球面上への制限であることからも容易に理解できる。）

量子力学で、球対称なポテンシャル V(r) に対する1粒子シュレーディンガー方程式（代表的なものは水素原子のシュレーディンガー方程式）

${\displaystyle \left\{-{\frac {\hbar ^{2}}{2m}}\nabla ^{2}+V(r)\right\}\psi ({\boldsymbol {r}})=E\psi ({\boldsymbol {r}})}$

を解いたときに、球面調和関数が現れる。量子力学では Y m
ℓ  の ℓ, m を量子数と呼び、それぞれ ℓ を方位量子数、m を磁気量子数という。

球面調和関数は軌道角運動量 ℓ と密接な関係がある。球面調和関数は ℓ² と ℓ_z の同時固有関数になっており、その固有値はそれぞれ ħ²ℓ(ℓ + 1), mħ である。すなわち

${\displaystyle {\boldsymbol {\ell }}^{2}Y_{\ell }^{m}=\hbar ^{2}\ell (\ell +1)Y_{\ell }^{m}}$

${\displaystyle \ell _{z}Y_{\ell }^{m}=m\hbar Y_{\ell }^{m}}$

となる。また、上昇下降演算子 ℓ₊, ℓ₋ を球面調和関数に作用させると

${\displaystyle \ell _{+}Y_{\ell }^{m}=\hbar {\sqrt {(\ell -m)(\ell +m+1)}}\,Y_{\ell }^{m+1}}$

${\displaystyle \ell _{-}Y_{\ell }^{m}=\hbar {\sqrt {(\ell +m)(\ell -m+1)}}\,Y_{\ell }^{m-1}}$

${\displaystyle \ell _{+}Y_{\ell }^{\ell }=0,\quad \ell _{-}Y_{\ell }^{-\ell }=0}$

となる。

• Jean Gallier (Department of Computer and Information Science University of Pennsylvania) (2013年). “Notes on Spherical Harmonics and Linear Representations of Lie Groups” (PDF). ペンシルバニア大学. 2017年8月29日閲覧。
• 野村隆昭 (2006年). “極座標・回転群・SL(2, R)” (PDF). 九州大学. 2017年1月4日閲覧。
• Brian C.Hall (July 1, 2013). Quantum Theory for Mathematicians. Graduate Texts in Mathematics 267. Springer 
• 高知大学自然科学系 田部井隆雄、神奈川県温泉地学研究所 里村幹夫、京都大学大学院理学研究科 福田洋一 (2004年). “4-4. ルジャンドルの多項式, 陪多項式”. 日本測地学会. 2017年1月4日閲覧。
• Watson, G. N.; Whittaker, E. T. (1927), A Course of Modern Analysis, Cambridge University Press, p. 392 
• Unsöld, Albrecht (1927). “Beiträge zur Quantenmechanik der Atome”. Annalen der Physik 387 (3): 355–393. Bibcode: 1927AnP...387..355U. doi:10.1002/andp.19273870304. ISSN 0003-3804. LCCN 50-13519. OCLC 5854993. 
• Stein, Elias; Weiss, Guido (November 1, 1971). Introduction to Fourier Analysis on Euclidean Spaces. Princeton mathematical series. Princeton, N.J.: Princeton University Press. ASIN 069108078X. ISBN 978-0-691-08078-9. NCID BA82681515. OCLC 919508312 

• 小出昭一郎『量子力学1』（改訂版）裳華房〈基礎物理学選書〉、1990年10月5日、89-96頁。ASIN 4785321326。ISBN 4-7853-2132-6。 NCID BN05389383。OCLC 835016094。全国書誌番号:91005557。 
• L. I. Schiff (1968) [1955]. Quantum Mechanics (3rd ed.). Singapore etc.: McGraw Hill. pp. 79-80. ASIN 0070856435. ISBN 0-07-085643-5. NCID BA1086214X. OCLC 632275975 
• Christian Helanow (2009年). “Spherical harmonics: a theoretical and graphical study” (PDF). 2017年1月4日閲覧。
• Joṥe Alvarado (2007年12月4日). “Group Theoretical Aspects of Quantum Mechanics” (PDF). 2016年12月1日閲覧。
• 野村隆昭:「球面調和函数と群の表現」、ISBN: 978-4535798182、日本評論社 (2018年7月)。
• Edmonds, A. R.: "Angular Momentum in Quantum Mechanics",　Princeton University Press, ISBN 978-0691025896 (1996). Reprint version.
• 山内恭彦：「回転群とその表現」、岩波書店、ISBN 978-4000051460 (1988年11月)。
• MartinJ. Mohlenkamp: "A FastTransform for Spherical Harmonics", The Journal of Fourier Analysis and Applications 5(2/3), pp.159-184 (1999).
• Kendall Atkinson and Weimin Han: "Spherical Harmonics and Approximations on the Unit Sphere: An Introduction", Springer, ISBN 978-3-642-25982-1 (2012).

ウィキメディア・コモンズには、球面調和関数に関連するカテゴリがあります。

• 調和関数

• ブリタニカ国際大百科事典 小項目事典『球関数』 - コトバンク
• Spherical harmonic - ブリタニカ百科事典