tecktom.com

出てくるもの:時間を分離した特性関数と1次元の運動

2026年10月15日6分で読めますAI生成教科書『虎の巻』第1部 第13章 Hamilton–Jacobi 方程式
目次
  1. 時間の分離
  2. 1次元の運動
  3. 変数が分離できる場合
  4. 木下の層

AI 生成・木下が未勉強

解析力学から様々な理論分野を組み立てる 第1部 第13章 §1.13.3 ★★★ HH が時刻を含まなければ S=W(q)−EtS=W(q)-Et。1次元では W(q)=∫p dqW(q)=\int p\,\mathrm{d} q で、運動は積分1つで求まる

時間の分離

Hamiltonian が時刻を陽に含まないとき、Hamilton–Jacobi 方程式 (1.13.2) で tt が現れるのは ∂S/∂t\partial S/\partial t だけである。そこで

S(q,α,t)=W(q,α)−Et,H(q,∂W∂q)=E(1.13.8)S(q,\alpha,t)=W(q,\alpha)-Et,\qquad H\Bigl(q,\frac{\partial W}{\partial q}\Bigr)=E\tag{1.13.8}

と置くと、tt を含まない方程式になる。後者を時間に依らない Hamilton–Jacobi 方程式、W(q)W(q) を特性関数と呼ぶ。EE はエネルギーの値で、完全解の定数の1つ(たとえば α1=E\alpha_1=E)として使う。特性関数は引数 qq を必ず付けて書き、仕事 WW や状態数 W(E)W(E) と区別する。

Jacobi の定理を α1=E\alpha_1=E に当てると β1=∂S/∂E=∂W/∂E−t\beta_1=\partial S/\partial E=\partial W/\partial E-t なので

∂W∂E=t+β1(1.13.9)\frac{\partial W}{\partial E}=t+\beta_1\tag{1.13.9}

となる。WW そのものも母関数とみなせて、F2=W(q,P)F_2=W(q,P) で P1=EP_1=E とすると、新しい Hamiltonian は K=H=E=P1K=H=E=P_1 で、新しい座標は Q˙1=∂K/∂P1=1\dot Q_1=\partial K/\partial P_1=1 より Q1=t+β1Q_1=t+\beta_1 と時刻に比例して進む。SS は新しい変数を全部止める変換、WW は時刻に比例して進む座標を1つ残す変換で、式 (1.13.9) はどちらの見方でも同じ式になる。

1次元の運動

自由度1、H=p2/(2m)+V(q)H=p^2/(2m)+V(q) なら、式 (1.13.8) は 12m(dW/dq)2+V(q)=E\frac{1}{2m}(\mathrm{d} W/\mathrm{d} q)^2+V(q)=E で、解は

W(q,E)=∫qp(q′,E) dq′,p(q,E)=±2m(E−V(q))(1.13.10)W(q,E)=\int^q p(q',E)\,\mathrm{d} q',\qquad p(q,E)=\pm\sqrt{2m\bigl(E-V(q)\bigr)}\tag{1.13.10}

である。符号は運動の向きで、転回点(E=V(q)E=V(q))を通るたびに入れ替わるので、WW は経路に沿って積分を続けた値として扱う。式 (1.13.9) は

t+β1=∂W∂E=∫q∂p∂E dq′=∫qmp(q′,E) dq′(1.13.11)t+\beta_1=\frac{\partial W}{\partial E}=\int^q\frac{\partial p}{\partial E}\,\mathrm{d} q'=\int^q\frac{m}{p(q',E)}\,\mathrm{d} q'\tag{1.13.11}

となり、時刻が位置の積分として求まる。これを qq について解けば運動 q(t)q(t) である。

調和振動子 V=12mω2q2V=\frac12m\omega^2q^2 では p=2mE−m2ω2q2p=\sqrt{2mE-m^2\omega^2q^2}(p>0p>0 の枝)で、

∂W∂E=∫0qm dq′2mE−m2ω2q′2=1ωarcsin⁡(qmω22E)\frac{\partial W}{\partial E}=\int_0^q\frac{m\,\mathrm{d} q'}{\sqrt{2mE-m^2\omega^2q'^2}}=\frac1\omega\arcsin\Bigl(q\sqrt{\frac{m\omega^2}{2E}}\Bigr)

なので、式 (1.13.9) から

