ラグランジュ力学

**ラグランジュ力学**（ラグランジュりきがく、英語：）は、一般化座標とその微分を基本変数として記述された古典力学である。フランスの物理学者ジョゼフ＝ルイ・ラグランジュが創始した。後のハミルトン力学と同様にニュートン力学を再定式化した解析力学の一形式である。

## 概要
ラグランジュ形式の解析力学は**[[最小作用の原理]]**によって構成される。
元々はニュートン的な力学の分野において成立したが、電磁気学や相対性理論でも応用することが出来て、これらの分野における基礎方程式（マクスウェル方程式、アインシュタイン方程式）を導き出すことが出来る。
また、量子力学においても、経路積分の方法は最小作用の原理に関連して考え出された方法である。

ラグランジュ形式では一般化座標によって記述されており、変数の取り方が任意である。
ニュートンの運動方程式はベクトルの方程式であり、デカルト座標以外では煩雑な座標変換が必要となるが、ラグランジュ形式においてはラグランジアンはスカラーであり座標変換が簡単である。

実際の計算上でも、例えば長さが一定の振り子などで円周上を運動する場合には、平面内の運動なのでニュートンの運動方程式では2つの方向の2変数が必要となるが、ラグランジュ形式では一般化座標として角度を選ぶことにより1変数の方程式が得られる。
もちろんニュートンの運動方程式はラグランジュ形式と等価なので適当な変換により同じ式が得られるが、ラグランジュ形式では直接得られる点で便利である。

## 定式化
ラグランジュ形式において、力学系の運動状態を指定する力学変数は**一般化座標** q(t)=(q_1(t),\ldots) である。力学系の性質は一般化座標とその微分（一般化速度）、および時間を変数とする関数 L(q(t), \dot{q}(t), t) によって記述される。この力学系の性質を記述する関数 L は**ラグランジュ関数**（**ラグランジアン**）と呼ばれる。

ラグランジュ形式において、作用汎関数はラグランジュ関数の時間積分
^{t_\text{F}} L(q(t), \dot{q}(t), t)\, dt
}}
として与えられる。
一般化座標は実際には起こらない運動の値も取りうるが、そこから実際の運動を導く方法が**[[最小作用の原理]]**である。すなわち、作用汎関数が最小となる運動が実際に起こる運動である。

作用の停留条件から、**ラグランジュの運動方程式**（オイラー＝ラグランジュ方程式）

が得られる。
これはニュートンの運動方程式と同等である。

### 運動量
一般化座標に共役な一般化運動量は、ラグランジアンの一般化速度による偏微分

によって定義される。
これは並進対称性から導かれる保存量である。

一般化運動量を用いると、ラグランジュの運動方程式は

となる。ニュートンの運動方程式との比較から、右辺は一般化された力と見ることも出来る。

ハミルトン形式では一般化座標と一般化運動量によって記述されている。
一般化運動量は正準共役量であり、共役運動量や正準運動量と呼ばれることもある。

## ラグランジュ関数
**ラグランジュ関数**（**ラグランジアン**、）は、物理的な力学系の動力学を記述するために用いられる関数である。
ラグランジアン L(q,\dot{q},t) は一般に運動エネルギー  とポテンシャル  の差

の形で書かれる。

ラグランジアンはエネルギーの次元を持つスカラーであるが、観測可能な物理量ではなく、その値自体に物理的な意味があるわけではない。特に、座標と時間の任意関数 f(q,t) の時間による全微分を加える変換

を行っても全く同じ力学系を表す。
この全微分は連鎖律により

となるので、この変換に対して、共役運動量は

と変換される。したがって、新たな共役運動量の時間微分は

となる。一方、一般化された力は

と変換される。任意関数  に作用する全微分  と座標の偏微分  が交換可能なので、この変換に対して運動方程式が保たれる。

## 座標変換
座標変換 q \mapsto Q が

で表されるとき、新たな座標の下でのラグランジアンは

で与えられ、新たなラグランジアンから導かれる運動方程式は
{\partial Q_I} -\frac{d}{dt} \frac{\partial\tilde{L}}{\partial\dot{Q}_I} =0
}}
である。このように写像の合成で座標変換を容易に行えることが一般化座標で表されるラグランジュ形式の利点の一つである。

座標変換の時間微分は連鎖律により

