tecktom.com

Hamilton–Jacobi と作用の波面

2026年10月22日17分で読めますAI生成教科書『基礎理論』第14章 シンプレクティック幾何と解析力学
目次
  1. 1本の解と解の族
  2. Lagrange 部分多様体
  3. Hamilton–Jacobi 方程式
  4. 波面・焦線・ほかの分野へのつながり
  5. 確認テスト
  6. 木下の層

AI 生成・木下が未勉強

物理のための微分幾何入門 第14章 §14.5 ★★★ 運動量が (q,t)(q,t) の関数(場)になるのは、解の族を1つ選んだときだけである

1本の解と解の族

Hamilton 形式では、位置と運動量は相空間の独立な座標で、どちらも相手の関数ではない。1本の解は相空間の曲線 t↦(q(t),p(t))t\mapsto(q(t),p(t)) で、各時刻に相空間の1点を指す。同じ位置を同じ時刻に通る解は、運動量の値を変えれば無数にある。したがって「位置 qq での運動量」という量は、運動方程式だけからは決まらない。

一方、Hamilton–Jacobi 方程式では p=∂S/∂qp=\partial S/\partial q と書き、運動量を位置と時刻の関数として扱う。これは配位空間の上の場である。この食い違いは、解を1本でなく族で考えると解消する。

nn 個のパラメータ a=(a1,…,an)a=(a^1,\dots,a^n) を持つ解の族をとる。各解を (q(t;a),p(t;a))(q(t;a),p(t;a)) と書く。時刻 tt での点の集まりは、相空間の nn 次元の部分多様体になる。

Lt={(q(t;a), p(t;a))∣a を動かす}L_t=\{(q(t;a),\,p(t;a))\mid a\text{ を動かす}\}

LtL_t が位置の空間の上のグラフになっているとする。すなわち、各位置の上に LtL_t の点がちょうど1つあるとする。このときその点の運動量を読むことで、関数 p=p(q,t)p=p(q,t) が決まる。

これが運動量の場の正体である。場は運動の法則から出るのではなく、どの解の族を考えるかという選択から出る。

例として、1次元の自由粒子 H=p2/2mH=p^2/2m で、時刻0に原点を出る解の族をとる。パラメータは初めの運動量 p0p_0 である。

q(t;p0)=p0mt,p(t;p0)=p0q(t;p_0)=\frac{p_0}{m}t,\qquad p(t;p_0)=p_0

t>0t>0 では p0=mq/tp_0=mq/t と解けるので、この族に沿った運動量の場は次のとおりである。

p(q,t)=mqt(14.21)p(q,t)=\frac{mq}{t}\tag{14.21}

時刻0では L0L_0 は運動量の軸で、位置の上のグラフになっていない。時間が進むと LtL_t は傾き m/tm/t の直線になってグラフになり、時刻が大きくなると位置の軸に近づく(図 14.4)。同じ自由粒子でも、時刻0に運動量をそろえて全位置から出る族を選ぶと、場は定数になる。族の選び方が変われば場も変わる。

図

図 14.4 自由粒子(m=1m=1)で時刻0に原点を出る解の族。左は (q,t)(q,t) 平面の軌道(紺)と、作用 S=mq2/2tS=mq^2/2t の等高線(朱。下から S=2,1,12S=2,1,\tfrac12)。右は相空間で、時刻 tt の族 LtL_t は直線 p=mq/tp=mq/t である。t=0t=0 の L0L_0 は pp 軸で、qq の上のグラフでない。t>0t>0 では LtL_t がグラフになり、運動量の場 p(q,t)p(q,t) が決まる。

Lagrange 部分多様体

場が関数 SS の勾配として書けるかどうかは、族のとり方による。これを判定するのが次の概念である。

定義 14.23(Lagrange 部分多様体)

2n2n 次元のシンプレクティック多様体の部分多様体 LL で、次元が半分のものを考える。包含写像について ιL∗ω=0\iota_L^*\omega=0 となるとき、これを Lagrange 部分多様体と呼ぶ。言い換えると、この部分多様体に接するベクトル U,VU,V については常に ω(U,V)=0\omega(U,V)=0 である。

命題 14.24(グラフが Lagrange である条件)

