古典力学の定式化
物理学において、 ウィリアム・ローワン・ハミルトン と カール・グスタフ・ヤコビ・ヤコビ にちなんで名付けられた ハミルトン・ヤコビ方程式は、 古典力学 の別の定式化であり、 ニュートンの運動の法則 、 ラグランジアン力学 、 ハミルトン力学 などの他の定式化と同等です 。
ハミルトン・ヤコビ方程式は、粒子の運動を波として表す力学の定式化である。この意味で、この方程式は、光の伝播と粒子の運動の間に類似点を見出すという、理論物理学における長年の目標(少なくとも18世紀の ヨハン・ベルヌーイ に遡る)を達成した。力学系が従う波動方程式は、 後述するように シュレーディンガー方程式と類似しているが、同一ではない。このため、ハミルトン・ヤコビ方程式は 古典力学から 量子力学 への「最も近いアプローチ」と考えられている 。 [1] [2] この関係の定性的形式は、 ハミルトンの光機械アナロジー と呼ばれる。
数学において、ハミルトン・ヤコビ方程式は、 変分法の 問題の一般化における極値 幾何学 を記述する 必要条件である。これは、 動的計画法 における ハミルトン・ヤコビ・ベルマン方程式 の特殊なケースとして理解することができる 。 [3]
概要 ハミルトン・ヤコビ方程式は、1階の 非線形 偏微分方程式である。
− ∂ 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] : 291 この解は、最小作用の原理 で使用される形式の不定積分によって
系ラグランジアンと関連付けることができる 。 [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}}\ \mathrm {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_{\mathcal {L}}(\mathbf {q} ,\mathbf {\dot {q}} ,t)=\left\{\partial ^{2}{\mathcal {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 {\mathcal {L}}}{\partial {\dot {q}}^{i}}}=\sum _{j=1}^{n}\left({\frac {\partial ^{2}{\mathcal {L}}}{\partial {\dot {q}}^{i}\partial {\dot {q}}^{j}}}{\ddot {q}}^{j}+{\frac {\partial ^{2}{\mathcal {L}}}{\partial {\dot {q}}^{i}\partial {q}^{j}}}{\dot {q}}^{j}\right)+{\frac {\partial ^{2}{\mathcal {L}}}{\partial {\dot {q}}^{i}\partial t}},\qquad i=1,\ldots ,n,} n × n {\displaystyle n\times n} H L {\displaystyle H_{\mathcal {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 {\mathcal {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 {\mathcal {L}}}{\partial {\dot {q}}^{i}}}\right|_{\mathbf {\dot {q}} =\mathbf {v} }\!\!\!\!\!\!\!,\quad i=1,\ldots ,n.}
証拠 以下の証明では、配置空間が の開集合であると仮定していますが、 基礎となる手法は任意の 空間 にも同様に当てはまります。この証明の文脈において、カリグラフィ文字は 作用関数、イタリック体は ハミルトンの主関数を表します。 R n , {\displaystyle \mathbb {R} ^{n},} S {\displaystyle {\mathcal {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 {\mathcal {S}}_{\delta \xi }[\gamma ,t_{1},t_{0}]} S {\displaystyle {\mathcal {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 {\mathcal {S}}_{\delta \xi }[\xi ,t_{1},t_{0}]=\int _{t_{0}}^{t_{1}}\left({\frac {\partial {\mathcal {L}}}{\partial \mathbf {q} }}-{\frac {d}{dt}}{\frac {\partial {\mathcal {L}}}{\partial \mathbf {\dot {q}} }}\right)\delta \xi \,dt+{\frac {\partial {\mathcal {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 {\mathcal {S}}_{\delta \xi }[\xi ,t;t_{0}]=\left.{\frac {\partial {\mathcal {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 {\mathcal {S}}_{\delta \gamma }[\gamma ,t]{\overset {\text{def}}{{}={}}}\left.{\frac {d{\mathcal {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 {\mathcal {S}}_{\delta \xi }[\xi ,t;t_{0}]} t > t 0 , {\displaystyle t>t_{0},} δ γ {\displaystyle \delta \gamma }
機械システムのハミルトニアン が与えられた場合、ハミルトン・ヤコビ方程式は ハミルトンの主関数に対する 1階の 非線形 偏微分方程式 である。 [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 {\mathcal {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}} -{\mathcal {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 {\mathcal {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}}={\mathcal {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).}
共役運動量は 一般化座標に関する
1次導関数に対応する。 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 }
ハミルトン・ヤコビ方程式は、 一般化座標 と時間 の関数 に関する単一の一次偏微分方程式 で ある。一般化運動量は、古典 作用 の微分としてのみ現れる 。 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 {\begin{aligned}&\mathbf {p} ={\frac {\partial G_{2}}{\partial \mathbf {q} }},\quad \mathbf {Q} ={\frac {\partial G_{2}}{\partial \mathbf {P} }},\quad \\&K(\mathbf {Q} ,\mathbf {P} ,t)=H(\mathbf {q} ,\mathbf {p} ,t)+{\frac {\partial G_{2}}{\partial t}}\end{aligned}}} 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 {\begin{aligned}&\mathbf {p} ={\frac {\partial G_{2}}{\partial \mathbf {q} }}={\frac {\partial S}{\partial \mathbf {q} }}\\[1ex]\rightarrow {}&H(\mathbf {q} ,\mathbf {p} ,t)+{\partial G_{2} \over \partial t}=0\\[1ex]\rightarrow {}&H{\left(\mathbf {q} ,{\frac {\partial S}{\partial \mathbf {q} }},t\right)}+{\partial S \over \partial t}=0.\end{aligned}}}
を解くと、 わかりやすいように成分で書かれた 便利な方程式
や 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 + 1 2 m a 2 ( sinh 2 μ + sin 2 ν ) [ ( d S μ d μ ) 2 + ( d S ν d ν ) 2 ] + U z ( z ) + 1 sinh 2 μ + sin 2 ν [ U μ ( μ ) + U ν ( ν ) ] = E . {\displaystyle {\begin{aligned}{\frac {1}{2m}}\left({\frac {dS_{z}}{dz}}\right)^{2}+{\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}\right]&\\{}+U_{z}(z)+{\frac {1}{\sinh ^{2}\mu +\sin ^{2}\nu }}\left[U_{\mu }(\mu )+U_{\nu }(\nu )\right]&=E.\end{aligned}}}
最初の 常微分方程式 を分離すると
、簡約されたハミルトン・ヤコビ方程式が得られます(並べ替えて両辺に分母を乗算した後)。
この方程式自体は、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 μ = Γ μ ( d S ν d ν ) 2 + 2 m a 2 U ν ( ν ) + 2 m a 2 ( Γ z − E ) sin 2 ν = Γ ν {\displaystyle {\begin{alignedat}{4}\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 }\\\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 }\end{alignedat}}} 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 + 1 2 m ( σ 2 + τ 2 ) [ ( d S σ d σ ) 2 + ( d S τ d τ ) 2 ] + U z ( z ) + 1 σ 2 + τ 2 [ U σ ( σ ) + U τ ( τ ) ] = E . {\displaystyle {\begin{aligned}{\frac {1}{2m}}\left({\frac {dS_{z}}{dz}}\right)^{2}+{\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}\right]&\\{}+U_{z}(z)+{\frac {1}{\sigma ^{2}+\tau ^{2}}}\left[U_{\sigma }(\sigma )+U_{\tau }(\tau )\right]&=E.\end{aligned}}}
最初の 常微分方程式 を分離すると
、簡約されたハミルトン・ヤコビ方程式が得られます(並べ替えて両辺に分母を乗算した後)。
この方程式自体は、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 σ ( σ ) + U τ ( τ ) ] = 2 m ( σ 2 + τ 2 ) ( E − Γ z ) {\displaystyle \left({\frac {dS_{\sigma }}{d\sigma }}\right)^{2}+\left({\frac {dS_{\tau }}{d\tau }}\right)^{2}+2m\left[U_{\sigma }(\sigma )+U_{\tau }(\tau )\right]=2m\left(\sigma ^{2}+\tau ^{2}\right)\left(E-\Gamma _{z}\right)} ( d S σ d σ ) 2 + 2 m U σ ( σ ) + 2 m σ 2 ( Γ z − E ) = Γ σ ( d S τ d τ ) 2 + 2 m U τ ( τ ) + 2 m τ 2 ( Γ z − E ) = Γ τ {\displaystyle {\begin{alignedat}{4}\left({\frac {dS_{\sigma }}{d\sigma }}\right)^{2}&+\,&2mU_{\sigma }(\sigma )&+\,&2m\sigma ^{2}\left(\Gamma _{z}-E\right)&=\,&\Gamma _{\sigma }\\\left({\frac {dS_{\tau }}{d\tau }}\right)^{2}&+\,&2mU_{\tau }(\tau )&+\,&2m\tau ^{2}\left(\Gamma _{z}-E\right)&=\,&\Gamma _{\tau }\end{alignedat}}} S {\displaystyle S}
波と粒子
光波面と軌跡 HJEは、軌跡と 波面 の二重性を確立する。 [10] 例えば、幾何光学では、光は「光線」と波のどちらとしても考えることができる。波面は、 時刻 に放射された光が 時刻 に到達した面として定義できる 。光線と波面は二重性を持ち、一方が分かれば、他方を推測することができる。 C t {\textstyle {\mathcal {C}}_{t}} t = 0 {\textstyle t=0} t {\textstyle t}
より正確には、幾何光学は変分問題であり、「作用」は経路に沿った 移動時間であり 、 媒質の 屈折率 、 無限小の弧長です。上記の定式化から、オイラー・ラグランジュの定式化を用いて光線の経路を計算することができます。あるいは、ハミルトン・ヤコビ方程式を解くことで波面を計算することもできます。どちらか一方を知ることで、もう一方を知ることができます。 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 {\mathcal {L}}}{\partial {\dot {\mathbf {q} }}}}={\boldsymbol {p}}} q ˙ {\textstyle {\dot {\mathbf {q} }}} L {\textstyle {\mathcal {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 ψ 0 ψ 0 . {\displaystyle {\frac {1}{2m}}\left(\nabla S\right)^{2}+U+{\frac {\partial S}{\partial t}}={\frac {i\hbar }{2m}}{\frac {\nabla ^{2}\psi _{0}}{\psi _{0}}}.}
上記のシュレーディンガー方程式の古典極限( )は、ハミルトン・ヤコビ方程式の次の変形と同一となる。 ℏ → 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 真空中の四元ポテンシャル を持つ電磁場中を運動する 静止質量 と電荷を持つ 粒子の場合 、計量テンソルによって決定される幾何学におけるハミルトン・ヤコビ方程式は 形を持ち
、ハミルトン主作用関数について解く ことで粒子の軌道と運動量のさらなる解を得ること ができる。 [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 ξ , y = − e c γ ∫ A y d ξ , z = − e 2 2 c 2 γ 2 ∫ ( A 2 − A 2 ¯ ) d ξ , ξ = c t − e 2 2 γ 2 c 2 ∫ ( A 2 − A 2 ¯ ) d ξ , p x = − e c A x , p y = − e c A y , p z = e 2 2 γ c ( A 2 − A 2 ¯ ) , E = c γ + e 2 2 γ c ( A 2 − A 2 ¯ ) , {\displaystyle {\begin{aligned}x&=-{\frac {e}{c\gamma }}\int A_{z}\,d\xi ,&y&=-{\frac {e}{c\gamma }}\int A_{y}\,d\xi ,\\[1ex]z&=-{\frac {e^{2}}{2c^{2}\gamma ^{2}}}\int \left(\mathrm {A} ^{2}-{\overline {\mathrm {A} ^{2}}}\right)\,d\xi ,&\xi &=ct-{\frac {e^{2}}{2\gamma ^{2}c^{2}}}\int \left(\mathrm {A} ^{2}-{\overline {\mathrm {A} ^{2}}}\right)\,d\xi ,\\[1ex]p_{x}&=-{\frac {e}{c}}A_{x},&p_{y}&=-{\frac {e}{c}}A_{y},\\[1ex]p_{z}&={\frac {e^{2}}{2\gamma c}}\left(\mathrm {A} ^{2}-{\overline {\mathrm {A} ^{2}}}\right),&{\mathcal {E}}&=c\gamma +{\frac {e^{2}}{2\gamma c}}\left(\mathrm {A} ^{2}-{\overline {\mathrm {A} ^{2}}}\right),\end{aligned}}} ξ = 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 , A x = c E 0 ω cos ω ξ 1 , A y = − c E 0 ω sin ω ξ 1 . {\displaystyle {\begin{aligned}E_{x}&=E_{0}\sin \omega \xi _{1},&E_{y}&=E_{0}\cos \omega \xi _{1},\\[1ex]A_{x}&={\frac {cE_{0}}{\omega }}\cos \omega \xi _{1},&A_{y}&=-{\frac {cE_{0}}{\omega }}\sin \omega \xi _{1}.\end{aligned}}}
したがって x = − e c E 0 ω sin ω ξ 1 , y = − e c E 0 ω cos ω ξ 1 , p x = − e E 0 ω cos ω ξ 1 , p y = e E 0 ω sin ω ξ 1 , {\displaystyle {\begin{aligned}x&=-{\frac {ecE_{0}}{\omega }}\sin \omega \xi _{1},&y&=-{\frac {ecE_{0}}{\omega }}\cos \omega \xi _{1},\\[1ex]p_{x}&=-{\frac {eE_{0}}{\omega }}\cos \omega \xi _{1},&p_{y}&={\frac {eE_{0}}{\omega }}\sin \omega \xi _{1},\end{aligned}}}
ここで、これは、粒子が 磁場ベクトルに沿った
方向の、一定の半径 と不変の運動量値を 持つ円軌道に沿って移動することを意味します。 ξ 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}}
単色直線偏波平面波 平坦で単色の直線偏波の場合、 軸に沿って電界が向くため 、 E {\displaystyle E} y {\displaystyle y} E y = E 0 cos ω ξ 1 , A y = − c E 0 ω sin ω ξ 1 , {\displaystyle {\begin{aligned}E_{y}&=E_{0}\cos \omega \xi _{1},&A_{y}&=-{\frac {cE_{0}}{\omega }}\sin \omega \xi _{1},\end{aligned}}} x = const , y = y 0 cos ω ξ 1 , y 0 = − e c E 0 γ ω 2 , z = C z y 0 sin 2 ω ξ 1 , C z = e E 0 8 γ ω , γ 2 = m 2 c 2 + e 2 E 0 2 2 ω 2 , {\displaystyle {\begin{aligned}x&={\text{const}},\\[1ex]y&=y_{0}\cos \omega \xi _{1},&y_{0}&=-{\frac {ecE_{0}}{\gamma \omega ^{2}}},\\[1ex]z&=C_{z}y_{0}\sin 2\omega \xi _{1},&C_{z}&={\frac {eE_{0}}{8\gamma \omega }},\\[1ex]\gamma ^{2}&=m^{2}c^{2}+{\frac {e^{2}E_{0}^{2}}{2\omega ^{2}}},\end{aligned}}} p x = 0 , p y = p y , 0 sin ω ξ 1 , p y , 0 = e E 0 ω , p z = − 2 C z p y , 0 cos 2 ω ξ 1 {\displaystyle {\begin{aligned}p_{x}&=0,\\[1ex]p_{y}&=p_{y,0}\sin \omega \xi _{1},&p_{y,0}&={\frac {eE_{0}}{\omega }},\\[1ex]p_{z}&=-2C_{z}p_{y,0}\cos 2\omega \xi _{1}\end{aligned}}}
これは、粒子の軌道が電場 ベクトルに沿って軸を向いた 8 の字を描くことを意味します。 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 , y = y 0 cos ω ξ 1 , y 0 = − e ρ 0 B 0 γ ω , z = C z y 0 sin 2 ω ξ 1 , C z = e ρ 0 B 0 8 c γ , γ 2 = m 2 c 2 + e 2 ρ 0 2 B 0 2 2 c 2 , {\displaystyle {\begin{aligned}x&={\text{constant}},\\y&=y_{0}\cos \omega \xi _{1},&y_{0}&=-{\frac {e\rho _{0}B_{0}}{\gamma \omega }},\\z&=C_{z}y_{0}\sin 2\omega \xi _{1},&C_{z}&={\frac {e\rho _{0}B_{0}}{8c\gamma }},\\\gamma ^{2}&=m^{2}c^{2}+{\frac {e^{2}\rho _{0}^{2}B_{0}^{2}}{2c^{2}}},\end{aligned}}} p x = 0 , p y = p y , 0 sin ω ξ 1 , p y , 0 = e ρ 0 B 0 c , p z = − 2 C z p y , 0 cos 2 ω ξ 1 , {\displaystyle {\begin{aligned}p_{x}&=0,\\p_{y}&=p_{y,0}\sin \omega \xi _{1},&p_{y,0}&={\frac {e\rho _{0}B_{0}}{c}},\\p_{z}&=-2C_{z}p_{y,0}\cos 2\omega \xi _{1},\end{aligned}}} 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 。 ^ ルドルフ・E・カルマン (1963). 「最適制御理論と変分法」. リチャード・ベルマン編著. 『数学的最適化技術』 . バークレー: カリフォルニア大学出版局. 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 eVから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 . ハミルトン, W. (1834). 「光学に応用された一般的な数学的手法の力学への応用について」 (PDF) . 英国協会報告書 : 513–518 . フェッター、A. & ワレッカ、J. (2003). 粒子と連続体の理論力学 . ドーバーブックス. 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。