tecktom.com

この章のまとめ

2026年10月15日10分で読めますAI生成教科書『虎の巻』第2部 第5章 場の Lagrangian 密度と場の Euler–Lagrange 方程式

AI 生成・木下が未勉強

解析力学から様々な理論分野を組み立てる 第2部 第5章 §2.5.7 ★★★ 粒子の LL を空間の各点に並べると L\mathcal L になり、部分積分を発散定理に替えるだけで第1部の道具がそのまま使える

この章で選んだものは Lagrangian 密度 L(ϕ,∂μϕ,x)\mathcal L(\phi,\partial_\mu\phi,x) で、作用は S[ϕ]=∫dt d3x L=1c∫L d4xS[\phi]=\int\mathrm{d} t\,\mathrm{d}^3x\,\mathcal L=\frac1c\int\mathcal L\,\mathrm{d}^4x である。場は、空間の点 x\bm{x} を番号とする連続無限個の一般化座標とみなすことができる。密度の形に書けるのは、近くどうしだけが相互作用する局所的な系である。

当てた道具は第1部の変分で、粒子と違うのは部分積分を発散定理に替えたことだけだった。その結果、作用の停留は各点ごとの偏微分方程式∂μ∂L∂(∂μϕ)−∂L∂ϕ=0\partial_\mu\frac{\partial \mathcal L}{\partial (\partial_\mu\phi)}-\frac{\partial \mathcal L}{\partial \phi}=0(定理 2.5.2)になり、L\mathcal L の選び方で波動方程式・Klein–Gordon 方程式・Schrödinger 方程式が出てきた。部分積分で出る境界の項からは、空間の端の条件(固定端・自由端・ばねつきの端)が決まる。また、全発散を足しても方程式は変わらない(命題 2.5.3)。各点で Legendre 変換すると、正準運動量の密度 π=∂L∂(∂tϕ)\pi=\frac{\partial \mathcal L}{\partial (\partial_t\phi)} と Hamiltonian 密度 H\mathcal H が作られ、Poisson 括弧 {ϕ(x),π(y)}=δ(3)(x−y)\{\phi(\bm{x}),\pi(\bm{y})\}=\delta^{(3)}(\bm{x}-\bm{y}) と局所的なエネルギーの保存則 ∂tH+∂xS=0\partial_t\mathcal H+\partial_x\mathcal S=0 が出る。

次の章では、ベクトルの場 AμA_\mu について Lorentz 不変な L=−14μ0FμνFμν−JμAμ\mathcal L=-\frac1{4\mu_0}F_{\mu\nu}F^{\mu\nu}-J^\mu A_\mu を選び、この章の定理 2.5.2 を4成分それぞれに当てて Maxwell 方程式を出す(第2部第6章)。そこでは A0A_0 の時間微分が L\mathcal L に現れないので、π\pi の1成分が恒等的に 00 になる。これがゲージの自由度(第2部第8章)の表れである。エネルギーの流れ S\mathcal S は場の Noether の定理(第2部第7章)で一般の形にする。

総合問題 ばねの床に載った原子の列から Klein–Gordon 型の場へ ★★★

オリジナル質量 mm の原子が間隔 aa で並び、隣どうしをばね定数 KK のばねでつなぐ。さらに各原子を、ばね定数 K0K_0 のばねでそれぞれの釣り合いの位置につなぐ(ばねの床)。nn 番目の原子の変位を qn(t)q_n(t) とすると、Lagrangian は

L=∑n[m2q˙n2−K2(qn+1−qn)2−K02qn2]L=\sum_n\Bigl[\frac m2\dot q_n^2-\frac K2(q_{n+1}-q_n)^2-\frac{K_0}{2}q_n^2\Bigr]

