理論物理学

球面調和関数

理論物理学

概要

量子力学で3次元空間における軌道角運動量演算子 $$\hbar\bm{L}=\bm{x}\times\bm{p}$$ に対し、$\bm{L}^2$ と $L_z$ の同時固有ケットを $\ket{l,m}$ とします。

球座標を $(r,\theta,\phi)$ で表し、同時固有ケットを $$ Y_{lm}(\theta,\phi) = \braket{\theta,\phi|l,m} $$ のように座標表示したものが、一般的に使われている球面調和関数です。

球面調和関数の具体的な表式は、 $$ Y_{lm}(\theta,\phi) = \frac{1}{\sqrt{2\pi}}\em^{im\phi}\sqrt{\frac{2l+1}{2}\frac{(l-m)!}{(l+m)!}}P_l^m(\cos\theta) $$ と与えられます。ここで $P_l^m(x)$ はLegendre陪関数です。
本記事では、Legendre陪関数の諸性質を用いた解法と、Legendre陪関数を知らない前提での解法の両方を実演します。

球面調和関数の表式には、いくつかの等価な書き換えが存在するので、どの流派を用いるのか注意が必要です。

解説

球座標表示

導入として、量子力学の角運動量演算子について、軽く触れておきましょう。

3次元空間における軌道角運動量演算子を $\hbar\bm{L}=\bm{x}\times\bm{p}$ とします。
具体的に3成分を具体的に書くと $$ \begin{align} L_x &= -i\left( y\pdv{z}-z\pdv{y} \right) \nonumber \\ L_y &= -i\left( z\pdv{x}-x\pdv{z} \right) \label{eq:Lxyz} \\ L_z &= -i\left( x\pdv{y}-y\pdv{x} \right) \nonumber \end{align} $$ となります。
正準交換関係から、$\bm{L}^2,\,L_3$ が同時対角化可能な演算子の組み合わせであることが分かります。同時固有ケットを $\ket{l,m}$ で表すと
$$ \begin{equation} \bm{L}^2\ket{l,m} = l(l+1)\ket{l,m} \;,\;\; L_z\ket{l,m} = m\ket{l,m} \label{eq:opeigen} \end{equation} $$ となりますが、固有値は $l =0,1,2,\cdots$ と $m=-l,-l+1,\cdots,l-1,l$ と量子化されます。
一般の角運動量代数では半整数の $l$ も許されますが、軌道角運動量なので空間の一価性により $l$ は整数に限られます。

ここまでの角運動量演算子の話は、以下の記事で詳しく説明しています:

球座標を次式で定義します。定義域は $r\geq0,\, 0\leq\theta\leq\pi,\, 0\leq\phi<2\pi$ です。 $$ \begin{equation} \bm{x} = \begin{pmatrix} x \\ y \\ z \end{pmatrix} = \begin{pmatrix} r\sin\theta\cos\phi \\ r\sin\theta\sin\phi \\ r\cos\theta \end{pmatrix} \label{eq:scoordinate} \end{equation} $$

球座標を対角的にする基底で、上記の演算子とケットを座標表示していきます。
角運動量演算子の $x,y,z$ 成分は、以下のようになります。

問1-1

$$ \begin{align} L_x &= i\left( \sin\phi\pdv{\theta}+\cot\theta\cos\phi\pdv{\phi} \right) \nonumber \\ L_y &= i\left( -\cos\phi\pdv{\theta}+\cot\theta\sin\phi\pdv{\phi} \right) \label{eq:Ls} \\ L_z &= -i\pdv{\phi} \nonumber \end{align} $$

※ $\cot\theta=1/\tan\theta$

直交座標 $x,y,z$ の微分を球座標に変換する。式\eqref{eq:scoordinate}に連鎖律 $$ \pdv{x_i} = \pdv{r}{x_i}\pdv{r} + \pdv{\theta}{x_i}\pdv{\theta} + \pdv{\phi}{x_i}\pdv{\phi} $$ を適用すると $$ \begin{align*} \pdv{x} &= \sin\theta\cos\phi\pdv{r}+\frac{1}{r}\cos\theta\cos\phi\pdv{\theta}-\frac{\sin\phi}{r\sin\theta}\pdv{\phi} \\ \pdv{y} &= \sin\theta\sin\phi\pdv{r}+\frac{1}{r}\cos\theta\sin\phi\pdv{\theta}+\frac{\cos\phi}{r\sin\theta}\pdv{\phi} \\ \pdv{z} &= \cos\theta\pdv{r}-\frac{1}{r}\sin\theta\pdv{\theta} \end{align*} $$ 式\eqref{eq:Lxyz}の直交座標 $x,y,z$ とその微分に、式\eqref{eq:scoordinate}と上式を代入する。

なお、ブラを略さずに正確に書くと $$ \bra{r, \theta,\phi}L_z = -i\pdv{\phi}\bra{r, \theta,\phi} $$ 等になります。

式\eqref{eq:Ls}より、角運動量の大きさ $\bm{L}^2=L_x^2+L_y^2+L_z^2$ および $L_z$ の昇降演算子 $L_\pm = L_x\pm iL_y$ も、球座標で表せるようになります。

問1-2

$$ \begin{equation} \bm{L}^2 = -\left[ \frac{1}{\sin\theta}\pdv{\theta}\left(\sin\theta\pdv{\theta}\right)+\frac{1}{\sin^2\theta}\pdv[2]{\phi} \right] \label{eq:L2s} \end{equation} $$

$$ \begin{align} L_+ &= \em^{i\phi}\left( \pdv{\theta}+i\cot\theta\pdv{\phi} \right) \label{eq:Lps} \\ L_- &= \em^{-i\phi}\left( -\pdv{\theta}+i\cot\theta\pdv{\phi} \right) \label{eq:Lms} \\ \end{align} $$

