tecktom.com

この章のまとめ

2026年10月15日10分で読めますAI生成教科書『虎の巻』第1部 第12章 正準変換:時間発展から入って母関数へ

AI 生成・木下が未勉強

解析力学から様々な理論分野を組み立てる 第1部 第12章 §1.12.7 ★★★ 時間発展 → 括弧を保つ変換 → 面積の保存 → 母関数 → K=H+∂F/∂tK=H+\partial F/\partial t の順に、正準変換を組み立てた

時間発展を出発点から行き先への写像(流れ)として読むと、自由粒子・調和振動子・振り子のどの例でも相空間の面積が保たれ、新しい変数の括弧が {Q,P}=1\{Q,P\}=1 のままだった。この性質、すなわち {Qi,Qj}=0\{Q_i,Q_j\}=0, {Pi,Pj}=0\{P_i,P_j\}=0, {Qi,Pj}=δij\{Q_i,P_j\}=\delta_{ij} を正準変換の定義に取った。行列で書くと Jacobi 行列がシンプレクティック MTJM=JM^{\mathsf T}\mathsf JM=\mathsf J であることで、1自由度では det⁡M=1\det M=1 と同じである。正準変換は Poisson 括弧を変えないので、時刻を含まない正準変換では HH を書き直しただけの KK について Hamilton 方程式の形が保たれる。無限小変換 δqi=ϵ ∂G/∂pi\delta q_i=\epsilon\,\partial G/\partial p_i, δpi=−ϵ ∂G/∂qi\delta p_i=-\epsilon\,\partial G/\partial q_i は1次まで正準で、その積み重ねである流れは、変分方程式から MTJMM^{\mathsf T}\mathsf JM が時刻によらないことを示して、厳密に正準と分かった。系として相空間の体積が保たれる(Liouville の定理)。時刻を刻んだ数値解法の1歩は正準とは限らず、正準でない Euler 法ではエネルギーが増え続ける。面積を保つことを線積分に直すと、∑ipi dqi−∑iPi dQi\sum_ip_i\,\mathrm{d} q_i-\sum_iP_i\,\mathrm{d} Q_i が全微分 dF\mathrm{d} F になることと同値で、この FF が母関数である。母関数の4つの型は互いに Legendre 変換で移り合い、恒等変換からのずれ F2=∑iqiPi+ϵGF_2=\sum_iq_iP_i+\epsilon G の1次の部分が生成子 GG である。変換が時刻を含むときは K=H+∂F/∂tK=H+\partial F/\partial t で、K=0K=0 にする母関数を探すことが次の章の Hamilton–Jacobi 方程式(第1部第13章)、運動を一様な回転にする変数を探すことが作用角変数(第1部第14章)になる。

この章の道具

  • 正準変換:{Qi,Qj}={Pi,Pj}=0\{Q_i,Q_j\}=\{P_i,P_j\}=0, {Qi,Pj}=δij\{Q_i,P_j\}=\delta_{ij}   ⟺  \iff MTJM=JM^{\mathsf T}\mathsf JM=\mathsf J(1自由度では det⁡M=1\det M=1)。括弧はどちらの変数で計算しても同じ。
  • 流れ φt\varphi_t は任意の HH で正準変換。系として det⁡M=1\det M=1(Liouville の定理)。
  • 母関数:∑ipi dqi−∑iPi dQi=dF1(q,Q)\sum_ip_i\,\mathrm{d} q_i-\sum_iP_i\,\mathrm{d} Q_i=\mathrm{d} F_1(q,Q)、pi=∂F1/∂qip_i=\partial F_1/\partial q_i, Pi=−∂F1/∂QiP_i=-\partial F_1/\partial Q_i。F2=F1+∑iQiPiF_2=F_1+\sum_iQ_iP_i で pi=∂F2/∂qip_i=\partial F_2/\partial q_i, Qi=∂F2/∂PiQ_i=\partial F_2/\partial P_i。
  • 時刻を含む変換では K=H+∂F/∂tK=H+\partial F/\partial t。生成子は F2=∑iqiPi+ϵGF_2=\sum_iq_iP_i+\epsilon G の1次の部分。

総合問題 調和振動子の流れを母関数で書く ★★★

