汎関数計算を使った解析力学
; draft
解析力学において「変分」を用いて導出される諸々を, 汎関数計算を用いて書き直してみたい. 飽くまで「試み」なので正確性は保証されない. 近い内容の文献がありましたら教えていただけると嬉しいです.
最小作用の原理
時間 t に依存する一般化座標と呼ばれるパラメータ qi(t) に対して, 作用 action と呼ばれる汎関数 S[qi] が存在し, 物理現象において座標 qi は作用 S[qi] が最小となるような経路が選ばれる. つまり,
δqj(t′)δS[qi(t)]=0.(t1<t′<t2)この古典的原理を最小作用の原理という.
変分を用いた最小作用の原理 δS[qi]=0 は 1 次の冪展開を用いて以下のように書き直される:
δS[qi]=∫t1t2dt′δqj(t′)δS[qi]δqj(t′)=0.
δqi(t′) は t1<t′<t2 で任意だから, 変分を用いた最小作用の原理は汎関数微分を用いた上の式と等価である.
Euler–Lagrange の運動方程式
作用は, 座標と時間に関する Lagrangian L(qi,q˙i,t) を用いて以下のように表される:
S[qi]=∫t1t2dtL(qi,q˙i,t).
最小作用の原理に対し, この Lagrangian が満たすべき条件を求めよう. t1<t<t2 における作用の汎関数微分は
δqj(t′)δS[qi(t)]=h→0limhS[qi(t)+hδjiδ(t−t′)]−S[qi(t)]=h→0limh1[∫t1t2dtL(qi(t)+hδjiδ(t−t′),q˙i(t)+hδjiδ˙(t−t′),t)−∫t1t2dtL(qi(t),q˙i(t),t)]=h→0limh1∫t1t2dt[∂qi(t)∂Lhδjiδ(t−t′)+∂q˙i(t)∂Lhδjiδ˙(t−t′)+o(ε2)]=∫t1t2dt[∂qj(t)∂Lδ(t−t′)+∂q˙j(t)∂Lδ˙(t−t′)]=∫t1t2dt[∂qj(t)∂L−dtd(∂q˙j(t)∂L)]δ(t−t′)+[∂q˙j(t)∂Lδ(t−t′)]t=t1t=t2.
したがって, t1<t′<t2 において,
δqi(t′)δS[qi(t)]=∂qi(t′)∂L−dt′d(∂q˙i(t′)∂L).
最小作用の原理を満たすとき, Lagrangian L(qi,q˙i,t) は以下の Euler–Lagrange の運動方程式を満たす:
∂qi∂L−dtd(∂q˙i∂L)=0.
または, 冪展開を用いて作用の変分を計算すると,
δS[qi(t)]=∫t1t2dt′δqj(t′)δS[qi(t)]δqj(t′)=∫t1t2dt′{∫t1t2dt[∂qj(t)∂L−dtd(∂q˙j(t)∂L)]δ(t−t′)+[∂q˙j(t)∂Lδ(t−t′)]t=t1t=t2}δqj(t′)=∫t1t2dt′∫t1t2dtδqj(t′)[∂qj(t)∂L−dtd(∂q˙j(t)∂L)]δ(t−t′)+[∫t1t2dt′δqj(t′)∂q˙j(t)∂Lδ(t−t′)]t=t1t=t2=∫t1t2dtδqj(t)[∂qj(t)∂L−dtd(∂q˙j(t)∂L)]+[δqj(t)∂q˙j(t)∂L]t=t1t=t2.
δqi(t) は t1<t<t2 で任意だから, 原理 δS[qi]=0 より, 同じ運動方程式が得られる.
Noether の定理
Noether の定理は, 時間と座標の微小変換 t↦t′=t+δt, qi(t)↦q′i(t′)=qi(t)+δqi(t) に対して保存量が存在することを主張する. これを汎関数微分で導出しよう. 混乱を防ぐため, この節では作用の積分変数を添字で書く. つまり, St[qi]≡S[qi(t)] である. 変分 δSt[qi] を t, q に関して1次で冪展開すると1,
δSt[qi]=∫t1t2dt0δqi(t0)δSt[qi]δLqi(t0)+∫t1t2dt0δτ(t0)δSτ(t)[qi]τ(t)=tδt(t0).
それぞれの汎関数微分を計算すると,
δqi(t0)δSt[qi]=∫t1t2dt[∂qi(t)∂L−dtd(∂q˙i(t)∂L)]δ(t−t0)+[∂q˙i(t)∂Lδ(t−t0)]t=t1t=t2,
δτ(t0)δSτ(t)[qi]=h→0limhSτ(t)+hδ(t−t0)[qi]−Sτ(t)[qi]=h→0limh1{∫τ(t1)τ(t2)d(τ(t)+hδ(t−t0))L(qi(τ(t)+hδ(t−t0)),q˙i(τ(t)+hδ(t−t0)),τ(t)+hδ(t−t0))−∫τ(t1)τ(t2)dτ(t)L(qi(τ(t)),q˙i(τ(t)),τ(t))}=h→0limh1{∫τ(t1)τ(t2)dτ(t)L(qi(τ(t)+hδ(t−t0)),q˙i(τ(t)+hδ(t−t0)),τ(t)+hδ(t−t0))+∫τ(t1)τ(t2)dτ(t)hL(qi(τ(t)+hδ(t−t0)),q˙i(τ(t)+hδ(t−t0)),τ(t)+hδ(t−t0))dτ(t)dδ(t−t0)−∫τ(t1)τ(t2)dτ(t)L(qi(τ(t)),q˙i(τ(t)),τ(t))}=h→0limh1{∫τ(t1)τ(t2)dτ(t)[∂qi(τ(t))∂Lq˙i(τ(t))hδ(t−t0)+∂q˙i(τ(t))∂Lq¨i(τ(t))hδ(t−t0)+∂τ(t)∂Lhδ(t−t0)]+∫τ(t1)τ(t2)dτ(t)hL(qi(τ(t)),q˙i(τ(t)),τ(t))dτ(t)dδ(t−t0)+O(h2)}=∫τ(t1)τ(t2)dτ(t){[∂qi(τ(t))∂Lq˙i(τ(t))+∂q˙i(τ(t))∂Lq¨i(τ(t))+∂τ(t)∂L]δ(t−t0)+L(qi(τ(t)),q˙i(τ(t)),τ(t))dτ(t)dδ(t−t0)}=∫τ(t1)τ(t2)dτ(t)dτ(t)d[L(qi(τ(t)),q˙i(τ(t)),τ(t))δ(t−t0)]=[L(qi(τ(t)),q˙i(τ(t)),τ(t))δ(t−t0)]τ(t)=τ(t1)τ(t)=τ(t2)
∴δτ(t0)δSτ(t)[qi]τ(t)=t=[L(qi(t),q˙i(t),t)δ(t−t0)]t=t1t=t2.
これらを変分の式に代入すると,
δSt[qi]=∫t1t2dt0δLqi(t0){∫t1t2dt[∂qi(t)∂L−dtd(∂q˙i(t)∂L)]δ(t−t0)+[∂q˙i(t)∂Lδ(t−t0)]t=t1t=t2}+∫t1t2dt0δt(t0)[L(qi(t),q˙i(t),t)δ(t−t0)]t=t1t=t2=∫t1t2dt0∫t1t2dtδLqi(t0)[∂qi(t)∂L−dtd(∂q˙i(t)∂L)]δ(t−t0)+[∫t1t2dt0δLqi(t0)∂q˙i(t)∂Lδ(t−t0)]t=t1t=t2+[∫t1t2dt0δt(t0)L(qi(t),q˙i(t),t)δ(t−t0)]t=t1t=t2=∫t1t2dtδLqi(t)[∂qi(t)∂L−dtd(∂q˙i(t)∂L)]+[δLqi(t)∂q˙i(t)∂L+δt(t)L(qi(t),q˙i(t),t)]t=t1t=t2=∫t1t2dtδLqi(t)[∂qi(t)∂L−dtd(∂q˙i(t)∂L)]+[δqi(t)∂q˙i(t)∂L−δt(q˙i(t)∂q˙i(t)∂L−L(qi(t),q˙i(t),t))]t=t1t=t2.
ここで, 最後の式の第一項は Euler–Lagrange の運動方程式より消え, 第二項の t1, t2 は任意である. したがって, δS=0 であるとすると, 対応する保存量が得られる:
時間の微小変換 t↦t′=t+δt に対し, 座標が qi(t)↦q′i(t′)=qi(t)+δqi(t) と変換されるとき, 作用が不変であるならば, 量
δQ≡δqipi−δtH≡δqi∂q˙i∂L−δt(∂q˙i∂Lq˙i−L)は保存する(Noether の定理 Noether’s theorem):
dtdδQ=0.
Hamilton の運動方程式
Hamiltonian H は
H(qi,∂q˙i∂L,t)=q˙i∂q˙i∂L−L
で定義される. Lagrangian の場合と同様に, 最小作用の原理に対し Hamiltonian が満たす条件を求めよう. 作用を Hamiltonian H(qi,∂q˙i∂L,t)=q˙i∂q˙i∂L−L で書き直すと,
S[qi(t)]=∫t1t2dt[q˙i(t)∂q˙i(t)∂L−H(qi(t),∂q˙i(t)∂L,t)].
これを汎関数微分して,
δqj(t′)δS[qi(t)]=h→0limh1{∫t1t2dt[(q˙i(t)+hδjiδ˙(t−t′))∂(q˙i(t)+hδjiδ˙(t−t′))∂L−H(qi(t)+hδjiδ(t−t′),∂(q˙i(t)+hδjiδ˙(t−t′))∂L,t)]−∫t1t2dt[q˙i(t)∂q˙i(t)∂L−H(qi(t),∂q˙i(t)∂L,t)]}pi(t)≡∂q˙i(t)∂L,∂(q˙i(t)+hδjiδ˙(t−t′))∂L=∂q˙i(t)∂L+∂q˙k(t)∂(∂q˙i(t)∂L)hδjkδ˙(t−t′)=pi(t)+∂q˙j(t)∂pihδ˙(t−t′)=h→0limh1{∫t1t2dt[(q˙i(t)+hδjiδ˙(t−t′))(pi(t)+∂q˙j(t)∂pihδ˙(t−t′))−H(qi(t)+hδjiδ(t−t′),pi(t)+∂q˙j(t)∂pihδ˙(t−t′),t)]−∫t1t2dt[q˙i(t)∂q˙i(t)∂L−H(qi(t),∂q˙i(t)∂L,t)]}=h→0limh1∫t1t2dt[q˙i(t)∂q˙j(t)∂pihδ˙(t−t′)+pj(t)hδ˙(t−t′)−∂qj(t)∂Hhδ(t−t′)−∂pi(t)∂H∂q˙j(t)∂pihδ˙(t−t′)+O(h2)]=∫t1t2dt[q˙i(t)∂q˙j(t)∂piδ˙(t−t′)+pj(t)δ˙(t−t′)−∂qj(t)∂Hδ(t−t′)−∂pi(t)∂H∂q˙j(t)∂piδ˙(t−t′)]=∫t1t2dt[−dtd{∂q˙j(t)∂pi(q˙i(t)−∂pi(t)∂H)}−(p˙j(t)+∂qj(t)∂H)]δ(t−t′)+[{∂q˙j(t)∂pi(q˙i(t)−∂pi(t)∂H)+pj(t)}δ(t−t′)]t=t1t=t2.
したがって, 作用の変分は,
δS[qi(t)]=∫t1t2dt′δqj(t′)δS[qi(t)]δqj(t′)=∫t1t2dt′δqj(t′)∫t1t2dt[−dtd{∂q˙j(t)∂pi(q˙i(t)−∂pi(t)∂H)}−(p˙j(t)+∂qj(t)∂H)]δ(t−t′)+∫t1t2dt′δqj(t′)[{∂q˙j(t)∂pi(q˙i(t)−∂pi(t)∂H)+pj(t)}δ(t−t′)]t=t1t=t2=∫t1t2dt′∫t1t2dtδqj(t′)[−dtd{∂q˙j(t)∂pi(q˙i(t)−∂pi(t)∂H)}−(p˙j(t)+∂qj(t)∂H)]δ(t−t′)+[∫t1t2dt′δqj(t′){∂q˙j(t)∂pi(q˙i(t)−∂pi(t)∂H)+pj(t)}δ(t−t′)]t=t1t=t2=∫t1t2dtδqj(t)[−dtd{∂q˙j(t)∂pi(q˙i(t)−∂pi(t)∂H)}−(p˙j(t)+∂qj(t)∂H)]+[δqj(t){∂q˙j(t)∂pi(q˙i(t)−∂pi(t)∂H)+pj(t)}]t=t1t=t2=∫t1t2dt[δq˙j(t)∂q˙j(t)∂pi(q˙i(t)−∂pi(t)∂H)−δqj(t)(p˙j(t)+∂qj(t)∂H)]+[δqj(t)pj(t)]t=t1t=t2(δpi(t)≡δq˙j(t)∂q˙j(t)∂pi)=∫t1t2dt[δpi(t)(q˙i(t)−∂pi(t)∂H)−δqj(t)(p˙j(t)+∂qj(t)∂H)]+[δqj(t)pj(t)]t=t1t=t2.
ここで, 第2項は両端固定の境界条件 δqj(t1)=δqj(t2)=0 より消え, 第1項の被積分関数の1項目は Hamiltonian の定義から ∂H/∂pi(t)=q˙i(t) となって消えるから, 結局
δS[qi(t)]=∫t1t2dtδqj(t)(p˙j(t)+∂qj(t)∂H)
となる. δqj(t) は任意であるから, 条件 δS[qi]=0 より, t1<t<t2 で次の運動方程式が得られる:
最小作用の原理を満たすとき, Hamiltonian は以下の Hamilton の運動方程式あるいは正準方程式 canonical equation を満たす:
p˙i=−∂qi∂H,q˙i=∂pi∂H.