⚙️

解析力学入門:ニュートン力学を超えて物理学の本質に迫る

に公開

はじめに

物理学を学ぶ理系学生にとって、解析力学は避けて通れない重要な分野です。高校や大学初年度で学ぶニュートン力学から一歩進んで、より洗練された数学的枠組みで力学現象を記述する——それが解析力学です。

「なぜわざわざニュートン力学を言い換える必要があるのか?」

この疑問を持つ学生は多いでしょう。しかし、解析力学を学ぶことで、物理学の見方が根本から変わります。量子力学、統計力学、場の理論、さらには一般相対性理論まで、現代物理学のほぼすべての分野が解析力学の言葉で書かれています。

本記事では、解析力学とは何か、なぜそれが必要なのか、そしてどのような概念から構成されているのかを、できる限り詳しく、かつ分かりやすく解説します。


第1章:ニュートン力学の限界と解析力学の必要性

1.1 ニュートン力学の偉大さ

まず、ニュートン力学の素晴らしさを確認しましょう。ニュートンの運動方程式

\mathbf{F} = m\mathbf{a}

は、驚くほどシンプルでありながら、惑星の運動から日常的な物体の運動まで、広範な現象を説明できます。17世紀にニュートンがこの法則を定式化して以来、物理学は飛躍的に発展しました。

1.2 ニュートン力学の困難

しかし、ニュートン力学にはいくつかの実用上の困難があります。

座標系への依存

ニュートンの運動方程式は、直交座標系(デカルト座標)で書くと最もシンプルです。

m\ddot{x} = F_x, \quad m\ddot{y} = F_y, \quad m\ddot{z} = F_z

しかし、問題によっては極座標や球座標を使った方が自然な場合があります。例えば、中心力場での運動を考えるとき、直交座標では方程式が複雑になります。

極座標 (r, \theta) での運動方程式は

m(\ddot{r} - r\dot{\theta}^2) = F_r
m(r\ddot{\theta} + 2\dot{r}\dot{\theta}) = F_\theta

となり、r\dot{\theta}^2 や 2\dot{r}\dot{\theta} といった「見かけの力」に対応する項が現れます。座標系を変えるたびに、このような補正項を手作業で導出しなければなりません。

拘束条件の扱い

多くの物理系には拘束条件があります。例えば:

  • 斜面上を滑る物体(斜面から離れない)
  • 振り子(紐の長さが一定)
  • 二重振り子(2つの紐の長さがそれぞれ一定)
  • 剛体(すべての点間の距離が一定)

ニュートン力学では、これらの拘束条件を抗力(垂直抗力、張力など)として扱います。しかし、抗力は運動の結果として決まるものであり、あらかじめ分からないことが多いのです。

例えば、単振り子の問題では、紐の張力 T を未知数として運動方程式を立て、拘束条件(紐の長さ一定)から T を消去する必要があります。問題が複雑になると、この作業は非常に煩雑になります。

多体系の複雑さ

N 個の質点からなる系を考えると、ニュートン力学では 3N 個の2階微分方程式を連立して解く必要があります。さらに、質点間に相互作用がある場合、各質点への力を個別に計算しなければなりません。

1.3 解析力学が解決すること

解析力学は、これらの困難を以下のように解決します:

  1. 座標系に依存しない定式化:どんな座標系でも同じ形の方程式が成り立つ
  2. 拘束条件の自然な取り込み:拘束力を陽に考えなくてよい
  3. 系全体を1つの量で特徴づける:ラグランジアンやハミルトニアンという1つの関数から、すべての運動方程式が導かれる
  4. 保存則の系統的な導出:対称性と保存則の関係が明確になる
  5. 他の分野への拡張:量子力学や場の理論への自然な橋渡し

第2章:変分原理——物理学の新しい視点

2.1 最小作用の原理

解析力学の根底には、変分原理という考え方があります。その中心となるのが最小作用の原理(より正確には停留作用の原理)です。

ニュートン力学では「各瞬間に力がどう働くか」を考えます。これに対し、変分原理では「運動全体を見て、ある量が極値をとるような経路が実現する」と考えます。

作用の定義

時刻 t_1 から t_2 まで運動する系を考えます。系の状態は位置 q と速度 \dot{q} で記述されるとします。このとき、作用(action)S を次のように定義します:

S[q] = \int_{t_1}^{t_2} L(q, \dot{q}, t) \, dt

ここで L はラグランジアンと呼ばれる関数です。作用 S は経路 q(t) の関数ではなく、経路全体を引数とする汎関数であることに注意してください。S[q] という記法でこれを表しています。

停留作用の原理

実際に実現する運動は、作用 S を停留させる(極値にする)経路である。

より数学的に言えば、始点 q(t_1) = q_1 と終点 q(t_2) = q_2 を固定したとき、経路を微小変化させても作用の1次の変化がゼロになる経路が物理的に実現します:

\delta S = 0

2.2 変分原理の歴史と哲学

変分原理のアイデアは古くからありました。

フェルマーの原理(1662年)

光学において、フェルマーは「光は2点間を最短時間で結ぶ経路を通る」という原理を提唱しました。これにより、反射の法則や屈折の法則(スネルの法則)が導かれます。

モーペルテュイの原理(1744年)

モーペルテュイは「自然は常に最小の作用で動く」という原理を提唱しました。彼の定式化は曖昧でしたが、この考え方が後の解析力学につながりました。

オイラーとラグランジュ(18世紀)

