logo

球面座標系におけるシュレーディンガー方程式 📂量子力学

球面座標系におけるシュレーディンガー方程式

方程式について

球座標系ではシュレディンガー方程式は次の通りだ。

22M[1r2r(r2ψr)+1r2sinθθ(sinθψθ)+1r2sin2θ2ψ2ϕ]+Vψ=Eψ(1) -\frac{\hbar^{2}}{2M}\left[\frac{1}{r^2}\frac{\partial}{\partial r} \left( r^2\frac{\partial \psi}{\partial r} \right) + \frac{1}{r^2\sin\theta}\frac{\partial}{\partial\theta}\left( \sin\theta \frac{\partial \psi}{\partial \theta} \right) + \frac{1}{r^2\sin^2\theta}\frac{\partial^2 \psi}{\partial^2 \phi} \right]+V\psi=E\psi \tag{1}

説明

3次元の時間に依存しないシュレディンガー方程式は次の通りだ。

22M2ψ+Vψ=Eψ -\frac{\hbar^{2}}{2M}\nabla^{2}\psi+V\psi=E\psi

ここで、MMは粒子の質量だ。ポテンシャルV=V(r)V = V(r)のように原点からの距離のみに依存するならば、問題を球座標系で解くのが良い。

球座標系でのラプラシアン

2f=1r2r(r2fr)+1r2sinθθ(sinθfθ)+1r2sin2θ2fϕ2 \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 \phi^{2} }

ラプラシアンが上記のようになっているので、シュレディンガー方程式は次の通りだ。

22M[1r2r(r2ψr)+1r2sinθθ(sinθψθ)+1r2sin2θ2ψ2ϕ]+Vψ=Eψ(1) -\frac{\hbar^{2}}{2M}\left[\frac{1}{r^2}\frac{\partial}{\partial r} \left( r^2\frac{\partial \psi}{\partial r} \right) + \frac{1}{r^2\sin\theta}\frac{\partial}{\partial\theta}\left( \sin\theta \frac{\partial \psi}{\partial \theta} \right) + \frac{1}{r^2\sin^2\theta}\frac{\partial^2 \psi}{\partial^2 \phi} \right] + V\psi = E\psi \tag{1}

波動関数ψ\psiが以下のように変数分離可能だと仮定しよう。

ψ(r,θ,ϕ)=R(r)Θ(θ)Φ(ϕ) \psi (r,\theta,\phi) = R(r) \Theta(\theta) \Phi( \phi)

すると、式(1)(1)は次の通りだ。

22M[ΘΦr2ddr(r2dRdr)+RΦr2sinθddθ(sinθdΘdθ)+RΘr2sin2θd2Φd2ϕ]+VRΘΦ=ERΘΦ -\frac{\hbar^{2}}{2M}\left[\frac{\Theta \Phi}{r^2}\frac{d}{d r} \left( r^2\frac{d R}{d r} \right) + \frac{R \Phi}{r^2\sin\theta}\frac{d}{d\theta}\left( \sin\theta \frac{d \Theta}{d \theta} \right) + \frac{R \Theta}{r^2\sin^2\theta}\frac{d^2 \Phi}{d^2 \phi} \right]+VR\Theta \Phi=ER\Theta \Phi

右辺の項を左辺に移項し、両辺に2Mr221RΘΦ-\dfrac{2Mr^{2}}{\hbar ^{2}}\dfrac{1}{R\Theta \Phi}を掛けて整理すると、以下のようになる。

[1Rddr(r2dRdr)2Mr22(V(r)E)]+[1Θsinθddθ(sinθdΘdθ)+1Φsin2θd2Φd2ϕ]=0 \left[ \frac{1}{R}\frac{d}{d r} \left( r^2\frac{d R}{d r} \right) -\frac{2Mr^{2}}{\hbar ^{2}}(V(r)-E) \right]+\left[ \frac{1}{\Theta \sin\theta}\frac{d}{d\theta}\left( \sin\theta \frac{d \Theta}{d \theta} \right) + \frac{1}{\Phi\sin^2\theta}\frac{d^2 \Phi}{d^2 \phi} \right]=0

ポテンシャル項と定数が追加されたが、基本的には球面座標系のラプラス方程式を解くことと大きな枠は同じだ。最初の中括弧で囲まれた項は変数rrのみに関する項で、第二項は変数θ\thetaおよびϕ\phiのみに影響を受ける項なので、それぞれの中括弧部分は定数だ。最初の項を(+1)\ell(\ell + 1)とする。すると第二項は(+1)-\ell(\ell + 1)だ。

