解析力学から様々な理論分野を組み立てる 第1部 第5章 §1.5.4 ★★★ 一様性・等方性で L = L ( v 2 ) L=L(v^2) L = L ( v 2 ) 、Galilei 変換で L L L が全微分だけ変わることを要求すると L = 1 2 m v 2 L=\frac12mv^2 L = 2 1 m v 2 。相対論では同じ手順で − m c 2 1 − v 2 / c 2 -mc^2\sqrt{1-v^2/c^2} − m c 2 1 − v 2 / c 2
一様性と等方性
1つの自由粒子の L L L を、運動の法則を前提にせず対称性だけから決める。慣性系で次の3つを要求する。
(i) 時間の一様性:L L L は時刻 t t t をあらわに含まない。
(ii) 空間の一様性:L L L は位置 r \bm{r} r をあらわに含まない。
(iii) 空間の等方性:L L L は速度 v = r ˙ \bm{v}=\dot{\bm{r}} v = r ˙ の向きによらない。
(i)(ii)(iii) から L L L は速さの2乗 v 2 = ∣ v ∣ 2 v^2=\lvert \bm{v}\rvert ^2 v 2 = ∣ v ∣ 2 だけの関数 L = L ( v 2 ) L=L(v^2) L = L ( v 2 ) になる。∂ L / ∂ r = 0 \partial L/\partial\bm{r}=0 ∂ L / ∂ r = 0 なので Euler–Lagrange 方程式は d d t ∂ L ∂ v = 0 \frac{\mathrm{d} }{\mathrm{d} t}\frac{\partial L}{\partial \bm{v}}=0 d t d ∂ v ∂ L = 0 で、∂ L / ∂ v = 2 L ′ ( v 2 ) v \partial L/\partial\bm{v}=2L'(v^2)\bm{v} ∂ L / ∂ v = 2 L ′ ( v 2 ) v が一定になる。L ′ L' L ′ が 0 0 0 でなければ v \bm{v} v は一定で、自由粒子は等速直線運動をする。慣性の法則は (i)〜(iii) から出る。しかし L ( v 2 ) L(v^2) L ( v 2 ) の関数形は決まらない。
Galilei 変換で全微分だけ変わること
ここで Galilei の相対性を要求する。速度 − u -\bm{u} − u で動く別の慣性系から見ると、粒子の速度は v ′ = v + u \bm{v}'=\bm{v}+\bm{u} v ′ = v + u になる。2つの系で同じ運動方程式が出ればよいので、命題 1.5.3 により、L ( v ′ 2 ) L(v'^2) L ( v ′2 ) と L ( v 2 ) L(v^2) L ( v 2 ) の差が時間の全微分 d d t Λ ( r , t ) \frac{\mathrm{d} }{\mathrm{d} t}\Lambda(\bm{r},t) d t d Λ ( r , t ) であればよい。u \bm{u} u が小さいとして1次まで展開すると
L ( ∣ v + u ∣ 2 ) = L ( v 2 ) + 2 L ′ ( v 2 ) v ⋅ u + O ( u 2 ) L\bigl(\lvert \bm{v}+\bm{u}\rvert ^2\bigr)=L(v^2)+2L'(v^2)\,\bm{v}\cdot\bm{u}+\mathcal{O}(u^2) L ( ∣ v + u ∣ 2 ) = L ( v 2 ) + 2 L ′ ( v 2 ) v ⋅ u + O ( u 2 )
となる。d d t Λ ( r , t ) = ∇ Λ ⋅ v + ∂ t Λ \frac{\mathrm{d} }{\mathrm{d} t}\Lambda(\bm{r},t)=\nabla\Lambda\cdot\bm{v}+\partial_t\Lambda d t d Λ ( r , t ) = ∇Λ ⋅ v + ∂ t Λ は v \bm{v} v の1次式で、係数 ∇ Λ \nabla\Lambda ∇Λ は r \bm{r} r と t t t だけの関数である。2 L ′ ( v 2 ) v ⋅ u 2L'(v^2)\,\bm{v}\cdot\bm{u} 2 L ′ ( v 2 ) v ⋅ u がこの形になるには、L ′ ( v 2 ) L'(v^2) L ′ ( v 2 ) が v \bm{v} v によらない定数でなければならない。
命題 1.5.4(自由粒子の Lagrangian ★★★)
(i)〜(iii) と、無限小の Galilei 変換で L L L が時間の全微分だけ変わることを要求すると、
L = 1 2 m v 2 ( m は定数 ) (1.5.7) L=\frac12m\,v^2\qquad(m\text{ は定数})\tag{1.5.7} L = 2 1 m v 2 ( m は定数 ) ( 1.5.7 )
となる(加える定数は運動に関係しないので落とした)。
係数を 1 2 m \frac12m 2 1 m と書いたのは、この定数を質量と呼ぶための約束である。有限の u \bm{u} u でも
1 2 m ∣ v + u ∣ 2 = 1 2 m v 2 + d d t ( m u ⋅ r + 1 2 m u 2 t ) \frac12m\lvert \bm{v}+\bm{u}\rvert ^2=\frac12mv^2+\frac{\mathrm{d} }{\mathrm{d} t}\Bigl(m\bm{u}\cdot\bm{r}+\frac12mu^2t\Bigr) 2 1 m ∣ v + u ∣ 2 = 2 1 m v 2 + d t d ( m u ⋅ r + 2 1 m u 2 t )
となり、要求は満たされている。m > 0 m>0 m > 0 は対称性からは出ない。m < 0 m<0 m < 0 だと、端を固定した経路を速く振動させるほど 1 2 m ∫ η ˙ 2 d t \frac12m\int\dot\eta^2\,\mathrm{d} t 2 1 m ∫ η ˙ 2 d t がいくらでも負になり、作用に下限が無くなる。1.5.5 節で見るように、真の経路は短い時間では作用の最小点であり、そのために m > 0 m>0 m > 0 が要る。
相互作用する N N N 個の粒子では、各粒子の 1 2 m a v a 2 \frac12m_av_a^2 2 1 m a v a 2 を足し、相互作用を表す項を加える。相互作用の項がどんな形かは対称性からは決まらない。速度によらない相互作用なら、それを − V ( r 1 , … , r N ) -V(\bm{r}_1,\dots,\bm{r}_N) − V ( r 1 , … , r N ) と書き、V V V が静力学で測る位置エネルギーと一致するように選ぶ。V V V が粒子の相対位置 r a − r b \bm{r}_a-\bm{r}_b r a − r b だけの関数なら、Galilei 変換で V V V は変わらない。こうして式 (1.5.1) が得られる。速度による相互作用(電磁場中の荷電粒子)では − V -V − V の代わりに − q ϕ + q v ⋅ A -q\phi+q\bm{v}\cdot\bm{A} − q ϕ + q v ⋅ A が入る(第2部第2章 )。
運動エネルギーの部分 ∑ a 1 2 m a v a 2 \sum_a\frac12m_av_a^2 ∑ a 2 1 m a v a 2 は、どの系でも同じ形をしている。これは自由な空間の対称性(時刻と位置の並進・回転・Galilei 変換)だけから決めたので、それらの対称性をすべて持っている。したがって、系ごとに違う物理はすべて相互作用の項 V V V に入り、どの保存量が残るかも V V V の対称性だけで決まる。
V V V の性質L L L を変えない変換保存する量 座標 y y y を含まない y y y 方向の並進p y p_y p y 原点からの距離 r r r だけによる 原点のまわりの回転 角運動量 L \bm{L} L 位置の差 r 1 − r 2 \bm{r}_1-\bm{r}_2 r 1 − r 2 だけによる 全体の並進 全運動量 p 1 + p 2 \bm{p}_1+\bm{p}_2 p 1 + p 2 時刻 t t t を含まない 時刻の並進 エネルギー
1行目の y y y のように L L L に現れない座標を循環座標 と呼ぶ。対称性から保存量を出す一般の定理(Noether の定理)は第1部第10章 で示す。
相対論では同じ手順が別の答えを出す
Galilei 変換の代わりに Lorentz 変換で L d t L\,\mathrm{d} t L d t が変わらないことを要求すると、自由粒子の L L L は L = − m c 2 1 − v 2 / c 2 L=-mc^2\sqrt{1-v^2/c^2} L = − m c 2 1 − v 2 / c 2 になる(第2部第3章 )。v ≪ c v\ll c v ≪ c で展開すると − m c 2 + 1 2 m v 2 + m v 4 8 c 2 + ⋯ -mc^2+\frac12mv^2+\frac{mv^4}{8c^2}+\cdots − m c 2 + 2 1 m v 2 + 8 c 2 m v 4 + ⋯ で、定数 − m c 2 -mc^2 − m c 2 は運動に関係せず、v 4 v^4 v 4 以降の項を落としたものが式 (1.5.7) である。
例題 1.5.4 Galilei 変換と運動エネルギーの形 ★★
オリジナル 1次元の自由粒子について次に答えよ。
(1) x ′ = x + u t x'=x+ut x ′ = x + u t (u u u は定数)として、1 2 m x ˙ ′ 2 \frac12m\dot x'^2 2 1 m x ˙ ′2 と 1 2 m x ˙ 2 \frac12m\dot x^2 2 1 m x ˙ 2 の差を、ある関数 Λ ( x , t ) \Lambda(x,t) Λ ( x , t ) の時間の全微分の形に書け。
(2) κ x ˙ 4 \kappa\dot x^4 κ x ˙ 4 を同じ変換で書き直したとき、u u u の1次の変化が Λ ( x , t ) \Lambda(x,t) Λ ( x , t ) の全微分の形に書けないことを示せ。
(3) L = − m c 2 1 − x ˙ 2 / c 2 L=-mc^2\sqrt{1-\dot x^2/c^2} L = − m c 2 1 − x ˙ 2 / c 2 を x ˙ \dot x x ˙ について4次まで展開せよ。地球の公転の速さ v = 3.0 × 10 4 m s − 1 v=3.0\times10^4\,\mathrm{m\,s^{-1}} v = 3.0 × 1 0 4 m s − 1 で、v 4 v^4 v 4 の項と v 2 v^2 v 2 の項の比を求めよ(c = 2.998 × 10 8 m s − 1 c=2.998\times10^8\,\mathrm{m\,s^{-1}} c = 2.998 × 1 0 8 m s − 1 )。
解答 (1) x ˙ ′ = x ˙ + u \dot x'=\dot x+u x ˙ ′ = x ˙ + u なので
1 2 m ( x ˙ + u ) 2 − 1 2 m x ˙ 2 = m u x ˙ + 1 2 m u 2 = d d t ( m u x + 1 2 m u 2 t ) . \frac12m(\dot x+u)^2-\frac12m\dot x^2=mu\dot x+\frac12mu^2=\frac{\mathrm{d} }{\mathrm{d} t}\Bigl(mux+\frac12mu^2t\Bigr). 2 1 m ( x ˙ + u ) 2 − 2 1 m x ˙ 2 = m u x ˙ + 2 1 m u 2 = d t d ( m ux + 2 1 m u 2 t ) .
Λ = m u x + 1 2 m u 2 t \Lambda=mux+\frac12mu^2t Λ = m ux + 2 1 m u 2 t である。
(2) κ ( x ˙ + u ) 4 − κ x ˙ 4 = 4 κ x ˙ 3 u + O ( u 2 ) \kappa(\dot x+u)^4-\kappa\dot x^4=4\kappa\dot x^3u+\mathcal{O}(u^2) κ ( x ˙ + u ) 4 − κ x ˙ 4 = 4 κ x ˙ 3 u + O ( u 2 ) 。一方 d d t Λ ( x , t ) = ∂ x Λ x ˙ + ∂ t Λ \frac{\mathrm{d} }{\mathrm{d} t}\Lambda(x,t)=\partial_x\Lambda\,\dot x+\partial_t\Lambda d t d Λ ( x , t ) = ∂ x Λ x ˙ + ∂ t Λ は x ˙ \dot x x ˙ の1次式である。4 κ x ˙ 3 u 4\kappa\dot x^3u 4 κ x ˙ 3 u は x ˙ \dot x x ˙ の3次式なので、どんな Λ ( x , t ) \Lambda(x,t) Λ ( x , t ) を選んでもこの形にならない。κ x ˙ 4 \kappa\dot x^4 κ x ˙ 4 を Lagrangian にすると、慣性系を取り替えたときに運動方程式の形が変わる。
(3) 1 − s = 1 − s 2 − s 2 8 − ⋯ \sqrt{1-s}=1-\frac s2-\frac{s^2}{8}-\cdots 1 − s = 1 − 2 s − 8 s 2 − ⋯ に s = x ˙ 2 / c 2 s=\dot x^2/c^2 s = x ˙ 2 / c 2 を入れて
L = − m c 2 + 1 2 m x ˙ 2 + m x ˙ 4 8 c 2 + O ( x ˙ 6 ) . L=-mc^2+\frac12m\dot x^2+\frac{m\dot x^4}{8c^2}+\mathcal{O}(\dot x^6). L = − m c 2 + 2 1 m x ˙ 2 + 8 c 2 m x ˙ 4 + O ( x ˙ 6 ) .
v 4 v^4 v 4 の項と v 2 v^2 v 2 の項の比は m v 4 / ( 8 c 2 ) m v 2 / 2 = v 2 4 c 2 = ( 3.0 × 10 4 ) 2 4 × ( 2.998 × 10 8 ) 2 = 2.5 × 10 − 9 \dfrac{mv^4/(8c^2)}{mv^2/2}=\dfrac{v^2}{4c^2}=\dfrac{(3.0\times10^4)^2}{4\times(2.998\times10^8)^2}=2.5\times10^{-9} m v 2 /2 m v 4 / ( 8 c 2 ) = 4 c 2 v 2 = 4 × ( 2.998 × 1 0 8 ) 2 ( 3.0 × 1 0 4 ) 2 = 2.5 × 1 0 − 9 である。地球の公転の速さでも、Galilei 変換から決めた 1 2 m v 2 \frac12mv^2 2 1 m v 2 との違いは 10 − 9 10^{-9} 1 0 − 9 の程度にとどまる。
近似の裏付け
型 B(極限) | 小さいもの v / c → 0 v/c\to0 v / c → 0 (速さが光速に比べて小さい) | 証明 証明済み(有限の時間区間で、力が位置について Lipschitz 連続なとき。運動方程式の解のパラメータへの連続依存) | 行き先 近似の数理 第2章
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com