tecktom.com

出てくるもの (2):HH の保存と循環座標

2026年10月15日10分で読めますAI生成教科書『虎の巻』第1部 第8章 Hamilton 形式と相空間
目次
  1. の時間変化
  2. が力学的エネルギーに等しい条件
  3. 循環座標
  4. 木下の層

AI 生成・木下が未勉強

解析力学から様々な理論分野を組み立てる 第1部 第8章 §1.8.4 ★★★ HH が tt を含まなければ HH が保存し、qkq_k を含まなければ pkp_k が保存して自由度が1つ減る

HH の時間変化

運動に沿った HH の変化を連鎖律で計算し、Hamilton 方程式 (1.8.6) を入れると

dHdt=∑i(∂H∂qiq˙i+∂H∂pip˙i)+∂H∂t=∑i(∂H∂qi∂H∂pi−∂H∂pi∂H∂qi)+∂H∂t=∂H∂t(1.8.14)\frac{\mathrm{d} H}{\mathrm{d} t}=\sum_i\Bigl(\frac{\partial H}{\partial q_i}\dot q_i+\frac{\partial H}{\partial p_i}\dot p_i\Bigr)+\frac{\partial H}{\partial t} =\sum_i\Bigl(\frac{\partial H}{\partial q_i}\frac{\partial H}{\partial p_i}-\frac{\partial H}{\partial p_i}\frac{\partial H}{\partial q_i}\Bigr)+\frac{\partial H}{\partial t}=\frac{\partial H}{\partial t}\tag{1.8.14}

を得る。和の中の2項は打ち消し合う。したがって HH が時刻 tt をあらわに含まなければ、HH の値は運動の間一定である。式 (1.8.5) の第3式から ∂H/∂t=−∂L/∂t\partial H/\partial t=-\partial L/\partial t なので、これは「LL が tt をあらわに含まなければ HH が保存する」と言っても同じである。時刻をずらす変換と HH の保存のつながりは、第1部第10章の Noether の定理で扱う(第1部第10章)。

HH が力学的エネルギーに等しい条件

保存する HH が力学的エネルギー T+VT+V に等しいとは限らない。等しくなる条件を調べる。運動エネルギーが速度の2次の同次式 T=12∑i,jaij(q)q˙iq˙jT=\frac12\sum_{i,j}a_{ij}(q)\dot q_i\dot q_j(aij=ajia_{ij}=a_{ji})で、VV が速度を含まないとする。このとき pi=∑jaijq˙jp_i=\sum_ja_{ij}\dot q_j なので ∑ipiq˙i=∑i,jaijq˙iq˙j=2T\sum_ip_i\dot q_i=\sum_{i,j}a_{ij}\dot q_i\dot q_j=2T となり、

H=∑ipiq˙i−L=2T−(T−V)=T+V(1.8.15)H=\sum_ip_i\dot q_i-L=2T-(T-V)=T+V\tag{1.8.15}

が成り立つ。直交座標から一般化座標への変換が時刻を含まない(拘束が時間によらない)なら、TT はこの形になる。

変換が時刻を含むと、TT に速度の1次の項と0次の項が入り、∑pq˙=2T\sum p\dot q=2T が成り立たなくなる。たとえば一定の角速度で回される枠の上で測った座標を使うと、HH は T+VT+V と違う関数になる。その HH が時刻を含まなければ HH は保存するが、T+VT+V は保存しない(枠を回し続ける外部の装置が系に仕事をするため)。

回る枠の上で時刻によらないポテンシャルの中を動く粒子(天体力学の制限三体問題で、一定の角速度で回る2天体の重力を受ける小天体)がその例である。枠の角速度を Ω\Omega とする。慣性系で測ったエネルギーと角運動量を EE と LzL_z と書くと、Hamiltonian は

H=E−ΩLzH=E-\Omega L_z

になる。EE も LzL_z も別々には保存しないが、この組み合わせは保存し、Jacobi 積分と呼ばれる。

章末の総合問題で、回転する輪にはめたビーズを使ってこれを確かめる。逆に、H=T+VH=T+V であっても、ポテンシャルが時刻を含めば HH は保存しない。HH が保存するかどうかと HH がエネルギーに等しいかどうかは、別の条件で決まる。

加速する台から測った座標では、回る枠と逆向きの食い違いが起こる。直線上の自由な粒子を、一定の加速度 aa で動く台から測った座標 qq で書く(第1部第2章の例と同じ座標)。このとき

