tecktom.com

授業で習った形との対応

2026年10月15日5分で読めますAI生成教科書『虎の巻』第1部 第8章 Hamilton 形式と相空間
目次
  1. 1階の連立方程式とエネルギー積分
  2. 木下の層

AI 生成・木下が未勉強

解析力学から様々な理論分野を組み立てる 第1部 第8章 §1.8.6 ★★ Newton の方程式を「位置と速度」の1階の連立方程式に直す書き方と、エネルギー積分の導き方は、Hamilton 方程式と dH/dt=∂H/∂t\mathrm{d} H/\mathrm{d} t=\partial H/\partial t の特別な場合

1階の連立方程式とエネルギー積分

力学の授業では、1次元の運動方程式 mx¨=F(x)=−V′(x)m\ddot x=F(x)=-V'(x) を、速度 v=x˙v=\dot x を新しい未知関数として

x˙=v,mv˙=−V′(x)(1.8.20)\dot x=v,\qquad m\dot v=-V'(x)\tag{1.8.20}

と1階の連立方程式に直すことが多い。微分方程式の授業では、(x,v)(x,v) の平面を「相平面」と呼んで解曲線を描く。また、運動方程式の両辺に x˙\dot x を掛けて

mx˙x¨+V′(x)x˙=ddt(12mx˙2+V(x))=0(1.8.21)m\dot x\ddot x+V'(x)\dot x=\frac{\mathrm{d} }{\mathrm{d} t}\Bigl(\frac12m\dot x^2+V(x)\Bigr)=0\tag{1.8.21}

とし、力学的エネルギーの保存を出す(エネルギー積分)。

この章の形と並べると、次の対応になる。

授業での形この章の形
変数位置 xx と速度 v=x˙v=\dot x一般化座標 qq と正準運動量 p=∂L/∂q˙p=\partial L/\partial\dot q
1階の方程式x˙=v\dot x=v, mv˙=Fm\dot v=Fq˙=∂H/∂p\dot q=\partial H/\partial p, p˙=−∂H/∂q\dot p=-\partial H/\partial q
エネルギー12mx˙2+V(x)\frac12m\dot x^2+V(x)H(q,p,t)=∑pq˙−LH(q,p,t)=\sum p\dot q-L
保存の導き方両辺に x˙\dot x を掛けるdHdt=∂H∂t\frac{\mathrm{d} H}{\mathrm{d} t}=\frac{\partial H}{\partial t}(式 (1.8.14))
平面の図相平面 (x,v)(x,v)相空間 (q,p)(q,p)

直交座標の粒子では p=mvp=mv なので、2つの形は質量を掛ける違いしかない。エネルギー積分で x˙\dot x を掛ける操作は、式 (1.8.14) の連鎖律の計算を具体的な HH について行ったものにあたる。一般化座標の運動量が速度に比例しないときは、2つの形が別物になる。極座標の pθ=mr2θ˙p_\theta=mr^2\dot\theta は角速度に rr に依存する係数が掛かったもので、電磁場の中では p=mv+qAp=m\bm{v}+q\bm{A} になる(第2部第2章)。「x˙\dot x を掛ける」操作は、一般の座標では何を掛ければよいかが分からないが、式 (1.8.14) はどの座標でも同じ計算で済む。

相平面 (x,v)(x,v) と相空間 (q,p)(q,p) の違いは、面積を保つかどうかに現れる。次の例題で、質量が位置によって変わる系を使って確かめる。量子力学で p^=−iℏ ∂/∂q\hat p=-\mathrm{i}\hbar\,\partial/\partial q に置き換える量も、速度でなく正準運動量である。Schrödinger 方程式は Hamiltonian H^(q^,p^)\hat H(\hat q,\hat p) から作られる(第4部第2章)。

例題 1.8.7 相平面 (q,q˙)(q,\dot q) と相空間 (q,p)(q,p) の面積 ★★

オリジナルM(q)>0M(q)>0 を滑らかな関数として、自由度1の Lagrangian L=12M(q)q˙2L=\frac12M(q)\dot q^2 を考える(たとえば曲がった針金に沿って動くビーズで、qq を針金に沿った長さ以外の座標にとった場合)。

  1. (1) pp と HH を求め、Hamilton 方程式を書け。
  2. (2) (q,p)(q,p) での速度場の発散が 00 であることを確かめよ。
  3. (3) v=q˙v=\dot q を使って運動方程式を q˙=v\dot q=v, v˙=G(q,v)\dot v=G(q,v) の形に直し、(q,v)(q,v) での速度場の発散を求めよ。
  4. (4) M(q)=1+q2M(q)=1+q^2, (q,v)=(1,1)(q,v)=(1,1) の点で、(3) の発散の値を求めよ。
解答

(1) p=M(q)q˙p=M(q)\dot q, H=p22M(q)H=\dfrac{p^2}{2M(q)}。Hamilton 方程式は q˙=pM\dot q=\dfrac{p}{M}, p˙=−∂H∂q=M′p22M2\dot p=-\frac{\partial H}{\partial q}=\dfrac{M'p^2}{2M^2}。

(2) ∂q˙∂q+∂p˙∂p=−M′pM2+M′pM2=0\frac{\partial \dot q}{\partial q}+\frac{\partial \dot p}{\partial p}=-\dfrac{M'p}{M^2}+\dfrac{M'p}{M^2}=0。

(3) Euler–Lagrange 方程式 ddt(Mq˙)=12M′q˙2\frac{\mathrm{d} }{\mathrm{d} t}(M\dot q)=\frac12M'\dot q^2 から Mv˙+M′v2=12M′v2M\dot v+M'v^2=\frac12M'v^2、すなわち v˙=−M′v22M\dot v=-\dfrac{M'v^2}{2M}。(q,v)(q,v) での発散は ∂v∂q+∂v˙∂v=0−M′vM=−M′(q) vM(q)\frac{\partial v}{\partial q}+\frac{\partial \dot v}{\partial v}=0-\dfrac{M'v}{M}=-\dfrac{M'(q)\,v}{M(q)} で、一般に 00 でない。

(4) M′=2qM'=2q なので発散は −2qv/(1+q2)=−2/2=−1-2qv/(1+q^2)=-2/2=-1。この点の近くの小さな領域は (q,v)(q,v) の平面では面積が単位時間あたり割合 11 で縮むが、(q,p)(q,p) の平面では面積が変わらない。面積を保つ性質は正準運動量 pp を使ったときに成り立つもので、速度 vv を使った平面では一般に成り立たない。

木下の勉強ステータス

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

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

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