tecktom.com

出てくるもの(3):電磁場のエネルギー運動量テンソル

2026年10月15日13分で読めますAI生成教科書『虎の巻』第2部 第7章 場の Noether の定理(電荷保存・エネルギー運動量テンソル)
目次
  1. 正準テンソルとその不具合
  2. 改良項を足す
  3. 源があるとき:物質と場を合わせた保存
  4. 木下の層

AI 生成・木下が未勉強

解析力学から様々な理論分野を組み立てる 第2部 第7章 §2.7.5 ★★★ 正準テンソルを対称に直すと、エネルギー密度・Poynting ベクトル・Maxwell の応力が1つのテンソルの成分になる

正準テンソルとその不具合

電磁場の Lagrangian 密度 L=−14μ0FμνFμν\mathcal L=-\frac{1}{4\mu_0}F_{\mu\nu}F^{\mu\nu}(源のない場合。第2部第6章)で、場は4本の AαA_\alpha である。FμνFμνF_{\mu\nu}F^{\mu\nu} を ∂μAα\partial_\mu A_\alpha で微分すると 4Fμα4F^{\mu\alpha} なので ∂L/∂(∂μAα)=−Fμα/μ0\partial\mathcal L/\partial(\partial_\mu A_\alpha)=-F^{\mu\alpha}/\mu_0 で、正準テンソル (2.7.16) は

Tcanμν=−1μ0Fμα∂νAα+14μ0δμνFαβFαβ(2.7.22)T_{\mathrm{can}}{}^\mu{}_\nu=-\frac1{\mu_0}F^{\mu\alpha}\partial_\nu A_\alpha+\frac{1}{4\mu_0}\delta^\mu{}_\nu F_{\alpha\beta}F^{\alpha\beta}\tag{2.7.22}

になる。この式には2つの不具合がある。1つめは AαA_\alpha そのもの(FF でなく)を含むので、ゲージ変換 Aα→Aα−∂αΛA_\alpha\to A_\alpha-\partial_\alpha\Lambda で値が変わること、2つめは μ,ν\mu,\nu について対称でないことである。エネルギー密度は測る量なので、ゲージの選び方によってはいけない。

改良項を足す

例題 2.7.2 の自由度を使う。Bμαν=FμαAν/μ0B^{\mu\alpha}{}_\nu=F^{\mu\alpha}A_\nu/\mu_0 は μ,α\mu,\alpha について反対称なので、その発散を足しても保存則と全体の4元運動量は変わらない。

1μ0∂α(FμαAν)=1μ0(∂αFμα)Aν+1μ0Fμα∂αAν\frac1{\mu_0}\partial_\alpha\bigl(F^{\mu\alpha}A_\nu\bigr)=\frac1{\mu_0}(\partial_\alpha F^{\mu\alpha})A_\nu+\frac1{\mu_0}F^{\mu\alpha}\partial_\alpha A_\nu

で、源のない Maxwell 方程式 ∂αFαμ=0\partial_\alpha F^{\alpha\mu}=0 から第1項は 00 である。これを式 (2.7.22) に足すと ∂νAα−∂αAν=Fνα\partial_\nu A_\alpha-\partial_\alpha A_\nu=F_{\nu\alpha} がまとまり、添字を上げて

Tμν=1μ0(FμαFαν+14ημνFαβFαβ)(2.7.23)T^{\mu\nu}=\frac{1}{\mu_0}\Bigl(F^{\mu\alpha}F_\alpha{}^\nu+\frac14\eta^{\mu\nu}F_{\alpha\beta}F^{\alpha\beta}\Bigr)\tag{2.7.23}

を得る。FF だけで書けているのでゲージ不変で、μ,ν\mu,\nu について対称である。この改良の手続きを Belinfante の方法と呼ぶ。トレースは Tμμ=1μ0(FμαFαμ+F2)=0T^\mu{}_\mu=\frac1{\mu_0}(F^{\mu\alpha}F_{\alpha\mu}+F^2)=0 で、これは光子に質量がないことと対応する。

成分を F0i=Ei/cF_{0i}=E_i/c、F12=−B3F_{12}=-B_3(記法の正本)で書き下すと

