tecktom.com

この章のまとめ

2026年10月17日12分で読めますAI生成教科書『虎の巻』第3部 第4章 停留位相と古典極限

AI 生成・木下が未勉強

量子力学の様々な構成方法 第3部 第4章 §3.4.7 ★★ 停留位相で経路積分から古典力学が出てくる。前因子は軌道の束の広がり、位相の飛びは焦点

経路積分の ℏ→0\hbar\to0 は、多重積分の停留位相である。停留条件は離散の Newton 方程式で、N→∞N\to\infty で作用の停留 δS=0\delta S=0 になり、古典軌道だけが残る。停留点のまわりの Gauss 積分の行列式は、Gel'fand–Yaglom の漸化式により Jacobi 方程式の解 f(t)=∂xf/∂vif(t)=\partial x_f/\partial v_i に移り、前因子は ∣∂2Scl/∂xf∂xi∣1/2\lvert \partial^2S_{\mathrm{cl}}/\partial x_f\partial x_i\rvert ^{1/2} になる(Van Vleck の公式 (3.4.14))。軌道の束が1点に集まる焦点を越えるたびに、Hessian の固有値が1つ負になり、位相が π/2\pi/2 遅れる。始点と終点を結ぶ古典軌道が複数あれば、寄与を位相つきで足し、∣K∣2\lvert K\rvert ^2 に干渉が現れる。

基準の側では、AeiS/ℏA\mathrm{e}^{\mathrm{i} S/\hbar} を Schrödinger 方程式に入れた ℏ0\hbar^0・ℏ1\hbar^1 の式が Hamilton–Jacobi 方程式と輸送方程式で、Van Vleck の公式は両方を満たす(命題 3.4.1)。点源の初期条件のもとで、2つの手続きは同じ関数を与える(命題 3.4.2)。時間を Fourier 変換すると WKB の振幅 1/p1/\sqrt p と位相 ∫p dx\int p\,\mathrm{d} x が出る。作用が2次式の系(自由粒子・一様な力・調和振動子・逆さの振動子)では、どの構成の古典極限も厳密になる(表 3.4.1)。

要点

  • 停留位相:∫g eiφ/ℏ≈∑2πℏ/∣φ′′∣ g eiφ/ℏ+iπsgn⁡φ′′/4\int g\,\mathrm{e}^{\mathrm{i}\varphi/\hbar}\approx\sum\sqrt{2\pi\hbar/\lvert \varphi''\rvert }\,g\,\mathrm{e}^{\mathrm{i}\varphi/\hbar+\mathrm{i}\pi\operatorname{sgn}\varphi''/4}。
  • 経路積分の停留点=古典軌道。揺らぎの行列式 ϵdet⁡M→f(t)\epsilon\det M\to f(t)(Jacobi 方程式 mf¨=−V′′fm\ddot f=-V''f、f(0)=0f(0)=0、f˙(0)=1\dot f(0)=1)。
  • Van Vleck:K≈∑(2πiℏ)−1/2∣∂2Scl/∂xf∂xi∣1/2eiScl/ℏ−iπν/2K\approx\sum(2\pi\mathrm{i}\hbar)^{-1/2}\lvert \partial^2S_{\mathrm{cl}}/\partial x_f\partial x_i\rvert ^{1/2}\mathrm{e}^{\mathrm{i} S_{\mathrm{cl}}/\hbar-\mathrm{i}\pi\nu/2}。
  • 位相 SclS_{\mathrm{cl}} は Hamilton–Jacobi 方程式、前因子の2乗は連続の式を満たす。残りは量子ポテンシャル QQ。
  • 時間の Fourier 変換で G≈miℏppiei∫p dx/ℏG\approx\frac{m}{\mathrm{i}\hbar\sqrt{pp_i}}\mathrm{e}^{\mathrm{i}\int p\,\mathrm{d} x/\hbar}(WKB)。

