tecktom.com

この章で作る道具:新しい Hamiltonian を 00 にする正準変換

2026年10月15日11分で読めますAI生成教科書『虎の巻』第1部 第13章 Hamilton–Jacobi 方程式
目次
  1. 母関数への条件
  2. 完全解と Jacobi の定理
  3. 木下の層

AI 生成・木下が未勉強

解析力学から様々な理論分野を組み立てる 第1部 第13章 §1.13.1 ★★★ K≡0K\equiv0 となる母関数を探す条件が Hamilton–Jacobi 方程式。解が見つかれば、運動は微分だけで取り出せる

前の章(第1部第12章)では、時間発展そのものが正準変換の連続であることから入り、母関数で正準変換を作る方法を見た。この章では、その母関数を1つの目的に使う。新しい変数 (Q,P)(Q,P) での Hamiltonian KK が恒等的に 00 になるように母関数を選べば、新しい変数の Hamilton 方程式は Q˙=P˙=0\dot Q=\dot P=0 となり、新しい変数はすべて定数になる。運動を解く問題は、そのような母関数を見つける問題に置き換わる。母関数が満たすべき条件を書き下すと、1本の偏微分方程式

H(q,∂S∂q,t)+∂S∂t=0H\Bigl(q,\frac{\partial S}{\partial q},t\Bigr)+\frac{\partial S}{\partial t}=0

が出る。これが Hamilton–Jacobi 方程式で、未知関数 S(q,t)S(q,t) を Hamilton の主関数と呼ぶ。

この章で作る道具と出てくるものは次のとおりである。1.13.1節で上の方程式を母関数の条件として導き、その解から運動を取り出す手順(Jacobi の定理)を示す。1.13.2節では同じ SS を別の側から作る。真の経路に沿った作用を終点の関数とみなすと、それが Hamilton–Jacobi 方程式を満たし、しかも時間発展という正準変換の母関数になっている。1.13.3節から1.13.5節が「出てくるもの」で、時間を分離した特性関数 W(q)W(q) による1次元の運動、中心力の変数分離による Kepler の軌道、そして SS を波の位相とみなす見方(作用の波面・特性曲線・Schrödinger 方程式の古典極限)を順に扱う。1.13.6節では、エネルギー保存から時間を積分で出す学部の解き方、WKB 近似の公式が、この章の式のどこにあたるかを並べる。

母関数への条件

自由度 ff の系を一般化座標 q=(q1,…,qf)q=(q_1,\dots,q_f) と正準運動量 p=(p1,…,pf)p=(p_1,\dots,p_f) で書き、Hamiltonian を H(q,p,t)H(q,p,t) とする。前の章で、旧い座標 qq と新しい運動量 PP を引数にとる母関数 F2(q,P,t)F_2(q,P,t) から

pi=∂F2∂qi,Qi=∂F2∂Pi,K=H+∂F2∂t(1.13.1)p_i=\frac{\partial F_2}{\partial q_i},\qquad Q_i=\frac{\partial F_2}{\partial P_i},\qquad K=H+\frac{\partial F_2}{\partial t}\tag{1.13.1}

によって正準変換 (q,p)→(Q,P)(q,p)\to(Q,P) が作れることを見た。新しい変数の時間発展は新しい Hamiltonian KK の Hamilton 方程式 Q˙i=∂K/∂Pi\dot Q_i=\partial K/\partial P_i, P˙i=−∂K/∂Qi\dot P_i=-\partial K/\partial Q_i に従う。

ここで K≡0K\equiv0 となる F2F_2 を選べたとする。右辺がどちらも 00 なので Q˙i=P˙i=0\dot Q_i=\dot P_i=0 で、2f2f 個の新しい変数はすべて定数になる。新しい運動量の値を Pi=αiP_i=\alpha_i、新しい座標の値を Qi=βiQ_i=\beta_i と書くと、式 (1.13.1) の2番目の式 βi=∂F2/∂αi\beta_i=\partial F_2/\partial\alpha_i はqq と tt と 2f2f 個の定数の関係になり、これを qq について解いたものが運動 q(t)q(t) である。K≡0K\equiv0 の条件は、式 (1.13.1) の1番目の式で pp を ∂F2/∂q\partial F_2/\partial q に置き換えると

H(q1,…,qf, ∂S∂q1,…,∂S∂qf, t)+∂S∂t=0(1.13.2)H\Bigl(q_1,\dots,q_f,\ \frac{\partial S}{\partial q_1},\dots,\frac{\partial S}{\partial q_f},\ t\Bigr)+\frac{\partial S}{\partial t}=0\tag{1.13.2}

