球面調和関数とは
この計算機は、複素球面調和関数 \(Y_n^m(\theta,\phi)\) を計算します。これはラプラス方程式の解の角度成分にあたる関数です。球面調和関数は物理学や応用数学のあらゆる分野に登場します。たとえば量子力学(原子軌道)、電磁気学、測地学、コンピュータグラフィックスの照明計算、地震学などです。本ツールでは、指定した天頂角(極角)θと方位角φにおける Y の実部・虚部・絶対値を求めます。内容は純粋な数学であり、国や地域を問わず同じ結果が得られます。
使い方
まず関数の定義を選びます。タイプAは量子力学で用いられる完全規格化の正規直交(コンドン=ショートレー位相)方式で、球面全体にわたる \(|Y|^2\) の積分が 1 になります。タイプBは非規格化の方式で、\(P_n^m(\cos\theta)\) に方位角の位相を掛けただけのシンプルな形です。次数 n(0, 1, 2, …)、位数 m(\(-n \le m \le n\))、天頂角 θ(度)、方位角 φ(度)を入力し、計算を実行すると実部と虚部が表示されます。
計算式の解説
$$Y_n^m(\theta,\phi) = N_{n,m}\cdot P_n^m(\cos\theta)\cdot e^{i\,m\phi}$$ここで \(e^{i\,m\phi} = \cos(m\phi) + i\sin(m\phi)\) です。ルジャンドル陪関数 \(P_n^m\) にはコンドン=ショートレー位相 \((-1)^m\) が含まれます。タイプAでは \(N = \sqrt{\frac{2n+1}{4\pi}\cdot\frac{(n-m)!}{(n+m)!}}\)、タイプBでは \(N = 1\) です。度で入力した角度は、計算前にラジアンへ変換されます。
計算例
タイプA、\(n = 2\)、\(m = 1\)、\(\theta = 30°\)、\(\phi = 30°\) の場合を見てみましょう。$$x = \cos 30° = 0.8660254$$ $$P_2^1(x) = -\sqrt{1-x^2}\cdot 3x = -1.2990381$$ $$N = \sqrt{\frac{5}{4\pi}\cdot\frac{1}{6}} = 0.2575162$$ となります。したがって \(N\cdot P = -0.3345283\)。\(\cos 30° = 0.8660254\)、\(\sin 30° = 0.5\) を用いると、実部は \(-0.2897113\)、虚部は \(-0.1672642\)、絶対値は \(0.3345283\) となります。
よくある質問
m の有効範囲は? 位数 m は \(-n \le m \le n\) を満たす整数でなければなりません。この範囲を外れると球面調和関数は定義されません。
なぜ極で Y がゼロになるの? \(\theta = 0°\) または \(180°\) では \(\sqrt{1-x^2} = 0\) となるため、\(m \ne 0\) のすべての場合で \(P_n^m = 0\) になります。有限の値を保つのは \(m = 0\) のときだけです。
どの符号規約を使っているの? コンドン=ショートレー位相 \((-1)^m\) を含めており、物理学で標準的に用いられる規約に一致します。