tecktom.com

当てる道具:特性関数を母関数にする

2026年10月15日7分で読めますAI生成教科書『虎の巻』第1部 第14章 作用角変数と断熱不変量
目次
  1. 母関数
  2. 角変数は1周で 増える
  3. 相空間の面積は変数を変えても同じ
  4. 木下の層

AI 生成・木下が未勉強

解析力学から様々な理論分野を組み立てる 第1部 第14章 §1.14.2 ★★★ 前の章の特性関数 W(q)W(q) の定数を EE から JJ に取り替えると、作用角変数への母関数になる

母関数 W(q,J)W(q,J)

前の章(第1部第13章)では、時刻を含まない HH に対して Hamilton–Jacobi 方程式を S=W(q)−EtS=W(q)-Et と分離し、特性関数 WW が H(q,∂W/∂q)=EH(q,\partial W/\partial q)=E を満たすことを見た。自由度1ではこれは ∂W/∂q=p(q,E)\partial W/\partial q=p(q,E) で、p(q,E)p(q,E) は等エネルギー線を pp について解いたものである。積分すれば

W(q,E)=∫q0qp(q′,E) dq′(1.14.7)W(q,E)=\int_{q_0}^{q}p(q',E)\,\mathrm{d} q'\tag{1.14.7}

となる。前の章ではこの WW を、定数 EE を新しい運動量とみなした母関数 F2F_2 として使い、新しい座標 ∂W/∂E=t+β1\partial W/\partial E=t+\beta_1 を得た。

ここでは定数の選び方を変える。EE は JJ の関数 E=H(J)E=H(J) なので、WW を qq と JJ の関数 W(q,J)=W(q,H(J))W(q,J)=W(q,H(J)) とみなし、JJ を新しい運動量にした母関数F2(q,J)=W(q,J)F_2(q,J)=W(q,J) を使う。母関数の規則(第1部第12章)p=∂F2/∂qp=\partial F_2/\partial q, θ=∂F2/∂J\theta=\partial F_2/\partial J から

p=∂W∂q=p(q,H(J)),θ=∂W∂J(1.14.8)p=\frac{\partial W}{\partial q}=p(q,H(J)),\qquad \theta=\frac{\partial W}{\partial J}\tag{1.14.8}

を得る。1つ目は「(q,p)(q,p) がエネルギー H(J)H(J) の等エネルギー線の上にある」ことを言っている。F2F_2 は時刻を含まないので、新しい Hamiltonian はもとの HH をそのまま新しい変数で書いたもので、その値は等エネルギー線の上で H(J)H(J) である。これで定理 1.14.2 の (1)(2) が示せた。θ\theta の運動はθ˙=∂H(J)/∂J\dot\theta=\partial H(J)/\partial J で、JJ が一定なので θ˙\dot\theta も一定である。

角変数は1周で 2π2\pi 増える

定理の (3) を示す。運動が1周するあいだに θ\theta がどれだけ増えるかは、θ=∂W/∂J\theta=\partial W/\partial J を qq について1周積分すればよい。JJ を固定して qq について微分してから積分すると

Δθ=∮∂θ∂q dq=∮∂∂q(∂W∂J)dq=∂∂J∮∂W∂q dq=∂∂J∮p dq=∂∂J(2πJ)=2π(1.14.9)\Delta\theta=\oint\frac{\partial \theta}{\partial q}\,\mathrm{d} q=\oint\frac{\partial }{\partial q}\Bigl(\frac{\partial W}{\partial J}\Bigr)\mathrm{d} q=\frac{\partial }{\partial J}\oint\frac{\partial W}{\partial q}\,\mathrm{d} q=\frac{\partial }{\partial J}\oint p\,\mathrm{d} q=\frac{\partial }{\partial J}(2\pi J)=2\pi\tag{1.14.9}

となる。3つ目の等号では qq についての微分と JJ についての微分の順を入れ替え、JJ についての微分を積分の外へ出した(積分の経路は JJ ごとに変わるが、閉曲線を1周するので端の寄与は無い)。最後の等号は作用変数の定義そのものである。JJ を ∮p dq/(2π)\oint p\,\mathrm{d} q/(2\pi) と定義したことが、そのまま角変数の周期を 2π2\pi にしている。JJ を別の定数倍にとれば、角変数の周期もその逆数倍になる。

1周で θ\theta が 2π2\pi 増え、その速さが一定値 ω=θ˙\omega=\dot\theta なので、1周にかかる時間 TT は ωT=2π\omega T=2\pi を満たす。これと式 (1.14.2) を合わせると

ω(J)=dHdJ=1dJ/dE=2πT(1.14.10)\omega(J)=\frac{\mathrm{d} H}{\mathrm{d} J}=\frac{1}{\mathrm{d} J/\mathrm{d} E}=\frac{2\pi}{T}\tag{1.14.10}

で、作用角変数の Hamilton 方程式から出る角速度と、前の章で面積の EE 微分から出した周期とが同じものであることが確かめられる。作用角変数を使うと、運動を解かずに周期が出る。 J(E)J(E) を面積として計算し、逆に解いて H(J)H(J) を出し、JJ で微分すればよい。

相空間の面積は変数を変えても同じ

正準変換は相空間の面積を保つ(第1部第8章・第1部第12章)。(θ,J)(\theta,J) の平面では、エネルギー EE の等エネルギー線はJ=J(E)J=J(E) の水平な線分 0≤θ<2π0\le\theta<2\pi で、その下の面積は 2πJ2\pi J である。(q,p)(q,p) の平面でこれに対応するのは等エネルギー線の囲む領域で、その面積はやはり 2πJ2\pi J である。作用角変数は、相空間の閉曲線を高さ JJ・幅 2π2\pi の長方形に「伸ばして並べる」変換になっている。調和振動子の式 (1.14.5) では、変換のヤコビアン ∂(q,p)/∂(θ,J)\partial(q,p)/\partial(\theta,J) を計算すると 11 になり、面積が保たれていることが直接確かめられる。

例題 1.14.2 箱の中の粒子の作用角変数 ★★★

オリジナル質量 mm の粒子が、0<q<a0<q<a の範囲を自由に動き、両端の壁で弾性的に跳ね返る(ポテンシャルは箱の中で 00、壁で無限大)。

  1. (1) エネルギー EE の運動の作用変数 J(E)J(E) を求めよ。
  2. (2) H(J)H(J) と ω(J)=dH/dJ\omega(J)=\mathrm{d} H/\mathrm{d} J を求め、粒子の速さ vv を使って ω=πv/a\omega=\pi v/a と書けることを示せ。これが往復の時間と合うことを確かめよ。
  3. (3) 右向きに動いているあいだ(p>0p>0)の角変数 θ\theta を qq で表せ。
  4. (4) m=1m=1, a=2a=2, E=1/2E=1/2 のとき JJ, ω\omega, 周期 TT の値を求めよ。
解答

(1) 箱の中では p=±2mEp=\pm\sqrt{2mE} で一定。軌道は図 1.14.1(b) の長方形で、面積は a×22mEa\times2\sqrt{2mE}。よって

J=12π⋅2a2mE=a2mEπ.J=\frac{1}{2\pi}\cdot2a\sqrt{2mE}=\frac{a\sqrt{2mE}}{\pi}.

(2) 逆に解いて E=H(J)=π2J22ma2E=H(J)=\dfrac{\pi^2J^2}{2ma^2}。ω=dHdJ=π2Jma2\omega=\frac{\mathrm{d} H}{\mathrm{d} J}=\dfrac{\pi^2J}{ma^2}。速さ vv の粒子では 2mE=mv\sqrt{2mE}=mv なので J=amv/πJ=amv/\pi を入れて ω=πv/a\omega=\pi v/a。周期は T=2π/ω=2a/vT=2\pi/\omega=2a/v で、距離 aa を往復する時間と一致する。調和振動子と違い、ω\omega が JJ(エネルギー)によって変わる。

(3) p>0p>0 の枝で W=∫0q2mE dq′=q2mE=πJq/aW=\int_0^q\sqrt{2mE}\,\mathrm{d} q'=q\sqrt{2mE}=\pi Jq/a。θ=∂W∂J=πq/a\theta=\frac{\partial W}{\partial J}=\pi q/a で、qq が 00 から aa まで進むあいだに θ\theta は 00 から π\pi まで一様に増える。帰り道(p<0p<0)では WW に −2mE (q−a)-\sqrt{2mE}\,(q-a) が足され、θ=π+π(a−q)/a\theta=\pi+\pi(a-q)/a が π\pi から 2π2\pi まで増える。1周で 2π2\pi で、粒子の位置は時間に比例して進むので θ\theta も時間に比例する。

(4) J=22⋅12/π=2/π≈0.637J=2\sqrt{2\cdot\frac12}/\pi=2/\pi\approx0.637、ω=π2(2/π)/(1⋅4)=π/2≈1.57\omega=\pi^2(2/\pi)/(1\cdot4)=\pi/2\approx1.57、T=2π/ω=4T=2\pi/\omega=4。速さは v=2E/m=1v=\sqrt{2E/m}=1 で、往復の距離 44 を速さ 11 で進む時間 44 と一致する。

木下の勉強ステータス

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

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

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