tecktom.com

この章のまとめ

2026年10月15日11分で読めますAI生成教科書『虎の巻』第1部 第13章 Hamilton–Jacobi 方程式

AI 生成・木下が未勉強

解析力学から様々な理論分野を組み立てる 第1部 第13章 §1.13.7 ★★★ Hamilton–Jacobi 方程式は、母関数の条件・終点の関数としての作用・波の位相の古典極限の3通りに読める

新しい Hamiltonian が恒等的に 00 になる正準変換を母関数 F2=S(q,α,t)F_2=S(q,\alpha,t) で作ろうとすると、SS は Hamilton–Jacobi 方程式 H(q,∂S/∂q,t)+∂S/∂t=0H(q,\partial S/\partial q,t)+\partial S/\partial t=0 を満たさなければならない。ff 個の独立な定数を含む完全解が見つかれば、Jacobi の定理により ∂S/∂αi=βi\partial S/\partial\alpha_i=\beta_i と pi=∂S/∂qip_i=\partial S/\partial q_i から運動が微分と代数だけで求まる。完全解の1つは、始点と終点を結ぶ真の経路の作用を終点の関数とみなしたもので、∂S/∂q=p\partial S/\partial q=p, ∂S/∂t=−H\partial S/\partial t=-H を満たし、時間発展という正準変換の母関数になる。始点を出た近くの経路が再び集まる共役点では、この SS は定義できない。HH が時刻を含まなければ S=W(q)−EtS=W(q)-Et と分離でき、1次元では W=∫p dqW=\int p\,\mathrm{d} q、時刻は ∂W/∂E\partial W/\partial E、周期は ∮p dq\oint p\,\mathrm{d} q の EE 微分になる。中心力では極座標で WW が分かれ、分離の定数として角運動量が現れ、Kepler 問題の軌道は積分1つで円錐曲線になる。SS の等値面は波面で、経路はそれに直交し、Hamilton–Jacobi 方程式の特性曲線が運動の経路である。Schrödinger 方程式に ψ=AeiS/ℏ\psi=A\mathrm{e}^{\mathrm{i} S/\hbar} を入れると、ℏ2\hbar^2 の項を除いて Hamilton–Jacobi 方程式と連続の式が出る。次の章(第1部第14章)では、周期運動の ∮p dq\oint p\,\mathrm{d} q を新しい運動量に選んで作用角変数を作り、ゆっくり変わる系でそれが保たれることを示す。

この章の道具

  • Hamilton–Jacobi 方程式:H(q,∂S/∂q,t)+∂S/∂t=0H(q,\partial S/\partial q,t)+\partial S/\partial t=0。F2=SF_2=S の正準変換で K≡0K\equiv0 となる条件。
  • Jacobi の定理:完全解 S(q,α,t)S(q,\alpha,t) から ∂S/∂αi=βi\partial S/\partial\alpha_i=\beta_i, pi=∂S/∂qip_i=\partial S/\partial q_i で運動が決まる。
  • 真の経路の作用 S(q,t;q0,t0)S(q,t;q_0,t_0):dS=∑p dq−H dt−∑p0 dq0+H0 dt0\mathrm{d} S=\sum p\,\mathrm{d} q-H\,\mathrm{d} t-\sum p_0\,\mathrm{d} q_0+H_0\,\mathrm{d} t_0。時間発展の母関数。
  • 時間の分離:S=W(q)−EtS=W(q)-Et、H(q,∂W/∂q)=EH(q,\partial W/\partial q)=E、∂W/∂E=t+β1\partial W/\partial E=t+\beta_1、周期 T=ddE∮p dqT=\frac{\mathrm{d} }{\mathrm{d} E}\oint p\,\mathrm{d} q。
  • 古典極限:ψ=AeiS/ℏ\psi=A\mathrm{e}^{\mathrm{i} S/\hbar} で ℏ→0\hbar\to0 とすると Hamilton–Jacobi 方程式と連続の式。条件は ∣d(ℏ/p)/dq∣≪1\lvert \mathrm{d}(\hbar/p)/\mathrm{d} q\rvert \ll1。