であるため、新たな座標に共役な運動量は
{\partial\dot{Q}_I} =\frac{\partial L}{\partial\dot{q}_i} \frac{\partial g_i}{\partial Q_I}
 =p_i \cdot \frac{\partial g_i}{\partial Q_I}
}}
となる。

### 母関数
座標変換は

で定義される母関数により生成される。
座標変換は

で与えられ、新たな運動量は

で与えられる。

先の任意関数によるラグランジュ関数の変換を伴う場合の母関数は

で与えられる。

## 拘束系
拘束条件が課された系にラグランジュ形式を用いる際に、一般座標を適当に選ぶことによって、拘束条件が常に満たされるようにすることができる。上で挙げた振り子の例であれば、座標変数に角度を選ぶことによって長さが一定という拘束条件が常に満たされるようにしている。
これの手法とは別に、ラグランジュの未定乗数法を用いて作用汎関数（ラグランジュ関数）に拘束条件を取り入れる方法がある。

一般化座標  に対して、拘束条件

が課されている場合を考える。
このとき、作用は
^{t_\text{F}} \beta(t)\, \varPhi(q,t)\, dt
}}
によって拘束条件が取り入れられる。ここで導入された  がラグランジュの未定乗数である。
拘束条件は全ての時間で成り立つので、未定乗数も各々の時間に対して導入される時間の関数である。

拘束条件が取り入れられた作用に対して最小作用の原理を適用して

が得られる。力学変数  に対応する運動方程式には「拘束力」 が加えられ、未定乗数に対応する運動方程式として拘束条件が導かれる。

## ハミルトン形式との関係
ハミルトン形式とラグランジュ形式はルジャンドル変換を通して等価である。ただし、ラグランジアンが退化している場合は、ルジャンドル変換が微分同相写像ではなくなり、ラグランジュ系からハミルトン系へ移行することができなくなる。この退化している場合の処方としてポール・ディラックの拘束理論が知られている。

## ラグランジュ形式による場の理論

特に相対論的な場の理論の場合では、ラグランジュ形式から出発するのが一般的である。その方が相対論的不変性などの対称性が見やすいからである。

力学変数としては場 \phi(x)を考える。作用積分は**ラグランジアン密度** \mathcal{L}(\phi, \partial\phi, x) により

で書かれる。その変分は
{\partial\phi} \delta\phi
 +\frac{\partial\mathcal{L}}{\partial(\partial_\mu\phi)}\partial_\mu\delta\phi \right) \sqrt{-g}\, d^dx \\
 &=\frac{1}{c} \int \left[ \frac{\partial\mathcal{L}}{\partial\phi}
 -\frac{1}{\sqrt{-g}} \partial_\mu \left( \frac{\partial\mathcal{L}}{\partial(\partial_\mu\phi)}
 \sqrt{-g} \right) \right] \delta\phi \sqrt{-g}\, d^dx
 +\frac{1}{c} \oint \frac{\partial\mathcal{L}}{\partial(\partial_\mu\phi)} \delta\phi \sqrt{-g}\, d\varSigma_\mu \\
\end{aligned}
}}
となり、ラグランジュの運動方程式として
 \frac{\delta S[\phi]}{\delta\phi(x)} =\frac{\partial\mathcal{L}}{\partial\phi}
 -\frac{1}{\sqrt{-g}} \partial_\mu \left( \frac{\partial\mathcal{L}}{\partial(\partial_\mu\phi)} \sqrt{-g} \right) =0
}}
が得られる。

