解析力学から様々な理論分野を組み立てる 第1部 第8章 §1.8.4 ★★★ H H H が t t t を含まなければ H H H が保存し、q k q_k q k を含まなければ p k p_k p k が保存して自由度が1つ減る
H H H の時間変化
運動に沿った H H H の変化を連鎖律で計算し、Hamilton 方程式 (1.8.6) を入れると
d H d t = ∑ i ( ∂ H ∂ q i q ˙ i + ∂ H ∂ p i p ˙ i ) + ∂ H ∂ t = ∑ i ( ∂ H ∂ q i ∂ H ∂ p i − ∂ H ∂ p i ∂ H ∂ q i ) + ∂ 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} d t d H = i ∑ ( ∂ q i ∂ H q ˙ i + ∂ p i ∂ H p ˙ i ) + ∂ t ∂ H = i ∑ ( ∂ q i ∂ H ∂ p i ∂ H − ∂ p i ∂ H ∂ q i ∂ H ) + ∂ t ∂ H = ∂ t ∂ H ( 1.8.14 )
を得る。和の中の2項は打ち消し合う。したがって H H H が時刻 t t t をあらわに含まなければ、H H H の値は運動の間一定である。式 (1.8.5) の第3式から ∂ H / ∂ t = − ∂ L / ∂ t \partial H/\partial t=-\partial L/\partial t ∂ H / ∂ t = − ∂ L / ∂ t なので、これは「L L L が t t t をあらわに含まなければ H H H が保存する」と言っても同じである。時刻をずらす変換と H H H の保存のつながりは、第1部第10章の Noether の定理で扱う(第1部第10章 )。
H H H が力学的エネルギーに等しい条件
保存する H H H が力学的エネルギー T + V T+V T + V に等しいとは限らない。等しくなる条件を調べる。運動エネルギーが速度の2次の同次式 T = 1 2 ∑ i , j a i j ( q ) q ˙ i q ˙ j T=\frac12\sum_{i,j}a_{ij}(q)\dot q_i\dot q_j T = 2 1 ∑ i , j a ij ( q ) q ˙ i q ˙ j (a i j = a j i a_{ij}=a_{ji} a ij = a j i )で、V V V が速度を含まないとする。このとき p i = ∑ j a i j q ˙ j p_i=\sum_ja_{ij}\dot q_j p i = ∑ j a ij q ˙ j なので ∑ i p i q ˙ i = ∑ i , j a i j q ˙ i q ˙ j = 2 T \sum_ip_i\dot q_i=\sum_{i,j}a_{ij}\dot q_i\dot q_j=2T ∑ i p i q ˙ i = ∑ i , j a ij q ˙ i q ˙ j = 2 T となり、
H = ∑ i p i q ˙ i − L = 2 T − ( T − V ) = T + V (1.8.15) H=\sum_ip_i\dot q_i-L=2T-(T-V)=T+V\tag{1.8.15} H = i ∑ p i q ˙ i − L = 2 T − ( T − V ) = T + V ( 1.8.15 )
が成り立つ。直交座標から一般化座標への変換が時刻を含まない(拘束が時間によらない)なら、T T T はこの形になる。
変換が時刻を含むと、T T T に速度の1次の項と0次の項が入り、∑ p q ˙ = 2 T \sum p\dot q=2T ∑ p q ˙ = 2 T が成り立たなくなる。たとえば一定の角速度で回される枠の上で測った座標を使うと、H H H は T + V T+V T + V と違う関数になる。その H H H が時刻を含まなければ H H H は保存するが、T + V T+V T + V は保存しない(枠を回し続ける外部の装置が系に仕事をするため)。
回る枠の上で時刻によらないポテンシャルの中を動く粒子(天体力学の制限三体問題で、一定の角速度で回る2天体の重力を受ける小天体)がその例である。枠の角速度を Ω \Omega Ω とする。慣性系で測ったエネルギーと角運動量を E E E と L z L_z L z と書くと、Hamiltonian は
H = E − Ω L z H=E-\Omega L_z H = E − Ω L z
になる。E E E も L z L_z L z も別々には保存しないが、この組み合わせは保存し、Jacobi 積分と呼ばれる。
章末の総合問題で、回転する輪にはめたビーズを使ってこれを確かめる。逆に、H = T + V H=T+V H = T + V であっても、ポテンシャルが時刻を含めば H H H は保存しない。H H H が保存するかどうかと H H H がエネルギーに等しいかどうかは、別の条件で決まる。
加速する台から測った座標では、回る枠と逆向きの食い違いが起こる。直線上の自由な粒子を、一定の加速度 a a a で動く台から測った座標 q q q で書く(第1部第2章 の例と同じ座標)。このとき
x = q + 1 2 a t 2 , L = m 2 ( q ˙ + a t ) 2 , p = m ( q ˙ + a t ) = m x ˙ , H = p q ˙ − L = p 2 2 m − a t 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} x = q + 2 1 a t 2 , L = 2 m ( q ˙ + a t ) 2 , p = m ( q ˙ + a t ) = m x ˙ , H = p q ˙ − L = 2 m p 2 − a t p
で、Hamiltonian は時刻を陽に含む。
この Hamiltonian は座標を含まないので、運動量 p p p は保存する。一方で式 (1.8.14) から d H d t = ∂ H / ∂ t = − a p \frac{\mathrm{d} H}{\mathrm{d} t}=\partial H/\partial t=-ap d t d H = ∂ H / ∂ t = − a p で、H H H は時刻に比例して変わる。慣性系で測ったエネルギー 1 2 m x ˙ 2 = p 2 / ( 2 m ) \frac12m\dot x^2=p^2/(2m) 2 1 m x ˙ 2 = p 2 / ( 2 m ) は一定のままなので、H H H の変化は台の座標を選んだことから来ている。回る枠では Hamiltonian が保存して T + V T+V T + V が保存せず、加速する台ではエネルギーが保存して Hamiltonian が保存しないので、2つの例で保存する量が入れ替わる。
循環座標
H H H がある座標 q k q_k q k を含まないとき、q k q_k q k を循環座標 と呼ぶ。Hamilton 方程式から p ˙ k = − ∂ H / ∂ q k = 0 \dot p_k=-\partial H/\partial q_k=0 p ˙ k = − ∂ H / ∂ q k = 0 なので、p k p_k p k は保存する。保存する p k p_k p k の値を定数 c k c_k c k として H H H に入れると、H H H は残りの f − 1 f-1 f − 1 個の座標と運動量、そして定数 c k c_k c k の関数になる。残りの 2 ( f − 1 ) 2(f-1) 2 ( f − 1 ) 本の方程式は q k q_k q k を含まないので、それだけで解くことができ、q k q_k q k は最後に q ˙ k = ∂ H / ∂ c k \dot q_k=\partial H/\partial c_k q ˙ k = ∂ H / ∂ c k を積分して求める。Lagrangian の形では、循環座標 q k q_k q k が L L L に入らなくても速度 q ˙ k \dot q_k q ˙ k は L L L に残るので、このように1変数を定数に置き換えることはできない。Hamilton の形では、循環座標に共役な運動量がそのまま定数になり、問題の自由度が1つ減る。
中心力の問題がその例である。例題 1.8.1 の H = p r 2 2 m + p θ 2 2 m r 2 + V ( r ) H=\frac{p_r^2}{2m}+\frac{p_\theta^2}{2mr^2}+V(r) H = 2 m p r 2 + 2 m r 2 p θ 2 + V ( r ) は θ \theta θ を含まないので p θ p_\theta p θ が保存し、動径方向の運動は自由度1の Hamiltonian
H r ( r , p r ) = p r 2 2 m + U e f f ( r ) , U e f f ( r ) = p θ 2 2 m r 2 + 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} H r ( r , p r ) = 2 m p r 2 + U eff ( r ) , U eff ( r ) = 2 m r 2 p θ 2 + V ( r ) ( 1.8.16 )
で決まる。U e f f U_{\mathrm{eff}} U eff を有効ポテンシャル と呼ぶ。第1項は、p θ p_\theta p θ を一定に保ったまま r r r を変えると回転の運動エネルギーが変わることを表し、例題 1.8.1 (3) の遠心力はこの項の − ∂ / ∂ r -\partial/\partial r − ∂ / ∂ r である。動径の運動が求まれば、角度は θ ˙ = ∂ H / ∂ p θ = p θ / ( m r 2 ) \dot\theta=\partial H/\partial p_\theta=p_\theta/(mr^2) θ ˙ = ∂ H / ∂ p θ = p θ / ( m r 2 ) を積分して求める。
図 1.8.2 Kepler 問題の有効ポテンシャル(r 0 = p θ 2 / ( m k ) r_0=p_\theta^2/(mk) r 0 = p θ 2 / ( mk ) を長さの単位とする)。極小の位置 r = r 0 r=r_0 r = r 0 が円軌道で、極小のまわりの2次の近似(実線)の曲がり具合から、動径方向の小振動の角振動数が決まる
近似の裏付け
型 B(極限) | 小さいもの 円軌道からのずれ ( r − r 0 ) / r 0 → 0 (r-r_0)/r_0\to0 ( r − r 0 ) / r 0 → 0 (U e f f U_{\mathrm{eff}} U eff を2次で打ち切る) | 証明 条件付き(U e f f U_{\mathrm{eff}} U eff の極小が非退化 U e f f ′ ′ ( r 0 ) > 0 U_{\mathrm{eff}}''(r_0)>0 U eff ′′ ( r 0 ) > 0 。軌道が極小の近くに留まることは H r H_r H r の保存から言える) | 行き先 近似の数理 第2章
例題 1.8.5 円軌道のまわりの振動と軌道の閉じ方 ★★★
オリジナル 質量 m m m の粒子が中心力のポテンシャル V ( r ) V(r) V ( r ) の中を動き、角運動量 p θ ≠ 0 p_\theta\neq0 p θ = 0 が与えられている。
(1) V = − k / r V=-k/r V = − k / r (k > 0 k>0 k > 0 )のとき、円軌道の半径 r 0 r_0 r 0 と、円軌道を回る角振動数 ω θ = θ ˙ \omega_\theta=\dot\theta ω θ = θ ˙ を求めよ。
(2) 同じ場合に、円軌道からわずかにずれた動径方向の振動の角振動数 ω r = U e f f ′ ′ ( r 0 ) / m \omega_r=\sqrt{U_{\mathrm{eff}}''(r_0)/m} ω r = U eff ′′ ( r 0 ) / m を求め、ω r = ω θ \omega_r=\omega_\theta ω r = ω θ を示せ。
(3) V = 1 2 κ r 2 V=\frac12\kappa r^2 V = 2 1 κ r 2 (κ > 0 \kappa>0 κ > 0 )のとき、同じ計算をして ω r = 2 ω θ \omega_r=2\omega_\theta ω r = 2 ω θ を示せ。
(4) (2)(3) の結果から、1周の間に r r r が何回極小になるかを答え、軌道の形を述べよ。
解答 (1) U e f f = p θ 2 2 m r 2 − k r U_{\mathrm{eff}}=\frac{p_\theta^2}{2mr^2}-\frac kr U eff = 2 m r 2 p θ 2 − r k 、U e f f ′ = − p θ 2 m r 3 + k r 2 = 0 U_{\mathrm{eff}}'=-\frac{p_\theta^2}{mr^3}+\frac k{r^2}=0 U eff ′ = − m r 3 p θ 2 + r 2 k = 0 から r 0 = p θ 2 m k r_0=\dfrac{p_\theta^2}{mk} r 0 = mk p θ 2 。ω θ = p θ m r 0 2 \omega_\theta=\dfrac{p_\theta}{mr_0^2} ω θ = m r 0 2 p θ 。
(2) U e f f ′ ′ = 3 p θ 2 m r 4 − 2 k r 3 U_{\mathrm{eff}}''=\frac{3p_\theta^2}{mr^4}-\frac{2k}{r^3} U eff ′′ = m r 4 3 p θ 2 − r 3 2 k に k = p θ 2 / ( m r 0 ) k=p_\theta^2/(mr_0) k = p θ 2 / ( m r 0 ) を入れると、r = r 0 r=r_0 r = r 0 で U e f f ′ ′ ( r 0 ) = 3 p θ 2 m r 0 4 − 2 p θ 2 m r 0 4 = p θ 2 m r 0 4 U_{\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} U eff ′′ ( r 0 ) = m r 0 4 3 p θ 2 − m r 0 4 2 p θ 2 = m r 0 4 p θ 2 。よって ω r = U e f f ′ ′ ( r 0 ) / m = p θ m r 0 2 = ω θ \omega_r=\sqrt{U''_{\mathrm{eff}}(r_0)/m}=\dfrac{p_\theta}{mr_0^2}=\omega_\theta ω r = U eff ′′ ( r 0 ) / m = m r 0 2 p θ = ω θ 。(無次元の単位 m = k = p θ = 1 m=k=p_\theta=1 m = k = p θ = 1 では r 0 = 1 r_0=1 r 0 = 1 , U e f f ( r 0 ) = − 1 / 2 U_{\mathrm{eff}}(r_0)=-1/2 U eff ( r 0 ) = − 1/2 , U e f f ′ ′ ( r 0 ) = 1 U''_{\mathrm{eff}}(r_0)=1 U eff ′′ ( r 0 ) = 1 で、図 1.8.2 の実線は − 1 2 + 1 2 ( r − 1 ) 2 -\frac12+\frac12(r-1)^2 − 2 1 + 2 1 ( r − 1 ) 2 である。)
(3) U e f f = p θ 2 2 m r 2 + 1 2 κ r 2 U_{\mathrm{eff}}=\frac{p_\theta^2}{2mr^2}+\frac12\kappa r^2 U eff = 2 m r 2 p θ 2 + 2 1 κ r 2 、U e f f ′ = 0 U'_{\mathrm{eff}}=0 U eff ′ = 0 から r 0 4 = p θ 2 / ( m κ ) r_0^4=p_\theta^2/(m\kappa) r 0 4 = p θ 2 / ( mκ ) 。U e f f ′ ′ ( r 0 ) = 3 p θ 2 m r 0 4 + κ = 4 κ U''_{\mathrm{eff}}(r_0)=\frac{3p_\theta^2}{mr_0^4}+\kappa=4\kappa U eff ′′ ( r 0 ) = m r 0 4 3 p θ 2 + κ = 4 κ なので ω r = 2 κ / m \omega_r=2\sqrt{\kappa/m} ω r = 2 κ / m 。一方 ω θ = p θ / ( m r 0 2 ) = κ / m \omega_\theta=p_\theta/(mr_0^2)=\sqrt{\kappa/m} ω θ = p θ / ( m r 0 2 ) = κ / m なので ω r = 2 ω θ \omega_r=2\omega_\theta ω r = 2 ω θ 。
(4) 角度が 2 π 2\pi 2 π 進む間に、動径の振動は Kepler 問題で1周期、調和振動子で2周期進む。したがって r r r が極小になる回数は1周あたり Kepler 問題で1回、調和振動子で2回で、どちらも1周で元の点に戻る閉じた軌道になる。Kepler 問題では力の中心が楕円の焦点(最も近い点が1つ)、調和振動子では力の中心が楕円の中心(最も近い点が2つ)にある。一般の V ( r ) V(r) V ( r ) では ω r / ω θ \omega_r/\omega_\theta ω r / ω θ が有理数とは限らず、近い円軌道は閉じない。
循環座標が見えるかどうかは座標の選び方で決まるが、量が保存するという事実は座標によらない。同じ中心力の問題をデカルト座標で書くと
H = p x 2 + p y 2 2 m + V ( x 2 + y 2 ) H=\frac{p_x^2+p_y^2}{2m}+V\bigl(\sqrt{x^2+y^2}\bigr) H = 2 m p x 2 + p y 2 + V ( x 2 + y 2 )
で、この式は2つの座標のどちらも含むので、循環座標は1つも無い。それでも L z = x p y − y p x L_z=xp_y-yp_x L z = x p y − y p x は運動の間一定で、これは極座標の p θ p_\theta p θ と同じ量である。保存量は座標を取り替えても保存量のままだが、それが運動量の1成分として表に出るのは、保存量に合わせて座標を選んだときだけである。
保存量(対称性)に合わせて循環座標を作る手順は、第1部第10章 で Noether の定理と合わせて扱う。詳しくは基礎理論『物理のための微分幾何入門』第13章 §13.5(流れの箱と循環座標)で扱う。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com