解析力学ノート

; draft

解析力学の基本事項をまとめる. 積極的に汎関数の計算も用いる.

最小作用の原理

古典力学において, 以下を原理として認める.

パラメータ ξiξ^i に依存する力学変数 Xα(ξi)X^α(ξ^i) に対して, 作用 action と呼ばれる汎関数 S[Xα]S[X^α] が存在し, 物理現象において力学変数 XαX^α は作用 S[Xα]S[X^α] が最小となるような関数が選ばれる. つまり, Xα(t)↦Xα(ξi)+δXα(ξi)X^α(t) ↦ X^α(ξ^i) + δX^α(ξ^i) (ただし境界固定 δXα=0δX^α=0) となる微小変換に対し, 作用が停留値を取る:

δS[Xα]≡S[Xα+δXα]−S[Xα]=0.δS[X^α] ≡ S[X^α + δX^α] - S[X^α] = 0.

この古典的原理を最小作用の原理という.

上の変分は 1 次の冪展開を用いて以下のように書き直される:

δS[Xα]=∫VdDXαδS[Xα]δXα(ξi)δXα(ξi).δS[X^α] = ∫_V \d{{}^D X^α} \fdv{S[X^α]}{X^α(ξ^i)} δX^α(ξ^i).

ただし 積分範囲 VV は力学変数 XαX^α のであり, DD はその次元である. δXα(ξi)δX^α(ξ^i) は {ξi}∈V−∂V\{ξ^i\} ∈ V - ∂V で任意だから, 最小作用の原理は汎関数微分を用いた以下の停留条件と等価である:

δS[qi(t)]δqj(t′)=0.(t1<t′<t2)\fdv{S[q^i(t)]}{q^j(t')} = 0. \quad (t_1<t'<t_2)

Euler–Lagrange の運動方程式

作用は, 座標と時間に関する Lagrangian L(qi,q˙i,t)L(q^i, \.q^i, t) を用いて以下のように表される:

S[qi]=∫t1t2dtL(qi,q˙i,t).S[q^i] = ∫_{t_1}^{t_2} \d{t} L(q^i, \.q^i, t).

最小作用の原理に対し, この Lagrangian が満たすべき条件を求めよう. qi↦qi+δqiq^i ↦ q^i + δq^i の変化に対し,

δS[qi]=∫t1t2dt[L(qi+δqi,q˙i+dδqidt,t)−L(qi,q˙i,t)]=∫t1t2dt[δqi∂L∂qi+dδqidt∂L∂q˙i]=∫t1t2dt[δqi∂L∂qi−δqiddt(∂L∂q˙i)+ddt(δqi∂L∂q˙i)]=∫t1t2dtδqi[∂L∂qi−ddt(∂L∂q˙i)]+[δqi∂L∂q˙i]t=t1t=t2.\begin{aligned} δS[q^i] =& ∫_{t_1}^{t_2} \d{t} \bqty{ L\pqty{q^i + δq^i, \.q^i + \dv{δq^i}{t}, t} - L(q^i, \.q^i, t) } \\ =& ∫_{t_1}^{t_2} \d{t} \bqty{ δq^i \pdv{L}{q^i} + \dv{δq^i}{t} \pdv{L}{\.q^i} } \\ =& ∫_{t_1}^{t_2} \d{t} \bqty{ δq^i \pdv{L}{q^i} - δq^i \dv{}{t} \pqty{\pdv{L}{\.q^i}} + \dv{}{t} \pqty{ δq^i \pdv{L}{\.q^i} } } \\ =& ∫_{t_1}^{t_2} \d{t} δq^i \bqty{ \pdv{L}{q^i} - \dv{}{t} \pqty{\pdv{L}{\.q^i}} } + \bqty{ δq^i \pdv{L}{\.q^i}}_{t=t_1}^{t=t_2}. \end{aligned}

ここで, 第2項は両端固定の境界条件 δqi(t1)=δqi(t2)=0δq^i(t_1)=δq^i(t_2)=0 より消える:

δS[qi]=∫t1t2dtδqi[∂L∂qi−ddt(∂L∂q˙i)].δS[q^i] = ∫_{t_1}^{t_2} \d{t} δq^i \bqty{ \pdv{L}{q^i} - \dv{}{t} \pqty{\pdv{L}{\.q^i}} }.

δqi(t)δq^i(t) は t1<t<t2t_1<t<t_2 で任意だから, 原理 δS[qi]=0δS[q^i] = 0 より, 運動方程式が得られる:

最小作用の原理を満たすとき, Lagrangian L(qi,q˙i,t)L(q^i,\.q^i,t) は以下の Euler–Lagrange の運動方程式を満たす:

∂L∂qi−ddt(∂L∂q˙i)=0.\pdv{L}{q^i} - \dv{}{t} \pqty{\pdv{L}{\.q^i}} = 0.

汎関数微分を用いても同様の結果が得られる. t1<t<t2t_1<t<t_2 における作用の汎関数微分は

δS[qi(t)]δqj(t′)=lim⁡h→0S[qi(t)+hδjiδ(t−t′)]−S[qi(t)]h=lim⁡h→01h[∫t1t2dtL(qi(t)+hδjiδ(t−t′),q˙i(t)+hδjiδ˙(t−t′),t)−∫t1t2dtL(qi(t),q˙i(t),t)]=lim⁡h→01h∫t1t2dt[∂L∂qi(t)hδjiδ(t−t′)+∂L∂q˙i(t)hδjiδ˙(t−t′)+o(ε2)]=∫t1t2dt[∂L∂qj(t)δ(t−t′)+∂L∂q˙j(t)δ˙(t−t′)]=∫t1t2dt[∂L∂qj(t)−ddt(∂L∂q˙j(t))]δ(t−t′)+[∂L∂q˙j(t)δ(t−t′)]t=t1t=t2.\begin{aligned} \fdv{S[q^i(t)]}{q^j(t')} &= \lim_{h → 0} \frac{S[q^i(t) + h δ^i_j δ(t - t')] - S[q^i(t)]}{h} \\ &= \lim_{h → 0} \frac1h \bqty{ ∫_{t_1}^{t_2} \d{t} L ( q^i(t) + h δ^i_j δ(t - t'), \.q^i(t) + h δ^i_j \.δ(t - t'), t) - ∫_{t_1}^{t_2} \d{t} L (q^i(t), \.q^i(t), t) } \\ &= \lim_{h → 0} \frac1h ∫_{t_1}^{t_2} \d{t} \bqty{ \pdv{L}{q^i(t)} h δ^i_j δ(t - t') + \pdv{L}{\.q^i(t)} h δ^i_j \.δ(t - t') + o(ε^2) } \\ &= ∫_{t_1}^{t_2} \d{t} \bqty{ \pdv{L}{q^j(t)} δ(t - t') + \pdv{L}{\.q^j(t)} \.δ(t - t') } \\ &= ∫_{t_1}^{t_2} \d{t} \bqty{\pdv{L}{q^j(t)} - \dv{}{t} \pqty{\pdv{L}{\.q^j(t)}}}δ(t - t') + \bqty{\pdv{L}{\.q^j(t)}δ(t-t')}_{t=t_1}^{t=t_2}. \\ \end{aligned}

したがって, t1<t′<t2t_1<t'<t_2 において,

δS[qi(t)]δqi(t′)=∂L∂qi(t′)−ddt′(∂L∂q˙i(t′)).\fdv{S[q^i(t)]}{q^i(t')} = \pdv{L}{q^i(t')} - \dv{}{t'} \pqty{\pdv{L}{\.q^i(t')}}.