演算子なので右側に任意の関数 $F(\bm{x})$ があると考えて計算するとよい。 $$ \begin{align*} L_x^2F(\bm{x}) &= i^2\left( \sin\phi\pdv{\theta}+\cot\theta\cos\phi\pdv{\phi} \right)^{\!2}F \\ &= -\sin\phi\left( \sin\phi\pdv[2]{F}{\theta}-\frac{\cos\phi}{\sin^2\theta}\pdv{F}{\phi}+\cot\theta\cos\phi\pdv{F}{\theta}{\phi} \right) \\ &\qquad -\cot\theta\cos\phi\left( \cos\phi\pdv{F}{\theta}+\sin\phi\pdv{F}{\theta}{\phi}-\cot\theta\sin\phi\pdv{F}{\phi}+\cot\theta\cos\phi\pdv[2]{F}{\phi} \right) \end{align*} $$ $$ \begin{align*} L_y^2F(\bm{x}) &= i^2\left( -\cos\phi\pdv{\theta}+\cot\theta\sin\phi\pdv{\phi} \right)^{\!2}F \\ &= -\cos\phi\left( \cos\phi\pdv[2]{F}{\theta}+\frac{\sin\phi}{\sin^2\theta}\pdv{F}{\phi}-\cot\theta\sin\phi\pdv{F}{\theta}{\phi} \right) \\ &\qquad -\cot\theta\sin\phi\left( \sin\phi\pdv{F}{\theta}-\cos\phi\pdv{F}{\theta}{\phi}+\cot\theta\cos\phi\pdv{F}{\phi}+\cot\theta\sin\phi\pdv[2]{F}{\phi} \right) \end{align*} $$ $$ L_z^2F(\bm{x}) = (-i)^2\left( \pdv{\phi} \right)^{\!2}F = -\pdv[2]{F}{\phi} $$ 以上を足すと $$ \begin{align*} \bm{L}^2F(\bm{x}) &= -\left[ \pdv[2]{F}{\theta}+\cot\theta\pdv{F}{\theta}+\cot^2\theta\pdv[2]{F}{\phi}+\pdv[2]{F}{\phi} \right] \\ &= -\left[ \frac{1}{\sin\theta}\pdv{\theta}\left(\sin\theta\pdv{\theta}\right)+\frac{1}{\sin^2\theta}\pdv[2]{\phi} \right]F \end{align*} $$

$L_\pm$ もまた定義通りに計算する。 $$ \begin{align*} L_\pm &= L_x \pm iL_y \\ &= i\left( \sin\phi\pdv{\theta}+\cot\theta\cos\phi\pdv{\phi} \right) \pm i\cdot i\left( -\cos\phi\pdv{\theta}+\cot\theta\sin\phi\pdv{\phi} \right) \\ &= \pm(\cos\phi\pm i\sin\phi)\pdv{\theta} +i\cot\theta(\cos\phi\pm i\sin\phi)\pdv{\phi} \\ &= \em^{\pm i\phi}\left( \pm\pdv{\theta}+i\cot\theta\pdv{\phi} \right) \end{align*} $$

式\eqref{eq:Lps}, \eqref{eq:Lms}を一見すると、$(L_\pm)^\dagger=L_\mp$ が成立していないように思われます。
昇降演算子のエルミート共役は単純に複素共役をとるだけでは不十分で、内積をきちんと計算すれば $(L_\pm)^\dagger=L_\mp$ が成立していることを確かめられます。

問1-3

発散しない任意の関数 $f(\theta,\phi), g(\theta,\phi)$ に対し、内積を $$ \braket{f|g} = \int_0^{2\pi}\!\mathrm{d}\phi\int_0^\pi\!\mathrm{d}\theta\sin\theta\, f^\ast (\theta,\phi)g(\theta,\phi) $$ と定義するとき、以下が成り立つ。($\braket{f|L_-g}$ についても同様) $$ \braket{f|L_+g} = \braket{L_-f|g} $$

$$ \begin{align*} \braket{f|L_+g} &= \int_0^{2\pi}\!\mathrm{d}\phi \int_0^\pi\!\mathrm{d}\theta\sin\theta\, f^\ast\cdot \em^{i\phi}\left( \pdv{g}{\theta}+i\cot\theta\pdv{g}{\phi} \right) \\ &= \int_0^{2\pi}\!\mathrm{d}\phi\, \em^{i\phi}\left( \Bigl[ \sin\theta f^\ast g \Bigr]_{\theta=0}^{\theta=\pi}- \int_0^\pi\!\mathrm{d}\theta\, \pdv{(\sin\theta f^\ast)}{\theta}g \right) \\ &\qquad +\int_0^\pi\!\mathrm{d}\theta\sin\theta\, i\cot\theta\left( \Bigl[ \em^{i\phi} f^\ast g \Bigr]_{\phi=0}^{\phi=2\pi}- \int_0^{2\pi}\!\mathrm{d}\phi\, \pdv{(\em^{i\phi}f^\ast)}{\phi}g \right) \end{align*} $$ 表面項は、$\sin0=\sin\pi=0$ および $\phi=0,2\pi$ の周期境界条件のため消える。 $$ \begin{align*} \braket{f|L_+g} &= \int_0^{2\pi}\!\mathrm{d}\phi \int_0^\pi\!\mathrm{d}\theta\, \left\{ -\em^{i\phi}\left( \cancel{\cos\theta f^\ast} +\sin\theta\pdv{f^\ast}{\theta} \right)g \right. \\ &\qquad\qquad \left. -i\cos\theta\left( \cancel{i\em^{i\phi} f^\ast} +\em^{i\phi}\pdv{f^\ast}{\phi} \right)g \right\} \\ &= \int_0^{2\pi}\!\mathrm{d}\phi \int_0^\pi\!\mathrm{d}\theta\sin\theta\, \em^{i\phi}\left( -\pdv{f^\ast}{\theta}-i\cot\theta\pdv{f^\ast}{\phi} \right) g \\ &= \int_0^{2\pi}\!\mathrm{d}\phi \int_0^\pi\!\mathrm{d}\theta\sin\theta\, \left\{ \em^{-i\phi}\left( -\pdv{\theta}+i\cot\theta\pdv{\phi} \right)f \right\}^\ast\cdot g \\ &= \braket{L_-f|g} \end{align*} $$

そして同時固有ケットを $$ Y_{lm}(\theta,\phi) = \braket{\theta,\phi|l,m} $$ のように座標表示したものを、球面調和関数と定義します。

