ハミルトン–ヤコビ方程式

物理学においてハミルトン–ヤコビ方程式 ( ) とは古典力学の再定式化であり、ニュートンの運動方程式、ラグランジュ力学、ハミルトン力学などの他の定式化と同値である。ハミルトン–ヤコビ方程式は力学系において保存される量を探し出す場合に特に便利であり、それはたとえ力学の問題それ自身が完全には解けない場合にでさえも可能である。
ハミルトン–ヤコビ方程式はまた、粒子の運動が波として表現される唯一の力学の定式化である。この視点から、ハミルトン–ヤコビ方程式は理論物理学の長らくの目標（少なくとも18世紀、ヨハン・ベルヌーイ以来）である、光の伝播と粒子の運動との類似性を見出す試みを達成したと見ることも出来る。力学系から得られる波動方程式は以下に示すとおり、シュレーディンガー方程式と、完全にではないがよく似ている。ハミルトン–ヤコビ方程式はこのような理由で、最も量子力学に近い古典力学の扱いであると考えられている。
数学的な定式化.
ハミルトン–ヤコビ方程式はハミルトンの主関数 () formula_1 に対する、一階の非線形偏微分方程式として以下のように表される。
,\dots,\frac{\partial S}{\partial q_{N}};t\right) + \frac{\partial S}{\partial t}=0.
後の節で示すように、この方程式はハミルトン力学において、formula_2 を古典的なハミルトニアン formula_3 の正準変換の母関数と見なすことにより導かれる。共役な運動量には一般化座標による formula_2 の一階の微分
が相当し、それは以下のように示される。
運動の経路をわずかに変化させた場合の作用の変化は以下により与えられる。
実際に起こる運動の経路はオイラー＝ラグランジュ方程式を満たすことから、formula_5 の積分の項はゼロである。最初の項で formula_6 とし、 formula_7 を簡単に formula_8 と書く。formula_9 を formula_10 と置き換え、最終的に
が得られる。この関係から、座標によるハミルトンの主関数 formula_11 の偏微分は、対応する運動量に等しいことが示された。Q.E.D.
同様に、一般化座標は下記のように、運動量の微分として得られる。式を逆に解いて、系の発展を得ることが出来る。すなわち、一般化座標が時間の関数として得られる。始状態での位置と速度は、formula_2 の積分の中で定数として現れ、それらは全エネルギー、角運動量、などの保存量（運動の積分）に対応する。
他の力学の記述との比較.
ハミルトン–ヤコビ方程式は"単一の" 、 formula_13 個の一般化座標 formula_14 と時間 formula_15 の関数 formula_2 に対する一階の偏微分方程式である。一般化運動量は formula_2 の微分としてしか現れない。顕著な特徴であるが、formula_2 は古典的な作用に等しい。
比較として、ラグランジュ力学での同値なオイラー＝ラグランジュ方程式にも、共役な運動量はやはり現れない。しかし、それは formula_13 個の"系" をなす、一般化座標の時間発展に関する一般には二階の微分方程式である。別の比較として、ハミルトンの正準方程式は同じように
formula_20 個の、一般化座標とそれに共役な formula_21 に対する一階の微分方程式の系である。
ハミルトン–ヤコビ方程式は、ハミルトンの原理の積分を最小化する問題と同値なので、ハミルトン–ヤコビ方程式は他の変分法の問題、あるいはさらに一般的な他の数学や物理学の領域、たとえば力学系、シンプレクティック幾何学、量子カオスの問題などにおいても便利である。例として、ハミルトン–ヤコビ方程式はリーマン多様体において測地線を求めるのに用いられるが、これはリーマン幾何学における重要な変分問題である。
記法.
以下では簡単のため、formula_22 のような太字の変数で formula_13 個の一般化座標を表す。
これらは回転操作でベクトルとしての変換を受ける必要はない。ドット積を、対応する成分の積の和として以下のように定義する。
導出.
第二種の母関数による正準変換 formula_24 は全て、以下のような関係を導く。
 = \mathbf{p}, \qquad