T00=ε02E2+B22μ0,T0i=Ti0=Sic,S=E×Bμ0,Tij=−σij,(2.7.24)T^{00}=\frac{\varepsilon_0}{2}\bm{E}^2+\frac{\bm{B}^2}{2\mu_0},\qquad T^{0i}=T^{i0}=\frac{S_i}{c},\quad \bm{S}=\frac{\bm{E}\times\bm{B}}{\mu_0},\qquad T^{ij}=-\sigma_{ij},\tag{2.7.24}

σij=ε0(EiEj−12δijE2)+1μ0(BiBj−12δijB2)\sigma_{ij}=\varepsilon_0\Bigl(E_iE_j-\frac12\delta_{ij}\bm{E}^2\Bigr)+\frac1{\mu_0}\Bigl(B_iB_j-\frac12\delta_{ij}\bm{B}^2\Bigr)

となる。σij\sigma_{ij} は Maxwell の応力テンソルで、負号は計量の符号の流儀 (+,−,−,−)(+,-,-,-) による。エネルギー密度・Poynting ベクトル・Maxwell の応力は、1つのテンソルの別々の成分だった。運動量密度は g=S/c2=ε0E×B\bm{g}=\bm{S}/c^2=\varepsilon_0\bm{E}\times\bm{B} である。

源があるとき:物質と場を合わせた保存

電流 JμJ^\mu があると、電磁場だけの TμνT^{\mu\nu} は保存しない。Maxwell 方程式 ∂μFμα=μ0Jα\partial_\mu F^{\mu\alpha}=\mu_0J^\alpha と Bianchi 恒等式 ∂λFμν+∂μFνλ+∂νFλμ=0\partial_\lambda F_{\mu\nu}+\partial_\mu F_{\nu\lambda}+\partial_\nu F_{\lambda\mu}=0 を使って式 (2.7.23) の発散を計算すると、

∂μTμν=1μ0[(∂μFμα)Fαν+Fμα∂μFαν+12Fαβ∂νFαβ].\partial_\mu T^\mu{}_\nu=\frac1{\mu_0}\Bigl[(\partial_\mu F^{\mu\alpha})F_{\alpha\nu}+F^{\mu\alpha}\partial_\mu F_{\alpha\nu}+\frac12F^{\alpha\beta}\partial_\nu F_{\alpha\beta}\Bigr].

第2項は μ,α\mu,\alpha の反対称性で 12Fμα(∂μFαν+∂αFνμ)\frac12F^{\mu\alpha}(\partial_\mu F_{\alpha\nu}+\partial_\alpha F_{\nu\mu}) と書け、Bianchi 恒等式でこれは −12Fμα∂νFμα-\frac12F^{\mu\alpha}\partial_\nu F_{\mu\alpha} になって第3項と打ち消し合う。残る第1項が

∂μTμν=−FνλJλ(2.7.25)\partial_\mu T^{\mu\nu}=-F^{\nu\lambda}J_\lambda\tag{2.7.25}

である。右辺を成分で書くと、ν=0\nu=0 では −E⋅J/c-\bm{E}\cdot\bm{J}/c、ν=i\nu=i では −(ρE+J×B)i-(\rho\bm{E}+\bm{J}\times\bm{B})_i で、

∂u∂t+∇⋅S=−J⋅E,∂gi∂t−∑j∂jσij=−(ρE+J×B)i(2.7.26)\frac{\partial u}{\partial t}+\nabla\cdot\bm{S}=-\bm{J}\cdot\bm{E},\qquad \frac{\partial g_i}{\partial t}-\sum_j\partial_j\sigma_{ij}=-(\rho\bm{E}+\bm{J}\times\bm{B})_i\tag{2.7.26}

になる。第1式は Poynting の定理で、場のエネルギーは電流にした仕事 J⋅E\bm{J}\cdot\bm{E} の分だけ減る。第2式の右辺は Lorentz 力の密度 f\bm{f} で、物質の運動量密度の増える速さに等しい。したがって g\bm{g} と物質の運動量密度の和は、応力 σij\sigma_{ij} の流れの分しか変わらない。場と物質を合わせた全運動量は保存する。