オリジナル調和振動子 H=p22m+12mω2q2H=\dfrac{p^2}{2m}+\dfrac12m\omega^2q^2 の時刻 00 から tt までの流れ (q,p)↦(Q,P)(q,p)\mapsto(Q,P) を考える(0<ωt<π0<\omega t<\pi)。

  1. (1) (使う道具:1.12.1節・1.12.2節/第1部第11章 の dAdt={A,H}\frac{\mathrm{d} A}{\mathrm{d} t}=\{A,H\})QQ, PP を qq, pp で表し、{Q,P}=1\{Q,P\}=1 を確かめよ。
  2. (2) (使う道具:1.12.4節の F1F_1)この流れの母関数 F1(q,Q,t)F_1(q,Q,t) を求めよ。
  3. (3) (使う道具:第1部第3章 の総合問題 (6) の部分積分)出発点 qq から時間 tt で QQ に着く実際の運動に沿った作用 Scl=∫0tL dτS_{\mathrm{cl}}=\int_0^tL\,\mathrm{d}\tau(L=12mq˙2−12mω2q2L=\frac12m\dot q^2-\frac12m\omega^2q^2)がScl=12(QP−qp)S_{\mathrm{cl}}=\frac12(QP-qp) と書けることを示し、F1=−SclF_1=-S_{\mathrm{cl}} を確かめよ。m=ω=1m=\omega=1, q=1q=1, Q=0Q=0, t=π/3t=\pi/3 での値を求めよ。
  4. (4) (使う道具:1.12.4節の Legendre 変換と無限小変換の母関数/第1部第9章)F2(q,P,t)=F1+QPF_2(q,P,t)=F_1+QP を qq, PP で表し、tt が小さいとき F2=qP+tH(q,P)+O(t2)F_2=qP+tH(q,P)+\mathcal{O}(t^2) となることを示せ。これは何を意味するか。
  5. (5) (使う道具:1.12.3節の Liouville の定理・1.12.5節の F1=12mωq2cot⁡QF_1=\frac12m\omega q^2\cot Q)H≤EH\le E の領域の面積を (q,p)(q,p) 平面で求めよ。次に、1.12.5節の変換 (q,p)↦(Q~,P~)(q,p)\mapsto(\tilde Q,\tilde P)(P~=H/ω\tilde P=H/\omega)で同じ領域が (Q~,P~)(\tilde Q,\tilde P) 平面のどんな領域に移るかを答え、面積が等しいことを確かめよ。また、この新しい変数で時刻 tt の流れがどんな写像になるかを答えよ。
解答

(1) 式 (1.12.2) のとおり Q=qcos⁡ωt+pmωsin⁡ωtQ=q\cos\omega t+\dfrac{p}{m\omega}\sin\omega t, P=pcos⁡ωt−mωqsin⁡ωtP=p\cos\omega t-m\omega q\sin\omega t。実際、dQdt={Q,H}\frac{\mathrm{d} Q}{\mathrm{d} t}=\{Q,H\} にあたる Q˙=P/m\dot Q=P/m, P˙=−mω2Q\dot P=-m\omega^2Q がこの式で成り立つ。{Q,P}=cos⁡ωt⋅cos⁡ωt−sin⁡ωtmω⋅(−mωsin⁡ωt)=1\{Q,P\}=\cos\omega t\cdot\cos\omega t-\dfrac{\sin\omega t}{m\omega}\cdot(-m\omega\sin\omega t)=1。

(2) QQ の式を pp について解くと p=mω(Q−qcos⁡ωt)sin⁡ωtp=\dfrac{m\omega(Q-q\cos\omega t)}{\sin\omega t}。これが ∂F1/∂q\partial F_1/\partial q なので、qq で積分して

F1(q,Q,t)=mωsin⁡ωt(qQ−q2+Q22cos⁡ωt)F_1(q,Q,t)=\frac{m\omega}{\sin\omega t}\Bigl(qQ-\frac{q^2+Q^2}{2}\cos\omega t\Bigr)

(QQ だけの関数の分の不定性は、次の PP の式で決まる)。確かめると P=−∂F1/∂Q=mω(Qcos⁡ωt−q)sin⁡ωtP=-\partial F_1/\partial Q=\dfrac{m\omega(Q\cos\omega t-q)}{\sin\omega t} で、pp の式を P=pcos⁡ωt−mωqsin⁡ωtP=p\cos\omega t-m\omega q\sin\omega t に入れたものと一致する(cos⁡2+sin⁡2=1\cos^2+\sin^2=1 を使う)。0<ωt<π0<\omega t<\pi では sin⁡ωt≠0\sin\omega t\neq0 なので qq と QQ が独立に動かせるので、F1F_1 を使うことができる。

(3) 第1部第3章の総合問題 (6) と同じく、q˙2=ddτ(qq˙)−qq¨\dot q^2=\frac{\mathrm{d} }{\mathrm{d} \tau}(q\dot q)-q\ddot q で書き換えると、運動方程式 q¨+ω2q=0\ddot q+\omega^2q=0 を満たす経路では

Scl=m2[qq˙]0t−m2∫0tq(q¨+ω2q)dτ=12(Q P−q p)S_{\mathrm{cl}}=\frac m2\bigl[q\dot q\bigr]_0^t-\frac m2\int_0^tq\bigl(\ddot q+\omega^2q\bigr)\mathrm{d}\tau=\frac12\bigl(Q\,P-q\,p\bigr)

(mq˙m\dot q は運動量で、終点で PP、始点で pp)。(2) の pp, PP を入れると