ハミルトン–ヤコビ方程式を導くためには、新しいハミルトニアン formula_25 が恒等的にゼロになるような母関数 formula_26 を取る。するとハミルトニアンの全ての微分はゼロになり、正準方程式は以下のように自明な関係になる。
すなわち、新しい一般化座標と運動量は運動の積分となる。新しい一般化運動量 formula_27 は通常 formula_28 ただし formula_29 と書かれる。
ハミルトン–ヤコビ方程式は変換後のハミルトニアン formula_25 に対する方程式として、
と導かれ、これは
,t\right) + {\partial S \over \partial t} = 0,
と、formula_31 とすれば同値である。
新しい一般化座標 formula_32 も同様に定数であり、formula_33 と書かれる。formula_34 について解けた場合、以下の便利な方程式が得られる。
あるいは明示的に成分で書くと
理想的に、これら formula_13 個の方程式は逆に解いて、元の一般化座標を定数 formula_36とformula_37 の関数として表せ、元の問題を解くことができる。
変数分離.
ハミルトン–ヤコビ方程式は変数分離によって解かれる場合に最も便利であり、その場合には保存量が直接的に求められる。例えば、ハミルトニアンが陽には時間 formula_15 に依っていない場合、formula_15 を分離する事が出来る。そのとき、時間微分 formula_40 は定数（通常 formula_41）となる必要があり、分離された解
を与える。時間に依存しない関数 formula_42 は時にハミルトンの特性関数と呼ばれる。簡約されたハミルトン–ヤコビ方程式は以下のようになる。
 \right) = E
他に変数分離が可能な状況として、ある一般化座標 formula_43 とその微分 formula_44 が一つの関数 formula_45 を通してのみハミルトニアンの中に現れるような場合を考える。
この場合、関数 formula_2 は二つの関数に分離でき、片方は formula_43 だけに依存して、他方は残りの一般化座標に依存する。
この形でハミルトン–ヤコビ方程式を置き換えると、関数 formula_48 は定数（以下 formula_49）となる事が示され、formula_50 に関する一階の常微分方程式 が得られる。
幸運な場合では、関数 formula_2 は formula_13 個の関数 formula_53 に完全に分離され以下のようになる。
この場合、問題は formula_13 個の常微分方程式に帰着する。
formula_2 が変数分離可能かどうかは、ハミルトニアンの形と一般化座標の選び方の両方に依存する。直交座標でハミルトニアンが時間に依存せず、一般化運動量について二次式である場合に、以下の条件を満たせば formula_2 は分離可能である。
すなわち、ポテンシャルエネルギーの項が加法的に各々の座標について分離可能で、各々の座標に対するポテンシャルエネルギーの項がハミルトニアンの対応する運動項と同じ座標依存の因子を掛けられている場合である（ステッケルの条件）。2自由度系（formula_57）の場合、系が直交座標、極座標、放物線座標、楕円座標のいずれかで変数分離可能であるとき、またそのときに限り、運動量について2次の運動の積分が存在し求積可能であることが知られている（ベルトラン・ダルブーの定理）。
直交曲線座標におけるいくつかの例を以下に示す。
球座標の例.
球座標におけるハミルトニアンは以下のように書かれる。
ハミルトン–ヤコビ方程式が完全に分離可能なのは、formula_58 が同じような以下の形式を持つ場合である。
 + \frac{U_{\phi}(\phi)}{r^{2}\sin^{2}\theta}.