x=q+12at2,L=m2(q˙+at)2,p=m(q˙+at)=mx˙,H=pq˙−L=p22m−at p\begin{gathered} x=q+\frac12at^2,\qquad L=\frac m2(\dot q+at)^2,\\ p=m(\dot q+at)=m\dot x,\qquad H=p\dot q-L=\frac{p^2}{2m}-at\,p \end{gathered}

で、Hamiltonian は時刻を陽に含む。

この Hamiltonian は座標を含まないので、運動量 pp は保存する。一方で式 (1.8.14) から dHdt=∂H/∂t=−ap\frac{\mathrm{d} H}{\mathrm{d} t}=\partial H/\partial t=-ap で、HH は時刻に比例して変わる。慣性系で測ったエネルギー 12mx˙2=p2/(2m)\frac12m\dot x^2=p^2/(2m) は一定のままなので、HH の変化は台の座標を選んだことから来ている。回る枠では Hamiltonian が保存して T+VT+V が保存せず、加速する台ではエネルギーが保存して Hamiltonian が保存しないので、2つの例で保存する量が入れ替わる。

循環座標

HH がある座標 qkq_k を含まないとき、qkq_k を循環座標と呼ぶ。Hamilton 方程式から p˙k=−∂H/∂qk=0\dot p_k=-\partial H/\partial q_k=0 なので、pkp_k は保存する。保存する pkp_k の値を定数 ckc_k として HH に入れると、HH は残りの f−1f-1 個の座標と運動量、そして定数 ckc_k の関数になる。残りの 2(f−1)2(f-1) 本の方程式は qkq_k を含まないので、それだけで解くことができ、qkq_k は最後に q˙k=∂H/∂ck\dot q_k=\partial H/\partial c_k を積分して求める。Lagrangian の形では、循環座標 qkq_k が LL に入らなくても速度 q˙k\dot q_k は LL に残るので、このように1変数を定数に置き換えることはできない。Hamilton の形では、循環座標に共役な運動量がそのまま定数になり、問題の自由度が1つ減る。

中心力の問題がその例である。例題 1.8.1 の H=pr22m+pθ22mr2+V(r)H=\frac{p_r^2}{2m}+\frac{p_\theta^2}{2mr^2}+V(r) は θ\theta を含まないので pθp_\theta が保存し、動径方向の運動は自由度1の Hamiltonian

Hr(r,pr)=pr22m+Ueff(r),Ueff(r)=pθ22mr2+V(r)(1.8.16)H_r(r,p_r)=\frac{p_r^2}{2m}+U_{\mathrm{eff}}(r),\qquad U_{\mathrm{eff}}(r)=\frac{p_\theta^2}{2mr^2}+V(r)\tag{1.8.16}

で決まる。UeffU_{\mathrm{eff}} を有効ポテンシャルと呼ぶ。第1項は、pθp_\theta を一定に保ったまま rr を変えると回転の運動エネルギーが変わることを表し、例題 1.8.1 (3) の遠心力はこの項の −∂/∂r-\partial/\partial r である。動径の運動が求まれば、角度は θ˙=∂H/∂pθ=pθ/(mr2)\dot\theta=\partial H/\partial p_\theta=p_\theta/(mr^2) を積分して求める。

図

図 1.8.2 Kepler 問題の有効ポテンシャル(r0=pθ2/(mk)r_0=p_\theta^2/(mk) を長さの単位とする)。極小の位置 r=r0r=r_0 が円軌道で、極小のまわりの2次の近似(実線)の曲がり具合から、動径方向の小振動の角振動数が決まる

近似の裏付け

型 B(極限) | 小さいもの 円軌道からのずれ (r−r0)/r0→0(r-r_0)/r_0\to0(UeffU_{\mathrm{eff}} を2次で打ち切る) | 証明 条件付き(UeffU_{\mathrm{eff}} の極小が非退化 Ueff′′(r0)>0U_{\mathrm{eff}}''(r_0)>0。軌道が極小の近くに留まることは HrH_r の保存から言える) | 行き先 近似の数理 第2章

例題 1.8.5 円軌道のまわりの振動と軌道の閉じ方 ★★★

オリジナル質量 mm の粒子が中心力のポテンシャル V(r)V(r) の中を動き、角運動量 pθ≠0p_\theta\neq0 が与えられている。

  1. (1) V=−k/rV=-k/r(k>0k>0)のとき、円軌道の半径 r0r_0 と、円軌道を回る角振動数 ωθ=θ˙\omega_\theta=\dot\theta を求めよ。
  2. (2) 同じ場合に、円軌道からわずかにずれた動径方向の振動の角振動数 ωr=Ueff′′(r0)/m\omega_r=\sqrt{U_{\mathrm{eff}}''(r_0)/m} を求め、ωr=ωθ\omega_r=\omega_\theta を示せ。
  3. (3) V=12κr2V=\frac12\kappa r^2(κ>0\kappa>0)のとき、同じ計算をして ωr=2ωθ\omega_r=2\omega_\theta を示せ。
  4. (4) (2)(3) の結果から、1周の間に rr が何回極小になるかを答え、軌道の形を述べよ。