または, 作用の変分を計算して,

δS[qi(t)]=∫t1t2dt′δS[qi(t)]δqj(t′)δqj(t′)=∫t1t2dt′{∫t1t2dt[∂L∂qj(t)−ddt(∂L∂q˙j(t))]δ(t−t′)+[∂L∂q˙j(t)δ(t−t′)]t=t1t=t2}δqj(t′)=∫t1t2dt′∫t1t2dtδqj(t′)[∂L∂qj(t)−ddt(∂L∂q˙j(t))]δ(t−t′)+[∫t1t2dt′δqj(t′)∂L∂q˙j(t)δ(t−t′)]t=t1t=t2=∫t1t2dtδqj(t)[∂L∂qj(t)−ddt(∂L∂q˙j(t))]+[δqj(t)∂L∂q˙j(t)]t=t1t=t2.\begin{aligned} δS[q^i(t)] &= ∫_{t_1}^{t_2} \d{t'} \fdv{S[q^i(t)]}{q^j(t')} δq^j(t') \\ &= ∫_{t_1}^{t_2} \d{t'} \qty{∫_{t_1}^{t_2} \d{t} \bqty{\pdv{L}{q^j(t)} - \dv{}{t} \pqty{\pdv{L}{\.q^j(t)}}}δ(t - t') + \bqty{\pdv{L}{\.q^j(t)}δ(t-t')}_{t=t_1}^{t=t_2}} δq^j(t') \\ &= ∫_{t_1}^{t_2} \d{t'} ∫_{t_1}^{t_2} \d{t} δq^j(t') \bqty{\pdv{L}{q^j(t)} - \dv{}{t} \pqty{\pdv{L}{\.q^j(t)}}}δ(t - t') + \bqty{∫_{t_1}^{t_2} \d{t'} δq^j(t') \pdv{L}{\.q^j(t)}δ(t-t')}_{t=t_1}^{t=t_2} \\ &= ∫_{t_1}^{t_2} \d{t} δq^j(t) \bqty{\pdv{L}{q^j(t)} - \dv{}{t} \pqty{\pdv{L}{\.q^j(t)}}} + \bqty{δq^j(t) \pdv{L}{\.q^j(t)}}_{t=t_1}^{t=t_2}. \\ \end{aligned}

このように汎関数微分を用いても同様に Euler–Lagrange の運動方程式が得られる.

例: 一次元一粒子系

一次元一粒子系の Lagrangian は,

L(q,q˙,t)=12mq˙2−V(q).L(q, \.q, t) = \frac12 m \.q^2 - V(q).

ここで,

∂L∂q=−∂V∂q,ddt(∂L∂q˙)=ddt(mq˙)=mq¨.\begin{aligned} \pdv{L}{q} &= - \pdv{V}{q}, & \dv{}{t}\pqty{\pdv{L}{\.q}} &= \dv{}{t} (m \.q) = m \"q. \end{aligned}

したがって, Euler–Lagrange の運動方程式は,

mq¨+∂V∂q=0.m\"q + \pdv{V}{q} = 0.

ポテンシャルが無い V=0V=0 ときの作用の表式を求める. 運動方程式 mq¨=0m\"q = 0 を解いて,

q˙(t)=q(t2)−q(t1)t2−t1.\.q(t) = \frac{q(t_2)-q(t_1)}{t_2-t_1}.

したがって, 作用は,

S[q]=∫t1t2dtm2(q(t2)−q(t1))2(t1−t2)2=m2(q(t2)−q(t1))2t2−t1.S[q] = ∫_{t_1}^{t_2} \d{t} \frac{m}{2} \frac{(q(t_2)-q(t_1))^2}{(t_1-t_2)^2} = \frac{m}{2} \frac{(q(t_2)-q(t_1))^2}{t_2-t_1}.

例: 調和振動子

調和振動子の Lagrangian は,

L(q,q˙,t)=12mq˙2−12mω2q2.L(q, \.q, t) = \frac12 m \.q^2 - \frac12 m ω^2 q^2.

ここで,

∂L∂q=−mω2q,ddt(∂L∂q˙)=ddt(mq˙)=mq¨.\begin{aligned} \pdv{L}{q} &= - m ω^2 q, & \dv{}{t}\pqty{\pdv{L}{\.q}} &= \dv{}{t} (m \.q) = m \"q. \end{aligned}

したがって, Euler–Lagrange の運動方程式は

mq¨+mω2q=0.m\"q + m ω^2 q = 0.

作用の表式を求める. 運動方程式を解いて,

q(t)=q1sin⁡ω(t−t2)−q2sin⁡ω(t−t1)sin⁡ω(t1−t2),q˙(t)=ωq1cos⁡ω(t−t2)−q2cos⁡ω(t−t1)sin⁡ω(t1−t2).\begin{aligned} q(t) &= \frac{q_1 \sin ω(t-t_2) - q_2 \sin ω(t-t_1)}{\sin ω(t_1-t_2)}, \\ \.q(t) &= ω \frac{q_1 \cos ω(t-t_2) - q_2 \cos ω(t-t_1)}{\sin ω(t_1-t_2)}. \\ \end{aligned}

したがって, 作用は,

S[q]=∫t1t2dtm2[{ωq1cos⁡ω(t−t2)−q2cos⁡ω(t−t1)sin⁡ω(t1−t2)}2−ω2{q1sin⁡ω(t−t2)−q2sin⁡ω(t−t1)sin⁡ω(t1−t2)}2]=∫t1t2dtmω22q12cos⁡2ω(t−t2)+q22cos⁡2ω(t−t1)−2q1q2cos⁡(2t−t1−t2)sin⁡2ω(t2−t1)=mω2sin⁡ω(t2−t1)[(q12+q22)cos⁡ω(t2−t1)−2q1q2].\begin{aligned} S[q] &= ∫_{t_1}^{t_2} \d{t} \frac{m}{2} \bqty{\qty{ω \frac{q_1 \cos ω(t-t_2) - q_2 \cos ω(t-t_1)}{\sin ω(t_1-t_2)}}^2 - ω^2 \qty{\frac{q_1 \sin ω(t-t_2) - q_2 \sin ω(t-t_1)}{\sin ω(t_1-t_2)}}^2} \\ &= ∫_{t_1}^{t_2} \d{t} \frac{mω^2}{2} \frac{q_1^2 \cos 2ω(t-t_2) + q_2^2 \cos 2ω(t-t_1) - 2q_1q_2\cos(2t-t_1-t_2)}{\sin^2 ω(t_2-t_1)} \\ &= \frac{mω}{2 \sin ω(t_2-t_1)} \bqty{(q_1^2+q_2^2) \cos ω(t_2-t_1) - 2q_1q_2}. \\ \end{aligned}

Noether の定理

時間の微小変換 t↦t′=t+δtt↦t'=t+δt に対し, 座標が qi(t)↦q′i(t′)=qi(t)+δqi(t)q^i(t)↦q'^i(t')=q^i(t)+δq^i(t) と変換されるとする. このとき t1<t<t2t_1<t<t_2 の作用は

δS[qi]=∫t1+δt(t1)t2+δt(t2)dt′L(q′i(t′),∂t′q′i(t′),t′)−∫t1t2dtL(qi(t),q˙i(t),t)(dt′=dt′dtdt=(1+δt˙)dt)=∫t1t2dt[(1+δt˙)L(q′i(t′),∂t′q′i(t′),t′)−L(qi(t),q˙i(t),t)](∂t′q′(t′)=dtdt′∂t(qi(t)+δqi(t))=(1−δt˙)(q˙i+δq˙i)=q˙i+δq˙i−q˙iδt˙)=∫t1t2dt[δt˙L+L(qi+δqi,q˙i+δq˙i−q˙iδt˙,t+δt)−L(qi,q˙i,t)]=∫t1t2dt[δt˙L+δqi∂L∂qi+(δq˙i−q˙iδt˙)∂L∂q˙i+δt∂L∂t](Lie 微分 δLqi(t)≡q′i(t)−qi(t)=δqi−q˙iδt)=∫t1t2dt[δt˙L+(δLqi+q˙iδt)∂L∂qi+(∂tδLqi+q¨iδt)∂L∂q˙i+δt∂L∂t]=∫t1t2dt[δt˙L+δLqi∂L∂qi+δtq˙i∂L∂qi+∂tδLqi∂L∂q˙i+δtq¨i∂L∂q˙i+δt∂L∂t]=∫t1t2dt[∂t(δtL)+δLqi∂L∂qi+∂t(δLqi∂L∂q˙i)−δLqiddt(∂L∂q˙i)]=∫t1t2dt{δLqi[∂L∂qi−ddt(∂L∂q˙i)]+ddt(δLqi∂L∂q˙i+δtL)}=∫t1t2dtδLqi[∂L∂qi−ddt(∂L∂q˙i)]+[δLqi∂L∂q˙i+δtL]t=t1t=t2=∫t1t2dtδLqi[∂L∂qi−ddt(∂L∂q˙i)]+[δqi∂L∂q˙i−δt(q˙i∂L∂q˙i−L)]t=t1t=t2.\begin{aligned} δS[q^i] &= ∫_{t_1+δt(t_1)}^{t_2+δt(t_2)} \d{t'} L(q'^i(t'),∂'_tq'^i(t'),t') - ∫_{t_1}^{t_2} \d{t} L(q^i(t),\.q^i(t),t) \\ & \quad \pqty{\d{t'} = \dv{t'}{t} \d{t} = (1+δ\.t) \d{t}} \\ &= ∫_{t_1}^{t_2} \d{t} \Big[ (1+δ\.t) L(q'^i(t'),∂'_tq'^i(t'),t') - L(q^i(t),\.q^i(t),t) \Big] \\ & \quad \pqty{ ∂'_tq'(t') = \dv{t}{t'} ∂_t (q^i(t)+δq^i(t)) = (1-δ\.t)(\.q^i+δ\.q^i) = \.q^i+δ\.q^i-\.q^iδ\.t } \\ &= ∫_{t_1}^{t_2} \d{t} \Big[ δ\.t L + L(q^i+δq^i,\.q^i+δ\.q^i-\.q^iδ\.t,t+δt) - L(q^i,\.q^i,t) \Big] \\ &= ∫_{t_1}^{t_2} \d{t} \bqty{δ\.t L + δq^i \pdv{L}{q^i} + (δ\.q^i-\.q^iδ\.t) \pdv{L}{\.q^i} + δt \pdv{L}{t}} \\ & \quad \pqty{\text{Lie 微分 $δ^Lq^i(t) ≡ q'^i(t) - q^i(t) = δq^i - \.q^i δt$}} \\ &= ∫_{t_1}^{t_2} \d{t} \bqty{δ\.t L + (δ^Lq^i + \.q^i δt) \pdv{L}{q^i} + (∂_tδ^Lq^i + \"q^i δt) \pdv{L}{\.q^i} + δt \pdv{L}{t}} \\ &= ∫_{t_1}^{t_2} \d{t} \bqty{δ\.t L + δ^Lq^i \pdv{L}{q^i} + δt \.q^i \pdv{L}{q^i} + ∂_tδ^Lq^i \pdv{L}{\.q^i} + δt \"q^i \pdv{L}{\.q^i} + δt \pdv{L}{t}} \\ &= ∫_{t_1}^{t_2} \d{t} \bqty{ ∂_t \pqty{δt L} + δ^Lq^i \pdv{L}{q^i} + ∂_t\pqty{δ^Lq^i \pdv{L}{\.q^i}} - δ^Lq^i \dv{}{t} \pqty{\pdv{L}{\.q^i}} } \\ &= ∫_{t_1}^{t_2} \d{t} \qty{ δ^Lq^i \bqty{\pdv{L}{q^i} - \dv{}{t} \pqty{\pdv{L}{\.q^i}}} + \dv{}{t} \pqty{δ^Lq^i \pdv{L}{\.q^i} + δt L} } \\ &= ∫_{t_1}^{t_2} \d{t} δ^Lq^i \bqty{\pdv{L}{q^i} - \dv{}{t} \pqty{\pdv{L}{\.q^i}}} + \bqty{δ^Lq^i \pdv{L}{\.q^i} + δt L}_{t=t_1}^{t=t_2} \\ &= ∫_{t_1}^{t_2} \d{t} δ^Lq^i \bqty{\pdv{L}{q^i} - \dv{}{t} \pqty{\pdv{L}{\.q^i}}} + \bqty{δq^i \pdv{L}{\.q^i} - δt \pqty{\.q^i \pdv{L}{\.q^i} - L}}_{t=t_1}^{t=t_2}. \\ \end{aligned}

ここで, 第一項は Euler–Lagrange の運動方程式より消えて, 第二項の t1t_1, t2t_2 は任意である1. したがって, この変換に対し作用が不変 δS=0δS=0 であるとすると, 対応する保存量が得られる:

時間の微小変換 t↦t′=t+δtt↦t'=t+δt に対し, 座標が qi(t)↦q′i(t′)=qi(t)+δqi(t)q^i(t)↦q'^i(t')=q^i(t)+δq^i(t) と変換されるとき, 作用が不変であるならば, 量

δQ≡δqi∂L∂q˙i−δt(∂L∂q˙iq˙i−L)δQ ≡ δq^i \pdv{L}{\.q^i} - δt \pqty{\pdv{L}{\.q^i} \.q^i - L}

は保存する(Noether の定理 Noether’s theorem):

dδQdt=0.\dv{δQ}{t} = 0.

同じ結果を汎関数微分を用いて導出しよう. 以下, 混乱を防ぐために作用の積分変数を添字で書く. つまり, St[qi]≡S[qi(t)]S_t[q^i]≡S[q^i(t)] である. 変分 δSt[qi]δS_t[q^i] を1次で展開すると2,

δSt[qi]=∫t1t2dt0δSt[qi]δqi(t0)δLqi(t0)+∫t1t2dt0δSτ(t)[qi]δτ(t0)∣τ(t)=tδt(t0).δS_t[q^i] = ∫_{t_1}^{t_2} \d{t_0} \fdv{S_t[q^i]}{q^i(t_0)} δ^Lq^i(t_0) + ∫_{t_1}^{t_2} \d{t_0} \left.\fdv{S_{τ(t)}[q^i]}{τ(t_0)}\right|_{τ(t)=t} δt(t_0).

それぞれの汎関数微分を計算すると,

δSt[qi]δqi(t0)=∫t1t2dt[∂L∂qi(t)−ddt(∂L∂q˙i(t))]δ(t−t0)+[∂L∂q˙i(t)δ(t−t0)]t=t1t=t2,\fdv{S_t[q^i]}{q^i(t_0)} = ∫_{t_1}^{t_2} \d{t} \bqty{\pdv{L}{q^i(t)} - \dv{}{t} \pqty{\pdv{L}{\.q^i(t)}}}δ(t - t_0) + \bqty{\pdv{L}{\.q^i(t)}δ(t-t_0)}_{t=t_1}^{t=t_2}, δSτ(t)[qi]δτ(t0)=lim⁡h→0Sτ(t)+hδ(t−t0)[qi]−Sτ(t)[qi]h=lim⁡h→01h{∫τ(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))}=lim⁡h→01h{∫τ(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−t0)dτ(t)−∫τ(t1)τ(t2)dτ(t)L(qi(τ(t)),q˙i(τ(t)),τ(t))}=lim⁡h→01h{∫τ(t1)τ(t2)dτ(t)[∂L∂qi(τ(t))q˙i(τ(t))hδ(t−t0)+∂L∂q˙i(τ(t))q¨i(τ(t))hδ(t−t0)+∂L∂τ(t)hδ(t−t0)]+∫τ(t1)τ(t2)dτ(t)hL(qi(τ(t)),q˙i(τ(t)),τ(t))dδ(t−t0)dτ(t)+O(h2)}=∫τ(t1)τ(t2)dτ(t){[∂L∂qi(τ(t))q˙i(τ(t))+∂L∂q˙i(τ(t))q¨i(τ(t))+∂L∂τ(t)]δ(t−t0)+L(qi(τ(t)),q˙i(τ(t)),τ(t))dδ(t−t0)dτ(t)}=∫τ(t1)τ(t2)dτ(t)ddτ(t)[L(qi(τ(t)),q˙i(τ(t)),τ(t))δ(t−t0)]=[L(qi(τ(t)),q˙i(τ(t)),τ(t))δ(t−t0)]τ(t)=τ(t1)τ(t)=τ(t2)\begin{aligned} \fdv{S_{τ(t)}[q^i]}{τ(t_0)} &= \lim_{h→0} \frac{S_{τ(t)+hδ(t-t_0)}[q^i] - S_{τ(t)}[q^i]}h \\ &= \lim_{h→0} \frac1h \left\{∫_{τ(t_1)}^{τ(t_2)} \d{(τ(t)+hδ(t-t_0))} L(q^i(τ(t)+hδ(t-t_0)),\.q^i(τ(t)+hδ(t-t_0)),τ(t)+hδ(t-t_0))\right. \\ &\qquad\qquad\quad - \left.∫_{τ(t_1)}^{τ(t_2)} \d{τ(t)} L(q^i(τ(t)),\.q^i(τ(t)),τ(t))\right\} \\ &= \lim_{h→0} \frac1h \left\{∫_{τ(t_1)}^{τ(t_2)} \d{τ(t)} L(q^i(τ(t)+hδ(t-t_0)),\.q^i(τ(t)+hδ(t-t_0)),τ(t)+hδ(t-t_0))\right. \\ &\qquad\qquad\quad + ∫_{τ(t_1)}^{τ(t_2)} \d{τ(t)} h L(q^i(τ(t)+hδ(t-t_0)),\.q^i(τ(t)+hδ(t-t_0)),τ(t)+hδ(t-t_0)) \dv{δ(t-t_0)}{τ(t)} \\ &\qquad\qquad\quad - \left.∫_{τ(t_1)}^{τ(t_2)} \d{τ(t)} L(q^i(τ(t)),\.q^i(τ(t)),τ(t))\right\} \\ &= \lim_{h→0} \frac1h \left\{∫_{τ(t_1)}^{τ(t_2)} \d{τ(t)} \bqty{\pdv{L}{q^i(τ(t))} \.q^i(τ(t)) hδ(t-t_0) + \pdv{L}{\.q^i(τ(t))} \"q^i(τ(t)) hδ(t-t_0) + \pdv{L}{τ(t)} hδ(t-t_0)}\right. \\ &\qquad\qquad\quad + \left.∫_{τ(t_1)}^{τ(t_2)} \d{τ(t)} h L(q^i(τ(t)),\.q^i(τ(t)),τ(t)) \dv{δ(t-t_0)}{τ(t)} + O(h^2) \right\} \\ &= ∫_{τ(t_1)}^{τ(t_2)} \d{τ(t)} \qty{\bqty{\pdv{L}{q^i(τ(t))} \.q^i(τ(t)) + \pdv{L}{\.q^i(τ(t))} \"q^i(τ(t)) + \pdv{L}{τ(t)}} δ(t-t_0) + L(q^i(τ(t)),\.q^i(τ(t)),τ(t)) \dv{δ(t-t_0)}{τ(t)}} \\ &= ∫_{τ(t_1)}^{τ(t_2)} \d{τ(t)} \dv{}{τ(t)} \bqty{L(q^i(τ(t)),\.q^i(τ(t)),τ(t)) δ(t-t_0)} \\ &= \bqty{L(q^i(τ(t)),\.q^i(τ(t)),τ(t)) δ(t-t_0)}_{τ(t)=τ(t_1)}^{τ(t)=τ(t_2)} \\ \end{aligned} ∴δSτ(t)[qi]δτ(t0)∣τ(t)=t=[L(qi(t),q˙i(t),t)δ(t−t0)]t=t1t=t2.∴ \left.\fdv{S_{τ(t)}[q^i]}{τ(t_0)}\right|_{τ(t)=t} = \bqty{L(q^i(t),\.q^i(t),t) δ(t-t_0)}_{t=t_1}^{t=t_2}. \\

これらを変分の式に代入すると,

δSt[qi]=∫t1t2dt0δLqi(t0){∫t1t2dt[∂L∂qi(t)−ddt(∂L∂q˙i(t))]δ(t−t0)+[∂L∂q˙i(t)δ(t−t0)]t=t1t=t2}+∫t1t2dt0δt(t0)[L(qi(t),q˙i(t),t)δ(t−t0)]t=t1t=t2=∫t1t2dt0∫t1t2dtδLqi(t0)[∂L∂qi(t)−ddt(∂L∂q˙i(t))]δ(t−t0)+[∫t1t2dt0δLqi(t0)∂L∂q˙i(t)δ(t−t0)]t=t1t=t2+[∫t1t2dt0δt(t0)L(qi(t),q˙i(t),t)δ(t−t0)]t=t1t=t2=∫t1t2dtδLqi(t)[∂L∂qi(t)−ddt(∂L∂q˙i(t))]+[δLqi(t)∂L∂q˙i(t)+δt(t)L(qi(t),q˙i(t),t)]t=t1t=t2=∫t1t2dtδLqi(t)[∂L∂qi(t)−ddt(∂L∂q˙i(t))]+[δqi(t)∂L∂q˙i(t)−δt(q˙i(t)∂L∂q˙i(t)−L(qi(t),q˙i(t),t))]t=t1t=t2.\begin{aligned} δS_t[q^i] &= ∫_{t_1}^{t_2} \d{t_0} δ^Lq^i(t_0) \qty{∫_{t_1}^{t_2} \d{t} \bqty{\pdv{L}{q^i(t)} - \dv{}{t} \pqty{\pdv{L}{\.q^i(t)}}}δ(t - t_0) + \bqty{\pdv{L}{\.q^i(t)}δ(t-t_0)}_{t=t_1}^{t=t_2}} \\ &\quad\quad+ ∫_{t_1}^{t_2} \d{t_0} δt(t_0) \bqty{L(q^i(t),\.q^i(t),t) δ(t-t_0)}_{t=t_1}^{t=t_2} \\ &= ∫_{t_1}^{t_2} \d{t_0} ∫_{t_1}^{t_2} \d{t} δ^Lq^i(t_0) \bqty{\pdv{L}{q^i(t)} - \dv{}{t} \pqty{\pdv{L}{\.q^i(t)}}}δ(t - t_0) + \bqty{∫_{t_1}^{t_2} \d{t_0} δ^Lq^i(t_0) \pdv{L}{\.q^i(t)}δ(t-t_0)}_{t=t_1}^{t=t_2} \\ &\quad\quad+ \bqty{∫_{t_1}^{t_2} \d{t_0} δt(t_0) L(q^i(t),\.q^i(t),t) δ(t-t_0)}_{t=t_1}^{t=t_2} \\ &= ∫_{t_1}^{t_2} \d{t} δ^Lq^i(t) \bqty{\pdv{L}{q^i(t)} - \dv{}{t} \pqty{\pdv{L}{\.q^i(t)}}} + \bqty{δ^Lq^i(t) \pdv{L}{\.q^i(t)} + δt(t) L(q^i(t),\.q^i(t),t)}_{t=t_1}^{t=t_2} \\ &= ∫_{t_1}^{t_2} \d{t} δ^Lq^i(t) \bqty{\pdv{L}{q^i(t)} - \dv{}{t} \pqty{\pdv{L}{\.q^i(t)}}} + \bqty{δq^i(t) \pdv{L}{\.q^i(t)} - δt \pqty{\.q^i(t) \pdv{L}{\.q^i(t)} - L(q^i(t),\.q^i(t),t)}}_{t=t_1}^{t=t_2}. \\ \end{aligned}

したがって同様に示された.

例: 空間並進に対する不変量

時間並進 t↦t′=t,qi(t)↦q′i(t′)=qi(t)+εit↦t'=t, q^i(t)↦q'^i(t')=q^i(t)+ε^i に対し, 作用が不変であるとき, 対応する保存量は

δQ=εi∂L∂q˙i=const.∴∂L∂q˙i=const.\begin{gathered} δQ = ε^i \pdv{L}{\.q^i} = \mathrm{const.} \\ ∴ \pdv{L}{\.q^i} = \mathrm{const.} \end{gathered}

この不変量 pi≡∂L∂q˙i\displaystyle p_i ≡ \pdv{L}{\.q^i} は一般化運動量と呼ばれる.

例: 時間並進に対する不変量

時間並進 t↦t′=t+εt↦t'=t+ε, qi(t)↦q′i(t′)=qi(t)q^i(t)↦q'^i(t')=q^i(t) に対し, 作用が不変であるとき, 対応する保存量は

δQ=−ε(q˙i∂L∂q˙i−L)=const.∴q˙i∂L∂q˙i−L=const.\begin{gathered} δQ = - ε \pqty{\.q^i \pdv{L}{\.q^i} - L} = \mathrm{const.} \\ ∴ \.q^i \pdv{L}{\.q^i} - L = \mathrm{const.} \end{gathered}

この不変量 H≡q˙i∂L∂q˙i−L=q˙ipi−L\displaystyle H ≡ \.q^i \pdv{L}{\.q^i} - L = \.q^i p_i - L は Hamiltonian と呼ばれる.

例: 空間回転に対する不変量

3 次元空間での一粒子 q=xq=\bm{x} を考える. 空間回転 t↦t′=t,x(t)↦x′(t′)=R(ε)x(t)=x(t)−ε×x(t)t↦t'=t, \bm{x}(t) ↦ \bm{x}'(t') = R(\bm{ε}) \bm{x}(t) = \bm{x}(t) - \bm{ε} × \bm{x}(t) に対し, 作用が不変であるとき, 対応する保存量は

δQ=(−ε×x)⋅∂L∂x˙=−ε⋅(x×∂L∂x˙)=const.δQ = (- \bm{ε} × \bm{x}) ⋅ \pdv{L}{\.{\bm{x}}} = - \bm{ε} ⋅ \pqty{\bm{x} × \pdv{L}{\.{\bm{x}}}} = \mathrm{const.} ∴x×∂L∂x˙=const.∴ \bm{x} × \pdv{L}{\.{\bm{x}}} = \mathrm{const.}

この不変量 l≡x×∂L∂x˙=x×p\displaystyle \bm{l} ≡ \bm{x} × \pdv{L}{\.{\bm{x}}} = \bm{x} × \bm{p} は角運動量と呼ばれる.

一般化運動量と Hamiltonian

Lagrangian LL が与えられたとき, qiq^i に対して

pi≡∂L∂q˙ip_i ≡ \pdv{L}{\.q^i}

を一般化運動量, または qiq^i に共役な運動量 conjugate momentum といい, 一般化座標とそれに共役な運動量の組 (qi,pi)(q^i, p_i) を正準変数 canonical variables という.

Lagrangian LL と正準変数 (qi,pi)(q^i, p_i) が与えられたとき,

H(qi,pi,t)≡q˙ipi−LH(q^i, p_i, t) ≡ \.q^i p_i - L

を Hamiltonian という.

一般化運動量と Hamiltonian は作用を端点で偏微分することで得ることもできる:

pi(t)=∂S∂qi(t),H(qi,pi,t)=−∂S∂t.p_i(t) = \pdv{S}{q^i(t)}, \quad H(q^i,p_i,t) = - \pdv{S}{t}.

ただし作用は S[q]=∫t0tdt′L(qi,q˙i,t′)S[q]=∫_{t_0}^{t} \d{t'} L(q^i,\.q^i,t') で与えられている. 実際, Norther の定理と同じ状況での変分は

δS[qi]=[δqi∂L∂q˙i−δt(q˙i∂L∂q˙i−L)]t′=t0t′=t=[δqipi−δtH]t′=t0t′=t.δS[q^i] = \bqty{δq^i \pdv{L}{\.q^i} - δt \pqty{\.q^i \pdv{L}{\.q^i} - L}}_{t'=t_0}^{t'=t} = \bqty{δq^i p_i - δt H}_{t'=t_0}^{t'=t}.

始点での変位を δt(t0)=δqi(t0)=0δt(t_0)=δq^i(t_0)=0 とすれば,

δS[qi]=δqipi−δtH.δS[q^i] = δq^i p_i - δt H.

この変分は経路の途中 t′∈(t0,t)t'∈(t_0,t) によらない形になっているから, 一点 tt での変位から求めたい全微分が得られる:

dS=dqipi−dtH.\d{S} = \d{q^i} p_i - \d{t} H.

これらの性質を組み合わせることで以下の方程式が得られる:

最小作用の原理を満たす作用 S[qi]=∫t0tdt′L(qi,q˙i,t′)S[q^i] = ∫_{t_0}^t \d{t'} L(q^i,\.q^i,t') に対し, 作用の端点 tt, q(t)q(t) での偏微分は Hamilton–Jacobi 方程式 Hamilton–Jacobi equation を満たす:

H(qi(t),∂S∂qi(t),t)+∂S∂t=0.H\pqty{q^i(t),\pdv{S}{q^i(t)},t}+\pdv{S}{t}=0.

Hamilton の運動方程式

Lagrangian の場合と同様に, 最小作用の原理に対し Hamiltonian が満たす条件を求めよう. Hamiltonian H(qi,pi,t)≡q˙ipi−LH(q^i, p_i, t) ≡ \.q^i p_i - L の全微分は,

dH=q˙idpi+pidq˙i−dL=q˙idpi+pidq˙i−∂L∂qidqi−pidq˙i−∂L∂tdt(∵dL=∂L∂qidqi+∂L∂q˙idq˙i+∂L∂tdt)=−∂L∂qidqi+q˙idpi−∂L∂tdt.\begin{aligned} \d{H} &= \.q^i \d{p_i} + p_i \d{\.q^i} - \d{L} \\ &= \.q^i \d{p_i} + p_i \d{\.q^i} - \pdv{L}{q^i} \d{q^i} - p_i \d{\.q^i} - \pdv{L}{t} \d{t} \\ & \quad \pqty{ ∵ \d{L} = \pdv{L}{q^i} \d{q^i} + \pdv{L}{\.q^i} \d{\.q^i} + \pdv{L}{t} \d{t} } \\ &= - \pdv{L}{q^i} \d{q^i} + \.q^i \d{p_i} - \pdv{L}{t} \d{t}. \end{aligned}

ここで, Euler-Lagrangian 方程式が成立するとき p˙i=∂L/∂qi\.p_i = ∂L / ∂q^i であることを用いると, Hamiltonian に関する運動方程式が得られる:

最小作用の原理を満たすとき, Hamiltonian は以下の Hamilton の運動方程式あるいは正準方程式 canonical equation を満たす:

p˙i=−∂H∂qi,q˙i=∂H∂pi.\.p_i = - \pdv{H}{q^i}, \quad \.q^i = \pdv{H}{p_i}.

Lagrangian が時間に陽に依存しないとき, Hamiltonian は保存する:

∂H∂t=−∂L∂t=0.\pdv{H}{t} = -\pdv{L}{t} = 0.

同じ結果を汎関数微分を用いて導出しよう. 作用を Hamiltonian H(qi,∂L∂q˙i,t)=q˙i∂L∂q˙i−L\displaystyle H\pqty{q^i,\pdv{L}{\.q^i},t} = \.q^i\pdv{L}{\.q^i} - L で書き直すと,

S[qi(t)]=∫t1t2dt[q˙i(t)∂L∂q˙i(t)−H(qi(t),∂L∂q˙i(t),t)].S[q^i(t)] = ∫_{t_1}^{t_2}\d{t} \bqty{\.q^i(t) \pdv{L}{\.q^i(t)} - H\pqty{q^i(t),\pdv{L}{\.q^i(t)},t}}.

これを汎関数微分して,

δS[qi(t)]δqj(t′)=lim⁡h→01h{∫t1t2dt[(q˙i(t)+hδjiδ˙(t−t′))∂L∂(q˙i(t)+hδjiδ˙(t−t′))−H(qi(t)+hδjiδ(t−t′),∂L∂(q˙i(t)+hδjiδ˙(t−t′)),t)]−∫t1t2dt[q˙i(t)∂L∂q˙i(t)−H(qi(t),∂L∂q˙i(t),t)]}(pi(t)≡∂L∂q˙i(t),∂L∂(q˙i(t)+hδjiδ˙(t−t′))=∂L∂q˙i(t)+∂∂q˙k(t)(∂L∂q˙i(t))hδjkδ˙(t−t′)=pi(t)+∂pi∂q˙j(t)hδ˙(t−t′))=lim⁡h→01h{∫t1t2dt[(q˙i(t)+hδjiδ˙(t−t′))(pi(t)+∂pi∂q˙j(t)hδ˙(t−t′))−H(qi(t)+hδjiδ(t−t′),pi(t)+∂pi∂q˙j(t)hδ˙(t−t′),t)]−∫t1t2dt[q˙i(t)∂L∂q˙i(t)−H(qi(t),∂L∂q˙i(t),t)]}=lim⁡h→01h∫t1t2dt[q˙i(t)∂pi∂q˙j(t)hδ˙(t−t′)+pj(t)hδ˙(t−t′)−∂H∂qj(t)hδ(t−t′)−∂H∂pi(t)∂pi∂q˙j(t)hδ˙(t−t′)+O(h2)]=∫t1t2dt[q˙i(t)∂pi∂q˙j(t)δ˙(t−t′)+pj(t)δ˙(t−t′)−∂H∂qj(t)δ(t−t′)−∂H∂pi(t)∂pi∂q˙j(t)δ˙(t−t′)]=∫t1t2dt[−ddt{∂pi∂q˙j(t)(q˙i(t)−∂H∂pi(t))}−(p˙j(t)+∂H∂qj(t))]δ(t−t′)+[{∂pi∂q˙j(t)(q˙i(t)−∂H∂pi(t))+pj(t)}δ(t−t′)]t=t1t=t2.\begin{aligned} \fdv{S[q^i(t)]}{q^j(t')} &= \lim_{h→0} \frac1h \left\{∫_{t_1}^{t_2}\d{t} \bqty{(\.q^i(t)+hδ^i_j\.δ(t-t')) \pdv{L}{(\.q^i(t)+hδ^i_j\.δ(t-t'))} - H\pqty{q^i(t)+hδ^i_jδ(t-t'),\pdv{L}{(\.q^i(t)+hδ^i_j\.δ(t-t'))},t}}\right. \\ &\qquad\qquad \left. - ∫_{t_1}^{t_2}\d{t} \bqty{\.q^i(t) \pdv{L}{\.q^i(t)} - H\pqty{q^i(t),\pdv{L}{\.q^i(t)},t}}\right\} \\ &\quad \pqty{\begin{aligned} p_i(t) ≡ \pdv{L}{\.q^i(t)}, \quad \pdv{L}{(\.q^i(t)+hδ^i_j\.δ(t-t'))} &= \pdv{L}{\.q^i(t)} + \pdv{}{\.q^{k}(t)}\pqty{\pdv{L}{\.q^{i}(t)}} hδ^k_j\.δ(t-t') \\ &= p_i(t)+\pdv{p_i}{\.q^{j}(t)} h\.δ(t-t') \end{aligned}} \\ &= \lim_{h→0} \frac1h \left\{∫_{t_1}^{t_2}\d{t} \bqty{(\.q^i(t)+hδ^i_j\.δ(t-t')) \pqty{p_i(t)+\pdv{p_i}{\.q^{j}(t)} h\.δ(t-t')} - H\pqty{q^i(t)+hδ^i_jδ(t-t'),p_i(t)+\pdv{p_i}{\.q^{j}(t)} h\.δ(t-t'),t}}\right. \\ &\qquad\qquad \left. - ∫_{t_1}^{t_2}\d{t} \bqty{\.q^i(t) \pdv{L}{\.q^i(t)} - H\pqty{q^i(t),\pdv{L}{\.q^i(t)},t}}\right\} \\ &= \lim_{h→0} \frac1h ∫_{t_1}^{t_2}\d{t} \bqty{\.q^i(t) \pdv{p_i}{\.q^{j}(t)} h\.δ(t-t') + p_j(t) h\.δ(t-t') - \pdv{H}{q^j(t)} hδ(t-t') - \pdv{H}{p_i(t)} \pdv{p_i}{\.q^{j}(t)} h\.δ(t-t') + O(h^2)} \\ &= ∫_{t_1}^{t_2}\d{t} \bqty{\.q^i(t) \pdv{p_i}{\.q^{j}(t)} \.δ(t-t') + p_j(t) \.δ(t-t') - \pdv{H}{q^j(t)} δ(t-t') - \pdv{H}{p_i(t)} \pdv{p_i}{\.q^{j}(t)} \.δ(t-t') } \\ &= ∫_{t_1}^{t_2}\d{t} \bqty{- \dv{}{t} \qty{\pdv{p_i}{\.q^{j}(t)} \pqty{\.q^i(t) - \pdv{H}{p_i(t)}}} - \pqty{\.p_j(t) + \pdv{H}{q^j(t)}}} δ(t-t') \\ &\quad + \bqty{\qty{\pdv{p_i}{\.q^{j}(t)} \pqty{\.q^i(t) - \pdv{H}{p_i(t)}} + p_j(t)} δ(t-t')}_{t=t_1}^{t=t_2}. \end{aligned}

したがって, 作用の変分は,

δS[qi(t)]=∫t1t2dt′δS[qi(t)]δqj(t′)δqj(t′)=∫t1t2dt′δqj(t′)∫t1t2dt[−ddt{∂pi∂q˙j(t)(q˙i(t)−∂H∂pi(t))}−(p˙j(t)+∂H∂qj(t))]δ(t−t′)+∫t1t2dt′δqj(t′)[{∂pi∂q˙j(t)(q˙i(t)−∂H∂pi(t))+pj(t)}δ(t−t′)]t=t1t=t2=∫t1t2dt′∫t1t2dtδqj(t′)[−ddt{∂pi∂q˙j(t)(q˙i(t)−∂H∂pi(t))}−(p˙j(t)+∂H∂qj(t))]δ(t−t′)+[∫t1t2dt′δqj(t′){∂pi∂q˙j(t)(q˙i(t)−∂H∂pi(t))+pj(t)}δ(t−t′)]t=t1t=t2=∫t1t2dtδqj(t)[−ddt{∂pi∂q˙j(t)(q˙i(t)−∂H∂pi(t))}−(p˙j(t)+∂H∂qj(t))]+[δqj(t){∂pi∂q˙j(t)(q˙i(t)−∂H∂pi(t))+pj(t)}]t=t1t=t2=∫t1t2dt[δq˙j(t)∂pi∂q˙j(t)(q˙i(t)−∂H∂pi(t))−δqj(t)(p˙j(t)+∂H∂qj(t))]+[δqj(t)pj(t)]t=t1t=t2(δpi(t)≡δq˙j(t)∂pi∂q˙j(t))=∫t1t2dt[δpi(t)(q˙i(t)−∂H∂pi(t))−δqj(t)(p˙j(t)+∂H∂qj(t))]+[δqj(t)pj(t)]t=t1t=t2.\begin{aligned} δS[q^i(t)] &= ∫_{t_1}^{t_2} \d{t'} \fdv{S[q^i(t)]}{q^j(t')} δq^j(t') \\ &= ∫_{t_1}^{t_2} \d{t'} δq^j(t') ∫_{t_1}^{t_2}\d{t} \bqty{- \dv{}{t} \qty{\pdv{p_i}{\.q^{j}(t)} \pqty{\.q^i(t) - \pdv{H}{p_i(t)}}} - \pqty{\.p_j(t) + \pdv{H}{q^j(t)}}} δ(t-t') \\ &\quad + ∫_{t_1}^{t_2} \d{t'} δq^j(t') \bqty{\qty{\pdv{p_i}{\.q^{j}(t)} \pqty{\.q^i(t) - \pdv{H}{p_i(t)}} + p_j(t)} δ(t-t')}_{t=t_1}^{t=t_2} \\ &= ∫_{t_1}^{t_2} \d{t'} ∫_{t_1}^{t_2}\d{t} δq^j(t') \bqty{- \dv{}{t} \qty{\pdv{p_i}{\.q^{j}(t)} \pqty{\.q^i(t) - \pdv{H}{p_i(t)}}} - \pqty{\.p_j(t) + \pdv{H}{q^j(t)}}} δ(t-t') \\ &\quad + \bqty{∫_{t_1}^{t_2} \d{t'} δq^j(t') \qty{\pdv{p_i}{\.q^{j}(t)} \pqty{\.q^i(t) - \pdv{H}{p_i(t)}} + p_j(t)} δ(t-t')}_{t=t_1}^{t=t_2} \\ &= ∫_{t_1}^{t_2}\d{t} δq^j(t) \bqty{- \dv{}{t} \qty{\pdv{p_i}{\.q^{j}(t)} \pqty{\.q^i(t) - \pdv{H}{p_i(t)}}} - \pqty{\.p_j(t) + \pdv{H}{q^j(t)}}} \\ &\quad + \bqty{δq^j(t) \qty{\pdv{p_i}{\.q^{j}(t)} \pqty{\.q^i(t) - \pdv{H}{p_i(t)}} + p_j(t)}}_{t=t_1}^{t=t_2} \\ &= ∫_{t_1}^{t_2}\d{t} \bqty{δ\.q^j(t) \pdv{p_i}{\.q^{j}(t)} \pqty{\.q^i(t) - \pdv{H}{p_i(t)}} - δq^j(t) \pqty{\.p_j(t) + \pdv{H}{q^j(t)}}} + \bqty{δq^j(t) p_j(t)}_{t=t_1}^{t=t_2} \\ &\quad \pqty{δp_i(t) ≡ δ\.q^j(t) \pdv{p_i}{\.q^{j}(t)}} \\ &= ∫_{t_1}^{t_2}\d{t} \bqty{δp_i(t) \pqty{\.q^i(t) - \pdv{H}{p_i(t)}} - δq^j(t) \pqty{\.p_j(t) + \pdv{H}{q^j(t)}}} + \bqty{δq^j(t) p_j(t)}_{t=t_1}^{t=t_2}. \\ \end{aligned}

ここで, 第2項は両端固定の境界条件 δqj(t1)=δqj(t2)=0δq^j(t_1)=δq^j(t_2)=0 より消える:

δS[qi(t)]=∫t1t2dt[δpi(t)(q˙i(t)−∂H∂pi(t))−δqj(t)(p˙j(t)+∂H∂qj(t))].δS[q^i(t)] = ∫_{t_1}^{t_2}\d{t} \bqty{δp_i(t) \pqty{\.q^i(t) - \pdv{H}{p_i(t)}} - δq^j(t) \pqty{\.p_j(t) + \pdv{H}{q^j(t)}}}.

また, δqj(t)δq^j(t), δq˙j(t)δ\.q^j(t) は任意であるから, δqj(t)δq^j(t) と δpj(t)δp_j(t) を独立に取ることができる. したがって, 条件 δS[qi]=0δS[q^i] = 0 より, t1<t<t2t_1<t<t_2 で正準方程式が得られる.

qiq^i と pip_i の独立性を陽に表した作用

S[qi(t),pi(t)]=∫t1t2dt[q˙i(t)pi(t)−H(qi(t),pi(t),t)].S[q^i(t), p_i(t)] = ∫_{t_1}^{t_2}\d{t} \bqty{\.q^i(t) p_i(t) - H\pqty{q^i(t),p_i(t),t}}.

も用いられる.

例: 一次元一粒子系

一次元一粒子系の Lagrangian は,

L(q,q˙,t)=12mq˙2−V(q).L(q, \.q, t) = \frac12 m \.q^2 - V(q).

ここで, 一般化運動量の定義より,

p=∂L∂q˙=mq˙.p = \pdv{L}{\.q} = m \.q.

したがって q˙=p/m\.q = p / m であるから, Hamiltonian の定義より,

H(q,p,t)=pmp−L(q,pm,t)=p22m+V(q).H(q, p, t) = \frac{p}{m} p - L\pqty{q, \frac{p}{m}, t} = \frac{p^2}{2m} + V(q).

ここで,

∂H∂q=dVdq,∂H∂p=pm.\begin{aligned} \pdv{H}{q} &= \dv{V}{q}, & \pdv{H}{p} &= \frac{p}{m}. \end{aligned}

したがって, Hamilton の運動方程式は,

p˙=−dVdq,q˙=pm.\begin{aligned} \.p &= - \dv{V}{q}, & \.q &= \frac{p}{m}. \end{aligned}

例: 調和振動子

調和振動子の Lagrangian は,

L(q,q˙,t)=12mq˙2−12mω2q2.L(q, \.q, t) = \frac12 m \.q^2 - \frac12 m ω^2 q^2.

ここで, 一般化運動量の定義より,

p=∂L∂q˙=mq˙.p = \pdv{L}{\.q} = m \.q.

したがって q˙=p/m\.q = p / m であるから, Hamiltonian の定義より,

H(q,p,t)=pmp−L(q,pm,t)=p22m+12mω2q2.H(q, p, t) = \frac{p}{m} p - L\pqty{q, \frac{p}{m}, t} = \frac{p^2}{2m} + \frac12 m ω^2 q^2.

ここで,

∂H∂q=mω2q,∂H∂p=pm.\begin{aligned} \pdv{H}{q} &= m ω^2 q, & \pdv{H}{p} &= \frac{p}{m}. \end{aligned}

したがって, Hamilton の運動方程式は,

p˙=−mω2q,q˙=pm.\begin{aligned} \.p &= - m ω^2 q, & \.q &= \frac{p}{m}. \end{aligned}

正準変換

正準変数の変換 (qi,pi)↦(Qj,Pj)=(Qj(qi,pi),Pj(qi,pi))(q^i, p_i) ↦ (Q^j, P_j) = (Q^j(q^i, p_i), P_j(q^i, p_i)) に対して Hamiltonian が H(qi,pi)↦K(Qj,Pj)H (q^i, p_i) ↦ K (Q^j, P_j) と変換されるとき, この正準変数の変換を正準変換 canonical transformation という. Hamiltonian の定義から, δS[qi,pi]=δ∫dt(q˙ipi−H)=0δS[q^i,p_i] = δ∫\d{t} (\.q^i p_i - H) = 0 かつ δS′[Qi,Pi]=δ∫dt(Q˙iPi−K)=0δS'[Q^i,P_i] = δ∫\d{t} (\.Q^i P_i - K) = 0. したがって, ある関数 WW が存在して,

(q˙ipi−H)−(Q˙iPi−K)=dWdt.∴dW=pidqi−PidQi−(H−K)dt.\begin{gathered} (\.q^i p_i - H) - (\.Q^i P_i - K) = \dv{W}{t}. \\ ∴\d{W} = p_i \d{q^i} - P_i \d{Q^i} - (H - K) \d{t}. \end{gathered}

または, 両辺に dQiPi/dt\d{Q^i P_i}/\d{t} を足して,

(q˙ipi−H)−(−QiP˙i−K)=ddt(W+QiPi)=:dW′dt.∴dW′=pidqi+QidPi−(H−K)dt.\begin{gathered} (\.q^i p_i - H) - (- Q^i \.P_i - K) = \dv{}{t} \pqty{W + Q^i P_i} =: \dv{W'}{t}. \\ ∴\d{W'} = p_i \d{q^i} + Q^i \d{P_i} - (H - K) \d{t}. \end{gathered}

これら W(qi,Qi,t)W(q^i, Q^i, t), W′(qi,Pi,t)W'(q^i, P_i, t) をどちらも母関数といい, 以下を満たす.

pi=∂W∂qi,Pi=−∂W∂Qi,K=H+∂W∂t,pi=∂W′∂qi,Qi=∂W′∂Pi,K=H+∂W′∂t.\begin{gathered} p_i = \pdv{W}{q^i}, \quad P_i = - \pdv{W}{Q^i}, \quad K = H + \pdv{W}{t}, \\ p_i = \pdv{W'}{q^i}, \quad Q^i = \pdv{W'}{P_i}, \quad K = H + \pdv{W'}{t}. \end{gathered}

Poisson 括弧

正準変数 (qi,pi)(q^i, p_i) に対し, Poisson 括弧 Poisson braket は以下で定義される演算である:

{A,B}P≡∂A∂qi∂B∂pi−∂B∂qi∂A∂pi.\{A, B\}_\mathrm{P} ≡ \pdv{A}{q^i}\pdv{B}{p_i} - \pdv{B}{q^i}\pdv{A}{p_i}.

正準変数自身は以下を満たす:

{qi,qj}P={pi,pj}P=0,{qi,pj}P=δji.\begin{aligned} \{q^i, q^j\}_\mathrm{P} = \{p_i, p_j\}_\mathrm{P} = 0, && \{q^i, p_j\}_\mathrm{P} = δ_j^i. \end{aligned}

また, Hamilton の運動方程式は以下のように書き換えられる:

{qi,H}P=dqidt,{pi,H}P=dpidt.\begin{aligned} \{q^i, H\}_\mathrm{P} = \dv{q^i}{t}, && \{p_i, H\}_\mathrm{P} = \dv{p_i}{t}. \end{aligned}

より一般に, 正準変数と時間に関する物理量 A(qi,pi,t)A(q^i, p_i, t) について, 時間微分に関して以下が成立する:

dAdt={A,H}P+∂A∂t.\dv{A}{t} = \{A, H\}_\mathrm{P} + \pdv{A}{t}.

実際, AA の時間による完全微分は,

dAdt=∂A∂qiq˙i+∂A∂pip˙i+∂A∂t=∂A∂qi∂H∂pi−∂A∂pi∂H∂qi+∂A∂t={A,H}P+∂A∂t.\begin{aligned} \dv{A}{t} &= \pdv{A}{q^i} \.q^i + \pdv{A}{p_i} \.p_i + \pdv{A}{t} \\ &= \pdv{A}{q^i} \pdv{H}{p_i} - \pdv{A}{p_i} \pdv{H}{q^i} + \pdv{A}{t} \\ &= \{A, H\}_\mathrm{P} + \pdv{A}{t}. \end{aligned}

参考文献

Footnotes

  1. 最小作用の原理の場合と違い, このときの δqiδq^i は両端固定でない. そのため, Euler-Lagrange の運動方程式の際に消えた発散項を, 今回の場合は消すことができない. ↩

  2. 積分変数を任意の関数 τ(t)τ(t) にしたときの作用 Sτ(t)[qi]S_{τ(t)}[q^i] は

    δSτ(t)[qi]=∫τ(t1)τ(t2)dτ(t0)δSτ(t)[qi]δqi(τ(t0))δLqi(τ(t0))+∫t1t2dt0δSτ(t)[qi]δτ(t0)δτ(t0).δS_τ(t)[q^i] = ∫_{τ(t_1)}^{τ(t_2)} \d{τ(t_0)} \fdv{S_{τ(t)}[q^i]}{q^i(τ(t_0))} δ^Lq^i(τ(t_0)) + ∫_{t_1}^{t_2} \d{t_0} \fdv{S_{τ(t)}[q^i]}{τ(t_0)} δτ(t_0).

    で与えられる. このとき τ(t)=tτ(t)=t, δτ(t)=δt(t)δτ(t) = δt(t) とすれば文中の式が得られる. ↩