物理学においてハミルトン–ヤコビ方程式(ハミルトン–ヤコビほうていしき、英語: Hamilton–Jacobi equation)とは古典力学の再定式化であり、ニュートンの運動方程式、ラグランジュ力学、ハミルトン力学などの他の定式化と同値である。ハミルトン–ヤコビ方程式は力学系において保存される量を探し出す場合に特に便利であり、それはたとえ力学の問題それ自身が完全には解けない場合にでさえも可能である。
ハミルトン–ヤコビ方程式はまた、粒子の運動が波として表現される唯一の力学の定式化である。この視点から、ハミルトン–ヤコビ方程式は理論物理学の長らくの目標(少なくとも18世紀、ヨハン・ベルヌーイ以来)である、光の伝播と粒子の運動との類似性を見出す試みを達成したと見ることも出来る。力学系から得られる波動方程式は以下に示すとおり、シュレーディンガー方程式と、完全にではないがよく似ている。ハミルトン–ヤコビ方程式はこのような理由で、最も量子力学に近い古典力学の扱いであると考えられている。[1][2]
数学的な定式化
ハミルトン–ヤコビ方程式はハミルトンの主関数 (英: Hamilton's principal function)
に対する、一階の非線形偏微分方程式として以下のように表される。[3]

後の節で示すように、この方程式はハミルトン力学において、
を古典的なハミルトニアン
の正準変換の母関数と見なすことにより導かれる。共役な運動量には一般化座標による
の一階の微分

が相当し、それは以下のように示される。
運動の経路をわずかに変化させた場合の作用の変化は以下により与えられる。
![{\displaystyle \delta S=\sum _{i=1}^{N}\left[{\frac {\partial L}{\partial {\dot {q}}_{k}}}\delta q_{k}\right]_{t_{1}}^{t_{2}}+\sum _{i=1}^{N}\int _{t_{1}}^{t_{2}}\left({\frac {\partial L}{\partial q_{k}}}-{\frac {d}{dt}}{\frac {\partial L}{\partial {\dot {q}}_{k}}}\right)\delta q_{k}\,dt.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/619d464778e4fc6ac5f177e6f5871ebb6b3e21d1)
実際に起こる運動の経路はオイラー=ラグランジュ方程式を満たすことから、
の積分の項はゼロである。最初の項で
とし、
を簡単に
と書く。
を
と置き換え、最終的に
.
が得られる。この関係から、座標によるハミルトンの主関数
の偏微分は、対応する運動量に等しいことが示された。Q.E.D.
同様に、一般化座標は下記のように、運動量の微分として得られる。式を逆に解いて、系の発展を得ることが出来る。すなわち、一般化座標が時間の関数として得られる。始状態での位置と速度は、
の積分の中で定数として現れ、それらは全エネルギー、角運動量、ラプラス–ルンゲ–レンツのベクトル(英語版)などの保存量(運動の積分)に対応する。[4]
記法
以下では簡単のため、
のような太字の変数で
個の一般化座標を表す。

これらは回転操作でベクトルとしての変換を受ける必要はない。ドット積を、対応する成分の積の和として以下のように定義する。

導出
第二種の母関数による正準変換
は全て、以下のような関係を導く。

ハミルトン–ヤコビ方程式を導くためには、新しいハミルトニアン
が恒等的にゼロになるような母関数
を取る。するとハミルトニアンの全ての微分はゼロになり、正準方程式は以下のように自明な関係になる。

すなわち、新しい一般化座標と運動量は運動の積分となる。新しい一般化運動量
は通常
ただし
と書かれる。
ハミルトン–ヤコビ方程式は変換後のハミルトニアン
に対する方程式として、

と導かれ、これは

と、
とすれば同値である。
新しい一般化座標
も同様に定数であり、
と書かれる。
について解けた場合、以下の便利な方程式が得られる。

あるいは明示的に成分で書くと

理想的に、これら
個の方程式は逆に解いて、元の一般化座標を定数
と
の関数として表せ、元の問題を解くことができる。
変数分離
ハミルトン–ヤコビ方程式は変数分離によって解かれる場合に最も便利であり、その場合には保存量が直接的に求められる。例えば、ハミルトニアンが陽には時間
に依っていない場合、
を分離する事が出来る。そのとき、時間微分
は定数(通常
)となる必要があり、分離された解

を与える。時間に依存しない関数
は時にハミルトンの特性関数と呼ばれる。簡約されたハミルトン–ヤコビ方程式は以下のようになる。

他に変数分離が可能な状況として、ある一般化座標
とその微分
が一つの関数
を通してのみハミルトニアンの中に現れるような場合を考える。

この場合、関数
は二つの関数に分離でき、片方は
だけに依存して、他方は残りの一般化座標に依存する。

この形でハミルトン–ヤコビ方程式を置き換えると、関数
は定数(以下
)となる事が示され、
に関する一階の常微分方程式 が得られる。

幸運な場合では、関数
は
個の関数
に完全に分離され以下のようになる。

この場合、問題は
個の常微分方程式に帰着する。
が変数分離可能かどうかは、ハミルトニアンの形と一般化座標の選び方の両方に依存する。直交座標でハミルトニアンが時間に依存せず、一般化運動量について二次式である場合に、以下の条件を満たせば
は分離可能である。
すなわち、ポテンシャルエネルギーの項が加法的に各々の座標について分離可能で、各々の座標に対するポテンシャルエネルギーの項がハミルトニアンの対応する運動項と同じ座標依存の因子を掛けられている場合である(ステッケルの条件)。2自由度系(
)の場合、系が直交座標、極座標、放物線座標、楕円座標のいずれかで変数分離可能であるとき、またそのときに限り、運動量について2次の運動の積分が存在し求積可能であることが知られている(ベルトラン・ダルブーの定理)[5]。
直交曲線座標におけるいくつかの例を以下に示す。
球座標の例
球座標におけるハミルトニアンは以下のように書かれる。
![{\displaystyle H={\frac {1}{2m}}\left[p_{r}^{2}+{\frac {p_{\theta }^{2}}{r^{2}}}+{\frac {p_{\phi }^{2}}{r^{2}\sin ^{2}\theta }}\right]+U(r,\theta ,\phi )}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8135e847f807057271308c42129f27a439ed25c2)
ハミルトン–ヤコビ方程式が完全に分離可能なのは、
が同じような以下の形式を持つ場合である。

ここで
,
,
は任意の関数とする。完全に分離された解
をハミルトン–ヤコビ方程式に代入すると以下が得られる。
![{\displaystyle {\frac {1}{2m}}\left({\frac {dS_{r}}{dr}}\right)^{2}+U_{r}(r)+{\frac {1}{2mr^{2}}}\left[\left({\frac {dS_{\theta }}{d\theta }}\right)^{2}+2mU_{\theta }(\theta )\right]+{\frac {1}{2mr^{2}\sin ^{2}\theta }}\left[\left({\frac {dS_{\phi }}{d\phi }}\right)^{2}+2mU_{\phi }(\phi )\right]=E}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0500b4017ff98b17688bbd5ce0f178b5cfff265c)
この式は常微分方程式の積分によって解け、最初に
に関する方程式は以下のようになる。

ただし
は運動の定数で、ハミルトン–ヤコビ方程式の
依存性は以下のように消去された。
![{\displaystyle {\frac {1}{2m}}\left({\frac {dS_{r}}{dr}}\right)^{2}+U_{r}(r)+{\frac {1}{2mr^{2}}}\left[\left({\frac {dS_{\theta }}{d\theta }}\right)^{2}+2mU_{\theta }(\theta )+{\frac {\Gamma _{\phi }}{\sin ^{2}\theta }}\right]=E}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0239284f86fc2a00497f5806ccbfc375d722385e)
次の常微分方程式は一般化座標
を含む。

再び
は運動の定数で、
は消去され、最後にハミルトン–ヤコビ方程式は常微分方程式

となり、これを積分すると
が求まる。
楕円柱座標の例
楕円柱座標(en:elliptic cylindrical coordinates)のハミルトニアンは以下のように書かれる。

ここで楕円の焦点は
軸上
の点にある。ハミルトン–ヤコビ方程式が完全に分離可能なのは、
が以下のように同じような形で与えられた場合である。

ただし
,
,
は任意の関数である。完全に分離された解
をハミルトン–ヤコビ方程式に代入することにより以下が得られる。
![{\displaystyle {\frac {1}{2m}}\left({\frac {dS_{z}}{dz}}\right)^{2}+U_{z}(z)+{\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}+2ma^{2}U_{\mu }(\mu )+2ma^{2}U_{\nu }(\nu )\right]=E}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ccc55248a1c7e94bb2bf451c2c2f707dd5bc7e35)
最初の常微分方程式、

を分離し、変形して両辺に分母を掛けると以下の簡約されたハミルトン–ヤコビ方程式が得られる。

さらにこれは独立な 2 つの常微分方程式


に分離でき、これらを解けば
の完全な解が得られる。
放物線柱座標の例
放物線柱座標(en:parabolic cylindrical coordinates)におけるハミルトニアンは

ハミルトン–ヤコビ方程式が完全に分離可能なのは、
が以下のように同じような形で与えられた場合である。

、
と
は任意の関数である。完全に分離された
をハミルトン–ヤコビ方程式に代入し、
![{\displaystyle {\frac {1}{2m}}\left({\frac {dS_{z}}{dz}}\right)^{2}+U_{z}(z)+{\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}+2mU_{\sigma }(\sigma )+2mU_{\tau }(\tau )\right]=E}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7eccdb6c70e861d3bbc68d7f82cfe80b61b4b715)
最初の常微分方程式

を分離し、変形して両辺に分母を掛けると以下の簡約されたハミルトン–ヤコビ方程式が得られる。

さらにこれは独立な 2 つの常微分方程式


に分離でき、これらを解けば
の完全な解が得られる。
参考文献
- W. Hamilton (1833). “On a General Method of Expressing the Paths of Light, and of the Planets, by the Coefficients of a Characteristic Function” (English). Dublin University Review: 795–826.
- W.Hamilton (1834). “On the Application to Dynamics of a General Mathematical Method previously Applied to Optics” (English). British Association Report: 513–518.
- H. Goldstein (2002). Classical Mechanics. Addison Wesley. ISBN 0-201-65702-3
- A. Fetter and J. Walecka (2003). Theoretical Mechanics of Particles and Continua. Dover Books. ISBN 0-486-43261-0
- L. D. Landau, L. M. Lifshitz (1975). Mechanics. Elsevier, Amsterdam ... Tokyo
- エリ・デ・ランダウ,イェ・エム・リフシッツ 著、恒藤敏彦,広重徹 訳『場の古典論』(原書第6版)東京図書、1978年。