tecktom.com

この章で選ぶ LL:運動エネルギー引くポテンシャル

2026年10月15日9分で読めますAI生成教科書『虎の巻』第1部 第5章 L=T-V の由来
目次
  1. 選ぶ
  2. 作用の値は何を測っているか
  3. 木下の層

AI 生成・木下が未勉強

解析力学から様々な理論分野を組み立てる 第1部 第5章 §1.5.1 ★★★ NN 個の質点の L=∑a12mar˙a2−VL=\sum_a\frac12m_a\dot{\bm{r}}_a^2-V を書き、自由落下で作用の値と中身を手で計算する

前の2章では、作用 S[q]=∫L dtS[q]=\int L\,\mathrm{d} t の LL を「与えられた関数」とだけ考え、停留の条件が Euler–Lagrange 方程式になることを数学として示した。この章で初めて LL に物理を入れる。粒子の力学で選ぶのは

L=T−V(運動エネルギー引くポテンシャル)L=T-V\qquad(\text{運動エネルギー引くポテンシャル})

で、この1行は第2部以降のすべての章の出発点になる。ただ、多くの教科書はこの式を天下りに与えるので、「なぜ足し算でなく引き算なのか」「なぜ運動エネルギーは速度の2乗なのか」「なぜ作用を最小にするのか」が分からないまま先へ進むことになる。この章ではこの3つの問いを別々に扱う。

答えを先に書くと、3つの問いの根拠はそれぞれ違う。運動エネルギーの形 12mv2\frac12mv^2 は、空間と時間の一様性・空間の等方性・Galilei 変換で物理が変わらないこと、という対称性の要請からほぼ一意に決まる。−V-V の符号と形は対称性からは出ず、Euler–Lagrange 方程式が Newton の運動方程式 mr¨=−∇Vm\ddot{\bm{r}}=-\nabla V と一致するように選んだものである。最後に、作用の停留という原理そのものは古典力学の中では説明されず、量子力学の経路積分で eiS/ℏ\mathrm{e}^{\mathrm{i} S/\hbar} の位相がそろう条件として現れる。このとき実際に要るのは停留であって最小ではないことも、調和振動子で確かめる。道具としては、Euler–Lagrange 方程式が座標の取り替えで形を変えないこと、LL に定数を掛けたり時間の全微分を足したりしても運動方程式が変わらないこと、の2つを示して使う。

選ぶ LL

質量 mam_a の質点が NN 個あり(a=1,…,Na=1,\dots,N)、位置を ra\bm{r}_a とする。力がポテンシャル V(r1,…,rN,t)V(\bm{r}_1,\dots,\bm{r}_N,t) から −∇aV-\nabla_aV として出る場合を考える。

定義 1.5.1(粒子の力学の Lagrangian)

L=T−V,T=∑a=1N12ma∣r˙a∣2,V=V(r1,…,rN,t)(1.5.1)L=T-V,\qquad T=\sum_{a=1}^N\frac12m_a\lvert \dot{\bm{r}}_a\rvert ^2,\qquad V=V(\bm{r}_1,\dots,\bm{r}_N,t)\tag{1.5.1}

を粒子の力学の Lagrangian とする。作用は S[q]=∫t1t2L dtS[q]=\int_{t_1}^{t_2}L\,\mathrm{d} t である。

TT は速度だけの関数、VV は位置(と時刻)だけの関数である。直交座標の成分 (x1,y1,z1,…,xN,yN,zN)(x_1,y_1,z_1,\dots,x_N,y_N,z_N) を一般化座標 q1,…,qfq_1,\dots,q_f(f=3Nf=3N)とすれば、第1部第4章の Euler–Lagrange 方程式

δSδqi(t)=∂L∂qi−ddt∂L∂q˙i=0(i=1,…,f)(1.5.2)\frac{\delta S}{\delta q_i(t)}=\frac{\partial L}{\partial q_i}-\frac{\mathrm{d} }{\mathrm{d} t}\frac{\partial L}{\partial \dot q_i}=0\qquad(i=1,\dots,f)\tag{1.5.2}