位置の空間の領域 UU の上のグラフ L={(q,P(q))∣q∈U}L=\{(q,P(q))\mid q\in U\} について、次は同値である。

  1. (a) LL は Lagrange 部分多様体である。
  2. (b) すべての添字の組で ∂Pi/∂qj=∂Pj/∂qi\partial P_i/\partial q^j=\partial P_j/\partial q^i が成り立つ。
  3. (c) (領域が可縮なら)ある関数 SS で Pi=∂S/∂qiP_i=\partial S/\partial q^i と書ける。
証明

包含写像を q↦(q,P(q))q\mapsto(q,P(q)) と書くと、ιL∗θ=Pi dqi\iota_L^*\theta=P_i\,\mathrm{d} q^i である。外微分をとると次のとおりである。

ιL∗ω=d(Pi dqi)=∂Pi∂qjdqj∧dqi=∑i<j(∂Pj∂qi−∂Pi∂qj)dqi∧dqj\iota_L^*\omega=\mathrm{d}(P_i\,\mathrm{d} q^i)=\frac{\partial P_i}{\partial q^j}\mathrm{d} q^j\wedge\mathrm{d} q^i=\sum_{i<j}\Bigl(\frac{\partial P_j}{\partial q^i}-\frac{\partial P_i}{\partial q^j}\Bigr)\mathrm{d} q^i\wedge\mathrm{d} q^j

これが0であることは (b) と同値である。(b) は1形式 Pi dqiP_i\,\mathrm{d} q^i が閉であることを意味し、可縮な領域では Poincaré の補題によって (c) と同値である。

ファイバー {q=q0}\{q=q_0\} は、その上で dqi=0\mathrm{d} q^i=0 なので Lagrange 部分多様体である。この場合はグラフではない。

命題 14.25(流れは Lagrange 部分多様体を保つ)

LL が Lagrange 部分多様体なら、正準変換による像 F(L)F(L) も Lagrange 部分多様体である。とくに Lt=Φt(L0)L_t=\Phi_t(L_0) は、どの時刻でも Lagrange である。

証明

像の包含写像は F∘ιLF\circ\iota_L である。

(F∘ιL)∗ω=ιL∗F∗ω=ιL∗ω=0(F\circ\iota_L)^*\omega=\iota_L^*F^*\omega=\iota_L^*\omega=0

原点から出る族では、L0L_0 はファイバーで Lagrange である。したがって LtL_t は常に Lagrange で、グラフになっている間は p=∂S/∂qp=\partial S/\partial q と書ける。逆に、時刻0に p=(−y,x)p=(-y,x) のような場を初期値に選ぶと、命題 14.24 の (b) が成り立たない。その場はどんな関数の勾配にもならず、Hamilton–Jacobi の理論が扱うのは Lagrange 部分多様体になる族である。

Hamilton–Jacobi 方程式

定理 14.26(Hamilton–Jacobi)

H(q,p,t)H(q,p,t) を Hamiltonian とする。

  1. (i) 時刻0の Lagrange 部分多様体 L0L_0 が、関数 S0S_0 の勾配のグラフであるとする。Lt=Φt(L0)L_t=\Phi_t(L_0) が時刻 [0,T][0,T] の間ずっと可縮な領域の上のグラフであるとする。このとき関数 S(q,t)S(q,t) が存在して、各時刻の族は p=∂S/∂qp=\partial S/\partial q のグラフになる。SS は初期条件 S(q,0)=S0(q)S(q,0)=S_0(q) と、次の Hamilton–Jacobi 方程式を満たす。

∂S∂t+H(q,∂S∂q,t)=0(14.22)\frac{\partial S}{\partial t}+H\Bigl(q,\frac{\partial S}{\partial q},t\Bigr)=0\tag{14.22} 2. (ii) 逆に、関数 S(q,t)S(q,t) が (14.22) を満たすとする。次の微分方程式の解 q(t)q(t) をとる。

q˙i=∂H∂pi(q,∂S∂q,t)\dot q^i=\frac{\partial H}{\partial p_i}\Bigl(q,\frac{\partial S}{\partial q},t\Bigr)

pi(t)=∂iS(q(t),t)p_i(t)=\partial_iS(q(t),t) と置くと、(q(t),p(t))(q(t),p(t)) は Hamilton 方程式の解である。

証明

