解析力学から様々な理論分野を組み立てる 第1部 第13章 §1.13.3 ★★★ H が時刻を含まなければ S=W(q)−Et。1次元では W(q)=∫pdq で、運動は積分1つで求まる
時間の分離
Hamiltonian が時刻を陽に含まないとき、Hamilton–Jacobi 方程式 (1.13.2) で t が現れるのは ∂S/∂t だけである。そこで
S(q,α,t)=W(q,α)−Et,H(q,∂q∂W)=E(1.13.8)
と置くと、t を含まない方程式になる。後者を時間に依らない Hamilton–Jacobi 方程式、W(q) を特性関数と呼ぶ。E はエネルギーの値で、完全解の定数の1つ(たとえば α1=E)として使う。特性関数は引数 q を必ず付けて書き、仕事 W や状態数 W(E) と区別する。
Jacobi の定理を α1=E に当てると β1=∂S/∂E=∂W/∂E−t なので
∂E∂W=t+β1(1.13.9)
となる。W そのものも母関数とみなせて、F2=W(q,P) で P1=E とすると、新しい Hamiltonian は K=H=E=P1 で、新しい座標は Q˙1=∂K/∂P1=1 より Q1=t+β1 と時刻に比例して進む。S は新しい変数を全部止める変換、W は時刻に比例して進む座標を1つ残す変換で、式 (1.13.9) はどちらの見方でも同じ式になる。
1次元の運動
自由度1、H=p2/(2m)+V(q) なら、式 (1.13.8) は 2m1(dW/dq)2+V(q)=E で、解は
W(q,E)=∫qp(q′,E)dq′,p(q,E)=±2m(E−V(q))(1.13.10)
である。符号は運動の向きで、転回点(E=V(q))を通るたびに入れ替わるので、W は経路に沿って積分を続けた値として扱う。式 (1.13.9) は
t+β1=∂E∂W=∫q∂E∂pdq′=∫qp(q′,E)mdq′(1.13.11)
となり、時刻が位置の積分として求まる。これを q について解けば運動 q(t) である。
調和振動子 V=21mω2q2 では p=2mE−m2ω2q2(p>0 の枝)で、
∂E∂W=∫0q2mE−m2ω2q′2mdq′=ω1arcsin(q2Emω2)
なので、式 (1.13.9) から
q(t)=mω22Esinω(t+β1)(1.13.12)
を得る。定数は振幅を決める E と位相を決める β1 で、運動方程式を一度も積分せずに一般解が出た。
W(q,E)=∫pdq は、相空間で等エネルギー線 H(q,p)=E の下の面積である。調和振動子の等エネルギー線は半軸 2E/mω2 と 2mE の楕円なので、1周の積分は楕円の面積
∮pdq=πmω22E2mE=ω2πE(1.13.13)
になる。この量を 2π で割った J=E/ω を新しい運動量に選ぶのが、次の章(第1部第14章)の作用角変数である。
変数が分離できる場合
自由度が2以上でも、Hamiltonian の形によっては W を各変数の関数の和 W=∑iWi(qi,α) に分けられる。最も簡単なのは循環座標で、qk が H に現れなければ pk は保存し、W は αkqk という項を持つ(∂W/∂qk=αk が保存する運動量の値)。一般に、方程式 (1.13.8) を変形して「q1 と ∂W/∂q1 だけを含む式」=「残りの変数だけを含む式」の形にできれば、両辺は定数でなければならず、W1 についての常微分方程式が1本切り出せる。この操作を繰り返して全部の変数が分かれるとき、系は変数分離できるという。どの座標系で分離できるかは Hamiltonian の対称性で決まり、中心力なら極座標で分離できる(次の節)。
例題 1.13.3 初期条件から定数を決める ★★
オリジナルm=1, ω=2 の調和振動子が、エネルギー E=8 を持ち、時刻 0 に q=1 を q が増える向きに通過した。
- (1) 式 (1.13.12) の振幅と β1 を求めよ。
- (2) 時刻 0 の運動量 p(0) を、Jacobi の定理の p=∂S/∂q から求めよ。
- (3) 相空間で1周する間の ∮pdq を求め、エネルギーと角振動数の比と比べよ。
解答
(1) 振幅は 2E/(mω2)=16/4=2。q(0)=2sin2β1=1 から sin2β1=21 で、q が増える向き(q˙(0)∝cos2β1>0)なので 2β1=π/6、β1=π/12。
(2) S=W−Et の q 微分は ∂W/∂q=p=+2mE−m2ω2q2(p>0 の枝)で、q=1 では p(0)=16−4=23。式 (1.13.12) から直接 mq˙(0)=2⋅2cos(π/6)=23 となり一致する。エネルギーは 21(23)2+21⋅4⋅12=6+2=8 で、与えられた値に戻る。
(3) 式 (1.13.13) から ∮pdq=2πE/ω=8π。E/ω=4 の 2π 倍である。エネルギーを変えずに β1 だけ変えると、同じ楕円の上で出発点がずれるだけなので、∮pdq は E だけで決まる。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com