[量子力学] 準古典的な系

April 25, 2020      

準古典的な系の波動関数

多数の粒子からなる準古典的な系について議論しよう.シュレーディンガー方程式を以下のように表す.

a2maΔaψ+(EU)ψ=0\begin{equation} \sum_a \frac{\hbar^2}{m_a} \Delta_a \psi + (E - U) \psi = 0 \end{equation}

これに波動関数ψ=eiS/\psi = e^{iS/\hbar}を代入する.

a1ma(aS)2ai2maΔaS=EU\begin{equation} % \label{eq:LandauLifshitz_Quantum_46.2} \sum_a \frac{1}{m_a}(\nabla_a S)^2 - \sum_a \frac{i\hbar}{2m_a} \Delta_a S = E - U \end{equation}

ただし波動関数のラプラシアンは以下のように整理した.以前の議論で「SSは作用(action)でラグランジアンの汎関数である」ということを仮定したが,今の時点ではそのことを前提としていないことに注意しよう.

Δψ=ΔeiS=eiS=[xeiSyeiSzeiS]=[iSxeiSiSyeiSiSzeiS]=i[2Sx2eiS+ i(Sx)2eiS]+i[2Sy2eiS+ i(Sy)2eiS]+i[2Sz2eiS+ i(Sz)2eiS]=[iΔS12(SS)]eiS\begin{align*} &\Delta \psi = \Delta e^{\frac{i}{\hbar}S} = \nabla \cdot \nabla e^{\frac{i}{\hbar}S} = \nabla \cdot \left[ \begin{array}{l} \frac{\partial}{\partial x} e^{\frac{i}{\hbar}S} \\ \frac{\partial}{\partial y} e^{\frac{i}{\hbar}S} \\ \frac{\partial}{\partial z} e^{\frac{i}{\hbar}S} \end{array} \right] = \nabla \cdot \left[ \begin{array}{l} \frac{i}{\hbar} \frac{\partial S}{\partial x} e^{\frac{i}{\hbar}S} \\ \frac{i}{\hbar} \frac{\partial S}{\partial y} e^{\frac{i}{\hbar}S} \\ \frac{i}{\hbar} \frac{\partial S}{\partial z} e^{\frac{i}{\hbar}S} \end{array} \right] \\ &= \frac{i}{\hbar} \left[ \frac{\partial^2 S}{\partial x^2} e^{\frac{i}{\hbar}S} +  \frac{i}{\hbar} \left( \frac{\partial S}{\partial x} \right)^2 e^{\frac{i}{\hbar}S} \right] + \frac{i}{\hbar} \left[ \frac{\partial^2 S}{\partial y^2} e^{\frac{i}{\hbar}S} +  \frac{i}{\hbar} \left( \frac{\partial S}{\partial y} \right)^2 e^{\frac{i}{\hbar}S} \right] + \frac{i}{\hbar} \left[ \frac{\partial^2 S}{\partial z^2} e^{\frac{i}{\hbar}S} +  \frac{i}{\hbar} \left( \frac{\partial S}{\partial z} \right)^2 e^{\frac{i}{\hbar}S} \right] \\ &= \left[ \frac{i}{\hbar} \Delta S - \frac{1}{\hbar^2} (\nabla S \cdot \nabla S) \right] e^{\frac{i}{\hbar}S} \end{align*}

ここでSSとして,\hbarより十分大きい項,\hbar程度の項,2\hbar^2程度の項,…の和という形を仮定する.

S=S0+iS1+(i)2S2+\begin{equation} % \label{eq:LandauLifshitz_Quantum_46.3} S = S_0 + \frac{\hbar}{i} S_1 + \left( \frac{\hbar}{i} \right)^2 S_2 + \cdots \end{equation}

1次元の場合に限れば以下のように書き直せる.