になる。ここで F2F_2 を S(q,α,t)S(q,\alpha,t) と書き直した。式 (1.13.2) を Hamilton–Jacobi 方程式、その解 SS を Hamilton の主関数と呼ぶ。SS は qq と tt の関数として引数つきで書き、作用汎関数 S[q]S[q](角括弧)と区別する。両者の関係は1.13.2節で明らかになる。

式 (1.13.2) は f+1f+1 個の変数 (q,t)(q,t) についての1階の偏微分方程式で、未知関数 SS の微分について一般に非線型である。2f2f 本の連立常微分方程式(Hamilton 方程式)を、1本の偏微分方程式に置き換えたことになる。偏微分方程式を解くほうが易しいとは限らないが、変数が分離できる系(1.13.3節・1.13.4節)では、積分を並べるだけで運動が求まる。

完全解と Jacobi の定理

偏微分方程式には解がいくらでもあるので、運動を取り出すにはどのような解が要るかを決めておく。

定義 1.13.1(完全解)

Hamilton–Jacobi 方程式 (1.13.2) の解 S(q,α,t)S(q,\alpha,t) が ff 個の定数 α=(α1,…,αf)\alpha=(\alpha_1,\dots,\alpha_f) を含み、

det⁡(∂2S∂qi ∂αj)≠0(1.13.3)\det\Bigl(\frac{\partial^2S}{\partial q_i\,\partial\alpha_j}\Bigr)\neq0\tag{1.13.3}

を満たすとき、SS を完全解と呼ぶ。SS に定数を足しても式 (1.13.2) は満たされるが、その足し算の定数は α\alpha に数えない。

条件 (1.13.3) は、ff 個の定数が互いに独立に SS の qq 依存性に入っていることを表す。これが成り立てば、βi=∂S/∂αi\beta_i=\partial S/\partial\alpha_i を qq について解くことができる(陰関数定理)。

定理 1.13.2(Jacobi の定理 ★★★)

S(q,α,t)S(q,\alpha,t) を完全解とし、β=(β1,…,βf)\beta=(\beta_1,\dots,\beta_f) を定数とする。

∂S∂αi(q,α,t)=βi,pi=∂S∂qi(q,α,t)(i=1,…,f)(1.13.4)\frac{\partial S}{\partial \alpha_i}(q,\alpha,t)=\beta_i,\qquad p_i=\frac{\partial S}{\partial q_i}(q,\alpha,t)\qquad(i=1,\dots,f)\tag{1.13.4}

の1番目の式を qq について解いて得た q(t)q(t) と、それを2番目の式に入れた p(t)p(t) は、Hamilton 方程式を満たす。2f2f 個の定数 (α,β)(\alpha,\beta) を動かすと、Hamilton 方程式の解がすべて得られる。

証明

K≡0K\equiv0 の正準変換として前の章の結果から従うが、ここでは式 (1.13.2) だけを使って直接確かめる。q(t)q(t) に沿って ∂S/∂αi\partial S/\partial\alpha_i が一定なので、その時間微分は

0=ddt∂S∂αi=∂2S∂t ∂αi+∑j∂2S∂qj ∂αiq˙j.0=\frac{\mathrm{d} }{\mathrm{d} t}\frac{\partial S}{\partial \alpha_i}=\frac{\partial^2S}{\partial t\,\partial\alpha_i}+\sum_j\frac{\partial^2S}{\partial q_j\,\partial\alpha_i}\dot q_j .

一方、式 (1.13.2) は α\alpha の値によらず恒等的に成り立つので、αi\alpha_i で微分してよく、連鎖律から

∂2S∂t ∂αi+∑j∂H∂pj∂2S∂qj ∂αi=0\frac{\partial^2S}{\partial t\,\partial\alpha_i}+\sum_j\frac{\partial H}{\partial p_j}\frac{\partial^2S}{\partial q_j\,\partial\alpha_i}=0

を得る(∂H/∂pj\partial H/\partial p_j は p=∂S/∂qp=\partial S/\partial q で評価する)。2式を引くと ∑j(∂2S/∂qj∂αi)(q˙j−∂H/∂pj)=0\sum_j\bigl(\partial^2S/\partial q_j\partial\alpha_i\bigr)\bigl(\dot q_j-\partial H/\partial p_j\bigr)=0 で、条件 (1.13.3) から q˙j=∂H/∂pj\dot q_j=\partial H/\partial p_j が従う。

次に pi(t)=∂S/∂qip_i(t)=\partial S/\partial q_i を時間で微分すると p˙i=∂2S/∂t∂qi+∑j(∂2S/∂qj∂qi)q˙j\dot p_i=\partial^2S/\partial t\partial q_i+\sum_j(\partial^2S/\partial q_j\partial q_i)\dot q_j である。式 (1.13.2) を qiq_i で微分すると、HH は qiq_i に直接と p=∂S/∂qp=\partial S/\partial q を通しての2通りで依存するので

