球面上で定義された特殊な数学関数
最初のいくつかの実球面調和関数の視覚的表現。青い部分は関数が正の領域を表し、黄色の部分は関数が負の領域を表します。原点からの表面の距離は、 角度方向の の絶対値を示します 。
はい
ℓ
メートル
(
θ
、
φ
)
{\displaystyle Y_{\ell }^{m}(\theta ,\varphi )}
(
θ
、
φ
)
{\displaystyle (\theta ,\varphi )}
数学 と 物理科学 において 、 球面調和関数は 球面 の表面上で定義される 特殊な関数 です。 多くの科学分野で 偏微分方程式を 解く際によく使用されます。 球面調和関数の表には、 一般的な球面調和関数のリストが含まれています。
球面調和関数は完全な直交関数 の集合 、つまり 正規直交基底 を形成するため、球面上で定義された各関数は、これらの球面調和関数の和として表すことができます。これは、円上で定義された 周期関数が、 フーリエ級数を介して 円関数 (正弦関数と余弦関数) の和として表されるのと似ています 。フーリエ級数の正弦関数と余弦関数のように、球面調和関数は、 右の図の関数の行に見られるように、 (空間) 角周波数 によって整理できます。さらに、球面調和関数は、3 次元の回転群 である SO(3) の 既約表現の 基底関数 であるため、SO(3) の 群論的 議論で中心的な役割を果たします 。
球面調和関数は、 球面領域で ラプラス方程式を解くことから生じる。ラプラス方程式の解である関数は 調和関数と呼ばれる。その名前にもかかわらず、球面調和関数は 直交座標 で最も単純な形を取り 、そこではラプラス方程式 に従う 次数 の 同次多項式 として定義できる。 同次性を使用して 前述の 次数の多項式から半径方向の依存性の因子を抽出すると、 球面座標 との関連がすぐに生じる。残りの因子は、球面角座標 と のみの関数、または同等にこれらの角度で指定される 方向 単位ベクトル の関数と見なすことができる。この設定では、それらは3次元のラプラス方程式の解の集合の角度部分と見なすことができ、この観点はしばしば代替定義として採用される。ただし、球面調和関数は、 球面上の標準的な円形計量に対する ラプラス-ベルトラミ演算子 に関して調和的な球面上の関数では ないことに注意してください。調和関数は 最大原理 を満たすため、この意味で球面上の唯一の調和関数は定数です。球面調和関数は、球面上の関数として、ラプラス-ベルトラミ演算子の固有関数です(高次元を参照)。
ℓ
{\displaystyle \ell}
(
x
、
ええ
、
ず
)
{\displaystyle (x,y,z)}
r
ℓ
{\displaystyle r^{\ell}}
ℓ
{\displaystyle \ell}
θ
{\displaystyle \theta}
φ
{\displaystyle \varphi}
r
{\displaystyle \mathbf {r} }
またはで 示される特定の球面調和関数の集合は、 1782年に ピエール・シモン・ド・ラプラス によって初めて導入されたため、ラプラスの球面調和関数として知られています。 [1]これらの関数は 直交 系を形成し 、したがって上で言及したように球面上の一般関数の展開の基礎となります。
はい
ℓ
メートル
(
θ
、
φ
)
{\displaystyle Y_{\ell }^{m}(\theta ,\varphi )}
はい
ℓ
メートル
(
r
)
{\displaystyle Y_{\ell }^{m}({\mathbf {r} })}
球面調和関数は、 多極静 電場や 電磁場 、 電子配置 、 重力場 、 ジオイド 、惑星や恒星の 磁場 、 宇宙マイクロ波背景放射などの表現を含む、多くの理論的および実際的な応用において重要です。3D コンピュータ グラフィックス では 、球面調和関数は、間接照明 ( アンビエント オクルージョン 、 グローバル イルミネーション 、 事前計算された放射伝達 など) や 3D 形状のモデリングなど、さまざまなトピックで役割を果たします。
歴史
ピエール=シモン・ラプラス 、1749–1827
球面調和関数は、ニュートンの万有引力の法則のニュートン力学ポテンシャルと関連して、3次元で初めて研究されました。1782年、ピエール・シモン・ド・ラプラスは著書『天体力学 』 で 、 点 x i に ある 質点 m i の 集合に関連する 点 x の重力ポテンシャルは 、次のように与えられると決定 しました。
R
3
→
R
{\displaystyle \mathbb {R} ^{3}\to \mathbb {R} }
五
(
x
)
=
∑
私
メートル
私
|
x
私
−
x
|
。
{\displaystyle V(\mathbf {x} )=\sum _{i}{\frac {m_{i}}{|\mathbf {x} _{i}-\mathbf {x} |}}.}
上記の合計の各項は、質点の個々のニュートンポテンシャルです。その直前に、 アドリアン・マリー・ルジャンドルは、ニュートンポテンシャルの r = | x | および r 1 = | x 1 | の累乗展開を調査しました 。彼は、 r ≤ r 1 の
場合、
1
|
x
1
−
x
|
=
ポ
0
(
コス
γ
)
1
r
1
+
ポ
1
(
コス
γ
)
r
r
1
2
+
ポ
2
(
コス
γ
)
r
2
r
1
3
+
⋯
{\displaystyle {\frac {1}{|\mathbf {x} _{1}-\mathbf {x} |}}=P_{0}(\cos \gamma ){\frac {1}{r_{1 }}}+P_{1}(\cos \gamma ){\frac {r}{r_{1}^{2}}}+P_{2}(\cos \gamma ){\frac {r^{2}}{r_{1}^{3}}}+\cdots }
ここで、 γ は ベクトル x と x 1 の 間の角度です。関数は ルジャンドル多項式 であり 、球面調和関数の特殊なケースとして導出できます。その後、ラプラスは 1782 年の回想録で、球面座標を使用してこれらの係数を調査し、 x 1 と x の 間の角度 γ を表しました。( 詳細については、
ルジャンドル多項式 § 応用を参照してください。)
ポ
私
:
[
−
1
、
1
]
→
R
{\displaystyle P_{i}:[-1,1]\to \mathbb {R} }
1867 年、 ウィリアム トムソン (ケルビン卿) と ピーター ガスリー テイトは、 自然哲学論文で 固体 球面調和関数 を導入し 、これらの関数に「球面調和関数」という名前を初めて導入しました。 固体調和関数は、 ラプラス方程式 の 同次 多項式解 でした。
ラプラス方程式を球座標で調べることで、トムソンとテイトはラプラスの球面調和関数を復元しました (調和多項式表現を参照)。「ラプラス係数」という用語は、ウィリアム ヒューウェル がこれらの線に沿って導入された特定の解のシステムを説明するために使用しましたが、他の人はこの呼称を、ラプラスとルジャンドルによって適切に導入された ゾーン球面調和関数 のために留保しました。
R
3
→
R
{\displaystyle \mathbb {R} ^{3}\to \mathbb {R} }
∂
2
あなた
∂
x
2
+
∂
2
あなた
∂
ええ
2
+
∂
2
あなた
∂
ず
2
=
0.
{\displaystyle {\frac {\partial ^{2}u}{\partial x^{2}}}+{\frac {\partial ^{2}u}{\partial y^{2}}}+{\frac {\partial ^{2}u}{\partial z^{2}}}=0.}
19 世紀に フーリエ級数が開発されたことで、 熱方程式 や 波動方程式 などのさまざまな物理的問題を長方形領域で解くことが可能になった。これは、関数を 三角関数の級数 に展開することで実現できる 。フーリエ級数の三角関数が弦の振動の基本モードを表すのに対し 、 球面調和関数はほぼ同じように 球の振動 の基本モードを表す。フーリエ級数の理論の多くの側面は、三角関数ではなく球面調和関数で展開することで一般化できる。さらに、三角関数が 複素指数関数 として同等に記述できるのと同様に、球面調和関数も複素数値関数として同等の形を持っていた。これは、もともとラプラスとルジャンドルが研究した天体力学などの
球対称性 を持つ問題にとって恩恵であった。
球面調和関数は物理学ですでに普及しており、20 世紀に量子力学が 誕生した際にその重要性が高まりました 。(複素数値の) 球面調和関数は、 軌道角運動量 演算子
の 2 乗の 固有関数 である
ため、 原子軌道 のさまざまな 量子化された 構成を表します。
S
2
→
C
{\displaystyle S^{2}\to \mathbb {C} }
−
私
ℏ
r
×
∇
、
{\displaystyle -i\hbar \mathbf {r} \times \nabla ,}
ラプラスの球面調和関数
(上から下) および(左から右) の 実数 (ラプラス) 球面調和関数 。ゾーン調和関数、セクター調和関数、テッセラル調和関数は、それぞれ左端の列、主対角線、およびその他の場所に描かれています。(負の次数の調和関数は、正の次数の調和関数に対して z 軸を中心 に回転して表示されます 。)
はい
ℓ
メートル
{\displaystyle Y_{\ell m}}
ℓ
=
0
、
…
、
4
{\displaystyle \ell =0,\dots ,4}
メートル
=
0
、
…
、
ℓ
{\displaystyle m=0,\dots ,\ell }
はい
ℓ
(
−
メートル
)
{\displaystyle Y_{\ell (-m)}}
90
∘
/
メートル
{\displaystyle 90^{\circ }/m}
実球面調和関数の別の図 。
はい
ℓ
メートル
{\displaystyle Y_{\ell m}}
ラプラス方程式は、 スカラー場 f のラプラシアン がゼロであることを規定している 。(ここで、スカラー場は複素数、すなわち(滑らかな)関数に対応するものと理解されている 。) 球面座標 では、これは次のようになる: [2]
ふ
:
R
3
→
C
{\displaystyle f:\mathbb {R} ^{3}\to \mathbb {C} }
∇
2
ふ
=
1
r
2
∂
∂
r
(
r
2
∂
ふ
∂
r
)
+
1
r
2
罪
θ
∂
∂
θ
(
罪
θ
∂
ふ
∂
θ
)
+
1
r
2
罪
2
θ
∂
2
ふ
∂
φ
2
=
0.
{\displaystyle \nabla^{2}f={\frac {1}{r^{2}}}{\frac {\partial }{\partial r}}\left(r^{2}{\frac {\partial f}{\partial r}}\right)+{\frac {1}{r^{2}\sin \theta }}{\frac {\partial }{\partial \theta }}\left(\sin \theta {\frac {\partial f}{\partial \theta }}\right)+{\frac {1}{r^{2}\sin ^{2}\theta }}{\frac {\partial ^{2}f}{\partial \varphi ^{2}}}=0.}
f ( r , θ , φ ) = R ( r ) Y ( θ , φ ) の形式の解を求める問題を考えてみましょう 。 変数分離 により、ラプラス方程式を適用することで2つの微分方程式が得られます。2番目の方程式は
、 Yが Y ( θ , φ ) = Θ( θ ) Φ( φ ) という形式である
と仮定して簡略化できます 。2番目の方程式に再び変数分離を適用すると、次の微分方程式のペアが得られます
。
1
R
d
d
r
(
r
2
d
R
d
r
)
=
λ
、
1
はい
1
罪
θ
∂
∂
θ
(
罪
θ
∂
はい
∂
θ
)
+
1
はい
1
罪
2
θ
∂
2
はい
∂
φ
2
=
−
λ
。
{\displaystyle {\frac {1}{R}}{\frac {d}{dr}}\left(r^{2}{\frac {dR}{dr}}\right)=\lambda ,\qquad {\frac {1}{Y}}{\frac {1}{\sin \theta }}{\frac {\partial }{\partial \theta }}\left(\sin \theta {\frac {\partial Y}{\partial \theta }}\right)+{\frac {1}{Y}}{\frac {1}{\sin ^{2}\theta }}{\frac {\partial ^{2}Y}{\partial \varphi ^{2}}}=-\lambda .}
1
Φ
d
2
Φ
d
φ
2
=
−
メートル
2
{\displaystyle {\frac {1}{\Phi }}{\frac {d^{2}\Phi }{d\varphi ^{2}}}=-m^{2}}
λ
罪
2
θ
+
罪
θ
Θ
d
d
θ
(
罪
θ
d
Θ
d
θ
)
=
メートル
2
{\displaystyle \lambda \sin ^{2}\theta +{\frac {\sin \theta }{\Theta }}{\frac {d}{d\theta }}\left(\sin \theta {\frac {d\Theta }{d\theta }}\right)=m^{2}}
となる。先験的 には、 m は 複素定数であるが、 Φ は 周期が 2 π を 均等に割り切る周期関数 でなければならないため 、 m は 必然的に整数となり、 Φ は 複素指数関数 e ± imφ の線形結合となる。解関数 Y ( θ , φ ) は、 θ = 0, π である球面の極で正則となる。 領域の境界点で第 2 方程式の 解 Θにこの正則性を課すことは、 シュトゥルム・リウヴィル問題 であり、パラメータ λは、 ℓ ≥ | m | となる負でない整数に対して λ = ℓ ( ℓ + 1) の形にならなければならない。これは、 軌道角運動量 の観点からも以下で説明される 。さらに、変数変換 t = cos θ により、この方程式はルジャンドル方程式 に変換され 、その解は、ルジャンドルの 随伴多項式 P メートル リットル (cos θ ) 。最後に、 R の方程式は R ( r ) = A r ℓ + B r − ℓ − 1 の形の解を持ちます。R 3 全体 にわたって解が正則であることが必要で、 B = 0 と なり ます。 [3]
ここで、解は特別な形式 Y ( θ , φ ) = Θ( θ ) Φ( φ )であると仮定しました。 ℓ の値が与えられた場合、この形式の独立した解は 2 ℓ + 1 個 あり、 − ℓ ≤ m ≤ ℓ となる整数 m ごとに 1 つずつあります 。 これらの角度解は 、ここでは 複素指数 関数として表される三角関数 と、ルジャンドルの随伴多項式
の積です。
はい
ℓ
メートル
:
S
2
→
C
{\displaystyle Y_{\ell}^{m}:S^{2}\to \mathbb {C} }
はい
ℓ
メートル
(
θ
、
φ
)
=
いいえ
e
私
メートル
φ
ポ
ℓ
メートル
(
コス
θ
)
{\displaystyle Y_{\ell }^{m}(\theta ,\varphi )=Ne^{im\varphi }P_{\ell }^{m}(\cos {\theta })}
満たす
r
2
∇
2
はい
ℓ
メートル
(
θ
、
φ
)
=
−
ℓ
(
ℓ
+
1
)
はい
ℓ
メートル
(
θ
、
φ
)
。
{\displaystyle r^{2}\nabla ^{2}Y_{\ell }^{m}(\theta ,\varphi )=-\ell (\ell +1)Y_{\ell }^{m}( \theta 、\varphi ).}
ここで は次数 ℓ 、次数 m の球面調和関数 と呼ばれ 、 は 準ルジャンドル多項式 、 N は正規化定数、 [4] θ と φ は それぞれ 余緯度と経度を表します。特に、 余緯度 θ 、つまり極角は、 北極の 0から赤道の π /2 、 南極の π までの範囲にあり、経度 φ 、つまり 方位角は、 0 ≤ φ < 2 π のすべての値を取ることができます。固定整数 ℓ に対して 、 固有値問題の
すべての解 Y ( θ , φ ) 、 、 は の 線形結合
です 。実際、そのような解のいずれについても、 r ℓ Y ( θ , φ )は、調和的な 同次多項式 (以下を参照)の球座標での表現であり、次元を数えると、そのような多項式が 2 ℓ + 1 線形独立で
あることがわかります。
はい
ℓ
メートル
:
S
2
→
C
{\displaystyle Y_{\ell}^{m}:S^{2}\to \mathbb {C} }
ポ
ℓ
メートル
:
[
−
1
、
1
]
→
R
{\displaystyle P_{\ell}^{m}:[-1,1]\to \mathbb {R} }
はい
:
S
2
→
C
{\displaystyle Y:S^{2}\to \mathbb {C} }
r
2
∇
2
はい
=
−
ℓ
(
ℓ
+
1
)
はい
{\displaystyle r^{2}\nabla ^{2}Y=-\ell (\ell +1)Y}
はい
ℓ
メートル
:
S
2
→
C
{\displaystyle Y_{\ell}^{m}:S^{2}\to \mathbb {C} }
R
3
→
C
{\displaystyle \mathbb {R} ^{3}\to \mathbb {C} }
原点を中心とする球における ラプラス方程式 の 一般解は、 球面調和関数の 線形結合に適切なスケール係数 r ℓ を乗じたものである。
ふ
:
R
3
→
C
{\displaystyle f:\mathbb {R} ^{3}\to \mathbb {C} }
Δ
ふ
=
0
{\displaystyle \Delta f=0}
ふ
(
r
、
θ
、
φ
)
=
∑
ℓ
=
0
∞
∑
メートル
=
−
ℓ
ℓ
ふ
ℓ
メートル
r
ℓ
はい
ℓ
メートル
(
θ
、
φ
)
、
{\displaystyle f(r,\theta ,\varphi )=\sum _{\ell =0}^{\infty }\sum _{m=-\ell }^{\ell }f_{\ell }^{m}r^{\ell }Y_{\ell }^{m}(\theta ,\varphi ),}
ここで、は 定数であり、因子 r ℓ Y ℓ m は (正則 ) 固体調和関数 として知られている 。このような展開は 球面上で有効である。
f
ℓ
m
∈
C
{\displaystyle f_{\ell }^{m}\in \mathbb {C} }
R
3
→
C
{\displaystyle \mathbb {R} ^{3}\to \mathbb {C} }
r
<
R
=
1
lim sup
ℓ
→
∞
|
f
ℓ
m
|
1
/
ℓ
.
{\displaystyle r<R={\frac {1}{\limsup _{\ell \to \infty }|f_{\ell }^{m}|^{{1}/{\ell }}}}.}
については 、代わりに の負のべきを持つ固体調和関数 ( 不規則 固体調和関数 )が選択されます。その場合、項を一致させて級数展開係数 を見つけるために、上で使用した テイラー級数 ( について ) の代わりに、 ローラン級数 ( について)で既知の領域の解を展開する必要があります 。
r
>
R
{\displaystyle r>R}
r
{\displaystyle r}
R
3
∖
{
0
}
→
C
{\displaystyle \mathbb {R} ^{3}\setminus \{\mathbf {0} \}\to \mathbb {C} }
r
=
∞
{\displaystyle r=\infty }
r
=
0
{\displaystyle r=0}
f
ℓ
m
∈
C
{\displaystyle f_{\ell }^{m}\in \mathbb {C} }
軌道角運動量
量子力学では、ラプラスの球面調和関数は 軌道角運動量の観点から理解されています [5]
ħ は 量子力学では慣例であり、 ħ = 1 の 単位で作業すると便利です 。球面調和関数は軌道角運動量の2乗の固有関数です
。ラプラスの球面調和関数は、軌道角運動量の2乗と方位軸の周りの回転の生成子の結合固有関数です。
L
=
−
i
ℏ
(
x
×
∇
)
=
L
x
i
+
L
y
j
+
L
z
k
.
{\displaystyle \mathbf {L} =-i\hbar (\mathbf {x} \times \mathbf {\nabla } )=L_{x}\mathbf {i} +L_{y}\mathbf {j} +L_{z}\mathbf {k} .}
L
2
=
−
r
2
∇
2
+
(
r
∂
∂
r
+
1
)
r
∂
∂
r
=
−
1
sin
θ
∂
∂
θ
sin
θ
∂
∂
θ
−
1
sin
2
θ
∂
2
∂
φ
2
.
{\displaystyle {\begin{aligned}\mathbf {L} ^{2}&=-r^{2}\nabla ^{2}+\left(r{\frac {\partial }{\partial r}}+1\right)r{\frac {\partial }{\partial r}}\\&=-{\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}}}.\end{aligned}}}
L
z
=
−
i
(
x
∂
∂
y
−
y
∂
∂
x
)
=
−
i
∂
∂
φ
.
{\displaystyle {\begin{aligned}L_{z}&=-i\left(x{\frac {\partial }{\partial y}}-y{\frac {\partial }{\partial x}}\right)\\&=-i{\frac {\partial }{\partial \varphi }}.\end{aligned}}}
これらの演算子は可換であり、 R 3 上の重み関数として 正規分布 に関して二乗積分可能な 関数 fの 重み付き ヒルベルト空間 上の 稠密に定義された 自己随伴演算子 です。
さらに、 L 2は 正の演算子 です 。
1
(
2
π
)
3
/
2
∫
R
3
|
f
(
x
)
|
2
e
−
|
x
|
2
/
2
d
x
<
∞
.
{\displaystyle {\frac {1}{(2\pi )^{3/2}}}\int _{\mathbb {R} ^{3}}|f(x)|^{2}e^{-|x|^{2}/2}\,dx<\infty .}
Y が L 2 と L z の結合固有関数である 場合 、定義により、
いくつかの実数 m および λ に対して成り立ちます。ここで、 m は 実際には整数でなければなりません。これは、 Y が 座標 φにおいて周期が 2 π を 均等に割り切る数でなければならないためです 。さらに、
および L x 、 L y 、 L z
のそれぞれが 自己随伴であるため、 λ ≥ m 2 となります。
L
2
Y
=
λ
Y
L
z
Y
=
m
Y
{\displaystyle {\begin{aligned}\mathbf {L} ^{2}Y&=\lambda Y\\L_{z}Y&=mY\end{aligned}}}
L
2
=
L
x
2
+
L
y
2
+
L
z
2
{\displaystyle \mathbf {L} ^{2}=L_{x}^{2}+L_{y}^{2}+L_{z}^{2}}
この結合固有空間をE λ , m で表し 、 上昇演算子と下降演算子 を
次 のように
定義します。すると、 L + と L − は L 2 と可換となり、 L + , L − , L z によって生成されるリー代数は 、交換関係を持つ
2 次 特殊線型リー代数 となります。
したがって、 L + : E λ , m → E λ , m +1 (これは「上昇演算子」です) および L − : E λ , m → E λ , m −1 (これは「下降演算子」です) となります。特に、 L
L
+
=
L
x
+
i
L
y
L
−
=
L
x
−
i
L
y
{\displaystyle {\begin{aligned}L_{+}&=L_{x}+iL_{y}\\L_{-}&=L_{x}-iL_{y}\end{aligned}}}
s
l
2
(
C
)
{\displaystyle {\mathfrak {sl}}_{2}(\mathbb {C} )}
[
L
z
,
L
+
]
=
L
+
,
[
L
z
,
L
−
]
=
−
L
−
,
[
L
+
,
L
−
]
=
2
L
z
.
{\displaystyle [L_{z},L_{+}]=L_{+},\quad [L_{z},L_{-}]=-L_{-},\quad [L_{+},L_{-}]=2L_{z}.}
か + : E λ , m → E λ , m + k は 、 k が十分に大きい場合、ゼロでなければなりません。 これ は、不等式 λ ≥ m 2 が各非自明な結合固有空間で成立する必要があるためです。Y ∈ E λ , m を 非ゼロの結合固有関数とし、 k を 次の式を満たす最小の整数と
します。
すると
、次の
式が成り立つ
ため 、正の整数 ℓ = m + kに対して λ = ℓ ( ℓ + 1) と なります。
L
+
k
Y
=
0.
{\displaystyle L_{+}^{k}Y=0.}
L
−
L
+
=
L
2
−
L
z
2
−
L
z
{\displaystyle L_{-}L_{+}=\mathbf {L} ^{2}-L_{z}^{2}-L_{z}}
0
=
L
−
L
+
k
Y
=
(
λ
−
(
m
+
k
)
2
−
(
m
+
k
)
)
Y
.
{\displaystyle 0=L_{-}L_{+}^{k}Y=(\lambda -(m+k)^{2}-(m+k))Y.}
これまではすべて球面座標表現で説明してきました が、完全な正規 直交球面ケット基底 でより抽象的に表現することもできます。
⟨
θ
,
φ
|
l
m
⟩
=
Y
l
m
(
θ
,
φ
)
{\displaystyle \langle \theta ,\varphi |lm\rangle =Y_{l}^{m}(\theta ,\varphi )}
調和多項式表現
球面調和関数は、特定の多項式関数 の単位球面への制限として表現できます。具体的には、
すべての実数 およびすべてのに対して で ある場合、
(複素数値)多項式関数は 次数の 同次で あるといえます。 が ラプラシアン である
場合 、
は 調和的で ある といえます 。次に、各 に対して 、次のように定義します。
R
3
→
C
{\displaystyle \mathbb {R} ^{3}\to \mathbb {C} }
p
:
R
3
→
C
{\displaystyle p:\mathbb {R} ^{3}\to \mathbb {C} }
ℓ
{\displaystyle \ell }
p
(
λ
x
)
=
λ
ℓ
p
(
x
)
{\displaystyle p(\lambda \mathbf {x} )=\lambda ^{\ell }p(\mathbf {x} )}
λ
∈
R
{\displaystyle \lambda \in \mathbb {R} }
x
∈
R
3
{\displaystyle x\in \mathbb {R} ^{3}}
p
{\displaystyle p}
Δ
p
=
0
,
{\displaystyle \Delta p=0,}
Δ
{\displaystyle \Delta }
ℓ
{\displaystyle \ell }
A
ℓ
=
{
harmonic polynomials
R
3
→
C
that are homogeneous of degree
ℓ
}
.
{\displaystyle \mathbf {A} _{\ell }=\left\{{\text{harmonic polynomials }}\mathbb {R} ^{3}\to \mathbb {C} {\text{ that are homogeneous of degree }}\ell \right\}.}
たとえば、 のとき 、 はすべての線形関数 の 3 次元空間になります 。これは、そのような関数はすべて自動的に調和関数になるからです。一方、 のときは 、5 次元空間になります。
ℓ
=
1
{\displaystyle \ell =1}
A
1
{\displaystyle \mathbf {A} _{1}}
R
3
→
C
{\displaystyle \mathbb {R} ^{3}\to \mathbb {C} }
ℓ
=
2
{\displaystyle \ell =2}
A
2
=
span
C
(
x
1
x
2
,
x
1
x
3
,
x
2
x
3
,
x
1
2
−
x
2
2
,
x
1
2
−
x
3
2
)
.
{\displaystyle \mathbf {A} _{2}=\operatorname {span} _{\mathbb {C} }(x_{1}x_{2},\,x_{1}x_{3},\,x_{2}x_{3},\,x_{1}^{2}-x_{2}^{2},\,x_{1}^{2}-x_{3}^{2}).}
任意の に対して、 次数の球面調和関数の 空間は、 の元の 球面への制限の空間に他なりません 。 [6]序論で示唆したように、この観点はおそらく「球面調和関数」(すなわち、 調和関数 の球面への制限 )という用語の起源です。
ℓ
{\displaystyle \ell }
H
ℓ
{\displaystyle \mathbf {H} _{\ell }}
ℓ
{\displaystyle \ell }
S
2
{\displaystyle S^{2}}
A
ℓ
{\displaystyle \mathbf {A} _{\ell }}
たとえば、任意の に対して、 公式は
の 定義域と共定義域を持つ
次数の同次多項式を定義します が、これは とは独立です 。この多項式は調和多項式であることが簡単にわかります。 を 球座標で書き 、 に制限すると、次
が
得られます。
これは と書き直すことができます。 ルジャンドル従属多項式
の公式を使用した後 、これを球面調和関数の公式として認識できます [7] (「特殊なケース」を参照)。
c
∈
C
{\displaystyle c\in \mathbb {C} }
p
(
x
1
,
x
2
,
x
3
)
=
c
(
x
1
+
i
x
2
)
ℓ
{\displaystyle p(x_{1},x_{2},x_{3})=c(x_{1}+ix_{2})^{\ell }}
ℓ
{\displaystyle \ell }
R
3
→
C
{\displaystyle \mathbb {R} ^{3}\to \mathbb {C} }
x
3
{\displaystyle x_{3}}
p
{\displaystyle p}
(
r
,
θ
,
φ
)
{\displaystyle (r,\theta ,\varphi )}
r
=
1
{\displaystyle r=1}
p
(
θ
,
φ
)
=
c
sin
(
θ
)
ℓ
(
cos
(
φ
)
+
i
sin
(
φ
)
)
ℓ
,
{\displaystyle p(\theta ,\varphi )=c\sin(\theta )^{\ell }(\cos(\varphi )+i\sin(\varphi ))^{\ell },}
p
(
θ
,
φ
)
=
c
(
1
−
cos
2
(
θ
)
)
ℓ
e
i
ℓ
φ
.
{\displaystyle p(\theta ,\varphi )=c\left({\sqrt {1-\cos ^{2}(\theta )}}\right)^{\ell }e^{i\ell \varphi }.}
P
ℓ
ℓ
{\displaystyle P_{\ell }^{\ell }}
Y
ℓ
ℓ
(
θ
,
φ
)
.
{\displaystyle Y_{\ell }^{\ell }(\theta ,\varphi ).}
コンベンション
直交性と正規化
ラプラス球面調和関数 では、いくつかの異なる正規化が一般的に使用さ れています。このセクション全体で、 に対して という標準的な規則を使用します ( ルジャンドル多項式 を 参照)
。これは、ロドリゲスの公式によって与えられる自然な正規化です。
S
2
→
C
{\displaystyle S^{2}\to \mathbb {C} }
m
>
0
{\displaystyle m>0}
P
ℓ
−
m
=
(
−
1
)
m
(
ℓ
−
m
)
!
(
ℓ
+
m
)
!
P
ℓ
m
{\displaystyle P_{\ell }^{-m}=(-1)^{m}{\frac {(\ell -m)!}{(\ell +m)!}}P_{\ell }^{m}}
から まで の 複素平面における 球面調和関数のプロット(色 付き ) (Mathematica 13.1 関数 ComplexPlot3D で作成)
Y
ℓ
m
(
θ
,
φ
)
{\displaystyle Y_{\ell }^{m}(\theta ,\varphi )}
ℓ
=
2
{\displaystyle \ell =2}
m
=
1
{\displaystyle m=1}
φ
=
π
{\displaystyle \varphi =\pi }
−
2
−
2
i
{\displaystyle -2-2i}
2
+
2
i
{\displaystyle 2+2i}
音響学 で は [8] ラプラス球面調和関数は一般的に次のように定義されます(これはこの記事で使用されている慣例です)。
一方、 量子力学
では次のように定義されます。 [9] [10]
Y
ℓ
m
(
θ
,
φ
)
=
(
2
ℓ
+
1
)
4
π
(
ℓ
−
m
)
!
(
ℓ
+
m
)
!
P
ℓ
m
(
cos
θ
)
e
i
m
φ
{\displaystyle Y_{\ell }^{m}(\theta ,\varphi )={\sqrt {{\frac {(2\ell +1)}{4\pi }}{\frac {(\ell -m)!}{(\ell +m)!}}}}\,P_{\ell }^{m}(\cos {\theta })\,e^{im\varphi }}
Y
ℓ
m
(
θ
,
φ
)
=
(
−
1
)
m
(
2
ℓ
+
1
)
4
π
(
ℓ
−
m
)
!
(
ℓ
+
m
)
!
P
ℓ
m
(
cos
θ
)
e
i
m
φ
{\displaystyle Y_{\ell }^{m}(\theta ,\varphi )=(-1)^{m}{\sqrt {{\frac {(2\ell +1)}{4\pi }}{\frac {(\ell -m)!}{(\ell +m)!}}}}\,P_{\ell }^{m}(\cos {\theta })\,e^{im\varphi }}
ここで、 コンドン・ショートリー位相を除いたルジャンドル多項式です(位相を2回カウントするのを避けるため)。
P
ℓ
m
{\displaystyle P_{\ell }^{m}}
どちらの定義でも、球面調和関数は正規直交であり
、 δ ij は クロネッカーデルタ 、 d Ω = sin( θ ) dφ dθ である。この正規化は確率が正規化されることを保証するため、量子力学で使用される。
∫
θ
=
0
π
∫
φ
=
0
2
π
Y
ℓ
m
Y
ℓ
′
m
′
∗
d
Ω
=
δ
ℓ
ℓ
′
δ
m
m
′
,
{\displaystyle \int _{\theta =0}^{\pi }\int _{\varphi =0}^{2\pi }Y_{\ell }^{m}\,Y_{\ell '}^{m'}{}^{*}\,d\Omega =\delta _{\ell \ell '}\,\delta _{mm'},}
∫
|
Y
ℓ
m
|
2
d
Ω
=
1.
{\displaystyle \int {|Y_{\ell }^{m}|^{2}d\Omega }=1.}
測地学 [11] とスペクトル解析
の分野では、
Y
ℓ
m
(
θ
,
φ
)
=
(
2
ℓ
+
1
)
(
ℓ
−
m
)
!
(
ℓ
+
m
)
!
P
ℓ
m
(
cos
θ
)
e
i
m
φ
{\displaystyle Y_{\ell }^{m}(\theta ,\varphi )={\sqrt {{(2\ell +1)}{\frac {(\ell -m)!}{(\ell +m)!}}}}\,P_{\ell }^{m}(\cos {\theta })\,e^{im\varphi }}
単位力を持つ
1
4
π
∫
θ
=
0
π
∫
φ
=
0
2
π
Y
ℓ
m
Y
ℓ
′
m
′
∗
d
Ω
=
δ
ℓ
ℓ
′
δ
m
m
′
.
{\displaystyle {\frac {1}{4\pi }}\int _{\theta =0}^{\pi }\int _{\varphi =0}^{2\pi }Y_{\ell }^{m}\,Y_{\ell '}^{m'}{}^{*}d\Omega =\delta _{\ell \ell '}\,\delta _{mm'}.}
対照的に、磁気 学 [11] の分野では、シュミットの半正規化高調波
Y
ℓ
m
(
θ
,
φ
)
=
(
ℓ
−
m
)
!
(
ℓ
+
m
)
!
P
ℓ
m
(
cos
θ
)
e
i
m
φ
{\displaystyle Y_{\ell }^{m}(\theta ,\varphi )={\sqrt {\frac {(\ell -m)!}{(\ell +m)!}}}\,P_{\ell }^{m}(\cos {\theta })\,e^{im\varphi }}
正規化されている
∫
θ
=
0
π
∫
φ
=
0
2
π
Y
ℓ
m
Y
ℓ
′
m
′
∗
d
Ω
=
4
π
(
2
ℓ
+
1
)
δ
ℓ
ℓ
′
δ
m
m
′
.
{\displaystyle \int _{\theta =0}^{\pi }\int _{\varphi =0}^{2\pi }Y_{\ell }^{m}\,Y_{\ell '}^{m'}{}^{*}d\Omega ={\frac {4\pi }{(2\ell +1)}}\delta _{\ell \ell '}\,\delta _{mm'}.}
量子力学でもこの正規化が時々使用され、 ジュリオ・ラカー にちなんでラカーの正規化と名付けられています。
上記の正規化された球面調和関数はすべて、
Y
ℓ
m
∗
(
θ
,
φ
)
=
(
−
1
)
−
m
Y
ℓ
−
m
(
θ
,
φ
)
,
{\displaystyle Y_{\ell }^{m}{}^{*}(\theta ,\varphi )=(-1)^{-m}Y_{\ell }^{-m}(\theta ,\varphi ),}
ここで、上付き文字 * は複素共役を表します。あるいは、この式は球面調和関数と ウィグナー D 行列 の関係から導かれます 。
コンドン・ショートリー相
球面調和関数の定義に関する混乱の原因の 1 つは、量子力学の文献で一般に コンドン ・ショートリー位相と呼ばれる位相係数 に関する ものです。量子力学のコミュニティでは、この 位相係数を ルジャンドル随伴多項式 の定義に含める か、球面調和関数の定義に追加するのが一般的な方法です。球面調和関数の定義にコンドン・ショートリー位相を使用する必要はありませんが、これを含めると、特に 昇順演算子 と降順演算子 の適用など、一部の量子力学演算が簡素化されます。測地学 [12] と磁気学のコミュニティでは、球面調和関数の定義にもルジャンドル随伴多項式の定義にもコンドン・ショートリー位相係数を含めることはありません。 [13]
(
−
1
)
m
{\displaystyle (-1)^{m}}
球面調和関数の実数基底は、複素数類似物 を用いて
次のように定義できる。ここでは一貫性を保つためにコンドン・ショートリー位相規則が用いられる。 実球面調和関数を用いて
複素数球面調和関数を定義する対応する逆方程式は 、
Y
ℓ
m
:
S
2
→
R
{\displaystyle Y_{\ell m}:S^{2}\to \mathbb {R} }
Y
ℓ
m
:
S
2
→
C
{\displaystyle Y_{\ell }^{m}:S^{2}\to \mathbb {C} }
Y
ℓ
m
=
{
i
2
(
Y
ℓ
m
−
(
−
1
)
m
Y
ℓ
−
m
)
if
m
<
0
Y
ℓ
0
if
m
=
0
1
2
(
Y
ℓ
−
m
+
(
−
1
)
m
Y
ℓ
m
)
if
m
>
0.
=
{
i
2
(
Y
ℓ
−
|
m
|
−
(
−
1
)
m
Y
ℓ
|
m
|
)
if
m
<
0
Y
ℓ
0
if
m
=
0
1
2
(
Y
ℓ
−
|
m
|
+
(
−
1
)
m
Y
ℓ
|
m
|
)
if
m
>
0.
=
{
2
(
−
1
)
m
ℑ
[
Y
ℓ
|
m
|
]
if
m
<
0
Y
ℓ
0
if
m
=
0
2
(
−
1
)
m
ℜ
[
Y
ℓ
m
]
if
m
>
0.
{\displaystyle {\begin{aligned}Y_{\ell m}&={\begin{cases}{\dfrac {i}{\sqrt {2}}}\left(Y_{\ell }^{m}-(-1)^{m}\,Y_{\ell }^{-m}\right)&{\text{if}}\ m<0\\Y_{\ell }^{0}&{\text{if}}\ m=0\\{\dfrac {1}{\sqrt {2}}}\left(Y_{\ell }^{-m}+(-1)^{m}\,Y_{\ell }^{m}\right)&{\text{if}}\ m>0.\end{cases}}\\&={\begin{cases}{\dfrac {i}{\sqrt {2}}}\left(Y_{\ell }^{-|m|}-(-1)^{m}\,Y_{\ell }^{|m|}\right)&{\text{if}}\ m<0\\Y_{\ell }^{0}&{\text{if}}\ m=0\\{\dfrac {1}{\sqrt {2}}}\left(Y_{\ell }^{-|m|}+(-1)^{m}\,Y_{\ell }^{|m|}\right)&{\text{if}}\ m>0.\end{cases}}\\&={\begin{cases}{\sqrt {2}}\,(-1)^{m}\,\Im [{Y_{\ell }^{|m|}}]&{\text{if}}\ m<0\\Y_{\ell }^{0}&{\text{if}}\ m=0\\{\sqrt {2}}\,(-1)^{m}\,\Re [{Y_{\ell }^{m}}]&{\text{if}}\ m>0.\end{cases}}\end{aligned}}}
Y
ℓ
m
:
S
2
→
C
{\displaystyle Y_{\ell }^{m}:S^{2}\to \mathbb {C} }
Y
ℓ
m
:
S
2
→
R
{\displaystyle Y_{\ell m}:S^{2}\to \mathbb {R} }
Y
ℓ
m
=
{
1
2
(
Y
ℓ
|
m
|
−
i
Y
ℓ
,
−
|
m
|
)
if
m
<
0
Y
ℓ
0
if
m
=
0
(
−
1
)
m
2
(
Y
ℓ
|
m
|
+
i
Y
ℓ
,
−
|
m
|
)
if
m
>
0.
{\displaystyle Y_{\ell }^{m}={\begin{cases}{\dfrac {1}{\sqrt {2}}}\left(Y_{\ell |m|}-iY_{\ell ,-|m|}\right)&{\text{if}}\ m<0\\[4pt]Y_{\ell 0}&{\text{if}}\ m=0\\[4pt]{\dfrac {(-1)^{m}}{\sqrt {2}}}\left(Y_{\ell |m|}+iY_{\ell ,-|m|}\right)&{\text{if}}\ m>0.\end{cases}}}
実球面調和関数は、 テッセラル球面調和関数 と呼ばれることもあります 。 [14] これらの関数は、上記の複素関数と同じ正規直交性を持っています。m > 0 の 実球面調和関数は コサイン型、 m < 0の実球面調和関数はサイン型と呼ばれます。この理由 は 、関数をルジャンドル多項式で次のように書くことでわかります。
Y
ℓ
m
:
S
2
→
R
{\displaystyle Y_{\ell m}:S^{2}\to \mathbb {R} }
Y
ℓ
m
:
S
2
→
C
{\displaystyle Y_{\ell }^{m}:S^{2}\to \mathbb {C} }
Y
ℓ
m
{\displaystyle Y_{\ell m}}
Y
ℓ
m
=
{
(
−
1
)
m
2
2
ℓ
+
1
4
π
(
ℓ
−
|
m
|
)
!
(
ℓ
+
|
m
|
)
!
P
ℓ
|
m
|
(
cos
θ
)
sin
(
|
m
|
φ
)
if
m
<
0
2
ℓ
+
1
4
π
P
ℓ
m
(
cos
θ
)
if
m
=
0
(
−
1
)
m
2
2
ℓ
+
1
4
π
(
ℓ
−
m
)
!
(
ℓ
+
m
)
!
P
ℓ
m
(
cos
θ
)
cos
(
m
φ
)
if
m
>
0
.
{\displaystyle Y_{\ell m}={\begin{cases}\left(-1\right)^{m}{\sqrt {2}}{\sqrt {{\dfrac {2\ell +1}{4\pi }}{\dfrac {(\ell -|m|)!}{(\ell +|m|)!}}}}\;P_{\ell }^{|m|}(\cos \theta )\ \sin(|m|\varphi )&{\text{if }}m<0\\[4pt]{\sqrt {\dfrac {2\ell +1}{4\pi }}}\ P_{\ell }^{m}(\cos \theta )&{\text{if }}m=0\\[4pt]\left(-1\right)^{m}{\sqrt {2}}{\sqrt {{\dfrac {2\ell +1}{4\pi }}{\dfrac {(\ell -m)!}{(\ell +m)!}}}}\;P_{\ell }^{m}(\cos \theta )\ \cos(m\varphi )&{\text{if }}m>0\,.\end{cases}}}
同じ正弦係数と余弦係数は、直交座標表現を扱う次のサブセクションでも見られます。
までの実球面調和関数のリストについては、 ここを 参照してください 。これは、上記の方程式の出力と一致することがわかります。
ℓ
=
4
{\displaystyle \ell =4}
量子化学における利用
水素原子の解析解からわかるように、波動関数の角度部分の固有関数は球面調和関数です。ただし、磁気項のない非相対論的シュレーディンガー方程式の解は実数にすることができます。これが、量子化学の基底関数で実数形式が広く使用されている理由です。プログラムで複素代数を使用する必要がないためです。ここで、実数関数は複素関数と同じ空間をカバーします。
たとえば、球面調和関数の表 からわかるように 、通常の p 関数 ( ) は複素数で軸方向が混在していますが、 実数バージョンは基本的に x 、 y 、 z だけです 。
ℓ
=
1
{\displaystyle \ell =1}
複素球面調和関数は、 から の 全体に を 次数 の 同次関数 として 拡張することによって 、つまり と設定すること
によって 、立体調和関数 を生じます。 は、 次数 の調和多項式と 同次多項式 の空間の基底である
ことがわかります。 より具体的には、これは回転群のこの表現の (正規化まで一意の) ゲルファント-ツェトリン基底であり 、その事実から、直交座標における の 明示的な式 を導くことができます。
Y
ℓ
m
{\displaystyle Y_{\ell }^{m}}
S
2
{\displaystyle S^{2}}
R
3
{\displaystyle \mathbb {R} ^{3}}
ℓ
{\displaystyle \ell }
R
ℓ
m
(
v
)
:=
‖
v
‖
ℓ
Y
ℓ
m
(
v
‖
v
‖
)
{\displaystyle R_{\ell }^{m}(v):=\|v\|^{\ell }Y_{\ell }^{m}\left({\frac {v}{\|v\|}}\right)}
R
ℓ
m
{\displaystyle R_{\ell }^{m}}
ℓ
{\displaystyle \ell }
S
O
(
3
)
{\displaystyle SO(3)}
R
ℓ
m
{\displaystyle R_{\ell }^{m}}
ヘルグロッツ生成関数
量子力学の慣例が に採用された場合 、 と
なります。
ここで、 は 、 、 を成分とするベクトルであり 、 は
複素座標を持つベクトルです。
Y
ℓ
m
:
S
2
→
C
{\displaystyle Y_{\ell }^{m}:S^{2}\to \mathbb {C} }
e
v
a
⋅
r
=
∑
ℓ
=
0
∞
∑
m
=
−
ℓ
ℓ
4
π
2
ℓ
+
1
r
ℓ
v
ℓ
λ
m
(
ℓ
+
m
)
!
(
ℓ
−
m
)
!
Y
ℓ
m
(
r
/
r
)
.
{\displaystyle e^{v{\mathbf {a} }\cdot {\mathbf {r} }}=\sum _{\ell =0}^{\infty }\sum _{m=-\ell }^{\ell }{\sqrt {\frac {4\pi }{2\ell +1}}}{\frac {r^{\ell }v^{\ell }{\lambda ^{m}}}{\sqrt {(\ell +m)!(\ell -m)!}}}Y_{\ell }^{m}(\mathbf {r} /r).}
r
{\displaystyle \mathbf {r} }
(
x
,
y
,
z
)
∈
R
3
{\displaystyle (x,y,z)\in \mathbb {R} ^{3}}
r
=
|
r
|
{\displaystyle r=|\mathbf {r} |}
a
=
z
^
−
λ
2
(
x
^
+
i
y
^
)
+
1
2
λ
(
x
^
−
i
y
^
)
.
{\displaystyle {\mathbf {a} }={\mathbf {\hat {z}} }-{\frac {\lambda }{2}}\left({\mathbf {\hat {x}} }+i{\mathbf {\hat {y}} }\right)+{\frac {1}{2\lambda }}\left({\mathbf {\hat {x}} }-i{\mathbf {\hat {y}} }\right).}
a
{\displaystyle \mathbf {a} }
a
=
[
1
2
(
1
λ
−
λ
)
,
−
i
2
(
1
λ
+
λ
)
,
1
]
.
{\displaystyle \mathbf {a} =[{\frac {1}{2}}({\frac {1}{\lambda }}-\lambda ),-{\frac {i}{2}}({\frac {1}{\lambda }}+\lambda ),1].}
の本質的な性質は、 それが null であることです。
a
{\displaystyle \mathbf {a} }
a
⋅
a
=
0.
{\displaystyle \mathbf {a} \cdot \mathbf {a} =0.}
と を 実パラメータとして とれば十分です。この生成関数を Herglotz にちなんで命名したのは、Courant & Hilbert 1962, §VII.7 に従っており、この発見は Herglotz による未発表のノートによるものです。
v
{\displaystyle v}
λ
{\displaystyle \lambda }
本質的に、球面調和関数のすべての特性は、この生成関数から導くことができる。 [15] この定義の直接的な利点は、ベクトルを 量子力学スピンベクトル演算子 に置き換えて 、 固体調和関数 の演算子アナログにすると 、 [16]標準化された 球面テンソル演算子 のセットの生成関数が得られることで ある
。
r
{\displaystyle \mathbf {r} }
J
{\displaystyle \mathbf {J} }
Y
ℓ
m
(
J
)
{\displaystyle {\mathcal {Y}}_{\ell }^{m}({\mathbf {J} })}
r
ℓ
Y
ℓ
m
(
r
/
r
)
{\displaystyle r^{\ell }Y_{\ell }^{m}(\mathbf {r} /r)}
Y
ℓ
m
(
J
)
{\displaystyle {\mathcal {Y}}_{\ell }^{m}({\mathbf {J} })}
e
v
a
⋅
J
=
∑
ℓ
=
0
∞
∑
m
=
−
ℓ
ℓ
4
π
2
ℓ
+
1
v
ℓ
λ
m
(
ℓ
+
m
)
!
(
ℓ
−
m
)
!
Y
ℓ
m
(
J
)
.
{\displaystyle e^{v{\mathbf {a} }\cdot {\mathbf {J} }}=\sum _{\ell =0}^{\infty }\sum _{m=-\ell }^{\ell }{\sqrt {\frac {4\pi }{2\ell +1}}}{\frac {v^{\ell }{\lambda ^{m}}}{\sqrt {(\ell +m)!(\ell -m)!}}}{\mathcal {Y}}_{\ell }^{m}({\mathbf {J} }).}
2 つの定義の平行性により、 は と同様に回転 (下記参照) によって変換されることが 保証され、これにより、 が球面テンソル演算子 であり 、 および となり、 クレプシュ-ゴルダン 合成定理や ウィグナー-エッカート定理 などのそのような演算子のすべての特性に従うことが保証されます 。さらに、これらは固定スケールまたは正規化による標準化されたセットです。
Y
ℓ
m
{\displaystyle {\mathcal {Y}}_{\ell }^{m}}
Y
ℓ
m
{\displaystyle Y_{\ell }^{m}}
T
q
(
k
)
{\displaystyle T_{q}^{(k)}}
k
=
ℓ
{\displaystyle k={\ell }}
q
=
m
{\displaystyle q=m}
ヘルグロッツの定義は、さらに因数分解して の多項式 と および の多項式にすることができる。 その場合 、次のようになる(コンドン・ショートリー位相)。
および m = 0
の場合 :
ここで
および
の場合、 これは次のように帰着する。
z
{\displaystyle z}
x
{\displaystyle x}
y
{\displaystyle y}
r
ℓ
(
Y
ℓ
m
Y
ℓ
−
m
)
=
[
2
ℓ
+
1
4
π
]
1
/
2
Π
¯
ℓ
m
(
z
)
(
(
−
1
)
m
(
A
m
+
i
B
m
)
(
A
m
−
i
B
m
)
)
,
m
>
0.
{\displaystyle r^{\ell }\,{\begin{pmatrix}Y_{\ell }^{m}\\Y_{\ell }^{-m}\end{pmatrix}}=\left[{\frac {2\ell +1}{4\pi }}\right]^{1/2}{\bar {\Pi }}_{\ell }^{m}(z){\begin{pmatrix}\left(-1\right)^{m}(A_{m}+iB_{m})\\(A_{m}-iB_{m})\end{pmatrix}},\qquad m>0.}
r
ℓ
Y
ℓ
0
≡
2
ℓ
+
1
4
π
Π
¯
ℓ
0
.
{\displaystyle r^{\ell }\,Y_{\ell }^{0}\equiv {\sqrt {\frac {2\ell +1}{4\pi }}}{\bar {\Pi }}_{\ell }^{0}.}
A
m
(
x
,
y
)
=
∑
p
=
0
m
(
m
p
)
x
p
y
m
−
p
cos
(
(
m
−
p
)
π
2
)
,
{\displaystyle A_{m}(x,y)=\sum _{p=0}^{m}{\binom {m}{p}}x^{p}y^{m-p}\cos \left((m-p){\frac {\pi }{2}}\right),}
B
m
(
x
,
y
)
=
∑
p
=
0
m
(
m
p
)
x
p
y
m
−
p
sin
(
(
m
−
p
)
π
2
)
,
{\displaystyle B_{m}(x,y)=\sum _{p=0}^{m}{\binom {m}{p}}x^{p}y^{m-p}\sin \left((m-p){\frac {\pi }{2}}\right),}
Π
¯
ℓ
m
(
z
)
=
[
(
ℓ
−
m
)
!
(
ℓ
+
m
)
!
]
1
/
2
∑
k
=
0
⌊
(
ℓ
−
m
)
/
2
⌋
(
−
1
)
k
2
−
ℓ
(
ℓ
k
)
(
2
ℓ
−
2
k
ℓ
)
(
ℓ
−
2
k
)
!
(
ℓ
−
2
k
−
m
)
!
r
2
k
z
ℓ
−
2
k
−
m
.
{\displaystyle {\bar {\Pi }}_{\ell }^{m}(z)=\left[{\frac {(\ell -m)!}{(\ell +m)!}}\right]^{1/2}\sum _{k=0}^{\left\lfloor (\ell -m)/2\right\rfloor }(-1)^{k}2^{-\ell }{\binom {\ell }{k}}{\binom {2\ell -2k}{\ell }}{\frac {(\ell -2k)!}{(\ell -2k-m)!}}\;r^{2k}\;z^{\ell -2k-m}.}
m
=
0
{\displaystyle m=0}
Π
¯
ℓ
0
(
z
)
=
∑
k
=
0
⌊
ℓ
/
2
⌋
(
−
1
)
k
2
−
ℓ
(
ℓ
k
)
(
2
ℓ
−
2
k
ℓ
)
r
2
k
z
ℓ
−
2
k
.
{\displaystyle {\bar {\Pi }}_{\ell }^{0}(z)=\sum _{k=0}^{\left\lfloor \ell /2\right\rfloor }(-1)^{k}2^{-\ell }{\binom {\ell }{k}}{\binom {2\ell -2k}{\ell }}\;r^{2k}\;z^{\ell -2k}.}
因数は 本質的にルジャンドルの随伴多項式 であり 、因数 は本質的に です 。
Π
¯
ℓ
m
(
z
)
{\displaystyle {\bar {\Pi }}_{\ell }^{m}(z)}
P
ℓ
m
(
cos
θ
)
{\displaystyle P_{\ell }^{m}(\cos \theta )}
(
A
m
±
i
B
m
)
{\displaystyle (A_{m}\pm iB_{m})}
e
±
i
m
φ
{\displaystyle e^{\pm im\varphi }}
例
上記に明示的に記載した、、 の 式を使用すると、 次のようになります。
Π
¯
ℓ
m
(
z
)
{\displaystyle {\bar {\Pi }}_{\ell }^{m}(z)}
A
m
(
x
,
y
)
{\displaystyle A_{m}(x,y)}
B
m
(
x
,
y
)
{\displaystyle B_{m}(x,y)}
Y
3
1
=
−
1
r
3
[
7
4
π
⋅
3
16
]
1
/
2
(
5
z
2
−
r
2
)
(
x
+
i
y
)
=
−
[
7
4
π
⋅
3
16
]
1
/
2
(
5
cos
2
θ
−
1
)
(
sin
θ
e
i
φ
)
{\displaystyle Y_{3}^{1}=-{\frac {1}{r^{3}}}\left[{\tfrac {7}{4\pi }}\cdot {\tfrac {3}{16}}\right]^{1/2}\left(5z^{2}-r^{2}\right)\left(x+iy\right)=-\left[{\tfrac {7}{4\pi }}\cdot {\tfrac {3}{16}}\right]^{1/2}\left(5\cos ^{2}\theta -1\right)\left(\sin \theta e^{i\varphi }\right)}
Y
4
−
2
=
1
r
4
[
9
4
π
⋅
5
32
]
1
/
2
(
7
z
2
−
r
2
)
(
x
−
i
y
)
2
=
[
9
4
π
⋅
5
32
]
1
/
2
(
7
cos
2
θ
−
1
)
(
sin
2
θ
e
−
2
i
φ
)
{\displaystyle Y_{4}^{-2}={\frac {1}{r^{4}}}\left[{\tfrac {9}{4\pi }}\cdot {\tfrac {5}{32}}\right]^{1/2}\left(7z^{2}-r^{2}\right)\left(x-iy\right)^{2}=\left[{\tfrac {9}{4\pi }}\cdot {\tfrac {5}{32}}\right]^{1/2}\left(7\cos ^{2}\theta -1\right)\left(\sin ^{2}\theta e^{-2i\varphi }\right)}
これは、ここ と ここ
に記載されている関数と一致することが確認できます 。
上記の式を使用して実球面調和関数を形成すると、項(コサイン)のみが 含まれ、項(サイン) のみが含まれる ことがわかります 。
m
>
0
{\displaystyle m>0}
A
m
{\displaystyle A_{m}}
m
<
0
{\displaystyle m<0}
B
m
{\displaystyle B_{m}}
r
ℓ
(
Y
ℓ
m
Y
ℓ
−
m
)
=
2
ℓ
+
1
2
π
Π
¯
ℓ
m
(
z
)
(
A
m
B
m
)
,
m
>
0.
{\displaystyle r^{\ell }\,{\begin{pmatrix}Y_{\ell m}\\Y_{\ell -m}\end{pmatrix}}={\sqrt {\frac {2\ell +1}{2\pi }}}{\bar {\Pi }}_{\ell }^{m}(z){\begin{pmatrix}A_{m}\\B_{m}\end{pmatrix}},\qquad m>0.}
m = 0
の場合 :
r
ℓ
Y
ℓ
0
≡
2
ℓ
+
1
4
π
Π
¯
ℓ
0
.
{\displaystyle r^{\ell }\,Y_{\ell 0}\equiv {\sqrt {\frac {2\ell +1}{4\pi }}}{\bar {\Pi }}_{\ell }^{0}.}
特殊なケースと値
のとき 、球面調和関数は通常の ルジャンドル多項式 に簡約されます 。
m
=
0
{\displaystyle m=0}
Y
ℓ
m
:
S
2
→
C
{\displaystyle Y_{\ell }^{m}:S^{2}\to \mathbb {C} }
Y
ℓ
0
(
θ
,
φ
)
=
2
ℓ
+
1
4
π
P
ℓ
(
cos
θ
)
.
{\displaystyle Y_{\ell }^{0}(\theta ,\varphi )={\sqrt {\frac {2\ell +1}{4\pi }}}P_{\ell }(\cos \theta ).}
、 またはもっと簡単に直交座標で 、
m
=
±
ℓ
{\displaystyle m=\pm \ell }
Y
ℓ
±
ℓ
(
θ
,
φ
)
=
(
∓
1
)
ℓ
2
ℓ
ℓ
!
(
2
ℓ
+
1
)
!
4
π
sin
ℓ
θ
e
±
i
ℓ
φ
,
{\displaystyle Y_{\ell }^{\pm \ell }(\theta ,\varphi )={\frac {(\mp 1)^{\ell }}{2^{\ell }\ell !}}{\sqrt {\frac {(2\ell +1)!}{4\pi }}}\sin ^{\ell }\theta \,e^{\pm i\ell \varphi },}
r
ℓ
Y
ℓ
±
ℓ
(
r
)
=
(
∓
1
)
ℓ
2
ℓ
ℓ
!
(
2
ℓ
+
1
)
!
4
π
(
x
±
i
y
)
ℓ
.
{\displaystyle r^{\ell }Y_{\ell }^{\pm \ell }({\mathbf {r} })={\frac {(\mp 1)^{\ell }}{2^{\ell }\ell !}}{\sqrt {\frac {(2\ell +1)!}{4\pi }}}(x\pm iy)^{\ell }.}
北極では、、 および が定義されていない場合、 を除くすべての球面調和関数は 消えます。
θ
=
0
{\displaystyle \theta =0}
φ
{\displaystyle \varphi }
m
=
0
{\displaystyle m=0}
Y
ℓ
m
(
0
,
φ
)
=
Y
ℓ
m
(
z
)
=
2
ℓ
+
1
4
π
δ
m
0
.
{\displaystyle Y_{\ell }^{m}(0,\varphi )=Y_{\ell }^{m}({\mathbf {z} })={\sqrt {\frac {2\ell +1}{4\pi }}}\delta _{m0}.}
対称性の特性
球面調和関数は、空間反転(パリティ)と回転の操作の下で、深く重要な特性を持ちます。
パリティ
球面調和関数は明確な偶奇性を持つ。つまり、原点の周りの反転に関して偶奇である。反転は演算子 で表される 。すると、さまざまな方法でわかるように(おそらく最も単純なのはヘルグロッツ生成関数から)、 は 単位ベクトルであり、
P
Ψ
(
r
)
=
Ψ
(
−
r
)
{\displaystyle P\Psi (\mathbf {r} )=\Psi (-\mathbf {r} )}
r
{\displaystyle \mathbf {r} }
Y
ℓ
m
(
−
r
)
=
(
−
1
)
ℓ
Y
ℓ
m
(
r
)
.
{\displaystyle Y_{\ell }^{m}(-\mathbf {r} )=(-1)^{\ell }Y_{\ell }^{m}(\mathbf {r} ).}
球面角に関して、パリティは座標を持つ点を に変換します 。 球面調和関数のパリティのステートメントは、次のようになります。
(これは次のように見ることができます。 ルジャンドル多項式は (−1) ℓ + m を与え 、指数関数から (−1) m を 与え、これらを合わせると球面調和関数のパリティは (−1) ℓ になります。)
{
θ
,
φ
}
{\displaystyle \{\theta ,\varphi \}}
{
π
−
θ
,
π
+
φ
}
{\displaystyle \{\pi -\theta ,\pi +\varphi \}}
Y
ℓ
m
(
θ
,
φ
)
→
Y
ℓ
m
(
π
−
θ
,
π
+
φ
)
=
(
−
1
)
ℓ
Y
ℓ
m
(
θ
,
φ
)
{\displaystyle Y_{\ell }^{m}(\theta ,\varphi )\to Y_{\ell }^{m}(\pi -\theta ,\pi +\varphi )=(-1)^{\ell }Y_{\ell }^{m}(\theta ,\varphi )}
パリティは実球面調和関数および高次元の球面調和関数に対しても成立します。つまり、 ℓ 次球面調和関数に 点反射を 適用すると、符号が (−1) ℓ 倍に変化します。
回転
m = 0 、 ℓ = 3 の実球関数の回転 。実関数が示されているため係数はウィグナーD行列と等しくないが、複素関数を再分解することで得られる。
単位ベクトルを に 送る原点周りの 回転を考えてみましょう 。この操作では、次数 と順序の球面調和関数は 、同じ次数の球面調和関数の線形結合に変換されます。つまり、
は 回転 に依存する 順序の行列です 。ただし、これはこの特性を表現する標準的な方法ではありません。標準的な方法では、次のように書きます。
R
{\displaystyle {\mathcal {R}}}
r
{\displaystyle \mathbf {r} }
r
′
{\displaystyle \mathbf {r} '}
ℓ
{\displaystyle \ell }
m
{\displaystyle m}
Y
ℓ
m
(
r
′
)
=
∑
m
′
=
−
ℓ
ℓ
A
m
m
′
Y
ℓ
m
′
(
r
)
,
{\displaystyle Y_{\ell }^{m}({\mathbf {r} }')=\sum _{m'=-\ell }^{\ell }A_{mm'}Y_{\ell }^{m'}({\mathbf {r} }),}
A
m
m
′
{\displaystyle A_{mm'}}
(
2
ℓ
+
1
)
{\displaystyle (2\ell +1)}
R
{\displaystyle {\mathcal {R}}}
Y
ℓ
m
(
r
′
)
=
∑
m
′
=
−
ℓ
ℓ
[
D
m
m
′
(
ℓ
)
(
R
)
]
∗
Y
ℓ
m
′
(
r
)
,
{\displaystyle Y_{\ell }^{m}({\mathbf {r} }')=\sum _{m'=-\ell }^{\ell }[D_{mm'}^{(\ell )}({\mathcal {R}})]^{*}Y_{\ell }^{m'}({\mathbf {r} }),}
ここで、は ウィグナーD行列 の要素の複素共役です 。特に、が 方位角の回転
であるとき、次の恒等式が得られます。
D
m
m
′
(
ℓ
)
(
R
)
∗
{\displaystyle D_{mm'}^{(\ell )}({\mathcal {R}})^{*}}
r
′
{\displaystyle \mathbf {r} '}
ϕ
0
{\displaystyle \phi _{0}}
Y
ℓ
m
(
r
′
)
=
Y
ℓ
m
(
r
)
e
i
m
ϕ
0
.
{\displaystyle Y_{\ell }^{m}({\mathbf {r} }')=Y_{\ell }^{m}({\mathbf {r} })e^{im\phi _{0}}.}
球面調和関数の回転挙動は、群論の観点からは、おそらくその典型的な特徴です。 次数の は、 次元 の群 SO(3) の既約表現の関数の基底セットを提供します 。解析法を使用して苦労して証明された球面調和関数に関する多くの事実 (加法定理など) は、対称性の方法を使用すると、より簡単な証明とより深い意味を獲得します。
Y
ℓ
m
{\displaystyle Y_{\ell }^{m}}
ℓ
{\displaystyle \ell }
(
2
ℓ
+
1
)
{\displaystyle (2\ell +1)}
球面調和関数展開
ラプラス球面調和関数は 正規直交関数の完全な集合を形成し、したがって 二乗可積分関数 の ヒルベルト空間 の 正規直交基底 を形成します。したがって、単位球面上では 、任意の二乗可積分関数は これらの線形結合として展開できます。
Y
ℓ
m
:
S
2
→
C
{\displaystyle Y_{\ell }^{m}:S^{2}\to \mathbb {C} }
L
C
2
(
S
2
)
{\displaystyle L_{\mathbb {C} }^{2}(S^{2})}
S
2
{\displaystyle S^{2}}
f
:
S
2
→
C
{\displaystyle f:S^{2}\to \mathbb {C} }
f
(
θ
,
φ
)
=
∑
ℓ
=
0
∞
∑
m
=
−
ℓ
ℓ
f
ℓ
m
Y
ℓ
m
(
θ
,
φ
)
.
{\displaystyle f(\theta ,\varphi )=\sum _{\ell =0}^{\infty }\sum _{m=-\ell }^{\ell }f_{\ell }^{m}\,Y_{\ell }^{m}(\theta ,\varphi ).}
この展開は平均二乗収束、つまり球面の L 2 における収束の意味で成り立ち、つまり
lim
N
→
∞
∫
0
2
π
∫
0
π
|
f
(
θ
,
φ
)
−
∑
ℓ
=
0
N
∑
m
=
−
ℓ
ℓ
f
ℓ
m
Y
ℓ
m
(
θ
,
φ
)
|
2
sin
θ
d
θ
d
φ
=
0.
{\displaystyle \lim _{N\to \infty }\int _{0}^{2\pi }\int _{0}^{\pi }\left|f(\theta ,\varphi )-\sum _{\ell =0}^{N}\sum _{m=-\ell }^{\ell }f_{\ell }^{m}Y_{\ell }^{m}(\theta ,\varphi )\right|^{2}\sin \theta \,d\theta \,d\varphi =0.}
展開係数はフーリエ係数 の類似物であり 、上記の式に球面調和関数の複素共役を掛け、立体角 Ω にわたって積分し、上記の直交関係を利用することで得られます。これは、基本的なヒルベルト空間理論によって厳密に正当化されます。正規直交化調和関数の場合、次の式が得られます。
f
ℓ
m
=
∫
Ω
f
(
θ
,
φ
)
Y
ℓ
m
∗
(
θ
,
φ
)
d
Ω
=
∫
0
2
π
d
φ
∫
0
π
d
θ
sin
θ
f
(
θ
,
φ
)
Y
ℓ
m
∗
(
θ
,
φ
)
.
{\displaystyle f_{\ell }^{m}=\int _{\Omega }f(\theta ,\varphi )\,Y_{\ell }^{m*}(\theta ,\varphi )\,d\Omega =\int _{0}^{2\pi }d\varphi \int _{0}^{\pi }\,d\theta \,\sin \theta f(\theta ,\varphi )Y_{\ell }^{m*}(\theta ,\varphi ).}
ℓ の係数が 十分急速に(たとえば 指数関数 的に)減少する場合、級数も f に 一様に収束します 。
二乗可積分関数は、 上記の
実数調和関数の和として展開することもできる。
f
:
S
2
→
R
{\displaystyle f:S^{2}\to \mathbb {R} }
Y
ℓ
m
:
S
2
→
R
{\displaystyle Y_{\ell m}:S^{2}\to \mathbb {R} }
f
(
θ
,
φ
)
=
∑
ℓ
=
0
∞
∑
m
=
−
ℓ
ℓ
f
ℓ
m
Y
ℓ
m
(
θ
,
φ
)
.
{\displaystyle f(\theta ,\varphi )=\sum _{\ell =0}^{\infty }\sum _{m=-\ell }^{\ell }f_{\ell m}\,Y_{\ell m}(\theta ,\varphi ).}
級数の収束は同じ意味で再び成り立ち、つまり、実球面調和関数は 正規直交関数の完全なセットを形成し、したがって 二乗可積分関数 の ヒルベルト空間 の 正規直交基底 を形成します。実調和関数による展開の利点は 、実関数の場合、 展開係数が 実数であることが保証されるのに対し、 (関数として考える) による展開の係数は その特性を持たないことです。
Y
ℓ
m
:
S
2
→
R
{\displaystyle Y_{\ell m}:S^{2}\to \mathbb {R} }
L
R
2
(
S
2
)
{\displaystyle L_{\mathbb {R} }^{2}(S^{2})}
Y
ℓ
m
{\displaystyle Y_{\ell m}}
f
:
S
2
→
R
{\displaystyle f:S^{2}\to \mathbb {R} }
f
ℓ
m
{\displaystyle f_{\ell m}}
f
ℓ
m
{\displaystyle f_{\ell }^{m}}
Y
ℓ
m
{\displaystyle Y_{\ell }^{m}}
f
:
S
2
→
C
⊃
R
{\displaystyle f:S^{2}\to \mathbb {C} \supset \mathbb {R} }
スペクトル分析
信号処理におけるパワースペクトル
関数 fの全パワーは、 信号処理の 文献では、 関数の積分の 2 乗をその定義域の面積で割ったものとして定義されています。 実数単位パワー球面調和関数の 正規直交性特性を使用すると、単位球上で定義された関数の全パワーが、 パーセバルの定理 の一般化によってそのスペクトル係数と関連していることを簡単に確認できます(ここで、定理はシュミット半正規化調和関数について述べられており、正規直交調和関数では関係が若干異なります)。
1
4
π
∫
Ω
|
f
(
Ω
)
|
2
d
Ω
=
∑
ℓ
=
0
∞
S
f
f
(
ℓ
)
,
{\displaystyle {\frac {1}{4\,\pi }}\int _{\Omega }|f(\Omega )|^{2}\,d\Omega =\sum _{\ell =0}^{\infty }S_{f\!f}(\ell ),}
どこ
S
f
f
(
ℓ
)
=
1
2
ℓ
+
1
∑
m
=
−
ℓ
ℓ
|
f
ℓ
m
|
2
{\displaystyle S_{f\!f}(\ell )={\frac {1}{2\ell +1}}\sum _{m=-\ell }^{\ell }|f_{\ell m}|^{2}}
は、角度パワースペクトル(シュミットの半正規化高調波の場合)として定義されます。同様に、2つの関数のクロスパワーを
次の
ように定義できます。
1
4
π
∫
Ω
f
(
Ω
)
g
∗
(
Ω
)
d
Ω
=
∑
ℓ
=
0
∞
S
f
g
(
ℓ
)
,
{\displaystyle {\frac {1}{4\,\pi }}\int _{\Omega }f(\Omega )\,g^{\ast }(\Omega )\,d\Omega =\sum _{\ell =0}^{\infty }S_{fg}(\ell ),}
S
f
g
(
ℓ
)
=
1
2
ℓ
+
1
∑
m
=
−
ℓ
ℓ
f
ℓ
m
g
ℓ
m
∗
{\displaystyle S_{fg}(\ell )={\frac {1}{2\ell +1}}\sum _{m=-\ell }^{\ell }f_{\ell m}g_{\ell m}^{\ast }}
はクロスパワースペクトルとして定義されます。関数 f と g の 平均がゼロの場合(つまり、スペクトル係数 f 00 と g 00 がゼロの場合)、 S ff ( ℓ ) と S fg ( ℓ )は、それぞれ次数 ℓ の関数の分散と共分散への寄与を表します 。(クロス)パワースペクトルは、次の形式のべき乗法則でよく近似されます。
S
f
f
(
ℓ
)
=
C
ℓ
β
.
{\displaystyle S_{f\!f}(\ell )=C\,\ell ^{\beta }.}
β = 0 の場合 、各次数が等しいパワーを持つため、スペクトルは「白」になります。 β < 0 の場合、長波長の低次数の方が高次数よりもパワーが大きいため、スペクトルは「赤」と呼ばれます。 最後に、 β > 0の場合、スペクトルは「青」と呼ばれます。 S ff ( ℓ ) の成長順序の条件は、次のセクションの f の微分可能性の順序に関連しています 。
微分可能性特性
元の関数 f の微分可能性の特性は、 S ff ( ℓ ) の漸近 特性 の観点からも 理解できます 。特に、 ℓ → ∞ のときに S ff ( ℓ ) が ℓ の任意の 有理関数 よりも速く減少する場合 、 f は 無限に微分可能 です 。さらに、 S ff ( ℓ ) が 指数的に減少する場合、 f は 球面上で
実際に 実解析的 です。
一般的な手法は、ソボレフ空間 の理論を使用することです。 S ff ( ℓ ) の成長と微分可能性の関係は、 フーリエ級数 の係数の成長に関する類似の結果に似ています 。具体的には、
の場合、 f は ソボレフ空間 H s ( S 2 ) 内にあります。特に、 ソボレフの埋め込み定理は、
すべての sに対して f が無限に微分可能である
ことを意味します 。
∑
ℓ
=
0
∞
(
1
+
ℓ
2
)
s
S
f
f
(
ℓ
)
<
∞
,
{\displaystyle \sum _{\ell =0}^{\infty }(1+\ell ^{2})^{s}S_{ff}(\ell )<\infty ,}
S
f
f
(
ℓ
)
=
O
(
ℓ
−
s
)
a
s
ℓ
→
∞
{\displaystyle S_{ff}(\ell )=O(\ell ^{-s})\quad {\rm {{as\ }\ell \to \infty }}}
代数的性質
加法定理
非常に興味深く、利用されている数学的結果は、球面調和関数の 加法定理 と呼ばれています。 球面座標がそれぞれ と である 2 つのベクトル r と r′ が 与えられた場合、 それらの間の
角度
は、右側に現れる三角関数の役割が球面調和関数によって果たされ、左側の役割がルジャンドル 多項式 によって果たされる関係によって与えられます。
(
r
,
θ
,
φ
)
{\displaystyle (r,\theta ,\varphi )}
(
r
,
θ
′
,
φ
′
)
{\displaystyle (r,\theta ',\varphi ')}
γ
{\displaystyle \gamma }
cos
γ
=
cos
θ
′
cos
θ
+
sin
θ
sin
θ
′
cos
(
φ
−
φ
′
)
{\displaystyle \cos \gamma =\cos \theta '\cos \theta +\sin \theta \sin \theta '\cos(\varphi -\varphi ')}
加法 定理 によれば [17]
ここで、 P ℓ は ℓ 次 ルジャンドル多項式 である 。この表現は実高調波と複素高調波の両方に有効である。 [18]この結果は、単位球の ポアソン核 の特性を使用して解析的に証明する ことも、ベクトル yを z 軸に沿うように回転させて 右辺を直接計算することで幾何学的に証明することもできる。 [19]
特に、 x = y のとき、これはUnsöldの定理 [20]を与え、これは
cos 2 θ + sin 2 θ = 1
の恒等式 を2次元に一般化する。
∑
m
=
−
ℓ
ℓ
Y
ℓ
m
∗
(
x
)
Y
ℓ
m
(
x
)
=
2
ℓ
+
1
4
π
{\displaystyle \sum _{m=-\ell }^{\ell }Y_{\ell m}^{*}(\mathbf {x} )\,Y_{\ell m}(\mathbf {x} )={\frac {2\ell +1}{4\pi }}}
展開( 1 )において、左辺は ℓ 次ゾーン球面調和関数 の定数倍である 。この観点から、高次元への次の一般化が得られる。Y j を n 球面上のℓ 次 球面調和関数の 空間 H ℓ の任意の正規直交基底と する。すると、 単位ベクトル xに対応する ℓ 次ゾーン調和関数は 、次のように分解される [21]。
P
ℓ
(
x
⋅
y
)
{\displaystyle P_{\ell }(\mathbf {x} \cdot \mathbf {y} )}
Z
x
(
ℓ
)
{\displaystyle Z_{\mathbf {x} }^{(\ell )}}
さらに、ゾーン調和関数は適切な ゲーゲンバウアー多項式 の定数倍として与えられます 。
Z
x
(
ℓ
)
(
y
)
{\displaystyle Z_{\mathbf {x} }^{(\ell )}({\mathbf {y} })}
( 2 )と( 3 )を組み合わせると、 x と y が球座標で表されているとき、 n = 2 次元で ( 1 )が得られる。最後に、 x = y で評価すると関数の恒等式が得られる
。
ここで、 ωn − 1 は( n −1)球
の体積である。
dim
H
ℓ
ω
n
−
1
=
∑
j
=
1
dim
(
H
ℓ
)
|
Y
j
(
x
)
|
2
{\displaystyle {\frac {\dim \mathbf {H} _{\ell }}{\omega _{n-1}}}=\sum _{j=1}^{\dim(\mathbf {H} _{\ell })}|Y_{j}({\mathbf {x} })|^{2}}
短縮ルール
もう一つの有用な恒等式は、2つの球面調和関数の積を球面調和関数の和として表すものである [22]。
この和の項の多くは自明にゼロである。この和の項がゼロでない値となるとの
値は、 3j-記号 の選択規則によって決定される 。
Y
a
,
α
(
θ
,
φ
)
Y
b
,
β
(
θ
,
φ
)
=
(
2
a
+
1
)
(
2
b
+
1
)
4
π
∑
c
=
0
∞
∑
γ
=
−
c
c
(
−
1
)
γ
2
c
+
1
(
a
b
c
α
β
−
γ
)
(
a
b
c
0
0
0
)
Y
c
,
γ
(
θ
,
φ
)
.
{\displaystyle Y_{a,\alpha }\left(\theta ,\varphi \right)Y_{b,\beta }\left(\theta ,\varphi \right)={\sqrt {\frac {\left(2a+1\right)\left(2b+1\right)}{4\pi }}}\sum _{c=0}^{\infty }\sum _{\gamma =-c}^{c}\left(-1\right)^{\gamma }{\sqrt {2c+1}}{\begin{pmatrix}a&b&c\\\alpha &\beta &-\gamma \end{pmatrix}}{\begin{pmatrix}a&b&c\\0&0&0\end{pmatrix}}Y_{c,\gamma }\left(\theta ,\varphi \right).}
c
{\displaystyle c}
γ
{\displaystyle \gamma }
クレプシュ・ゴルダン係数
クレプシュ・ゴルダン係数は、2 つの球面調和関数の積を球面調和関数自体で展開したときに現れる係数です。本質的に同じ計算を行うために、ウィグナーの 3-jm 記号 、 ラカー係数 、 スレーター積分 など、さまざまな手法が利用できます。抽象的には、クレプシュ・ゴルダン係数は、 回転群 の 2 つの 既約表現の テンソル積 を既約表現の和として表します。適切に正規化されると、係数は重複度になります。
球面調和関数の可視化
単位球面とその節線上の の概略図。 は 、極を通る m 個の 大円と、等しい緯度の ℓ − m 個の円に沿って 0 に等しくなります 。関数は、これらの線を横切るたびに符号が変わります。
Y
ℓ
m
{\displaystyle Y_{\ell m}}
ℜ
[
Y
ℓ
m
]
{\displaystyle \Re [Y_{\ell m}]}
n = 5 次の球面調和関数の 3D カラー プロット 。n = ℓ であることに注意してください 。
ラプラス球面調和関数は、 その「 節線 」、つまり、または となる 球面上の点の集合を考えることで視覚化できます 。 の節線は ℓ 円で構成されます。経度に沿って | m | 円、緯度に沿って ℓ −| m | 円があります。 および方向 の のゼロの数を それぞれ数えることで、各タイプの節線の数を判断できます。 を の関数として考えると 、ルジャンドル陪多項式の実数部と虚数部はそれぞれ ℓ −| m | ゼロを持ち、それぞれが節の「緯度線」を生み出します。 一方、 を の関数として考えると 、三角関数の sin 関数と cos 関数は 2| m | ゼロを持ち、それぞれが節の「経度線」を生み出します。
Y
ℓ
m
{\displaystyle Y_{\ell }^{m}}
ℜ
[
Y
ℓ
m
]
=
0
{\displaystyle \Re [Y_{\ell }^{m}]=0}
ℑ
[
Y
ℓ
m
]
=
0
{\displaystyle \Im [Y_{\ell }^{m}]=0}
Y
ℓ
m
{\displaystyle Y_{\ell }^{m}}
Y
ℓ
m
{\displaystyle Y_{\ell }^{m}}
θ
{\displaystyle \theta }
φ
{\displaystyle \varphi }
Y
ℓ
m
{\displaystyle Y_{\ell }^{m}}
θ
{\displaystyle \theta }
Y
ℓ
m
{\displaystyle Y_{\ell }^{m}}
φ
{\displaystyle \varphi }
球面調和関数の次数 m が 0 の場合 (図の左上)、球面調和関数は経度に依存せず、帯状 と呼ばれます 。 この ような球面調和関数は、帯状球面関数 の特殊なケースです 。ℓ = | m | の場合(図の右下)、緯度にゼロ交差はなく、関数は 扇状 と呼ばれます。その他の場合、関数は 球面を チェックし、 テッセラル と呼ばれます。
より一般的な ℓ 次の球面調和関数は必ずしもラプラス基底の球面調和関数とは限らず 、その節点集合はかなり一般的な種類のものになり得る。 [23]
Y
ℓ
m
{\displaystyle Y_{\ell }^{m}}
球面調和関数のリスト
コンドン・ショートリー位相規則を使用する
最初のいくつかの正規直交化ラプラス球面調和関数の解析表現:
Y
ℓ
m
:
S
2
→
C
{\displaystyle Y_{\ell }^{m}:S^{2}\to \mathbb {C} }
Y
0
0
(
θ
,
φ
)
=
1
2
1
π
{\displaystyle Y_{0}^{0}(\theta ,\varphi )={\frac {1}{2}}{\sqrt {\frac {1}{\pi }}}}
Y
1
−
1
(
θ
,
φ
)
=
1
2
3
2
π
sin
θ
e
−
i
φ
Y
1
0
(
θ
,
φ
)
=
1
2
3
π
cos
θ
Y
1
1
(
θ
,
φ
)
=
−
1
2
3
2
π
sin
θ
e
i
φ
{\displaystyle {\begin{aligned}Y_{1}^{-1}(\theta ,\varphi )&={\frac {1}{2}}{\sqrt {\frac {3}{2\pi }}}\,\sin \theta \,e^{-i\varphi }\\Y_{1}^{0}(\theta ,\varphi )&={\frac {1}{2}}{\sqrt {\frac {3}{\pi }}}\,\cos \theta \\Y_{1}^{1}(\theta ,\varphi )&={\frac {-1}{2}}{\sqrt {\frac {3}{2\pi }}}\,\sin \theta \,e^{i\varphi }\end{aligned}}}
Y
2
−
2
(
θ
,
φ
)
=
1
4
15
2
π
sin
2
θ
e
−
2
i
φ
Y
2
−
1
(
θ
,
φ
)
=
1
2
15
2
π
sin
θ
cos
θ
e
−
i
φ
Y
2
0
(
θ
,
φ
)
=
1
4
5
π
(
3
cos
2
θ
−
1
)
Y
2
1
(
θ
,
φ
)
=
−
1
2
15
2
π
sin
θ
cos
θ
e
i
φ
Y
2
2
(
θ
,
φ
)
=
1
4
15
2
π
sin
2
θ
e
2
i
φ
{\displaystyle {\begin{aligned}Y_{2}^{-2}(\theta ,\varphi )&={\frac {1}{4}}{\sqrt {\frac {15}{2\pi }}}\,\sin ^{2}\theta \,e^{-2i\varphi }\\Y_{2}^{-1}(\theta ,\varphi )&={\frac {1}{2}}{\sqrt {\frac {15}{2\pi }}}\,\sin \theta \,\cos \theta \,e^{-i\varphi }\\Y_{2}^{0}(\theta ,\varphi )&={\frac {1}{4}}{\sqrt {\frac {5}{\pi }}}\,(3\cos ^{2}\theta -1)\\Y_{2}^{1}(\theta ,\varphi )&={\frac {-1}{2}}{\sqrt {\frac {15}{2\pi }}}\,\sin \theta \,\cos \theta \,e^{i\varphi }\\Y_{2}^{2}(\theta ,\varphi )&={\frac {1}{4}}{\sqrt {\frac {15}{2\pi }}}\,\sin ^{2}\theta \,e^{2i\varphi }\end{aligned}}}
高次元
古典的な球面調和関数は、 3次元ユークリッド空間内の 単位球面上の複素数値関数として定義されます。球面調和関数は 、次のように 高次元ユークリッド空間に一般化することができ、関数 につながります 。 [24] P ℓ を 、 n 個 の実変数を持つ 次数 ℓの 複素数値同次多項式 の 空間 を表すものと し 、ここでは関数 として考えます 。つまり、多項式 p は、 任意の 実数 に対して 、
S
2
{\displaystyle S^{2}}
R
3
{\displaystyle \mathbb {R} ^{3}}
R
n
{\displaystyle \mathbb {R} ^{n}}
S
n
−
1
→
C
{\displaystyle S^{n-1}\to \mathbb {C} }
R
n
→
C
{\displaystyle \mathbb {R} ^{n}\to \mathbb {C} }
λ
∈
R
{\displaystyle \lambda \in \mathbb {R} }
p
(
λ
x
)
=
λ
ℓ
p
(
x
)
.
{\displaystyle p(\lambda \mathbf {x} )=\lambda ^{\ell }p(\mathbf {x} ).}
A ℓ を すべての調和多項式 からなる P ℓ の部分空間とします 。 これら
は(正則) 立体球面調和関数 です。H ℓ を A ℓ から
の制限によって得られる 単位球面上の関数の空間とします
。
A
ℓ
:=
{
p
∈
P
ℓ
∣
Δ
p
=
0
}
.
{\displaystyle \mathbf {A} _{\ell }:=\{p\in \mathbf {P} _{\ell }\,\mid \,\Delta p=0\}\,.}
S
n
−
1
:=
{
x
∈
R
n
∣
|
x
|
=
1
}
{\displaystyle S^{n-1}:=\{\mathbf {x} \in \mathbb {R} ^{n}\,\mid \,\left|x\right|=1\}}
H
ℓ
:=
{
f
:
S
n
−
1
→
C
∣
for some
p
∈
A
ℓ
,
f
(
x
)
=
p
(
x
)
for all
x
∈
S
n
−
1
}
.
{\displaystyle \mathbf {H} _{\ell }:=\left\{f:S^{n-1}\to \mathbb {C} \,\mid \,{\text{ for some }}p\in \mathbf {A} _{\ell },\,f(\mathbf {x} )=p(\mathbf {x} ){\text{ for all }}\mathbf {x} \in S^{n-1}\right\}.}
次のプロパティが保持されます。
ストーン・ワイエルシュトラス の定理により、 空間 H ℓ の和は、一様位相 に関する 上の連続関数の 集合内で 稠密で ある 。結果として、これらの空間の和は、球面上の二乗可積分関数の空間 L 2 ( S n −1 ) でも稠密である。したがって、球面上のすべての二乗可積分関数は、球面調和関数の級数に一意に分解され、級数は L 2 の 意味で収束する。
C
(
S
n
−
1
)
{\displaystyle C(S^{n-1})}
S
n
−
1
{\displaystyle S^{n-1}}
全てのf ∈ H ℓ に対して 、 Δ S n −1 は S n −1 上の ラプラス・ベルトラミ演算子 である。この演算子は3次元のラプラス演算子の角度部分の類似物である。すなわち、 n 次元 のラプラス演算子は次のように分解される。
Δ
S
n
−
1
f
=
−
ℓ
(
ℓ
+
n
−
2
)
f
.
{\displaystyle \Delta _{S^{n-1}}f=-\ell (\ell +n-2)f.}
∇
2
=
r
1
−
n
∂
∂
r
r
n
−
1
∂
∂
r
+
r
−
2
Δ
S
n
−
1
=
∂
2
∂
r
2
+
n
−
1
r
∂
∂
r
+
r
−
2
Δ
S
n
−
1
{\displaystyle \nabla ^{2}=r^{1-n}{\frac {\partial }{\partial r}}r^{n-1}{\frac {\partial }{\partial r}}+r^{-2}\Delta _{S^{n-1}}={\frac {\partial ^{2}}{\partial r^{2}}}+{\frac {n-1}{r}}{\frac {\partial }{\partial r}}+r^{-2}\Delta _{S^{n-1}}}
ストークスの定理 と前述の性質から 、空間 H ℓ は L 2 ( S n −1 ) の内積に関して直交すること がわかります 。つまり、 f ∈ H ℓ および k ≠ ℓ に対して g ∈ H k です。
∫
S
n
−
1
f
g
¯
d
Ω
=
0
{\displaystyle \int _{S^{n-1}}f{\bar {g}}\,\mathrm {d} \Omega =0}
逆に、空間 H ℓ は Δ S n −1 の固有空間そのものである 。特に、 スペクトル定理を リースポテンシャル に適用すると 、空間 H ℓ がL 2 ( S n −1 ) において対直交かつ完全であることが証明される 。
Δ
S
n
−
1
−
1
{\displaystyle \Delta _{S^{n-1}}^{-1}}
あらゆる同次多項式p∈Pℓは、pj∈Aj の 形 で 一意 に 表さ れる [ 25 ] 。 特に 、
p
(
x
)
=
p
ℓ
(
x
)
+
|
x
|
2
p
ℓ
−
2
+
⋯
+
{
|
x
|
ℓ
p
0
ℓ
e
v
e
n
|
x
|
ℓ
−
1
p
1
(
x
)
ℓ
o
d
d
{\displaystyle p(x)=p_{\ell }(x)+|x|^{2}p_{\ell -2}+\cdots +{\begin{cases}|x|^{\ell }p_{0}&\ell {\rm {\ even}}\\|x|^{\ell -1}p_{1}(x)&\ell {\rm {\ odd}}\end{cases}}}
dim
H
ℓ
=
(
n
+
ℓ
−
1
n
−
1
)
−
(
n
+
ℓ
−
3
n
−
1
)
=
(
n
+
ℓ
−
2
n
−
2
)
+
(
n
+
ℓ
−
3
n
−
2
)
.
{\displaystyle \dim \mathbf {H} _{\ell }={\binom {n+\ell -1}{n-1}}-{\binom {n+\ell -3}{n-1}}={\binom {n+\ell -2}{n-2}}+{\binom {n+\ell -3}{n-2}}.}
高次元の球面調和関数の直交基底は、 変数分離 法によって 、球面ラプラシアンに対するシュトゥルム・リウヴィル問題を解くことによって
帰納的に 構築することができる。
ここで、 φは S n −1 上の球面座標系の軸座標である 。このような手順の最終結果は [26] であり
、ここでインデックスは | ℓ 1 | ≤ ℓ 2 ≤ ⋯ ≤ ℓ n −1 を満たし、固有値は − ℓ n −1 ( ℓ n −1 + n −2)である。積の関数は、 ルジャンドル関数 によって定義される。
Δ
S
n
−
1
=
sin
2
−
n
φ
∂
∂
φ
sin
n
−
2
φ
∂
∂
φ
+
sin
−
2
φ
Δ
S
n
−
2
{\displaystyle \Delta _{S^{n-1}}=\sin ^{2-n}\varphi {\frac {\partial }{\partial \varphi }}\sin ^{n-2}\varphi {\frac {\partial }{\partial \varphi }}+\sin ^{-2}\varphi \Delta _{S^{n-2}}}
Y
ℓ
1
,
…
ℓ
n
−
1
(
θ
1
,
…
θ
n
−
1
)
=
1
2
π
e
i
ℓ
1
θ
1
∏
j
=
2
n
−
1
j
P
¯
ℓ
j
ℓ
j
−
1
(
θ
j
)
{\displaystyle Y_{\ell _{1},\dots \ell _{n-1}}(\theta _{1},\dots \theta _{n-1})={\frac {1}{\sqrt {2\pi }}}e^{i\ell _{1}\theta _{1}}\prod _{j=2}^{n-1}{}_{j}{\bar {P}}_{\ell _{j}}^{\ell _{j-1}}(\theta _{j})}
j
P
¯
L
ℓ
(
θ
)
=
2
L
+
j
−
1
2
(
L
+
ℓ
+
j
−
2
)
!
(
L
−
ℓ
)
!
sin
2
−
j
2
(
θ
)
P
L
+
j
−
2
2
−
(
ℓ
+
j
−
2
2
)
(
cos
θ
)
.
{\displaystyle {}_{j}{\bar {P}}_{L}^{\ell }(\theta )={\sqrt {{\frac {2L+j-1}{2}}{\frac {(L+\ell +j-2)!}{(L-\ell )!}}}}\sin ^{\frac {2-j}{2}}(\theta )P_{L+{\frac {j-2}{2}}}^{-\left(\ell +{\frac {j-2}{2}}\right)}(\cos \theta )\,.}
表現理論との関連
ℓ 次の球面調和関数の 空間 H ℓは 、点( SO(3) )の周りの回転の 対称 群とその二重被覆 SU(2)の 表現 である 。実際、回転は2次元 球面 に作用し、したがって、 ψを
球面 調和関数、 ρを 回転とする関数合成によって
H ℓ にも作用する。表現 H ℓは SO(3)の 既約表現 である。 [27]
ψ
↦
ψ
∘
ρ
−
1
{\displaystyle \psi \mapsto \psi \circ \rho ^{-1}}
H ℓ の元は、 A ℓ の元の球面、すなわち 3次元ユークリッド空間 R 3 上の次数 ℓ の同次調和多項式への制限として生じます。 ψ ∈ A ℓ の 分極 により、インデックス上で対称な
係数が存在し、これは要件によって一意に決定されます。 ψ が調和である
という条件は、 テンソルがすべてのインデックスのペア上で トレース フリーでなければならない という主張と同等です。したがって、 SO(3) の既約表現として 、 H ℓ は 次数ℓ の トレースレス対称テンソル の空間に同型です 。
ψ
i
1
…
i
ℓ
{\displaystyle \psi _{i_{1}\dots i_{\ell }}}
ψ
(
x
1
,
…
,
x
n
)
=
∑
i
1
…
i
ℓ
ψ
i
1
…
i
ℓ
x
i
1
⋯
x
i
ℓ
.
{\displaystyle \psi (x_{1},\dots ,x_{n})=\sum _{i_{1}\dots i_{\ell }}\psi _{i_{1}\dots i_{\ell }}x_{i_{1}}\cdots x_{i_{\ell }}.}
ψ
i
1
…
i
ℓ
{\displaystyle \psi _{i_{1}\dots i_{\ell }}}
より一般的には、高次元でも同様のことが言えます。n球面 上 の球面調和関数の 空間 H ℓ は 、トレースレス対称 ℓテンソルに対応する SO( n +1) の既約表現です。しかし、 SO(2) と SO(3) の既約テンソル表現はすべてこの種類のものですが 、高次元の特殊直交群には、このようには生じない追加の既約表現があります。
特殊直交群には、テンソル表現ではない追加の スピン表現 があり、 通常は 球面調和関数ではありません 。例外はSO(3)の スピン表現です。厳密に言えば、これらはSO(3)の 二重被覆 SU(2)の表現です。次に、SU(2)は単位 四元数 の群と同一視され、したがって 3次元球面 と一致します。3次元球面上の球面調和関数の空間は、四元数の乗算による作用に関して、SO(3)の特定のスピン表現です。
半球調和関数との関連
球面調和関数は2つの関数セットに分けることができます。 [28] 1つは半球関数(HSH)で、半球上で直交かつ完全です。もう1つは補完半球調和関数(CHSH)です。
一般化
2次元球面 の 角度 保存対称性は、 メビウス変換 群 PSL(2, C ) によって記述される。この群に関して、球面は通常の リーマン球面 と同等である。群 PSL(2, C ) は(適切な) ローレンツ群 と同型であり 、2次元球面へのその作用は、 ミンコフスキー空間 におけるローレンツ群の 天球 への作用と一致する。ローレンツ群に対する球面調和関数の類似物は、 超幾何級数 によって与えられる。さらに、球面調和関数は、超幾何級数で再表現することができ、 SO(3) = PSU(2)は PSL(2, C ) の サブグループ である 。
より一般的には、超幾何級数は任意 の対称空間 の対称性を記述するように一般化することができ、特に、超幾何級数は任意の リー群 に対して展開することができる。 [29] [30] [31] [32]
参照
ウィキメディア・コモンズには、球面調和関数 に関連するメディアがあります 。
注記
^ 3次元の球面調和関数に対するさまざまなアプローチの歴史的説明は、MacRobert 1967の第4章に記載されています。「ラプラス球面調和関数」という用語は一般的に使用されています。Courant & Hilbert 1962およびMeijer & Bauer 2004を参照してください。
^ ここで採用されている球面調和関数へのアプローチは、(Courant & Hilbert 1962、§V.8、§VII.5)に記載されています。
^ 物理的な応用では、無限遠で消える解が採用されることが多く、 A = 0 となります。これは球面調和関数の角度部分には影響しません。
^ Weisstein, Eric W. 「球面調和関数」 。mathworld.wolfram.com . 2023年5月10日 閲覧 。
^ エドモンズ 1957、§2.5
^ ホール 2013 セクション 17.6
^ ホール 2013 補題 17.16
^ ウィリアムズ、アール G. (1999)。 フーリエ音響:音響放射と近接場音響ホログラフィー 。カリフォルニア州サンディエゴ:アカデミックプレス 。ISBN 0080506909 . OCLC 181010993.
^ メサイア、アルバート (1999)。 量子力学:2巻を1冊にまとめた (2巻を1冊にまとめた完全版復刻版)。ニューヨーク州ミネオラ:ドーバー。pp . 520– 523。ISBN 0486409244 。
^ クロード・コーエン=タンヌージ;バーナード・ディウ;フランク・ラロエ (1996)。 量子力学 。スーザン・リード・ヘムリー訳。他。ワイリー・インターサイエンス: ワイリー。 ISBN 9780471569527 。
^ ab ブレイクリー、リチャード(1995)。 重力と磁気の応用におけるポテンシャル理論 。ケンブリッジ、イングランド、ニューヨーク:ケンブリッジ大学出版局。p. 113。ISBN 978-0521415088 。
^ Heiskanen と Moritz、物理測地学、1967 年、式 1-62
^ Weisstein, Eric W. 「Condon-Shortley Phase」 。mathworld.wolfram.com 。 2022年11月2日 閲覧 。
^ ホイッタカー&ワトソン 1927年、392ページ。
^ 例えば、Garg, A., Classical Electrodynamics in a Nutshell (Princeton University Press、2012) の付録 A を参照。
^ Li, Feifei; Braun, Carol; Garg, Anupam (2013)、「スピンに対する Weyl-Wigner-Moyal 形式論」、 Europhysics Letters 、 102 (6): 60006、 arXiv : 1210.4075 、 Bibcode :2013EL....10260006L、 doi :10.1209/0295-5075/102/60006、 S2CID 119610178
^ Edmonds, AR (1996). 量子力学における角運動量 . プリンストン大学出版局. p. 63.
^これは ℓ 次球面調和関数の任意の正規直交基底に対して有効です。単位べき調和関数の場合、 4 π の係数を除去する必要があります 。
^ ホイッタカー&ワトソン 1927年、395ページ
^ 未回収 1927
^ スタイン&ワイス 1971、§IV.2
^ ブリンク, DM; サッチラー, GR 角運動量 。オックスフォード大学出版局。146 ページ。
^ エレメンコ、ヤコブソン、ナディラシビリ 2007
^ ソロメンツェフ 2001; スタイン&ワイス 1971, §Iv.2
^参照: Axler, Sheldon、Ramey, Wade (1995) 『調和多項式とディリクレ型問題』 の系 1.8
^ 樋口 篤 (1987). 「N球面上の対称テンソル球面調和関数とド・シッター群SO(N,1)への応用」. 数理物理学ジャーナル . 28 (7): 1553– 1566. Bibcode :1987JMP....28.1553H. doi :10.1063/1.527513.
^ ホール 2013 補論 17.17
^ Zheng Y, Wei K, Liang B, Li Y, Chu X (2019-12-23). 「球面キャップ上のゼルニケのような関数:光学表面フィッティングとグラフィックスレンダリングにおける原理と応用」. Optics Express . 27 (26): 37180– 37195. Bibcode :2019OExpr..2737180Z. doi : 10.1364/OE.27.037180 . ISSN 1094-4087. PMID 31878503.
^ N. Vilenkin、 「特殊関数と群表現の理論」 、Am. Math. Soc. Transl.、vol. 22、(1968)。
^ JD Talman、 「特殊関数、群論的アプローチ 」(EP Wignerの講義に基づく)、WA Benjamin、ニューヨーク(1968年)。
^ W.ミラー、 「対称性と変数の分離」、 アディソン・ウェズリー、リーディング(1977年)。
^ A. Wawrzyńczyk、 「Group Representations and Special Functions」 、Polish Scientific Publishers。Warszawa (1984)。
参考文献
引用文献
クーラント、リチャード 、 ヒルベルト、デイヴィッド (1962)、 数理物理学の方法、第1巻 、ワイリー・インターサイエンス 。
エドモンズ、AR(1957)、 量子力学における角運動量 、プリンストン大学出版、 ISBN 0-691-07912-9
エレメンコ、アレクサンドル; ヤコブソン、ドミトリー; ナディラシビリ、ニコライ (2007)、「S² と R² 上の節点集合と節点領域について」、 Annales de l'Institut Fourier 、 57 (7): 2345– 2360、 doi : 10.5802/aif.2335 、 ISSN 0373-0956、 MR 2394544
ホール、ブライアン C. (2013)、 「数学者のための量子理論」 、Graduate Texts in Mathematics、vol. 267、Springer、 Bibcode :2013qtm..book.....H、 ISBN 978-1461471158
MacRobert, TM (1967)、 球面調和関数: 調和関数に関する基礎論文とその応用 、Pergamon Press 。
マイヤー、ポール・ハーマン・エルンスト、バウアー、エドモンド(2004)、 群論:量子力学への応用 、ドーバー、 ISBN 978-0-486-43798-9 。
Solomentsev, ED (2001) [1994]、「球面調和関数」、 数学百科事典 、 EMS Press 。
スタイン、エリアス 、 ワイス、グイド (1971)、 ユークリッド空間のフーリエ解析入門 、プリンストン、ニュージャージー:プリンストン大学出版、 ISBN 978-0-691-08078-9 。
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 。
ウィテカー、ET 、 ワトソン、GN (1927)、 現代分析講座 、 ケンブリッジ大学出版局 、p. 392 。
一般的な参考文献
EW ホブソン、 「球面調和関数と楕円調和関数の理論 」、(1955) チェルシー出版、 ISBN 978-0-8284-0104-3 。
C. Müller, Spherical Harmonics 、(1966) Springer、Lecture Notes in Mathematics、Vol. 17、 ISBN 978-3-540-03600-5 。
EUコンドンとGHショートリー、 「原子スペクトルの理論 」(1970年)ケンブリッジ大学出版局、 ISBN 0-521-09209-4 、 第3章を参照 。
JD ジャクソン 『古典電気力学』 、 ISBN 0-471-30932-X
アルバート・メサイア 『量子力学 』 第2巻(2000年)ドーバー。ISBN 0-486-40924-4 。
Press, WH; Teukolsky, SA; Vetterling, WT; Flannery, BP (2007)、「セクション 6.7. 球面調和関数」、Numerical Recipes: The Art of Scientific Computing (第 3 版)、ニューヨーク: Cambridge University Press、 ISBN 978-0-521-88068-8
DA Varshalovich、AN Moskalev、VK Khersonskii 角運動量の量子理論 、(1988) World Scientific Publishing Co.、シンガポール、 ISBN 9971-5-0107-4
Weisstein、Eric W. 「球面調和関数」 。MathWorld 。
Maddock, John、Boost.Math における球面調和関数
外部リンク