tecktom.com

この章のまとめ

2026年10月15日10分で読めますAI生成教科書『虎の巻』第4部 第7章 WKB と Hamilton–Jacobi

AI 生成・木下が未勉強

解析力学から様々な理論分野を組み立てる 第4部 第7章 §4.7.7 ★★★ WKB は Hamilton–Jacobi 方程式の解に振幅と位相の数え方を足したもの。量子化は作用変数を ℏ\hbar 単位で数えること、トンネルは虚数の時間の作用

この章では、1粒子の Hamiltonian H=p2/(2m)+V(x)H=p^2/(2m)+V(x) をそのまま量子化した Schrödinger 方程式を、位相 σ\sigma の厳密な方程式 (4.7.4) に書き直し、ℏ/i\hbar/\mathrm{i} の冪で展開した。最低次 (σ0′)2=p2(\sigma_0')^2=p^2 は時間に依らない Hamilton–Jacobi 方程式で、位相は特性関数 W(x,E)W(x,E)、時間を含めれば主関数 S=W−EtS=W-Et である。次の次数は確率の流れの保存で、振幅の2乗は古典の粒子の滞在時間の密度 1/p1/p になった。WKB の波は、1本の軌道でなく、等エネルギーの曲線 ΛE\Lambda_E 全体に対応している。

転回点での振幅の発散は、滑らかな曲線 ΛE\Lambda_E を xx 軸に射影したときだけの特異点で、運動量表示へ移ると消えた。一様な力では運動量表示の WKB が厳密で、それを Fourier 変換すると Airy 関数になる。Fourier 変換を停留位相で評価すると、位相の関数は Legendre 変換 W=W~+pxW=\tilde W+px で移り、2つの枝の寄与に ±π/4\pm\pi/4 の位相が付いた。転回点を通るごとに π/2\pi/2 ずつ遅れるこの規則を1周分数えたものが Maslov の指数 μM\mu_{\mathrm M} で、Bohr–Sommerfeld の条件は作用変数の量子化 J=ℏ(n+μM/4)J=\hbar(n+\mu_{\mathrm M}/4) になった。準位の間隔 ℏω(J)\hbar\omega(J) は作用角変数の ω=∂H/∂J\omega=\partial H/\partial J そのものである。古典的に入れない領域では運動量が虚数になり、トンネルの確率の指数 2Θ2\Theta は、逆さにしたポテンシャルの中の虚数の時間の運動の作用で決まった。

第1部で用意した道具との対応をまとめると、WKB の位相は Hamilton–Jacobi 方程式(第13章)の解であり、量子化の条件は作用変数(第14章)で書け、位置表示と運動量表示のあいだの移り変わりは Legendre 変換(第7章)と母関数(第12章)の規則に従う。次の第4部第8章では、位相 eiS/ℏ\mathrm{e}^{\mathrm{i} S/\hbar} を1本の軌道でなくすべての経路について足し合わせる経路積分を扱い、この章の停留位相がその ℏ→0\hbar\to0 の評価として、虚数の時間の作用がトンネルの経路として、もう一度現れる。

総合問題 重力の中で跳ねる中性子 ★★★

オリジナル質量 mm の粒子が、高さ x>0x>0 で重力のポテンシャル V=mgxV=mgx を受け、x=0x=0 の水平な鏡(硬い壁、ψ(0)=0\psi(0)=0)で跳ね返る。エネルギーを EE とする。

  1. (1) 作用変数 J(E)J(E) を求め、周期 T(E)=2π dJ/dET(E)=2\pi\,\mathrm{d} J/\mathrm{d} E が、速さ v0=2E/mv_0=\sqrt{2E/m} で鏡から跳ね上がって戻るまでの時間 2v0/g2v_0/g に等しいことを確かめよ。(使う道具:第1部第13章の T=d/dE∮p dqT=\mathrm{d}/\mathrm{d} E\oint p\,\mathrm{d} q・第1部第14章の作用変数)
  2. (2) Maslov の指数を数え、WKB の準位を EnWKB=(3π2(n+34))2/3EgE_n^{\mathrm{WKB}}=\bigl(\frac{3\pi}{2}(n+\frac34)\bigr)^{2/3}E_{\mathrm g}、Eg=(ℏ2mg2/2)1/3E_{\mathrm g}=(\hbar^2mg^2/2)^{1/3} の形に求めよ。(使う道具:4.7.3節・4.7.4節)
  3. (3) 厳密な準位を求めよ。x>0x>0 で x→∞x\to\infty に減衰する解は式 (4.7.15) の Ai\mathrm{Ai} であることを使い、Ai\mathrm{Ai} の零点 a1=−2.3381a_1=-2.3381, a2=−4.0879a_2=-4.0879, a3=−5.5206a_3=-5.5206, a4=−6.7867a_4=-6.7867 で表して、n=0,1,2,3n=0,1,2,3 で (2) と比べよ。(使う道具:4.7.3節の運動量表示の厳密解)
  4. (4) 中性子(m=1.675×10−27 kgm=1.675\times10^{-27}\,\mathrm{kg}、g=9.81 m s−2g=9.81\,\mathrm{m\,s^{-2}}、ℏ=1.055×10−34 J s\hbar=1.055\times10^{-34}\,\mathrm{J\,s})で、EgE_{\mathrm g}、基底状態のエネルギー、基底状態のエネルギーの古典の粒子が上がれる最高点の高さを求めよ。第1部第13章の総合問題で求めた「古典極限が破れる高さ」v∗2/(2g)≈3.7 μmv_*^2/(2g)\approx3.7\,\mathrm{\mu m}(v∗=(ℏg/m)1/3v_*=(\hbar g/m)^{1/3})と比べ、この状態を古典的に扱えない理由を述べよ。(使う道具:第1部第13章の総合問題・4.7.1節)
  5. (5) 厳密な準位の n=10n=10 と n=11n=11 の間隔を、式 (4.7.19) の 2πℏ/T(E)2\pi\hbar/T(E)(EE は2つの準位の中点)と比べよ。Ai\mathrm{Ai} の零点は a11=−13.6914a_{11}=-13.6914, a12=−14.5278a_{12}=-14.5278 とする。(使う道具:4.7.4節)
解答

(1) 粒子は 0<x<E/(mg)0<x<E/(mg) を往復する。∮p dx=2∫0E/(mg)2m(E−mgx) dx=423E3/2gm\oint p\,\mathrm{d} x=2\int_0^{E/(mg)}\sqrt{2m(E-mgx)}\,\mathrm{d} x=\dfrac{4\sqrt2}{3}\dfrac{E^{3/2}}{g\sqrt m} なので

J(E)=223π E3/2gm,T(E)=2πdJdE=22 E1/2gm=2g2Em=2v0g.J(E)=\frac{2\sqrt2}{3\pi}\,\frac{E^{3/2}}{g\sqrt m},\qquad T(E)=2\pi\frac{\mathrm{d} J}{\mathrm{d} E}=\frac{2\sqrt2\,E^{1/2}}{g\sqrt m}=\frac{2}{g}\sqrt{\frac{2E}{m}}=\frac{2v_0}{g}.

上がる時間 v0/gv_0/g と落ちる時間 v0/gv_0/g の和に一致する(検算 scripts/verify_ch4-07.py)。

(2) 1周で滑らかな転回点(最高点 x=E/(mg)x=E/(mg))を1回、壁(鏡)を1回通るので μM=1+2=3\mu_{\mathrm M}=1+2=3。J(E)=ℏ(n+34)J(E)=\hbar(n+\frac34) を EE について解くと

EnWKB=(3πℏgm22(n+34))2/3=(3π2(n+34))2/3(ℏ2mg22)1/3.E_n^{\mathrm{WKB}}=\Bigl(\frac{3\pi\hbar g\sqrt m}{2\sqrt2}\Bigl(n+\frac34\Bigr)\Bigr)^{2/3}=\Bigl(\frac{3\pi}{2}\Bigl(n+\frac34\Bigr)\Bigr)^{2/3}\Bigl(\frac{\hbar^2mg^2}{2}\Bigr)^{1/3}.

(3) 一様な力 F=mgF=mg の厳密解は、式 (4.7.15) で z=(x−E/(mg))/ℓgz=(x-E/(mg))/\ell_{\mathrm g}、ℓg=(2m2g/ℏ2)−1/3=(ℏ2/(2m2g))1/3\ell_{\mathrm g}=(2m^2g/\hbar^2)^{-1/3}=(\hbar^2/(2m^2g))^{1/3} と置いた Ai(z)\mathrm{Ai}(z) である。鏡の条件 ψ(0)=0\psi(0)=0 は Ai(−E/(mgℓg))=0\mathrm{Ai}(-E/(mg\ell_{\mathrm g}))=0 で、mgℓg=Egmg\ell_{\mathrm g}=E_{\mathrm g} なので、厳密な準位は

En=−an+1 Eg,n=0,1,2,…E_n=-a_{n+1}\,E_{\mathrm g},\qquad n=0,1,2,\dots

である。数値で比べると(単位 EgE_{\mathrm g})、

nn0123
厳密 −an+1-a_{n+1}2.33814.08795.52066.7867
WKB2.32034.08185.51726.7845
相対誤差−0.76%-0.76\%−0.15%-0.15\%−0.062%-0.062\%−0.033%-0.033\%

基底状態でも誤差は 1% 未満である(検算 scripts/verify_ch4-07.py)。鏡の側には転回点が無く、WKB が破れるのは最高点のまわりだけで、しかもそこは一様な力なので停留位相の −π/4-\pi/4 の精度が高い。壁を滑らかな転回点と同じに数えて n+12n+\frac12 を使うと、基底状態は (3π/4)2/3=1.77(3\pi/4)^{2/3}=1.77 となり、24% もずれる。

(4) Eg=((1.055×10−34)2×1.675×10−27×9.812/2)1/3=9.64×10−32 J=0.602 peVE_{\mathrm g}=\bigl((1.055\times10^{-34})^2\times1.675\times10^{-27}\times9.81^2/2\bigr)^{1/3}=9.64\times10^{-32}\,\mathrm J=0.602\,\mathrm{peV}(1 peV=10−12 eV1\,\mathrm{peV}=10^{-12}\,\mathrm{eV})。基底状態は E0=2.338 Eg=1.41 peVE_0=2.338\,E_{\mathrm g}=1.41\,\mathrm{peV} で、古典の粒子が上がれる最高点は E0/(mg)=13.7 μmE_0/(mg)=13.7\,\mathrm{\mu m} である。長さの単位 ℓg=Eg/(mg)=5.87 μm\ell_{\mathrm g}=E_{\mathrm g}/(mg)=5.87\,\mathrm{\mu m} は、第1部第13章の v∗2/(2g)=3.7 μmv_*^2/(2g)=3.7\,\mathrm{\mu m} と同じ組み合わせ ℏ2/3m−2/3g−1/3\hbar^{2/3}m^{-2/3}g^{-1/3} でできていて、比は ℓg/(v∗2/(2g))=22/3=1.59\ell_{\mathrm g}/(v_*^2/(2g))=2^{2/3}=1.59 である(検算 scripts/verify_ch4-07.py)。最高点 13.7 μm13.7\,\mathrm{\mu m} は 3.7 μm3.7\,\mathrm{\mu m} の 3.73.7 倍しかなく、鏡から最高点までの全体が古典極限の破れる長さの数倍に収まっている。4.7.1節の例題の言葉では、補正の目安 ℏg/(mv3)\hbar g/(mv^3) は鏡の位置(速さ v0=16 mm s−1v_0=16\,\mathrm{mm\,s^{-1}})でも 0.140.14 あり、最高点から 3.7 μm3.7\,\mathrm{\mu m} 下より上(高さの約 27%)では 1 を超える。したがって、この状態は「鏡の上で跳ねる小さな粒」としては扱えず、波として扱うほかない。このような中性子の量子状態は、実際に鏡の上の中性子の透過を測る実験で観測されている1V. V. Nesvizhevsky et al. (2002), Nature 415, 297.下の脚注へ。WKB の準位が (3) のように 1% 以内で合うのは、この状況でも位相の数え方(Maslov の指数)が正しく、誤差が転回点の近くの層にしか生じないためである。

(5) 厳密な間隔は (−a12)−(−a11)=14.5278−13.6914=0.8363 Eg(-a_{12})-(-a_{11})=14.5278-13.6914=0.8363\,E_{\mathrm g}(零点の値は5桁目で丸めてある)。一方、(1) の T(E)T(E) から 2πℏ/T(E)=πℏgm/2E2\pi\hbar/T(E)=\pi\hbar g\sqrt m/\sqrt{2E} で、ℏgm=2 Eg3/2\hbar g\sqrt m=\sqrt2\,E_{\mathrm g}^{3/2} を使うと2πℏ/T(E)=πEg/E/Eg2\pi\hbar/T(E)=\pi E_{\mathrm g}/\sqrt{E/E_{\mathrm g}}。中点 E=14.110 EgE=14.110\,E_{\mathrm g} を入れると 0.8364 Eg0.8364\,E_{\mathrm g} で、厳密な間隔と4桁目まで一致する(検算 scripts/verify_ch4-07.py)。準位の間隔はエネルギーとともに E−1/2E^{-1/2} で狭くなり、これは高く跳ねるほど周期 T∝E1/2T\propto E^{1/2} が長くなることの量子版である。

1 V. V. Nesvizhevsky et al. (2002), Nature 415, 297.

木下の勉強ステータス

木下はまだこの記事を勉強していません。

AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com

← AI生成教科書『虎の巻』の記事一覧へ