古典力学の定式化
物理学において、 ウィリアム・ローワン・ハミルトン と カール・グスタフ・ヤコブ・ヤコビ にちなんで名付けられた ハミルトン・ヤコビ方程式は、 古典力学 の別の定式化であり、 ニュートンの運動の法則 、 ラグランジアン力学 、 ハミルトン力学 などの他の定式化と同等です 。
ハミルトン・ヤコビ方程式は、粒子の運動を波として表すことができる力学の定式化である。この意味で、この方程式は、光の伝播と粒子の運動の類似性を見つけるという理論物理学の長年の目標(少なくとも18世紀の ヨハン・ベルヌーイにまで遡る)を達成した。機械システムが従う波動方程式は、以下に説明するように、 シュレーディンガー方程式 に似ているが同一ではない 。このため、ハミルトン・ヤコビ方程式は、 古典力学から 量子力学 への「最も近いアプローチ」であると考えられている 。 [1] [2] この関係の質的形式は、 ハミルトンの光機械類似性 と呼ばれている。
数学において、ハミルトン・ヤコビ方程式は 変分法の 問題の一般化における 極限 幾何学を記述する 必要条件である。これは 動的計画法 の ハミルトン・ヤコビ・ベルマン方程式 の特殊なケースとして理解することができる 。 [3]
概要
ハミルトン・ヤコビ方程式は、一階の 非線形 偏微分方程式である。
−
∂
S
∂
t
=
H
(
q
,
∂
S
∂
q
,
t
)
.
{\displaystyle -{\frac {\partial S}{\partial t}}=H\!\!\left(\mathbf {q} ,{\frac {\partial S}{\partial \mathbf {q} }},t\right).}
座標
q
{\displaystyle \mathbf {q} }
にある粒子のシステムについて 。関数は システムのエネルギーを与える システムの ハミルトニアンです。方程式の解は 作用汎関数 であり 、 [4]古い教科書では ハミルトンの主関数 と呼ばれています。解は、 最小作用の原理 で使用される形式の不定積分によって システムラグランジアンに関連付けることができます 。 [5] : 431
定作用の幾何学的表面はシステムの軌跡に垂直であり、システムダイナミクスの波面のようなビューを作成します。ハミルトン-ヤコビ方程式のこの特性は、古典力学と量子力学を結び付けます。 [6] : 175
H
{\displaystyle H}
S
{\displaystyle S}
L
{\displaystyle \ {\mathcal {L}}\ }
S
=
∫
L
d
t
+
s
o
m
e
c
o
n
s
t
a
n
t
{\displaystyle \ S=\int {\mathcal {L}}\ \operatorname {d} t+~{\mathsf {some\ constant}}~}
表記
太字の変数は 一般化座標 のリストを表します 。
q
{\displaystyle \mathbf {q} }
N
{\displaystyle N}
q
=
(
q
1
,
q
2
,
…
,
q
N
−
1
,
q
N
)
{\displaystyle \mathbf {q} =(q_{1},q_{2},\ldots ,q_{N-1},q_{N})}
変数またはリスト上のドットは時間微分を表します( ニュートン記法を 参照)。たとえば、
q
˙
=
d
q
d
t
.
{\displaystyle {\dot {\mathbf {q} }}={\frac {d\mathbf {q} }{dt}}.}
同じ数の座標を持つ2つのリスト間のドット積 表記 は、対応する成分の積の合計を表す略記法です。たとえば、
p
⋅
q
=
∑
k
=
1
N
p
k
q
k
.
{\displaystyle \mathbf {p} \cdot \mathbf {q} =\sum _{k=1}^{N}p_{k}q_{k}.}
作用関数(別名ハミルトンの主関数)
意味
ヘッセ行列は 逆行列であるとする 。この関係は
、 オイラー・ラグランジュ方程式が 2階常微分方程式の系を形成すること
を示している 。行列を逆行列にすると、 この系は次のように変換される
。
H
L
(
q
,
q
˙
,
t
)
=
{
∂
2
L
/
∂
q
˙
i
∂
q
˙
j
}
i
j
{\textstyle H_{\cal {L}}(\mathbf {q} ,\mathbf {\dot {q}} ,t)=\left\{\partial ^{2}{\cal {L}}/\partial {\dot {q}}^{i}\partial {\dot {q}}^{j}\right\}_{ij}}
d
d
t
∂
L
∂
q
˙
i
=
∑
j
=
1
n
(
∂
2
L
∂
q
˙
i
∂
q
˙
j
q
¨
j
+
∂
2
L
∂
q
˙
i
∂
q
j
q
˙
j
)
+
∂
2
L
∂
q
˙
i
∂
t
,
i
=
1
,
…
,
n
,
{\displaystyle {\frac {d}{dt}}{\frac {\partial {\cal {L}}}{\partial {\dot {q}}^{i}}}=\sum _{j=1}^{n}\left({\frac {\partial ^{2}{\cal {L}}}{\partial {\dot {q}}^{i}\partial {\dot {q}}^{j}}}{\ddot {q}}^{j}+{\frac {\partial ^{2}{\cal {L}}}{\partial {\dot {q}}^{i}\partial {q}^{j}}}{\dot {q}}^{j}\right)+{\frac {\partial ^{2}{\cal {L}}}{\partial {\dot {q}}^{i}\partial t}},\qquad i=1,\ldots ,n,}
n
×
n
{\displaystyle n\times n}
H
L
{\displaystyle H_{\cal {L}}}
q
¨
i
=
F
i
(
q
,
q
˙
,
t
)
,
i
=
1
,
…
,
n
.
{\displaystyle {\ddot {q}}^{i}=F_{i}(\mathbf {q} ,\mathbf {\dot {q}} ,t),\ i=1,\ldots ,n.}
時刻と 配置空間の 点を固定します。存在定理と一意性定理は、 条件とを持つすべての 初期 値問題 に対して、 局所的に一意の解を持つことを保証します。さらに、 異なる初期速度を持つ極値が で交差しないほど 十分に小さい時間間隔があるとします。 後者は、任意の と任意の に対して、 となる 極値は最大でも 1 つしかないことを意味します 。を 作用 関数に 代入する と、 ハミルトンの主関数 (HPF) が得られます。
t
0
{\displaystyle t_{0}}
q
0
∈
M
{\displaystyle \mathbf {q} _{0}\in M}
v
0
,
{\displaystyle \mathbf {v} _{0},}
γ
|
τ
=
t
0
=
q
0
{\displaystyle \gamma |_{\tau =t_{0}}=\mathbf {q} _{0}}
γ
˙
|
τ
=
t
0
=
v
0
{\displaystyle {\dot {\gamma }}|_{\tau =t_{0}}=\mathbf {v} _{0}}
γ
=
γ
(
τ
;
t
0
,
q
0
,
v
0
)
.
{\displaystyle \gamma =\gamma (\tau ;t_{0},\mathbf {q} _{0},\mathbf {v} _{0}).}
(
t
0
,
t
1
)
{\displaystyle (t_{0},t_{1})}
v
0
{\displaystyle \mathbf {v} _{0}}
M
×
(
t
0
,
t
1
)
.
{\displaystyle M\times (t_{0},t_{1}).}
q
∈
M
{\displaystyle \mathbf {q} \in M}
t
∈
(
t
0
,
t
1
)
,
{\displaystyle t\in (t_{0},t_{1}),}
γ
=
γ
(
τ
;
t
,
t
0
,
q
,
q
0
)
{\displaystyle \gamma =\gamma (\tau ;t,t_{0},\mathbf {q} ,\mathbf {q} _{0})}
γ
|
τ
=
t
0
=
q
0
{\displaystyle \gamma |_{\tau =t_{0}}=\mathbf {q} _{0}}
γ
|
τ
=
t
=
q
.
{\displaystyle \gamma |_{\tau =t}=\mathbf {q} .}
γ
=
γ
(
τ
;
t
,
t
0
,
q
,
q
0
)
{\displaystyle \gamma =\gamma (\tau ;t,t_{0},\mathbf {q} ,\mathbf {q} _{0})}
S
(
q
,
t
;
q
0
,
t
0
)
=
def
∫
t
0
t
L
(
γ
(
τ
;
⋅
)
,
γ
˙
(
τ
;
⋅
)
,
τ
)
d
τ
,
{\displaystyle S(\mathbf {q} ,t;\mathbf {q} _{0},t_{0})\ {\stackrel {\text{def}}{=}}\int _{t_{0}}^{t}{\mathcal {L}}(\gamma (\tau ;\cdot ),{\dot {\gamma }}(\tau ;\cdot ),\tau )\,d\tau ,}
どこ
γ
=
γ
(
τ
;
t
,
t
0
,
q
,
q
0
)
,
{\displaystyle \gamma =\gamma (\tau ;t,t_{0},\mathbf {q} ,\mathbf {q} _{0}),}
γ
|
τ
=
t
0
=
q
0
,
{\displaystyle \gamma |_{\tau =t_{0}}=\mathbf {q} _{0},}
γ
|
τ
=
t
=
q
.
{\displaystyle \gamma |_{\tau =t}=\mathbf {q} .}
運動 量 は次のように定義されます。このセクションでは 、HPF がわかると、の依存性がなくなること
を 示します。
p
i
(
q
,
q
˙
,
t
)
=
∂
L
/
∂
q
˙
i
.
{\textstyle p_{i}(\mathbf {q} ,\mathbf {\dot {q}} ,t)=\partial {\cal {L}}/\partial {\dot {q}}^{i}.}
p
i
{\displaystyle p_{i}}
q
˙
{\displaystyle \mathbf {\dot {q}} }
実際、時刻と 配置空間の 点を固定します。すべての時刻 と点について、 ハミルトン の主関数の定義からの(唯一の)極値を とします 。 で の速度 を とします 。すると
t
0
{\displaystyle t_{0}}
q
0
{\displaystyle \mathbf {q} _{0}}
t
{\displaystyle t}
q
,
{\displaystyle \mathbf {q} ,}
γ
=
γ
(
τ
;
t
,
t
0
,
q
,
q
0
)
{\displaystyle \gamma =\gamma (\tau ;t,t_{0},\mathbf {q} ,\mathbf {q} _{0})}
S
{\displaystyle S}
v
=
def
γ
˙
(
τ
;
t
,
t
0
,
q
,
q
0
)
|
τ
=
t
{\displaystyle \mathbf {v} \,{\stackrel {\text{def}}{=}}\,{\dot {\gamma }}(\tau ;t,t_{0},\mathbf {q} ,\mathbf {q} _{0})|_{\tau =t}}
τ
=
t
{\displaystyle \tau =t}
∂
S
∂
q
i
=
∂
L
∂
q
˙
i
|
q
˙
=
v
,
i
=
1
,
…
,
n
.
{\displaystyle {\frac {\partial S}{\partial q^{i}}}=\left.{\frac {\partial {\cal {L}}}{\partial {\dot {q}}^{i}}}\right|_{\mathbf {\dot {q}} =\mathbf {v} }\!\!\!\!\!\!\!,\quad i=1,\ldots ,n.}
証拠
以下の証明では、配置空間が の開集合であると仮定していますが、基礎となる手法は任意の 空間 にも同様に適用されます 。この証明の文脈では、カリグラフィ文字は 作用関数を示し、イタリック体は ハミルトンの主関数を示します。
R
n
,
{\displaystyle \mathbb {R} ^{n},}
S
{\displaystyle {\cal {S}}}
S
{\displaystyle S}
ステップ 1。 を配置空間内のパスとし、 に沿ったベクトル場と します 。(各 ベクトルは、 摂動 、 点 における機械システムの 微小変化 または 仮想変位 と呼ばれます )。 点 における 作用の 方向の 変化は 、式で与えられること
を思い出してください。ここで、右側の偏導関数を計算した後 、 と
を代入する必要があります 。(この式は、部分積分によるガトー微分の定義に従います)。
ξ
=
ξ
(
t
)
{\displaystyle \xi =\xi (t)}
δ
ξ
=
δ
ξ
(
t
)
{\displaystyle \delta \xi =\delta \xi (t)}
ξ
{\displaystyle \xi }
t
,
{\displaystyle t,}
δ
ξ
(
t
)
{\displaystyle \delta \xi (t)}
ξ
(
t
)
{\displaystyle \xi (t)}
δ
S
δ
ξ
[
γ
,
t
1
,
t
0
]
{\displaystyle \delta {\cal {S}}_{\delta \xi }[\gamma ,t_{1},t_{0}]}
S
{\displaystyle {\cal {S}}}
ξ
{\displaystyle \xi }
δ
ξ
{\displaystyle \delta \xi }
δ
S
δ
ξ
[
ξ
,
t
1
,
t
0
]
=
∫
t
0
t
1
(
∂
L
∂
q
−
d
d
t
∂
L
∂
q
˙
)
δ
ξ
d
t
+
∂
L
∂
q
˙
δ
ξ
|
t
0
t
1
,
{\displaystyle \delta {\cal {S}}_{\delta \xi }[\xi ,t_{1},t_{0}]=\int _{t_{0}}^{t_{1}}\left({\frac {\partial {\cal {L}}}{\partial \mathbf {q} }}-{\frac {d}{dt}}{\frac {\partial {\cal {L}}}{\partial \mathbf {\dot {q}} }}\right)\delta \xi \,dt+{\frac {\partial {\cal {L}}}{\partial \mathbf {\dot {q}} }}\,\delta \xi {\Biggl |}_{t_{0}}^{t_{1}},}
q
i
=
ξ
i
(
t
)
{\displaystyle q^{i}=\xi ^{i}(t)}
q
˙
i
=
ξ
˙
i
(
t
)
{\displaystyle {\dot {q}}^{i}={\dot {\xi }}^{i}(t)}
が極値であると 仮定します。 はオイラー・ラグランジュ方程式を満たすので、積分項は消えます。 の開始点 が固定されている場合、オイラー・ラグランジュ方程式を導出するのに使用されたのと同じ論理により、 したがって、
ξ
{\displaystyle \xi }
ξ
{\displaystyle \xi }
ξ
{\displaystyle \xi }
q
0
{\displaystyle \mathbf {q} _{0}}
δ
ξ
(
t
0
)
=
0.
{\displaystyle \delta \xi (t_{0})=0.}
δ
S
δ
ξ
[
ξ
,
t
;
t
0
]
=
∂
L
∂
q
˙
|
q
˙
=
ξ
˙
(
t
)
q
=
ξ
(
t
)
δ
ξ
(
t
)
.
{\displaystyle \delta {\cal {S}}_{\delta \xi }[\xi ,t;t_{0}]=\left.{\frac {\partial {\cal {L}}}{\partial \mathbf {\dot {q}} }}\right|_{\mathbf {\dot {q}} ={\dot {\xi }}(t)}^{\mathbf {q} =\xi (t)}\,\delta \xi (t).}
ステップ2。 をHPFの定義からの(唯一の)極値、 に沿ったベクトル場 、および と 「互換性のある」 の変形と します 。正確には、
γ
=
γ
(
τ
;
q
,
q
0
,
t
,
t
0
)
{\displaystyle \gamma =\gamma (\tau ;\mathbf {q} ,\mathbf {q} _{0},t,t_{0})}
δ
γ
=
δ
γ
(
τ
)
{\displaystyle \delta \gamma =\delta \gamma (\tau )}
γ
,
{\displaystyle \gamma ,}
γ
ε
=
γ
ε
(
τ
;
q
ε
,
q
0
,
t
,
t
0
)
{\displaystyle \gamma _{\varepsilon }=\gamma _{\varepsilon }(\tau ;\mathbf {q} _{\varepsilon },\mathbf {q} _{0},t,t_{0})}
γ
{\displaystyle \gamma }
δ
γ
.
{\displaystyle \delta \gamma .}
γ
ε
|
ε
=
0
=
γ
,
{\displaystyle \gamma _{\varepsilon }|_{\varepsilon =0}=\gamma ,}
γ
˙
ε
|
ε
=
0
=
δ
γ
,
{\displaystyle {\dot {\gamma }}_{\varepsilon }|_{\varepsilon =0}=\delta \gamma ,}
γ
ε
|
τ
=
t
0
=
γ
|
τ
=
t
0
=
q
0
.
{\displaystyle \gamma _{\varepsilon }|_{\tau =t_{0}}=\gamma |_{\tau =t_{0}}=\mathbf {q} _{0}.}
HPFとガトー微分の定義により、
δ
S
δ
γ
[
γ
,
t
]
=
def
d
S
[
γ
ε
,
t
]
d
ε
|
ε
=
0
=
d
S
(
γ
ε
(
t
)
,
t
)
d
ε
|
ε
=
0
=
∂
S
∂
q
δ
γ
(
t
)
.
{\displaystyle \delta {\cal {S}}_{\delta \gamma }[\gamma ,t]{\overset {\text{def}}{{}={}}}\left.{\frac {d{\cal {S}}[\gamma _{\varepsilon },t]}{d\varepsilon }}\right|_{\varepsilon =0}=\left.{\frac {dS(\gamma _{\varepsilon }(t),t)}{d\varepsilon }}\right|_{\varepsilon =0}={\frac {\partial S}{\mathbf {\partial q} }}\,\delta \gamma (t).}
ここでは、それを考慮してコンパクトさを 重視しました 。
q
=
γ
(
t
;
q
,
q
0
,
t
,
t
0
)
{\displaystyle \mathbf {q} =\gamma (t;\mathbf {q} ,\mathbf {q} _{0},t,t_{0})}
t
0
{\displaystyle t_{0}}
ステップ 3. ここで、ステップ 1 の 式に とを代入し 、その結果をステップ 2 で導出した式と比較します。 ベクトル場に対して が任意に選択されたという事実により、証明が完了します。
ξ
=
γ
{\displaystyle \xi =\gamma }
δ
ξ
=
δ
γ
{\displaystyle \delta \xi =\delta \gamma }
δ
S
δ
ξ
[
ξ
,
t
;
t
0
]
{\displaystyle \delta {\cal {S}}_{\delta \xi }[\xi ,t;t_{0}]}
t
>
t
0
,
{\displaystyle t>t_{0},}
δ
γ
{\displaystyle \delta \gamma }
機械システムのハミルトニアン が与えられた場合 、ハミルトン・ヤコビ方程式は ハミルトンの主関数に対する一階 非線形 偏微分方程式 である。 [7]
H
(
q
,
p
,
t
)
{\displaystyle H(\mathbf {q} ,\mathbf {p} ,t)}
S
{\displaystyle S}
−
∂
S
∂
t
=
H
(
q
,
∂
S
∂
q
,
t
)
.
{\displaystyle -{\frac {\partial S}{\partial t}}=H\left(\mathbf {q} ,{\frac {\partial S}{\partial \mathbf {q} }},t\right).}
導出
初期速度が 極値の場合 (HPFの定義の前の議論を参照)、
ξ
=
ξ
(
t
;
t
0
,
q
0
,
v
0
)
,
{\displaystyle \xi =\xi (t;t_{0},\mathbf {q} _{0},\mathbf {v} _{0}),}
v
0
=
ξ
˙
|
t
=
t
0
{\displaystyle \mathbf {v} _{0}={\dot {\xi }}|_{t=t_{0}}}
L
(
ξ
(
t
)
,
ξ
˙
(
t
)
,
t
)
=
d
S
(
ξ
(
t
)
,
t
)
d
t
=
[
∂
S
∂
q
q
˙
+
∂
S
∂
t
]
q
˙
=
ξ
˙
(
t
)
q
=
ξ
(
t
)
.
{\displaystyle {\cal {L}}(\xi (t),{\dot {\xi }}(t),t)={\frac {dS(\xi (t),t)}{dt}}=\left[{\frac {\partial S}{\partial \mathbf {q} }}\mathbf {\dot {q}} +{\frac {\partial S}{\partial t}}\right]_{\mathbf {\dot {q}} ={\dot {\xi }}(t)}^{\mathbf {q} =\xi (t)}.}
の公式と、
を満たす ハミルトニアン座標ベースの定義
から、 (一意に解ける) 方程式 が得られる。
ここで 、
p
i
=
p
i
(
q
,
t
)
{\displaystyle p_{i}=p_{i}(\mathbf {q} ,t)}
H
(
q
,
p
,
t
)
=
p
q
˙
−
L
(
q
,
q
˙
,
t
)
,
{\displaystyle H(\mathbf {q} ,\mathbf {p} ,t)=\mathbf {p} \mathbf {\dot {q}} -{\cal {L}}(\mathbf {q} ,\mathbf {\dot {q}} ,t),}
q
˙
(
p
,
q
,
t
)
{\displaystyle \mathbf {\dot {q}} (\mathbf {p} ,\mathbf {q} ,t)}
q
˙
)
{\displaystyle \mathbf {\dot {q}} )}
p
=
∂
L
(
q
,
q
˙
,
t
)
∂
q
˙
,
{\textstyle \mathbf {p} ={\frac {\partial {\cal {L}}(\mathbf {q} ,\mathbf {\dot {q}} ,t)}{\partial \mathbf {\dot {q}} }},}
∂
S
∂
t
=
L
(
q
,
q
˙
,
t
)
−
∂
S
∂
q
q
˙
=
−
H
(
q
,
∂
S
∂
q
,
t
)
,
{\displaystyle {\frac {\partial S}{\partial t}}={\cal {L}}(\mathbf {q} ,\mathbf {\dot {q}} ,t)-{\frac {\partial S}{\mathbf {\partial q} }}\mathbf {\dot {q}} =-H\left(\mathbf {q} ,{\frac {\partial S}{\partial \mathbf {q} }},t\right),}
q
=
ξ
(
t
)
{\displaystyle \mathbf {q} =\xi (t)}
q
˙
=
ξ
˙
(
t
)
.
{\displaystyle \mathbf {\dot {q}} ={\dot {\xi }}(t).}
あるいは、以下に述べるように、ハミルトン力学 からハミルトン・ヤコビ方程式を導くこともできる。これは、 古典ハミルトン方程式の
標準変換 の 生成関数 として 扱うことによって行われる。
S
{\displaystyle S}
H
=
H
(
q
1
,
q
2
,
…
,
q
N
;
p
1
,
p
2
,
…
,
p
N
;
t
)
.
{\displaystyle H=H(q_{1},q_{2},\ldots ,q_{N};p_{1},p_{2},\ldots ,p_{N};t).}
共役運動量は 一般化座標に関する
の一次導関数に対応する。
S
{\displaystyle S}
p
k
=
∂
S
∂
q
k
.
{\displaystyle p_{k}={\frac {\partial S}{\partial q_{k}}}.}
ハミルトン・ヤコビ方程式の解として、主関数には 未定定数が含まれます。最初 の未定定数は と表され 、最後の未定定数は の積分から得られます 。
N
+
1
{\displaystyle N+1}
N
{\displaystyle N}
α
1
,
α
2
,
…
,
α
N
{\displaystyle \alpha _{1},\,\alpha _{2},\dots ,\alpha _{N}}
∂
S
∂
t
{\displaystyle {\frac {\partial S}{\partial t}}}
と の関係は 、これら の運動定数 を用いて 位相空間 における軌道を記述します 。さらに、これらの量
も運動定数であり、これらの方程式を逆にすると、すべての および定数と時間 の関数として見つけることができます 。 [8]
p
{\displaystyle \mathbf {p} }
q
{\displaystyle \mathbf {q} }
β
k
=
∂
S
∂
α
k
,
k
=
1
,
2
,
…
,
N
{\displaystyle \beta _{k}={\frac {\partial S}{\partial \alpha _{k}}},\quad k=1,2,\ldots ,N}
q
{\displaystyle \mathbf {q} }
α
{\displaystyle \alpha }
β
{\displaystyle \beta }
ハミルトン・ヤコビ方程式は、 一般化座標 と時間 の関数に関する 単一 の 1 階偏微分方程式です。一般化運動量は、 古典作用 の導関数としてのみ現れます 。
N
{\displaystyle N}
q
1
,
q
2
,
…
,
q
N
{\displaystyle q_{1},\,q_{2},\dots ,q_{N}}
t
{\displaystyle t}
S
{\displaystyle S}
比較すると、 ラグランジュ力学 の同等の オイラー・ラグランジュの運動方程式 にも共役運動量は現れませんが、これらの方程式は の システム であり、一般に一般化座標の時間発展に関する2次方程式です。同様に、 ハミルトンの運動方程式 は、一般化座標とその共役運動量 の時間発展に関する2 N の1次方程式の別の システム です 。
N
{\displaystyle N}
p
1
,
p
2
,
…
,
p
N
{\displaystyle p_{1},\,p_{2},\dots ,p_{N}}
HJE はハミルトンの原理 のような積分最小化問題と等価な表現であるため 、 変分法の 他の問題や、より一般的には 力学系 、 シンプレクティック幾何学 、 量子カオス などの 数学 や 物理学 の他の分野でも有用である。例えば、ハミルトン・ヤコビ方程式は、 リーマン多様体 上の 測地線を決定するために使用できる。これは リーマン幾何学 における 重要な 変分問題 である。しかし、計算ツールとしては、偏微分方程式を解くのは非常に複雑であることが知られているが、独立変数を分離できる場合は別である。この場合、HJE は計算上有用となる。 [5] : 444
タイプ 2 生成関数 を含む任意の標準 変換は 関係式につながり
、新しい変数 と新しいハミルトニアン
に関するハミルトン方程式 は同じ形式になります。
G
2
(
q
,
P
,
t
)
{\displaystyle G_{2}(\mathbf {q} ,\mathbf {P} ,t)}
p
=
∂
G
2
∂
q
,
Q
=
∂
G
2
∂
P
,
K
(
Q
,
P
,
t
)
=
H
(
q
,
p
,
t
)
+
∂
G
2
∂
t
{\displaystyle \mathbf {p} ={\partial G_{2} \over \partial \mathbf {q} },\quad \mathbf {Q} ={\partial G_{2} \over \partial \mathbf {P} },\quad K(\mathbf {Q} ,\mathbf {P} ,t)=H(\mathbf {q} ,\mathbf {p} ,t)+{\partial G_{2} \over \partial t}}
P
,
Q
{\displaystyle \mathbf {P} ,\,\mathbf {Q} }
K
{\displaystyle K}
P
˙
=
−
∂
K
∂
Q
,
Q
˙
=
+
∂
K
∂
P
.
{\displaystyle {\dot {\mathbf {P} }}=-{\partial K \over \partial \mathbf {Q} },\quad {\dot {\mathbf {Q} }}=+{\partial K \over \partial \mathbf {P} }.}
HJE を導出するには、 新しいハミルトニアン を作成するような生成関数を選択します 。したがって、その導関数もすべてゼロになり、変換されたハミルトン方程式は自明になる
ため、新しい一般化座標 と運動量は 運動 定数 になります。これらは定数であるため、このコンテキストでは、新しい一般化運動量 は通常 、つまり と表され 、新しい 一般化座標 は通常 と表されるため です 。
G
2
(
q
,
P
,
t
)
{\displaystyle G_{2}(\mathbf {q} ,\mathbf {P} ,t)}
K
=
0
{\displaystyle K=0}
P
˙
=
Q
˙
=
0
{\displaystyle {\dot {\mathbf {P} }}={\dot {\mathbf {Q} }}=0}
P
{\displaystyle \mathbf {P} }
α
1
,
α
2
,
…
,
α
N
{\displaystyle \alpha _{1},\,\alpha _{2},\dots ,\alpha _{N}}
P
m
=
α
m
{\displaystyle P_{m}=\alpha _{m}}
Q
{\displaystyle \mathbf {Q} }
β
1
,
β
2
,
…
,
β
N
{\displaystyle \beta _{1},\,\beta _{2},\dots ,\beta _{N}}
Q
m
=
β
m
{\displaystyle Q_{m}=\beta _{m}}
生成関数をハミルトンの主関数に任意の定数を加えた値に設定すると 、
HJEが自動的に発生する。
A
{\displaystyle A}
G
2
(
q
,
α
,
t
)
=
S
(
q
,
t
)
+
A
,
{\displaystyle G_{2}(\mathbf {q} ,{\boldsymbol {\alpha }},t)=S(\mathbf {q} ,t)+A,}
p
=
∂
G
2
∂
q
=
∂
S
∂
q
→
H
(
q
,
p
,
t
)
+
∂
G
2
∂
t
=
0
→
H
(
q
,
∂
S
∂
q
,
t
)
+
∂
S
∂
t
=
0.
{\displaystyle \mathbf {p} ={\frac {\partial G_{2}}{\partial \mathbf {q} }}={\frac {\partial S}{\partial \mathbf {q} }}\,\rightarrow \,H(\mathbf {q} ,\mathbf {p} ,t)+{\partial G_{2} \over \partial t}=0\,\rightarrow \,H\left(\mathbf {q} ,{\frac {\partial S}{\partial \mathbf {q} }},t\right)+{\partial S \over \partial t}=0.}
これらをについて解くと、
わかりやすいように成分で書かれた
便利な方程式や が得られる。
S
(
q
,
α
,
t
)
{\displaystyle S(\mathbf {q} ,{\boldsymbol {\alpha }},t)}
Q
=
β
=
∂
S
∂
α
,
{\displaystyle \mathbf {Q} ={\boldsymbol {\beta }}={\partial S \over \partial {\boldsymbol {\alpha }}},}
Q
m
=
β
m
=
∂
S
(
q
,
α
,
t
)
∂
α
m
.
{\displaystyle Q_{m}=\beta _{m}={\frac {\partial S(\mathbf {q} ,{\boldsymbol {\alpha }},t)}{\partial \alpha _{m}}}.}
理想的には、これらの N 方程式を逆転させて、 定数およびの関数として元の 一般化座標 を見つけ、 元の問題を解決できます。
q
{\displaystyle \mathbf {q} }
α
,
β
,
{\displaystyle {\boldsymbol {\alpha }},\,{\boldsymbol {\beta }},}
t
{\displaystyle t}
変数の分離
問題が 変数の加法的分離を 許す場合、HJEは直接 運動定数 を導く。例えば、 ハミルトニアンが時間に明示的に依存しない場合、時間 tは 分離できる。その場合、 HJEの時間微分は定数でなければならず、通常は( )で表され、分離された解を与える。
ここで、
時間に依存しない関数は、 省略された作用 または ハミルトンの特性関数 [5] :434 と呼ばれることがあり、 [9] :607 と 表記される こともある ( 作用原理名を 参照)。簡約されたハミルトン・ヤコビ方程式は次のように書ける。
∂
S
∂
t
{\displaystyle {\frac {\partial S}{\partial t}}}
−
E
{\displaystyle -E}
S
=
W
(
q
1
,
q
2
,
…
,
q
N
)
−
E
t
{\displaystyle S=W(q_{1},q_{2},\ldots ,q_{N})-Et}
W
(
q
)
{\displaystyle W(\mathbf {q} )}
S
0
{\displaystyle S_{0}}
H
(
q
,
∂
S
∂
q
)
=
E
.
{\displaystyle H\left(\mathbf {q} ,{\frac {\partial S}{\partial \mathbf {q} }}\right)=E.}
他の変数の分離可能性を説明するために、ある 一般化座標 とその導関数が
ハミルトニアンにおいて
単一の関数として一緒に現れると仮定する。
q
k
{\displaystyle q_{k}}
∂
S
∂
q
k
{\displaystyle {\frac {\partial S}{\partial q_{k}}}}
ψ
(
q
k
,
∂
S
∂
q
k
)
{\displaystyle \psi \left(q_{k},{\frac {\partial S}{\partial q_{k}}}\right)}
H
=
H
(
q
1
,
q
2
,
…
,
q
k
−
1
,
q
k
+
1
,
…
,
q
N
;
p
1
,
p
2
,
…
,
p
k
−
1
,
p
k
+
1
,
…
,
p
N
;
ψ
;
t
)
.
{\displaystyle H=H(q_{1},q_{2},\ldots ,q_{k-1},q_{k+1},\ldots ,q_{N};p_{1},p_{2},\ldots ,p_{k-1},p_{k+1},\ldots ,p_{N};\psi ;t).}
その場合、関数 Sは2つの関数に分割することができ、1つは q k のみに依存し、もう1つは残りの 一般化座標 のみに依存する。
S
=
S
k
(
q
k
)
+
S
rem
(
q
1
,
…
,
q
k
−
1
,
q
k
+
1
,
…
,
q
N
,
t
)
.
{\displaystyle S=S_{k}(q_{k})+S_{\text{rem}}(q_{1},\ldots ,q_{k-1},q_{k+1},\ldots ,q_{N},t).}
これらの式をハミルトン・ヤコビ方程式に代入すると、関数 ψは 定数(ここでは と表記 )でなければならないことがわかり、次の1階常 微分方程式 が得られる。
Γ
k
{\displaystyle \Gamma _{k}}
S
k
(
q
k
)
,
{\displaystyle S_{k}(q_{k}),}
ψ
(
q
k
,
d
S
k
d
q
k
)
=
Γ
k
.
{\displaystyle \psi \left(q_{k},{\frac {dS_{k}}{dq_{k}}}\right)=\Gamma _{k}.}
幸運なケースでは、関数は完全に 関数 に分離できる。
S
{\displaystyle S}
N
{\displaystyle N}
S
m
(
q
m
)
,
{\displaystyle S_{m}(q_{m}),}
S
=
S
1
(
q
1
)
+
S
2
(
q
2
)
+
⋯
+
S
N
(
q
N
)
−
E
t
.
{\displaystyle S=S_{1}(q_{1})+S_{2}(q_{2})+\cdots +S_{N}(q_{N})-Et.}
このような場合、問題は 常微分方程式 に移行します。
N
{\displaystyle N}
S の分離可能性は、ハミルトニアンと 一般化座標 の選択の両方に依存します 。 直交座標 と、時間依存性がなく 一般化運動量に関して 2 次 であるハミルトニアンの場合、各座標の位置エネルギー項にハミルトニアンの対応する運動量項の座標依存係数を乗じると、各座標の位置エネルギー項が加法的に分離可能であれば、S は完全に分離可能になります ( シュテッケル条件)。説明のために、次のセクションでは 直交座標 の例をいくつか 示します。
S
{\displaystyle S}
さまざまな座標系の例
球座標
球座標 では、保存ポテンシャル U 内を運動する自由粒子のハミルトニアンは 次のように表される。
H
=
1
2
m
[
p
r
2
+
p
θ
2
r
2
+
p
ϕ
2
r
2
sin
2
θ
]
+
U
(
r
,
θ
,
ϕ
)
.
{\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 ).}
ハミルトン・ヤコビ方程式は、 類似の形式で記述できる
関数が存在する 限り、これらの座標で完全に分離可能である。
U
r
(
r
)
,
U
θ
(
θ
)
,
U
ϕ
(
ϕ
)
{\displaystyle U_{r}(r),U_{\theta }(\theta ),U_{\phi }(\phi )}
U
{\displaystyle U}
U
(
r
,
θ
,
ϕ
)
=
U
r
(
r
)
+
U
θ
(
θ
)
r
2
+
U
ϕ
(
ϕ
)
r
2
sin
2
θ
.
{\displaystyle U(r,\theta ,\phi )=U_{r}(r)+{\frac {U_{\theta }(\theta )}{r^{2}}}+{\frac {U_{\phi }(\phi )}{r^{2}\sin ^{2}\theta }}.}
完全に分離された溶液を
HJEに代入すると、
S
=
S
r
(
r
)
+
S
θ
(
θ
)
+
S
ϕ
(
ϕ
)
−
E
t
{\displaystyle S=S_{r}(r)+S_{\theta }(\theta )+S_{\phi }(\phi )-Et}
1
2
m
(
d
S
r
d
r
)
2
+
U
r
(
r
)
+
1
2
m
r
2
[
(
d
S
θ
d
θ
)
2
+
2
m
U
θ
(
θ
)
]
+
1
2
m
r
2
sin
2
θ
[
(
d
S
ϕ
d
ϕ
)
2
+
2
m
U
ϕ
(
ϕ
)
]
=
E
.
{\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.}
この方程式は、ハミルトン・ヤコビ方程式から依存性
を排除する 運動定数 である
、 の方程式から始めて、 常微分方程式 を連続的に積分することによって解くことができる。
ϕ
{\displaystyle \phi }
(
d
S
ϕ
d
ϕ
)
2
+
2
m
U
ϕ
(
ϕ
)
=
Γ
ϕ
{\displaystyle \left({\frac {dS_{\phi }}{d\phi }}\right)^{2}+2mU_{\phi }(\phi )=\Gamma _{\phi }}
Γ
ϕ
{\displaystyle \Gamma _{\phi }}
ϕ
{\displaystyle \phi }
1
2
m
(
d
S
r
d
r
)
2
+
U
r
(
r
)
+
1
2
m
r
2
[
1
sin
2
θ
(
d
S
θ
d
θ
)
2
+
2
m
sin
2
θ
U
θ
(
θ
)
+
Γ
ϕ
]
=
E
.
{\displaystyle {\frac {1}{2m}}\left({\frac {dS_{r}}{dr}}\right)^{2}+U_{r}(r)+{\frac {1}{2mr^{2}}}\left[{\frac {1}{\sin ^{2}\theta }}\left({\frac {dS_{\theta }}{d\theta }}\right)^{2}+{\frac {2m}{\sin ^{2}\theta }}U_{\theta }(\theta )+\Gamma _{\phi }\right]=E.}
次の 常微分方程式には 一般化座標 が含まれます。
ここで は 再び 運動の定数 であり、依存性を排除して HJE を最終的な 常微分方程式 に簡約します。
この方程式を積分すると の解が完成します 。
θ
{\displaystyle \theta }
1
sin
2
θ
(
d
S
θ
d
θ
)
2
+
2
m
sin
2
θ
U
θ
(
θ
)
+
Γ
ϕ
=
Γ
θ
{\displaystyle {\frac {1}{\sin ^{2}\theta }}\left({\frac {dS_{\theta }}{d\theta }}\right)^{2}+{\frac {2m}{\sin ^{2}\theta }}U_{\theta }(\theta )+\Gamma _{\phi }=\Gamma _{\theta }}
Γ
θ
{\displaystyle \Gamma _{\theta }}
θ
{\displaystyle \theta }
1
2
m
(
d
S
r
d
r
)
2
+
U
r
(
r
)
+
Γ
θ
2
m
r
2
=
E
{\displaystyle {\frac {1}{2m}}\left({\frac {dS_{r}}{dr}}\right)^{2}+U_{r}(r)+{\frac {\Gamma _{\theta }}{2mr^{2}}}=E}
S
{\displaystyle S}
楕円円筒座標
楕円柱座標 におけるハミルトニアンは、 楕円 の 焦点 が -軸上 の にある
場合に 書き表される
。ハミルトン・ヤコビ方程式は、 が類似した形を持つという条件で、これらの座標において完全に分離可能である。
ここで 、 、 は 任意の関数である。完全に分離された解を
HJEに代入すると、次の式が得られる
。
H
=
p
μ
2
+
p
ν
2
2
m
a
2
(
sinh
2
μ
+
sin
2
ν
)
+
p
z
2
2
m
+
U
(
μ
,
ν
,
z
)
{\displaystyle H={\frac {p_{\mu }^{2}+p_{\nu }^{2}}{2ma^{2}\left(\sinh ^{2}\mu +\sin ^{2}\nu \right)}}+{\frac {p_{z}^{2}}{2m}}+U(\mu ,\nu ,z)}
±
a
{\displaystyle \pm a}
x
{\displaystyle x}
U
{\displaystyle U}
U
(
μ
,
ν
,
z
)
=
U
μ
(
μ
)
+
U
ν
(
ν
)
sinh
2
μ
+
sin
2
ν
+
U
z
(
z
)
{\displaystyle U(\mu ,\nu ,z)={\frac {U_{\mu }(\mu )+U_{\nu }(\nu )}{\sinh ^{2}\mu +\sin ^{2}\nu }}+U_{z}(z)}
U
μ
(
μ
)
{\displaystyle U_{\mu }(\mu )}
U
ν
(
ν
)
{\displaystyle U_{\nu }(\nu )}
U
z
(
z
)
{\displaystyle U_{z}(z)}
S
=
S
μ
(
μ
)
+
S
ν
(
ν
)
+
S
z
(
z
)
−
E
t
{\displaystyle S=S_{\mu }(\mu )+S_{\nu }(\nu )+S_{z}(z)-Et}
1
2
m
(
d
S
z
d
z
)
2
+
U
z
(
z
)
+
1
2
m
a
2
(
sinh
2
μ
+
sin
2
ν
)
[
(
d
S
μ
d
μ
)
2
+
(
d
S
ν
d
ν
)
2
+
2
m
a
2
U
μ
(
μ
)
+
2
m
a
2
U
ν
(
ν
)
]
=
E
.
{\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.}
最初の 常微分方程式
を分離すると
、簡約されたハミルトン・ヤコビ方程式が得られます(並べ替えと両辺の分母の乗算後)。
この方程式自体は、2 つの独立した 常微分方程式
に分離でき
、これを解くと、 の完全な解が得られます 。
1
2
m
(
d
S
z
d
z
)
2
+
U
z
(
z
)
=
Γ
z
{\displaystyle {\frac {1}{2m}}\left({\frac {dS_{z}}{dz}}\right)^{2}+U_{z}(z)=\Gamma _{z}}
(
d
S
μ
d
μ
)
2
+
(
d
S
ν
d
ν
)
2
+
2
m
a
2
U
μ
(
μ
)
+
2
m
a
2
U
ν
(
ν
)
=
2
m
a
2
(
sinh
2
μ
+
sin
2
ν
)
(
E
−
Γ
z
)
{\displaystyle \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 )=2ma^{2}\left(\sinh ^{2}\mu +\sin ^{2}\nu \right)\left(E-\Gamma _{z}\right)}
(
d
S
μ
d
μ
)
2
+
2
m
a
2
U
μ
(
μ
)
+
2
m
a
2
(
Γ
z
−
E
)
sinh
2
μ
=
Γ
μ
{\displaystyle \left({\frac {dS_{\mu }}{d\mu }}\right)^{2}+2ma^{2}U_{\mu }(\mu )+2ma^{2}\left(\Gamma _{z}-E\right)\sinh ^{2}\mu =\Gamma _{\mu }}
(
d
S
ν
d
ν
)
2
+
2
m
a
2
U
ν
(
ν
)
+
2
m
a
2
(
Γ
z
−
E
)
sin
2
ν
=
Γ
ν
{\displaystyle \left({\frac {dS_{\nu }}{d\nu }}\right)^{2}+2ma^{2}U_{\nu }(\nu )+2ma^{2}\left(\Gamma _{z}-E\right)\sin ^{2}\nu =\Gamma _{\nu }}
S
{\displaystyle S}
放物線円筒座標
放物線円筒座標 におけるハミルトニアンは 次のように書ける。
H
=
p
σ
2
+
p
τ
2
2
m
(
σ
2
+
τ
2
)
+
p
z
2
2
m
+
U
(
σ
,
τ
,
z
)
.
{\displaystyle H={\frac {p_{\sigma }^{2}+p_{\tau }^{2}}{2m\left(\sigma ^{2}+\tau ^{2}\right)}}+{\frac {p_{z}^{2}}{2m}}+U(\sigma ,\tau ,z).}
ハミルトン・ヤコビ方程式は 、 、 、 が任意の関数
で ある類似の形を持つ
限り、これらの座標で完全に分離可能である。完全に分離された解を
HJEに
代入すると、次の式が得られる。
U
{\displaystyle U}
U
(
σ
,
τ
,
z
)
=
U
σ
(
σ
)
+
U
τ
(
τ
)
σ
2
+
τ
2
+
U
z
(
z
)
{\displaystyle U(\sigma ,\tau ,z)={\frac {U_{\sigma }(\sigma )+U_{\tau }(\tau )}{\sigma ^{2}+\tau ^{2}}}+U_{z}(z)}
U
σ
(
σ
)
{\displaystyle U_{\sigma }(\sigma )}
U
τ
(
τ
)
{\displaystyle U_{\tau }(\tau )}
U
z
(
z
)
{\displaystyle U_{z}(z)}
S
=
S
σ
(
σ
)
+
S
τ
(
τ
)
+
S
z
(
z
)
−
E
t
+
constant
{\displaystyle S=S_{\sigma }(\sigma )+S_{\tau }(\tau )+S_{z}(z)-Et+{\text{constant}}}
1
2
m
(
d
S
z
d
z
)
2
+
U
z
(
z
)
+
1
2
m
(
σ
2
+
τ
2
)
[
(
d
S
σ
d
σ
)
2
+
(
d
S
τ
d
τ
)
2
+
2
m
U
σ
(
σ
)
+
2
m
U
τ
(
τ
)
]
=
E
.
{\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.}
最初の 常微分方程式
を分離すると
、簡約されたハミルトン・ヤコビ方程式が得られます(並べ替えと両辺の分母の乗算後)。
この方程式自体は、2 つの独立した 常微分方程式
に分離でき
、これを解くと、 の完全な解が得られます 。
1
2
m
(
d
S
z
d
z
)
2
+
U
z
(
z
)
=
Γ
z
{\displaystyle {\frac {1}{2m}}\left({\frac {dS_{z}}{dz}}\right)^{2}+U_{z}(z)=\Gamma _{z}}
(
d
S
σ
d
σ
)
2
+
(
d
S
τ
d
τ
)
2
+
2
m
U
σ
(
σ
)
+
2
m
U
τ
(
τ
)
=
2
m
(
σ
2
+
τ
2
)
(
E
−
Γ
z
)
{\displaystyle \left({\frac {dS_{\sigma }}{d\sigma }}\right)^{2}+\left({\frac {dS_{\tau }}{d\tau }}\right)^{2}+2mU_{\sigma }(\sigma )+2mU_{\tau }(\tau )=2m\left(\sigma ^{2}+\tau ^{2}\right)\left(E-\Gamma _{z}\right)}
(
d
S
σ
d
σ
)
2
+
2
m
U
σ
(
σ
)
+
2
m
σ
2
(
Γ
z
−
E
)
=
Γ
σ
{\displaystyle \left({\frac {dS_{\sigma }}{d\sigma }}\right)^{2}+2mU_{\sigma }(\sigma )+2m\sigma ^{2}\left(\Gamma _{z}-E\right)=\Gamma _{\sigma }}
(
d
S
τ
d
τ
)
2
+
2
m
U
τ
(
τ
)
+
2
m
τ
2
(
Γ
z
−
E
)
=
Γ
τ
{\displaystyle \left({\frac {dS_{\tau }}{d\tau }}\right)^{2}+2mU_{\tau }(\tau )+2m\tau ^{2}\left(\Gamma _{z}-E\right)=\Gamma _{\tau }}
S
{\displaystyle S}
波と粒子
光波面と軌跡
HJEは、軌跡と 波面 の二重性を確立します。 [10] たとえば、幾何光学では、光は「光線」としても波としても考えられます。波面は、 時刻に放射された光が 時刻に到達した表面として定義できます 。光線と波面は二重であり、一方がわかれば、もう一方は推測できます。
C
t
{\textstyle {\cal {C}}_{t}}
t
=
0
{\textstyle t=0}
t
{\textstyle t}
より正確には、幾何光学は変分問題であり、「作用」は 経路に沿った移動時間であり、 は 媒体の 屈折率 、は 無限小の弧の長さです。上記の定式化から、オイラー・ラグランジュ定式化を使用して光線の経路を計算できます。または、ハミルトン・ヤコビ方程式を解くことで波面を計算することもできます。1 つを知ることは、もう 1 つを知ることにつながります。
T
{\textstyle T}
T
=
1
c
∫
A
B
n
d
s
{\displaystyle T={\frac {1}{c}}\int _{A}^{B}n\,ds}
n
{\textstyle n}
d
s
{\textstyle ds}
上記の双対性は非常に一般的であり、変分原理から派生する すべての システムに適用されます。つまり、オイラー・ラグランジュ方程式を使用して軌道を計算するか、ハミルトン・ヤコビ方程式を使用して波面を計算します。
時刻において最初に であったシステムの、 時刻 における波面は、 となる 点の集合として定義されます 。 が 既知であれば、運動量は直ちに推定されます。
t
{\textstyle t}
q
0
{\textstyle \mathbf {q} _{0}}
t
0
{\textstyle t_{0}}
q
{\textstyle \mathbf {q} }
S
(
q
,
t
)
=
const
{\textstyle S(\mathbf {q} ,t)={\text{const}}}
S
(
q
,
t
)
{\textstyle S(\mathbf {q} ,t)}
p
=
∂
S
∂
q
.
{\displaystyle \mathbf {p} ={\frac {\partial S}{\partial \mathbf {q} }}.}
がわかれば、 の方程式を解くことで 軌道 の接線が計算されます。 ここで は ラグランジアンです。軌道は の情報から復元されます 。
p
{\textstyle \mathbf {p} }
q
˙
{\textstyle {\dot {\mathbf {q} }}}
∂
L
∂
q
˙
=
p
{\displaystyle {\frac {\partial {\cal {L}}}{\partial {\dot {\mathbf {q} }}}}={\boldsymbol {p}}}
q
˙
{\textstyle {\dot {\mathbf {q} }}}
L
{\textstyle {\cal {L}}}
q
˙
{\textstyle {\dot {\mathbf {q} }}}
シュレーディンガー方程式との関係
関数の等 値 面は、任意の時刻 t で決定できます。 時間の関数としての -等値 面の運動は、等値面上の点から始まる粒子の運動によって定義されます。このような等値面の運動は、 波動方程式に 厳密には従いませんが、 -空間 を移動する 波 として考えることができます 。これを示すために、 S が 波の
位相 を表すものとします。
ここで、は 指数引数を無次元にするために導入された 定数 ( プランク定数) です。 波 の 振幅の変化は、 を 複素数 と することで表すことができます 。ハミルトン・ヤコビ方程式は と書き直され、
これが シュレーディンガー方程式 となります。
S
(
q
,
t
)
{\displaystyle S(\mathbf {q} ,t)}
S
{\displaystyle S}
q
{\displaystyle \mathbf {q} }
q
{\displaystyle \mathbf {q} }
ψ
=
ψ
0
e
i
S
/
ℏ
{\displaystyle \psi =\psi _{0}e^{iS/\hbar }}
ℏ
{\displaystyle \hbar }
S
{\displaystyle S}
ℏ
2
2
m
∇
2
ψ
−
U
ψ
=
ℏ
i
∂
ψ
∂
t
{\displaystyle {\frac {\hbar ^{2}}{2m}}\nabla ^{2}\psi -U\psi ={\frac {\hbar }{i}}{\frac {\partial \psi }{\partial t}}}
逆に、シュレーディンガー方程式と の仮定から始めると 、 次 のように推論できる。 [11]
ψ
{\displaystyle \psi }
1
2
m
(
∇
S
)
2
+
U
+
∂
S
∂
t
=
i
ℏ
2
m
∇
2
S
.
{\displaystyle {\frac {1}{2m}}\left(\nabla S\right)^{2}+U+{\frac {\partial S}{\partial t}}={\frac {i\hbar }{2m}}\nabla ^{2}S.}
上記のシュレーディンガー方程式の古典的極限( )は、ハミルトン・ヤコビ方程式の次の変形と同一になる。
ℏ
→
0
{\displaystyle \hbar \rightarrow 0}
1
2
m
(
∇
S
)
2
+
U
+
∂
S
∂
t
=
0.
{\displaystyle {\frac {1}{2m}}\left(\nabla S\right)^{2}+U+{\frac {\partial S}{\partial t}}=0.}
アプリケーション
重力場におけるHJE
曲がった空間を移動する 静止質量
の粒子に対する 形式 [12]の エネルギー-運動量関係 を用いると、 は アインシュタインの場の方程式から解かれた 計量テンソル (すなわち、 逆計量 ) の 反変 座標であり 、 は 光速 で ある 。4 元運動量を 作用の 4元勾配 に等しく設定すると 、
計量 によって決定される幾何学 、
すなわち 重力場
におけるハミルトン-ヤコビ方程式が得られる 。
g
α
β
P
α
P
β
−
(
m
c
)
2
=
0
{\displaystyle g^{\alpha \beta }P_{\alpha }P_{\beta }-(mc)^{2}=0}
m
{\displaystyle m}
g
α
β
{\displaystyle g^{\alpha \beta }}
c
{\displaystyle c}
P
α
{\displaystyle P_{\alpha }}
S
{\displaystyle S}
P
α
=
−
∂
S
∂
x
α
{\displaystyle P_{\alpha }=-{\frac {\partial S}{\partial x^{\alpha }}}}
g
{\displaystyle g}
g
α
β
∂
S
∂
x
α
∂
S
∂
x
β
−
(
m
c
)
2
=
0
,
{\displaystyle g^{\alpha \beta }{\frac {\partial S}{\partial x^{\alpha }}}{\frac {\partial S}{\partial x^{\beta }}}-(mc)^{2}=0,}
電磁場におけるHJE
真空中の4元ポテンシャル を持つ電磁場中を運動する 静止質量 と電荷を持つ 粒子の場合 、計量テンソルによって決定される幾何学におけるハミルトン・ヤコビ方程式は次のような 形を持ち
、ハミルトン主作用関数について解くこと で粒子の軌道と運動量のさらなる解を得ることができる。 [13]
ここで 、およびは ベクトルポテンシャルのサイクル平均
である。
m
{\displaystyle m}
e
{\displaystyle e}
A
i
=
(
ϕ
,
A
)
{\displaystyle A_{i}=(\phi ,\mathrm {A} )}
g
i
k
=
g
i
k
{\displaystyle g^{ik}=g_{ik}}
g
i
k
(
∂
S
∂
x
i
+
e
c
A
i
)
(
∂
S
∂
x
k
+
e
c
A
k
)
=
m
2
c
2
{\displaystyle g^{ik}\left({\frac {\partial S}{\partial x^{i}}}+{\frac {e}{c}}A_{i}\right)\left({\frac {\partial S}{\partial x^{k}}}+{\frac {e}{c}}A_{k}\right)=m^{2}c^{2}}
S
{\displaystyle S}
x
=
−
e
c
γ
∫
A
z
d
ξ
,
{\displaystyle x=-{\frac {e}{c\gamma }}\int A_{z}\,d\xi ,}
y
=
−
e
c
γ
∫
A
y
d
ξ
,
{\displaystyle y=-{\frac {e}{c\gamma }}\int A_{y}\,d\xi ,}
z
=
−
e
2
2
c
2
γ
2
∫
(
A
2
−
A
2
¯
)
d
ξ
,
{\displaystyle z=-{\frac {e^{2}}{2c^{2}\gamma ^{2}}}\int (\mathrm {A} ^{2}-{\overline {\mathrm {A} ^{2}}})\,d\xi ,}
ξ
=
c
t
−
e
2
2
γ
2
c
2
∫
(
A
2
−
A
2
¯
)
d
ξ
,
{\displaystyle \xi =ct-{\frac {e^{2}}{2\gamma ^{2}c^{2}}}\int (\mathrm {A} ^{2}-{\overline {\mathrm {A} ^{2}}})\,d\xi ,}
p
x
=
−
e
c
A
x
,
p
y
=
−
e
c
A
y
,
{\displaystyle p_{x}=-{\frac {e}{c}}A_{x},\quad p_{y}=-{\frac {e}{c}}A_{y},}
p
z
=
e
2
2
γ
c
(
A
2
−
A
2
¯
)
,
{\displaystyle p_{z}={\frac {e^{2}}{2\gamma c}}(\mathrm {A} ^{2}-{\overline {\mathrm {A} ^{2}}}),}
E
=
c
γ
+
e
2
2
γ
c
(
A
2
−
A
2
¯
)
,
{\displaystyle {\mathcal {E}}=c\gamma +{\frac {e^{2}}{2\gamma c}}(\mathrm {A} ^{2}-{\overline {\mathrm {A} ^{2}}}),}
ξ
=
c
t
−
z
{\displaystyle \xi =ct-z}
γ
2
=
m
2
c
2
+
e
2
c
2
A
¯
2
{\displaystyle \gamma ^{2}=m^{2}c^{2}+{\frac {e^{2}}{c^{2}}}{\overline {A}}^{2}}
A
¯
{\displaystyle {\overline {\mathbf {A} }}}
円偏波
円偏波 の場合 、
E
x
=
E
0
sin
ω
ξ
1
,
E
y
=
E
0
cos
ω
ξ
1
,
{\displaystyle E_{x}=E_{0}\sin \omega \xi _{1},\quad E_{y}=E_{0}\cos \omega \xi _{1},}
A
x
=
c
E
0
ω
cos
ω
ξ
1
,
A
y
=
−
c
E
0
ω
sin
ω
ξ
1
.
{\displaystyle A_{x}={\frac {cE_{0}}{\omega }}\cos \omega \xi _{1},\quad A_{y}=-{\frac {cE_{0}}{\omega }}\sin \omega \xi _{1}.}
したがって
、
ここで は、粒子が 磁場ベクトルに沿った方向の
永久半径 と不変の運動量値を持つ 円軌道に沿って移動することを意味します。
x
=
−
e
c
E
0
ω
sin
ω
ξ
1
,
{\displaystyle x=-{\frac {ecE_{0}}{\omega }}\sin \omega \xi _{1},}
y
=
−
e
c
E
0
ω
cos
ω
ξ
1
,
{\displaystyle y=-{\frac {ecE_{0}}{\omega }}\cos \omega \xi _{1},}
p
x
=
−
e
E
0
ω
cos
ω
ξ
1
,
{\displaystyle p_{x}=-{\frac {eE_{0}}{\omega }}\cos \omega \xi _{1},}
p
y
=
e
E
0
ω
sin
ω
ξ
1
,
{\displaystyle p_{y}={\frac {eE_{0}}{\omega }}\sin \omega \xi _{1},}
ξ
1
=
ξ
/
c
{\displaystyle \xi _{1}=\xi /c}
e
c
E
0
/
γ
ω
2
{\displaystyle ecE_{0}/\gamma \omega ^{2}}
e
E
0
/
ω
2
{\displaystyle eE_{0}/\omega ^{2}}
単色直線偏波平面波
軸に沿った方向を向いた 電界を持つ平坦で単色の直線偏波の場合
、
粒子の軌跡は 8 の字型になり、その軸は電界 ベクトルに沿って配向されます。
E
{\displaystyle E}
y
{\displaystyle y}
E
y
=
E
0
cos
ω
ξ
1
,
{\displaystyle E_{y}=E_{0}\cos \omega \xi _{1},}
A
y
=
−
c
E
0
ω
sin
ω
ξ
1
,
{\displaystyle A_{y}=-{\frac {cE_{0}}{\omega }}\sin \omega \xi _{1},}
x
=
const
,
{\displaystyle x={\text{const}},}
y
0
=
−
e
c
E
0
γ
ω
2
,
{\displaystyle y_{0}=-{\frac {ecE_{0}}{\gamma \omega ^{2}}},}
y
=
y
0
cos
ω
ξ
1
,
z
=
C
z
y
0
sin
2
ω
ξ
1
,
{\displaystyle y=y_{0}\cos \omega \xi _{1},\quad z=C_{z}y_{0}\sin 2\omega \xi _{1},}
C
z
=
e
E
0
8
γ
ω
,
γ
2
=
m
2
c
2
+
e
2
E
0
2
2
ω
2
,
{\displaystyle C_{z}={\frac {eE_{0}}{8\gamma \omega }},\quad \gamma ^{2}=m^{2}c^{2}+{\frac {e^{2}E_{0}^{2}}{2\omega ^{2}}},}
p
x
=
0
,
{\displaystyle p_{x}=0,}
p
y
,
0
=
e
E
0
ω
,
{\displaystyle p_{y,0}={\frac {eE_{0}}{\omega }},}
p
y
=
p
y
,
0
sin
ω
ξ
1
,
{\displaystyle p_{y}=p_{y,0}\sin \omega \xi _{1},}
p
z
=
−
2
C
z
p
y
,
0
cos
2
ω
ξ
1
{\displaystyle p_{z}=-2C_{z}p_{y,0}\cos 2\omega \xi _{1}}
E
{\displaystyle E}
ソレノイド磁場を持つ電磁波
軸方向(ソレノイド)磁場を持つ電磁波の場合: [14]
したがって、
有効半径 、誘導性 、巻き数 、ソレノイド巻き線を流れる電流の大きさを持つソレノイド内の磁場の大きさは
です 。ソレノイド磁場の軸対称性により、粒子の運動は、 任意の方位角でソレノイド軸に垂直な平面内で8の字軌道に沿って発生します 。
E
=
E
ϕ
=
ω
ρ
0
c
B
0
cos
ω
ξ
1
,
{\displaystyle E=E_{\phi }={\frac {\omega \rho _{0}}{c}}B_{0}\cos \omega \xi _{1},}
A
ϕ
=
−
ρ
0
B
0
sin
ω
ξ
1
=
−
L
s
π
ρ
0
N
s
I
0
sin
ω
ξ
1
,
{\displaystyle A_{\phi }=-\rho _{0}B_{0}\sin \omega \xi _{1}=-{\frac {L_{s}}{\pi \rho _{0}N_{s}}}I_{0}\sin \omega \xi _{1},}
x
=
constant
,
{\displaystyle x={\text{constant}},}
y
0
=
−
e
ρ
0
B
0
γ
ω
,
{\displaystyle y_{0}=-{\frac {e\rho _{0}B_{0}}{\gamma \omega }},}
y
=
y
0
cos
ω
ξ
1
,
{\displaystyle y=y_{0}\cos \omega \xi _{1},}
z
=
C
z
y
0
sin
2
ω
ξ
1
,
{\displaystyle z=C_{z}y_{0}\sin 2\omega \xi _{1},}
C
z
=
e
ρ
0
B
0
8
c
γ
,
{\displaystyle C_{z}={\frac {e\rho _{0}B_{0}}{8c\gamma }},}
γ
2
=
m
2
c
2
+
e
2
ρ
0
2
B
0
2
2
c
2
,
{\displaystyle \gamma ^{2}=m^{2}c^{2}+{\frac {e^{2}\rho _{0}^{2}B_{0}^{2}}{2c^{2}}},}
p
x
=
0
,
{\displaystyle p_{x}=0,}
p
y
,
0
=
e
ρ
0
B
0
c
,
{\displaystyle p_{y,0}={\frac {e\rho _{0}B_{0}}{c}},}
p
y
=
p
y
,
0
sin
ω
ξ
1
,
{\displaystyle p_{y}=p_{y,0}\sin \omega \xi _{1},}
p
z
=
−
2
C
z
p
y
,
0
cos
2
ω
ξ
1
,
{\displaystyle p_{z}=-2C_{z}p_{y,0}\cos 2\omega \xi _{1},}
B
0
{\displaystyle B_{0}}
ρ
0
{\displaystyle \rho _{0}}
L
s
{\displaystyle L_{s}}
N
s
{\displaystyle N_{s}}
I
0
{\displaystyle I_{0}}
y
z
{\displaystyle yz}
φ
{\displaystyle \varphi }
参照
数学ポータル 物理学ポータル
参考文献
^ ゴールドスタイン、ハーバート (1980年)。 古典力学 ( 第2版)。マサチューセッツ州レディング:アディソンウェスレー。pp.484-492。ISBN 978-0-201-02918-5 。 (特に491ページの最後の段落から始まる議論)
^ 桜井, JJ (1994). 現代量子力学 (改訂版). マサチューセッツ州レディング: アディソンウェスレー. pp. 103–107. ISBN 0-201-53929-2 。
^ Kálmán, Rudolf E. (1963). 「最適制御理論と変分法」。 Bellman, Richard (編)。 数学的最適化手法 。 バークレー: カリフォルニア大学出版局。 pp. 309–331。 OCLC 1033974。
^ Hand, LN; Finch, JD (2008). 解析力学 . ケンブリッジ大学出版局. ISBN 978-0-521-57572-0 。
^ abc ゴールドスタイン、ハーバート; プール、チャールズ P.; サフコ、ジョン L. (2008). 古典力学 (3、[後編] 版). サンフランシスコ ミュンヘン: アディソン ウェズリー. ISBN 978-0-201-65702-9 。
^ クーパースミス、ジェニファー(2017)。怠惰な宇宙:最小作用の原理入門。オックスフォード、イギリス/ニューヨーク、ニューヨーク:オックスフォード大学出版局 。ISBN 978-0-19-874304-0 。
^ Hand, LN; Finch, JD (2008). 解析力学 . ケンブリッジ大学出版局. ISBN 978-0-521-57572-0 。
^ ゴールドスタイン、ハーバート (1980年)。 古典力学 ( 第2版)。マサチューセッツ州レディング:アディソンウェスレー。p.440。ISBN 978-0-201-02918-5 。
^ Hanc, Jozef; Taylor, Edwin F.; Tuleja, Slavomir (2005-07-01). 「1次元および2次元の変分力学」. American Journal of Physics . 73 (7): 603–610. Bibcode :2005AmJPh..73..603H. doi :10.1119/1.1848516. ISSN 0002-9505.
^ Houchmandzadeh, Bahram (2020). 「ハミルトン・ヤコビ方程式:代替アプローチ」. American Journal of Physics . 85 (5): 10.1119/10.0000781. arXiv : 1910.09414 . Bibcode :2020AmJPh..88..353H. doi :10.1119/10.0000781. S2CID 204800598.
^ ゴールドスタイン、ハーバート (1980年)。 古典力学 ( 第2版)。マサチューセッツ州レディング:アディソンウェスレー。pp.490-491。ISBN 978-0-201-02918-5 。
^ ウィーラー、ジョン; ミスナー、チャールズ; ソーン、キップ (1973)。 重力 。WHフリーマン&カンパニー 。pp.649、1188。ISBN 978-0-7167-0344-0 。
^ Landau, L. ; Lifshitz, E. (1959). 場の古典理論 . マサチューセッツ州レディング: Addison-Wesley. OCLC 17966515.
^ EV Shun'ko; DE Stevenson; VS Belkin (2014). 「6~100 eVの範囲でプラズマ電子エネルギーを制御可能な誘導結合プラズマリアクター」. IEEE Transactions on Plasma Science . 42, part II (3): 774–785. Bibcode :2014ITPS...42..774S. doi :10.1109/TPS.2014.2299954. S2CID 34765246.
さらに読む
アーノルド、VI (1989)。 古典力学の数学的手法 (第2版)。ニューヨーク:シュプリンガー 。ISBN 0-387-96890-3 。
ハミルトン、W. (1833)。「光と惑星の軌道を特性関数の係数で表す一般的な方法について」 (PDF) 。 ダブリン大学レビュー : 795–826。
Hamilton, W. (1834)。「光学に以前適用されていた一般的な数学的手法の力学への応用について」 (PDF) 。 英国協会レポート : 513–518。
Fetter, A. & Walecka, J. (2003). 粒子と連続体の理論力学 . Dover Books. ISBN 978-0-486-43261-8 。
Landau, LD ; Lifshitz, EM (1975). 力学 . アムステルダム: Elsevier.
桜井, JJ (1985). 現代量子力学 . ベンジャミン/カミングス出版. ISBN 978-0-8053-7501-5 。
Jacobi、CGJ (1884)、 Vorlesungen über Dynamik 、CGJ Jacobi's Gesammelte Werke (ドイツ語)、ベルリン: G. Reimer、 OL 14009561M
中根美千代フレイザー、クレイグ G. (2002)。 「ハミルトン・ヤコビ・ダイナミクスの初期の歴史」。 ケンタウルス座 。 44 (3–4): 161–227。 土井 :10.1111/j.1600-0498.2002.tb00613.x。 PMID 17357243。