## ラグランジュ関数の存在条件
座標の2階微分  について高々1次である次の運動方程式
:A_{ij}(q,\dot{q},t)\ddot{q}^j + B_i(q,\dot{q},t)=0, \quad (i,j=1,\dots,N)
（ただし ）を導くラグランジュ関数が局所的に存在する必要十分条件は以下であることがヘルムホルツにより調べられている：
:\begin{align}
A_{ij} &= A_{ji},\\
\frac{\partial A_{jk}}{\partial \dot{q}^i} &= \frac{\partial A_{ki}}{\partial \dot{q}^j},\\
\frac{\partial B_i}{\partial \dot{q}^j}+\frac{\partial B_j}{\partial \dot{q}^i} &= 2\left(\dot{q}^k\frac{\partial }{\partial q^k}+\frac{\partial }{\partial t}\right)A_{ij},\\
2\left(\frac{\partial B_i}{\partial q^j}-\frac{\partial B_j}{\partial q^i}\right) &= \left(\dot{q}^k\frac{\partial }{\partial q^k}+\frac{\partial }{\partial t}\right) \left(\frac{\partial B_i}{\partial \dot{q}^j}-\frac{\partial B_j}{\partial \dot{q}^i}\right).
\end{align}
このとき、ラグランジュ関数は以下で与えられる：
:\begin{align}
L(q,\dot{q},t) &= K(q,\dot{q},t)+D_i(q,t)\dot{q}^i+C(q,t),\\
K(q,\dot{q},t) &= \int_0^1 \dot{q}^iH_i(q,\dot{q}\tau,t)\mathrm{d}\tau,\\
D_i(q,t) &= q^j\int_0^1 \mathrm{d}\tau\ \tau Z_{[i,j]}(\tau q,t)+\frac{\partial G}{\partial q^i},\\
C(q,t) &= q^i\int_0^1 \mathrm{d}\tau\ Y_i(q\tau,t),\\
H_i(q,\dot{q},t) &:= \int_0^1 \mathrm{d}\tau\ A_{ij}(q,\tau\dot{q},t)\dot{q}^j,\\
Z_{[i,j]}(q,t) &\equiv -Z_{[j,i]}(q,t) := \frac{1}{2}\left(\frac{\partial B_i}{\partial \dot{q}^j}-\frac{\partial B_j}{\partial \dot{q}^i}\right)+\frac{\partial^2 K}{\partial q^i \partial \dot{q}^j}-\frac{\partial^2 K}{\partial \dot{q}^i \partial q^j},\\
Y_i(q,t) &:= \left\{\frac{\partial^2 K}{\partial q^i \partial \dot{q}^j}+\frac{1}{2}\left(\frac{\partial B_i}{\partial \dot{q}^j}-\frac{\partial B_j}{\partial \dot{q}^i}\right)\right\}\dot{q}^j - \frac{\partial K}{\partial q^i} - B_i + \frac{\partial^2 K}{\partial \dot{q}^i \partial t} + \frac{\partial D_i}{\partial t}.
\end{align}
ここで  は  の任意関数である。

## 具体例
### 相対論的な粒子系
相対論的な系では、時間は位置と共に4元ベクトルとなるので、時間は力学変数となり、運動のパラメータではなくなる。パラメータを  として、力学変数を

とする。ここで  は時空の添え字で、 は粒子を区別する添え字である。
自由粒子系を考えると、作用積分は

である。ここで  は平坦な時空の計量で \eta = \mathrm{diag}(-1, 1, \ldots, 1) である。
平方根の中が正である為に、作用積分の段階で運動は時間的なものに限定されている。

ラグランジュの運動方程式は

となる。
ここで、一般化運動量は

}}

}}
である。
固有時間 c^2d\tau_i^2 = \eta_{\rho\sigma}dX_i^\rho dX_i^\sigma を使うと

となる。

#### 補助変数の導入
この作用は平方根の中に微分を含む形のため扱いが困難である。
補助変数  を導入して別の形に書くことが出来る。

としてマクスウェルの方程式が導かれる。

#### 電磁場中の粒子系
物質場として相対論的な粒子系を考え、相互作用項として

を考える。

このとき、物質  に関する運動方程式は

となり、ローレンツ力を再現する。

また、4元電流密度は
 \frac{\delta S_\text{int}[X,A]}{\delta A_\mu(x)}
 =\sum_i \frac{q_ic}{\sqrt{-g}} \int \dot{X}_i^\mu(\lambda)\, \delta^4(X_i(\lambda)-x)\, d\lambda
}}
となる。

### 一般相対性理論
一般相対性理論においては、平坦な時空の計量は曲がった時空の計量  に置き換えられ、これが力学変数となる。
作用積分は

と書かれる。
重力場の項は

である。
ここで  はスカラー曲率である。
アインシュタイン方程式は時空の計量  の運動方程式として導かれる。

## 脚注

### 注釈

### 出典


## 参考文献
* 
* 
* 
* 

## 関連項目
* オイラー＝ラグランジュ方程式
* 最小作用の原理
* 解析力学 - ハミルトン力学
* ネーターの定理

## 外部リンク
*  - スカラーペディア百科事典「ラグランジュ力学」の項目。
* 大貫義郎:「拘束系の経路積分」数理解析研究所講究録、1260巻、2002年、pp.62-76.
* 須藤大樹：拘束系の正準量子化と経路積分量子化」(平成18年)