この保存が成り立つためには、場の運動量 ε0E×B\varepsilon_0\bm{E}\times\bm{B} を数に入れなければならない。動く2つの電荷の間の磁気力は作用反作用の法則を満たさないことがあり(例題 2.7.6)、そのとき粒子だけの運動量は保存しない。第1部第10章で全運動量の保存を並進の対称性から出したとき、系は粒子だけだった。電磁気では並進の対称性を持つのは粒子と場を合わせた作用なので、保存するのも粒子と場を合わせた運動量である。

例題 2.7.5 平面波の放射圧 ★★

オリジナルxx の向きに進む平面波 E=E0cos⁡(kx−ωt)ey\bm{E}=E_0\cos(kx-\omega t)\bm{e}_y、B=(E0/c)cos⁡(kx−ωt)ez\bm{B}=(E_0/c)\cos(kx-\omega t)\bm{e}_z を考える。

  1. (1) T00=uT^{00}=u, SxS_x, TxxT^{xx} を求め、Sx=cuS_x=cu、Txx=uT^{xx}=u を示せ。
  2. (2) この波が x=0x=0 の面に置いた板に垂直に当たり、すべて吸収されるとき、板が単位面積あたり受ける力を時間平均の強度 I=⟨Sx⟩I=\langle S_x \rangle で書け。完全に反射されるときはどうか。
  3. (3) 地球の軌道での太陽光の強度 I=1360.8 W m−2I=1360.8\,\mathrm{W\,m^{-2}} 1G. Kopp and J. L. Lean (2011), Geophys. Res. Lett. 38, L01706(全太陽放射照度 1360.8±0.5 W m−21360.8\pm0.5\,\mathrm{W\,m^{-2}})。下の脚注へで、吸収するときの圧力と、1 m21\,\mathrm{m^2} の板が1日(86400 s86400\,\mathrm s)に受け取る運動量を求めよ。
解答

(1) c2=1/(ε0μ0)c^2=1/(\varepsilon_0\mu_0) から B22μ0=E02cos⁡22μ0c2=ε0E02cos⁡22\frac{B^2}{2\mu_0}=\frac{E_0^2\cos^2}{2\mu_0c^2}=\frac{\varepsilon_0E_0^2\cos^2}{2} で、電場と磁場のエネルギーは等しく u=ε0E02cos⁡2(kx−ωt)u=\varepsilon_0E_0^2\cos^2(kx-\omega t)。E×B=EyBzex\bm{E}\times\bm{B}=E_yB_z\bm{e}_x から Sx=E02cos⁡2/(μ0c)=c ε0E02cos⁡2=cuS_x=E_0^2\cos^2/(\mu_0c)=c\,\varepsilon_0E_0^2\cos^2=cu。応力は Ex=Bx=0E_x=B_x=0 なので σxx=−ε02Ey2−12μ0Bz2=−u\sigma_{xx}=-\frac{\varepsilon_0}{2}E_y^2-\frac{1}{2\mu_0}B_z^2=-u で、Txx=−σxx=uT^{xx}=-\sigma_{xx}=u。

(2) TxxT^{xx} は「xx 成分の運動量が xx の向きに流れる密度」なので、板に毎秒・単位面積あたり流れ込む運動量は ⟨Txx⟩=⟨u⟩=I/c\langle T^{xx} \rangle =\langle u \rangle =I/c である。吸収されるとこの運動量がすべて板に移り、圧力は I/cI/c になる。反射されると、入ってきた運動量 I/cI/c に加えて逆向きに出ていく波の運動量 I/cI/c の分も板が受け持つので、圧力は 2I/c2I/c になる。

(3) I/c=1360.8/(2.998×108)=4.54×10−6 PaI/c=1360.8/(2.998\times10^8)=4.54\times10^{-6}\,\mathrm{Pa}(反射なら 9.08×10−6 Pa9.08\times10^{-6}\,\mathrm{Pa})。1 m21\,\mathrm{m^2} の板が受ける力は 4.54×10−6 N4.54\times10^{-6}\,\mathrm N で、1日に受け取る運動量は 4.54×10−6×86400=0.392 kg m s−14.54\times10^{-6}\times86400=0.392\,\mathrm{kg\,m\,s^{-1}} である。

例題 2.7.6 作用反作用を満たさない磁気力 ★★