(ii) から示す。pi(t)=∂iS(q(t),t)p_i(t)=\partial_iS(q(t),t) を時刻で微分すると、連鎖律から

p˙i=∂i∂tS+∂j∂iS q˙j\dot p_i=\partial_i\partial_tS+\partial_j\partial_iS\,\dot q^j

である。(14.22) を qiq^i で偏微分すると、HH の第2引数にも位置が入っているので次のようになる。

∂i∂tS=−∂H∂qi−∂H∂pj ∂i∂jS\partial_i\partial_tS=-\frac{\partial H}{\partial q^i}-\frac{\partial H}{\partial p_j}\,\partial_i\partial_jS

これを代入し、q˙j=∂H/∂pj\dot q^j=\partial H/\partial p_j を使う。2階微分を含む2項が打ち消し合う。

p˙i=−∂H∂qi−∂H∂pj∂i∂jS+∂j∂iS ∂H∂pj=−∂H∂qi\dot p_i=-\frac{\partial H}{\partial q^i}-\frac{\partial H}{\partial p_j}\partial_i\partial_jS+\partial_j\partial_iS\,\frac{\partial H}{\partial p_j}=-\frac{\partial H}{\partial q^i}

位置の方程式は仮定そのものなので、Hamilton 方程式の2式がそろう。

(i) を示す。

命題 14.25 から、LtL_t は Lagrange である。グラフなので、命題 14.24 から各時刻で LtL_t は p=∂S~/∂qp=\partial\tilde S/\partial q と書ける。ここで S~(q,t)\tilde S(q,t) は、各時刻で定数を除いて決まる関数である。場を Pi(q,t)=∂iS~P_i(q,t)=\partial_i\tilde S と書く。

LtL_t の点は Hamilton 方程式に従うので、p˙i=−∂H/∂qi\dot p_i=-\partial H/\partial q^i である。一方、点 (q(t),P(q(t),t))(q(t),P(q(t),t)) の運動量の変化を連鎖律で書くと、次の式になる。

p˙i=∂tPi+∂jPi q˙j=∂tPi+∂iPj ∂H∂pj\dot p_i=\partial_tP_i+\partial_jP_i\,\dot q^j=\partial_tP_i+\partial_iP_j\,\frac{\partial H}{\partial p_j}

2つ目の等号では、位置の方程式と命題 14.24 (b) の対称性を使った。

場を代入した関数を h(q,t)=H(q,P(q,t),t)h(q,t)=H(q,P(q,t),t) と置く。2つの表し方を等しいと置くと、次のように書ける。

∂tPi=−∂H∂qi−∂H∂pj∂iPj=−∂ih\partial_tP_i=-\frac{\partial H}{\partial q^i}-\frac{\partial H}{\partial p_j}\partial_iP_j=-\partial_ih

Pi=∂iS~P_i=\partial_i\tilde S なので、次が成り立つ。

∂i(∂tS~+h)=0\partial_i\bigl(\partial_t\tilde S+h\bigr)=0

したがって ∂tS~+h\partial_t\tilde S+h は、位置によらない関数 c(t)c(t) である。

その時間積分を引いたものを、新しく SS と置く。

S(q,t)=S~(q,t)−∫0tc(s) dsS(q,t)=\tilde S(q,t)-\int_0^tc(s)\,\mathrm{d} s

勾配は変わらず、(14.22) が成り立つ。S(q,0)S(q,0) は定数のずれを合わせて S0S_0 にとる。

この定理が、運動量が場になる理由の答えである。(14.22) の ∂S/∂q\partial S/\partial q は1本の軌道の運動量ではない。Lagrange 部分多様体になる解の族 LtL_t を、位置の上のグラフとして読んだものである。時空の各点を族の解がちょうど1本通るので、その解の運動量がその点に割り当てられる。

命題 14.27(関数 S は作用である)

定理 14.26 の族の解に沿って、SS の変化率は Lagrangian に等しい。

dSdt=piq˙i−H=L\frac{d S}{d t}=p_i\dot q^i-H=L

したがって S(q,t)S(q,t) は、初期値 S0(q0)S_0(q_0) に族の解に沿った作用積分を足したものである。

S(q,t)=S0(q0)+∫0tL dtS(q,t)=S_0(q_0)+\int_0^tL\,\mathrm{d} t