式\eqref{eq:opeigen}を式\eqref{eq:Ls}, \eqref{eq:L2s}で座標表示することにより、$Y_{lm}(\theta,\phi)$ が満たすべき固有値方程式は以下の2つとなります。 $$ \begin{align} -\left[ \frac{1}{\sin\theta}\pdv{\theta}\left(\sin\theta\pdv{\theta}\right)+\frac{1}{\sin^2\theta}\pdv[2]{\phi} \right] Y_{lm}(\theta,\phi) &= l(l+1)Y_{lm}(\theta,\phi) \label{eq:L2eigeneq} \\ -i\pdv{\phi}Y_{lm}(\theta,\phi) &= mY_{lm}(\theta,\phi) \label{eq:Lzeigeneq} \end{align} $$ 演算子 $\bm{L}^2, L_z$ は $\theta,\phi$ のみの関数であり、$r$ については何も制約を与えません。そのため、球面調和関数に $r$ 依存性を含めることはせず、球座標のブラも $\bra{\theta,\phi}$ と記載しています。

本記事では量子力学の文脈で、球面調和関数を導入しました。
しかし本来は量子力学の仮定を用いることなく、球座標のラプラシアンから直接定義される概念です。

3次元空間のLaplace方程式を考えます。 $$ \triangle F(\bm{x}) = \left[ \pdv[2]{x}+\pdv[2]{y}+\pdv[2]{z} \right]F(\bm{x}) =0 $$ ラプラシアンを球座標で書くと、以下のようになります。 $$ \left[ \frac{1}{r}\pdv[2]{r}r+\frac{1}{r^2}\left\{ \frac{1}{\sin\theta}\pdv{\theta}\left(\sin\theta\pdv{\theta}\right)+\frac{1}{\sin^2\theta}\pdv[2]{\phi}\right\} \right]F(r,\theta,\phi) =0 $$ 解を $F(r,\theta,\phi)=R(r)Y(\theta,\phi)$ と変数分離形で書き表して、上式の $\{\cdots\}Y(\theta,\phi)$ に関する固有値方程式を調べることで、式\eqref{eq:L2eigeneq}, \eqref{eq:Lzeigeneq}に帰着します。

$\phi$ 部分の固有関数

$\phi$ に関する微分方程式\eqref{eq:Lzeigeneq}の解は、よく知られているように指数関数となります。 $$ Y_{lm}(\theta,\phi) \propto \em^{im\phi} $$ $m$ は整数なので、$Y_{lm}(\theta,\phi+2\pi)=Y_{lm}(\theta,\phi)$ が保証されます。

球面調和関数は、$\theta$ と $\phi$ の変数分離形となります。
$\phi$ の部分だけで正規直交性を満たすように、規格化係数を取り出して定義しておきます。 $$ \begin{equation} Y_{lm}(\theta,\phi) = \frac{1}{\sqrt{2\pi}}\Theta_{lm}(\theta)\em^{im\phi} \label{eq:YTheta} \end{equation} $$

問2

$$ \begin{equation} \int_0^{2\pi} \!\mathrm{d}\phi \left(\frac{1}{\sqrt{2\pi}}\em^{im\phi}\right)^{\!\!\ast} \cdot \left(\frac{1}{\sqrt{2\pi}}\em^{im’\phi}\right) = \delta_{mm’} \label{eq:m_norm} \end{equation} $$

