tecktom.com

当てる道具:作用を終点の関数として見る

2026年10月15日11分で読めますAI生成教科書『虎の巻』第1部 第13章 Hamilton–Jacobi 方程式
目次
  1. 終点を動かしたときの作用の変化
  2. 時間発展の母関数
  3. 真の経路が1本でなくなる所
  4. 木下の層

AI 生成・木下が未勉強

解析力学から様々な理論分野を組み立てる 第1部 第13章 §1.13.2 ★★★ 真の経路に沿った作用を終点の関数とみなすと Hamilton–Jacobi 方程式の解になり、時間発展という正準変換の母関数になる

終点を動かしたときの作用の変化

前の節では Hamilton–Jacobi 方程式を母関数の条件として導いたが、その解 SS が何を表すかはまだ分からない。ここでは第1部第3章・第4章の変分の道具を当てて、SS が作用であることを示す。

時刻 t0t_0 に点 q0q_0 を出て時刻 tt に点 qq に着く真の経路(Euler–Lagrange 方程式の解)が1本だけあるとし、その経路に沿った作用

S(q,t;q0,t0)≔∫t0tL(q(t′),q˙(t′),t′) dt′(q(t′) は真の経路)(1.13.5)S(q,t;q_0,t_0)\coloneqq\int_{t_0}^{t}L\bigl(q(t'),\dot q(t'),t'\bigr)\,\mathrm{d} t'\qquad(q(t')\text{ は真の経路})\tag{1.13.5}

を、終点 (q,t)(q,t) と始点 (q0,t0)(q_0,t_0) の普通の関数とみなす。作用汎関数 S[q]S[q] が「経路を入れて数を返す」のに対して、この SS は「端点を入れて、その端点を結ぶ真の経路の作用を返す」。

終点の位置を qq から q+δqq+\delta q へ動かすと、真の経路も少し変わる。第1部第4章で示したとおり、作用の第1変分は内部の Euler–Lagrange 方程式の項と端の項に分かれ、真の経路では内部の項が消えて端の項だけが残る。始点を固定し終点だけを動かすので

δS=∑ipi δqi∣t,すなわち∂S∂qi=pi(1.13.6)\delta S=\sum_ip_i\,\delta q_i\Big|_{t},\qquad\text{すなわち}\qquad \frac{\partial S}{\partial q_i}=p_i\tag{1.13.6}

を得る。pi=∂L/∂q˙ip_i=\partial L/\partial\dot q_i は終点に着いたときの正準運動量である。

時刻 tt についての微分は次のように出す。終点を真の経路そのものに沿って動かすと、SS は被積分関数の値 LL の割合で増える。一方、そのとき (q,t)(q,t) は dqi=q˙i dt\mathrm{d} q_i=\dot q_i\,\mathrm{d} t だけ動くので、連鎖律から

L=dSdt∣経路に沿って=∂S∂t+∑i∂S∂qiq˙i=∂S∂t+∑ipiq˙i.L=\frac{\mathrm{d} S}{\mathrm{d} t}\Big|_{\text{経路に沿って}}=\frac{\partial S}{\partial t}+\sum_i\frac{\partial S}{\partial q_i}\dot q_i=\frac{\partial S}{\partial t}+\sum_ip_i\dot q_i .

よって ∂S/∂t=L−∑ipiq˙i=−H\partial S/\partial t=L-\sum_ip_i\dot q_i=-H である。これと式 (1.13.6) を合わせると H(q,∂S/∂q,t)+∂S/∂t=0H(q,\partial S/\partial q,t)+\partial S/\partial t=0 となり、真の経路の作用 (1.13.5) は Hamilton–Jacobi 方程式 (1.13.2) を満たす。始点についても同じ計算で、端の項の符号が逆になるので ∂S/∂q0i=−p0i\partial S/\partial q_{0i}=-p_{0i}, ∂S/∂t0=H(q0,p0,t0)\partial S/\partial t_0=H(q_0,p_0,t_0) となる。まとめると

dS=∑ipi dqi−H dt−∑ip0i dq0i+H0 dt0.(1.13.7)\mathrm{d} S=\sum_ip_i\,\mathrm{d} q_i-H\,\mathrm{d} t-\sum_ip_{0i}\,\mathrm{d} q_{0i}+H_0\,\mathrm{d} t_0 .\tag{1.13.7}

時間発展の母関数

式 (1.13.7) を前の章の母関数の言葉で読む。始点 t0t_0 を固定し、始点の位置 q0q_0 を新しい座標 QQ とみなすと、S(q,t;Q)S(q,t;Q) は旧い座標 qq と新しい座標 QQ を引数にとる母関数 F1(q,Q,t)F_1(q,Q,t) の形をしている。F1F_1 の規則 p=∂F1/∂qp=\partial F_1/\partial q, P=−∂F1/∂QP=-\partial F_1/\partial Q に式 (1.13.7) を当てると、新しい運動量は P=−∂S/∂q0=p0P=-\partial S/\partial q_0=p_0 になり、新しい Hamiltonian は K=H+∂S/∂t=0K=H+\partial S/\partial t=0 になる。つまり、真の経路の作用は「時刻 tt の状態 (q,p)(q,p) を時刻 t0t_0 の状態 (q0,p0)(q_0,p_0) へ戻す正準変換」の母関数である。前の章で「時間発展は正準変換」と述べたが、その変換の母関数を具体的に書くと作用になる。始点の位置 q0q_0 を定数 α\alpha とみなせば、S(q,t;q0,t0)S(q,t;q_0,t_0) は1.13.1節の意味での完全解の1つで、Jacobi の定理の β=∂S/∂α=−p0\beta=\partial S/\partial\alpha=-p_0 は始点の運動量(の符号を変えたもの)である。

真の経路が1本でなくなる所

この構成は、端点を結ぶ真の経路が1本だけ、という仮定に頼っている。調和振動子 L=12mq˙2−12mω2q2L=\frac12m\dot q^2-\frac12m\omega^2q^2 では、q0q_0 を出る真の経路は q(t′)=q0cos⁡ωt′+(p0/mω)sin⁡ωt′q(t')=q_0\cos\omega t'+(p_0/m\omega)\sin\omega t' で、ωt=π\omega t=\pi では p0p_0 によらず q(t)=−q0q(t)=-q_0 になる(図 1.13.1)。この時刻には q=−q0q=-q_0 を通る真の経路が無数にあり、それ以外の qq を通る真の経路は1本もないので、S(q,t;q0)S(q,t;q_0) は定義できない。完全解の条件 (1.13.3) で言えば、例題 1.13.2 で扱う主関数では ∂2S/∂q ∂q0=−mω/sin⁡ωt\partial^2S/\partial q\,\partial q_0=-m\omega/\sin\omega t が ωt→π\omega t\to\pi で発散する。同じ始点を出た近くの経路が再び集まるこの点は、第1部第5章で焦点と呼んだもの(図 1.5.2)と同じで、ここでは始点の共役点と呼ぶ。第1部第3章では、端を固定した ∫(12q˙2−12ω2q2)dt\int(\frac12\dot q^2-\frac12\omega^2q^2)\mathrm{d} t の停留点が、区間の長さ TT が π/ω\pi/\omega を超えると最小点でなくなることを2次の変分で見た。その境目の T=π/ωT=\pi/\omega が、ここで経路が再び集まる時刻と一致している。

図

図 1.13.1 調和振動子で、同じ点 q0q_0 を時刻 00 に出る真の経路(初めの運動量 p0/(mωq0)=−1.5,−0.75,0,0.75,1.5p_0/(m\omega q_0)=-1.5,-0.75,0,0.75,1.5。厳密な解なので点線)。すべての経路が ωt=π\omega t=\pi で −q0-q_0 に集まる。この時刻には、(q0,0)(q_0,0) と (q,t)(q,t) を結ぶ真の経路が1本に決まらず、主関数 S(q,t;q0)S(q,t;q_0) が定義できない

例題 1.13.2 調和振動子の主関数が Hamilton–Jacobi 方程式を満たすことを確かめる ★★★

オリジナル調和振動子 L=12mq˙2−12mω2q2L=\frac12m\dot q^2-\frac12m\omega^2q^2 の真の経路の作用は、第1部第5章の式 (1.5.8) でqa→q0q_{\mathrm a}\to q_0, qb→qq_{\mathrm b}\to q, T→tT\to t と読み替えた

S(q,t;q0)=mω2sin⁡ωt[(q2+q02)cos⁡ωt−2qq0](0<ωt<π)S(q,t;q_0)=\frac{m\omega}{2\sin\omega t}\bigl[(q^2+q_0^2)\cos\omega t-2qq_0\bigr]\qquad(0<\omega t<\pi)

である(∂S/∂q=p\partial S/\partial q=p と ∂S/∂t=−E\partial S/\partial t=-E は第5章の例題 1.5.5 で q0=0q_0=0 の場合を確かめた)。

  1. (1) ∂S/∂q\partial S/\partial q と ∂S/∂t\partial S/\partial t を計算し、SS が Hamilton–Jacobi 方程式 (1.13.2) を満たすことを代入で確かめよ。
  2. (2) m=ω=1m=\omega=1, q0=1q_0=1, q=12q=\frac12, t=π/3t=\pi/3 のとき、SS・p=∂S/∂qp=\partial S/\partial q・−∂S/∂t-\partial S/\partial t の値を求め、(1) の式が数値でも成り立つことを確かめよ。
  3. (3) q0=1q_0=1, q=12q=\frac12 を止めて ωt→π\omega t\to\pi とすると、p0=−∂S/∂q0p_0=-\partial S/\partial q_0 と ∂2S/∂q ∂q0\partial^2S/\partial q\,\partial q_0 はどうなるか。その結果を共役点(図 1.13.1)と結びつけて読め。
解答

(1) qq で微分すると

∂S∂q=mω(qcos⁡ωt−q0)sin⁡ωt.\frac{\partial S}{\partial q}=\frac{m\omega(q\cos\omega t-q_0)}{\sin\omega t}.

tt で微分するには S=mω2[(q2+q02)cot⁡ωt−2qq0/sin⁡ωt]S=\frac{m\omega}{2}\bigl[(q^2+q_0^2)\cot\omega t-2qq_0/\sin\omega t\bigr] と書き、(cot⁡x)′=−1/sin⁡2x(\cot x)'=-1/\sin^2x, (1/sin⁡x)′=−cos⁡x/sin⁡2x(1/\sin x)'=-\cos x/\sin^2x を使って

∂S∂t=−mω22sin⁡2ωt[q2+q02−2qq0cos⁡ωt].\frac{\partial S}{\partial t}=-\frac{m\omega^2}{2\sin^2\omega t}\bigl[q^2+q_0^2-2qq_0\cos\omega t\bigr].

一方 12m(∂S∂q)2=mω22sin⁡2ωt[q2cos⁡2ωt−2qq0cos⁡ωt+q02]\frac{1}{2m}\bigl(\frac{\partial S}{\partial q}\bigr)^2=\frac{m\omega^2}{2\sin^2\omega t}\bigl[q^2\cos^2\omega t-2qq_0\cos\omega t+q_0^2\bigr] なので、2つを足すと qq0qq_0 と q02q_0^2 の項が消え、

∂S∂t+12m(∂S∂q)2=mω2q22sin⁡2ωt(cos⁡2ωt−1)=−12mω2q2.\frac{\partial S}{\partial t}+\frac{1}{2m}\Bigl(\frac{\partial S}{\partial q}\Bigr)^2=\frac{m\omega^2q^2}{2\sin^2\omega t}\bigl(\cos^2\omega t-1\bigr)=-\frac12m\omega^2q^2 .

よって ∂S∂t+12m(∂S∂q)2+12mω2q2=0\frac{\partial S}{\partial t}+\frac{1}{2m}\bigl(\frac{\partial S}{\partial q}\bigr)^2+\frac12m\omega^2q^2=0 で、Hamilton–Jacobi 方程式が成り立つ。

(2) cos⁡(π/3)=12\cos(\pi/3)=\frac12, sin⁡(π/3)=32\sin(\pi/3)=\frac{\sqrt3}{2} を入れると

S=13[54⋅12−1]=−38≈−0.2165,p=1/4−13/2=−32,−∂S∂t=12⋅34[14+1−12]=12.\begin{gathered}S=\frac{1}{\sqrt3}\Bigl[\frac54\cdot\frac12-1\Bigr]=-\frac{\sqrt3}{8}\approx-0.2165,\\ p=\frac{1/4-1}{\sqrt3/2}=-\frac{\sqrt3}{2},\qquad -\frac{\partial S}{\partial t}=\frac{1}{2\cdot\frac34}\Bigl[\frac14+1-\frac12\Bigr]=\frac12 .\end{gathered}

12p2+12q2=38+18=12\frac12p^2+\frac12q^2=\frac38+\frac18=\frac12 なので ∂S∂t+H=−12+12=0\frac{\partial S}{\partial t}+H=-\frac12+\frac12=0 となり、(1) の式が数値でも成り立つ。このとき p0=(q−q0cos⁡ωt)/sin⁡ωt=0p_0=(q-q_0\cos\omega t)/\sin\omega t=0 で、経路は q(t′)=cos⁡t′q(t')=\cos t'、エネルギーは 12q02=12\frac12q_0^2=\frac12 である。

(3) p0=−∂S∂q0=mω(q−q0cos⁡ωt)sin⁡ωtp_0=-\frac{\partial S}{\partial q_0}=\dfrac{m\omega(q-q_0\cos\omega t)}{\sin\omega t} で、ωt→π\omega t\to\pi では分子が mω(q+q0)=32mω≠0m\omega(q+q_0)=\frac32m\omega\neq0、分母が +0+0 に近づくので p0→+∞p_0\to+\infty である。また ∂2S/∂q ∂q0=−mω/sin⁡ωt→−∞\partial^2S/\partial q\,\partial q_0=-m\omega/\sin\omega t\to-\infty で、完全解の条件 (1.13.3) の行列式が発散する。ωt=π\omega t=\pi には q0q_0 を出たすべての経路が −q0=−1-q_0=-1 に集まるので、q=12q=\frac12 に着くには初めの運動量を限りなく大きくするしかなく、ωt=π\omega t=\pi ちょうどではそのような経路が無い。主関数が ωt→π\omega t\to\pi で発散するのは、終点 qq と始点 q0q_0 を結ぶ真の経路が1本に決まらなくなる共役点を、式の側から見たものである。

木下の勉強ステータス

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

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

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