S22miS2m=EU(x)\begin{equation} % \label{eq:LandauLifshitz_Quantum_46.4} \frac{S'^2}{2m} - i \hbar \frac{S''}{2m} = E - U(x) \end{equation}

最初の近似として\hbarについて1次以上の項を無視すると,次の関係が得られる.

S022m=EU(x)S0=±2m[EU(x)]\begin{gather} \frac{{S'_0}^2}{2m} = E - U(x) \\ S'_0 = \pm \sqrt{2m [E - U(x)]} \end{gather}

積分の中身は運動量p(x)p(x)になっている(2mK=mv\sqrt{2m K} = mv).古典的な系であればEU(x)E \le U(x)の範囲での運動しか実現しないので,ルートの中は正である.

S0=±2m[EU(x)]dx=±pdx,  p=2m[EU(x)]\begin{align} S_0 &= \pm \int \sqrt{2m [E - U(x)]} dx \\ &= \pm \int p dx, ~~ p = \sqrt{2m[E - U(x)]} \end{align}

この関係が成り立つためには,(4)の2項目が1項目に比べて十分小さくなければならない.

ddx(S) 1,  where ddx(S)=SS2\begin{gather} % \label{eq:LandauLifshitz_Quantum_46.6} \left| \frac{d}{dx} \left( \frac{\hbar}{S'} \right) \right|  \ll 1, ~~ \mathrm{where}~ \frac{d}{dx} \left( \frac{\hbar}{S'} \right) = - \frac{\hbar S''}{S'^2} \end{gather}

S=pS' = pと代入してやれば(9)は次のように書き換えられる.

mFp31\begin{equation} % \label{eq:LandauLifshitz_Quantum_46.7} \frac{m \hbar |F|}{p^3} \ll 1 \end{equation}

ただし変形には以下の関係を用いた.

dpdx=ddx2m[EU(x)]=122m[EU(x)]ddx(2m[EU(x)])=mpdU(x)dx=mFp\begin{align*} \frac{dp}{dx} &= \frac{d}{dx} \sqrt{2m[E - U(x)]} = \frac{1}{2\sqrt{2m[E - U(x)]}} \frac{d}{dx} \left( 2m[E - U(x)] \right) \\ &= - \frac{m}{p} \frac{dU(x)}{dx} = \frac{mF}{p} \end{align*}

(10)から,運動量が非常に小さい時は近似が成り立たないことがわかる.特にp(x)=0p(x)=0あるいは,E=U(x)E=U(x)となる点においては明らかに不成立である.そこで次に,\hbarの1次の項まで含めた場合について考える.S=S0+iS1S = S_0 + \frac{\hbar}{i} S_1を(4)に代入し,2次以上の項は無視して整理すると以下の通り.

(S0+iS1)2i(S0+iS1)=2m[EU(x)]S02+2iS0S1iS0=2m[EU(x)]2S0S1+S0=0\begin{gather} \left( S'_0 + \frac{\hbar}{i} S'_1 \right)^2 - i\hbar \left( S''_0 + \frac{\hbar}{i} S''_1 \right) = 2m [E - U(x)] \\ {S'_0}^2 + \frac{2\hbar}{i} S'_0 S'_1 - i\hbar S''_0 = 2m [E - U(x)] \\ 2 S'_0 S'_1 + S''_0 = 0 \end{gather}
S1=S02S0S1=12log(S0)=12logp\begin{gather} S'_1 = -\frac{S''_0}{2 S'_0} \\ S_1 = -\frac{1}{2} \log (S'_0) = -\frac{1}{2} \log p \end{gather}

よって\hbarの1次の項まで考慮したSSは以下の通り.

S=±pdx12logp\begin{equation} S = \pm \int p dx - \frac{1}{2} \log p \end{equation}

これをもともと仮定した波動関数ψ=eiS/\psi = e^{iS/\hbar}に代入すると,以下の形で一般解を表すことができる.

ψ=c1peipdx+c2peipdx\begin{equation} % \label{eq:LandauLifshitz_Quantum_46.9} \psi = \frac{c_1}{\sqrt{p}} e^{\frac{i}{\hbar} \int p dx} + \frac{c_2}{\sqrt{p}} e^{-\frac{i}{\hbar} \int p dx} \end{equation}
eiS=exp[i(±pdx12logp ×i)]=exp[±ipdx]×exp[12logp]\begin{equation} e^{\frac{i}{\hbar}S} = \exp \left[ \frac{i}{\hbar} \left( \pm \int p dx - \frac{1}{2} \log p  \times \frac{\hbar}{i} \right) \right] = \exp \left[ \pm \frac{i}{\hbar} \int p dx \right] \times \exp \left[ - \frac{1}{2} \log p \right] \end{equation}

一方で古典的には生じないE<U(x)E < U(x)の範囲については,以下のように表すことができる.

ψ=c1pe1pdx+c2pe1pdx\begin{equation} % \label{eq:LandauLifshitz_Quantum_46.10} \psi = \frac{c_1}{\sqrt{|p|}} e^{\frac{1}{\hbar} \int |p| dx} + \frac{c_2}{\sqrt{|p|}} e^{-\frac{1}{\hbar} \int |p| dx} \end{equation}

準古典的な系の境界条件

x=ax=aU(a)=EU(a)=Eとなるポテンシャルの境界とする.このとき境界の右側と左側での解は次のように表される.

ψ=cpexp[1xapdx],  for x>a\begin{equation} % \label{eq:LandauLifshitz_Quantum_47.1} \psi = \frac{c}{\sqrt{|p|}} \exp\left[ -\frac{1}{\hbar} \int_x^a |p| dx \right], ~~ \mathrm{for}~ x>a \end{equation}
ψ=c1pexp[iaxpdx]+c2pexp[iaxpdx],   for x<a\begin{equation} % \label{eq:LandauLifshitz_Quantum_47.2} \psi = \frac{c_1}{\sqrt{p}} \exp \left[\frac{i}{\hbar} \int_a^x p dx \right] + \frac{c_2}{\sqrt{p}} \exp \left[-\frac{i}{\hbar} \int_a^x p dx \right],  ~~ \mathrm{for}~ x<a \end{equation}

これらの係数を決定するために,境界x=ax=aで接続するための条件を考えよう.まず,x=ax=a付近でのポテンシャルを次のように近似する.

EU(x)=F0(xa),  where F0=dUdxx=a<0\begin{equation} % \label{eq:LandauLifshitz_Quantum_47.3} E - U(x) = F_0(x-a),~~ \mathrm{where}~ F_0 = - \left.\frac{dU}{dx}\right|_{x=a} < 0 \end{equation}

これをp=2m[EU(x)]p = \sqrt{2m[E -U(x)]}に代入すると,運動量の積分は次のように計算できる.

1axp dx=1ax2m[EU(x)]dx=1ax2mF0(xa)12dx=2mF023[(xa)32]ax=232mF0(xa)32\begin{align} &\frac{1}{\hbar} \int_a^x p~dx = \frac{1}{\hbar} \int_a^x \sqrt{2m[E -U(x)]} dx = \frac{1}{\hbar} \int_a^x \sqrt{2mF_0} (x-a)^\frac{1}{2} dx \\ &= \frac{\sqrt{2mF_0}}{\hbar} \frac{2}{3} \left[ (x-a)^\frac{3}{2} \right]^x_a = \frac{2}{3\hbar} \sqrt{2mF_0} (x-a)^\frac{3}{2} \end{align}

ボーア・ゾンマーフェルトの量子化条件

ψ=cpcos[1bxpdxπ4]right side of x=bψ=cpcos[1xapdxπ4] left side of x=a\begin{align} % \label{eq:LandauLifshitz_Quantum_48.1} &\psi = \frac{c}{\sqrt{p}} \cos \left[ \frac{1}{\hbar} \int_b^x p dx - \frac{\pi}{4} \right] && \mathrm{right~side~of~} x=b \\ &\psi = \frac{c'}{\sqrt{p}} \cos \left[ \frac{1}{\hbar} \int_x^a p dx - \frac{\pi}{4} \right] && \mathrm{left~side~of~} x=a \end{align}

左側を基準に考えれば,x=bx=bで位相がπ4-\frac{\pi}{4}で,そこから1bxpdx\frac{1}{\hbar} \int_b^x pdxに従って位相が増えていく.c=cc=c'の場合,x=ax=aでの位相が±π4+2nπ\pm \frac{\pi}{4} + 2n\piになっていればx=ax=aでのcos\cosの値は一致するが,その他の点でも一致するためにはπ4+2nπ\frac{\pi}{4} + 2n\piでなければならない.結局(25, 26)の2式が一致するためには,位相について次の関係が成り立つ必要がある.

if c=c:1bapdxπ4=π4+2nπif c=c:1bapdxπ4=π4+π+2nπ\begin{align} &\mathrm{if~} c=c': &&\frac{1}{\hbar} \int^a_b p dx - \frac{\pi}{4} = \frac{\pi}{4} + 2n \pi \\ &\mathrm{if~} c=-c': &&\frac{1}{\hbar} \int^a_b p dx - \frac{\pi}{4} = \frac{\pi}{4}+ \pi + 2n \pi \end{align}

まとめて書けば次のように表される.

1bapdxπ2=nπ,  with c=(1)nc\begin{equation} \frac{1}{\hbar} \int^a_b p dx - \frac{\pi}{2} = n \pi, ~~\mathrm{with}~ c = (-1)^n c' \end{equation}

周期的な運動の1周分の積分をpdx=2bapdx\oint p dx = 2 \int_b^a p dxとして書き直すと,(30)の関係が得られる.これはBohr-Sommerfeldの量子化条件に対応する.

12πpdx=n+12\begin{equation} % \label{eq:LandauLifshitz_Quantum_48.2} \frac{1}{2\pi \hbar} \oint p dx = n + \frac{1}{2} \end{equation}

さてpdx\oint pdxという積分は,xp平面上の面積と考えることもできる.この領域をn分割すると,分割された領域ひとつひとつの面積は(およそ)2π2\pi\hbarとなる.一方でnという数は系のエネルギー状態に番号を付けたもの(量子数)に対応する.つまりひとつのエネルギー状態はxp平面上の2π2\pi\hbarの面積に対応している.さらに言い換えれば,xp平面上のΔpΔx\Delta p \Delta xの面積には(31)だけのエネルギー状態が対応している,と考えられる.

ΔpΔx2π\begin{equation} % \label{eq:LandauLifshitz_Quantum_48.4} \frac{\Delta p \Delta x}{2\pi \hbar} \end{equation}