最初のいくつかの実球面調和関数の視覚的表現。青色の部分は関数が正である領域を表し、黄色の部分は負である領域を表します。表面の原点からの距離は、絶対値を示します。Y ℓ m ( θ 、 φ ) {\displaystyle Y_{\ell }^{m}(\theta ,\varphi )} 角度方向( θ 、 φ ) {\displaystyle (\theta ,\varphi )} 。 数学 や物理学 において、球面調和関数は 球面 上で定義される特殊な関数 です。これらは多くの科学分野で偏微分方程式 を解く際によく用いられます。球面調和関数表には、 よく用いられる球面調和関数の一覧が掲載されています。
球面調和関数は直交関数 の完全な集合、すなわち正規直交基底 を形成するため、球面上で定義された特定の関数は、これらの球面調和関数の和として表すことができます。これは、 フーリエ級数によって 円関数 (正弦関数と余弦関数)の和として表せる円上で定義された周期関数 と似ています。フーリエ級数の正弦関数と余弦関数と同様に、球面調和関数は(空間)角周波数 によって整理することができ、右の図の関数の行に示されています。さらに、球面調和関数は、 3次元回転群である SO(3) の既約表現 の基底関数 であり、 SO(3)の群論的 議論において中心的な役割を果たします。
球面調和関数は、球面領域におけるラプラス方程式 の解から得られます。ラプラス方程式の解となる関数は調和関数と呼ばれます。その名前にもかかわらず、球面調和関数は デカルト座標系 において最も単純な形をとり、次数 の同次多項式 として定義できます。ℓ {\displaystyle \ell } で( x 、 y 、 z ) {\displaystyle (x,y,z)} ラプラス方程式に従う。同次性を用いて半径依存性の因子を抽出すると、球面座標との関連がすぐに生じる。 r ℓ \displaystyle r^{\ell }} 上記の次数多項式からℓ {\displaystyle \ell } 残りの因子は球面角度座標の関数とみなすことができる。θ {\displaystyle \theta } そしてφ {\displaystyle \varphi } のみ、または同等に方向 単位ベクトル r {\displaystyle \mathbf {r} } specified by these angles. In this setting, they may be viewed as the angular portion of a set of solutions to Laplace's equation in three dimensions, and this viewpoint is often taken as an alternative definition. Notice, however, that spherical harmonics are not functions on the sphere which are harmonic with respect to the Laplace-Beltrami operator for the standard round metric on the sphere: the only harmonic functions in this sense on the sphere are the constants, since harmonic functions satisfy the Maximum principle . Spherical harmonics, as functions on the sphere, are eigenfunctions of the Laplace-Beltrami operator (see Higher dimensions ).
A specific set of spherical harmonics, denoted Y ℓ m ( θ , φ ) {\displaystyle Y_{\ell }^{m}(\theta ,\varphi )} or Y ℓ m ( r ) {\displaystyle Y_{\ell }^{m}({\mathbf {r} })} , are known as Laplace's spherical harmonics, as they were first introduced by Pierre Simon de Laplace in 1782.[ 1] These functions form an orthogonal system, and are thus basic to the expansion of a general function on the sphere as alluded to above.
Spherical harmonics are important in many theoretical and practical applications, including the representation of multipole electrostatic and electromagnetic fields , electron configurations , gravitational fields , geoids , the magnetic fields of planetary bodies and stars, and the cosmic microwave background radiation . In 3D computer graphics , spherical harmonics play a role in a wide variety of topics including indirect lighting (ambient occlusion , global illumination , precomputed radiance transfer , etc.) and modelling of 3D shapes.
History Pierre-Simon Laplace , 1749–1827Spherical harmonics were first investigated in connection with the Newtonian potential of Newton's law of universal gravitation in three dimensions. In 1782, Pierre-Simon de Laplace had, in his Mécanique Céleste , determined that the gravitational potential R 3 → R \displaystyle \mathbb {R} ^{3}\to \mathbb {R} } at a point x associated with a set of point masses m i located at points x i was given by
V ( x ) = ∑ i m i | x i − x | . {\displaystyle V(\mathbf {x} )=\sum _{i}{\frac {m_{i}}{|\mathbf {x} _{i}-\mathbf {x} |}}.}
Each term in the above summation is an individual Newtonian potential for a point mass. Just prior to that time, Adrien-Marie Legendre had investigated the expansion of the Newtonian potential in powers of r = | x | and r 1 = | x 1 | . He discovered that if r ≤ r 1 then
1 | x 1 − x | = P 0 ( cos γ ) 1 r 1 + P 1 ( cos γ ) r r 1 2 + P 2 ( cos γ ) 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 }
where γ is the angle between the vectors x and x 1 . The functions P i : [ − 1 , 1 ] → R {\displaystyle P_{i}:[-1,1]\to \mathbb {R} } are the Legendre polynomials , and they can be derived as a special case of spherical harmonics. Subsequently, in his 1782 memoir, Laplace investigated these coefficients using spherical coordinates to represent the angle γ between x 1 and x . (See Legendre polynomials § Applications for more detail.)
In 1867, William Thomson (Lord Kelvin) and Peter Guthrie Tait introduced the solid spherical harmonics in their Treatise on Natural Philosophy , and also first introduced the name of "spherical harmonics" for these functions. The solid harmonics were homogeneous polynomial solutions R 3 → R \displaystyle \mathbb {R} ^{3}\to \mathbb {R} } of Laplace's equation ∂ 2 u ∂ x 2 + ∂ 2 u ∂ y 2 + ∂ 2 u ∂ z 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.} By examining Laplace's equation in spherical coordinates, Thomson and Tait recovered Laplace's spherical harmonics. (See Harmonic polynomial representation .) The term "Laplace's coefficients" was employed by William Whewell to describe the particular system of solutions introduced along these lines, whereas others reserved this designation for the zonal spherical harmonics that had properly been introduced by Laplace and Legendre.
19世紀に発展したフーリエ級数は、 熱方程式 や波動方程式 の解法など、矩形領域におけるさまざまな物理問題の解決を可能にした。これは、関数を三角関数の級数に展開することによって実現できた。フーリエ級数における三角関数は 弦 の基本振動モードを表すのに対し、球面調和関数はほぼ同じように球の 基本振動モードを表す。フーリエ級数の理論の多くの側面は、三角関数ではなく球面調和関数による展開によって一般化できた。さらに、三角関数が複素指数関数 として等価に表せるのと同様に、球面調和関数も複素数値関数として等価な形式を持つことができた。これは、ラプラスとルジャンドルが最初に研究した天体力学 のような、球対称性 を持つ問題にとって大きな利点となった。
物理学における球面調和関数の普及は、20世紀に誕生した量子力学 におけるその重要性の基礎を築いた。(複素数値の)球面調和関数S 2 → C \displaystyle S^{2}\to \mathbb {C} } は軌道角運動量 演算子 の二乗の固有関数 である。− 私 ℏ r × ∇ 、 {\displaystyle -i\hbar \mathbf {r} \times \nabla ,} したがって、それらは原子軌道 の 異なる量子化された 配置を表している。
ラプラスの球面調和関数実数(ラプラス)球面調和関数Y ℓ m {\displaystyle Y_{\ell m}} のためにℓ = 0 、 … 、 4 {\displaystyle \ell =0,\dots ,4} (上から下へ)そしてm = 0 、 … 、 ℓ {\displaystyle m=0,\dots ,\ell } (左から右へ)。ゾーン調和、セクター調和、テッセラル調和は、それぞれ左端の列、主対角線、その他の場所に示されています。(負の次数調和)Y ℓ ( − m ) {\displaystyle Y_{\ell (-m)}} z 軸を中心に回転させた状態で表示されます。90 ∘ / m 90°/m (正の次数のものに関して)。高調波の視覚的な分かりやすさを向上させるために回転を追加しました。 実数球面調和関数の別図Y ℓ m {\displaystyle Y_{\ell m}} 。 ラプラス方程式は 、スカラー場f のラプラシアン がゼロであることを規定する。(ここでスカラー場は複素数、すなわち(滑らかな)関数に対応するものと理解される。)f : R 3 → C {\displaystyle f:\mathbb {R} ^{3}\to \mathbb {C} } .)球面座標 では次のようになります。[ 2 ]
∇ 2 f = 1 r 2 ∂ ∂ r ( r 2 ∂ f ∂ r ) + 1 r 2 罪 θ ∂ ∂ θ ( 罪 θ ∂ f ∂ θ ) + 1 r 2 罪 2 θ ∂ 2 f ∂ φ 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 つの微分方程式が得られます。 1 R d d r ( r 2 d R d r ) = λ 、 1 Y 1 罪 θ ∂ ∂ θ ( 罪 θ ∂ Y ∂ θ ) + 1 Y 1 罪 2 θ ∂ 2 Y ∂ φ 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 .} 2 番目の式は、 Y が Y ( θ , φ ) = Θ( θ ) Φ( φ ) の形である と仮定すると、簡略化できます。2 番目の式に再び変数分離法を適用すると、次の微分方程式の組が得られます。
1 Φ d 2 Φ d φ 2 = − m 2 {\displaystyle {\frac {1}{\Phi }}{\frac {d^{2}\Phi }{d\varphi ^{2}}}=-m^{2}} λ 罪 2 θ + 罪 θ Θ d d θ ( 罪 θ d Θ d θ ) = m 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 に対して。先験的に 、m は複素定数ですが、Φ は周期が 2π を 割り切る周期関数 でなければならないため、m は 必然的に整数であり、Φは複素指数関数 e ± imφ の線形結合です。解関数Y ( θ , φ ) は、 θ = 0, π の球の極で正則です。この正則性を領域の境界点で第 2 式の解Θ に課すことは、パラメータλ を ℓ ≥ | m |を 満たす 非負の整数に対してλ = ℓ ( ℓ + 1) の形に強制するシュトゥルム・リウヴィル問題です 。これは、軌道角運動量 の観点からも以下で 説明されています。さらに、変数変換t = cos θ により、この方程式はルジャンドル方程式に変換され、その解は 関連するルジャンドル多項式 P m ℓ (cos θ ) の倍数です。最後に、 R の方程式はR ( r ) = A r ℓ + B r − ℓ − 1 の形の解を持ち、 R 3 全体で解が正則である必要があるため、B = 0 と なります。[ 3 ]
ここでは、解は特別な形式Y ( θ , φ ) = Θ( θ ) Φ( φ )を持つと仮定した。 ℓ の値が与えられた場合、この形式の独立した解は2 ℓ + 1 個あり、それぞれ− ℓ ≤ m ≤ ℓ を満たす整数m ごとに 1 つずつ存在する。 これらの角度解Y ℓ m : S 2 → C {\displaystyle Y_{\ell }^{m}:S^{2}\to \mathbb {C} } これらは、ここでは複素指数関数 として表される三角関数 と、関連するルジャンドル多項式の積です。
Y ℓ m ( θ 、 φ ) = N e 私 m φ P ℓ m ( コス θ ) {\displaystyle Y_{\ell }^{m}(\theta ,\varphi )=Ne^{im\varphi }P_{\ell }^{m}(\cos {\theta })}
満たす r 2 ∇ 2 Y ℓ m ( θ 、 φ ) = − ℓ ( ℓ + 1 ) Y ℓ m ( θ 、 φ ) 。 {\displaystyle r^{2}\nabla ^{2}Y_{\ell }^{m}(\theta ,\varphi )=-\ell (\ell +1)Y_{\ell }^{m}(\theta ,\varphi ).}
ここY ℓ m : S 2 → C {\displaystyle Y_{\ell }^{m}:S^{2}\to \mathbb {C} } は、次数ℓ 、階数 m の球面調和関数 と呼ばれます。P ℓ m : [ − 1 、 1 ] → R {\displaystyle P_{\ell }^{m}:[-1,1]\to \mathbb {R} } は関連ルジャンドル多項式 、N は正規化定数、[ 4 ] であり、θ とφ は それぞれ余緯度と経度を表します。特に、余緯度 θ 、つまり極角は、北極で0 、赤道で π /2 、南極で π まで変化し、経度φ 、 つまり 方位角は、 0 ≤ φ < 2 π のすべての値をとることができます。固定された整数ℓ に対して、すべての解Y ( θ , φ ) 、Y : S 2 → C {\displaystyle Y:S^{2}\to \mathbb {C} } 固有値問題の r 2 ∇ 2 Y = − ℓ ( ℓ + 1 ) Y {\displaystyle r^{2}\nabla ^{2}Y=-\ell (\ell +1)Y} は、 Y ℓ m : S 2 → C {\displaystyle Y_{\ell }^{m}:S^{2}\to \mathbb {C} } 実際、そのような解のいずれに対しても、r ℓ Y ( θ , φ )は球面座標における同次多項式 の表現である。R 3 → C {\displaystyle \mathbb {R} ^{3}\to \mathbb {C} } これは調和的であるため(下記参照)、次元を数えると、そのような線形独立な多項式 が 2ℓ + 1 個存在することがわかります。
一般的な解決策f : R 3 → C {\displaystyle f:\mathbb {R} ^{3}\to \mathbb {C} } ラプラス方程式 へΔ f = 0 {\displaystyle \Delta f=0} 原点を中心とする球面では、適切なスケール係数r ℓ を乗じた球面調和関数の線形結合 である。
f ( r 、 θ 、 φ ) = ∑ ℓ = 0 ∞ ∑ m = − ℓ ℓ f ℓ m r ℓ Y ℓ m ( θ 、 φ ) 、 {\displaystyle f(r,\theta ,\varphi )=\sum _{\ell =0}^{\infty }\sum _{m=-\ell }^{\ell }f_{\ell }^{m}r^{\ell }Y_{\ell }^{m}(\theta ,\varphi ),}
どこでf ℓ m ∈ C {\displaystyle f_{\ell }^{m}\in \mathbb {C} } は定数であり、因子r ℓ Y ℓ m は (正則 )固体高調波として知られています。 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 ] L = − 私 ℏ ( x × ∇ ) = L x 私 + 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} .} ħは 量子力学において慣例的に用いられるものであり、ħ = 1 となる単位系で作業するのが便利である。球面調和関数は軌道角運動量の二乗の固有関数である。 L 2 = − r 2 ∇ 2 + ( r ∂ ∂ r + 1 ) r ∂ ∂ r = − 1 罪 θ ∂ ∂ θ 罪 θ ∂ ∂ θ − 1 罪 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 = − 私 ( x ∂ ∂ y − y ∂ ∂ x ) = − 私 ∂ ∂ φ 。 {\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の 重み付き ヒルベルト空間 上の密に定義された 自己共役演算子 です。 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 .} さらに、L 2は 正演算子 である。
Y が L2 とLz の 同時固有関数である場合、定義により L 2 Y = λ Y L z Y = m Y {\displaystyle {\begin{aligned}\mathbf {L} ^{2}Y&=\lambda Y\\L_{z}Y&=mY\end{aligned}}} ある実数m とλ に対して。ここで、m は 実際には整数でなければなりません。なぜなら、Y は 2π を割り切る周期を持つ座標φ で周期的になければならないからです。さらに、 L 2 = L x 2 + L y 2 + L z 2 {\displaystyle \mathbf {L} ^{2}=L_{x}^{2}+L_{y}^{2}+L_{z}^{2}} また、L x 、L y 、L zはそれぞれ自己共役であるため、 λ ≥ m 2 が成り立つ。
この結合固有空間をE λ , m と表記し、昇降演算子を 次のように 定義する。L + = L x + 私 L y L − = L x − 私 L y {\displaystyle {\begin{aligned}L_{+}&=L_{x}+iL_{y}\\L_{-}&=L_{x}-iL_{y}\end{aligned}}} すると、L + とL − は L 2 と可換になり、 L + 、L − 、L z によって生成されるリー代数は、2 位の特別な線形リー代数 になります。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}.} したがって、L + : E λ , m → E λ , m +1 (これは「上昇演算子」です) およびL − : E λ , m → E λ , m −1 (これは「下降演算子」です) となります。特に、L k + : E λ , m → E λ , m + k は 、 k が十分に大きい場合、ゼロでなければなりません。なぜなら、 λ ≥ m 2 という不等式が各非自明な結合固有空間で成り立たなければならないからです。Y ∈ E λ , m を ゼロでない結合固有関数とし、k を 最小の整数とします 。 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.} したがって、λ = ℓ ( ℓ + 1) は 正の整数ℓ = m + k に対して成り立つ。
上記はすべて球面座標表現で計算されたものであり、⟨ θ 、 φ | 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 \mathbf {x} \in \mathbb {R} ^{3}} 私たちはこう言いますp {\displaystyle p} 調和的 である場合 Δ p = 0 、 {\displaystyle \Delta p=0,} どこΔ {\displaystyle \Delta } はラプラシアン です。次に、各についてℓ {\displaystyle \ell } 定義する A ℓ = { 調和多項式 R 3 → C 次数が同質である ℓ } 。 {\displaystyle \mathbf {A} _{\ell }=\left\{{\text{harmonic polynomials }}\mathbb {R} ^{3}\to \mathbb {C} {\text{ that are homogeneous of degree }}\ell \right\}.}
例えば、ℓ = 1 {\displaystyle \ell =1} 、A 1 {\displaystyle \mathbf {A} _{1}} これは、すべての線形関数の3次元空間です。R 3 → C {\displaystyle \mathbb {R} ^{3}\to \mathbb {C} } なぜなら、そのような関数はすべて自動的に調和的になるからである。一方、ℓ = 2 {\displaystyle \ell =2} つまり、5次元空間が存在するということです。 A 2 = スパン C ( x 1 x 2 、 x 1 x 3 、 x 2 x 3 、 x 1 2 − x 2 2 、 2 x 3 2 − x 1 2 − x 2 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},\,2x_{3}^{2}-x_{1}^{2}-x_{2}^{2}).}
いかなる場合でもℓ {\displaystyle \ell } 空間H ℓ {\displaystyle \mathbf {H} _{\ell }} 次数 の球面調和関数ℓ {\displaystyle \ell } それは球体に対する制約の空間にすぎないS 2 {\displaystyle S^{2}} 要素のA ℓ {\displaystyle \mathbf {A} _{\ell }} [ 6 ] 序論で示唆したように、この視点はおそらく「球面調和関数」(つまり、調和関数の球面への制限) という 用語の起源である。
例えば、任意のc ∈ C {\displaystyle c\in \mathbb {C} } 式 p ( x 1 、 x 2 、 x 3 ) = c ( x 1 + 私 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 罪 ( θ ) ℓ ( コス ( φ ) + 私 罪 ( φ ) ) ℓ 、 {\displaystyle p(\theta ,\varphi )=c\sin(\theta )^{\ell }(\cos(\varphi )+i\sin(\varphi ))^{\ell },} これは次のように書き換えることができます p ( θ 、 φ ) = c ( 1 − コス 2 ( θ ) ) ℓ e 私 ℓ φ 。 {\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 ).} [ 7 ] (特殊なケース 。)
慣例
直交性と正規化 ラプラス球面調和関数には、いくつかの異なる正規化方法が一般的に用いられている。S 2 → C {\displaystyle S^{2}\to \mathbb {C} } これらの正規化の選択は、コンドン・ショートリー位相の位置とは無関係である。
球面調和関数のグラフY ℓ m ( θ 、 φ ) {\displaystyle Y_{\ell }^{m}(\theta ,\varphi )} とℓ = 2 {\displaystyle \ell =2} そしてm = 1 {\displaystyle m=1} そしてφ = π {\displaystyle \varphi =\pi } 複素平面において− 2 − 2 私 {\displaystyle -2-2i} に2 + 2 私 {\displaystyle 2+2i} Mathematica 13.1の関数ComplexPlot3Dで作成された色を使用 コンドン・ショートリー段階を含めるとP ℓ m {\displaystyle P_{\ell }^{m}} 正規直交複素球面調和関数は Y ℓ m ( θ 、 φ ) = ( 2 ℓ + 1 ) 4 π ( ℓ − m ) ! ( ℓ + m ) ! P ℓ m ( コス θ ) e 私 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 }.} これは本稿で使用されている慣例であり、音響学では一般的である。[ 11 ]
同じ正規直交関数は量子力学ではしばしば[ 12 ] [ 13 ]のように記述される。 Y ℓ m ( θ 、 φ ) = ( − 1 ) m ( 2 ℓ + 1 ) 4 π ( ℓ − m ) ! ( ℓ + m ) ! P ~ ℓ m ( コス θ ) e 私 m φ 、 {\displaystyle Y_{\ell }^{m}(\theta ,\varphi )=(-1)^{m}{\sqrt {{\frac {(2\ell +1)}{4\pi }}{\frac {(\ell -m)!}{(\ell +m)!}}}}\,{\widetilde {P}}_{\ell }^{m}(\cos {\theta })\,e^{im\varphi },} どこP ~ ℓ m {\displaystyle {\widetilde {P}}_{\ell }^{m}} これは、コンドン・ショートリー位相を含まない関連ルジャンドル関数を表します。これにより、位相を二重に数えることを防ぎます。
コンドン・ショートリー位相のどちらの配置でも、正規直交球面調和関数は以下を満たす。 ∫ θ = 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'},} ここでδ ij はクロネッカーのデルタ であり、d Ω = sin( θ ) dφ dθ で ある。[ 9 ] この正規化は、確率が正規化されることを保証するため、量子力学で使用される。 ∫ | Y ℓ m | 2 d Ω = 1. {\displaystyle \int {|Y_{\ell }^{m}|^{2}d\Omega }=1.}
測地学 [ 14 ] とスペクトル解析の分野では、対応する4 π {\displaystyle 4\pi } -正規化された高調波。上記と同じ位相規約で記述すると、これらは Y ℓ m ( θ 、 φ ) = ( 2 ℓ + 1 ) ( ℓ − m ) ! ( ℓ + m ) ! P ℓ m ( コス θ ) e 私 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'}.} 測地学の応用では、この正規化は一般的にコンドン・ショートリー位相の省略と組み合わせて行われる。[ 15 ]
一方、磁気学 [ 14 ] の分野では、シュミット半正規化高調波がよく用いられる。上記と同じ位相規約で記述すると、これらは次のようになる。 Y ℓ m ( θ 、 φ ) = ( ℓ − m ) ! ( ℓ + m ) ! P ℓ m ( コス θ ) e 私 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 ),} ここで、上付き文字*は 複素共役 を表します。同様に、( − 1 ) m = ( − 1 ) − m {\displaystyle (-1)^{m}=(-1)^{-m}} 整数m {\displaystyle m} これは次のように書くこともできます( − 1 ) − m {\displaystyle (-1)^{-m}} この方程式は、球面調和関数とウィグナーD行列 の関係からも導かれる。[ 9 ]
複素球面調和関数が正規化され位相が調整されていると仮定すると、 Y ℓ m * = ( − 1 ) m Y ℓ − m 、 {\displaystyle Y_{\ell }^{m}{}^{*}=(-1)^{m}Y_{\ell }^{-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 = { 私 2 ( Y ℓ m − ( − 1 ) m Y ℓ − m ) もし m < 0 Y ℓ 0 もし m = 0 1 2 ( Y ℓ − m + ( − 1 ) m Y ℓ m ) もし m > 0. = { 私 2 ( Y ℓ − | m | − ( − 1 ) m Y ℓ | m | ) もし m < 0 Y ℓ 0 もし m = 0 1 2 ( Y ℓ − | m | + ( − 1 ) m Y ℓ | m | ) もし m > 0. = { 2 ( − 1 ) m ℑ [ Y ℓ | m | ] もし m < 0 Y ℓ 0 もし m = 0 2 ( − 1 ) m ℜ [ Y ℓ m ] もし 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 | − 私 Y ℓ 、 − | m | ) もし m < 0 Y ℓ 0 もし m = 0 ( − 1 ) m 2 ( Y ℓ | m | + 私 Y ℓ 、 − | m | ) もし 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}}}
実球面調和関数Y ℓ m : S 2 → R {\displaystyle Y_{\ell m}:S^{2}\to \mathbb {R} } これらはテッセラル球面調和関数 とも呼ばれることがある。[ 16 ] これらの関数は複素関数と同じ正規直交性を持つ。Y ℓ m : S 2 → C {\displaystyle Y_{\ell }^{m}:S^{2}\to \mathbb {C} } 上記。実球面調和関数Y ℓ m {\displaystyle Y_{\ell m}} m > 0 の場合はコサイン型、 m < 0 の場合はサイン型と呼ばれます。その理由は、関数をルジャンドル多項式で表すとわかります。 Y ℓ m = { ( − 1 ) m 2 2 ℓ + 1 4 π ( ℓ − | m | ) ! ( ℓ + | m | ) ! P ℓ | m | ( コス θ ) 罪 ( | m | φ ) もし m < 0 2 ℓ + 1 4 π P ℓ 0 ( コス θ ) もし m = 0 ( − 1 ) m 2 2 ℓ + 1 4 π ( ℓ − m ) ! ( ℓ + m ) ! P ℓ m ( コス θ ) コス ( m φ ) もし 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 }^{0}(\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}}} ここP ℓ m {\displaystyle P_{\ell }^{m}} は、コンドン・ショートリー位相を含む関連ルジャンドル関数を表します。位相の位置が異なると、表示される符号も一貫して変更する必要があります。[ 8 ]
同じ正弦と余弦の係数は、デカルト座標系を扱う次の小節でも見られます。
実数球面調和関数の一覧については、こちらを ご覧ください。ℓ = 4 {\displaystyle \ell =4} これは、上記の方程式の出力と一致していることがわかる。
対称性 球面調和関数は、空間反転(パリティ)と回転の操作において、深く重要な特性を持つ。
パリティ 球面調和関数は明確なパリティを持ちます。つまり、原点に関する反転に関して偶関数か奇関数かのいずれかになります。反転は演算子で表されます。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} ).}
球面角の観点から見ると、パリティは座標を持つ点を変換します。{ θ 、 φ } {\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 を与え、指数関数からは(−1) m が得られ、これらを合わせると球面調和関数のパリティは(−1) ℓ となります。)
実球面調和関数、および高次元の球面調和関数ではパリティが引き続き成り立ちます。次数ℓ の球面調和関数に点反射を適用すると、符号が (−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 ).}
この展開は平均二乗収束の意味で成り立つ。つまり球面のL2における収束であり、
リム N → ∞ ∫ 0 2 π ∫ 0 π | f ( θ 、 φ ) − ∑ ℓ = 0 N ∑ m = − ℓ ℓ f ℓ m Y ℓ m ( θ 、 φ ) | 2 罪 θ 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 θ 罪 θ 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} } )その特性は持っていません。
完全性関係 完全性関係は、すべてのインデックスの和とディラックのデルタ関数 との間の分布等式です。球面調和関数の場合、ディラックのデルタは、球面上の方位角と極角の2つのディラックのデルタ関数のテンソル積です。関係は[ 19 ] で示されます。
∑ l = 0 ∞ ∑ m = − l l Y l m * ( θ ′ 、 ϕ ′ ) Y l m ( θ 、 ϕ ) = δ ( ϕ − ϕ ′ ) δ ( コス θ − コス θ ′ ) 。 {\displaystyle \sum _{l=0}^{\infty }\sum _{m=-l}^{l}Y_{lm}^{*}(\theta \,',\,\phi \,')Y_{lm}(\theta ,\,\phi )=\delta (\phi -\phi \,')\,\delta (\cos \theta -\cos \theta \,').}
スペクトル分析
信号処理におけるパワースペクトル 信号処理の 文献では、関数f の全パワーは、関数の二乗の積分をその定義域の面積で割ったものとして定義されています。実数単位パワー球面調和関数の直交性を利用すると、単位球上で定義された関数の全パワーが、 パーセバルの定理 の一般化によってそのスペクトル係数と関連付けられることが容易に確認できます(ここでは、定理はシュミット半正規化調和関数について述べられており、直交調和関数では関係が若干異なります)。[ 20 ] [ 21 ]
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 の微分可能性の次数と関連しています。
代数的性質
加法定理 非常に興味深く、有用な数学的結果の一つに、球面調和関数の加法定理 がある。球面座標を持つ2つのベクトルr とr′が与えられた場合、 ( r 、 θ 、 φ ) {\displaystyle (r,\theta ,\varphi )} そして( r 、 θ ′ 、 φ ′ ) {\displaystyle (r,\theta ',\varphi ')} それぞれ、角度γ {\displaystyle \gamma } 両者の関係は、 コス γ = コス θ ′ コス θ + 罪 θ 罪 θ ′ コス ( φ − φ ′ ) {\displaystyle \cos \gamma =\cos \theta '\cos \theta +\sin \theta \sin \theta '\cos(\varphi -\varphi ')} ここで、右辺に現れる三角関数の役割は球面調和関数によって担われ、左辺の三角関数の役割はルジャンドル多項式 によって担われます。
加法定理は [ 23 ] と述べている。
ここで、Pℓ は 次数ℓ のルジャンドル多項式 です。この式は、実数および複素数の調和関数の両方に有効です。[ 24 ] この結果は、単位球内のポアソン核 の性質を使用して解析的に証明することも、ベクトルyを z 軸に沿って回転させてから右辺を直接計算することによって幾何学的に証明することもできます。[ 25 ]
特に、x = y の場合、これはウンソルドの定理[ 26 ]を与える。 ∑ 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 }}} これは、恒等式cos 2 θ + sin 2 θ = 1 を 2 次元に一般化したものである。
展開式(1 )において、左辺P ℓ ( x ⋅ y ) {\displaystyle P_{\ell }(\mathbf {x} \cdot \mathbf {y} )} は、次数ℓ のゾーン球面調和関数 の定数倍である。この観点から、より高次元への以下の一般化が得られる。Y j を n 次元 球面上の次数ℓ の球面調和関数の空間H ℓ の任意の正規直交基底とする。すると、Z x ( ℓ ) {\displaystyle Z_{\mathbf {x} }^{(\ell )}} 、単位ベクトルx に対応する次数ℓ のゾーン調和は、[ 27 ]のように分解される。
さらに、ゾーン調和Z x ( ℓ ) ( y ) {\displaystyle Z_{\mathbf {x} }^{(\ell )}({\mathbf {y} })} 適切なゲーゲンバウアー多項式 の定数倍として与えられる。
( 2 )と( 3 )を組み合わせると、 x とyが 球面座標で表されている場合、次元n =2における( 1 )が得られる。最後に、 x = y で評価すると、関数恒等式が得られる 。薄暗い H ℓ ω n − 1 = ∑ j = 1 薄暗い ( 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}} ここで、ω n −1は ( n −1) 球の体積である。
収縮ルール もう1つの便利な恒等式は、2つの球面調和関数の積を球面調和関数の和として表します[ 28 ]。 Y 1 α ( θ 、 φ ) Y b β ( θ 、 φ ) = ( 2 1 + 1 ) ( 2 b + 1 ) 4 π ∑ c = 0 ∞ ∑ γ = − c c ( − 1 ) γ 2 c + 1 ( 1 b c α β − γ ) ( 1 b c 0 0 0 ) Y c γ ( θ 、 φ ) = ( 2 1 + 1 ) ( 2 b + 1 ) 4 π ∑ c = | 1 − b | 1 + b 1 2 c + 1 ⟨ 1 α b β | c α + β ⟩ ⟨ 1 0 b 0 | c 0 ⟩ Y c γ ( θ 、 φ ) {\displaystyle {\begin{aligned}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)\\&={\sqrt {\frac {\left(2a+1\right)\left(2b+1\right)}{4\pi }}}\sum _{c=|a-b|}^{a+b}{\frac {1}{\sqrt {2c+1}}}\left\langle a\,\alpha \,b\,\beta |c\,\alpha +\beta \right\rangle \left\langle a\,0\,b\,0|c\,0\right\rangle Y_{c}^{\gamma }\left(\theta ,\varphi \right)\end{aligned}}} この和の項の多くは自明にゼロです。c {\displaystyle c} そしてγ {\displaystyle \gamma } この和にゼロでない項をもたらすものは、 3j シンボル の選択規則によって決定されます。
クレブシュ・ゴルダン係数クレブシュ・ゴルダン係数は、2 つの球面調和関数の積を球面調和関数自体で展開した際に現れる係数です。本質的に同じ計算を行うためのさまざまな手法があり、ウィグナー3-jm 記号 、ラカー係数 、スレーター積分 などが含まれます。抽象的に言えば、クレブシュ・ゴルダン係数は、回転群 の2 つの既約表現の テンソル積 を既約表現の和として表します。適切に正規化すれば、係数は重複度になります。
球面調和関数の一覧 最初のいくつかの正規直交ラプラス球面調和関数の解析的表現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 π 罪 θ e − 私 φ Y 1 0 ( θ 、 φ ) = 1 2 3 π コス θ Y 1 1 ( θ 、 φ ) = − 1 2 3 2 π 罪 θ e 私 φ {\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 π 罪 2 θ e − 2 私 φ Y 2 − 1 ( θ 、 φ ) = 1 2 15 2 π 罪 θ コス θ e − 私 φ Y 2 0 ( θ 、 φ ) = 1 4 5 π ( 3 コス 2 θ − 1 ) Y 2 1 ( θ 、 φ ) = − 1 2 15 2 π 罪 θ コス θ e 私 φ Y 2 2 ( θ 、 φ ) = 1 4 15 2 π 罪 2 θ e 2 私 φ {\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}}}
高次元 古典的な球面調和関数は、単位球面上の複素数値関数として定義される。S 2 {\displaystyle S^{2}} 3次元ユークリッド空間内 R 3 {\displaystyle \mathbb {R} ^{3}} 球面調和関数は、より高次元のユークリッド空間に一般化することができる。R n {\displaystyle \mathbb {R} ^{n}} 以下のように、関数につながるS n − 1 → C {\displaystyle S^{n-1}\to \mathbb {C} } [ 31 ] Pℓ を n 個の実変数に関する次数ℓ の複素数値同次多項式 の空間とし、 ここでは関数として考える。R n → C {\displaystyle \mathbb {R} ^{n}\to \mathbb {C} } つまり、任意の実数に対して、多項式pが P ℓ に属するのは、次の条件が満たされる場合である。λ ∈ R {\displaystyle \lambda \in \mathbb {R} } 1つは
p ( λ x ) = λ ℓ p ( x ) 。 {\displaystyle p(\lambda \mathbf {x} )=\lambda ^{\ell }p(\mathbf {x} ).}
Aℓ を、 Pℓ のすべての 調和多項式 からなる部分空間とする。 A ℓ := { p ∈ P ℓ ∣ Δ p = 0 } 。 {\displaystyle \mathbf {A} _{\ell }:=\{p\in \mathbf {P} _{\ell }\,\mid \,\Delta p=0\}\,.} これらは(通常の)固体球面調和関数 である。Hℓを 単位球面上の関数の空間とする 。 S n − 1 := { x ∈ R n ∣ | x | = 1 } {\displaystyle S^{n-1}:=\{\mathbf {x} \in \mathbb {R} ^{n}\,\mid \,\left|x\right|=1\}} A ℓ からの制限によって得られたH ℓ := { f : S n − 1 → C ∣ 一部の人にとって p ∈ A ℓ 、 f ( x ) = p ( x ) すべての人々のために 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 ℓ の和は集合内で稠密 であるC ( S n − 1 ) {\displaystyle C(S^{n-1})} 連続関数のS n − 1 {\displaystyle S^{n-1}} 一様 位相に関して、ストーン・ワイエルシュトラス の定理により、これらの空間の和は球面上の二乗可積分関数の空間L 2 ( S n −1 )においても稠密である。したがって、球面上のすべての二乗可積分関数は、 L 2 の 意味で収束する球面調和関数の級数に一意に分解される。すべてのf ∈ H ℓ に対して、Δ S n − 1 f = − ℓ ( ℓ + n − 2 ) f 。 {\displaystyle \Delta _{S^{n-1}}f=-\ell (\ell +n-2)f.} ここで、Δ S n −1 はS n −1 上のラプラス・ベルトラミ演算子 である。この演算子は、3 次元のラプラシアンの角度部分のアナログである。すなわち、n 次元のラプラシアンは次のように分解される。∇ 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 ) の内積に関して直交することがわかる。つまり、∫ S n − 1 f g ¯ d Ω = 0 {\displaystyle \int _{S^{n-1}}f{\bar {g}}\,\mathrm {d} \Omega =0} f ∈ H ℓ およびg ∈ H k ( k ≠ ℓ の場合)について。逆に、空間H ℓ はΔ S n −1 の固有空間に正確に一致する。特に、スペクトル定理を リースポテンシャル に適用すると、Δ S n − 1 − 1 {\displaystyle \Delta _{S^{n-1}}^{-1}} 空間H ℓ がL 2 ( S n −1 ) において互いに直交し完備であることの別の証明を与える。 すべての同次多項式p ∈ P ℓ は一意に次の形式で記述できます[ 32 ] 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}}} ここでp j ∈ A j である。特に、薄暗い 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 = 罪 2 − n φ ∂ ∂ φ 罪 n − 2 φ ∂ ∂ φ + 罪 − 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}}} ここで、φは S n −1 上の球座標系における軸方向座標である。このような手順の最終結果は[ 33 ]である。 Y ℓ 1 、 … ℓ n − 1 ( θ 1 、 … θ n − 1 ) = 1 2 π e 私 ℓ 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})} ここで、添え字は| ℓ1 | ≤ ℓ2 ≤ ⋯ ≤ ℓn - 1 を 満たし、固有値は −ℓn-1 (ℓn-1 + n-2) である。積 に 含まれる 関数 は 、ルジャンドル 関数 を用いて定義さ れる 。 j P ¯ L ℓ ( θ ) = 2 L + j − 1 2 ( L + ℓ + j − 2 ) ! ( L − ℓ ) ! 罪 2 − j 2 ( θ ) P L + j − 2 2 − ( ℓ + j − 2 2 ) ( コス θ ) 。 {\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ℓにも作用する 。 ψ ↦ ψ ∘ ρ − 1 {\displaystyle \psi \mapsto \psi \circ \rho ^{-1}} ψ は球面調和関数、ρ は回転である。表現Hℓ は SO(3)の既約表現 である。 [ 34 ]
H ℓ の要素は、 A ℓ の要素の球面への制限として現れます。A ℓ は、3 次元ユークリッド空間R 3 上の次数ℓ の調和多項式です。ψ ∈ A ℓ の偏極 により、係数が存在します。 ψ 私 1 … 私 ℓ {\displaystyle \psi _{i_{1}\dots i_{\ell }}} インデックスに関して対称であり、要件によって一意に決定される ψ ( x 1 、 … 、 x n ) = ∑ 私 1 … 私 ℓ ψ 私 1 … 私 ℓ x 私 1 ⋯ x 私 ℓ 。 {\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 }}.} ψ が調和的である という条件は、テンソルが ψ 私 1 … 私 ℓ {\displaystyle \psi _{i_{1}\dots i_{\ell }}} 添え字のペアごとにトレースが ゼロでなければならない。したがって、 SO(3) の既約表現として、H ℓ は 次数ℓのトレース ゼロ対称テンソル の空間と同型である。
より一般的には、同様の記述は高次元でも成り立ちます。n次元 球面 上の球面調和関数の空間H ℓ は、トレースレス対称ℓテンソルに対応する SO( n +1) の既約表現です。しかし、SO(2) およびSO(3) のすべての既約テンソル表現はこのタイプであるのに対し、高次元の特殊直交群には、この方法では生じない追加の既約表現があります。
特殊直交群には、テンソル表現ではない追加のスピン表現が あり、通常は 球面調和関数ではありません。例外はSO(3)のスピン表現 です。厳密に言えば、これらはSO(3)の二重被覆SU(2)の表現です。SU(2)は単位 四元数 の群と同一視され、したがって3次元球面 と一致します。3次元球面上の球面調和関数の空間は、四元数乗算による作用に関して、SO(3)の特定のスピン表現です。
半球調和関数との関連性 球面調和関数は、2 つの関数セットに分けられます。[ 35 ] 1 つは半球調和関数 (HSH) で、半球上で直交かつ完全です。もう 1 つは相補半球調和関数 (CHSH) です。
注記 ↑ 3次元における球面調和関数のさまざまなアプローチに関する歴史的説明は、 MacRobert 1967の第 4 章に記載されています。「 ラプラス 球面調和関数」という用語は一般的に使用されています。Courant & Hilbert 1962 およびMeijer & Bauer 2004 を参照してください。 ↑ ここで採用されている球面調和関数へのアプローチは、( Courant & Hilbert 1962 、§V.8、§VII.5) に記載されている。 ↑ 物理的な応用では、無限遠でゼロになる解を採用することが多く、その結果A = 0 となります。これは球面調和関数の角度部分には影響しません。 ↑ Weisstein, Eric W. "Spherical Harmonic" . mathworld.wolfram.com . 2023-05-10 に取得. ↑ エドモンズ 1957 、§2.5 ↑ ホール 2013 セクション 17.6↑ Hall 2013 補題 17.161 2 3 Weisstein, Eric W. 「Condon-Shortley Phase」 . MathWorld—A Wolfram Web Resource . Wolfram Research . 2026-05-02 取得 。 1 2 3 「DLMF: §14.30 球面調和関数と回転楕円調和関数」 。NIST 数学関数デジタルライブラリ 。米国国立標準技術研究所。 2026年5月2日 取得 。 ↑ Heiskanen and Moritz, Physical Geodesy , 1967, eq. 1-62. ↑ ウィリアムズ、アール・G. (1999). フーリエ音響学:音響放射と近接場音響ホログラフィー . サンディエゴ、カリフォルニア州:アカデミック・プレス. ISBN 0-08-050690-9 OCLC 181010993 ↑ メサイア、アルバート( 1999). 量子力学:2巻を1冊にまとめたもの (2巻を1冊にまとめた完全復刻 版)。ニューヨーク州ミネオラ:ドーバー出版。520–523 頁 。ISBN 0-486-40924-4 。↑ クロード・コーエン=タンヌージ;バーナード・ディウ;フランク・ラロエ (1996)。 量子力学 。スーザン・リード・ヘムリー訳。他 。ワイリー・インターサイエンス: ワイリー。 ISBN 978-0-471-56952-7 。1 2 ブレイクリー、リチャード (1995)。 重力と磁気応用におけるポテンシャル理論 。ケンブリッジ、イギリス、 ニューヨーク:ケンブリッジ大学出版局。p . 113。ISBN 978-0-521-41508-8 。↑ 引用エラー: 指定された名前の参照SHTOOLS complexが呼び出されましたが、定義されていません (ヘルプ ページを参照してください)。 ↑ Whittaker & Watson 1927 、p. 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 ↑ ジャクソン、ジョン・デヴィッド (1999)。 古典電気力学 (第 3 版、[Nachdr.] 版)。ニューヨーク州ホーボーケン:ワイリー。 p. 108.ISBN 978-0-471-30932-1 。↑ Dahlen, FA; Simons, Frederik J. (2008). "地球物理学および宇宙論における球面上のスペクトル推定". Geophysical Journal International . 174 (3): 774– 807. arXiv : 0705.3083 . doi : 10.1111/j.1365-246X.2008.03854.x . ↑ 引用エラー: 指定された名前の参照SteinWeissが呼び出されましたが、定義されていません (ヘルプ ページを参照してください)。 ↑ 引用エラー: 指定された名前の参照DaiXuが呼び出されましたが、定義されていません (ヘルプ ページを参照してください)。 ↑ Edmonds, AR (1996). Angular Momentum In Quantum Mechanics . Princeton University Press. p. 63 . ↑ これは、次数ℓ の球面調和関数の任意の正規直交基底に対して有効です。単位べき乗調和関数の場合は、 4π の係数を取り除く必要があります。 ↑ Whittaker & Watson 1927 、p. 395 ↑ 未販売 1927 ↑ スタイン& ワイス 1971 、§IV.2 ↑ Brink, DM; Satchler, GR Angular Momentum . Oxford University Press. p. 146. ↑ 「球面振動 – 球面調和関数」 。 ↑ エレメンコ、ヤコブソン、 ナディラシビリ 2007 ↑ Solomentsev 2001 ; Stein & Weiss 1971 、§Iv.2 ↑ 参照: Axler, Sheldon; Ramey, Wade (1995)『 Harmonic Polynomials and Dirichlet-Type Problems』の系1.8 ↑ 樋口篤(1987) 「N球面上の対称テンソル球面調和関数とド・ジッター群SO(N,1)への応用」 . Journal of Mathematical Physics . 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.、第22巻、(1968年)。 ↑ JD タルマン、『特殊関数、群論的アプローチ』( EP ウィグナー の講義に基づく)、WA ベンジャミン、ニューヨーク(1968年)。 ↑ W. ミラー、『対称性と変数分離』、 アディソン・ウェスリー、レディング(1977年)。 ↑ A. Wawrzyńczyk、『群表現と特殊関数』 、ポーランド科学出版社、ワルシャワ(1984年)。
参考文献
引用文献 クーラント、リチャード ;ヒルベルト、デイヴィッド (1962)、『数理物理学の方法、第1巻』 、ワイリー・インターサイエンス 。Edmonds, AR (1957)、『量子力学における角運動量』 、プリンストン大学出版局、Bibcode : 1957amqm.book.....E、ISBN 0-691-07912-9 Eremenko, Alexandre; Jakobson, Dmitry; Nadirashvili, Nikolai (2007)、「S² および R² 上の節点集合と節点領域について」、Annales de l'Institut Fourier 、57 (7): 2345–2360 、doi : 10.5802/aif.2335 、ISSN 0373-0956、MR 2394544 Hall, Brian C. (2013), 『数学者のための量子論』 、大学院数学テキスト、第 267巻、Springer、Bibcode : 2013qtm..book.....H、ISBN 978-1-4614-7115-8 MacRobert, TM (1967), Spherical harmonics: An elementary treatise on harmonic functions, with applications , Pergamon Press 。Meijer, Paul Herman Ernst; Bauer, Edmond (2004), Group theory: The application to quantum mechanics , Dover, ISBN 978-0-486-43798-9 。Solomentsev, ED (2001) [1994]、「球面調和関数」、数学百科事典 、EMS Press 。Stein, Elias ; Weiss, Guido (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)『現代解析学講義 』ケンブリッジ大学出版局 、 392ページ 。
一般的な参考文献 EW ホブソン著『球面調和関数と楕円体調和関数の理論』 (1955年)、チェルシー出版、ISBN 978-0-8284-0104-3 。 C. ミュラー著『球面調和関数 』(1966年)、シュプリンガー社、数学講義録、第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 版)、ニューヨーク:ケンブリッジ大学出版局、ISBN 978-0-521-88068-8 DA Varshalovich、AN Moskalev、VK Khersonskii角運動量の量子理論 、(1988) World Scientific Publishing Co.、シンガポール、ISBN 9971-5-0107-4 ワイススタイン、エリック・W. 「球面調和関数」 . MathWorld .マドック、ジョン、「Boost.Mathにおける球面調和関数」