q(t)=2Emω2 sin⁡ω(t+β1)(1.13.12)q(t)=\sqrt{\frac{2E}{m\omega^2}}\,\sin\omega(t+\beta_1)\tag{1.13.12}

を得る。定数は振幅を決める EE と位相を決める β1\beta_1 で、運動方程式を一度も積分せずに一般解が出た。

W(q,E)=∫p dqW(q,E)=\int p\,\mathrm{d} q は、相空間で等エネルギー線 H(q,p)=EH(q,p)=E の下の面積である。調和振動子の等エネルギー線は半軸 2E/mω2\sqrt{2E/m\omega^2} と 2mE\sqrt{2mE} の楕円なので、1周の積分は楕円の面積

∮p dq=π2Emω22mE=2πEω(1.13.13)\oint p\,\mathrm{d} q=\pi\sqrt{\frac{2E}{m\omega^2}}\sqrt{2mE}=\frac{2\pi E}{\omega}\tag{1.13.13}

になる。この量を 2π2\pi で割った J=E/ωJ=E/\omega を新しい運動量に選ぶのが、次の章(第1部第14章)の作用角変数である。

変数が分離できる場合

自由度が2以上でも、Hamiltonian の形によっては WW を各変数の関数の和 W=∑iWi(qi,α)W=\sum_iW_i(q_i,\alpha) に分けられる。最も簡単なのは循環座標で、qkq_k が HH に現れなければ pkp_k は保存し、WW は αkqk\alpha_kq_k という項を持つ(∂W/∂qk=αk\partial W/\partial q_k=\alpha_k が保存する運動量の値)。一般に、方程式 (1.13.8) を変形して「q1q_1 と ∂W/∂q1\partial W/\partial q_1 だけを含む式」=「残りの変数だけを含む式」の形にできれば、両辺は定数でなければならず、W1W_1 についての常微分方程式が1本切り出せる。この操作を繰り返して全部の変数が分かれるとき、系は変数分離できるという。どの座標系で分離できるかは Hamiltonian の対称性で決まり、中心力なら極座標で分離できる(次の節)。

例題 1.13.3 初期条件から定数を決める ★★

オリジナルm=1m=1, ω=2\omega=2 の調和振動子が、エネルギー E=8E=8 を持ち、時刻 00 に q=1q=1 を qq が増える向きに通過した。

  1. (1) 式 (1.13.12) の振幅と β1\beta_1 を求めよ。
  2. (2) 時刻 00 の運動量 p(0)p(0) を、Jacobi の定理の p=∂S/∂qp=\partial S/\partial q から求めよ。
  3. (3) 相空間で1周する間の ∮p dq\oint p\,\mathrm{d} q を求め、エネルギーと角振動数の比と比べよ。
解答

(1) 振幅は 2E/(mω2)=16/4=2\sqrt{2E/(m\omega^2)}=\sqrt{16/4}=2。q(0)=2sin⁡2β1=1q(0)=2\sin2\beta_1=1 から sin⁡2β1=12\sin2\beta_1=\frac12 で、qq が増える向き(q˙(0)∝cos⁡2β1>0\dot q(0)\propto\cos2\beta_1>0)なので 2β1=π/62\beta_1=\pi/6、β1=π/12\beta_1=\pi/12。

(2) S=W−EtS=W-Et の qq 微分は ∂W/∂q=p=+2mE−m2ω2q2\partial W/\partial q=p=+\sqrt{2mE-m^2\omega^2q^2}(p>0p>0 の枝)で、q=1q=1 では p(0)=16−4=23p(0)=\sqrt{16-4}=2\sqrt3。式 (1.13.12) から直接 mq˙(0)=2⋅2cos⁡(π/6)=23m\dot q(0)=2\cdot2\cos(\pi/6)=2\sqrt3 となり一致する。エネルギーは 12(23)2+12⋅4⋅12=6+2=8\frac12(2\sqrt3)^2+\frac12\cdot4\cdot1^2=6+2=8 で、与えられた値に戻る。

(3) 式 (1.13.13) から ∮p dq=2πE/ω=8π\oint p\,\mathrm{d} q=2\pi E/\omega=8\pi。E/ω=4E/\omega=4 の 2π2\pi 倍である。エネルギーを変えずに β1\beta_1 だけ変えると、同じ楕円の上で出発点がずれるだけなので、∮p dq\oint p\,\mathrm{d} q は EE だけで決まる。

木下の勉強ステータス

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

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

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