がそのまま当てはまる。∂L/∂x˙a=max˙a\partial L/\partial\dot x_a=m_a\dot x_a、∂L/∂xa=−∂V/∂xa\partial L/\partial x_a=-\partial V/\partial x_a なので、式 (1.5.2) は

mar¨a=−∇aV(a=1,…,N)(1.5.3)m_a\ddot{\bm{r}}_a=-\nabla_aV\qquad(a=1,\dots,N)\tag{1.5.3}

となる。これが Newton の運動方程式で、式 (1.5.1) の LL は、停留させると Newton の運動方程式が出るように選んだものになっている。ただしこの確認だけでは「同じ方程式を出す LL はほかに無いのか」「なぜ TT と VV の差なのか」には答えていない。それを1.5.3節と1.5.4節で調べる。

作用の値は何を測っているか

T−VT-V の意味を、いちばん簡単な例で見ておく。鉛直上向きを zz とし、一様な重力 V=mgzV=mgz の中の質点を考える。L=12mz˙2−mgzL=\frac12m\dot z^2-mgz から z¨=−g\ddot z=-g が出る。時刻 00 に z=0z=0 から投げ上げ、時刻 TT にふたたび z=0z=0 に戻る経路は

zcl(t)=g2 t(T−t)(1.5.4)z_{\mathrm{cl}}(t)=\frac g2\,t(T-t)\tag{1.5.4}

で(添字 cl は古典的な経路の意)、最高点の高さは gT2/8gT^2/8 である。同じ投げ上げは第1部第1章の例題 1.1.1 で、真の経路をずらしても作用が1次で変わらないこと(停留)を見るのに使い、下の作用の値 S[zcl]=−mg2T3/24S[z_{\mathrm{cl}}]=-mg^2T^3/24 もそこで出ている。ここでは同じ経路で TT と VV の兼ね合いを見る。この経路で運動エネルギーと位置エネルギーの時間平均を取ると

⟨T⟩=mg2T224,⟨V⟩=mg2T212,S[zcl]=∫0T(T−V) dt=−mg2T324\langle T \rangle =\frac{mg^2T^2}{24},\qquad \langle V \rangle =\frac{mg^2T^2}{12},\qquad S[z_{\mathrm{cl}}]=\int_0^T(T-V)\,\mathrm{d} t=-\frac{mg^2T^3}{24}

となる(この節では TT が運動エネルギーと時間の長さの2つの意味で出る。時間の長さは積分の上端と、結果の式の T2T^2, T3T^3 にだけ現れる)。この値は 第1部第3章の例題で数学として求めた ∫(12q˙2+cq) dt\int(\frac12\dot q^2+cq)\,\mathrm{d} t の停留値 −c2/24-c^2/24 と同じ形で、c=gc=g、m=1m=1、T=1T=1 とおけば一致する。

作用 S=∫(T−V) dtS=\int(T-V)\,\mathrm{d} t が小さくなるのは、VV が大きい所(高い所)に長くいて、しかも TT が小さい(ゆっくり動く)経路である。ところが、決まった時間 TT のうちに高い所まで行って戻ってくるには速く動く必要があり、TT が大きくなる。真の経路 (1.5.4) は、高く上がることと速く動かないことの兼ね合いで決まっている。次の例題でこれを数で確かめる。

例題 1.5.1 放物線の族で作用を比べる ★★★

オリジナル質量 mm の質点が一様な重力の中で、時刻 00 と時刻 TT に高さ 00 にある。高さ AA をパラメータとする経路の族

zA(t)=4A t(T−t)T2(中点 t=T/2 で高さ A)z_A(t)=\frac{4A\,t(T-t)}{T^2}\qquad(\text{中点 }t=T/2\text{ で高さ }A)

を考える。

  1. (1) 作用 S(A)=∫0T(12mz˙A2−mgzA)dtS(A)=\int_0^T\bigl(\frac12m\dot z_A^2-mgz_A\bigr)\mathrm{d} t を AA の関数として求めよ。
  2. (2) S(A)S(A) を最小にする AA と、そのときの SS を求め、真の経路 (1.5.4) と比べよ。
  3. (3) m=0.1 kgm=0.1\,\mathrm{kg}、g=9.81 m s−2g=9.81\,\mathrm{m\,s^{-2}}、T=2 sT=2\,\mathrm{s} で、(2) の AA と SS の値を求めよ。
