解析力学から様々な理論分野を組み立てる 第1部 第13章 §1.13.6 ★★ エネルギー保存から時間を積分で出す解き方、Binet の方程式、WKB の公式は、どれも特性関数 W ( q ) W(q) W ( q ) を経由した式の一部
エネルギー保存から時間を出す解き方
1次元の運動を学部の力学で解くときは、まずエネルギー保存 1 2 m q ˙ 2 + V ( q ) = E \frac12m\dot q^2+V(q)=E 2 1 m q ˙ 2 + V ( q ) = E を書き、q ˙ = 2 ( E − V ) / m \dot q=\sqrt{2(E-V)/m} q ˙ = 2 ( E − V ) / m を変数分離して
t − t 0 = ∫ q 0 q d q ′ 2 ( E − V ( q ′ ) ) / m , 周期 T = ∮ d q q ˙ (1.13.23) t-t_0=\int_{q_0}^{q}\frac{\mathrm{d} q'}{\sqrt{2(E-V(q'))/m}},\qquad \text{周期}\ T=\oint\frac{\mathrm{d} q}{\dot q}\tag{1.13.23} t − t 0 = ∫ q 0 q 2 ( E − V ( q ′ )) / m d q ′ , 周期 T = ∮ q ˙ d q ( 1.13.23 )
とする。この章の式 (1.13.11) では、時刻は特性関数をエネルギーで微分した ∂ W / ∂ E = ∫ ( ∂ p / ∂ E ) d q \partial W/\partial E=\int(\partial p/\partial E)\,\mathrm{d} q ∂ W / ∂ E = ∫ ( ∂ p / ∂ E ) d q として出た。∂ p / ∂ E = m / p = 1 / q ˙ \partial p/\partial E=m/p=1/\dot q ∂ p / ∂ E = m / p = 1/ q ˙ なので、被積分関数は同じものである。違いは、この章の形では時刻が「W W W の E E E 微分」という1つの量にまとまっていることで、1周の積分をとると
T = ∮ d q q ˙ = d d E ∮ p d q (1.13.24) T=\oint\frac{\mathrm{d} q}{\dot q}=\frac{\mathrm{d} }{\mathrm{d} E}\oint p\,\mathrm{d} q\tag{1.13.24} T = ∮ q ˙ d q = d E d ∮ p d q ( 1.13.24 )
となる。周期は、相空間で等エネルギー線が囲む面積をエネルギーで微分したものである。調和振動子では式 (1.13.13) を微分して T = 2 π / ω T=2\pi/\omega T = 2 π / ω を得る。この関係は次の章で作用変数 J = 1 2 π ∮ p d q J=\frac{1}{2\pi}\oint p\,\mathrm{d} q J = 2 π 1 ∮ p d q を新しい運動量に選ぶ理由になる。
Kepler 問題
学部の力学では、Kepler 問題は角運動量の保存 m r 2 φ ˙ = α φ mr^2\dot\varphi=\alpha_\varphi m r 2 φ ˙ = α φ で時間微分を φ \varphi φ 微分に替え、u = 1 / r u=1/r u = 1/ r について Binet の方程式
d 2 u d φ 2 + u = m k α φ 2 \frac{\mathrm{d}^2u}{\mathrm{d}\varphi^2}+u=\frac{mk}{\alpha_\varphi^2} d φ 2 d 2 u + u = α φ 2 mk
を導き、それを解いて u = 1 λ ( 1 + e cos ( φ − φ 0 ) ) u=\frac{1}{\lambda}\bigl(1+e\cos(\varphi-\varphi_0)\bigr) u = λ 1 ( 1 + e cos ( φ − φ 0 ) ) を得ることが多い。この章では、同じ結果を1.13.4 節で ∂ W / ∂ α φ = β φ \partial W/\partial\alpha_\varphi=\beta_\varphi ∂ W / ∂ α φ = β φ の積分1つから出した。Binet の方程式は2階の常微分方程式で、積分定数 e e e , φ 0 \varphi_0 φ 0 は解いた後で初期条件から決める。Hamilton–Jacobi の側では、分離の定数 E E E , α φ \alpha_\varphi α φ が最初から式の中にあり、離心率は e = 1 + 2 E α φ 2 / ( m k 2 ) e=\sqrt{1+2E\alpha_\varphi^2/(mk^2)} e = 1 + 2 E α φ 2 / ( m k 2 ) としてそれらで書かれる。
WKB の公式
量子力学の教科書では、時間に依らない1次元の Schrödinger 方程式の近似解として
ψ ( q ) ≈ C p ( q ) exp ( ± i ℏ ∫ q p ( q ′ ) d q ′ ) , p ( q ) = 2 m ( E − V ( q ) ) \psi(q)\approx\frac{C}{\sqrt{p(q)}}\exp\Bigl(\pm\frac{\mathrm{i}}{\hbar}\int^qp(q')\,\mathrm{d} q'\Bigr),\qquad p(q)=\sqrt{2m(E-V(q))} ψ ( q ) ≈ p ( q ) C exp ( ± ℏ i ∫ q p ( q ′ ) d q ′ ) , p ( q ) = 2 m ( E − V ( q ))
を WKB の公式として習う。指数の中の ∫ p d q \int p\,\mathrm{d} q ∫ p d q は特性関数 (1.13.10) そのもので、符号 ± \pm ± は W W W の2つの枝(右向きと左向きの運動)にあたる。前の因子 1 / p 1/\sqrt p 1/ p は、時間に依らない場合の連続の式 (1.13.21) 、すなわち d ( A 2 p / m ) / d q = 0 \mathrm{d}(A^2p/m)/\mathrm{d} q=0 d ( A 2 p / m ) / d q = 0 から出る A ∝ p − 1 / 2 A\propto p^{-1/2} A ∝ p − 1/2 である。古典的には、粒子の速さが遅い所ほど長く留まるので見つかる確率が高い、という事情を表している。
学部の授業での形 この章の形 1次元の時間 t = ∫ d q / 2 ( E − V ) / m t=\int\mathrm{d} q/\sqrt{2(E-V)/m} t = ∫ d q / 2 ( E − V ) / m t + β 1 = ∂ W / ∂ E t+\beta_1=\partial W/\partial E t + β 1 = ∂ W / ∂ E 周期 T = ∮ d q / q ˙ T=\oint\mathrm{d} q/\dot q T = ∮ d q / q ˙ T = d d E ∮ p d q T=\frac{\mathrm{d} }{\mathrm{d} E}\oint p\,\mathrm{d} q T = d E d ∮ p d q Kepler の軌道 Binet の方程式を解く ∂ W / ∂ α φ = β φ \partial W/\partial\alpha_\varphi=\beta_\varphi ∂ W / ∂ α φ = β φ を積分保存量 運動方程式から探す 分離の定数 E E E , α φ \alpha_\varphi α φ , α θ \alpha_\theta α θ WKB の位相と振幅 ∫ p d q / ℏ \int p\,\mathrm{d} q/\hbar ∫ p d q /ℏ と 1 / p 1/\sqrt p 1/ p 特性関数 W W W と連続の式 光線 アイコナール方程式 ∣ ∇ S ∣ 2 = n 2 \lvert \nabla\mathcal S\rvert ^2=n^2 ∣ ∇ S ∣ 2 = n 2 ∣ ∇ W ∣ 2 = 2 m ( E − V ) \lvert \nabla W\rvert ^2=2m(E-V) ∣ ∇ W ∣ 2 = 2 m ( E − V )
例題 1.13.6 4次のポテンシャルの周期を面積から出す ★★
オリジナル 質量 m m m の粒子が V ( q ) = a q 4 V(q)=aq^4 V ( q ) = a q 4 (a > 0 a>0 a > 0 )の中で振動する。
(1) ∮ p d q \oint p\,\mathrm{d} q ∮ p d q を E E E の冪と定数 I = ∫ 0 1 1 − u 4 d u I=\int_0^1\sqrt{1-u^4}\,\mathrm{d} u I = ∫ 0 1 1 − u 4 d u で表せ。
(2) 式 (1.13.24) を使って周期 T T T を求め、T T T が E − 1 / 4 E^{-1/4} E − 1/4 に比例することを示せ。
(3) m = a = E = 1 m=a=E=1 m = a = E = 1 のときの周期の値を求めよ(I ≈ 0.87402 I\approx0.87402 I ≈ 0.87402 )。エネルギーを16倍にすると周期は何倍になるか。
解答 (1) 転回点は q t = ( E / a ) 1 / 4 q_{\mathrm t}=(E/a)^{1/4} q t = ( E / a ) 1/4 で、対称性から ∮ p d q = 4 ∫ 0 q t 2 m ( E − a q 4 ) d q \oint p\,\mathrm{d} q=4\int_0^{q_{\mathrm t}}\sqrt{2m(E-aq^4)}\,\mathrm{d} q ∮ p d q = 4 ∫ 0 q t 2 m ( E − a q 4 ) d q 。q = q t u q=q_{\mathrm t}u q = q t u と置くと a q 4 = E u 4 aq^4=Eu^4 a q 4 = E u 4 なので
∮ p d q = 4 2 m E q t ∫ 0 1 1 − u 4 d u = 4 2 m a − 1 / 4 I E 3 / 4 . \oint p\,\mathrm{d} q=4\sqrt{2mE}\,q_{\mathrm t}\int_0^1\sqrt{1-u^4}\,\mathrm{d} u=4\sqrt{2m}\,a^{-1/4}I\,E^{3/4}. ∮ p d q = 4 2 m E q t ∫ 0 1 1 − u 4 d u = 4 2 m a − 1/4 I E 3/4 .
(2) T = d d E ∮ p d q = 3 4 ⋅ 4 2 m a − 1 / 4 I E − 1 / 4 = 3 2 m a − 1 / 4 I E − 1 / 4 T=\frac{\mathrm{d} }{\mathrm{d} E}\oint p\,\mathrm{d} q=\frac34\cdot4\sqrt{2m}\,a^{-1/4}I\,E^{-1/4}=3\sqrt{2m}\,a^{-1/4}I\,E^{-1/4} T = d E d ∮ p d q = 4 3 ⋅ 4 2 m a − 1/4 I E − 1/4 = 3 2 m a − 1/4 I E − 1/4 。エネルギーが大きいほど周期が短くなる。調和振動子では面積が E E E の1乗に比例し、周期は E E E によらなかったが、4次のポテンシャルでは面積の増え方が E 3 / 4 E^{3/4} E 3/4 と遅く、その差が周期の E E E 依存性になる。
(3) T = 3 2 × 0.87402 ≈ 3.708 T=3\sqrt2\times0.87402\approx3.708 T = 3 2 × 0.87402 ≈ 3.708 。Hamilton 方程式 q ¨ = − 4 q 3 \ddot q=-4q^3 q ¨ = − 4 q 3 を q ( 0 ) = 1 q(0)=1 q ( 0 ) = 1 , q ˙ ( 0 ) = 0 \dot q(0)=0 q ˙ ( 0 ) = 0 から数値で積分して q ˙ \dot q q ˙ が2度目に 0 0 0 になる時刻の2倍を測っても、3.708 3.708 3.708 になる。エネルギーを16倍にすると周期は 16 − 1 / 4 = 1 2 16^{-1/4}=\frac12 1 6 − 1/4 = 2 1 倍になる。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com