tecktom.com

この章で選ぶもの:Hamiltonian H(q,p,t)H(q,p,t)

2026年10月15日8分で読めますAI生成教科書『虎の巻』第1部 第8章 Hamilton 形式と相空間
目次
  1. 正準運動量と Hamiltonian
  2. 例:1次元の粒子と調和振動子
  3. 木下の層

AI 生成・木下が未勉強

解析力学から様々な理論分野を組み立てる 第1部 第8章 §1.8.1 ★★★ LL の速度 q˙\dot q を正準運動量 pp に取り替えた関数 H(q,p,t)=∑ipiq˙i−LH(q,p,t)=\sum_ip_i\dot q_i-L を選ぶ

第1部第4章までで、運動は作用 S[q]=∫L(q,q˙,t) dtS[q]=\int L(q,\dot q,t)\,\mathrm{d} t の停留として書けるようになった。停留の条件は Euler–Lagrange 方程式で、一般化座標 qiq_i についての2階の微分方程式が ff 本並ぶ。この章では同じ運動を、一般化座標 qiq_i と正準運動量 pi=∂L/∂q˙ip_i=\partial L/\partial\dot q_i の 2f2f 個の変数で書き直す。使う関数は Lagrangian LL を Legendre 変換して作る Hamiltonian H(q,p,t)=∑ipiq˙i−LH(q,p,t)=\sum_ip_i\dot q_i-L で、運動方程式はq˙i=∂H/∂pi\dot q_i=\partial H/\partial p_i, p˙i=−∂H/∂qi\dot p_i=-\partial H/\partial q_i という1階の方程式 2f2f 本になる(Hamilton 方程式)。

書き直して増えるものは3つある。1つ目は相空間で、(q,p)(q,p) を座標とする 2f2f 次元の空間の各点から軌道が1本ずつ出ていき、自由度1の系なら運動の全体を平面上の図(相図)として見渡すことができる。2つ目は保存量の読み方で、HH が時刻を含まなければ HH が保存し、HH がある座標 qkq_k を含まなければ pkp_k が保存して、問題の自由度が1つ減る。3つ目は流れの性質で、Hamilton 方程式の流れは相空間の面積(一般には 2f2f 次元の体積)を変えない。この性質は第3部の統計力学(Liouville の定理)の出発点になり、数値計算では長時間の積分の方法を選ぶ基準になる。最後に、Newton の運動方程式を1階の連立方程式に直す授業での書き方と、この章の形を並べる。

正準運動量と Hamiltonian

自由度 ff の系の Lagrangian L(q,q˙,t)L(q,\dot q,t) が与えられているとする。第1部第7章の式 (1.7.8) で作った正準運動量を、各 ii について

pi=∂L∂q˙i(q,q˙,t)(i=1,…,f)(1.8.1)p_i=\frac{\partial L}{\partial \dot q_i}(q,\dot q,t)\qquad(i=1,\dots,f)\tag{1.8.1}

と書く(qiq_i に共役な正準運動量)。直交座標の粒子 L=12mx˙2−V(x)L=\frac12m\dot x^2-V(x) では p=mx˙p=m\dot x で、Newton 力学の運動量と同じものになる。極座標の角度 θ\theta に共役な pθp_\theta は角運動量になり(下の例題)、電磁場の中の荷電粒子では pp は mx˙m\dot x と一致しない(第2部第2章)。正準運動量は「LL を速度で偏微分したもの」として定義される量で、速度に質量を掛けたものとは限らない。

式 (1.8.1) は (q,q˙,t)(q,\dot q,t) から pp への対応を与える。qq と tt を止めたとき、この対応が q˙\dot q について解けてq˙i=q˙i(q,p,t)\dot q_i=\dot q_i(q,p,t) と書けることを仮定する。十分条件は、f×ff\times f 行列 ∂2L/∂q˙i∂q˙j\partial^2L/\partial\dot q_i\partial\dot q_j が正則(行列式が 00 でない)であることで、運動エネルギーが速度の正定値2次式である普通の力学系はこれを満たす。この条件を満たす LL を正則な Lagrangian と呼ぶ。

定義 1.8.1(相空間(Hamiltonian は第7章のもの))

正則な Lagrangian L(q,q˙,t)L(q,\dot q,t) に対して、第1部第7章の式 (1.7.9) で作った Hamiltonian を ff 自由度で書くと

H(q,p,t)=∑i=1fpi q˙i(q,p,t)−L(q,q˙(q,p,t),t)(1.8.2)H(q,p,t)=\sum_{i=1}^fp_i\,\dot q_i(q,p,t)-L\bigl(q,\dot q(q,p,t),t\bigr)\tag{1.8.2}

である。右辺の q˙\dot q はすべて (q,p,t)(q,p,t) の関数として書き換えてから代入する。(q1,…,qf,p1,…,pf)(q_1,\dots,q_f,p_1,\dots,p_f) を座標とする 2f2f 次元の空間を相空間と呼び、その点を (q,p)(q,p) と書く。

HH は LL の速度の変数についての Legendre 変換である(第1部第7章)。第1部第7章で示したとおり、Legendre 変換は関数を「接線の傾き pp と切片」で表し直す操作で、LL が q˙\dot q について凸なら情報を失わずに逆に戻すことができる。この章では、第7章の取り替えない変数の規則(式 (1.7.10))を ff 自由度で使う。変換した変数については ∂H/∂pi=q˙i\partial H/\partial p_i=\dot q_i、変換に加わらない変数(qq と tt)については ∂H/∂qi=−∂L/∂qi\partial H/\partial q_i=-\partial L/\partial q_i と ∂H/∂t=−∂L/∂t\partial H/\partial t=-\partial L/\partial t が成り立つ。