解答

(1) z˙A=4A(T−2t)/T2\dot z_A=4A(T-2t)/T^2 なので

∫0Tz˙A2 dt=16A2T4∫0T(T−2t)2 dt=16A2T4⋅T33=16A23T,∫0TzA dt=4AT2⋅T36=2AT3.\int_0^T\dot z_A^2\,\mathrm{d} t=\frac{16A^2}{T^4}\int_0^T(T-2t)^2\,\mathrm{d} t=\frac{16A^2}{T^4}\cdot\frac{T^3}{3}=\frac{16A^2}{3T},\qquad \int_0^Tz_A\,\mathrm{d} t=\frac{4A}{T^2}\cdot\frac{T^3}{6}=\frac{2AT}{3}.

よって

S(A)=8mA23T−2mgAT3.S(A)=\frac{8mA^2}{3T}-\frac{2mgAT}{3}.

第1項(運動エネルギーの寄与)は AA の2乗で増え、第2項(位置エネルギーの寄与、符号は負)は AA に比例して減る。

(2) dSdA=16mA3T−2mgT3=0\frac{\mathrm{d} S}{\mathrm{d} A}=\frac{16mA}{3T}-\frac{2mgT}{3}=0 から A=gT28A=\frac{gT^2}{8}。このとき

S=8m3T⋅g2T464−2mgT3⋅gT28=mg2T324−mg2T312=−mg2T324.S=\frac{8m}{3T}\cdot\frac{g^2T^4}{64}-\frac{2mgT}{3}\cdot\frac{gT^2}{8}=\frac{mg^2T^3}{24}-\frac{mg^2T^3}{12}=-\frac{mg^2T^3}{24}.

A=gT2/8A=gT^2/8 は真の経路の最高点の高さで、zAz_A はこのとき式 (1.5.4) と一致する。族が真の経路を含んでいるので、族の中の最小は真の作用と同じ値になる。

(3) A=9.81×4/8=4.905 mA=9.81\times4/8=4.905\,\mathrm{m}、S=−0.1×9.812×8/24=−3.21 J sS=-0.1\times9.81^2\times8/24=-3.21\,\mathrm{J\,s}。AA をこれより大きくすると運動エネルギーの増え方が、小さくすると位置エネルギーの減り方の不足が上回り、どちらでも SS は大きくなる(図 1.5.1)。

図

図 1.5.1 一様な重力の中で端を z(0)=z(T)=0z(0)=z(T)=0 に固定した経路の族の作用(m=0.1 kgm=0.1\,\mathrm{kg}、T=2 sT=2\,\mathrm{s})。点線の放物線の族は真の経路を含むので、最小値 −3.21 J s-3.21\,\mathrm{J\,s}(黒丸)が真の作用である。実線の正弦の族は真の経路を含まない近似で、最小値 −3.16 J s-3.16\,\mathrm{J\,s} は真の値より上にある

図 1.5.1 の実線は、経路を Asin⁡(πt/T)A\sin(\pi t/T) の形に限って最小を探したもので、最小は A=4gT2/π3A=4gT^2/\pi^3、S=−4mg2T3/π4S=-4mg^2T^3/\pi^4 になる。真の値との比は (4/π4)/(1/24)=96/π4≈0.986(4/\pi^4)/(1/24)=96/\pi^4\approx0.986 で、近似の作用は真の作用より 1.4%1.4\% だけ上にある。この LL では真の経路のまわりの2次の変分が 12m∫η˙2 dt≥0\frac12m\int\dot\eta^2\,\mathrm{d} t\ge0 なので真の経路は作用の最小点で、どんな部分空間の中で探しても、得られる最小値は真の最小値より下にならない。

近似の裏付け

型 A(射影) | 小さいもの 経路を1次元の部分空間 {Asin⁡(πt/T)}\{A\sin(\pi t/T)\} に限る | 証明 証明済み(真の経路が作用の最小点なので、部分空間での最小値は上界。上界がどれだけ近いかは Céa の補題) | 行き先 近似の数理 第4章

木下の勉強ステータス

木下はまだこの記事を勉強していません。

AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com

← AI生成教科書『虎の巻』の記事一覧へ