ここで formula_59, formula_60, formula_61 は任意の関数とする。完全に分離された解 formula_62 をハミルトン–ヤコビ方程式に代入すると以下が得られる。
\frac{1}{2m r^{2}} \left[ \left( \frac{dS_{\theta}}{d\theta} \right)^{2} + 2m U_{\theta}(\theta) \right] + 
\frac{1}{2m r^{2}\sin^{2}\theta} \left[ \left( \frac{dS_{\phi}}{d\phi} \right)^{2} + 2m U_{\phi}(\phi) \right] = E
この式は常微分方程式の積分によって解け、最初に formula_63 に関する方程式は以下のようになる。
ただし formula_64 は運動の定数で、ハミルトン–ヤコビ方程式の formula_63 依存性は以下のように消去された。
\frac{1}{2m r^{2}} \left[ \left( \frac{dS_{\theta}}{d\theta} \right)^{2} + 2m U_{\theta}(\theta) + \frac{\Gamma_{\phi}}{\sin^{2}\theta} \right] = E
次の常微分方程式は一般化座標 formula_66 を含む。
再び formula_67 は運動の定数で、formula_66 は消去され、最後にハミルトン–ヤコビ方程式は常微分方程式
となり、これを積分すると formula_2 が求まる。
楕円柱座標の例.
楕円柱座標（）のハミルトニアンは以下のように書かれる。
\frac{p_{z}^{2}}{2m} + U(\mu, \nu, z)
ここで楕円の焦点は formula_70 軸上 formula_71 の点にある。ハミルトン–ヤコビ方程式が完全に分離可能なのは、 formula_58 が以下のように同じような形で与えられた場合である。
ただし formula_73, formula_74, formula_75 は任意の関数である。完全に分離された解 formula_76 をハミルトン–ヤコビ方程式に代入することにより以下が得られる。
\frac{1}{2ma^{2} \left( \sinh^{2} \mu + \sin^{2} \nu\right)} \left[ \left( \frac{dS_{\mu}}{d\mu} \right)^{2} + \left( \frac{dS_{\nu}}{d\nu} \right)^{2} + 2m a^{2} U_{\mu}(\mu) + 2m a^{2} U_{\nu}(\nu)\right] = E
最初の常微分方程式、
を分離し、変形して両辺に分母を掛けると以下の簡約されたハミルトン–ヤコビ方程式が得られる。
さらにこれは独立な 2 つの常微分方程式
</math>
に分離でき、これらを解けば formula_2 の完全な解が得られる。
放物線柱座標の例.
放物線柱座標（）におけるハミルトニアンは
\frac{p_{z}^{2}}{2m} + U(\sigma, \tau, z)
ハミルトン–ヤコビ方程式が完全に分離可能なのは、formula_58 が以下のように同じような形で与えられた場合である。
 + U_{z}(z)
formula_80、 formula_81 と formula_75 は任意の関数である。完全に分離された formula_83 をハミルトン–ヤコビ方程式に代入し、
\frac{1}{2m \left( \sigma^{2} + \tau^{2} \right)} \left[ \left( \frac{dS_{\sigma}}{d\sigma} \right)^{2} + \left( \frac{dS_{\tau}}{d\tau} \right)^{2} + 2m U_{\sigma}(\sigma) + 2m U_{\tau}(\tau)\right] = E
最初の常微分方程式
を分離し、変形して両辺に分母を掛けると以下の簡約されたハミルトン–ヤコビ方程式が得られる。
さらにこれは独立な 2 つの常微分方程式
</math>
に分離でき、これらを解けばformula_2の完全な解が得られる。
シュレーディンガー方程式との関係.
関数 formula_86 の、三次元空間上のは、すべての時間 formula_15 について定められる。ある formula_2 の等値面の、時間の関数としての運動は、等値面上のある点 formula_22 から始まる粒子の運動により定義される。そのような等値面の運動は formula_22 空間を運動する"波動"と考えることができるが、その運動は完全に波動方程式に従うわけではない。これを示すため、formula_2 で波の位相を表すようにすると
ここで formula_92 は指数関数の引数を無次元化するために導入した定数である。波の振幅は formula_2 を複素数にすることによって表現する。そうしてハミルトン–ヤコビ方程式を書き直すと
これはシュレーディンガー方程式の"非線形"な変種である。
逆に、シュレーディンガー方程式と formula_48 に関する仮設からスタートすると以下のようになる。
上のシュレーディンガー方程式の古典極限 formula_95 が、以下のようなハミルトン–ヤコビ方程式の変種と等しいことが分かった。
具体例.
以下、formula_96 は計量テンソルの共変な成分であり、formula_97 は粒子の静止質量、formula_98 は光速である。