総合問題 一様な重力の中の粒子を Hamilton–Jacobi 方程式で6通りに見る ★★★

オリジナル質量 mm の粒子が、鉛直上向きの座標 qq の上で一様な重力 V(q)=mgqV(q)=mgq を受けて運動する(H=p2/(2m)+mgqH=p^2/(2m)+mgq)。時刻 00 に q0q_0 を出て時刻 tt に qq に着く運動を考える。

  1. (1) (使う道具:第1部第4章の Euler–Lagrange 方程式と1.13.2節)端点を結ぶ真の経路を求め、その作用がS(q,t;q0)=m(q−q0)22t−mgt(q+q0)2−mg2t324S(q,t;q_0)=\dfrac{m(q-q_0)^2}{2t}-\dfrac{mgt(q+q_0)}{2}-\dfrac{mg^2t^3}{24} になることを示せ。
  2. (2) (使う道具:1.13.1節・1.13.2節)(1) の SS が Hamilton–Jacobi 方程式を満たし、∂S/∂q=p\partial S/\partial q=p, ∂S/∂q0=−p0\partial S/\partial q_0=-p_0 となることを確かめよ。
  3. (3) (使う道具:1.13.3節)時間に依らない Hamilton–Jacobi 方程式の解 W(q,E)W(q,E) を、∂W/∂q=+2m(E−mgq)\partial W/\partial q=+\sqrt{2m(E-mgq)} の枝で求めよ。式 (1.13.9) から q(t)q(t) を求めよ。
  4. (4) (使う道具:第1部第11章の Poisson 括弧・第1部第12章の母関数)(3) の新しい変数 Q=∂W/∂E−tQ=\partial W/\partial E-t, P=EP=E を (q,p)(q,p) で表し、{Q,P}=1\{Q,P\}=1 を確かめよ。F2=W(q,P)−PtF_2=W(q,P)-Pt の新しい Hamiltonian が K=0K=0 であることも確かめよ。
  5. (5) (使う道具:1.13.5節と近似の数理 第7章)古典極限の条件 (1.13.22) の左辺を速さ vv で表せ。中性子(m=1.675×10−27 kgm=1.675\times10^{-27}\,\mathrm{kg})が地上の重力(g=9.81 m s−2g=9.81\,\mathrm{m\,s^{-2}})の中を v=1 m s−1v=1\,\mathrm{m\,s^{-1}} で動くときの値と、左辺が 11 になる速さ v∗v_*、最高点から速さ v∗v_* になるまでに落ちる高さ v∗2/(2g)v_*^2/(2g) を求めよ(ℏ=1.055×10−34 J s\hbar=1.055\times10^{-34}\,\mathrm{J\,s})。
  6. (6) (使う道具:第1部第7章の Legendre 変換)上向きに動き続ける運動(p>0p>0)について、Φ(E)=W(q,E)−W(q0,E)−Et\Phi(E)=W(q,E)-W(q_0,E)-Et を EE で停留させた値が S(q,t;q0)S(q,t;q_0) に等しいことを、m=g=1m=g=1, q0=0q_0=0, q=38q=\frac38, t=12t=\frac12 の場合に確かめよ。
解答

(1) Euler–Lagrange 方程式は mq¨=−mgm\ddot q=-mg。q(0)=q0q(0)=q_0, q(t)=qq(t)=q を満たす解は

q(τ)=q0+(q−q0)τt+g2τ(t−τ)(0≤τ≤t)q(\tau)=q_0+(q-q_0)\frac{\tau}{t}+\frac{g}{2}\tau(t-\tau)\qquad(0\le\tau\le t)

である(q¨=−g\ddot q=-g で、τ=0,t\tau=0,t で第3項が消える)。q˙(τ)=(q−q0)/t+g2(t−2τ)\dot q(\tau)=(q-q_0)/t+\frac g2(t-2\tau) を L=12mq˙2−mgqL=\frac12m\dot q^2-mgq に入れて τ\tau について 00 から tt まで積分すると、∫0t12mq˙2 dτ=m(q−q0)22t+mg2t324\int_0^t\frac12m\dot q^2\,\mathrm{d}\tau=\frac{m(q-q_0)^2}{2t}+\frac{mg^2t^3}{24}(交差項は ∫0t(t−2τ)dτ=0\int_0^t(t-2\tau)\mathrm{d}\tau=0 で消える)、∫0tmgq dτ=mg[(q+q0)t2+gt312]\int_0^tmgq\,\mathrm{d}\tau=mg\bigl[\frac{(q+q_0)t}{2}+\frac{gt^3}{12}\bigr] なので