例:1次元の粒子と調和振動子

L=12mq˙2−V(q)L=\frac12m\dot q^2-V(q) のとき、p=mq˙p=m\dot q を q˙=p/m\dot q=p/m と解いて式 (1.8.2) に入れると

H=p⋅pm−[12m(pm)2−V(q)]=p22m+V(q)(1.8.3)H=p\cdot\frac pm-\Bigl[\frac12m\Bigl(\frac pm\Bigr)^2-V(q)\Bigr]=\frac{p^2}{2m}+V(q)\tag{1.8.3}

となる。この場合の HH は運動エネルギーとポテンシャルの和で、力学的エネルギーと同じ値をとる。ただし HH は速度でなく運動量の関数として書かれている。調和振動子 V=12mω2q2V=\frac12m\omega^2q^2 では H=p2/(2m)+12mω2q2H=p^2/(2m)+\frac12m\omega^2q^2 で、qq と pp を対等に含む。

HH がいつでも力学的エネルギーに等しいわけではない。H=T+VH=T+V が成り立つ条件と、成り立たないのに HH が保存する例は1.8.4節で扱う。

HH は pp の関数として書く

HH の値だけを見れば、pp で書いても q˙\dot q で書いても同じ数になる。しかし1.8.2節の Hamilton 方程式は、qq についての偏微分を「pp を止めて」とることを要求する。q˙\dot q で書いた式を「q˙\dot q を止めて」偏微分すると、違う量が出る。次の例題の (3) がその例である。

例題 1.8.1 平面の極座標の Hamiltonian ★★★

オリジナル質量 mm の粒子が平面内でポテンシャル V(r)V(r) の中を動く。極座標 (r,θ)(r,\theta) で L=12m(r˙2+r2θ˙2)−V(r)L=\frac12m(\dot r^2+r^2\dot\theta^2)-V(r) である。

  1. (1) 正準運動量 prp_r, pθp_\theta を求め、pθp_\theta が原点のまわりの角運動量であることを確かめよ。
  2. (2) Hamiltonian H(r,θ,pr,pθ)H(r,\theta,p_r,p_\theta) を求めよ。
  3. (3) HH の値を速度で書いた式 Ev(r,r˙,θ˙)=12m(r˙2+r2θ˙2)+V(r)E_{\mathrm v}(r,\dot r,\dot\theta)=\frac12m(\dot r^2+r^2\dot\theta^2)+V(r) について、「r˙,θ˙\dot r,\dot\theta を止めた rr 偏微分」と「pr,pθp_r,p_\theta を止めた rr 偏微分 ∂H/∂r\partial H/\partial r」を比べよ。
解答

(1) pr=∂L/∂r˙=mr˙p_r=\partial L/\partial\dot r=m\dot r、pθ=∂L/∂θ˙=mr2θ˙p_\theta=\partial L/\partial\dot\theta=mr^2\dot\theta。直交座標の位置 r=(rcos⁡θ,rsin⁡θ)\bm{r}=(r\cos\theta,r\sin\theta) と速度から角運動量の zz 成分を計算するとm(xy˙−yx˙)=m rcos⁡θ (r˙sin⁡θ+rθ˙cos⁡θ)−m rsin⁡θ (r˙cos⁡θ−rθ˙sin⁡θ)=mr2θ˙m(x\dot y-y\dot x)=m\,r\cos\theta\,(\dot r\sin\theta+r\dot\theta\cos\theta)-m\,r\sin\theta\,(\dot r\cos\theta-r\dot\theta\sin\theta)=mr^2\dot\theta で、pθp_\theta と一致する。

(2) r˙=pr/m\dot r=p_r/m, θ˙=pθ/(mr2)\dot\theta=p_\theta/(mr^2) を式 (1.8.2) に入れると、∑pq˙=pr2/m+pθ2/(mr2)\sum p\dot q=p_r^2/m+p_\theta^2/(mr^2) とL=pr2/(2m)+pθ2/(2mr2)−V(r)L=p_r^2/(2m)+p_\theta^2/(2mr^2)-V(r) の差から

H=pr22m+pθ22mr2+V(r).H=\frac{p_r^2}{2m}+\frac{p_\theta^2}{2mr^2}+V(r).

(3) 速度を止めると ∂Ev/∂r=mrθ˙2+V′(r)\partial E_{\mathrm v}/\partial r=mr\dot\theta^2+V'(r)。運動量を止めると

∂H∂r=−pθ2mr3+V′(r)=−mrθ˙2+V′(r).\frac{\partial H}{\partial r}=-\frac{p_\theta^2}{mr^3}+V'(r)=-mr\dot\theta^2+V'(r).

第1項の符号が逆になる。pθp_\theta を止めて rr を大きくすると θ˙=pθ/(mr2)\dot\theta=p_\theta/(mr^2) が小さくなり、回転の運動エネルギー pθ2/(2mr2)p_\theta^2/(2mr^2) が減るためである。Hamilton 方程式で使う正しい方は後者で、p˙r=−∂H/∂r=mrθ˙2−V′(r)\dot p_r=-\partial H/\partial r=mr\dot\theta^2-V'(r) の第1項が遠心力になる。前者を使うと遠心力の向きが逆になる。

木下の勉強ステータス

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

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

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