オリジナル電荷 q1q_1 が原点を +x+x の向きに速さ vv で、電荷 q2q_2 が点 (0,d,0)(0,d,0) を +y+y の向きに速さ uu で通過する瞬間を考える。v,u≪cv,u\ll c とし、点電荷 qq が速度 v\bm{v} で動くときにつくる磁場を、最低次の B=μ04πq v×err2\bm{B}=\frac{\mu_0}{4\pi}\frac{q\,\bm{v}\times\bm{e}_r}{r^2}(er\bm{e}_r は電荷から測る点への単位ベクトル)で近似する。

  1. (1) q2q_2 が受ける磁気力と、q1q_1 が受ける磁気力を求めよ。
  2. (2) q1=q2=1.0 μCq_1=q_2=1.0\,\mu\mathrm C, d=0.10 md=0.10\,\mathrm m, v=u=1.0×103 m/sv=u=1.0\times10^3\,\mathrm{m/s} で、磁気力と Coulomb 力の大きさを比べよ。
  3. (3) 粒子の運動量の和はこの瞬間にどう変わるか。保存則とどう折り合うかを述べよ。
解答

(1) q1q_1 から q2q_2 への単位ベクトルは ey\bm{e}_y で、q2q_2 の位置での磁場は B1=μ0q1v4πd2ex×ey=μ0q1v4πd2ez\bm{B}_1=\frac{\mu_0q_1v}{4\pi d^2}\bm{e}_x\times\bm{e}_y=\frac{\mu_0q_1v}{4\pi d^2}\bm{e}_z。q2q_2 が受ける力は q2u ey×B1=μ0q1q2vu4πd2exq_2u\,\bm{e}_y\times\bm{B}_1=\frac{\mu_0q_1q_2vu}{4\pi d^2}\bm{e}_x である。q2q_2 から q1q_1 への単位ベクトルは −ey-\bm{e}_y で、q2q_2 の速度 u eyu\,\bm{e}_y と平行なので B2=0\bm{B}_2=\bm{0}、q1q_1 が受ける磁気力は 00 である(図 2.7.3)。

(2) 磁気力は μ04πq1q2vud2=10−7×10−12×10610−2=1.0×10−11 N\frac{\mu_0}{4\pi}\frac{q_1q_2vu}{d^2}=10^{-7}\times\frac{10^{-12}\times10^6}{10^{-2}}=1.0\times10^{-11}\,\mathrm N、Coulomb 力は 14πε0q1q2d2=8.99×109×10−12/10−2=0.899 N\frac{1}{4\pi\varepsilon_0}\frac{q_1q_2}{d^2}=8.99\times10^9\times10^{-12}/10^{-2}=0.899\,\mathrm N。比は vu/c2=1.1×10−11vu/c^2=1.1\times10^{-11} で、磁気力は電気力に比べて (v/c)2(v/c)^2 の程度に小さい。

(3) Coulomb 力は(この近似で)大きさが等しく向きが逆なので、粒子の運動量の和の変化は磁気力だけで決まり、μ0q1q2vu4πd2ex\frac{\mu_0q_1q_2vu}{4\pi d^2}\bm{e}_x の速さで増える。粒子だけの運動量は保存しない。式 (2.7.26) によれば、その分だけ場の運動量 ∫ε0E×B d3x\int\varepsilon_0\bm{E}\times\bm{B}\,\mathrm{d}^3x が減り、粒子と場を合わせた運動量は保存する。

図

図 2.7.3 例題 2.7.6 の配置。q1q_1 の磁場は q2q_2 の位置で紙面の手前向き(⊙\odot)で、q2q_2 は +x+x の向きの磁気力を受ける。q1q_1 は q2q_2 の速度の延長線上にあるので q2q_2 の磁場を受けず、磁気力は片方にしか働かない

近似の裏付け

型 B(極限) | 小さいもの v/c→0v/c\to0(電荷の速さが光速に比べて小さい。磁場は Biot–Savart の最低次) | 証明 要確認 | 行き先 近似の数理 第1章

ほかの巻では

1 G. Kopp and J. L. Lean (2011), Geophys. Res. Lett. 38, L01706(全太陽放射照度 1360.8±0.5 W m−21360.8\pm0.5\,\mathrm{W\,m^{-2}})。

木下の勉強ステータス

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

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

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