解析力学から様々な理論分野を組み立てる 第1部 第11章 §1.11.3 ★★★ d A d t = { A , H } + ∂ A ∂ t \frac{\mathrm{d} A}{\mathrm{d} t}=\{A,H\}+\frac{\partial A}{\partial t} d t d A = { A , H } + ∂ t ∂ A 。{ G , H } = 0 \{G,H\}=0 { G , H } = 0 は「G G G が保存する」と「G G G の変換で H H H が変わらない」を同時に言う
運動方程式を括弧で書く
時間を陽に含まない関数 A ( q , p ) A(q,p) A ( q , p ) については、第1部第9章 の式 (1.9.19) で d A d t = { A , H } \frac{\mathrm{d} A}{\mathrm{d} t}=\{A,H\} d t d A = { A , H } が出ている。ここで足すのは、A A A が時刻を陽に含むときの ∂ A ∂ t \frac{\partial A}{\partial t} ∂ t ∂ A の項だけである。相空間の関数 A ( q , p , t ) A(q,p,t) A ( q , p , t ) を運動の軌道 ( q ( t ) , p ( t ) ) (q(t),p(t)) ( q ( t ) , p ( t )) に沿って時間微分すると、連鎖律と Hamilton 方程式から
d A d t = ∑ i = 1 f ( ∂ A ∂ q i q ˙ i + ∂ A ∂ p i p ˙ i ) + ∂ A ∂ t = { A , H } + ∂ A ∂ t (1.11.10) \frac{\mathrm{d} A}{\mathrm{d} t}=\sum_{i=1}^f\Bigl(\frac{\partial A}{\partial q_i}\dot q_i+\frac{\partial A}{\partial p_i}\dot p_i\Bigr)+\frac{\partial A}{\partial t}
=\{A,H\}+\frac{\partial A}{\partial t}\tag{1.11.10} d t d A = i = 1 ∑ f ( ∂ q i ∂ A q ˙ i + ∂ p i ∂ A p ˙ i ) + ∂ t ∂ A = { A , H } + ∂ t ∂ A ( 1.11.10 )
を得る。第1項は「時刻 d t \mathrm{d} t d t の間に H H H が点を動かした分の変化」で、式 (1.11.3) の G = H G=H G = H , ϵ = d t \epsilon=\mathrm{d} t ϵ = d t の場合にあたる。第2項は A A A の形そのものが時刻とともに変わる分である。Hamilton 方程式 q ˙ i = { q i , H } \dot q_i=\{q_i,H\} q ˙ i = { q i , H } , p ˙ i = { p i , H } \dot p_i=\{p_i,H\} p ˙ i = { p i , H } は、式 (1.11.10) を座標と運動量に当てたものになる。
この式から、時間を陽に含まない関数 A ( q , p ) A(q,p) A ( q , p ) が保存する条件は { A , H } = 0 \{A,H\}=0 { A , H } = 0 である。H H H 自身が時間を陽に含まなければ、d H d t = { H , H } = 0 \frac{\mathrm{d} H}{\mathrm{d} t}=\{H,H\}=0 d t d H = { H , H } = 0 でエネルギーが保存する。
Noether の定理を括弧で書く
G ( q , p ) G(q,p) G ( q , p ) を時間を陽に含まない関数とする。反対称性 (1.11.5) から、次の2つの式は同じ1つの式である。
d G d t = { G , H } = − { H , G } = − δ H ϵ ( δ H は G が生成する変換による H の変化 ) . (1.11.11) \frac{\mathrm{d} G}{\mathrm{d} t}=\{G,H\}=-\{H,G\}=-\frac{\delta H}{\epsilon}\qquad(\delta H\text{ は }G\text{ が生成する変換による }H\text{ の変化}).\tag{1.11.11} d t d G = { G , H } = − { H , G } = − ϵ δ H ( δ H は G が生成する変換による H の変化 ) . ( 1.11.11 )
定理 1.11.3(Hamilton 形式の Noether の定理 ★★★)
時間を陽に含まない Hamiltonian H H H と関数 G ( q , p ) G(q,p) G ( q , p ) について、次の2つは同値である。
(i) G G G は運動の保存量である(d G d t = 0 \frac{\mathrm{d} G}{\mathrm{d} t}=0 d t d G = 0 )。
(ii) G G G が生成する無限小変換で H H H は変わらない(δ H = ϵ { H , G } = 0 \delta H=\epsilon\{H,G\}=0 δ H = ϵ { H , G } = 0 )。
どちらも { G , H } = 0 \{G,H\}=0 { G , H } = 0 と同じである。
この定理は、第1部第10章 の命題 1.10.3 を括弧で書き直したものである。証明は式 (1.11.11) の1行で足りる。
第1部第10章では、Lagrangian L L L を変えない変換 δ q i = ϵ K i ( q ) \delta q_i=\epsilon K_i(q) δ q i = ϵ K i ( q ) から保存量 Q = ∑ i p i K i Q=\sum_ip_iK_i Q = ∑ i p i K i を作った。定理 1.11.3 はその Hamilton 形式の版で、2つの点で形が変わっている。1つ目は、変換が q q q だけでなく p p p も混ぜてよいことである(G G G は q q q と p p p の任意の関数でよい)。2つ目は、向きが両方向になったことである。保存量 G G G が見つかれば、それを生成子にした変換は必ず H H H の対称性になる。Lagrangian 形式で変換を q q q だけの変換に限っている限り、保存量から対称性へ戻る向きは一般には成り立たない。
調和振動子 H = p 2 2 m + 1 2 m ω 2 q 2 H=\frac{p^2}{2m}+\frac12m\omega^2q^2 H = 2 m p 2 + 2 1 m ω 2 q 2 を例にとる。H H H 自身は保存するので、H H H が生成する変換(時間発展)で H H H は変わらない。2次元の等方的な調和振動子では、L z L_z L z のほかに H x − H y H_x-H_y H x − H y (x x x 方向と y y y 方向のエネルギーの差)も保存し、それが生成する変換は q q q と p p p を混ぜる。この変換は座標空間の変換としては書けないので、Lagrangian 形式の Noether の定理(δ q = ϵ K ( q ) \delta q=\epsilon K(q) δ q = ϵK ( q ) の形)からは出てこない。章末の総合問題で扱う。
H H H だけが特別に見える理由
3つの生成子との括弧を並べると、どれも「その変換でずらした向きの微分」になる。
{ F , p x } = ∂ F ∂ x , { F , L z } = ∂ F ∂ θ , { F , H } = d F d t (1.11.12) \{F,p_x\}=\frac{\partial F}{\partial x},\qquad\{F,L_z\}=\frac{\partial F}{\partial \theta},\qquad\{F,H\}=\frac{\mathrm{d} F}{\mathrm{d} t}\tag{1.11.12} { F , p x } = ∂ x ∂ F , { F , L z } = ∂ θ ∂ F , { F , H } = d t d F ( 1.11.12 )
ここで ∂ F / ∂ θ \partial F/\partial\theta ∂ F / ∂ θ は、位置と運動量を一緒に z z z 軸のまわりに回す角度での微分である。最後の式は F F F が時刻を陽に含まない場合に限る。3つの式は同じ形をしていて、括弧の上では H H H だけが特別な位置にあるわけではない。
H H H が特別に見えるのは、運動を追うパラメータに時刻を選んだからである。時刻を選ぶと、時刻をずらす生成子が発展を進める量になり、残りの生成子は発展の間に一定かどうかを調べる量になる。
定理 1.11.3 の式 { p x , H } = 0 \{p_x,H\}=0 { p x , H } = 0 は、この非対称を取り払うと2通りに読める。1つは、時間発展で p x p_x p x が変わらないこと(運動量の保存)である。もう1つは、x x x 方向にずらしても H H H が変わらないこと(並進の対称性)である。
発展のパラメータは時刻でなくてもよい。この巻のほかの章に、次の3つの例がある。
光線 :第2部第10章 では、光軸の座標 z z z を発展のパラメータに選んだ。z z z をずらす生成子は光線の Hamiltonian
H 光 ( x , y , p x , p y ; z ) = − n ( x , y , z ) 2 − p x 2 − p y 2 H_{\text{光}}(x,y,p_x,p_y;z)=-\sqrt{n(x,y,z)^2-p_x^2-p_y^2} H 光 ( x , y , p x , p y ; z ) = − n ( x , y , z ) 2 − p x 2 − p y 2
で、− H 光 -H_{\text{光}} − H 光 が z z z 方向の運動量の役をする。2. 相対論 :第2部第3章 では、エネルギーと運動量が4元運動量 p μ = ( E / c , p ) p^\mu=(E/c,\bm{p}) p μ = ( E / c , p ) にまとまる。時刻とエネルギーの組は、位置と運動量の組と同じ関係にある。3. 量子力学 :時間発展の e − i H ^ t / ℏ \mathrm{e}^{-\mathrm{i}\hat Ht/\hbar} e − i H ^ t /ℏ と並進の e − i p ^ a / ℏ \mathrm{e}^{-\mathrm{i}\hat pa/\hbar} e − i p ^ a /ℏ は同じ形のユニタリ演算子である(第1部第9章 ・第4部第2章 )。
時間を陽に含む保存量と Poisson の定理
G G G が時間を陽に含む場合は、式 (1.11.10) の2つの項が打ち消し合えば保存する。{ G , H } ≠ 0 \{G,H\}\neq0 { G , H } = 0 でも保存量になりうる。自由粒子 H = p 2 / ( 2 m ) H=p^2/(2m) H = p 2 / ( 2 m ) の K = m q − p t K=mq-pt K = m q − pt がその例で、{ K , H } = m { q , H } = p \{K,H\}=m\{q,H\}=p { K , H } = m { q , H } = p 、∂ K ∂ t = − p \frac{\partial K}{\partial t}=-p ∂ t ∂ K = − p なので d K d t = 0 \frac{\mathrm{d} K}{\mathrm{d} t}=0 d t d K = 0 である。K / m = q − p t / m K/m=q-pt/m K / m = q − pt / m は「時刻 0 0 0 での位置」で、等速運動ではこれが一定になる。K K K は Galilei 変換(一定の速度で動く座標系への乗り移り)の生成子で、δ q = ϵ { q , K } = − ϵ t \delta q=\epsilon\{q,K\}=-\epsilon t δ q = ϵ { q , K } = − ϵ t , δ p = ϵ { p , K } = − ϵ m \delta p=\epsilon\{p,K\}=-\epsilon m δ p = ϵ { p , K } = − ϵ m は、速度を − ϵ -\epsilon − ϵ だけずらす変換である。
保存量を2つ持っていれば、その括弧も保存量になる。
定理 1.11.4(Poisson の定理 ★★)
F ( q , p , t ) F(q,p,t) F ( q , p , t ) と G ( q , p , t ) G(q,p,t) G ( q , p , t ) がともに保存量なら、{ F , G } \{F,G\} { F , G } も保存量である。
証明 式 (1.11.10) を { F , G } \{F,G\} { F , G } に当てる。Jacobi 恒等式 (1.11.8) を ( H , F , G ) (H,F,G) ( H , F , G ) に使い反対称性で並べ替えると{ { F , G } , H } = { F , { G , H } } + { { F , H } , G } \{\{F,G\},H\}=\{F,\{G,H\}\}+\{\{F,H\},G\} {{ F , G } , H } = { F , { G , H }} + {{ F , H } , G } となる。時間の偏微分は積の微分と同じ規則で括弧の中に入るので∂ ∂ t { F , G } = { ∂ F ∂ t , G } + { F , ∂ G ∂ t } \frac{\partial }{\partial t}\{F,G\}=\{\frac{\partial F}{\partial t},G\}+\{F,\frac{\partial G}{\partial t}\} ∂ t ∂ { F , G } = { ∂ t ∂ F , G } + { F , ∂ t ∂ G } 。両方を足すと
d d t { F , G } = { F , d G d t } + { d F d t , G } (1.11.13) \frac{\mathrm{d} }{\mathrm{d} t}\{F,G\}=\Bigl\{F,\frac{\mathrm{d} G}{\mathrm{d} t}\Bigr\}+\Bigl\{\frac{\mathrm{d} F}{\mathrm{d} t},G\Bigr\}\tag{1.11.13} d t d { F , G } = { F , d t d G } + { d t d F , G } ( 1.11.13 )
を得る。仮定から右辺は 0 0 0 である。
式 (1.11.13) は、時間微分が括弧に対して積の微分と同じ規則で当たることを言っている。Poisson の定理はその特別な場合である。たとえば中心力の中の粒子では L x L_x L x と L y L_y L y が保存するので、式 (1.11.9) から L z = { L x , L y } L_z=\{L_x,L_y\} L z = { L x , L y } も保存する。ただし、括弧から新しい保存量が出るとは限らない。{ F , G } \{F,G\} { F , G } が 0 0 0 や定数、あるいはすでに知っている保存量の関数になることも多い(自由粒子の K K K と p p p では { K , p } = m \{K,p\}=m { K , p } = m で定数になる)。
相空間に点の集まりを広げたときの密度 ρ ( q , p , t ) \rho(q,p,t) ρ ( q , p , t ) は、点と一緒に流れるので d ρ d t = 0 \frac{\mathrm{d} \rho}{\mathrm{d} t}=0 d t d ρ = 0 、すなわち ∂ ρ ∂ t = − { ρ , H } = { H , ρ } \frac{\partial \rho}{\partial t}=-\{\rho,H\}=\{H,\rho\} ∂ t ∂ ρ = − { ρ , H } = { H , ρ } を満たす。これが Liouville 方程式で、統計力学の出発点になる(第3部第1章 )。
例題 1.11.3 一様な重力の中の Galilei 変換の生成子 ★★
オリジナル 一様な重力の中の1次元の運動 H = p 2 2 m + m g q H=\dfrac{p^2}{2m}+mgq H = 2 m p 2 + m g q (q q q は鉛直上向き)を考える。
(1) 自由粒子の Galilei 変換の生成子 K = m q − p t K=mq-pt K = m q − pt について、d K d t \frac{\mathrm{d} K}{\mathrm{d} t} d t d K を式 (1.11.10) で求めよ。
(2) K ′ = m q − p t − 1 2 m g t 2 K'=mq-pt-\frac12mgt^2 K ′ = m q − pt − 2 1 m g t 2 が保存量であることを示せ。
(3) K ′ K' K ′ の値は何を表すか。運動方程式の解 q ( t ) = q 0 + v 0 t − 1 2 g t 2 q(t)=q_0+v_0t-\frac12gt^2 q ( t ) = q 0 + v 0 t − 2 1 g t 2 を代入して答えよ。
解答 (1) { q , H } = p / m \{q,H\}=p/m { q , H } = p / m , { p , H } = − m g \{p,H\}=-mg { p , H } = − m g から { K , H } = m ⋅ p m − t ( − m g ) = p + m g t \{K,H\}=m\cdot\frac pm-t(-mg)=p+mgt { K , H } = m ⋅ m p − t ( − m g ) = p + m g t 。∂ K ∂ t = − p \frac{\partial K}{\partial t}=-p ∂ t ∂ K = − p なので
d K d t = p + m g t − p = m g t ≠ 0. \frac{\mathrm{d} K}{\mathrm{d} t}=p+mgt-p=mgt\neq0 . d t d K = p + m g t − p = m g t = 0.
重力があると K K K は保存しない。
(2) { K ′ , H } \{K',H\} { K ′ , H } は K K K と同じで p + m g t p+mgt p + m g t (− 1 2 m g t 2 -\frac12mgt^2 − 2 1 m g t 2 は q , p q,p q , p を含まないので括弧に寄与しない)。∂ K ′ ∂ t = − p − m g t \frac{\partial K'}{\partial t}=-p-mgt ∂ t ∂ K ′ = − p − m g t なので d K ′ d t = 0 \frac{\mathrm{d} K'}{\mathrm{d} t}=0 d t d K ′ = 0 。
(3) p = m q ˙ = m ( v 0 − g t ) p=m\dot q=m(v_0-gt) p = m q ˙ = m ( v 0 − g t ) を入れるとK ′ = m ( q 0 + v 0 t − 1 2 g t 2 ) − m ( v 0 − g t ) t − 1 2 m g t 2 = m q 0 K'=m\bigl(q_0+v_0t-\tfrac12gt^2\bigr)-m(v_0-gt)t-\tfrac12mgt^2=mq_0 K ′ = m ( q 0 + v 0 t − 2 1 g t 2 ) − m ( v 0 − g t ) t − 2 1 m g t 2 = m q 0 。K ′ / m K'/m K ′ / m は時刻 0 0 0 の位置 q 0 q_0 q 0 で、どの時刻に計算しても同じ値になる。重力の中でも、軌道を時刻 0 0 0 までさかのぼった位置が保存量になっている。時間を陽に含む保存量は、こうした「初期値を今の値で表したもの」として作ることができる。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com