[量子力学] 中心場での粒子の運動

April 24, 2020      

重心運動と相対運動

古典力学において2つの粒子が相互作用する運動は,重心運動と相対運動の形に書き換えることができた.まずはこの変形について確認しておこう.2つの粒子が相互作用するような系のハミルトニアンは次のように表される.

H=m1r˙122+m1r˙222+U(r2r1)\begin{equation} H = \frac{m_1 \dot{\boldsymbol{r}}_1^2}{2} + \frac{m_1 \dot{\boldsymbol{r}}_2^2}{2} + U(|\boldsymbol{r}_2 - \boldsymbol{r}_1|) \end{equation}

ここで,2つの粒子の重心ベクトルR\boldsymbol{R}と相対位置のベクトルr\boldsymbol{r}を導入する.

R=m1r1+m2r2m1+m2r=r2r1\begin{align} \boldsymbol{R} &= \frac{m_1 \boldsymbol{r}_1 + m_2 \boldsymbol{r}_2}{m_1 + m_2} \\ \boldsymbol{r} &= \boldsymbol{r}_2 - \boldsymbol{r}_1 \end{align}

r1,r2\boldsymbol{r}_1, \boldsymbol{r}_2について解いてやれば以下の通り.

r1=Rm2m1+m2rr2=R+m1m1+m2r\begin{align} \boldsymbol{r}_1 &= \boldsymbol{R} - \frac{m_2}{m_1 + m_2} \boldsymbol{r} \\ \boldsymbol{r}_2 &= \boldsymbol{R} + \frac{m_1}{m_1 + m_2} \boldsymbol{r} \end{align}

これを用いてハミルトニアンを書き直すと,重心位置と相対位置に関する項に分離することができる.ここでm=m1m2m1+m2m = \frac{m_1 m_2}{m_1 + m_2}とおいて,この量を換算質量(reduced mass)と呼ぶ.

H=m12(R˙22m2(R˙r˙)m1+m2 +m22r˙2(m1+m2)2 )+m22(R˙2+2m1(R˙r˙)m1+m2+m12r˙2(m1+m2)2)+U(r)=m1+m22R˙2+12m1m2m1+m2r˙2+U(r)\begin{align*} H &= \frac{m_1}{2} \left(\dot{\boldsymbol{R}}^2 -\frac{2m_2 (\dot{\boldsymbol{R}} \cdot \dot{\boldsymbol{r}})}{m_1+m_2}  + \frac{m_2^2 \dot{\boldsymbol{r}}^2}{(m_1+m_2)^2} \right) + \frac{m_2}{2} \left(\dot{\boldsymbol{R}}^2 +\frac{2m_1 (\dot{\boldsymbol{R}} \cdot \dot{ \boldsymbol{r}})}{m_1+m_2} + \frac{m_1^2 \dot{\boldsymbol{r}}^2}{(m_1+m_2)^2} \right) + U(|\boldsymbol{r}|) \\ &= \frac{m_1 + m_2}{2} \dot{\boldsymbol{R}}^2 + \frac{1}{2} \frac{m_1 m_2}{m_1 + m_2} \dot{\boldsymbol{r}}^2 + U(|\boldsymbol{r}|) \end{align*}

同様に量子的な2粒子系についても,ハミルトニアン(演算子)は(6)から(7)の形に変形することができる.

H^=22m1Δ122m2Δ2+U(r)\begin{equation} % \label{eq:LandauLifshitz_Quantum_32.1} \hat{H} = - \frac{\hbar^2}{2m_1} \Delta_1 - \frac{\hbar^2}{2m_2} \Delta_2 + U(r) \end{equation}
H^=22(m1+m2)ΔR22mΔr+U(r)\begin{equation} % \label{eq:LandauLifshitz_Quantum_32.3} \hat{H} = - \frac{\hbar^2}{2(m_1+m_2)} \Delta_R - \frac{\hbar^2}{2m} \Delta_r + U(r) \end{equation}

このとき波動関数は,ψ(r)φ(R)\psi(r)\varphi(R)の形で表すことができる.念のため確認しておこう.