この章の限界も書いておく。焦点の近くと転回点の近くでは式 (3.4.14) が発散し、一様な近似が要る。古典軌道が無い所(エネルギーの足りない障壁の向こう)では、停留点が実の経路でなく虚数の時間の経路になり、トンネル効果とインスタントンとして第3部第7章(第3部第7章)で扱う。虚時間の経路積分の停留位相(鞍点法)は第3部第6章(第3部第6章)で統計力学に使う。

総合問題 調和振動子を Van Vleck の公式だけで解き直す ★★★

オリジナル質量 mm、角周波数 ω\omega の調和振動子 V=12mω2x2V=\frac12m\omega^2x^2(共通例題①)を、停留位相の道具だけで扱う。

  1. (1) xix_i から時間 tt(ωt\omega t は π\pi の整数倍でない)で xfx_f に着く古典軌道と作用 SclS_{\mathrm{cl}} を求めよ。
    (使う道具:3.4.2節・第3部第3章)
  2. (2) Jacobi 方程式を解いて前因子を求め、焦点の数 ν\nu を ωt\omega t で書け。Van Vleck の公式による KK が、0<ωt<π0<\omega t<\pi の Mehler 核を t→t−i0t\to t-\mathrm{i}0 で連続につないだものに一致することを示せ。
    (使う道具:3.4.3節・3.4.4節・第3部第3章)
  3. (3) SclS_{\mathrm{cl}} が Hamilton–Jacobi 方程式を、前因子の2乗が連続の式を満たすことを確かめ、(2) の KK が Schrödinger 方程式を厳密に満たす理由を述べよ。
    (使う道具:3.4.5節・3.4.6節)
  4. (4) 0<ωt<2π0<\omega t<2\pi で跡 Z(t)=∫dx K(x,t;x,0)Z(t)=\int\mathrm{d} x\,K(x,t;x,0) を求め、Z(t)=∑ne−iEnt/ℏZ(t)=\sum_n\mathrm{e}^{-\mathrm{i} E_nt/\hbar} と見て準位 EnE_n を読み取れ。
    (使う道具:(2)・3.4.1節の Fresnel 積分・第2部第4章のスペクトル)
  5. (5) (4) の零点エネルギー ℏω/2\hbar\omega/2 が焦点の位相とどう関わるかを述べよ。さらに 87{}^{87}Rb(m=86.909 um=86.909\,\mathrm u)を ω=2π×50 Hz\omega=2\pi\times50\,\mathrm{Hz} の調和型の罠に入れたとき、最初の焦点の時刻と基底状態の幅 ℏ/(2mω)\sqrt{\hbar/(2m\omega)} を求めよ。
    (使う道具:3.4.4節・(4))
解答

(1) 運動方程式 x¨=−ω2x\ddot x=-\omega^2x の、端の条件を満たす解は xcl(s)=[xisin⁡ω(t−s)+xfsin⁡ωs]/sin⁡ωtx_{\mathrm{cl}}(s)=[x_i\sin\omega(t-s)+x_f\sin\omega s]/\sin\omega t である(sin⁡ωt≠0\sin\omega t\neq0 なので一意)。作用は、運動方程式を使って Scl=m2[xx˙]0tS_{\mathrm{cl}}=\frac m2[x\dot x]_0^t と書け、

Scl(xf,t;xi)=mω2sin⁡ωt[(xf2+xi2)cos⁡ωt−2xfxi]S_{\mathrm{cl}}(x_f,t;x_i)=\frac{m\omega}{2\sin\omega t}\bigl[(x_f^2+x_i^2)\cos\omega t-2x_fx_i\bigr]

となる。

(2) V′′=mω2V''=m\omega^2 なので Jacobi 方程式は f¨=−ω2f\ddot f=-\omega^2f で、f(t)=sin⁡ωt/ωf(t)=\sin\omega t/\omega である。−∂2Scl/∂xf∂xi=mω/sin⁡ωt=m/f(t)-\partial^2S_{\mathrm{cl}}/\partial x_f\partial x_i=m\omega/\sin\omega t=m/f(t) と合う。ff は ωs=π,2π,…\omega s=\pi,2\pi,\dots で 0 になるので ν=⌊ωt/π⌋\nu=\lfloor\omega t/\pi\rfloor である。Van Vleck の公式は