Scl=mω2sin⁡ωt(Q2cos⁡ωt−qQ−qQ+q2cos⁡ωt)=mω2sin⁡ωt((q2+Q2)cos⁡ωt−2qQ)=−F1.S_{\mathrm{cl}}=\frac{m\omega}{2\sin\omega t}\Bigl(Q^2\cos\omega t-qQ-qQ+q^2\cos\omega t\Bigr)=\frac{m\omega}{2\sin\omega t}\Bigl((q^2+Q^2)\cos\omega t-2qQ\Bigr)=-F_1 .

流れの母関数は、実際の運動に沿った作用の符号を変えたものである。m=ω=1m=\omega=1, q=1q=1, Q=0Q=0, t=π/3t=\pi/3 ではF1=1sin⁡(π/3)(0−12cos⁡π3)=23⋅(−14)=−36≈−0.289F_1=\dfrac{1}{\sin(\pi/3)}\Bigl(0-\dfrac12\cos\dfrac\pi3\Bigr)=\dfrac{2}{\sqrt3}\cdot\Bigl(-\dfrac14\Bigr)=-\dfrac{\sqrt3}{6}\approx-0.289、Scl≈0.289S_{\mathrm{cl}}\approx0.289。作用を端の点の関数と見ることは第1部第13章の Hamilton の主関数 S(q,t)S(q,t) につながる。

(4) P=mω(Qcos⁡ωt−q)sin⁡ωtP=\dfrac{m\omega(Q\cos\omega t-q)}{\sin\omega t} を QQ について解くと Q=qcos⁡ωt+Ptan⁡ωtmωQ=\dfrac{q}{\cos\omega t}+\dfrac{P\tan\omega t}{m\omega}。これを F1+QPF_1+QP に入れて整理すると

F2(q,P,t)=qPcos⁡ωt+tan⁡ωt2(mωq2+P2mω)F_2(q,P,t)=\frac{qP}{\cos\omega t}+\frac{\tan\omega t}{2}\Bigl(m\omega q^2+\frac{P^2}{m\omega}\Bigr)

(∂F2/∂q=P+mωqsin⁡ωtcos⁡ωt=p\partial F_2/\partial q=\dfrac{P+m\omega q\sin\omega t}{\cos\omega t}=p, ∂F2/∂P=Q\partial F_2/\partial P=Q で確かめることができる。ωt=π/2\omega t=\pi/2 では cos⁡ωt=0\cos\omega t=0 で F2F_2 を使うことができない。Q=p/(mω)Q=p/(m\omega) と P=−mωqP=-m\omega q になり qq と PP が独立でなくなるからである)。t→0t\to0 で 1/cos⁡ωt=1+O(t2)1/\cos\omega t=1+\mathcal{O}(t^2), tan⁡ωt=ωt+O(t3)\tan\omega t=\omega t+\mathcal{O}(t^3) なので

F2=qP+t(12mω2q2+P22m)+O(t2)=qP+t H(q,P)+O(t2).F_2=qP+t\Bigl(\frac12m\omega^2q^2+\frac{P^2}{2m}\Bigr)+\mathcal{O}(t^2)=qP+t\,H(q,P)+\mathcal{O}(t^2).

式 (1.12.16) と比べると、ϵ=t\epsilon=t, G=HG=H の無限小変換の母関数である。時間発展の生成子が HH であること(第1部第9章)を、有限の時間の流れの母関数から読み取ることができる。

(5) H≤EH\le E は半軸 2E/(mω2)\sqrt{2E/(m\omega^2)}(qq 方向)と 2mE\sqrt{2mE}(pp 方向)の楕円の内部で、面積は π2E/(mω2)2mE=2πE/ω\pi\sqrt{2E/(m\omega^2)}\sqrt{2mE}=2\pi E/\omega。1.12.5節の変換では H=ωP~H=\omega\tilde P なので、H≤EH\le E は 0≤P~≤E/ω0\le\tilde P\le E/\omega になり、Q~\tilde Q は 0≤Q~<2π0\le\tilde Q<2\pi を1周する。(Q~,P~)(\tilde Q,\tilde P) 平面では幅 2π2\pi・高さ E/ωE/\omega の帯で、面積は 2πE/ω2\pi E/\omega。楕円の面積と等しく、正準変換で面積が保たれることと合う。新しい変数では K=ωP~K=\omega\tilde P なので、時刻 tt の流れは (Q~,P~)↦(Q~+ωt,P~)(\tilde Q,\tilde P)\mapsto(\tilde Q+\omega t,\tilde P) で、帯を横に ωt\omega t だけずらす平行移動になる。(1) で (q,p)(q,p) 平面の回転だった流れが、変数を取り替えると平行移動になる。面積を保つことは Q~\tilde Q 方向の平行移動では明らかで、正準変換で変数を移しても保たれることが Liouville の定理の中身である。

木下の勉強ステータス

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

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

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