$$ \begin{align*} & \int_0^{2\pi} \!\mathrm{d}\phi \left(\frac{1}{\sqrt{2\pi}}\em^{im\phi}\right)^{\!\!\ast} \cdot \left(\frac{1}{\sqrt{2\pi}}\em^{im’\phi}\right) \\ &= \frac{1}{2\pi}\int_0^{2\pi} \!\mathrm{d}\phi\; \em^{i(m’-m)\phi} \\ &= \left\{ \begin{array}{lc} \displaystyle\frac{1}{2\pi}\int_0^{2\pi} \!\mathrm{d}\phi\; 1 & (m=m’) \\ \displaystyle\frac{1}{2\pi}\left[ \frac{1}{i(m’-m)}\em^{i(m’-m)\phi} \right]_0^{2\pi} & (m\neq m’) \end{array} \right. \\ &= \delta_{mm’} \end{align*} $$

$\theta$ 部分の固有関数(Legendre陪関数の方法)

この節では、Legendre陪関数の性質を既知として、$\theta$ 部分の固有関数を導出していきます
以下の公式を証明無しに使います。

Legendre陪関数 $P_l^m(x)$ は、Legendreの陪微分方程式 $$ \begin{equation} \left[\dv{x}(1-x^2)\dv{x}+l(l+1)-\frac{m^2}{1-x^2}\right]P_l^m(x)=0 \label{eq:Legendre_diffeq} \end{equation} $$ の解であり、$l=0,1,2,\cdots,\,m=-l,-l+1,\cdots,l-1,l$ のとき $x=\pm1$ で正則である。
具体的な表式は、Legendre多項式 $P_l(x)$ のRodriguesの公式を用いて、以下で与えられる。 $$ \begin{align} P_l^m(x) &= (-1)^m(1-x^2)^\frac{m}{2}\dv[m]{x}P_l(x) \nonumber \\ &= \frac{(-1)^m}{2^l l!}(1-x^2)^\frac{m}{2}\dv[l+m]{x}(x^2-1)^l \label{eq:rodrigues} \end{align} $$ 公式 $$ \begin{equation} P_l^{-m}(x) = (-1)^m\frac{(l-m)!}{(l+m)!}P_l^m(x) \label{eq:Legendre_-m} \end{equation} $$ $$ \begin{equation} \int_{-1}^{1}\!\mathrm{d}x\,P_l^m(x)P_{l’}^{m}(x)=\frac{2}{2l+1}\frac{(l+m)!}{(l-m)!}\delta_{ll’} \label{eq:Legendre_orth} \end{equation} $$ が成り立つ。

ハーポルホード高橋
ハーポルホード高橋

$Y,\Theta,P$ について、 添字 $l,m$ を上付きにするか下付きにするかは、本記事では特にこだわっていません。
位相因子を付けるかどうかで使い分けている流派もあるようです。

$\theta$ に関する微分方程式\eqref{eq:L2eigeneq}は、$\phi$ の部分の解を代入すると、以下のようになります。 $$ \begin{equation} -\left[ \frac{1}{\sin\theta}\dv{\theta}\left(\sin\theta\dv{\theta}\right)-\frac{m^2}{\sin^2\theta} \right]\Theta_{lm}(\theta) = l(l+1)\Theta_{lm}(\theta) \label{eq:Thetadiffeq} \end{equation} $$

$x=\cos\theta$ に置換します。極角の定義域は$0\leq\theta\leq\pi$なので、$-1\leq x\leq1$とは1対1で対応します。 $$ \begin{equation} \dv{x} = \dv{\theta}{(\cos\theta)}\dv{\theta} = -\frac{1}{\sin\theta}\dv{\theta} \label{eq:dvx_dvcos} \end{equation} $$ を用いて式\eqref{eq:Thetadiffeq}を変形すると、 $$ \begin{align*} -\left[ \frac{-1}{\sin\theta}\dv{\theta}\left((1-\cos^2\theta)\frac{-1}{\sin\theta}\dv{\theta}\right)+l(l+1)-\frac{m^2}{1-\cos^2\theta} \right]\Theta_{lm}(\theta) &= 0 \\ \therefore \left[ \dv{x}\left((1-x^2)\dv{x}\right)+l(l+1)-\frac{m^2}{1-x^2} \right]\left.\Theta_{lm}(\theta)\right|_{\theta\,\mapsto\,x} &= 0 \end{align*} $$ これはLegendreの陪微分方程式\eqref{eq:Legendre_diffeq}と一致しており、解はLegendre陪関数 $P_l^m(x)$ で与えられます。
意味のある解になるためには波動関数が $\theta=0,\pi$ で発散してはいけませんが、正則である条件は角運動量演算子の性質から導かれた $l,m$ の制約と矛盾していません。 $$ \Theta_{lm}(\theta) \propto P_l^m(\cos\theta) $$ なお、微分方程式は $m\to -m$ で不変なので $P_l^{-m}(\cos\theta)$ も解ですが、式\eqref{eq:Legendre_-m}の性質から、$P_l^m(\cos\theta)$ と独立な解ではありません。

$\Theta_{lm}(\theta)$ だけについても正規直交性を満たすように、規格化係数を付けておきます。 $$ \Theta_{lm}(\theta) = \sqrt{\frac{2l+1}{2}\frac{(l-m)!}{(l+m)!}}P_l^m(\cos\theta) $$

問3

$$ \begin{equation} \int_0^{\pi} \!\sin\theta \mathrm{d}\theta\; \Theta_{lm}(\theta)^\ast \cdot \Theta_{l’m}(\theta) = \delta_{ll’} \label{eq:l_norm} \end{equation} $$

$\cos\theta=x$ に置換して、Legendre陪関数の直交性\eqref{eq:Legendre_orth}を用いる。 $$ \begin{align*} & \int_0^{\pi} \!\sin\theta \mathrm{d}\theta\; \Theta_{lm}(\theta)^\ast \cdot \Theta_{l’m}(\theta) \\ &= \sqrt{\frac{(2l+1)(2l’+1)}{2^2}\frac{(l-m)!(l’-m)!}{(l+m)!(l’+m)!}} \int_0^{\pi} \!\sin\theta \mathrm{d}\theta P_l^m(\cos\theta)P_{l’}^m(\cos\theta) \\ &= \sqrt{\frac{(2l+1)(2l’+1)}{2^2}\frac{(l-m)!(l’-m)!}{(l+m)!(l’+m)!}} \int_{-1}^{1} \!\mathrm{d}x P_l^m(x)P_{l’}^m(x) \\ &= \sqrt{\frac{(2l+1)(2l’+1)}{2^2}\frac{(l-m)!(l’-m)!}{(l+m)!(l’+m)!}} \cdot \frac{2}{2l+1}\frac{(l+m)!}{(l-m)!}\delta_{ll’} \\ &= \delta_{ll’} \end{align*} $$

まとめると、球面調和関数は $$ \begin{equation} Y_{lm}(\theta,\phi) = \frac{1}{\sqrt{2\pi}}\em^{im\phi}\sqrt{\frac{2l+1}{2}\frac{(l-m)!}{(l+m)!}}P_l^m(\cos\theta) \label{eq:Y_main} \end{equation} $$ となります。式\eqref{eq:m_norm}, \eqref{eq:l_norm}より、立体角積分 $\int\mathrm{d}\Omega = \int\sin\theta\mathrm{d}\theta\mathrm{d}\phi$ に関して正規直交性を満たします。 $$ \int\mathrm{d}\Omega\,Y_{lm}^\ast(\theta,\phi)Y_{l’m’}(\theta,\phi) = \delta_{ll’}\delta_{mm’} $$

$\theta$ 部分の固有関数(昇降演算子の方法)

今度は、Legendre陪関数の知識を使わずに $\Theta_{lm}(\theta)$ を導出することを試みます
そのためには、角運動量演算子の一般論で証明した、固有値 $m$ の昇降演算子に注目します。

昇降演算子 $L_\pm$ について $$ \begin{align*} L_+\ket{l,m} = \sqrt{(l-m)(l+m+1)}\ket{l,m+1} \\ L_-\ket{l,m} = \sqrt{(l+m)(l-m+1)}\ket{l,m-1} \end{align*} $$ 特に $$ L_+\ket{l,l} = L_-\ket{l,-l} = 0 $$ が成り立つ。

導出は以下の記事をご参照ください。

これらの式は、左から球座標のブラ $\bra{\theta,\phi}$ を作用させることで、球面調和関数 $Y_{lm}(\theta,\phi)$ に対する方程式になります。
方針としては、$m$ が最大の状態 $Y_{ll}(\theta,\phi)$ に対して下降演算子 $L_-$ を作用させることで、任意の $Y_{lm}(\theta,\phi)$ が得られることを利用します。使う式は以下の2つです。 $$ L_+ Y_{ll}(\theta,\phi) = 0 $$ $$ L_- Y_{lm}(\theta,\phi) = \sqrt{(l+m)(l-m+1)}Y_{l,m-1}(\theta,\phi) $$

まずは $\Theta_{ll}(\theta)$ の具体形を求めます。
$L_+$ を式\eqref{eq:Lps}で座標表示に書き直すと、以下の微分方程式が得られます。

問4

$$ \left( \dv{\theta}-l\frac{\cos\theta}{\sin\theta} \right)\Theta_{ll}(\theta) = 0 $$

式\eqref{eq:Lps}より $$ \em^{i\phi}\left( \pdv{\theta}+i\cot\theta\pdv{\phi} \right)Y_{ll}(\theta) = 0 $$ $\phi$ の部分は式\eqref{eq:YTheta}で既に得られているので、代入すると $$ \begin{align*} \em^{i\phi}\left( \pdv{\theta}+i\frac{\cos\theta}{\sin\theta}\pdv{\phi} \right)\frac{1}{\sqrt{2\pi}}\Theta_{ll}(\theta)\em^{il\phi} &= 0 \\ \therefore\;\; \frac{1}{\sqrt{2\pi}}\em^{i(l+1)\phi} \left( \dv{\theta}-l\frac{\cos\theta}{\sin\theta} \right)\Theta_{ll}(\theta) &= 0 \end{align*} $$

よく観察すると、微分方程式の解は $$ \Theta_{ll}(\theta) \propto \sin^l\theta $$ であることが分かります。規格化条件から係数を決めると、 $$ \begin{equation} \Theta_{ll}(\theta) = \frac{1}{2^ll!}\sqrt{\frac{(2l+1)!}{2}} \sin^l\theta \label{eq:Theta_ll} \end{equation} $$ となります。

問5

式\eqref{eq:Theta_ll}の $\Theta_{ll}(\theta)$ は規格化$$ \int_0^\pi \!\sin\theta\mathrm{d}\theta\; |\Theta_{ll}(\theta)|^2 = 1 $$を満たす。

$u=\cos\theta$ とおくと、$\dv{u}{\theta}=-\sin\theta$ より $$ \begin{align*} \int_0^\pi \!\mathrm{d}\theta\, \sin^{2l+1}\theta &= \int_0^\pi \!\mathrm{d}\theta\, \sin\theta(1-\cos^2\theta)^l \\ &= \int_{-1}^1\!\mathrm{d}u\, (1-u^2)^l \\ &= \int_{-1}^1\!\mathrm{d}u\, (1-u)^l(1+u)^l \\ &= \int_{-1}^1\!\mathrm{d}u\, (1-u)^l \cdot \frac{l!}{(2l)!}\dv[l]{u}(1+u)^{2l} \\ &= \frac{l!}{(2l)!} \left( \left[ (1-u)^l\dv[l-1]{u}(1+u)^{2l} \right]_{-1}^1 + l\int_{-1}^1\!\mathrm{d}u\, (1-u)^{l-1} \cdot \dv[l-1]{u}(1+u)^{2l} \right) \\ &= \cdots \\ &= 0 + \frac{l!}{(2l)!}\cdot l! \int_{-1}^1\!\mathrm{d}u\, (1+u)^{2l} \\ &= \frac{(l!)^2}{(2l)!}\left[ \frac{1}{2l+1}(1+u)^{2l+1} \right]_{-1}^1 \\ &= \frac{(l!)^2}{(2l)!}\frac{2^{2l+1}}{2l+1} \\ &= (2^ll!)^2\frac{2}{(2l+1)!} \end{align*} $$ よって $$ \begin{align*} \int_0^\pi \!\sin\theta\mathrm{d}\theta\, \left|\frac{1}{2^ll!}\sqrt{\frac{(2l+1)!}{2}} \sin^l\theta \right|^2 &= \frac{1}{(2^ll!)^2}\frac{(2l+1)!}{2} \int_0^\pi \!\mathrm{d}\theta\, \sin^{2l+1}\theta \\ &= \frac{1}{(2^ll!)^2}\frac{(2l+1)!}{2} \cdot (2^ll!)^2\frac{2}{(2l+1)!} \\ &= 1 \end{align*} $$

次に、下降演算子の作用を考えます。
$L_-$ を式\eqref{eq:Lms}で座標表示に書き直すと、$\Theta_{lm}(\theta)$ の漸化式が得られます。

問6

$$ \sin^m\theta\cdot\Theta_{l,m-1}(\theta) = \gamma_{lm}\dv{(\cos\theta)}\Bigl( \sin^m\theta\cdot\Theta_{lm}(\theta) \Bigr) $$ ただし係数を $\gamma_{lm}$ を以下で定める。 $$ \gamma_{lm} = \frac{1}{\sqrt{(l+m)(l-m+1)}} $$

$Y_{lm}(\theta)$ に下降演算子 $L_-$ を作用させると $$ \begin{align*} L_- Y_{lm}(\theta,\phi) &= \sqrt{(l+m)(l-m+1)}Y_{l,m-1}(\theta,\phi) \\ &= \sqrt{\frac{(l+m)(l-m+1)}{2\pi}}\em^{i(m-1)\phi}\Theta_{l,m-1}(\theta) \end{align*} $$ 一方で $L_-$ の座標表示\eqref{eq:Lms}を考えると $$ \begin{align*} L_- Y_{lm}(\theta,\phi) & = \em^{-i\phi}\left( -\pdv{\theta}+i\cot\theta\pdv{\phi} \right) Y_{lm}(\theta,\phi) \\ &= \frac{1}{\sqrt{2\pi}}\em^{-i\phi}\em^{im\phi}\left( -\dv{\theta}-m\cot\theta \right) \Theta_{lm}(\theta) \\ &= \frac{1}{\sqrt{2\pi}}\em^{i(m-1)\phi}\left( \sin\theta\dv{\Theta_{lm}(\theta)}{(\cos\theta)}+m\dv{\sin\theta}{(\cos\theta)}\Theta_{lm}(\theta) \right) \\ &= \frac{1}{\sqrt{2\pi}}\em^{i(m-1)\phi} \frac{\sin\theta}{\sin^m\theta}\dv{(\cos\theta)}\Bigl( \sin^m\theta\cdot\Theta_{lm}(\theta) \Bigr) \end{align*} $$ 2行目から3行目の変形では、 $$ \dv{\theta} = \dv{\cos\theta}{\theta}\dv{(\cos\theta)} = -\sin\theta\dv{(\cos\theta)} $$ および、$\cos^2\theta+\sin^2\theta=1$ の両辺を $\cos\theta$ で微分して $$ 2\cos\theta+2\sin\theta\dv{\sin\theta}{(\cos\theta)}=0 \;\;\Longrightarrow\;\; \cot\theta = -\dv{\sin\theta}{(\cos\theta)} $$ を用いた。
右辺どうしを等号で結んで整理すると、与式を得る。

この漸化式を $(l-m)$ 回用いて $Y_{lm}(\theta,\phi)$ を変形していくと、式\eqref{eq:YTheta}右辺の $\Theta_{lm}(\theta)$ は、具体形が式\eqref{eq:Theta_ll}で分かっている $\Theta_{ll}(\theta)$ で表せるようになります。
$\Theta_{ll}(\theta)$ を代入することで、以下が答えとなります。

問7

$$ \begin{equation} Y_{lm}(\theta,\phi) = \frac{1}{\sqrt{2\pi}}\em^{im\phi}\cdot \sqrt{\frac{2l+1}{2}\frac{(l-m)!}{(l+m)!}}\frac{1}{2^ll!}\frac{1}{\sin^m\theta}\dv[l-m]{(\cos\theta)}\sin^{2l}\theta \label{eq:Y_annihilation} \end{equation} $$

式\eqref{eq:YTheta}に対して変形する。 $$ \begin{align*} Y_{lm}(\theta,\phi) &= \frac{1}{\sqrt{2\pi}}\em^{im\phi}\Theta_{lm}(\theta) \\ &= \frac{1}{\sqrt{2\pi}}\em^{im\phi}\frac{1}{\sin^m\theta}\cdot \Bigl( \sin^m\theta\cdot\Theta_{lm}(\theta) \Bigr) \\ &= \gamma_{l,m+1} \frac{1}{\sqrt{2\pi}}\em^{im\phi}\frac{1}{\sin^m\theta}\cdot \dv{(\cos\theta)}\Bigl( \sin^{m+1}\theta\cdot\Theta_{l,m+1}(\theta) \Bigr) \\ &= \gamma_{l,m+1}\gamma_{l,m+2} \frac{1}{\sqrt{2\pi}}\em^{im\phi}\frac{1}{\sin^m\theta}\cdot \dv[2]{(\cos\theta)}\Bigl( \sin^{m+2}\theta\cdot\Theta_{l,m+2}(\theta) \Bigr) \\ &= \;\cdots \\ &= (\gamma_{l,m+1}\gamma_{l,m+2}\cdots\gamma_{ll}) \frac{1}{\sqrt{2\pi}}\em^{im\phi}\frac{1}{\sin^m\theta}\cdot \dv[l-m]{(\cos\theta)}\Bigl( \sin^{l}\theta\cdot\Theta_{ll}(\theta) \Bigr) \end{align*} $$ $\Theta_{ll}(\theta)$ に式\eqref{eq:Theta_ll}を代入すると $$ \begin{align*} Y_{lm}(\theta,\phi) &= (\gamma_{l,m+1}\gamma_{l,m+2}\cdots\gamma_{ll}) \frac{1}{\sqrt{2\pi}}\em^{im\phi}\frac{1}{\sin^m\theta}\cdot \dv[l-m]{(\cos\theta)}\left(\frac{1}{2^ll!}\sqrt{\frac{(2l+1)!}{2}}\sin^{2l}\theta \right) \\ &= \frac{1}{\sqrt{2\pi}}\em^{im\phi} \cdot (\gamma_{l,m+1}\gamma_{l,m+2}\cdots\gamma_{ll})\sqrt{\frac{(2l+1)!}{2}} \cdot \frac{1}{2^ll!}\frac{1}{\sin^m\theta}\dv[l-m]{(\cos\theta)}\sin^{2l}\theta \end{align*} $$ 係数部分を整理する。 $$ \begin{align*} & (\gamma_{l,m+1}\gamma_{l,m+2}\cdots\gamma_{ll})\sqrt{\frac{(2l+1)!}{2}} \\ &= \frac{1}{\sqrt{(l+m+1)(l-m)}}\cdot\frac{1}{\sqrt{(l+m+2)(l-m-1)}}\cdots\frac{1}{\sqrt{2l\cdot1}}\cdot\sqrt{\frac{(2l+1)!}{2}} \\ &= \sqrt{\frac{(l+m)!}{(2l)!(l-m)!}}\cdot\sqrt{\frac{(2l+1)!}{2}} \\ &= \sqrt{\frac{2l+1}{2}\frac{(l+m)!}{(l-m)!}} \end{align*} $$

$\ket{l,m}$ の時点で規格化まで考慮されているため、$\Theta_{ll}(\theta)$ の係数を式\eqref{eq:Theta_ll}のように合わせておけば、正規直交性 $$ \int\mathrm{d}\Omega\,Y_{lm}^\ast(\theta,\phi)Y_{l’m’}(\theta,\phi) = \delta_{ll’}\delta_{mm’} $$ が自動的に成り立ちます。直交性は、演算子 $\bm{L}^2, L_z$ に対して異なる固有値を持つことから自明です。

導出過程は省略しますが、$\Theta_{l,-l}(\theta)$ に上昇演算子を作用させる方法でも、本質的には同じ答えが得られます。

表記の違い

球面調和関数は、導出方法や流派によってさまざまな表記が存在します。
典型的な6種類の表記を、以下にまとめます。

$$ \begin{equation} Y_{lm}(\theta,\phi) = \frac{1}{\sqrt{2\pi}}\em^{im\phi}\sqrt{\frac{2l+1}{2}\frac{(l-m)!}{(l+m)!}}P_l^m(\cos\theta) \label{eq:Y_A} \tag{A} \end{equation} $$ $$ \begin{equation} Y_{lm}(\theta,\phi) = \frac{1}{\sqrt{2\pi}}\em^{im\phi}\sqrt{\frac{2l+1}{2}\frac{(l+m)!}{(l-m)!}}\frac{1}{2^ll!}\frac{1}{\sin^m\theta}\dv[l-m]{(\cos\theta)}\sin^{2l}\theta \label{eq:Y_B} \tag{B} \end{equation} $$ $$ \begin{equation} Y_{lm}(\theta,\phi) = (-1)^{l+m}\frac{1}{\sqrt{2\pi}}\em^{im\phi}\sqrt{\frac{2l+1}{2}\frac{(l-m)!}{(l+m)!}}\frac{1}{2^ll!}\sin^m\theta\dv[l+m]{(\cos\theta)}\sin^{2l}\theta \label{eq:Y_C} \tag{C} \end{equation} $$ $$ \begin{equation} Y_{lm}(\theta,\phi) = (-1)^\frac{|m|-m}{2}\frac{1}{\sqrt{2\pi}}\em^{im\phi}\sqrt{\frac{2l+1}{2}\frac{(l-|m|)!}{(l+|m|)!}}P_l^{|m|}(\cos\theta) \label{eq:Y_D} \tag{D} \end{equation} $$ $$ \begin{equation} Y_{lm}(\theta,\phi) = (-1)^\frac{|m|-m}{2}\frac{1}{\sqrt{2\pi}}\em^{im\phi}\sqrt{\frac{2l+1}{2}\frac{(l+|m|)!}{(l-|m|)!}}\frac{1}{2^ll!}\frac{1}{\sin^{|m|}\theta}\dv[l-|m|]{(\cos\theta)}\sin^{2l}\theta \label{eq:Y_E} \tag{E} \end{equation} $$ $$ \begin{equation} Y_{lm}(\theta,\phi) = (-1)^\frac{|m|+m}{2}\frac{1}{\sqrt{2\pi}}\em^{im\phi}\sqrt{\frac{2l+1}{2}\frac{(l-|m|)!}{(l+|m|)!}}\frac{1}{2^ll!}\sin^{|m|}\theta\dv[l+|m|]{(\cos\theta)}\sin^{2l}\theta \label{eq:Y_F} \tag{F} \end{equation} $$

これらには符号の違いがあります。

(A) = (C) = (D)
$\updownarrow \;\times(-1)^l$
(B) = (E) = (F)

$\theta$ 部分の固有関数をどのように求めるかによって、得られる表式が変わります。
それぞれの表式は、以下の方法で自然に導出されます。

  • Legendre陪関数を利用する方法 → (A)(式\eqref{eq:Y_main}と同じ)
  • $\Theta_{ll}(\theta)$ に下降演算子を作用させる方法 → (B)(式\eqref{eq:Y_annihilation}と同じ)
  • $\Theta_{l,-l}(\theta)$ に上昇演算子を作用させる方法 → (C)
  • $m$ で場合分けして (A) を組み合わせたもの → (D)
  • $m$ で場合分けして (B), (C) を組み合わせたもの → (E)
  • $m$ で場合分けして (B), (C) を組み合わせたもの → (F)

符号の違いは微分方程式からは決まらない任意の位相因子であり、どの定義を採用するかという慣習の問題です。
物理量の計算では、内積 $\int\mathrm{d}\Omega\,Y_{lm}^* Y_{l’m’}$ の形で現れるので、観測可能な物理現象に影響は出ません。しかし学習時や数値計算の際には、文献がどの定義を用いているか意識する必要があります。

ハーポルホード高橋
ハーポルホード高橋

導出を一度理解した後は、(A) の定義が一番シンプルで実用的です。
Legendre陪関数の性質を、数学公式として使うことができて便利です。

表式 (A) ~ (F) の等価性と符号の違いの証明は、以下の通りです。

問8-1

(A) $\Rightarrow$ (C)

(A) の $P_l^m(\cos\theta)$ に、Rodriguesの公式\eqref{eq:rodrigues}を適用する。 $$ \begin{align*} P_l^m(\cos\theta) &= \frac{(-1)^m}{2^ll!}(1-\cos^2\theta)^\frac{m}{2}\dv[l+m]{(\cos\theta)}(\cos^2\theta-1)^l \\ &= (-1)^{l+m}\frac{1}{2^ll!}\sin^m\theta\dv[l+m]{(\cos\theta)}\sin^{2l}\theta \end{align*} $$ よって (A) は (C) に一致する。

問8-2

(A) $\Rightarrow$ $(-1)^l\times$ (B)

公式\eqref{eq:rodrigues}, \eqref{eq:Legendre_-m}より $$ \begin{align*} P_l^m(x) &= (-1)^m\frac{(l+m)!}{(l-m)!}P_l^{-m}(x) \\ &= (-1)^m\frac{(l+m)!}{(l-m)!} \cdot \frac{(-1)^{-m}}{2^ll!}(1-x^2)^{-\frac{m}{2}}\dv[l-m]{x}(x^2-1)^l \end{align*} $$ したがって $$ \begin{align*} \sqrt{\frac{(l-m)!}{(l+m)!}}P_l^m(\cos\theta) &= \sqrt{\frac{(l-m)!}{(l+m)!}} \cdot \frac{(l+m)!}{(l-m)!}\frac{1}{2^ll!}(1-\cos^2\theta)^{-\frac{m}{2}}\dv[l-m]{(\cos\theta)}(\cos^2\theta-1)^l \\ &= (-1)^l\sqrt{\frac{(l+m)!}{(l-m)!}}\frac{1}{2^ll!}\frac{1}{\sin^m\theta}\dv[l-m]{(\cos\theta)}\sin^{2l}\theta \end{align*} $$ 右辺先頭の $(-1)^l$ が余分に掛かるため、(A) は $(-1)^l\times$ (B) に一致する。

問8-3

(D) $\Rightarrow$ (A)

$m\geq0$ のとき、単に絶対値を外せば (D) は (A) に一致する。
$m<0$ のとき、公式\eqref{eq:Legendre_-m}より $$ \begin{align*} (-1)^\frac{|m|-m}{2}\sqrt{\frac{(l-|m|)!}{(l+|m|)!}}P_l^{|m|}(\cos\theta) &= (-1)^{-m}\sqrt{\frac{(l+m)!}{(l-m)!}}P_l^{-m}(\cos\theta) \\ &= (-1)^m\sqrt{\frac{(l+m)!}{(l-m)!}} \cdot (-1)^m\frac{(l-m)!}{(l+m)!}P_l^m(\cos\theta) \\ &= \sqrt{\frac{(l-m)!}{(l+m)!}}P_l^m(\cos\theta) \end{align*} $$ となるので、(D) は (A) に一致する。

問8-4

$m\geq0$ のとき (E) $\Rightarrow$ (B)
$m<0$ のとき (E) $\Rightarrow$ $(-1)^l\times$ (C)

$m\geq0$ のとき、単に絶対値を外せば (E) は (B) に一致する。
$m<0$ のとき、 $$ \begin{align*} & (-1)^\frac{|m|-m}{2}\sqrt{\frac{(l+|m|)!}{(l-|m|)!}}\frac{1}{\sin^{|m|}\theta}\dv[l-|m|]{(\cos\theta)}\sin^{2l}\theta \\ =&\; (-1)^{-m}\sqrt{\frac{(l-m)!}{(l+m)!}}\frac{1}{\sin^{-m}\theta}\dv[l+m]{(\cos\theta)}\sin^{2l}\theta \\ =&\; (-1)^l\cdot (-1)^{l+m}\sqrt{\frac{(l-m)!}{(l+m)!}}\sin^m\theta\dv[l+m]{(\cos\theta)}\sin^{2l}\theta \end{align*} $$ となるので、(E) は $(-1)^l\times$ (C) に一致する。

問8-5

$m\geq0$ のとき (F) $\Rightarrow$ $(-1)^l\times$ (C)
$m<0$ のとき (F) $\Rightarrow$ (B)

$m\geq0$ のとき、単に絶対値を外して $$ (-1)^\frac{|m|+m}{2} = (-1)^m = (-1)^l \cdot (-1)^{l+m} $$ に注意すれば、 (F) は $(-1)^l\times$ (C) に一致する。
$m<0$ のとき、 $$ \begin{align*} & (-1)^\frac{|m|+m}{2}\sqrt{\frac{(l-|m|)!}{(l+|m|)!}}\sin^{|m|}\theta\dv[l+|m|]{(\cos\theta)}\sin^{2l}\theta \\ =&\; \sqrt{\frac{(l+m)!}{(l-m)!}}\sin^{-m}\theta\dv[l-m]{(\cos\theta)}\sin^{2l}\theta \\ =&\; \sqrt{\frac{(l+m)!}{(l-m)!}}\frac{1}{\sin^m\theta}\dv[l-m]{(\cos\theta)}\sin^{2l}\theta \end{align*} $$ となるので、(F) は (B) に一致する。

具体例

$l=0,1,2,3$ について、Legendre陪関数 $P_l^m(x)$ と球面調和関数 $Y_{lm}(\theta,\phi)$ の一覧表を示します。
Legendre陪関数は式\eqref{eq:rodrigues}、球面調和関数は式\eqref{eq:Y_main}((A)と同じ)の定義に従います。

$l$$m$$P_l^m(x)$$Y_{lm}(\theta,\phi)$
$0$$0$$1$$\frac{1}{4\pi}$
$1$$-1$$\frac{1}{2}(1-x^2)^{1/2}$$\sqrt{\frac{3}{8\pi}}\sin\theta\,\em^{-i\phi}$
$1$$0$$x$$\sqrt{\frac{3}{4\pi}}\cos\theta$
$1$$1$$-(1-x^2)^{1/2}$$-\sqrt{\frac{3}{8\pi}}\sin\theta\,\em^{i\phi}$
$2$$-2$$\frac{1}{8}(1-x^2)$$\sqrt{\frac{15}{32\pi}}\sin^2\theta\,\em^{-2i\phi}$
$2$$-1$$\frac{1}{2}x(1-x^2)^{1/2}$$\sqrt{\frac{15}{8\pi}}\sin\theta\cos\theta\,\em^{-i\phi}$
$2$$0$$\frac{1}{2}(3x^2-1)$$\sqrt{\frac{5}{16\pi}}(3\cos^2\theta-1)$
$2$$1$$-3x(1-x^2)^{1/2}$$-\sqrt{\frac{15}{8\pi}}\sin\theta\cos\theta\,\em^{i\phi}$
$2$$2$$3(1-x^2)$$\sqrt{\frac{15}{32\pi}}\sin^2\theta\,\em^{2i\phi}$
$3$$-3$$\frac{1}{48}(1-x^2)^{3/2}$$\sqrt{\frac{35}{64\pi}}\sin^3\theta\,\em^{-3i\phi}$
$3$$-2$$\frac{1}{8}x(1-x^2)$$\sqrt{\frac{105}{32\pi}}\sin^2\theta\cos\theta\,\em^{-2i\phi}$
$3$$-1$$\frac{1}{8}(5x^2-1)(1-x^2)^{1/2}$$\sqrt{\frac{21}{64\pi}}\sin\theta(5\cos^2\theta-1)\,\em^{-i\phi}$
$3$$0$$\frac{1}{2}(5x^3-3x)$$\sqrt{\frac{7}{16\pi}}(5\cos^3\theta-3\cos\theta)$
$3$$1$$-\frac{3}{2}(5x^2-1)(1-x^2)^{1/2}$$-\sqrt{\frac{21}{64\pi}}\sin\theta(5\cos^2\theta-1)\,\em^{i\phi}$
$3$$2$$15x(1-x^2)$$\sqrt{\frac{105}{32\pi}}\sin^2\theta\cos\theta\,\em^{2i\phi}$
$3$$3$$-15(1-x^2)^{3/2}$$-\sqrt{\frac{35}{64\pi}}\sin^3\theta\,\em^{3i\phi}$

参考文献

タイトルとURLをコピーしました