K=12πiℏ∣mωsin⁡ωt∣1/2exp⁡(iℏScl−iπ2⌊ωtπ⌋)K=\frac{1}{\sqrt{2\pi\mathrm{i}\hbar}}\Bigl\lvert\frac{m\omega}{\sin\omega t}\Bigr\rvert^{1/2}\exp\Bigl(\frac{\mathrm{i}}\hbar S_{\mathrm{cl}}-\frac{\mathrm{i}\pi}{2}\Bigl\lfloor\frac{\omega t}{\pi}\Bigr\rfloor\Bigr)

である。一方 Mehler 核の前因子 mω/(2πiℏsin⁡ωt)\sqrt{m\omega/(2\pi\mathrm{i}\hbar\sin\omega t)} を t→t−i0t\to t-\mathrm{i}0 でつなぐと、ωt=nπ\omega t=n\pi の近くで sin⁡ω(t−i0)=sin⁡ωt−i0⋅ωcos⁡nπ\sin\omega(t-\mathrm{i}0)=\sin\omega t-\mathrm{i}0\cdot\omega\cos n\pi となる。nn が奇数なら sin⁡ωt\sin\omega t は正から負へ変わりながら虚部 +0+0 の側(上半平面)を通り、nn が偶数なら負から正へ変わりながら虚部 −0-0 の側(下半平面)を通る。どちらの場合も sin⁡ω(t−i0)\sin\omega(t-\mathrm{i}0) は 0 のまわりを反時計回りに半周するので、偏角が π\pi 増え、1/sin⁡1/\sin の偏角は π\pi 減る。平方根の位相はその半分の −π/2-\pi/2 だけ変わり、nn 回越えると e−iπn/2\mathrm{e}^{-\mathrm{i}\pi n/2} が付く。これは上の式の e−iπν/2\mathrm{e}^{-\mathrm{i}\pi\nu/2} に等しい(ωt=1,4,7,11\omega t=1,4,7,11 で数値でも確かめた)。

(3) ∂tScl+(∂xfScl)2/(2m)+12mω2xf2=0\partial_tS_{\mathrm{cl}}+(\partial_{x_f}S_{\mathrm{cl}})^2/(2m)+\frac12m\omega^2x_f^2=0 が成り立つことは直接の計算で確かめられる(verify)。前因子の2乗 D=mω/sin⁡ωt\mathcal D=m\omega/\sin\omega t は xfx_f によらないので、連続の式は ∂tD+D ∂xfv=0\partial_t\mathcal D+\mathcal D\,\partial_{x_f}v=0 となる。v=∂xfScl/m=ω(xfcos⁡ωt−xi)/sin⁡ωtv=\partial_{x_f}S_{\mathrm{cl}}/m=\omega(x_f\cos\omega t-x_i)/\sin\omega t から ∂xfv=ωcot⁡ωt\partial_{x_f}v=\omega\cot\omega t、∂tD=−mω2cos⁡ωt/sin⁡2ωt=−D ωcot⁡ωt\partial_t\mathcal D=-m\omega^2\cos\omega t/\sin^2\omega t=-\mathcal D\,\omega\cot\omega t で、成り立つ。前因子 AA が xfx_f によらないので、式 (3.4.19) の右辺 ℏ22m(∂x2A/A)K\frac{\hbar^2}{2m}(\partial_x^2A/A)K は 0 で、KK は Schrödinger 方程式を厳密に満たす。作用が2次式であることの帰結である。

(4) 対角 xf=xi=xx_f=x_i=x では Scl=mωx2(cos⁡ωt−1)sin⁡ωt=−mωx2tan⁡ωt2S_{\mathrm{cl}}=\frac{m\omega x^2(\cos\omega t-1)}{\sin\omega t}=-m\omega x^2\tan\frac{\omega t}{2} である。0<ωt<π0<\omega t<\pi では