ここで右辺の積分は、q0q_0 から出て点 (q,t)(q,t) に着く族の解に沿ってとる。

証明

連鎖律と (14.22) から、次のとおりである。

dSdt=∂tS+∂iS q˙i=−H+piq˙i\frac{d S}{d t}=\partial_tS+\partial_iS\,\dot q^i=-H+p_i\dot q^i

右辺は、第3章 §3.4 の Legendre 写像 FL\mathbb FL で結ばれた Lagrangian に等しい。

時刻を加えた拡大相空間の1形式 θ−H dt\theta-H\,\mathrm{d} t を考える。dS\mathrm{d} S は、この1形式を族の解が掃く (n+1)(n+1) 次元の部分多様体へ引き戻したものに等しい。SS が1価に決まるのは、この引き戻しが閉形式だからである。

波面・焦線・ほかの分野へのつながり

SS の等高面を作用の波面と呼ぶ。Hamiltonian が H=gijpipj/2m+V(q)H=g^{ij}p_ip_j/2m+V(q) の形なら、速度は q˙i=gijpj/m\dot q^i=g^{ij}p_j/m である。これは、計量で添字を上げた ∂S\partial S の勾配ベクトルに比例する(第5章 §5.2)。したがって族の軌道は、計量の意味で波面に直交する。

図 14.4 の左は (q,t)(q,t) 平面の図なので、軌道と等高線は直交して見えない。各時刻の位置の空間で見ると、軌道は波面に直交している。3次元なら、波面は原点を中心とする球面である。

光学では、屈折率 n(x)n(\bm{x}) の媒質の中の光線の族について、位相が eikonal 方程式 ∣∇S∣2=n2\left\vert \nabla S\right\vert ^2=n^2 を満たす。光線は波面に垂直で、第10章 §10.5 の Fermat の原理の光線と一致する。量子力学では、波動関数を ψ=A eiS/ℏ\psi=A\,\mathrm{e}^{\mathrm{i} S/\hbar} と置き、Planck 定数を0に近づけた主要項をとる。すると位相が (14.22) を満たし(WKB 近似)、古典的な作用の波面は量子力学の波の位相の極限になる。

族の解が交わると、LtL_t は位置の上のグラフでなくなる。1次元の自由粒子で、時刻0に位置 q0q_0 から運動量 −kq0-kq_0 で出る族をとる。各解は q(t)=q0(1−kt/m)q(t)=q_0(1-kt/m) である(k>0k>0)。時刻 m/km/k ですべての解が原点に集まり、その時刻の LtL_t は運動量の軸になる。

このように解が集まる点の集まりを焦線(コースティック)と呼ぶ。焦線では運動量の場も SS も特異になるが、LtL_t 自体は相空間の滑らかな直線のままである。特異になるのは LtL_t から位置の空間への射影であって、相空間の幾何ではない。焦線を越えると族の解は再び1本ずつに分かれ、LtL_t はまたグラフになる。

例題 14.7 原点から出る自由粒子の作用 ★★★

オリジナル1次元の自由粒子 H=p2/2mH=p^2/2m で、時刻0に原点を出る解の族を考える(t>0t>0)。

  1. (1) S(q,t)=mq2/2tS(q,t)=mq^2/2t が (14.22) を満たすことを示せ。その勾配が (14.21) に一致することも示せ。
  2. (2) 原点を出て点 (q,t)(q,t) に着く解に沿った作用積分を計算し、SS に一致することを確かめよ。
  3. (3) 波面 S=S= 一定は、時刻とともにどう動くか。その速さを粒子の速さと比べよ。
解答

(1) 偏微分は次のとおりである。

∂S∂q=mqt,∂S∂t=−mq22t2\frac{\partial S}{\partial q}=\frac{mq}{t},\qquad \frac{\partial S}{\partial t}=-\frac{mq^2}{2t^2}

1つ目は (14.21) に一致する。これを Hamiltonian に入れると H=(mq/t)2/2m=mq2/2t2H=(mq/t)^2/2m=mq^2/2t^2 で、(14.22) が成り立つ。

(2) 原点を出て点 (q,t)(q,t) に着く解は、速さ v=q/tv=q/t の等速運動である。Lagrangian L=mv2/2L=mv^2/2 は一定なので、次のとおりである。

