0

多様体上のEL方程式

21
0
$$$$

定義

多様体を$M$,曲線を$c: (\alpha, \beta) \to M,\, t \mapsto c(t)$,ラグランジアンを接束$TM$上の関数$L: TM \to \mathbb{R},\, (x, v) \mapsto L(x,v)$,作用を$S: C^{\infty}((\alpha, \beta), M) \to \mathbb{R},\, c \mapsto S[c]$とする。$C^{\infty}((\alpha, \beta), M)$は区間$(\alpha, \beta)$から多様体$M$への$C^{\infty}$級写像全体の集合,つまり曲線全体の集合である。閉区間$[a,b]$$[a,b]\subset(\alpha, \beta)$となるようにとる。作用の値$S[c]$は曲線$c$の速度ベクトルを$c'$として,
\begin{equation} S[c] = \int_{a}^{b}L(c(t), c'(t))dt \end{equation}
で与えられる。今は時間に陽に依存しないラグランジアンを考えているが依存していても以降の議論は変わらない。

多様体上の変分法

曲線の変形を$\Phi: (-\epsilon, \epsilon)\times M\to M,\,(h,q) \mapsto\Phi_{h}(q)$とする。これは$h\in(-\epsilon,\epsilon)$をひとつ決めると$M$上の点$q$を動かした点$\Phi_{h}(q)$を与える。$\Phi_{h}(c)$は変形された曲線を表す。この定義にはユークリッド空間の座標を使っていないので多様体上の変分法を表現できる
変形の条件として以下を課す。

  • $h=0$はもとの点の位置を変えない:$\Phi_{0}(q) =q$
  • 端点の固定:任意の$h$に対して$\Phi_{h}(c(a)) = c(a), \Phi_{h}(c(b)) = c(b)$
  • $\Phi$$C^{\infty}$級である。変形された曲線は自己交差を起こさない。

変形された曲線に対する作用の値は次のようになる。
\begin{equation} S[\Phi_{h}(c)] = \int_{a}^{b}L(\Phi_{h}(c(t)), \frac{\partial}{\partial t}\Phi_{h}(c(t)))dt \end{equation}

EL方程式の導出

曲線$c$を固定すると$S[\Phi_{(-)}(c)]$$h$の関数である。状況としてはひとつの決まった曲線をパラメータに応じてぐにゃぐにゃ曲げている。曲線$c$が作用の停留を与えているならば,$h=0$における微分係数が$0$である($h$$1$変数関数だから例えば2次関数を思い浮かべてもらえればよい)。よって停留作用の原理は以下で与えられる。

停留作用の原理

\begin{equation} \left.\frac{d}{dh}S[\Phi_{h}(c)] \right|_{h=0} =0 \end{equation}

多様体$M$のチャートを$(U,\phi)$,局所座標を$x^{\>i}$とする。接束$TM$に自然に誘導されるチャートの局所座標を$(x^{\>i},v^{\>i})$とする。速度ベクトルの第$i$成分$\dot{c}^{\>i}$には$\displaystyle \dot{c}^{\>i} = \frac{d}{dt}(x^{\>i} \circ c)$の関係がある。合成関数の微分を行うと,
\begin{align} &\left.\frac{d}{dh}S[\Phi_{h}(c)] \right|_{h=0} \notag \\ &= \int_{a}^{b} \sum_{i} \left.\frac{\partial L}{\partial x^{\>i}} \right|_{\substack{x = c(t) \\v = c'(t)}} \left.\frac{\partial}{\partial h}\left(x^{\>i} \circ \Phi_{h}(c(t))\right)\right|_{h=0} \notag \\ &\quad+ \left.\frac{\partial L}{\partial v^{\>i}} \right|_{\substack{x = c(t) \\v = c'(t)}} \left.\frac{\partial}{\partial h}\frac{\partial}{\partial t}\left(x^{\>i} \circ \Phi_{h}(c(t))\right)\right|_{h=0} dt \end{align}
となる。$x^{\>i}, v^{\>i}$の微分係数をとっている点$x =c(t), v=c'(t)$は,$\Phi_{0}(c(t)) = c(t)$から来ている。計算の途中ではこの表示を省く。$\sum$の中の第2項は変形$\Phi$$C^{\infty}$級であることから$h$の微分と$t$の微分を交換することができる[1]
$t$の微分を前に持ってきて部分積分を行う。もうこうなれば普通のEL方程式の導出と同じである。
\begin{align} &\int_{a}^{b}\frac{\partial L}{\partial v^{\>i}} \frac{\partial}{\partial t} \left.\frac{\partial}{\partial h}\left(x^{\>i} \circ \Phi_{h}(c(t))\right)\right|_{h=0} dt \notag \\ &= \left[ \frac{\partial L}{\partial v^{\>i}} \left.\frac{\partial}{\partial h}\left(x^{\>i} \circ \Phi_{h}(c(t))\right)\right|_{h=0} \right]_{a}^{b} \notag \\ &\quad- \int_{a}^{b}\left(\frac{d}{dt}\frac{\partial L}{\partial v^{\>i}}\right) \left.\frac{\partial}{\partial h}\left(x^{\>i} \circ \Phi_{h}(c(t))\right)\right|_{h=0}dt \notag \\ &=\int_{a}^{b} - \left(\frac{d}{dt}\frac{\partial L}{\partial v^{\>i}}\right) \left.\frac{\partial}{\partial h}\left(x^{\>i} \circ \Phi_{h}(c(t))\right)\right|_{h=0}dt \end{align}
境界項は端点においては任意の$h$で変形がない(定値写像である)ので微分係数がゼロになり消える。ここで,$V^{\>i} =\displaystyle \left.\frac{\partial}{\partial h}\left(x^{\>i} \circ \Phi_{h}(c(t))\right)\right|_{h=0} $とおく。$V^{\>i}$は変分ベクトル場というベクトル場の第$i$成分を定める(これは点$c(t)$における接ベクトルの成分である)。部分積分で得た第2項を元の積分の中に戻すと,
\begin{equation} \int_{a}^{b}\sum_{i} \left( \frac{\partial L}{\partial x^{\>i}} - \frac{d}{dt}\frac{\partial L}{\partial v^{\>i}} \right)V^{\>i}dt = 0 \end{equation}
変分法の基本補題[2]を適用することで,多様体上のEL方程式を得る。

多様体上のEL方程式

\begin{align} &\left.\frac{d}{dh}S[\Phi_{h}(c)] \right|_{h=0} =0 \notag \\ &\iff \left. \frac{d}{dt}\frac{\partial L}{\partial v^{\>i}} \right|_{\substack{x = c(t) \\v = c'(t)}} = \left. \frac{\partial L}{\partial x^{\>i}} \right|_{\substack{x = c(t) \\v = c'(t)}} \end{align}


[1]: 微分の順序交換は$C^{2}$級で十分であるがそう限る必要はないので$C^{\infty}$級を課した。

[2]: 証明したことないです。すみません。

投稿日:22日前
更新日:22日前
OptHub AI Competition

この記事を高評価した人

高評価したユーザはいません

この記事に送られたバッジ

バッジはありません。

投稿者

物理学徒です。幾何学だいすき!

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中