22(m1+m2)ΔR(φψ)22mΔr(φψ)+U(r)(φψ)=E(φψ)2ψ2(m1+m2)ΔRφ2φ2mΔrψ+U(r)(φψ)=(ER+Er)(φψ){2ψ2(m1+m2)ΔRφ=ER(φψ)2φ2mΔrψ+U(r)(φψ)=Er(φψ)\begin{gather} - \frac{\hbar^2}{2(m_1+m_2)} \Delta_R (\varphi \psi) - \frac{\hbar^2}{2m} \Delta_r (\varphi \psi) + U(r) (\varphi \psi) = E (\varphi \psi) \\ - \frac{\hbar^2 \psi}{2(m_1+m_2)} \Delta_R \varphi - \frac{\hbar^2 \varphi}{2m} \Delta_r \psi + U(r) (\varphi \psi) = (E_R + E_r) (\varphi \psi) \\ \to \left\{ \begin{array}{l} - \frac{\hbar^2 \psi}{2(m_1 + m_2)} \Delta_R \varphi = E_R (\varphi \psi) \\ - \frac{\hbar^2 \varphi}{2m} \Delta_r \psi + U(r) (\varphi \psi) = E_r (\varphi \psi) \end{array} \right. \end{gather}

ラプラシアンに関係ないφ,ψ\varphi, \psiを掃ってやれば,重心位置と相対位置ばらばらにシュレーディンガー方程式を解いてやればよいことがわかる.相対位置に関しては,1粒子が原点からの距離に応じたポテンシャルの中を運動する式に対応する.

{22(m1+m2)ΔRφ=ERφ22mΔrψ+U(r)ψ=Erψ\begin{equation} \left\{ \begin{array}{l} - \frac{\hbar^2}{2(m_1 + m_2)} \Delta_R \varphi = E_R \varphi \\ - \frac{\hbar^2}{2m} \Delta_r \psi + U(r) \psi= E_r \psi \end{array} \right. \end{equation}

中心場での粒子の運動

中心場での粒子の運動は以下の式で表される.

Δψ+2m2[EU(r)]ψ=0\begin{equation} % \label{eq:LandauLifshitz_Quantum_32.4} \Delta \psi + \frac{2m}{\hbar^2} [E - U(r)] \psi = 0 \end{equation}

原点に関して球対称な式なので,ここからは極座標で考えよう.極座標系のラプラシアンは次のように表された([流体力学] 円筒座標・極座標のナブラとラプラシアン).

Δ=2r2+2rr+1r22θ2+1r2tanθθ+1r2sin2θ2ϕ2=1r2r(r2r)+1r2[1sinθθ(sinθθ)+1sin2θ22ϕ]\begin{align} \Delta &= \frac{\partial^2}{\partial r^2} + \frac{2}{r} \frac{\partial}{\partial r} + \frac{1}{r^2} \frac{\partial^2}{\partial \theta^2} + \frac{1}{r^2 \tan \theta} \frac{\partial}{\partial \theta} + \frac{1}{r^2 \sin^2 \theta} \frac{\partial^2}{\partial \phi^2} \\ &= \frac{1}{r^2} \frac{\partial}{\partial r} \left( r^2 \frac{\partial}{\partial r} \right) + \frac{1}{r^2} \left[ \frac{1}{\sin \theta} \frac{\partial}{\partial \theta} \left( \sin \theta \frac{\partial}{\partial \theta} \right) + \frac{1}{\sin^2 \theta} \frac{\partial^2}{\partial^2 \phi} \right] \end{align}

また,角運動量の2乗の演算子は次のように表された([量子力学] 角運動量).

l^2=1sin2θ22ϕcosθsinθθ22θ\begin{equation} \hat{l}_2 = - \frac{1}{\sin^2 \theta} \frac{\partial^2}{\partial^2 \phi}- \frac{\cos \theta}{\sin \theta} \frac{\partial}{\partial \theta} - \frac{\partial^2}{\partial^2 \theta} \end{equation}

これらを用いると(12)は次のように書き直せる.

22m[1r2r(r2ψr)+l^2r2ψ]+U(r)ψ=Eψ\begin{equation} % \label{eq:LandauLifshitz_Quantum_32.6} \frac{\hbar^2}{2m} \left[ - \frac{1}{r^2} \frac{\partial}{\partial r} \left( r^2 \frac{\partial \psi}{\partial r} \right) + \frac{\hat{l}_2}{r^2} \psi \right] + U(r) \psi = E \psi \end{equation}

この式の解を(17)の形と仮定しよう.ただし,Ylm(θ,ϕ)Y_{lm}(\theta, \phi)は角運動量の2乗の固有関数になっていた([量子力学] 角運動量).

ψ=R(r)Ylm(θ,ϕ)\begin{equation} % \label{eq:LandauLifshitz_Quantum_32.7} \psi = R(r) Y_{lm} (\theta, \phi) \end{equation}
l^2Ylm=l(l+1)Ylm\begin{equation} \hat{l}_2 Y_{lm} = l(l+1) Y_{lm} \end{equation}

これを実際に(16)に代入すると,以下のように整理できる.

22m[Ylm(θ,ϕ)r2r(r2R(r)r)+R(r)l(l+1)r2Ylm(θ,ϕ)]+U(r)R(r)Ylm(θ,ϕ)=ER(r)Ylm(θ,ϕ)1r2ddr(r2dR(r)dr)R(r)l(l+1)r2+2m2[EU(r)]R(r)=0\begin{gather*} % \label{eq:LandauLifshitz_Quantum_32.8} \frac{\hbar^2}{2m} \left[ - \frac{Y_{lm}(\theta, \phi)}{r^2} \frac{\partial}{\partial r} \left( r^2 \frac{\partial R(r)}{\partial r} \right) + R(r) \frac{l(l+1)}{r^2} Y_{lm}(\theta, \phi) \right] + U(r) R(r)Y_{lm}(\theta, \phi) = E R(r)Y_{lm}(\theta, \phi) \\ \frac{1}{r^2} \frac{d}{d r} \left( r^2 \frac{d R(r)}{d r} \right) - R(r) \frac{l(l+1)}{r^2} + \frac{2m}{\hbar^2} [E - U(r)] R(r) = 0 \end{gather*}

ここで,次のような置換を行う.

R(r)=χ(r)rdR(r)dr=1rdχdrχr2ddr(rdχdrχ)=dχdr+rd2rdr2dχdr=rd2rdr2\begin{gather} % \label{eq:LandauLifshitz_Quantum_32.9} R(r) = \frac{\chi(r)}{r} \\ \frac{d R(r)}{d r} = \frac{1}{r} \frac{d\chi}{dr} - \frac{\chi}{r^2} \\ \frac{d}{d r}\left( r \frac{d\chi}{dr} - \chi \right) = \frac{d\chi}{dr} + r \frac{d^2 r}{dr^2} - \frac{d\chi}{dr} = r \frac{d^2 r}{dr^2} \end{gather}
1rd2χdr2χ(r)rl(l+1)r2+2m2[EU(r)]χ(r)r=0d2χdr2+[2m2(EU(r))l(l+1)r2]χ=0\begin{gather} % \label{eq:LandauLifshitz_Quantum_32.10} \frac{1}{r} \frac{d^2 \chi}{dr^2} - \frac{\chi(r)}{r} \frac{l(l+1)}{r^2} + \frac{2m}{\hbar^2} [E-U(r)] \frac{\chi(r)}{r} = 0 \\ \frac{d^2 \chi}{dr^2} + \left[ \frac{2m}{\hbar^2} (E-U(r)) - \frac{l(l+1)}{r^2} \right] \chi = 0 \end{gather}

ポテンシャルを次のように置きなおせば,(23)はポテンシャルUl(r)U_l(r)の中の1次元運動を表すシュレーディンガー方程式と考えられる.

d2χdr2+2m2[EUl(r)]χ=0,  where Ul(r)=U(r)+22ml(l+1)r2\begin{equation} % \label{eq:LandauLifshitz_Quantum_32.12} \frac{d^2 \chi}{dr^2} + \frac{2m}{\hbar^2} [E-U_l(r)] \chi = 0, ~~\mathrm{where}~U_l(r) = U(r) + \frac{\hbar^2}{2m} \frac{l(l+1)}{r^2} \end{equation}

1次元運動で離散スペクトルをもつとき,エネルギー固有値に対して波動関数は縮退していないので,(24)でllが決まっていて,エネルギーが与えられれば波動関数が決定される.周方向の波動関数はパラメタllmmで決定されたので,最終的に,中心場での粒子の運動を表す波動関数はE,l,mE, l, mによって完全に決定されることがわかる.

動径方向のシュレーディンガー方程式が解けたとして,離散スペクトルの場合についてとりうるエネルギー固有値を並べて,小さいほうから番号をnr=0,1,2,n_r =0,1,2, \cdotsと振ってやろう.この波動関数を決定する3つのパラメタは,nrn_rが動径量子数(radial quantum number),llは方位量子数(azimuthal quantum number),mmは磁気量子数(magnetic quantum number)と呼ばれる.