解答

(1) Ueff=pθ22mr2−krU_{\mathrm{eff}}=\frac{p_\theta^2}{2mr^2}-\frac kr、Ueff′=−pθ2mr3+kr2=0U_{\mathrm{eff}}'=-\frac{p_\theta^2}{mr^3}+\frac k{r^2}=0 から r0=pθ2mkr_0=\dfrac{p_\theta^2}{mk}。ωθ=pθmr02\omega_\theta=\dfrac{p_\theta}{mr_0^2}。

(2) Ueff′′=3pθ2mr4−2kr3U_{\mathrm{eff}}''=\frac{3p_\theta^2}{mr^4}-\frac{2k}{r^3} に k=pθ2/(mr0)k=p_\theta^2/(mr_0) を入れると、r=r0r=r_0 で Ueff′′(r0)=3pθ2mr04−2pθ2mr04=pθ2mr04U_{\mathrm{eff}}''(r_0)=\frac{3p_\theta^2}{mr_0^4}-\frac{2p_\theta^2}{mr_0^4}=\frac{p_\theta^2}{mr_0^4}。よって ωr=Ueff′′(r0)/m=pθmr02=ωθ\omega_r=\sqrt{U''_{\mathrm{eff}}(r_0)/m}=\dfrac{p_\theta}{mr_0^2}=\omega_\theta。(無次元の単位 m=k=pθ=1m=k=p_\theta=1 では r0=1r_0=1, Ueff(r0)=−1/2U_{\mathrm{eff}}(r_0)=-1/2, Ueff′′(r0)=1U''_{\mathrm{eff}}(r_0)=1 で、図 1.8.2 の実線は −12+12(r−1)2-\frac12+\frac12(r-1)^2 である。)

(3) Ueff=pθ22mr2+12κr2U_{\mathrm{eff}}=\frac{p_\theta^2}{2mr^2}+\frac12\kappa r^2、Ueff′=0U'_{\mathrm{eff}}=0 から r04=pθ2/(mκ)r_0^4=p_\theta^2/(m\kappa)。Ueff′′(r0)=3pθ2mr04+κ=4κU''_{\mathrm{eff}}(r_0)=\frac{3p_\theta^2}{mr_0^4}+\kappa=4\kappa なので ωr=2κ/m\omega_r=2\sqrt{\kappa/m}。一方 ωθ=pθ/(mr02)=κ/m\omega_\theta=p_\theta/(mr_0^2)=\sqrt{\kappa/m} なので ωr=2ωθ\omega_r=2\omega_\theta。

(4) 角度が 2π2\pi 進む間に、動径の振動は Kepler 問題で1周期、調和振動子で2周期進む。したがって rr が極小になる回数は1周あたり Kepler 問題で1回、調和振動子で2回で、どちらも1周で元の点に戻る閉じた軌道になる。Kepler 問題では力の中心が楕円の焦点(最も近い点が1つ)、調和振動子では力の中心が楕円の中心(最も近い点が2つ)にある。一般の V(r)V(r) では ωr/ωθ\omega_r/\omega_\theta が有理数とは限らず、近い円軌道は閉じない。

循環座標が見えるかどうかは座標の選び方で決まるが、量が保存するという事実は座標によらない。同じ中心力の問題をデカルト座標で書くと

H=px2+py22m+V(x2+y2)H=\frac{p_x^2+p_y^2}{2m}+V\bigl(\sqrt{x^2+y^2}\bigr)

で、この式は2つの座標のどちらも含むので、循環座標は1つも無い。それでも Lz=xpy−ypxL_z=xp_y-yp_x は運動の間一定で、これは極座標の pθp_\theta と同じ量である。保存量は座標を取り替えても保存量のままだが、それが運動量の1成分として表に出るのは、保存量に合わせて座標を選んだときだけである。

保存量(対称性)に合わせて循環座標を作る手順は、第1部第10章で Noether の定理と合わせて扱う。詳しくは基礎理論『物理のための微分幾何入門』第13章 §13.5(流れの箱と循環座標)で扱う。

木下の勉強ステータス

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

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

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