∫0tL dt=12m(qt)2t=mq22t=S(q,t)\int_0^tL\,\mathrm{d} t=\frac12m\Bigl(\frac{q}{t}\Bigr)^2t=\frac{mq^2}{2t}=S(q,t)

命題 14.27 で S0=0S_0=0 とした場合に当たる。

(3) S=cS=c の点は q=±2ct/mq=\pm\sqrt{2ct/m} で、波面は時刻の平方根に比例して広がる。波面の速さを計算する。

dqdt=122cmt=q2t\frac{d q}{d t}=\frac12\sqrt{\frac{2c}{mt}}=\frac{q}{2t}

これは、その位置を通る粒子の速さ q/tq/t の半分である。量子力学の言葉では、位相速度が群速度の半分になることに対応する。

例題 14.8 完全解と包絡線 ★★

オリジナル1次元の自由粒子で、定数 α\alpha を含む関数 Sα(q,t)=αq−α2t/2mS_\alpha(q,t)=\alpha q-\alpha^2t/2m を考える。

  1. (1) どの α\alpha でも SαS_\alpha が (14.22) を満たすことを示し、対応する解の族を答えよ。
  2. (2) ∂Sα/∂α=0\partial S_\alpha/\partial\alpha=0 を解いて得た α\alpha を代入すると、例題 14.7 の作用になることを示せ。
解答

(1) 偏微分は ∂Sα/∂q=α\partial S_\alpha/\partial q=\alpha と ∂Sα/∂t=−α2/2m\partial S_\alpha/\partial t=-\alpha^2/2m である。したがって ∂tSα+(∂qSα)2/2m=0\partial_tS_\alpha+(\partial_qS_\alpha)^2/2m=0 が成り立つ。運動量の場は p=αp=\alpha で一定で、族は「運動量をそろえて全位置から出る解」である。この族の LtL_t は水平な直線 p=αp=\alpha である。

(2) ∂Sα/∂α=q−αt/m=0\partial S_\alpha/\partial\alpha=q-\alpha t/m=0 から α=mq/t\alpha=mq/t である。代入すると次のとおりである。

S=mqtq−12m(mqt)2t=mq2t−mq22t=mq22tS=\frac{mq}{t}q-\frac{1}{2m}\Bigl(\frac{mq}{t}\Bigr)^2t=\frac{mq^2}{t}-\frac{mq^2}{2t}=\frac{mq^2}{2t}

運動量ごとの平面波の族を重ね、各点で停留する運動量を選ぶと、点源から出る波の作用が得られる。完全解から軌道を得る Jacobi の方法は、虎の巻『解析力学から様々な理論分野を組み立てる』第1部第13章で扱う。

確認テスト

オリジナル

  1. 14.5-1 平面の相空間の2つのグラフ p=(y,x)p=(y,x) と p=(−y,x)p=(-y,x) のうち、Lagrange 部分多様体はどちらか。Lagrange である方について、場を勾配に持つ関数を求めよ。
  2. 14.5-2 調和振動子 H=p2/2m+mω02q2/2H=p^2/2m+m\omega_0^2q^2/2 で、時刻0に原点を出る解の族を考える。0<t<π/ω00<t<\pi/\omega_0 での運動量の場を求めよ。次の関数が (14.22) を満たすことを確かめ、t=π/ω0t=\pi/\omega_0 で何が起こるかを答えよ。

S=12mω0q2cot⁡ω0tS=\frac12m\omega_0q^2\cot\omega_0t 3. 14.5-3 一様な重力 H=p2/2m+mgqH=p^2/2m+mgq で、時刻0に原点を出る解の族を考える。運動量の場と作用 S(q,t)S(q,t) を求めよ。作用は、原点から出る解に沿った作用積分にとる。4. 14.5-4 本文の焦点を結ぶ族を考える。自由粒子で、初期値は q(0)=q0q(0)=q_0、p(0)=−kq0p(0)=-kq_0 である。0<t<m/k0<t<m/k での運動量の場を求めよ。S(q,0)=−kq2/2S(q,0)=-kq^2/2 となる作用を求め、(14.22) を満たすことも確かめよ。

木下の勉強ステータス

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

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

← AI生成教科書『基礎理論』の記事一覧へ