1Rddr(r2dRdr)2Mr22(V(r)E)=(+1)1Θsinθddθ(sinθdΘdθ)+1Φsin2θd2Φd2ϕ=(+1) \begin{align*} \frac{1}{R}\frac{d}{d r} \left( r^2\frac{d R}{d r} \right) -\frac{2Mr^{2}}{\hbar ^{2}}(V(r)-E) &= \ell(\ell+1) \tag{2} \\ \frac{1}{\Theta \sin\theta}\frac{d}{d\theta}\left( \sin\theta \frac{d \Theta}{d \theta} \right) + \frac{1}{\Phi\sin^2\theta}\frac{d^2 \Phi}{d^2 \phi} &= -\ell(\ell+1) \tag{3} \end{align*}

角度に関する方程式(3)(3)の解は特に球面調和関数と呼ばれ、次の通りだ。 Ym(θ,ϕ)=eimϕPm(cosθ) Y_{\ell}^{m}(\theta, \phi) = e^{\i m\phi}P_{\ell}^{m}(\cos \theta)

(3)(3)からmmが直接的に見えないが、変数分離の過程でθ\thetaϕ\phiを分離定数として現れる。Pm(cosθ)P_{\ell}^{m}(\cos\theta)連立レジャンドル多項式だ。規格化すると以下のようになる。

Ym(θ,ϕ)=2+14π(m)!(+m)!Pm(cosθ)eimϕ Y_{\ell}^{m} (\theta,\phi)=\sqrt{\frac{2\ell+1}{4\pi}\frac{(\ell-m)!}{(\ell+m)!}}P_{\ell}^{m}(\cos\theta)e^{\i m\phi}

次に、残っているのは半径成分の方程式(2)(2)だ。両辺にRRを掛けて、右辺にエネルギーEEに関する項だけが残るように整理する。

22Mr2ddr(r2dRdr)+(V22Ml(l+1)r2)R=ER(4) -\frac{ \hbar^{2} }{ 2Mr^{2} }\frac{d}{d r} \left( r^2\frac{d R}{d r} \right) +\left(V- \frac{\hbar ^{2}}{2M}\frac{l(l+1)}{r^{2}} \right)R =ER \tag{4}

この時にrR(r)=u(r)rR(r)=u(r)で置換すると、方程式が簡単になる。

R=urdRdr=1rdudr1r2uddr(r2dRdr)=ddr(rdudru)=rd2ud2 \begin{align*} R&=\frac{u}{r} \\ \frac{ dR }{ dr }&=\frac{1}{r}\frac{ d u }{ d r}-\frac{1}{r^{2}}u \\ \frac{ d }{ d r }\left(r^{2} \frac{ d R}{ d r } \right)&=\frac{ d }{ dr }\left( r \frac{ d u }{ dr }-u \right)=r\frac{ d ^{2}u}{ d^{2} } \end{align*}

よって(4)(4)は以下の通りだ。

22Mr2rd2udr2+(V22M(+1)r2)ur=Eur    22Md2udr2+(V22M(+1)r2)u=Eu \begin{align*} &&-\frac{ \hbar^{2} }{ 2Mr^{2} }r\frac{ d ^{2}u}{ dr^{2} } +\left(V- \frac{\hbar ^{2}}{2M}\frac{\ell(\ell+1)}{r^{2}} \right)\frac{u}{r} &= E\frac{u}{r} \\ \implies &&-\frac{ \hbar^{2} }{ 2M }\frac{ d ^{2}u}{ dr^{2} } +\left(V- \frac{\hbar ^{2}}{2M}\frac{\ell(\ell+1)}{r^{2}} \right)u &= Eu \end{align*}

この形は以下の1次元シュレディンガー方程式と非常に似ている。

22md2ψdx2+Vψ=Eψ -\frac{\hbar^{2}}{2m}\frac{ d ^{2} \psi}{ d x^{2} }+V\psi = E\psi

違いはポテンシャルがV=V22M(+1)r2V = V- \dfrac{\hbar ^{2}}{2M}\dfrac{\ell(\ell+1)}{r^{2}}に変わったことだけだ。1r2\dfrac{1}{r^{2}}に比例する第二項を遠心力項と呼ぶ。u(r)u(r)Ylm(cosθ)Y_{l}^{m}(\cos\theta)と同じく規格化条件を満たさなければならない。

0R(r)2r2dr=0u(r)2dr \int_{0}^{\infty}|R(r)|^{2}r^{2}dr=\int _{0}^{\infty} |u(r)|^{2}dr

V(r)V(r)の一般的な解はここまでだ。あとは与えられた問題でポテンシャル関数V(r)V(r)が正確に与えられていれば、それに従って解決できる。