である。ρ=m/a\rho=m/a、Y=KaY=Ka、v=Y/ρv=\sqrt{Y/\rho}、ω0=K0/m\omega_0=\sqrt{K_0/m} と置く。

  1. (1) (使う道具:第2部第4章・2.5.1節)qn(t)=ϕ(t,na)q_n(t)=\phi(t,na) と置いて連続体の極限を取り、Lagrangian 密度を ρ,v,ω0\rho,v,\omega_0 で書け。ρ=0.010 kg m−1\rho=0.010\,\mathrm{kg\,m^{-1}}、Y=100 NY=100\,\mathrm N、ω0=300 s−1\omega_0=300\,\mathrm s^{-1}、a=0.050 ma=0.050\,\mathrm m のとき、mm、KK、K0K_0 はいくらか。
  2. (2) (使う道具:2.5.3節)場の Euler–Lagrange 方程式を書け。
  3. (3) (使う道具:2.5.3節の例題 2.5.3・第2部第4章)分散関係を求めよ。原子の列のまま平面波 qn=ei(kna−ωt)q_n=\mathrm{e}^{\mathrm{i}(kna-\omega t)} を入れた分散関係と比べ、(1) の数値で kℓ=π/2, 3π/2, 5π/2k\ell=\pi/2,\ 3\pi/2,\ 5\pi/2(ℓ=1.0 m\ell=1.0\,\mathrm m)のとき、連続体の ω\omega が格子の ω\omega よりどれだけの割合で大きいかを求めよ。
  4. (4) (使う道具:2.5.4節)長さ ℓ=1.0 m\ell=1.0\,\mathrm m の列の端 x=0x=0 を固定し、x=ℓx=\ell を自由端にする。固有振動の角周波数 ωn\omega_n を求め、下から3つの振動数を数値で答えよ。
  5. (5) (使う道具:2.5.5節・第1部第8章)正準運動量の密度と Hamiltonian 密度を求め、エネルギーの流れ S\mathcal S を書け。(4) の端の条件のもとで全エネルギーが保存することを示し、基本モード ϕ=Asin⁡k1xcos⁡ω1t\phi=A\sin k_1x\cos\omega_1t のエネルギーを求めよ。
解答

(1) 運動の項と隣どうしのばねの項は 2.5.1節と同じで、ρ2(∂tϕ)2−Y2(∂xϕ)2\frac\rho2(\partial_t\phi)^2-\frac Y2(\partial_x\phi)^2 になる。床のばねの項は ∑nK02qn2→∫dxaK02ϕ2\sum_n\frac{K_0}{2}q_n^2\to\int\frac{\mathrm{d} x}{a}\frac{K_0}{2}\phi^2 で、K0a=ma⋅K0m=ρ ω02\frac{K_0}{a}=\frac ma\cdot\frac{K_0}{m}=\rho\,\omega_0^2 だから

L=ρ2[(∂tϕ)2−v2(∂xϕ)2−ω02ϕ2].\mathcal L=\frac\rho2\Bigl[(\partial_t\phi)^2-v^2(\partial_x\phi)^2-\omega_0^2\phi^2\Bigr].

数値は m=ρa=5.0×10−4 kgm=\rho a=5.0\times10^{-4}\,\mathrm{kg}、K=Y/a=2.0×103 N m−1K=Y/a=2.0\times10^3\,\mathrm{N\,m^{-1}}、K0=mω02=45 N m−1K_0=m\omega_0^2=45\,\mathrm{N\,m^{-1}}。

(2) ∂L∂(∂tϕ)=ρ ∂tϕ\frac{\partial \mathcal L}{\partial (\partial_t\phi)}=\rho\,\partial_t\phi、∂L∂(∂xϕ)=−ρv2∂xϕ\frac{\partial \mathcal L}{\partial (\partial_x\phi)}=-\rho v^2\partial_x\phi、∂L∂ϕ=−ρω02ϕ\frac{\partial \mathcal L}{\partial \phi}=-\rho\omega_0^2\phi を式 (2.5.12) に入れて

∂t2ϕ−v2∂x2ϕ+ω02ϕ=0.\partial_t^2\phi-v^2\partial_x^2\phi+\omega_0^2\phi=0 .

c→vc\to v、κ→ω0/v\kappa\to\omega_0/v と読み替えると Klein–Gordon 方程式 (2.5.17) と同じ形で、床のばねが「質量」の項の役をする。

(3) 平面波 ei(kx−ωt)\mathrm{e}^{\mathrm{i}(kx-\omega t)} を入れて ω2=ω02+v2k2\omega^2=\omega_0^2+v^2k^2。ω0\omega_0 より低い角周波数の波は伝わらない。原子の列では mq¨n=K(qn+1−2qn+qn−1)−K0qnm\ddot q_n=K(q_{n+1}-2q_n+q_{n-1})-K_0q_n に平面波を入れ、eika−2+e−ika=−4sin⁡2(ka/2)\mathrm{e}^{\mathrm{i} ka}-2+\mathrm{e}^{-\mathrm{i} ka}=-4\sin^2(ka/2) を使って

ωlat2=ω02+4Kmsin⁡2ka2.\omega_{\mathrm{lat}}^2=\omega_0^2+\frac{4K}{m}\sin^2\frac{ka}{2}.