S=m(q−q0)22t−mgt(q+q0)2−mg2t324.S=\frac{m(q-q_0)^2}{2t}-\frac{mgt(q+q_0)}{2}-\frac{mg^2t^3}{24}.

q0=q=0q_0=q=0, t=Tt=T(高さ 00 から投げ上げて高さ 00 に戻る運動)とおくと S=−mg2T3/24S=-mg^2T^3/24 で、第1部第1章の例題 1.1.1 の S(0)S(0) と一致する。

(2) ∂S/∂q=m(q−q0)/t−mgt/2\partial S/\partial q=m(q-q_0)/t-mgt/2 は経路の終点の運動量 mq˙(t)=m(q−q0)/t−mgt/2m\dot q(t)=m(q-q_0)/t-mgt/2 に等しい。∂S/∂q0=−m(q−q0)/t−mgt/2=−mq˙(0)=−p0\partial S/\partial q_0=-m(q-q_0)/t-mgt/2=-m\dot q(0)=-p_0 である。∂S/∂t=−m(q−q0)22t2−mg(q+q0)2−mg2t28\partial S/\partial t=-\frac{m(q-q_0)^2}{2t^2}-\frac{mg(q+q_0)}{2}-\frac{mg^2t^2}{8} と、12m(∂S/∂q)2=m(q−q0)22t2−mg(q−q0)2+mg2t28\frac{1}{2m}(\partial S/\partial q)^2=\frac{m(q-q_0)^2}{2t^2}-\frac{mg(q-q_0)}{2}+\frac{mg^2t^2}{8} を足すと−mgq-mgq が残り、mgqmgq を足して 00 になる。Hamilton–Jacobi 方程式が成り立つ。

(3) 12m(dW/dq)2+mgq=E\frac{1}{2m}(\mathrm{d} W/\mathrm{d} q)^2+mgq=E から dW/dq=2m(E−mgq)\mathrm{d} W/\mathrm{d} q=\sqrt{2m(E-mgq)} で、積分して

W(q,E)=−22m3mg (E−mgq)3/2W(q,E)=-\frac{2\sqrt{2m}}{3mg}\,(E-mgq)^{3/2}

(qq で微分すると 2m (E−mgq)1/2\sqrt{2m}\,(E-mgq)^{1/2} に戻る)。∂W/∂E=−2mmg(E−mgq)1/2=−pmg\partial W/\partial E=-\frac{\sqrt{2m}}{mg}(E-mgq)^{1/2}=-\frac{p}{mg} なので、式 (1.13.9) は −p/(mg)=t+β1-p/(mg)=t+\beta_1、すなわち p=−mg(t+β1)p=-mg(t+\beta_1) である(p>0p>0 の枝なので t<−β1t<-\beta_1)。E=p2/(2m)+mgqE=p^2/(2m)+mgq を qq について解いて

q(t)=Emg−g2(t+β1)2.q(t)=\frac{E}{mg}-\frac g2(t+\beta_1)^2 .

最高点の高さ E/(mg)E/(mg) に時刻 −β1-\beta_1 に達する放物運動である。

(4) Q=−p/(mg)−tQ=-p/(mg)-t, P=p2/(2m)+mgqP=p^2/(2m)+mgq なので

{Q,P}=∂Q∂q∂P∂p−∂Q∂p∂P∂q=0⋅pm−(−1mg)mg=1.\{Q,P\}=\frac{\partial Q}{\partial q}\frac{\partial P}{\partial p}-\frac{\partial Q}{\partial p}\frac{\partial P}{\partial q}=0\cdot\frac pm-\Bigl(-\frac{1}{mg}\Bigr)mg=1 .