∂2S∂t ∂qi+∂H∂qi+∑j∂H∂pj∂2S∂qj ∂qi=0.\frac{\partial^2S}{\partial t\,\partial q_i}+\frac{\partial H}{\partial q_i}+\sum_j\frac{\partial H}{\partial p_j}\frac{\partial^2S}{\partial q_j\,\partial q_i}=0 .

いま示した q˙j=∂H/∂pj\dot q_j=\partial H/\partial p_j を使って2式を比べると p˙i=−∂H/∂qi\dot p_i=-\partial H/\partial q_i を得る。最後に、任意の初期値 (q(t0),p(t0))(q(t_0),p(t_0)) に対して、p(t0)=∂S/∂qp(t_0)=\partial S/\partial q を満たす α\alpha が選べれば(条件 (1.13.3) がその点の近くでこれを保証する)、β\beta は βi=∂S/∂αi\beta_i=\partial S/\partial\alpha_i で決まり、その初期値を通る解が得られる。

定理の手順は、微分方程式を解く操作を一切含まない。完全解が手に入れば、あとは偏微分と代数方程式だけで運動が決まる。手間はすべて、完全解を見つける段階にある。

自由粒子の完全解

1次元の自由粒子 H=p2/(2m)H=p^2/(2m) の Hamilton–Jacobi 方程式は 12m(∂S/∂q)2+∂S/∂t=0\frac{1}{2m}(\partial S/\partial q)^2+\partial S/\partial t=0 である。S=αq+g(t)S=\alpha q+g(t) と置くと α2/(2m)+g′(t)=0\alpha^2/(2m)+g'(t)=0 なので、S(q,α,t)=αq−α2t/(2m)S(q,\alpha,t)=\alpha q-\alpha^2t/(2m) は完全解になる(∂2S/∂q∂α=1≠0\partial^2S/\partial q\partial\alpha=1\neq0)。Jacobi の定理から β=∂S/∂α=q−αt/m\beta=\partial S/\partial\alpha=q-\alpha t/m で q(t)=β+αt/mq(t)=\beta+\alpha t/m、p=∂S/∂q=αp=\partial S/\partial q=\alpha となる。定数 α\alpha は運動量、β\beta は時刻 00 での位置である。

例題 1.13.1 自由粒子の別の完全解 ★★

オリジナル1次元の自由粒子 H=p2/(2m)H=p^2/(2m) について、aa を定数として S(q,a,t)=m(q−a)22tS(q,a,t)=\dfrac{m(q-a)^2}{2t}(t>0t>0)を考える。

  1. (1) SS が Hamilton–Jacobi 方程式を満たし、完全解の条件 (1.13.3) を満たすことを示せ。
  2. (2) Jacobi の定理で q(t)q(t) と p(t)p(t) を求めよ。
  3. (3) (2) の定数 aa, β\beta は物理的に何を表すか。本文の例の (α,β)(\alpha,\beta) と比べよ。
解答

(1) ∂S/∂q=m(q−a)/t\partial S/\partial q=m(q-a)/t, ∂S/∂t=−m(q−a)2/(2t2)\partial S/\partial t=-m(q-a)^2/(2t^2) なので

12m(∂S∂q)2+∂S∂t=m(q−a)22t2−m(q−a)22t2=0.\frac{1}{2m}\Bigl(\frac{\partial S}{\partial q}\Bigr)^2+\frac{\partial S}{\partial t}=\frac{m(q-a)^2}{2t^2}-\frac{m(q-a)^2}{2t^2}=0 .

∂2S/∂q ∂a=−m/t≠0\partial^2S/\partial q\,\partial a=-m/t\neq0 なので完全解である。

(2) β=∂S/∂a=−m(q−a)/t\beta=\partial S/\partial a=-m(q-a)/t を qq について解くと q(t)=a−βt/mq(t)=a-\beta t/m。p=∂S/∂q=m(q−a)/t=−βp=\partial S/\partial q=m(q-a)/t=-\beta で、運動量は一定である。

(3) t→0t\to0 で q→aq\to a なので、aa は時刻 00 での位置、−β-\beta は運動量である。本文の例では α\alpha が運動量、β\beta が時刻 00 の位置で、役割が入れ替わっている。どちらからも同じ自由粒子の運動が得られるが、完全解の選び方によって、定数のうちどれを新しい運動量 PP と呼び、どれを新しい座標 QQ と呼ぶかが変わる。完全解は1つに決まらず、どの完全解を使っても、Jacobi の定理で得られる運動の全体は同じである。

木下の勉強ステータス

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

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

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