第1章―ラグランジュ力学を復習する
こんにちは、管理人のturtleです。
本稿はあくまで管理人の備忘録として書いておりますが、解析力学を学びたい読者ならだれでも理解できる形を目指しております。本稿はこちらの子記事となっておりますので、よければ合わせてご覧ください。
では、親記事で説明したラグランジュ力学の復習からしていきます。
1 最小作用の原理
まずは、次のように物理量として「ラグランジアン\(L\)」と「作用\(S\)」を定義して、次のようなことを原理としましょう。
定義1.1(ラグランジアンの定義)
系の一般化座標を\(\boldsymbol{q}=(q_1,q_2,…,q_N)\)としたとき、各時刻における系の状態\((\boldsymbol{q}, \dot{\boldsymbol{q}}, t)\)を引数として、その系の自由度や対称性、相互作用を正しく表すように選ばれる物理量\(L(\boldsymbol{q}, \dot{\boldsymbol{q}}, t)\)を、ラグランジアンと定義する。
定義1.2(作用の定義)
系の一般化座標を\(\boldsymbol{q}=(q_1,q_2,…,q_N)\)、系に対応するラグランジアンを\(L(\boldsymbol{q}, \dot{\boldsymbol{q}}, t)\)としたとき、
$$
S=\int_{t_i}^{t_f} L(\boldsymbol{q}, \dot{\boldsymbol{q}}, t) dt
$$なる物理量\(S\)を、時刻\(t_i\)から時刻\(t_f\)までの作用と定義する。
原理1.3(最小作用の原理)
現実に物体が通過する経路\(\boldsymbol{q}(t)\)においては作用\(S\)が停留する、すなわち、始点\(\boldsymbol{q}(t_i)\)と終点 \(\boldsymbol{q}(t_f)\)を固定したときの作用の変分\(\delta S\)がゼロになる。
かなり突飛にいろんなことが出てきて、何を言っているのかわからないように思えますが、くわしくはこちらで説明しています。
大事なことだけ書いておくと、原理というのは、証明はできないけれどたしかに現実の運動がそれにしたがっている、という主張のことで、古典力学におけるニュートンの運動方程式、古典電磁気学におけるマクスウェル方程式みたいなものですね。
今回、解析力学をはじめるにあたって、原理1.3、すなわち「現実に物体が通過する経路\(\boldsymbol{q}(t)\)においては、作用\(S=\int_{t_i}^{t_f} L(\boldsymbol{q}, \dot{\boldsymbol{q}}, t) dt\)が停留する」という主張が、現実に起こる運動と非常によく整合しているから、証明なしに受け入れてしまおうということです。
2 オイラー=ラグランジュ方程式
先ほどの原理1.3から、以下の定理1.4が導かれます。
定理1.4(オイラー=ラグランジュ方程式)
系の一般化座標を\(\boldsymbol{q}=(q_1,q_2,…,q_N)\)、系に対応するラグランジアンを\(L(\boldsymbol{q}, \dot{\boldsymbol{q}}, t)\)としたとき、\(i=1,2,…,N\)に対して以下が成り立つ。
$$
\frac{\partial L}{\partial q_i} – \frac{d}{dt} \frac{\partial L}{\partial \dot{q}_i} =0
$$
さて、ここまでがラグランジュ力学の出発点です。本稿で行う議論は、ここまでわかっていれば大丈夫です。
第2章―ラグランジュ力学を連続体に拡張する
さて、ここまでのラグランジュ力学についての議論は、系の自由度が有限であって、かつ離散的であることを前提としているため、空間に連続的に分布する弦の振動や流体の運動、あるいは空間そのものの性質である電磁場のように、無限個の自由度を持つ連続的な物理系を扱うときには、この枠組みをそのまま適用することができないという限界があります。
そこで本稿では、このような空間的に連続な系もラグランジュ形式で扱えるように、理論を拡張することを目指していきます。
1 連続体への翻訳
有限の自由度を持つ離散的な系から、無限の自由度を持つ連続的な系へ、理論を拡張するために、次のように離散的な量の連続的な極限をとることを考えます。
これまで、離散的な物理系における\(N\)個の質点は、離散的な添字\(i=1, 2, …, N\)によって区別され、それぞれの質点の状態は、時間\(t\)のみに依存する一般化座標\(q_i(t)\)によって記述されていました。これに対して、連続的な物理系においては、対象となる物理量が空間全体に隙間なく連続して分布しているとみなすので、個々の要素を区別するために用いられていた離散的な添字\(i\)は、連続的に変化する空間座標\(\boldsymbol{x}\)へと置き換えられます。
この添字の連続化に伴って、時間\(t\)のみに依存していた一般化座標\(q_i(t)\)を、空間座標 \(\boldsymbol{x}\)と時間\(t\)の独立変数双方に依存する関数\(\phi(\boldsymbol{x}, t)\)へと拡張することにします(この新しい変数は一般に「場」と呼ばれます)。
以上をまとめると、有限自由度の系から無限自由度の系への移行にあたって、\(\{q_i(t)\}_{i=1}^N \rightarrow \phi(\boldsymbol{x}, t)\)のように、状態変数を変化させるということです。
また、離散的な物理系における系全体のラグランジアン\(L\)は、各粒子\(i\)のラグランジアン\(L_i\)の和として表現されるため、連続的な極限においては、ラグランジアンを空間の各点に「ラグランジアン密度\(\mathcal{L}\)」として分配して、\(L = \int \mathcal{L} d\boldsymbol{x}\)とかけます。このラグランジアン密度\(\mathcal{L}\)がどのような変数に依存するかを考えましょう。
離散的な物理系のラグランジアン\(L(q_i, \dot{q}_i, t)\)から類推すると、\(\mathcal{L}\)は場の値 \(\phi(\boldsymbol{x}, t)\)およびその時間微分である場の変化率\(\dot{\phi}(\boldsymbol{x}, t)\)に依存するでしょう。
さらに連続体特有の重要な性質として、\(\mathcal{L}\)は場の空間微分\(\nabla\phi\)にも依存します。というのは、離散的な物理系において、その極限が連続体であるような系(バネで繋がれた隣り合う無数の質点など)を考えると、そのラグランジアンは隣接する座標の差\(q_{i+1} – q_i\)の項を含み、この微小な座標の差は極限で場の空間微分\(\nabla\phi\)になるためです。
したがって、連続的な物理系を記述するラグランジアン密度は、一般に\(\mathcal{L}(\phi, \nabla\phi, \dot{\phi}, t)\)という関数形で表現されるのが自然です。
2 場におけるオイラー=ラグランジュ方程式
以上より、系全体の作用\(S\)は、時間\(t\)および空間体積\(V\)に関する積分として、
$$
S = \int_{t_1}^{t_2} L dt = \int_{t_1}^{t_2} \int_V \mathcal{L}\left(\phi, \nabla\phi, \dot{\phi}, t\right) d\boldsymbol{x} dt
$$のように定義されます。これに対してラグランジュ力学は、現実にあらわれる運動が作用\(S\)を停留させるような場の変動を選択するという”最小作用の原理”を、基本的な枠組みとして要請するので、定理として以下が成り立ちます。
定理2.1(\(N\)次元の場におけるオイラー=ラグランジュ方程式)
\(N\)次元の場\(\phi(\boldsymbol{x},t)\)において、ラグランジアン密度を\(\mathcal{L}(\phi, \nabla\phi, \dot{\phi},t)\)とすると、以下が成り立つ。
$$
\frac{\partial}{\partial t}\left( \frac{\partial \mathcal{L}}{\partial \dot{\phi}} \right) + \sum_{j=1}^N \frac{\partial}{\partial x_j}\left( \frac{\partial \mathcal{L}}{\partial (\partial_{x_j}\phi)} \right) – \frac{\partial \mathcal{L}}{\partial \phi} = 0
$$
定理2.1の証明
先に述べたように、ラグランジアン密度\(\mathcal{L}\)を用いた作用\(S\)は、時間\(t\)および空間体積\(V\)に関する二重積分として、
$$
S = \int_{t_1}^{t_2} L dt = \int_{t_1}^{t_2} \int_V \mathcal{L}\left(\phi, \nabla\phi, \dot{\phi}, t\right) d\boldsymbol{x} dt
$$と定義される。ここで、現実の場の変動\(\phi(\boldsymbol{x},t)\)に対して、時間と空間の境界を固定した仮想的な微小変化\(\delta \phi(\boldsymbol{x},t)\)をさせた経路\(\phi+\delta \phi\)を考えると、対応する作用の変分\(\delta S\)がゼロになる。
一方で、この変分\(\delta S\)を二重積分の表式から計算すると、ラグランジアン密度\(\mathcal{L}\)の、変数\(\phi, \nabla\phi, \dot{\phi}=\frac{\partial \phi}{\partial t}\)についての多変数関数のテイラー展開を用いて、
\begin{align}
\delta S &= \int_{t_1}^{t_2} \int_V \left\{ \frac{\partial \mathcal{L}}{\partial \phi}\delta\phi + \sum_{j=1}^N \frac{\partial \mathcal{L}}{\partial \left( \frac{\partial \phi}{\partial x_j} \right)}\delta \left( \frac{\partial \phi}{\partial x_j} \right) + \frac{\partial \mathcal{L}}{\left( \frac{\partial \phi}{\partial t} \right)} \delta\left( \frac{\partial \phi}{\partial t} \right) \right\} \, d\boldsymbol{x} \, dt \\
&= \int_{t_1}^{t_2} \int_V \left\{ \frac{\partial \mathcal{L}}{\partial \Phi} \delta \phi + \sum_{j=1}^N \frac{\partial \mathcal{L}}{\partial \left( \frac{\partial \Phi}{\partial x_j} \right)} \frac{\partial}{\partial x_j} (\delta \phi) + \frac{\partial \mathcal{L}}{\partial \left( \frac{\partial \phi}{\partial t} \right)} \frac{\partial}{\partial t} (\delta \phi) \right\} \, d\boldsymbol{x} \, dt \\
&= \int_{t_1}^{t_2} \int_V \left\{ \frac{\partial \mathcal{L}}{\partial \phi} – \sum_{j=1}^N \frac{\partial}{\partial x_j} \left( \frac{\partial \mathcal{L}}{\partial \left( \frac{\partial \phi}{\partial x_j} \right)} \right) – \frac{\partial}{\partial t} \left( \frac{\partial \mathcal{L}}{\partial \left( \frac{\partial \phi}{\partial t} \right)} \right) \right\} \delta \phi \, d\boldsymbol{x} \, dt \\
\end{align}とかける。ただし、2つ目の等号は、極限の定義を考えると、
$$
\delta \left( \frac{\partial}{\partial a} \Phi(a) \right) = \lim_{a_2 \to a_1} \left( \frac{\left( \Phi(a_2)+\delta \Phi(a_2) \right) – \left( \Phi(a_1)+\delta \Phi(a_1) \right)}{a_2-a_1} – \frac{\Phi(a_2) – \Phi(a_1)}{a_2-a_1} \right) = \frac{\partial}{\partial a} \delta \Phi(a)
$$より、変分の微分と微分の変分が等しいこと、つまり変分と偏微分の順序を入れ替えることができること、3つ目の等号は、第2項の各\(j\)については空間について部分積分を行うことで、
\begin{align*}
\int_{x_{j,1}}^{x_{j,2}} dx_j \left[ \frac{\partial \mathcal{L}}{\partial \left( \frac{\partial \phi}{\partial x_j} \right)} \frac{\partial}{\partial x_j} (\delta \phi) \right] &= \left[ \frac{\partial \mathcal{L}}{\partial \left( \frac{\partial \phi}{\partial x_j} \right)} \delta \phi \right]_{x_{j,1}}^{x_{j,2}} – \int_{x_{j,1}}^{x_{j,2}} dx_j \left[ \frac{\partial}{\partial x_j} \left( \frac{\partial \mathcal{L}}{\partial \left( \frac{\partial \Phi}{\partial x_j} \right)} \right) \delta \phi \right] \\
&= – \int_{x_{j,1}}^{x_{j,2}} dx_j \left[ \frac{\partial}{\partial x_j} \left( \frac{\partial \mathcal{L}}{\partial \left( \frac{\partial \Phi}{\partial x_j} \right)} \right) \delta \phi \right]
\end{align*}となること、第3項については時間について部分積分を行うことで、
\begin{align*}
\int_{t_1}^{t_2} dt \left[ \frac{\partial \mathcal{L}}{\partial \left( \frac{\partial \Phi}{\partial t} \right)} \frac{\partial}{\partial t} (\delta \phi) \right] &= \left[ \frac{\partial \mathcal{L}}{\partial \left( \frac{\partial \Phi}{\partial t} \right)} \delta \Phi \right]_{t_1}^{t_2} – \int_{t_1}^{t_2} dt \left[ \frac{\partial}{\partial t} \left( \frac{\partial \mathcal{L}}{\partial \left( \frac{\partial \Phi}{\partial t} \right)} \right) \delta \Phi \right] \\
&= – \int_{t_1}^{t_2} dt \left[ \frac{\partial}{\partial t} \left( \frac{\partial \mathcal{L}}{\partial \left( \frac{\partial \Phi}{\partial t} \right)} \right) \delta \Phi \right]
\end{align*}となることを用いた。さて、任意の変分\(\delta \phi\)に対して、\(\delta S = 0\)が成り立つことから、変分法の基本補題より、中括弧\(\{\quad\}\)の中身自体がすべての\(\boldsymbol{x}, t\)において恒等的にゼロでなければならないので、以下が成り立つことが示される。
$$
\frac{\partial}{\partial t}\left( \frac{\partial \mathcal{L}}{\partial \dot{\phi}} \right) + \sum_{j=1}^N \frac{\partial}{\partial x_j}\left( \frac{\partial \mathcal{L}}{\partial (\partial_{x_j}\phi)} \right) – \frac{\partial \mathcal{L}}{\partial \phi} = 0
$$□
さて、この方程式によって、場の力学も解析力学の枠組みに落とし込んで考えられるようになりました。弦や電磁場など具体的な例への適用を見たい方は、こちらの第2章の7節をご覧ください。
以上、ラグランジュ力学の連続体への拡張を説明してみました。お疲れさまでした

コメント