解析力学から様々な理論分野を組み立てる 第2部 第5章 §2.5.1 ★★★ 場の Lagrangian は空間積分 L = ∫ L d 3 x L=\int\mathcal L\,\mathrm{d}^3x L = ∫ L d 3 x で、中身 L \mathcal L L は各点の値と1階微分だけで決まる
前の章(第2部第4章 )では、ばねでつないだ原子の列の Lagrangian から出発し、格子間隔 a a a を 0 0 0 に送る極限で、原子の変位 q n ( t ) q_n(t) q n ( t ) の並びが連続な関数 ϕ ( t , x ) \phi(t,x) ϕ ( t , x ) に置き換わり、運動方程式が波動方程式になることを見た。そのとき Lagrangian は原子についての和から x x x についての積分 L = ∫ L d x L=\int\mathcal L\,\mathrm{d} x L = ∫ L d x に変わり、積分の中身 L \mathcal L L は各点の ϕ \phi ϕ の時間微分と空間微分だけで書けていた。この章では、この積分の中身を主役にする。
この章で選ぶものはLagrangian 密度 L ( ϕ , ∂ t ϕ , ∇ ϕ , t , x ) \mathcal L(\phi,\partial_t\phi,\nabla\phi,t,\bm{x}) L ( ϕ , ∂ t ϕ , ∇ ϕ , t , x ) である。作用は L \mathcal L L を時間と空間の両方について積分したもの
S [ ϕ ] = ∫ t 1 t 2 d t ∫ d 3 x L S[\phi]=\int_{t_1}^{t_2}\mathrm{d} t\int\mathrm{d}^3x\;\mathcal L S [ ϕ ] = ∫ t 1 t 2 d t ∫ d 3 x L
になる。当てる道具は第1部の汎関数と変分(第1部第3章 )と Euler–Lagrange 方程式の導き方(第1部第4章 )で、変えるのは1か所だけである。粒子では時刻 t t t の1変数で部分積分したが、場では時刻と空間の4変数で部分積分する。そのための道具が発散定理で、これを使うと場の Euler–Lagrange 方程式
∂ μ ∂ L ∂ ( ∂ μ ϕ ) − ∂ L ∂ ϕ = 0 \partial_\mu\frac{\partial \mathcal L}{\partial (\partial_\mu\phi)}-\frac{\partial \mathcal L}{\partial \phi}=0 ∂ μ ∂ ( ∂ μ ϕ ) ∂ L − ∂ ϕ ∂ L = 0
が出てくる。この1本の式から、弦の波動方程式・Klein–Gordon 方程式・Schrödinger 方程式が、L \mathcal L L の選び方の違いとして出てくる。後半では、部分積分で出る境界の項から弦の端の条件(固定端・自由端・ばねつきの端)を導き、粒子の Legendre 変換を空間の各点で行って Hamiltonian 密度とエネルギーの流れを作る。電磁場の Lagrangian は次の章(第2部第6章 )、場の Noether の定理とエネルギー運動量テンソルはその次の章(第2部第7章 )で扱う。
原子の列から弦へ(前の章の結果)
前の章(第2部第4章 )の出発点は、質量 m m m の原子が間隔 a a a で1列に並び、隣どうしがばね定数 K K K のばねでつながった系だった。線密度 ρ = m / a \rho=m/a ρ = m / a と1次元の弾性率 Y = K a Y=Ka Y = K a を固定して格子間隔 a a a を 0 0 0 に送ると、列の Lagrangian は第4章の式 (2.4.14) の形になる。この章では変位を ϕ ( t , x ) \phi(t,x) ϕ ( t , x ) と書いて
L = ∫ d x [ ρ 2 ( ∂ t ϕ ) 2 − Y 2 ( ∂ x ϕ ) 2 ] , ρ = m a , Y = K a (2.5.1) L=\int\mathrm{d} x\;\Bigl[\frac\rho2(\partial_t\phi)^2-\frac Y2(\partial_x\phi)^2\Bigr],\qquad \rho=\frac ma,\quad Y=Ka\tag{2.5.1} L = ∫ d x [ 2 ρ ( ∂ t ϕ ) 2 − 2 Y ( ∂ x ϕ ) 2 ] , ρ = a m , Y = K a ( 2.5.1 )
と書く。ρ \rho ρ は単位長さあたりの質量(線密度)、Y Y Y は1次元の弾性率である。棒の縦振動なら Y Y Y は Young 率と断面積の積、張った弦の横振動なら張力に当たる(2.5.6 節)。前の章では、この L L L から波動方程式 ρ ∂ t 2 ϕ = Y ∂ x 2 ϕ \rho\,\partial_t^2\phi=Y\partial_x^2\phi ρ ∂ t 2 ϕ = Y ∂ x 2 ϕ と波の速さ v = Y / ρ v=\sqrt{Y/\rho} v = Y / ρ を出した。
Lagrangian 密度
式 (2.5.1) の被積分関数を取り出して L \mathcal L L と書く。
定義 2.5.1(Lagrangian 密度)
空間の各点 x \bm{x} x に実数の値 ϕ ( t , x ) \phi(t,\bm{x}) ϕ ( t , x ) を持つ量を場 と呼ぶ。場の Lagrangian が
L ( t ) = ∫ d 3 x L ( ϕ ( t , x ) , ∂ t ϕ ( t , x ) , ∇ ϕ ( t , x ) , t , x ) (2.5.2) L(t)=\int\mathrm{d}^3x\;\mathcal L\bigl(\phi(t,\bm{x}),\,\partial_t\phi(t,\bm{x}),\,\nabla\phi(t,\bm{x}),\,t,\,\bm{x}\bigr)\tag{2.5.2} L ( t ) = ∫ d 3 x L ( ϕ ( t , x ) , ∂ t ϕ ( t , x ) , ∇ ϕ ( t , x ) , t , x ) ( 2.5.2 )
の形に書けるとき、L \mathcal L L をLagrangian 密度 と呼ぶ。作用は
S [ ϕ ] = ∫ t 1 t 2 L d t = ∫ t 1 t 2 d t ∫ d 3 x L (2.5.3) S[\phi]=\int_{t_1}^{t_2}L\,\mathrm{d} t=\int_{t_1}^{t_2}\mathrm{d} t\int\mathrm{d}^3x\;\mathcal L\tag{2.5.3} S [ ϕ ] = ∫ t 1 t 2 L d t = ∫ t 1 t 2 d t ∫ d 3 x L ( 2.5.3 )
である。L \mathcal L L の次元は「エネルギー÷体積」(1次元の場ではエネルギー÷長さ)になる。
粒子の力学との対応はこうなる。粒子では一般化座標 q i q_i q i に番号 i = 1 , … , f i=1,\dots,f i = 1 , … , f が付いていた。場では番号の代わりに空間の点 x \bm{x} x が付き、ϕ ( t , x ) \phi(t,\bm{x}) ϕ ( t , x ) は「番号 x \bm{x} x の一般化座標の時刻 t t t での値」である。自由度は連続無限個になり、和 ∑ i \sum_i ∑ i は積分 ∫ d 3 x \int\mathrm{d}^3x ∫ d 3 x に変わる。同じ記号 x \bm{x} x が出てくるが、x \bm{x} x は力学変数ではなく番号であり、時刻 t t t と同じく「どこで値を読むか」を指定する引数である。
Lagrangian が密度の積分の形 (2.5.2) に書けることは、どの場でも成り立つわけではない。原子の列では、各原子が隣の原子とだけばねでつながっていたので、L L L は「点 n n n の付近の量」の和になっていた。隣どうしの差 q n + 1 − q n q_{n+1}-q_n q n + 1 − q n は連続の極限で微分 ∂ x ϕ \partial_x\phi ∂ x ϕ になる。こうして、近くどうしだけが相互作用する系(局所的な系)では、L \mathcal L L が各点の ϕ \phi ϕ とその1階微分で書くことができる。遠く離れた点どうしが直接相互作用する系(たとえば重力や Coulomb 力で直接結ばれた粒子の集まり)は、このままでは密度の形にならない。そうした遠隔の相互作用も、間に場(電磁場)を立ててその場を局所的な L \mathcal L L で書けば、密度の形に収まる。次の章で電磁場の L \mathcal L L を書くのはこのためである。
この巻では L \mathcal L L が ϕ \phi ϕ の2階以上の微分を含まない場合だけを扱う。粒子の L ( q , q ˙ , t ) L(q,\dot q,t) L ( q , q ˙ , t ) が加速度を含まなかったのと同じ制限で、物理で使う場の理論(波動・電磁場・Klein–Gordon・Dirac)はすべてこの範囲に入る。
相対論の記法
時刻と空間をまとめて x μ = ( c t , x ) x^\mu=(ct,\bm{x}) x μ = ( c t , x ) と書くと、∂ μ = ( 1 c ∂ t , ∇ ) \partial_\mu=(\frac1c\partial_t,\nabla) ∂ μ = ( c 1 ∂ t , ∇ ) 、体積要素は d 4 x = c d t d 3 x \mathrm{d}^4x=c\,\mathrm{d} t\,\mathrm{d}^3x d 4 x = c d t d 3 x だから、作用は
S [ ϕ ] = 1 c ∫ L ( ϕ , ∂ μ ϕ , x ) d 4 x (2.5.4) S[\phi]=\frac1c\int\mathcal L\bigl(\phi,\partial_\mu\phi,x\bigr)\,\mathrm{d}^4x\tag{2.5.4} S [ ϕ ] = c 1 ∫ L ( ϕ , ∂ μ ϕ , x ) d 4 x ( 2.5.4 )
と書ける(記法は 第2部第3章 と同じで、計量は η μ ν = d i a g ( + 1 , − 1 , − 1 , − 1 ) \eta_{\mu\nu}=\mathrm{diag}(+1,-1,-1,-1) η μν = diag ( + 1 , − 1 , − 1 , − 1 ) )。L \mathcal L L が ∂ μ ϕ \partial_\mu\phi ∂ μ ϕ を通してだけ微分を含むことを、引数を ∂ μ ϕ \partial_\mu\phi ∂ μ ϕ とまとめて書いて表す。前の係数 1 / c 1/c 1/ c は、L \mathcal L L をエネルギー密度の次元に保つための約束である。弦のように相対論と関係しない場では、d t d 3 x \mathrm{d} t\,\mathrm{d}^3x d t d 3 x (1次元なら d t d x \mathrm{d} t\,\mathrm{d} x d t d x )の形 (2.5.3) をそのまま使う。
図 2.5.1 原子の列の分散関係(点線)と、その連続体の極限である弦の分散関係(実線)。k a → 0 ka\to0 k a → 0 で2本は接し、ずれの割合は ( k a ) 2 / 24 (ka)^2/24 ( k a ) 2 /24 から始まる。k a = π / 4 ka=\pi/4 k a = π /4 でも差は 2.6 % にとどまる
式 (2.5.1) は格子間隔を 0 0 0 に送った極限なので、格子の運動方程式と比べてどれだけずれるかは分散関係で測ることができる。第4章の式 (2.4.16) で、格子の分散関係 ω = 2 K / m sin ( k a / 2 ) \omega=2\sqrt{K/m}\,\sin(ka/2) ω = 2 K / m sin ( k a /2 ) を k a ka k a について展開して
ω = v k ( 1 − ( k a ) 2 24 + O ( ( k a ) 4 ) ) , v = a K / m = Y / ρ (2.5.5) \omega=vk\Bigl(1-\frac{(ka)^2}{24}+\mathcal{O}\bigl((ka)^4\bigr)\Bigr),\qquad v=a\sqrt{K/m}=\sqrt{Y/\rho}\tag{2.5.5} ω = v k ( 1 − 24 ( k a ) 2 + O ( ( k a ) 4 ) ) , v = a K / m = Y / ρ ( 2.5.5 )
を得ている(図 2.5.1 )。連続体の L \mathcal L L は、k a ka k a の小さい(波長が格子間隔より十分長い)成分について、相対誤差 ( k a ) 2 / 24 (ka)^2/24 ( k a ) 2 /24 で格子を再現する。
近似の裏付け
型 B(極限) | 小さいもの k a → 0 ka\to0 k a → 0 (波長が格子間隔 a a a より十分長い) | 証明 条件付き(波長の長い成分だけ。分散関係のずれ ( k a ) 2 / 24 (ka)^2/24 ( k a ) 2 /24 は式 (2.5.5) (第4章の式 (2.4.16) )で評価できる。解そのものの収束の評価は要確認) | 行き先 近似の数理 第9章
例題 2.5.1 原子の列から作った弦の上を進む波 ★★
オリジナル 質量 m = 1.0 × 10 − 3 k g m=1.0\times10^{-3}\,\mathrm{kg} m = 1.0 × 1 0 − 3 kg の原子が間隔 a = 1.0 × 10 − 2 m a=1.0\times10^{-2}\,\mathrm m a = 1.0 × 1 0 − 2 m で並び、ばね定数 K = 1.0 × 10 3 N m − 1 K=1.0\times10^{3}\,\mathrm{N\,m^{-1}} K = 1.0 × 1 0 3 N m − 1 のばねでつながっている。
(1) 連続体の線密度 ρ \rho ρ 、弾性率 Y Y Y 、波の速さ v v v を求めよ。
(2) 進む波 ϕ = A cos ( k x − ω t ) \phi=A\cos(kx-\omega t) ϕ = A cos ( k x − ω t ) (ω = v k \omega=vk ω = v k )について、Lagrangian 密度 L = ρ 2 ( ∂ t ϕ ) 2 − Y 2 ( ∂ x ϕ ) 2 \mathcal L=\frac\rho2(\partial_t\phi)^2-\frac Y2(\partial_x\phi)^2 L = 2 ρ ( ∂ t ϕ ) 2 − 2 Y ( ∂ x ϕ ) 2 が各点で 0 0 0 になることを示せ。
(3) 同じ波のエネルギー密度 E = ρ 2 ( ∂ t ϕ ) 2 + Y 2 ( ∂ x ϕ ) 2 \mathcal E=\frac\rho2(\partial_t\phi)^2+\frac Y2(\partial_x\phi)^2 E = 2 ρ ( ∂ t ϕ ) 2 + 2 Y ( ∂ x ϕ ) 2 を求めよ。A = 1.0 × 10 − 3 m A=1.0\times10^{-3}\,\mathrm m A = 1.0 × 1 0 − 3 m 、波長 0.50 m 0.50\,\mathrm m 0.50 m のとき、時間平均の値はいくらか。
(4) 定在波 ϕ = A sin k x cos ω t \phi=A\sin kx\cos\omega t ϕ = A sin k x cos ω t では L \mathcal L L は各点では 0 0 0 にならないが、1波長と1周期にわたる平均は 0 0 0 になることを示せ。
解答 (1) ρ = m / a = 0.10 k g m − 1 \rho=m/a=0.10\,\mathrm{kg\,m^{-1}} ρ = m / a = 0.10 kg m − 1 、Y = K a = 10 N Y=Ka=10\,\mathrm N Y = K a = 10 N 、v = Y / ρ = 10 m s − 1 v=\sqrt{Y/\rho}=10\,\mathrm{m\,s^{-1}} v = Y / ρ = 10 m s − 1 。
(2) ∂ t ϕ = A ω sin ( k x − ω t ) \partial_t\phi=A\omega\sin(kx-\omega t) ∂ t ϕ = A ω sin ( k x − ω t ) 、∂ x ϕ = − A k sin ( k x − ω t ) \partial_x\phi=-Ak\sin(kx-\omega t) ∂ x ϕ = − A k sin ( k x − ω t ) だから
L = A 2 2 ( ρ ω 2 − Y k 2 ) sin 2 ( k x − ω t ) . \mathcal L=\frac{A^2}{2}\bigl(\rho\omega^2-Yk^2\bigr)\sin^2(kx-\omega t). L = 2 A 2 ( ρ ω 2 − Y k 2 ) sin 2 ( k x − ω t ) .
ω = v k \omega=vk ω = v k なら ρ ω 2 = ρ v 2 k 2 = Y k 2 \rho\omega^2=\rho v^2k^2=Yk^2 ρ ω 2 = ρ v 2 k 2 = Y k 2 なので括弧が 0 0 0 になり、L = 0 \mathcal L=0 L = 0 が各点・各時刻で成り立つ。進む波では、運動エネルギーの密度と位置エネルギーの密度が各点でつねに等しい。
(3) 同じ計算で E = A 2 2 ( ρ ω 2 + Y k 2 ) sin 2 ( k x − ω t ) = ρ ω 2 A 2 sin 2 ( k x − ω t ) \mathcal E=\frac{A^2}{2}(\rho\omega^2+Yk^2)\sin^2(kx-\omega t)=\rho\omega^2A^2\sin^2(kx-\omega t) E = 2 A 2 ( ρ ω 2 + Y k 2 ) sin 2 ( k x − ω t ) = ρ ω 2 A 2 sin 2 ( k x − ω t ) 。sin 2 \sin^2 sin 2 の時間平均は 1 / 2 1/2 1/2 だから、平均は 1 2 ρ ω 2 A 2 \frac12\rho\omega^2A^2 2 1 ρ ω 2 A 2 である。k = 2 π / 0.50 = 12.6 m − 1 k=2\pi/0.50=12.6\,\mathrm m^{-1} k = 2 π /0.50 = 12.6 m − 1 、ω = v k = 126 s − 1 \omega=vk=126\,\mathrm s^{-1} ω = v k = 126 s − 1 を入れると
E ‾ = 1 2 × 0.10 × ( 126 ) 2 × ( 1.0 × 10 − 3 ) 2 = 7.90 × 10 − 4 J m − 1 . \overline{\mathcal E}=\tfrac12\times0.10\times(126)^2\times(1.0\times10^{-3})^2=7.90\times10^{-4}\,\mathrm{J\,m^{-1}}. E = 2 1 × 0.10 × ( 126 ) 2 × ( 1.0 × 1 0 − 3 ) 2 = 7.90 × 1 0 − 4 J m − 1 .
(4) ∂ t ϕ = − A ω sin k x sin ω t \partial_t\phi=-A\omega\sin kx\sin\omega t ∂ t ϕ = − A ω sin k x sin ω t 、∂ x ϕ = A k cos k x cos ω t \partial_x\phi=Ak\cos kx\cos\omega t ∂ x ϕ = A k cos k x cos ω t 、Y k 2 = ρ ω 2 Yk^2=\rho\omega^2 Y k 2 = ρ ω 2 を使うと
L = ρ ω 2 A 2 2 ( sin 2 k x sin 2 ω t − cos 2 k x cos 2 ω t ) . \mathcal L=\frac{\rho\omega^2A^2}{2}\bigl(\sin^2kx\,\sin^2\omega t-\cos^2kx\,\cos^2\omega t\bigr). L = 2 ρ ω 2 A 2 ( sin 2 k x sin 2 ω t − cos 2 k x cos 2 ω t ) .
これは節(sin k x = 0 \sin kx=0 sin k x = 0 )の上では − 1 2 ρ ω 2 A 2 cos 2 ω t ≠ 0 -\frac12\rho\omega^2A^2\cos^2\omega t\neq0 − 2 1 ρ ω 2 A 2 cos 2 ω t = 0 である。1波長の平均で sin 2 k x \sin^2kx sin 2 k x と cos 2 k x \cos^2kx cos 2 k x はともに 1 / 2 1/2 1/2 、1周期の平均で sin 2 ω t \sin^2\omega t sin 2 ω t と cos 2 ω t \cos^2\omega t cos 2 ω t もともに 1 / 2 1/2 1/2 なので、平均は ρ ω 2 A 2 2 ( 1 4 − 1 4 ) = 0 \frac{\rho\omega^2A^2}{2}\bigl(\frac14-\frac14\bigr)=0 2 ρ ω 2 A 2 ( 4 1 − 4 1 ) = 0 になる。運動方程式を満たす波では、L \mathcal L L の平均が 0 0 0 になる(運動と位置のエネルギーが平均として等分される)。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com