オイラーとラグランジュが変分法を数学的に厳密化し、解析力学の基礎を築きました。特にラグランジュの著作『解析力学』(Mécanique analytique, 1788年)は、力学を純粋に解析的な方法で記述した画期的な著作でした。

哲学的な意味

変分原理は、しばしば目的論的(teleological)な解釈を誘います。「自然は最も効率的な道を選ぶ」という見方です。しかし、現代物理学では、変分原理は因果的な法則(各瞬間での法則)と数学的に等価であり、どちらが「本質的」かという問いには明確な答えはないと考えられています。

むしろ、変分原理の強みは、対称性と保存則の関係を明らかにすることや、量子力学への拡張(ファインマンの経路積分)に自然につながることにあります。

2.3 変分法の数学

変分原理を理解するために、変分法の基本を学びましょう。

汎関数とその変分

通常の関数は数を数に対応させます:f: \mathbb{R} \to \mathbb{R}

汎関数は関数を数に対応させます:S: \{q(t)\} \to \mathbb{R}

汎関数の例:

  • 曲線の長さ:L[y] = \int_a^b \sqrt{1 + (y')^2} \, dx
  • 作用:S[q] = \int_{t_1}^{t_2} L(q, \dot{q}, t) \, dt

汎関数 S[q] の変分 \delta S は、経路を微小変化させたときの S の変化です。経路を q(t) \to q(t) + \epsilon \eta(t) と変化させると(\eta(t_1) = \eta(t_2) = 0 として端点を固定)、

\delta S = \lim_{\epsilon \to 0} \frac{S[q + \epsilon\eta] - S[q]}{\epsilon}

オイラー=ラグランジュ方程式の導出

作用の変分を計算しましょう:

S[q + \epsilon\eta] = \int_{t_1}^{t_2} L(q + \epsilon\eta, \dot{q} + \epsilon\dot{\eta}, t) \, dt

\epsilon について1次まで展開すると:

S[q + \epsilon\eta] \approx S[q] + \epsilon \int_{t_1}^{t_2} \left( \frac{\partial L}{\partial q}\eta + \frac{\partial L}{\partial \dot{q}}\dot{\eta} \right) dt

第2項を部分積分します。\dot{\eta} を含む項について:

\int_{t_1}^{t_2} \frac{\partial L}{\partial \dot{q}}\dot{\eta} \, dt = \left[ \frac{\partial L}{\partial \dot{q}}\eta \right]_{t_1}^{t_2} - \int_{t_1}^{t_2} \frac{d}{dt}\left(\frac{\partial L}{\partial \dot{q}}\right)\eta \, dt

\eta(t_1) = \eta(t_2) = 0 より、境界項は消えます。したがって:

\delta S = \int_{t_1}^{t_2} \left( \frac{\partial L}{\partial q} - \frac{d}{dt}\frac{\partial L}{\partial \dot{q}} \right)\eta \, dt

\delta S = 0 がすべての \eta(t) に対して成り立つためには:

\boxed{\frac{d}{dt}\frac{\partial L}{\partial \dot{q}} - \frac{\partial L}{\partial q} = 0}

これがオイラー=ラグランジュ方程式(Euler-Lagrange equation)です。


第3章:ラグランジュ形式の力学

3.1 一般化座標

解析力学の強力さの源は、一般化座標(generalized coordinates)という概念にあります。

自由度

系の自由度とは、系の状態を完全に指定するのに必要な独立なパラメータの数です。

  • 空間中の1質点:自由度3(x, y, z)
  • 平面上の1質点:自由度2(x, y または r, \theta)
  • 単振り子:自由度1(角度 \theta)
  • 二重振り子:自由度2(2つの角度 \theta_1, \theta_2)
  • 剛体:自由度6(重心位置3 + 回転3)

N 個の質点からなる系で、k 個の拘束条件がある場合、自由度は f = 3N - k です。

一般化座標の定義

系の自由度が f のとき、f 個の独立な変数 q_1, q_2, \ldots, q_f で系の配位(状態)を完全に記述できます。これらを一般化座標と呼びます。

一般化座標は必ずしも長さの次元を持つ必要はありません。角度、面積比、無次元量など、何でも構いません。重要なのは:

  1. 系の配位を完全に決定できること
  2. 拘束条件を自動的に満たすこと
  3. 互いに独立であること

例:二重振り子

二重振り子を考えましょう。2つの質点の位置は直交座標では4つの座標 (x_1, y_1, x_2, y_2) で記述されますが、2つの紐の長さ一定という拘束条件があるため、自由度は2です。

自然な一般化座標は、2つの振り子の角度 \theta_1, \theta_2 です。直交座標との関係は:

x_1 = l_1 \sin\theta_1, \quad y_1 = -l_1 \cos\theta_1
x_2 = l_1 \sin\theta_1 + l_2 \sin\theta_2, \quad y_2 = -l_1 \cos\theta_1 - l_2 \cos\theta_2

この変換により、拘束条件(紐の長さ一定)は自動的に満たされます。

3.2 ラグランジアン

ラグランジアンの定義

ラグランジアン L は次のように定義されます:

\boxed{L = T - V}

ここで T は運動エネルギー、V はポテンシャルエネルギーです。

この定義は一見すると天下り的ですが、以下で見るように、この定義によってオイラー=ラグランジュ方程式がニュートンの運動方程式と等価になります。

1次元の場合

1次元での自由粒子を考えます。

T = \frac{1}{2}m\dot{x}^2, \quad V = V(x)
L = \frac{1}{2}m\dot{x}^2 - V(x)

オイラー=ラグランジュ方程式を適用します:

\frac{\partial L}{\partial x} = -\frac{dV}{dx}, \quad \frac{\partial L}{\partial \dot{x}} = m\dot{x}
\frac{d}{dt}(m\dot{x}) - \left(-\frac{dV}{dx}\right) = 0
m\ddot{x} = -\frac{dV}{dx} = F

これはまさにニュートンの運動方程式です!

多自由度系の場合

f 個の一般化座標 q_1, \ldots, q_f を持つ系では、f 個のオイラー=ラグランジュ方程式が得られます:

\frac{d}{dt}\frac{\partial L}{\partial \dot{q}_i} - \frac{\partial L}{\partial q_i} = 0 \quad (i = 1, 2, \ldots, f)

3.3 具体例:様々な系のラグランジアン

例1:単振り子

長さ l、質量 m の単振り子を考えます。一般化座標として角度 \theta を選びます。

位置と速度:

x = l\sin\theta, \quad y = -l\cos\theta
\dot{x} = l\dot{\theta}\cos\theta, \quad \dot{y} = l\dot{\theta}\sin\theta

運動エネルギー:

T = \frac{1}{2}m(\dot{x}^2 + \dot{y}^2) = \frac{1}{2}ml^2\dot{\theta}^2

ポテンシャルエネルギー(y = 0 を基準):

V = mgy = -mgl\cos\theta

ラグランジアン:

L = \frac{1}{2}ml^2\dot{\theta}^2 + mgl\cos\theta

オイラー=ラグランジュ方程式:

\frac{\partial L}{\partial \theta} = -mgl\sin\theta, \quad \frac{\partial L}{\partial \dot{\theta}} = ml^2\dot{\theta}
ml^2\ddot{\theta} + mgl\sin\theta = 0
\ddot{\theta} + \frac{g}{l}\sin\theta = 0

これは振り子の運動方程式そのものです。

例2:中心力場での運動

質量 m の粒子が中心力 F(r) = -dV/dr を受けて運動する場合を考えます。

極座標 (r, \theta) を使うと:

\dot{x} = \dot{r}\cos\theta - r\dot{\theta}\sin\theta
\dot{y} = \dot{r}\sin\theta + r\dot{\theta}\cos\theta

運動エネルギー:

T = \frac{1}{2}m(\dot{x}^2 + \dot{y}^2) = \frac{1}{2}m(\dot{r}^2 + r^2\dot{\theta}^2)

ラグランジアン:

L = \frac{1}{2}m(\dot{r}^2 + r^2\dot{\theta}^2) - V(r)

r についての方程式:

m\ddot{r} - mr\dot{\theta}^2 = -\frac{dV}{dr}

\theta についての方程式:

\frac{d}{dt}(mr^2\dot{\theta}) = 0

2番目の式は角運動量 L_z = mr^2\dot{\theta} の保存を表しています!ラグランジアンが \theta に依存しない(\theta について対称)ことから、この保存則が自動的に導かれます。

例3:二重振り子

二重振り子は、解析力学の威力を示す良い例です。

質量 m_1, m_2、長さ l_1, l_2 の二重振り子を考えます。

位置:

x_1 = l_1\sin\theta_1, \quad y_1 = -l_1\cos\theta_1
x_2 = l_1\sin\theta_1 + l_2\sin\theta_2, \quad y_2 = -l_1\cos\theta_1 - l_2\cos\theta_2

速度:

\dot{x}_1 = l_1\dot{\theta}_1\cos\theta_1, \quad \dot{y}_1 = l_1\dot{\theta}_1\sin\theta_1
\dot{x}_2 = l_1\dot{\theta}_1\cos\theta_1 + l_2\dot{\theta}_2\cos\theta_2
\dot{y}_2 = l_1\dot{\theta}_1\sin\theta_1 + l_2\dot{\theta}_2\sin\theta_2

運動エネルギー(計算を進めると):

T = \frac{1}{2}(m_1 + m_2)l_1^2\dot{\theta}_1^2 + \frac{1}{2}m_2l_2^2\dot{\theta}_2^2 + m_2l_1l_2\dot{\theta}_1\dot{\theta}_2\cos(\theta_1 - \theta_2)

ポテンシャルエネルギー:

V = -(m_1 + m_2)gl_1\cos\theta_1 - m_2gl_2\cos\theta_2

ラグランジアン L = T - V からオイラー=ラグランジュ方程式を導くと、連立した非線形微分方程式が得られます。この系はカオス的な振る舞いを示すことで有名です。

ニュートン力学で同じ問題を解こうとすると、各質点への力(重力と張力)を個別に書き下し、拘束条件から張力を消去する必要があり、非常に煩雑です。一方、ラグランジュ形式では、座標変換と微分だけで運動方程式が得られます。

3.4 一般化運動量と一般化力

一般化運動量

一般化座標 q_i に共役な運動量(一般化運動量)p_i を次のように定義します:

\boxed{p_i = \frac{\partial L}{\partial \dot{q}_i}}

直交座標では p_x = m\dot{x} となり、通常の運動量と一致します。

極座標では:

  • p_r = m\dot{r}(動径方向の運動量)
  • p_\theta = mr^2\dot{\theta}(角運動量)

一般化運動量は必ずしも「質量×速度」の次元を持ちません。角度が座標の場合、対応する運動量は角運動量の次元を持ちます。

一般化力

一般化座標 q_i に対応する一般化力 Q_i は次のように定義されます:

Q_i = \frac{\partial L}{\partial q_i}

ポテンシャル V(q) がある場合:

Q_i = -\frac{\partial V}{\partial q_i}

オイラー=ラグランジュ方程式の別の形

これらの定義を使うと、オイラー=ラグランジュ方程式は:

\dot{p}_i = \frac{\partial L}{\partial q_i}

または、ポテンシャル力の場合:

\dot{p}_i = -\frac{\partial V}{\partial q_i}

これは「運動量の時間変化 = 力」というニュートンの第2法則の一般化です。

3.5 循環座標と保存則

循環座標(可無視座標)

ラグランジアン L がある座標 q_j に陽に依存しないとき、その座標を循環座標(cyclic coordinate)または可無視座標(ignorable coordinate)と呼びます。

\frac{\partial L}{\partial q_j} = 0

このとき、オイラー=ラグランジュ方程式より:

\frac{d}{dt}\frac{\partial L}{\partial \dot{q}_j} = 0

すなわち、一般化運動量 p_j が保存します。

例:角運動量の保存

中心力場 V = V(r) では、ラグランジアンは \theta に依存しません:

L = \frac{1}{2}m(\dot{r}^2 + r^2\dot{\theta}^2) - V(r)

したがって \theta は循環座標であり、対応する運動量

p_\theta = mr^2\dot{\theta}

が保存します。これは角運動量の保存則です。

対称性と保存則

循環座標の存在は、系の対称性と深く関係しています。

  • \theta が循環座標 ⟺ 回転対称性 ⟺ 角運動量保存
  • 空間の一様性(並進対称性) ⟺ 運動量保存
  • 時間の一様性(L が t に陽に依存しない) ⟺ エネルギー保存

このように、「ラグランジアンがある座標 θ に依存しない」ということは、物理的には「系を回転させても物理法則が変わらない(回転対称性がある)」ことを意味します。この「対称性と保存則の対応」は、後に解説するネーターの定理として一般化されます。


第4章:ハミルトン形式の力学

4.1 ルジャンドル変換

ラグランジュ形式からハミルトン形式への移行には、ルジャンドル変換という数学的操作を用います。

熱力学での例

熱力学では、内部エネルギー U(S, V)(S: エントロピー、V: 体積)から、ルジャンドル変換によって自由エネルギーを定義します:

F(T, V) = U - TS

ここで T = \partial U/\partial S です。これにより、独立変数が (S, V) から (T, V) に変わります。

力学での適用

ラグランジアン L(q, \dot{q}, t) を、変数 \dot{q} について ルジャンドル変換します。

まず、一般化運動量を思い出しましょう:

p = \frac{\partial L}{\partial \dot{q}}

この関係から \dot{q} を p の関数として解きます:\dot{q} = \dot{q}(q, p, t)

ハミルトニアン H を次のように定義します:

\boxed{H(q, p, t) = p\dot{q} - L(q, \dot{q}, t)}

右辺の \dot{q} は \dot{q}(q, p, t) で置き換えます。

4.2 ハミルトニアンの物理的意味

自然な力学系の場合

「自然な」力学系(運動エネルギーが \dot{q} の2次同次式で、ポテンシャルが位置のみに依存)では、ハミルトニアンは全エネルギーに等しくなります。

運動エネルギーが

T = \frac{1}{2}\sum_{i,j} a_{ij}(q)\dot{q}_i\dot{q}_j

の形の場合、オイラーの同次関数定理より

\sum_i \dot{q}_i \frac{\partial T}{\partial \dot{q}_i} = 2T

したがって

H = \sum_i p_i\dot{q}_i - L = \sum_i \dot{q}_i\frac{\partial L}{\partial \dot{q}_i} - L = 2T - (T - V) = T + V

つまり H = T + V = E(全エネルギー)です。

一般の場合

しかし、時間に依存する座標変換を含む系や、磁場中の荷電粒子など、ハミルトニアンが全エネルギーと一致しない場合もあります。ハミルトニアンはあくまで「ルジャンドル変換で定義される量」であり、エネルギーとの同一視は常に成り立つわけではありません。

4.3 ハミルトンの正準方程式

ハミルトニアン H(q, p, t) の全微分を計算します:

dH = \frac{\partial H}{\partial q}dq + \frac{\partial H}{\partial p}dp + \frac{\partial H}{\partial t}dt

一方、H = p\dot{q} - L の定義から:

dH = \dot{q}dp + pd\dot{q} - \frac{\partial L}{\partial q}dq - \frac{\partial L}{\partial \dot{q}}d\dot{q} - \frac{\partial L}{\partial t}dt

p = \partial L/\partial \dot{q} より、第2項と第4項は打ち消し合います。また、オイラー=ラグランジュ方程式より \partial L/\partial q = \dot{p} なので:

dH = \dot{q}dp - \dot{p}dq - \frac{\partial L}{\partial t}dt

両者を比較すると、ハミルトンの正準方程式が得られます:

\boxed{\dot{q} = \frac{\partial H}{\partial p}, \quad \dot{p} = -\frac{\partial H}{\partial q}}

また:

\frac{\partial H}{\partial t} = -\frac{\partial L}{\partial t}

方程式の特徴

  1. 1階の連立微分方程式:ラグランジュ形式では2階微分方程式でしたが、ハミルトン形式では1階の方程式が2倍の数だけあります。

  2. 対称的な構造:q と p の方程式は符号を除いて対称です。これは正準構造の表れです。

  3. 相空間:(q, p) を座標とする空間を相空間(phase space)と呼びます。系の状態は相空間の1点で表され、時間発展は相空間内の軌道として表現されます。

4.4 具体例:調和振動子

ラグランジアンからハミルトニアンへ

1次元調和振動子(バネ定数 k、質量 m)を考えます。

ラグランジアン:

L = \frac{1}{2}m\dot{x}^2 - \frac{1}{2}kx^2

一般化運動量:

p = \frac{\partial L}{\partial \dot{x}} = m\dot{x}

\dot{x} を p で表す:

\dot{x} = \frac{p}{m}

ハミルトニアン:

H = p\dot{x} - L = p \cdot \frac{p}{m} - \frac{1}{2}m\left(\frac{p}{m}\right)^2 + \frac{1}{2}kx^2
H = \frac{p^2}{m} - \frac{p^2}{2m} + \frac{1}{2}kx^2 = \frac{p^2}{2m} + \frac{1}{2}kx^2

これは確かに全エネルギー T + V です。

正準方程式

\dot{x} = \frac{\partial H}{\partial p} = \frac{p}{m}
\dot{p} = -\frac{\partial H}{\partial x} = -kx

第1式を時間微分して第2式を代入すると:

m\ddot{x} = \dot{p} = -kx

これは調和振動子の運動方程式です。

相空間での軌道

調和振動子の解は:

x(t) = A\cos(\omega t + \phi), \quad p(t) = -m\omega A\sin(\omega t + \phi)

ここで \omega = \sqrt{k/m} です。

相空間 (x, p) での軌道を見ると:

\frac{x^2}{A^2} + \frac{p^2}{(m\omega A)^2} = 1

これは楕円であり、エネルギー一定の軌道を表しています。

4.5 ポアソン括弧

ハミルトン力学の美しい構造を表す重要な概念がポアソン括弧(Poisson bracket)です。

定義

相空間上の2つの関数 f(q, p, t) と g(q, p, t) に対し、ポアソン括弧を次のように定義します:

\boxed{\{f, g\} = \frac{\partial f}{\partial q}\frac{\partial g}{\partial p} - \frac{\partial f}{\partial p}\frac{\partial g}{\partial q}}

多自由度系では:

\{f, g\} = \sum_i \left( \frac{\partial f}{\partial q_i}\frac{\partial g}{\partial p_i} - \frac{\partial f}{\partial p_i}\frac{\partial g}{\partial q_i} \right)

基本的な性質

  1. 反対称性:\{f, g\} = -\{g, f\}

  2. 双線形性:\{af + bg, h\} = a\{f, h\} + b\{g, h\}

  3. ライプニッツ則:\{fg, h\} = f\{g, h\} + g\{f, h\}

  4. ヤコビ恒等式:\{f, \{g, h\}\} + \{g, \{h, f\}\} + \{h, \{f, g\}\} = 0

基本ポアソン括弧

座標と運動量のポアソン括弧は:

\{q_i, q_j\} = 0, \quad \{p_i, p_j\} = 0, \quad \{q_i, p_j\} = \delta_{ij}

運動方程式のポアソン括弧表示

任意の物理量 f(q, p, t) の時間発展は:

\frac{df}{dt} = \{f, H\} + \frac{\partial f}{\partial t}

特に、q と p について:

\dot{q} = \{q, H\} = \frac{\partial H}{\partial p}
\dot{p} = \{p, H\} = -\frac{\partial H}{\partial q}

これはハミルトンの正準方程式そのものです!

保存量の条件

物理量 f が保存する(df/dt = 0)条件は:

\{f, H\} + \frac{\partial f}{\partial t} = 0

f が時間に陽に依存しない場合:

\{f, H\} = 0

つまり、f が保存量であることと、f が H とポアソン可換であることは同値です。

量子力学への橋渡し

ポアソン括弧の構造は、量子力学の交換子と深い関係があります。

量子力学では、物理量は演算子(オペレーター)で表され、交換子 [\hat{A}, \hat{B}] = \hat{A}\hat{B} - \hat{B}\hat{A} が定義されます。

ディラックは、古典力学と量子力学の対応関係(正準量子化)を:

\{f, g\}_{\text{classical}} \longrightarrow \frac{1}{i\hbar}[\hat{f}, \hat{g}]_{\text{quantum}}

という置き換えで表しました。特に基本交換関係:

[\hat{x}, \hat{p}] = i\hbar

は、基本ポアソン括弧 \{x, p\} = 1 に対応しています。


第5章:ネーターの定理——対称性と保存則

5.1 対称性と物理法則

物理学において、対称性は根本的な役割を果たします。対称性とは、ある変換に対して系(または法則)が不変であることです。

対称性の例

  • 空間並進対称性:物理法則は空間のどこでも同じ
  • 時間並進対称性:物理法則は時間によらない
  • 回転対称性:物理法則は方向によらない
  • ローレンツ対称性:物理法則は慣性系によらない(特殊相対性理論)
  • ゲージ対称性:電磁気学や素粒子物理学の基礎

エミー・ネーターは1918年に、対称性と保存則の深い関係を定理として証明しました。

5.2 ネーターの定理

定理の主張

連続的な対称性があれば、対応する保存量が存在する。

より正確には:

ラグランジアン L が連続変換 q \to q + \epsilon\delta q(\epsilon は無限小パラメータ)の下で不変(または全微分だけ変化)ならば、次の量 Q が保存する:

\boxed{Q = \frac{\partial L}{\partial \dot{q}}\delta q}

導出

ラグランジアン L(q, \dot{q}) を考えます。ここでは簡単のため、変換に対してラグランジアンそのものが不変(\delta L = 0)である場合を示します。変換 q \to q' = q + \epsilon\delta q(q) の下で、L の変化は:

\delta L = \frac{\partial L}{\partial q}\delta q + \frac{\partial L}{\partial \dot{q}}\delta \dot{q}

オイラー=ラグランジュ方程式より \partial L/\partial q = d/dt(\partial L/\partial \dot{q}) なので:

\delta L = \frac{d}{dt}\left(\frac{\partial L}{\partial \dot{q}}\right)\delta q + \frac{\partial L}{\partial \dot{q}}\frac{d(\delta q)}{dt} = \frac{d}{dt}\left(\frac{\partial L}{\partial \dot{q}}\delta q\right)

L が対称性を持つ、すなわち \delta L = 0 なら:

\frac{d}{dt}\left(\frac{\partial L}{\partial \dot{q}}\delta q\right) = 0

したがって Q = (\partial L/\partial \dot{q})\delta q = p\,\delta q が保存量です。

5.3 具体例

空間並進対称性 → 運動量保存

自由粒子を考えます:L = \frac{1}{2}m\dot{x}^2

並進変換 x \to x + \epsilon(\delta x = 1)の下で L は不変です(L は x に依存しない)。

保存量:

Q = \frac{\partial L}{\partial \dot{x}} \cdot 1 = m\dot{x} = p

運動量が保存します。

回転対称性 → 角運動量保存

中心力場:L = \frac{1}{2}m(\dot{r}^2 + r^2\dot{\theta}^2) - V(r)

回転変換 \theta \to \theta + \epsilon(\delta\theta = 1)の下で L は不変です。

保存量:

Q = \frac{\partial L}{\partial \dot{\theta}} \cdot 1 = mr^2\dot{\theta} = L_z

角運動量が保存します。

時間並進対称性 → エネルギー保存

L が時間 t に陽に依存しない場合を考えます。

厳密には、時間並進はラグランジアンの不変性ではなく作用の不変性として扱う必要がありますが、結果としてハミルトニアン H が保存量となります。

\frac{dH}{dt} = \frac{\partial H}{\partial t} = -\frac{\partial L}{\partial t} = 0

自然な系では H = E(全エネルギー)なので、エネルギー保存則が導かれます。

対称性と保存則の対応表

対称性 変換 保存量
空間の一様性 並進 \mathbf{r} \to \mathbf{r} + \mathbf{a} 運動量 \mathbf{p}
空間の等方性 回転 角運動量 \mathbf{L}
時間の一様性 時間並進 t \to t + s エネルギー E
ゲージ対称性 \psi \to e^{i\alpha}\psi 電荷

5.4 ネーターの定理の深い意味

ネーターの定理は、物理学の根本原理を明らかにしています:

  1. 保存則は対称性の帰結:エネルギー、運動量、角運動量などの保存則は、時空の対称性から自然に導かれます。

  2. 対称性が物理法則を決定する:現代物理学では、基礎理論を構築する際に、まず対称性を指定し、その対称性を満たすラグランジアンを構成するというアプローチが取られます。

  3. 破れた対称性の物理:対称性が「破れる」現象(自発的対称性の破れ)は、素粒子物理学で重要な役割を果たします(ヒッグス機構など)。


第6章:正準変換と発展した話題

6.1 正準変換

動機

ハミルトン力学では、適切な変数を選ぶことで問題が劇的に簡単になることがあります。では、どのような変数変換が「許される」のでしょうか?

(q, p) \to (Q, P) という変換を考えます。新しい変数でもハミルトンの正準方程式

\dot{Q} = \frac{\partial K}{\partial P}, \quad \dot{P} = -\frac{\partial K}{\partial Q}

が成り立つような変換を正準変換(canonical transformation)と呼びます。ここで K は新しいハミルトニアンです。

正準変換の条件

変換 (q, p) \to (Q, P) が正準変換である条件は、基本ポアソン括弧が保たれることです:

\{Q_i, Q_j\}_{q,p} = 0, \quad \{P_i, P_j\}_{q,p} = 0, \quad \{Q_i, P_j\}_{q,p} = \delta_{ij}

母関数

正準変換は母関数(generating function)を用いて系統的に構成できます。

例えば、母関数 F_2(q, P, t) を使うと:

p = \frac{\partial F_2}{\partial q}, \quad Q = \frac{\partial F_2}{\partial P}, \quad K = H + \frac{\partial F_2}{\partial t}

という関係から正準変換が定まります。

例:調和振動子の作用-角変数

調和振動子 H = \frac{p^2}{2m} + \frac{1}{2}kx^2 を考えます。

正準変換により、新しい変数 (J, \phi)(作用変数と角変数)に移ると:

K = \omega J

という非常に簡単な形になります。ここで J はエネルギーを角振動数で割ったもの、\phi は位相角です。

正準方程式は:

\dot{J} = 0, \quad \dot{\phi} = \omega

J は保存し(これは当然、エネルギー保存の言い換え)、\phi は等速で増加します。

6.2 ハミルトン=ヤコビ理論

ハミルトン=ヤコビ方程式

解析力学の究極の目標の1つは、運動を積分(求積)することです。ハミルトン=ヤコビ理論は、この目標に向けた強力な方法を提供します。

母関数 S(q, P, t) を使って、新しいハミルトニアン K = 0 となる正準変換を見つけることを考えます。このとき:

H\left(q, \frac{\partial S}{\partial q}, t\right) + \frac{\partial S}{\partial t} = 0

これがハミルトン=ヤコビ方程式です。

物理的意味

ハミルトン=ヤコビ方程式の解 S は「主関数」と呼ばれ、作用に相当します。また、S の等位面は波面を定義し、粒子の軌道はこの波面に直交します。

この構造は、量子力学の波動関数の位相と深く関係しています。実際、シュレディンガーはハミルトン=ヤコビ理論を手がかりにして波動方程式を発見しました。

6.3 リウヴィルの定理

相空間の体積保存

ハミルトン力学には、統計力学につながる重要な性質があります。

リウヴィルの定理:正準方程式に従う時間発展では、相空間の体積が保存される。

より正確には、相空間内の任意の領域を考え、その領域内の点が正準方程式に従って運動するとき、領域の体積(\int dq\,dp)は時間によらず一定です。

証明の概略

1自由度系で考えます。相空間の「速度場」は (\dot{q}, \dot{p}) = (\partial H/\partial p, -\partial H/\partial q) です。

この速度場の発散を計算すると:

\frac{\partial \dot{q}}{\partial q} + \frac{\partial \dot{p}}{\partial p} = \frac{\partial^2 H}{\partial q \partial p} - \frac{\partial^2 H}{\partial p \partial q} = 0

発散がゼロということは、流れが非圧縮性であること、つまり体積が保存されることを意味します。

統計力学への影響

リウヴィルの定理は、統計力学の基礎となります。多数の粒子からなる系の状態を、相空間内の点の集合(アンサンブル)として考えるとき、この集合の「密度」がどう時間発展するかを記述するリウヴィル方程式が導かれます。


第7章:拡張と応用

7.1 拘束系の力学

ホロノミック拘束とノンホロノミック拘束

ホロノミック拘束は、座標のみの関係式で表される拘束です:

f(q_1, q_2, \ldots, q_n, t) = 0

例:振り子の紐の長さ一定 x^2 + y^2 = l^2

ノンホロノミック拘束は、速度も含む関係式で、積分して座標だけの関係にできない拘束です:

\sum_i a_i(q, t)\dot{q}_i + b(q, t) = 0

例:滑らずに転がる球(接点での速度条件)

ホロノミック拘束は一般化座標の選び方で自動的に満たせますが、ノンホロノミック拘束はラグランジュの未定乗数法などで扱う必要があります。

ダランベールの原理

仮想仕事の原理:系が平衡にあるとき、拘束に適合する仮想変位に対して、力がなす仮想仕事の総和はゼロです。

\sum_i \mathbf{F}_i \cdot \delta\mathbf{r}_i = 0

ダランベールの原理:これを動的な場合に拡張したものです。

\sum_i (\mathbf{F}_i - m_i\ddot{\mathbf{r}}_i) \cdot \delta\mathbf{r}_i = 0

この原理から、オイラー=ラグランジュ方程式を導出できます。

7.2 剛体の力学

剛体の自由度

空間内の剛体は6自由度を持ちます(重心位置3 + 姿勢3)。

姿勢を表すには、オイラー角 (\phi, \theta, \psi) がよく使われます。

剛体のラグランジアン

剛体の運動エネルギーは:

T = \frac{1}{2}M|\dot{\mathbf{R}}|^2 + \frac{1}{2}\boldsymbol{\omega} \cdot \mathbf{I} \cdot \boldsymbol{\omega}

第1項は重心の並進運動、第2項は重心まわりの回転運動です。\mathbf{I} は慣性テンソルです。

オイラー角を使ってラグランジアンを書くと、ジャイロスコープや歳差運動などの複雑な現象を系統的に扱えます。

7.3 連続体と場の理論

離散系から連続系へ

N 個の質点からなるバネ鎖を考えます。N \to \infty、間隔 \to 0 の極限で、系は連続体(弦)になります。

離散系のラグランジアン

L = \sum_i \left(\frac{1}{2}m\dot{q}_i^2 - \frac{1}{2}k(q_{i+1} - q_i)^2\right)

は、連続極限で

L = \int dx \, \mathcal{L}\left(\phi, \frac{\partial\phi}{\partial t}, \frac{\partial\phi}{\partial x}\right)

というラグランジアン密度 \mathcal{L} の積分になります。

場の方程式

オイラー=ラグランジュ方程式も一般化されます。場 \phi(x, t) に対して:

\frac{\partial}{\partial t}\frac{\partial\mathcal{L}}{\partial(\partial_t\phi)} + \frac{\partial}{\partial x}\frac{\partial\mathcal{L}}{\partial(\partial_x\phi)} - \frac{\partial\mathcal{L}}{\partial\phi} = 0

弦の場合、これは波動方程式

\frac{\partial^2\phi}{\partial t^2} = c^2\frac{\partial^2\phi}{\partial x^2}

になります。

電磁気学と相対論

マクスウェル方程式も、ラグランジアン密度

\mathcal{L} = -\frac{1}{4}F_{\mu\nu}F^{\mu\nu}

から導かれます。ここで F_{\mu\nu} は電磁場テンソルです。

一般相対性理論も、アインシュタイン=ヒルベルト作用

S = \int d^4x \sqrt{-g} R

から導かれます。R は時空の曲率を表すリッチスカラー、g は時空の計量(距離の測り方)を表す計量テンソルの行列式です。

7.4 量子力学への発展

正準量子化

ハミルトン形式は、量子力学への自然な橋渡しを提供します。

正準量子化の手順:

  1. 座標と運動量を演算子 \hat{q}, \hat{p} に置き換える
  2. ポアソン括弧を交換子に置き換える:\{q, p\} = 1 \to [\hat{q}, \hat{p}] = i\hbar
  3. ハミルトニアンを演算子 \hat{H} に置き換える

シュレディンガー方程式

i\hbar\frac{\partial}{\partial t}|\psi\rangle = \hat{H}|\psi\rangle

は、古典的なハミルトンの運動方程式の量子版と見なせます。

経路積分

ファインマンは、量子力学の別の定式化として経路積分を提唱しました。

粒子が点 A から点 B に行く確率振幅は:

\langle B|A \rangle = \int \mathcal{D}[q(t)] \, e^{iS[q]/\hbar}

すべての可能な経路 q(t) について「作用の位相」を足し合わせます。古典極限(\hbar \to 0)では、作用を停留させる経路(古典的軌道)が主に寄与し、最小作用の原理が再現されます。

この定式化は、作用 S と変分原理が量子力学でも本質的な役割を果たすことを示しています。


第8章:解析力学を学ぶ意義

8.1 なぜ解析力学を学ぶのか

計算の効率化

複雑な拘束を持つ系や、曲線座標を使う問題では、解析力学は圧倒的に効率的です。ラグランジアンを書き下ろし、微分するだけで運動方程式が得られます。

概念的な理解

「力」という直感的だが座標に依存する概念から離れ、「エネルギー」という座標に依存しないスカラー量を基礎とすることで、物理現象の本質に迫れます。

現代物理学の言語

  • 量子力学:正準量子化、経路積分
  • 統計力学:相空間、リウヴィル方程式、エルゴード理論
  • 相対論的場の理論:ラグランジアン密度、ゲージ対称性
  • 一般相対性理論:作用原理、測地線方程式

現代物理学のほぼすべての分野が、解析力学の言葉で記述されています。

対称性の視点

ネーターの定理を通じて、対称性と保存則の関係を理解できます。これは現代物理学の根幹をなす考え方です。

8.2 学習のアドバイス

段階的な学習

  1. まずラグランジュ形式:変分原理、一般化座標、オイラー=ラグランジュ方程式
  2. 次にハミルトン形式:相空間、正準方程式、ポアソン括弧
  3. さらに正準変換とネーターの定理:対称性と保存則

具体例を大切に

抽象的な定式化だけでなく、具体的な問題(振り子、中心力場、剛体など)を自分の手で解くことが重要です。

数学的な準備

  • 多変数の微積分
  • 線形代数
  • 常微分方程式
  • (できれば)変分法の基礎

おすすめの教科書

  • ランダウ=リフシッツ『力学』:簡潔で本質的
  • ゴールドスタイン『古典力学』:詳細で網羅的
  • アーノルド『古典力学の数学的方法』:数学的に厳密

まとめ

解析力学は、ニュートン力学を再定式化したものですが、単なる「言い換え」ではありません。

  • 変分原理という新しい視点を提供し
  • 一般化座標で拘束を自然に扱い
  • 対称性と保存則の関係(ネーターの定理)を明らかにし
  • ポアソン括弧という代数構造を通じて量子力学への橋渡しをし
  • 相空間の概念で統計力学の基礎を与え
  • 場の理論への拡張を可能にする

現代物理学を学ぶ上で、解析力学は避けて通れない基礎です。最初は抽象的に感じるかもしれませんが、具体的な問題を解きながら学ぶことで、その美しさと有用性が実感できるはずです。

物理学の深い構造を理解する旅の第一歩として、解析力学をぜひ楽しんでください。


付録:主要な公式集

ラグランジュ形式

ラグランジアン

L = T - V

オイラー=ラグランジュ方程式

\frac{d}{dt}\frac{\partial L}{\partial \dot{q}_i} - \frac{\partial L}{\partial q_i} = 0

一般化運動量

p_i = \frac{\partial L}{\partial \dot{q}_i}

ハミルトン形式

ハミルトニアン

H = \sum_i p_i\dot{q}_i - L

正準方程式

\dot{q}_i = \frac{\partial H}{\partial p_i}, \quad \dot{p}_i = -\frac{\partial H}{\partial q_i}

ポアソン括弧

\{f, g\} = \sum_i \left(\frac{\partial f}{\partial q_i}\frac{\partial g}{\partial p_i} - \frac{\partial f}{\partial p_i}\frac{\partial g}{\partial q_i}\right)

基本ポアソン括弧

\{q_i, p_j\} = \delta_{ij}, \quad \{q_i, q_j\} = 0, \quad \{p_i, p_j\} = 0

運動方程式(ポアソン括弧形式)

\frac{df}{dt} = \{f, H\} + \frac{\partial f}{\partial t}

ネーターの定理

対称性変換 q \to q + \epsilon\delta q でラグランジアンが不変なら、

Q = \frac{\partial L}{\partial \dot{q}}\delta q

が保存量。

Discussion