Z(t)=mω2πiℏsin⁡ωt∫dx e−imωx2tan⁡(ωt/2)/ℏ=mω2πiℏsin⁡ωtπℏimωtan⁡(ωt/2)Z(t)=\sqrt{\frac{m\omega}{2\pi\mathrm{i}\hbar\sin\omega t}}\int\mathrm{d} x\,\mathrm{e}^{-\mathrm{i} m\omega x^2\tan(\omega t/2)/\hbar} =\sqrt{\frac{m\omega}{2\pi\mathrm{i}\hbar\sin\omega t}}\sqrt{\frac{\pi\hbar}{\mathrm{i} m\omega\tan(\omega t/2)}}

で、sin⁡ωttan⁡(ωt/2)=2sin⁡2(ωt/2)\sin\omega t\tan(\omega t/2)=2\sin^2(\omega t/2) を使うと Z(t)=12isin⁡(ωt/2)Z(t)=\frac{1}{2\mathrm{i}\sin(\omega t/2)} を得る。xx の積分は tan⁡(ωt/2)>0\tan(\omega t/2)>0 の Fresnel 積分で、位相 e−iπ/4\mathrm{e}^{-\mathrm{i}\pi/4} を持つ。π<ωt<2π\pi<\omega t<2\pi では KK の位相 e−iπ/2\mathrm{e}^{-\mathrm{i}\pi/2} と、tan⁡(ωt/2)<0\tan(\omega t/2)<0 になった Fresnel 積分の位相の変化 e+iπ/2\mathrm{e}^{+\mathrm{i}\pi/2} が打ち消し、同じ式がそのまま成り立つ(tt に負の虚部を与えた数値積分で確かめた)。q=e−iωtq=\mathrm{e}^{-\mathrm{i}\omega t}(tt に負の小さな虚部を与えて ∣q∣<1\lvert q\rvert <1)とおくと

Z(t)=e−iωt/21−e−iωt=∑n=0∞e−iω(n+12)tZ(t)=\frac{\mathrm{e}^{-\mathrm{i}\omega t/2}}{1-\mathrm{e}^{-\mathrm{i}\omega t}}=\sum_{n=0}^\infty\mathrm{e}^{-\mathrm{i}\omega(n+\frac12)t}

なので、準位は En=ℏω(n+12)E_n=\hbar\omega(n+\frac12) で、どれも縮退していない。共通例題①の正答と同じである。

(5) 跡の xx の積分の停留点は x=0x=0(原点に静止した周期軌道)で、そのまわりの軌道は1周期 2π/ω2\pi/\omega の間に焦点を2つ通り、位相 e−iπ=−1\mathrm{e}^{-\mathrm{i}\pi}=-1 を受け取る。(4) の Z(t)Z(t) は tt が1周期進むと sin⁡(ωt/2)\sin(\omega t/2) の符号が変わって −1-1 倍になり、分子 e−iωt/2\mathrm{e}^{-\mathrm{i}\omega t/2} のこの −1-1 が焦点の位相そのものである。これがなければ準位は ℏωn\hbar\omega n になり、零点エネルギー ℏω/2\hbar\omega/2 は現れない。WKB の Bohr–Sommerfeld の条件 J=(n+12)ℏJ=(n+\frac12)\hbar の 12\frac12 と同じものを、経路の側の焦点の数え方から得たことになる。

87{}^{87}Rb では、最初の焦点は π/ω=1/(2×50) s=10 ms\pi/\omega=1/(2\times50)\,\mathrm s=10\,\mathrm{ms} である。m=86.909×1.6605×10−27 kg=1.443×10−25 kgm=86.909\times1.6605\times10^{-27}\,\mathrm{kg}=1.443\times10^{-25}\,\mathrm{kg} なので、基底状態の幅は 1.055×10−34/(2×1.443×10−25×314.2) m=1.08 μm\sqrt{1.055\times10^{-34}/(2\times1.443\times10^{-25}\times314.2)}\,\mathrm m=1.08\,\mu\mathrm m である。罠の中のある点から出た原子の波は 10 ms 後に反対側の鏡像の点に集まり、その時刻の伝播関数は e−iπ/2δ(xf+xi)\mathrm{e}^{-\mathrm{i}\pi/2}\delta(x_f+x_i) になる。

木下の勉強ステータス

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

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

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