4Kmsin⁡2ka2=v2k2(1−(ka)212+⋯ )\frac{4K}{m}\sin^2\frac{ka}2=v^2k^2\bigl(1-\frac{(ka)^2}{12}+\cdots\bigr) なので、連続体の方がわずかに高い。数値ではkℓ=π/2k\ell=\pi/2(ka=0.079ka=0.079)で 5.5×10−55.5\times10^{-5}、kℓ=3π/2k\ell=3\pi/2(ka=0.24ka=0.24)で 1.6×10−31.6\times10^{-3}、kℓ=5π/2k\ell=5\pi/2(ka=0.39ka=0.39)で 5.6×10−35.6\times10^{-3} の割合だけ大きく、kk が大きいほどずれが大きい(図 2.5.1 の点線と実線の開き方と同じ)。

(4) ϕ=sin⁡kxcos⁡ωt\phi=\sin kx\cos\omega t は x=0x=0 の固定端を満たす。自由端の条件 ∂L∂(∂xϕ)∣x=ℓ=0\frac{\partial \mathcal L}{\partial (\partial_x\phi)}\big|_{x=\ell}=0 から cos⁡kℓ=0\cos k\ell=0、kn=(2n−1)π2ℓk_n=\frac{(2n-1)\pi}{2\ell} で、

ωn=ω02+v2kn2.\omega_n=\sqrt{\omega_0^2+v^2k_n^2}.

v=100 m s−1v=100\,\mathrm{m\,s^{-1}} だから ω1=339 s−1\omega_1=339\,\mathrm s^{-1}、ω2=559 s−1\omega_2=559\,\mathrm s^{-1}、ω3=841 s−1\omega_3=841\,\mathrm s^{-1}、振動数は fn=ωn/(2π)=53.9, 88.9, 134 Hzf_n=\omega_n/(2\pi)=53.9,\ 88.9,\ 134\,\mathrm{Hz}。床のばねが無ければ fn=25, 75, 125 Hzf_n=25,\ 75,\ 125\,\mathrm{Hz} の奇数倍の列だったが、ω0\omega_0 のために低いモードほど持ち上がり、等間隔でなくなる。

(5) π=ρ ∂tϕ\pi=\rho\,\partial_t\phi、

H=π ∂tϕ−L=π22ρ+ρv22(∂xϕ)2+ρω022ϕ2.\mathcal H=\pi\,\partial_t\phi-\mathcal L=\frac{\pi^2}{2\rho}+\frac{\rho v^2}{2}(\partial_x\phi)^2+\frac{\rho\omega_0^2}{2}\phi^2 .

∂tH\partial_t\mathcal H を計算して運動方程式 (2) を使うと ∂tH+∂xS=0\partial_t\mathcal H+\partial_x\mathcal S=0、S=−ρv2∂xϕ ∂tϕ\mathcal S=-\rho v^2\partial_x\phi\,\partial_t\phi となる(床のばねの項は ∂t(ρω022ϕ2)=ρω02ϕ ∂tϕ\partial_t(\frac{\rho\omega_0^2}{2}\phi^2)=\rho\omega_0^2\phi\,\partial_t\phi で、運動方程式の ω02ϕ\omega_0^2\phi の項と打ち消し合うので、流れには入らない)。xx で積分すると dHdt=−[S]0ℓ\frac{\mathrm{d} H}{\mathrm{d} t}=-[\mathcal S]_0^\ell で、x=0x=0 では ∂tϕ=0\partial_t\phi=0、x=ℓx=\ell では ∂xϕ=0\partial_x\phi=0 なので dHdt=0\frac{\mathrm{d} H}{\mathrm{d} t}=0。基本モードでは、∫0ℓsin⁡2k1x dx=∫0ℓcos⁡2k1x dx=ℓ/2\int_0^\ell\sin^2k_1x\,\mathrm{d} x=\int_0^\ell\cos^2k_1x\,\mathrm{d} x=\ell/2 を使って

H=ρA2ℓ4[ω12sin⁡2ω1t+(v2k12+ω02)cos⁡2ω1t]=ρ ω12A2ℓ4H=\frac{\rho A^2\ell}{4}\Bigl[\omega_1^2\sin^2\omega_1t+(v^2k_1^2+\omega_0^2)\cos^2\omega_1t\Bigr]=\frac{\rho\,\omega_1^2A^2\ell}{4}

となり、例題 2.5.5 の ω\omega を ω1\omega_1 に替えた形になる。

木下の勉強ステータス

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

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

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