F2=W(q,P)−PtF_2=W(q,P)-Pt では p=∂F2/∂q=∂W/∂qp=\partial F_2/\partial q=\partial W/\partial q を HH に入れると H=PH=P となり、K=H+∂F2/∂t=P−P=0K=H+\partial F_2/\partial t=P-P=0 である。QQ と PP はどちらも運動の定数で、−Q=β1-Q=\beta_1 が最高点に達する時刻にあたる。

(5) p=mvp=mv で dp/dq=−m2g/p\mathrm{d} p/\mathrm{d} q=-m^2g/p((3) の p2=2m(E−mgq)p^2=2m(E-mgq) を微分した)なので

∣ddqℏp∣=ℏp2∣dpdq∣=ℏm2gp3=ℏgmv3.\Bigl|\frac{\mathrm{d} }{\mathrm{d} q}\frac{\hbar}{p}\Bigr|=\frac{\hbar}{p^2}\Bigl|\frac{\mathrm{d} p}{\mathrm{d} q}\Bigr|=\frac{\hbar m^2g}{p^3}=\frac{\hbar g}{mv^3}.

中性子で v=1 m s−1v=1\,\mathrm{m\,s^{-1}} なら 1.055×10−34×9.811.675×10−27×1≈6.2×10−7\frac{1.055\times10^{-34}\times9.81}{1.675\times10^{-27}\times1}\approx6.2\times10^{-7} で、条件は十分に満たされる。左辺が 11 になる速さは v∗=(ℏg/m)1/3≈8.5×10−3 m s−1v_*=(\hbar g/m)^{1/3}\approx8.5\times10^{-3}\,\mathrm{m\,s^{-1}} で、最高点からの高さに直すと v∗2/(2g)≈3.7×10−6 mv_*^2/(2g)\approx3.7\times10^{-6}\,\mathrm m である。最高点(転回点)から数 μm\mu\mathrm m の範囲では Hamilton–Jacobi 方程式による記述が使えず、そこでは波動関数を Schrödinger 方程式で直接扱う必要がある(第4部第7章)。

(6) Φ′(E)=∂W/∂E(q)−∂W/∂E(q0)−t=(p0−p)/(mg)−t\Phi'(E)=\partial W/\partial E(q)-\partial W/\partial E(q_0)-t=\bigl(p_0-p\bigr)/(mg)-t なので、停留の条件は p0−p=mgtp_0-p=mgt で、これは (1) の経路で運動量が時間 tt の間に mgtmgt 減ることと同じである。数値では m=g=1m=g=1 で p=2(E−38)p=\sqrt{2(E-\frac38)}, p0=2Ep_0=\sqrt{2E}、2E−2E−34=12\sqrt{2E}-\sqrt{2E-\frac34}=\frac12 を解いて E=12E=\frac12(p0=1p_0=1, p=12p=\frac12)。そのとき

Φ=W(38,12)−W(0,12)−12⋅12=−124+13−14=124\Phi=W\Bigl(\tfrac38,\tfrac12\Bigr)-W\Bigl(0,\tfrac12\Bigr)-\tfrac12\cdot\tfrac12=-\frac{1}{24}+\frac13-\frac14=\frac{1}{24}

(W(q,E)=−223(E−q)3/2W(q,E)=-\frac{2\sqrt2}{3}(E-q)^{3/2} で、(1/8)3/2=1/(162)(1/8)^{3/2}=1/(16\sqrt2), (1/2)3/2=1/(22)(1/2)^{3/2}=1/(2\sqrt2) を使った)。(1) の式からも S=(3/8)21−(1/2)(3/8)2−(1/2)324=964−332−1192=124S=\frac{(3/8)^2}{1}-\frac{(1/2)(3/8)}{2}-\frac{(1/2)^3}{24}=\frac{9}{64}-\frac{3}{32}-\frac{1}{192}=\frac{1}{24} で、一致する。S=W−EtS=W-Et のうち EE を消して時刻 tt を変数にする操作は、−∂Φ/∂E=0-\partial\Phi/\partial E=0 の停留点で値をとる Legendre 変換で、EE と tt が互いに共役な変数の組になっている。

木下の勉強ステータス

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

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

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