解析力学から様々な理論分野を組み立てる 第4部 第8章 §4.8.4 ★★★ 停留位相の評価を経路の空間で行うと、古典の経路だけが残り、前因子は Van Vleck の行列式、共役点を越えるごとに位相 − π / 2 -\pi/2 − π /2 が付く。最小作用の原理が「最小」でなく「停留」である理由がここにある
1変数の停留位相
大きなパラメータ λ \lambda λ を含む積分 I ( λ ) = ∫ d x g ( x ) e i λ φ ( x ) I(\lambda)=\int\mathrm{d} x\,g(x)\,\mathrm{e}^{\mathrm{i}\lambda\varphi(x)} I ( λ ) = ∫ d x g ( x ) e i λ φ ( x ) を考える。λ \lambda λ が大きいと、φ \varphi φ が変化する所では被積分関数が速く振動し、隣り合う寄与が打ち消し合う。打ち消し合わずに残るのは、φ ′ ( x ∗ ) = 0 \varphi'(x_*)=0 φ ′ ( x ∗ ) = 0 となる停留点 x ∗ x_* x ∗ の近くだけである。そこで φ \varphi φ を2次まで展開し、Fresnel 積分 (4.8.7) を当てると
I ( λ ) ≈ ∑ x ∗ g ( x ∗ ) 2 π λ ∣ φ ′ ′ ( x ∗ ) ∣ exp [ i λ φ ( x ∗ ) + i π 4 sgn φ ′ ′ ( x ∗ ) ] (4.8.24) I(\lambda)\approx\sum_{x_*}g(x_*)\sqrt{\frac{2\pi}{\lambda\lvert \varphi''(x_*)\rvert }}\;\exp\Bigl[\mathrm{i}\lambda\varphi(x_*)+\frac{\mathrm{i}\pi}{4}\operatorname{sgn}\varphi''(x_*)\Bigr]\tag{4.8.24} I ( λ ) ≈ x ∗ ∑ g ( x ∗ ) λ ∣ φ ′′ ( x ∗ )∣ 2 π exp [ i λ φ ( x ∗ ) + 4 i π sgn φ ′′ ( x ∗ ) ] ( 4.8.24 )
を得る。第4部第7章では、この公式を運動量表示から位置表示への Fourier 変換に当てて、転回点の接続の公式を出した(第4部第7章 )。各停留点の寄与は λ − 1 / 2 \lambda^{-1/2} λ − 1/2 の大きさで、停留点から離れた所の寄与はそれより速く小さくなる。
経路の空間での停留位相
経路積分 (4.8.16) は、N − 1 N-1 N − 1 個の変数 q 1 , … , q N − 1 q_1,\dots,q_{N-1} q 1 , … , q N − 1 についての積分で、大きなパラメータは 1 / ℏ 1/\hbar 1/ℏ 、位相の関数は S N ( q 1 , … , q N − 1 ) S_N(q_1,\dots,q_{N-1}) S N ( q 1 , … , q N − 1 ) である。ℏ \hbar ℏ が作用に比べて小さいとき、式 (4.8.24) を多変数で使うと、寄与は ∂ S N / ∂ q k = 0 \partial S_N/\partial q_k=0 ∂ S N / ∂ q k = 0 (k = 1 , … , N − 1 k=1,\dots,N-1 k = 1 , … , N − 1 )を満たす経路の近くから来る。N → ∞ N\to\infty N → ∞ でこの条件は δ S = 0 \delta S=0 δ S = 0 、すなわち端を固定した Euler–Lagrange 方程式になる。最小作用の原理は、量子力学の経路の和で位相がそろう条件として出てくる。第1部第5章で、古典力学の中には作用を停留させる理由が無く、理由は量子力学にあると述べたのは、この意味である(第1部第5章 )。停留点が最小点か鞍点かは問われない。式 (4.8.24) は φ ′ ′ \varphi'' φ ′′ の符号によらず使え、符号は位相 ± π / 4 \pm\pi/4 ± π /4 にだけ現れる。光学で同じ議論をすると、光路長 ∫ n d s \int n\,\mathrm{d} s ∫ n d s の停留する光線だけが残る Fermat の原理になる(第2部第10章 )。
停留点のまわりの2次の積分は、第3節の揺らぎの積分そのものである。一般の V V V について、古典の経路のまわりで作用を2次まで展開し、その揺らぎの積分を行った結果は
K ( q f , q i ; T ) ≈ ∑ 古典の経路 1 2 π i ℏ ∣ ∂ 2 S c l ∂ q f ∂ q i ∣ exp [ i ℏ S c l ( q f , q i ; T ) − i π 2 ν ] (4.8.25) K(q_f,q_i;T)\approx\sum_{\text{古典の経路}}\sqrt{\frac{1}{2\pi\mathrm{i}\hbar}\Bigl|\frac{\partial^2S_{\mathrm{cl}}}{\partial q_f\,\partial q_i}\Bigr|}\;
\exp\Bigl[\frac{\mathrm{i}}{\hbar}S_{\mathrm{cl}}(q_f,q_i;T)-\frac{\mathrm{i}\pi}{2}\nu\Bigr]\tag{4.8.25} K ( q f , q i ; T ) ≈ 古典の経路 ∑ 2 π i ℏ 1 ∂ q f ∂ q i ∂ 2 S cl exp [ ℏ i S cl ( q f , q i ; T ) − 2 i π ν ] ( 4.8.25 )
になる。和は q i q_i q i から時間 T T T で q f q_f q f に着く古典の経路すべてについて取り、ν \nu ν はその経路が途中で通り過ぎた共役点の数である。この式を Van Vleck の公式と呼ぶ1 1 J. H. Van Vleck (1928), Proceedings of the National Academy of Sciences 14 , 178–188。経路積分からの導出は C. Morette (1951), Physical Review 81 , 848–852。 。2次の L L L では式 (4.8.23) がこの形の厳密な場合で、近似の記号 ≈ \approx ≈ が等号になる。
前因子の意味を見る。第1部第13章で、主関数の微分が端の運動量になることを示した。始点では p i = − ∂ S c l / ∂ q i p_i=-\partial S_{\mathrm{cl}}/\partial q_i p i = − ∂ S cl / ∂ q i である(第1部第13章 )。これを q f q_f q f で微分すると − ∂ 2 S c l / ∂ q f ∂ q i = ∂ p i / ∂ q f = ( ∂ q f / ∂ p i ) − 1 -\partial^2S_{\mathrm{cl}}/\partial q_f\partial q_i=\partial p_i/\partial q_f=(\partial q_f/\partial p_i)^{-1} − ∂ 2 S cl / ∂ q f ∂ q i = ∂ p i / ∂ q f = ( ∂ q f / ∂ p i ) − 1 となる。∂ q f / ∂ p i \partial q_f/\partial p_i ∂ q f / ∂ p i は、始点の運動量を少し変えたときに終点がどれだけ動くかを表し、古典の経路の束がどれだけ広がるかの目安である。経路の束が広がっている所では振幅が小さく、束が1点に集まる共役点(∂ q f / ∂ p i = 0 \partial q_f/\partial p_i=0 ∂ q f / ∂ p i = 0 )では前因子が発散する。これは第4部第7章で、転回点で WKB の振幅が発散したことの、経路の側での姿である。調和振動子では q f = q i cos ω T + p i sin ω T / ( m ω ) q_f=q_i\cos\omega T+p_i\sin\omega T/(m\omega) q f = q i cos ω T + p i sin ω T / ( mω ) なので ∂ q f / ∂ p i = sin ω T / ( m ω ) \partial q_f/\partial p_i=\sin\omega T/(m\omega) ∂ q f / ∂ p i = sin ω T / ( mω ) で、式 (4.8.23) の前因子と合う。
共役点を越えると位相 − π / 2 -\pi/2 − π /2
調和振動子の核を ω T > π \omega T>\pi ω T > π へ延ばす。ω T \omega T ω T が π \pi π を越えると式 (4.8.21) の n = 1 n=1 n = 1 のモードの係数が負になり、そのモードの Fresnel 積分の位相が e i π / 4 \mathrm{e}^{\mathrm{i}\pi/4} e i π /4 から e − i π / 4 \mathrm{e}^{-\mathrm{i}\pi/4} e − i π /4 に変わる。比 F / F 0 F/F_0 F / F 0 には e − i π / 2 \mathrm{e}^{-\mathrm{i}\pi/2} e − i π /2 が掛かる。したがって π < ω T < 2 π \pi<\omega T<2\pi π < ω T < 2 π では
K ( q f , q i ; T ) = m ω 2 π i ℏ ∣ sin ω T ∣ exp [ i ℏ S c l − i π 2 ] (4.8.26) K(q_f,q_i;T)=\sqrt{\frac{m\omega}{2\pi\mathrm{i}\hbar\lvert \sin\omega T\rvert }}\,\exp\Bigl[\frac{\mathrm{i}}{\hbar}S_{\mathrm{cl}}-\frac{\mathrm{i}\pi}{2}\Bigr]\tag{4.8.26} K ( q f , q i ; T ) = 2 π i ℏ ∣ sin ω T ∣ mω exp [ ℏ i S cl − 2 i π ] ( 4.8.26 )
で、式 (4.8.25) の ν = 1 \nu=1 ν = 1 の場合である。一般に ν \nu ν は ω T / π \omega T/\pi ω T / π の整数部分で、図 4.8.2 で行列式が 0 を横切った回数に等しい。
この位相は、ちょうど共役点の時刻 T = π / ω T=\pi/\omega T = π / ω で何が起きるかを見ると確かめることができる。ω T \omega T ω T が π \pi π に下から近づくと sin ω T = s → + 0 \sin\omega T=s\to+0 sin ω T = s → + 0 、cos ω T → − 1 \cos\omega T\to-1 cos ω T → − 1 でS c l → − m ω ( q f + q i ) 2 / ( 2 s ) S_{\mathrm{cl}}\to-m\omega(q_f+q_i)^2/(2s) S cl → − mω ( q f + q i ) 2 / ( 2 s ) となり、核は q f = − q i q_f=-q_i q f = − q i のまわりに幅 ℏ s / ( m ω ) \sqrt{\hbar s/(m\omega)} ℏ s / ( mω ) で鋭く集まる Gauss 型の関数になる。その q f q_f q f についての積分は、式 (4.8.7) で a = − m ω / ( ℏ s ) < 0 a=-m\omega/(\hbar s)<0 a = − mω / ( ℏ s ) < 0 として m ω / ( 2 π i ℏ s ) 2 π ℏ s / ( m ω ) e − i π / 4 = e − i π / 2 \sqrt{m\omega/(2\pi\mathrm{i}\hbar s)}\sqrt{2\pi\hbar s/(m\omega)}\,\mathrm{e}^{-\mathrm{i}\pi/4}=\mathrm{e}^{-\mathrm{i}\pi/2} mω / ( 2 π i ℏ s ) 2 π ℏ s / ( mω ) e − i π /4 = e − i π /2 である。したがって
K ( q f , q i ; π ω ) = e − i π / 2 δ ( q f + q i ) . K\Bigl(q_f,q_i;\frac{\pi}{\omega}\Bigr)=\mathrm{e}^{-\mathrm{i}\pi/2}\,\delta(q_f+q_i). K ( q f , q i ; ω π ) = e − i π /2 δ ( q f + q i ) .
半周期後の状態は始めの状態を q → − q q\to-q q → − q と裏返したものに − i -\mathrm{i} − i を掛けたものになる。固有状態の言葉では、e − i H ^ π / ( ω ℏ ) \mathrm{e}^{-\mathrm{i}\hat H\pi/(\omega\hbar)} e − i H ^ π / ( ω ℏ ) の固有値はe − i π ( n + 1 / 2 ) = − i ( − 1 ) n \mathrm{e}^{-\mathrm{i}\pi(n+1/2)}=-\mathrm{i}(-1)^n e − i π ( n + 1/2 ) = − i ( − 1 ) n で、( − 1 ) n (-1)^n ( − 1 ) n は ψ n \psi_n ψ n の偶奇なので、裏返しに − i -\mathrm{i} − i を掛けた演算子と一致する。共役点で古典の経路が1点に集まるとき、核は δ \delta δ 関数になり、通り過ぎるたびに位相 − π / 2 -\pi/2 − π /2 を残す。第4部第7章で転回点ごとに数えた Maslov の位相 − π / 2 -\pi/2 − π /2 は、1周期の閉じた軌道が転回点の近くで通る共役点の数を数えたものに当たる。
図 4.8.3 例題 4.8.4 の積分。厳密な値 2 π J 0 ( λ ) 2\pi J_0(\lambda) 2 π J 0 ( λ ) (黒の点線)と、2つの停留点 x = 0 , π x=0,\pi x = 0 , π の寄与を足した停留位相の近似 2 2 π / λ cos ( λ − π / 4 ) 2\sqrt{2\pi/\lambda}\cos(\lambda-\pi/4) 2 2 π / λ cos ( λ − π /4 ) (青の実線、λ ≥ 2 \lambda\ge2 λ ≥ 2 )。λ = 3 \lambda=3 λ = 3 ですでに曲線はほとんど重なる。停留点が1つだけでは振動は出ず、2つの停留点の位相 ± ( λ − π / 4 ) \pm(\lambda-\pi/4) ± ( λ − π /4 ) の干渉が振動になる
例題 4.8.4 2つの停留点の干渉 ★★
オリジナル I ( λ ) = ∫ 0 2 π e i λ cos x d x I(\lambda)=\int_0^{2\pi}\mathrm{e}^{\mathrm{i}\lambda\cos x}\,\mathrm{d} x I ( λ ) = ∫ 0 2 π e i λ c o s x d x を考える。厳密には I ( λ ) = 2 π J 0 ( λ ) I(\lambda)=2\pi J_0(\lambda) I ( λ ) = 2 π J 0 ( λ ) (J 0 J_0 J 0 は0次の Bessel 関数)である。
(1) 停留点を求め、式 (4.8.24) で I ( λ ) ≈ 2 2 π / λ cos ( λ − π / 4 ) I(\lambda)\approx2\sqrt{2\pi/\lambda}\,\cos(\lambda-\pi/4) I ( λ ) ≈ 2 2 π / λ cos ( λ − π /4 ) を導け。
(2) λ = 10 \lambda=10 λ = 10 と λ = 30 \lambda=30 λ = 30 で、厳密な値 2 π J 0 ( 10 ) = − 1.5453 2\pi J_0(10)=-1.5453 2 π J 0 ( 10 ) = − 1.5453 、2 π J 0 ( 30 ) = − 0.5427 2\pi J_0(30)=-0.5427 2 π J 0 ( 30 ) = − 0.5427 と比べよ。
解答 (1) φ ( x ) = cos x \varphi(x)=\cos x φ ( x ) = cos x の停留点は φ ′ = − sin x = 0 \varphi'=-\sin x=0 φ ′ = − sin x = 0 から x = 0 x=0 x = 0 と x = π x=\pi x = π である(x = 2 π x=2\pi x = 2 π は x = 0 x=0 x = 0 と同じ点で、周期関数の積分なので1回だけ数える)。x = 0 x=0 x = 0 では φ = 1 \varphi=1 φ = 1 、φ ′ ′ = − 1 \varphi''=-1 φ ′′ = − 1 、x = π x=\pi x = π では φ = − 1 \varphi=-1 φ = − 1 、φ ′ ′ = + 1 \varphi''=+1 φ ′′ = + 1 である。式 (4.8.24) で g = 1 g=1 g = 1 として
I ( λ ) ≈ 2 π λ ( e i λ − i π / 4 + e − i λ + i π / 4 ) = 2 2 π λ cos ( λ − π 4 ) . I(\lambda)\approx\sqrt{\frac{2\pi}{\lambda}}\Bigl(\mathrm{e}^{\mathrm{i}\lambda-\mathrm{i}\pi/4}+\mathrm{e}^{-\mathrm{i}\lambda+\mathrm{i}\pi/4}\Bigr)=2\sqrt{\frac{2\pi}{\lambda}}\cos\Bigl(\lambda-\frac\pi4\Bigr). I ( λ ) ≈ λ 2 π ( e i λ − i π /4 + e − i λ + i π /4 ) = 2 λ 2 π cos ( λ − 4 π ) .
2つの停留点は、経路積分で言えば2本の古典の経路に当たる。それぞれの寄与の大きさは同じで、位相 λ φ ( x ∗ ) \lambda\varphi(x_*) λ φ ( x ∗ ) と ± π / 4 \pm\pi/4 ± π /4 が違うので干渉する。
(2) λ = 10 \lambda=10 λ = 10 では 2 2 π / 10 = 1.5853 2\sqrt{2\pi/10}=1.5853 2 2 π /10 = 1.5853 、cos ( 10 − π / 4 ) = − 0.9780 \cos(10-\pi/4)=-0.9780 cos ( 10 − π /4 ) = − 0.9780 で、近似は − 1.5504 -1.5504 − 1.5504 である。厳密な値 − 1.5453 -1.5453 − 1.5453 との差は 0.3 % 0.3\% 0.3% である。λ = 30 \lambda=30 λ = 30 では近似が − 0.5396 -0.5396 − 0.5396 、厳密な値が − 0.5427 -0.5427 − 0.5427 で、差は 0.003 0.003 0.003 (0.6 % 0.6\% 0.6% )である。差の絶対値は λ \lambda λ とともに小さくなる。λ = 30 \lambda=30 λ = 30 で割合が大きく見えるのは、そこで cos ( λ − π / 4 ) \cos(\lambda-\pi/4) cos ( λ − π /4 ) が小さく、値そのものが小さいためである(図 4.8.3 )。
近似の裏付け
型 B(極限) | 小さいもの ℏ / S c l → 0 \hbar/S_{\mathrm{cl}}\to0 ℏ/ S cl → 0 (作用が ℏ \hbar ℏ に比べて大きい) | 証明 要確認(有限次元の積分の停留位相は、停留点が孤立して非退化なら証明済み。経路積分そのものに当てたときの前提は原論文で未確認。共役点のごく近くでは2次の展開が使えない) | 行き先 近似の数理 第7章
1 J. H. Van Vleck (1928), Proceedings of the National Academy of Sciences 14 , 178–188。経路積分からの導出は C. Morette (1951), Physical Review 81 , 848–852。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com