解析力学から様々な理論分野を組み立てる 第2部 第7章 §2.7.5 ★★★ 正準テンソルを対称に直すと、エネルギー密度・Poynting ベクトル・Maxwell の応力が1つのテンソルの成分になる
正準テンソルとその不具合
電磁場の Lagrangian 密度 L = − 1 4 μ 0 F μ ν F μ ν \mathcal L=-\frac{1}{4\mu_0}F_{\mu\nu}F^{\mu\nu} L = − 4 μ 0 1 F μν F μν (源のない場合。第2部第6章 )で、場は4本の A α A_\alpha A α である。F μ ν F μ ν F_{\mu\nu}F^{\mu\nu} F μν F μν を ∂ μ A α \partial_\mu A_\alpha ∂ μ A α で微分すると 4 F μ α 4F^{\mu\alpha} 4 F μα なので ∂ L / ∂ ( ∂ μ A α ) = − F μ α / μ 0 \partial\mathcal L/\partial(\partial_\mu A_\alpha)=-F^{\mu\alpha}/\mu_0 ∂ L / ∂ ( ∂ μ A α ) = − F μα / μ 0 で、正準テンソル (2.7.16) は
T c a n μ ν = − 1 μ 0 F μ α ∂ ν A α + 1 4 μ 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} T can μ ν = − μ 0 1 F μα ∂ ν A α + 4 μ 0 1 δ μ ν F α β F α β ( 2.7.22 )
になる。この式には2つの不具合がある。1つめは A α A_\alpha A α そのもの(F F F でなく)を含むので、ゲージ変換 A α → A α − ∂ α Λ A_\alpha\to A_\alpha-\partial_\alpha\Lambda A α → A α − ∂ α Λ で値が変わること、2つめは μ , ν \mu,\nu μ , ν について対称でないことである。エネルギー密度は測る量なので、ゲージの選び方によってはいけない。
改良項を足す
例題 2.7.2 の自由度を使う。B μ α ν = F μ α A ν / μ 0 B^{\mu\alpha}{}_\nu=F^{\mu\alpha}A_\nu/\mu_0 B μα ν = F μα A ν / μ 0 は μ , α \mu,\alpha μ , α について反対称なので、その発散を足しても保存則と全体の4元運動量は変わらない。
1 μ 0 ∂ α ( F μ α A ν ) = 1 μ 0 ( ∂ α F μ α ) A ν + 1 μ 0 F μ α ∂ α 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 μ 0 1 ∂ α ( F μα A ν ) = μ 0 1 ( ∂ α F μα ) A ν + μ 0 1 F μα ∂ α A ν
で、源のない Maxwell 方程式 ∂ α F α μ = 0 \partial_\alpha F^{\alpha\mu}=0 ∂ α F α μ = 0 から第1項は 0 0 0 である。これを式 (2.7.22) に足すと ∂ ν A α − ∂ α A ν = F ν α \partial_\nu A_\alpha-\partial_\alpha A_\nu=F_{\nu\alpha} ∂ ν A α − ∂ α A ν = F ν α がまとまり、添字を上げて
T μ ν = 1 μ 0 ( F μ α F α ν + 1 4 η μ ν 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} T μν = μ 0 1 ( F μα F α ν + 4 1 η μν F α β F α β ) ( 2.7.23 )
を得る。F F F だけで書けているのでゲージ不変で、μ , ν \mu,\nu μ , ν について対称である。この改良の手続きを Belinfante の方法と呼ぶ。トレースは T μ μ = 1 μ 0 ( F μ α F α μ + F 2 ) = 0 T^\mu{}_\mu=\frac1{\mu_0}(F^{\mu\alpha}F_{\alpha\mu}+F^2)=0 T μ μ = μ 0 1 ( F μα F α μ + F 2 ) = 0 で、これは光子に質量がないことと対応する。
成分を F 0 i = E i / c F_{0i}=E_i/c F 0 i = E i / c 、F 12 = − B 3 F_{12}=-B_3 F 12 = − B 3 (記法の正本)で書き下すと
T 00 = ε 0 2 E 2 + B 2 2 μ 0 , T 0 i = T i 0 = S i c , S = E × B μ 0 , T i j = − σ i j , (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} T 00 = 2 ε 0 E 2 + 2 μ 0 B 2 , T 0 i = T i 0 = c S i , S = μ 0 E × B , T ij = − σ ij , ( 2.7.24 )
σ i j = ε 0 ( E i E j − 1 2 δ i j E 2 ) + 1 μ 0 ( B i B j − 1 2 δ i j B 2 ) \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 = ε 0 ( E i E j − 2 1 δ ij E 2 ) + μ 0 1 ( B i B j − 2 1 δ ij B 2 )
となる。σ i j \sigma_{ij} σ ij は Maxwell の応力テンソルで、負号は計量の符号の流儀 ( + , − , − , − ) (+,-,-,-) ( + , − , − , − ) による。エネルギー密度・Poynting ベクトル・Maxwell の応力は、1つのテンソルの別々の成分だった。運動量密度は g = S / c 2 = ε 0 E × B \bm{g}=\bm{S}/c^2=\varepsilon_0\bm{E}\times\bm{B} g = S / c 2 = ε 0 E × B である。
源があるとき:物質と場を合わせた保存
電流 J μ J^\mu J μ があると、電磁場だけの T μ ν T^{\mu\nu} T μν は保存しない。Maxwell 方程式 ∂ μ F μ α = μ 0 J α \partial_\mu F^{\mu\alpha}=\mu_0J^\alpha ∂ μ F μα = μ 0 J α と Bianchi 恒等式 ∂ λ F μ ν + ∂ μ F ν λ + ∂ ν F λ μ = 0 \partial_\lambda F_{\mu\nu}+\partial_\mu F_{\nu\lambda}+\partial_\nu F_{\lambda\mu}=0 ∂ λ F μν + ∂ μ F ν λ + ∂ ν F λ μ = 0 を使って式 (2.7.23) の発散を計算すると、
∂ μ T μ ν = 1 μ 0 [ ( ∂ μ F μ α ) F α ν + F μ α ∂ μ F α ν + 1 2 F α β ∂ ν 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]. ∂ μ T μ ν = μ 0 1 [ ( ∂ μ F μα ) F α ν + F μα ∂ μ F α ν + 2 1 F α β ∂ ν F α β ] .
第2項は μ , α \mu,\alpha μ , α の反対称性で 1 2 F μ α ( ∂ μ F α ν + ∂ α F ν μ ) \frac12F^{\mu\alpha}(\partial_\mu F_{\alpha\nu}+\partial_\alpha F_{\nu\mu}) 2 1 F μα ( ∂ μ F α ν + ∂ α F ν μ ) と書け、Bianchi 恒等式でこれは − 1 2 F μ α ∂ ν F μ α -\frac12F^{\mu\alpha}\partial_\nu F_{\mu\alpha} − 2 1 F μα ∂ ν F μα になって第3項と打ち消し合う。残る第1項が
∂ μ T μ ν = − F ν λ J λ (2.7.25) \partial_\mu T^{\mu\nu}=-F^{\nu\lambda}J_\lambda\tag{2.7.25} ∂ μ T μν = − F ν λ J λ ( 2.7.25 )
である。右辺を成分で書くと、ν = 0 \nu=0 ν = 0 では − E ⋅ J / c -\bm{E}\cdot\bm{J}/c − E ⋅ J / c 、ν = i \nu=i ν = i では − ( ρ E + J × B ) i -(\rho\bm{E}+\bm{J}\times\bm{B})_i − ( ρ E + J × B ) i で、
∂ u ∂ t + ∇ ⋅ S = − J ⋅ E , ∂ g i ∂ t − ∑ j ∂ j σ i j = − ( ρ 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} ∂ t ∂ u + ∇ ⋅ S = − J ⋅ E , ∂ t ∂ g i − j ∑ ∂ j σ ij = − ( ρ E + J × B ) i ( 2.7.26 )
になる。第1式は Poynting の定理で、場のエネルギーは電流にした仕事 J ⋅ E \bm{J}\cdot\bm{E} J ⋅ E の分だけ減る。第2式の右辺は Lorentz 力の密度 f \bm{f} f で、物質の運動量密度の増える速さに等しい。したがって g \bm{g} g と物質の運動量密度の和は、応力 σ i j \sigma_{ij} σ ij の流れの分しか変わらない。場と物質を合わせた全運動量は保存する。
この保存が成り立つためには、場の運動量 ε 0 E × B \varepsilon_0\bm{E}\times\bm{B} ε 0 E × B を数に入れなければならない。動く2つの電荷の間の磁気力は作用反作用の法則を満たさないことがあり(例題 2.7.6 )、そのとき粒子だけの運動量は保存しない。第1部第10章 で全運動量の保存を並進の対称性から出したとき、系は粒子だけだった。電磁気では並進の対称性を持つのは粒子と場を合わせた作用なので、保存するのも粒子と場を合わせた運動量である。
例題 2.7.5 平面波の放射圧 ★★
オリジナル x x x の向きに進む平面波 E = E 0 cos ( k x − ω t ) e y \bm{E}=E_0\cos(kx-\omega t)\bm{e}_y E = E 0 cos ( k x − ω t ) e y 、B = ( E 0 / c ) cos ( k x − ω t ) e z \bm{B}=(E_0/c)\cos(kx-\omega t)\bm{e}_z B = ( E 0 / c ) cos ( k x − ω t ) e z を考える。
(1) T 00 = u T^{00}=u T 00 = u , S x S_x S x , T x x T^{xx} T xx を求め、S x = c u S_x=cu S x = c u 、T x x = u T^{xx}=u T xx = u を示せ。
(2) この波が x = 0 x=0 x = 0 の面に置いた板に垂直に当たり、すべて吸収されるとき、板が単位面積あたり受ける力を時間平均の強度 I = ⟨ S x ⟩ I=\langle S_x \rangle I = ⟨ S x ⟩ で書け。完全に反射されるときはどうか。
(3) 地球の軌道での太陽光の強度 I = 1360.8 W m − 2 I=1360.8\,\mathrm{W\,m^{-2}} I = 1360.8 W m − 2 1 1 G. Kopp and J. L. Lean (2011), Geophys. Res. Lett. 38, L01706(全太陽放射照度 1360.8 ± 0.5 W m − 2 1360.8\pm0.5\,\mathrm{W\,m^{-2}} 1360.8 ± 0.5 W m − 2 )。 で、吸収するときの圧力と、1 m 2 1\,\mathrm{m^2} 1 m 2 の板が1日(86400 s 86400\,\mathrm s 86400 s )に受け取る運動量を求めよ。
解答 (1) c 2 = 1 / ( ε 0 μ 0 ) c^2=1/(\varepsilon_0\mu_0) c 2 = 1/ ( ε 0 μ 0 ) から B 2 2 μ 0 = E 0 2 cos 2 2 μ 0 c 2 = ε 0 E 0 2 cos 2 2 \frac{B^2}{2\mu_0}=\frac{E_0^2\cos^2}{2\mu_0c^2}=\frac{\varepsilon_0E_0^2\cos^2}{2} 2 μ 0 B 2 = 2 μ 0 c 2 E 0 2 c o s 2 = 2 ε 0 E 0 2 c o s 2 で、電場と磁場のエネルギーは等しく u = ε 0 E 0 2 cos 2 ( k x − ω t ) u=\varepsilon_0E_0^2\cos^2(kx-\omega t) u = ε 0 E 0 2 cos 2 ( k x − ω t ) 。E × B = E y B z e x \bm{E}\times\bm{B}=E_yB_z\bm{e}_x E × B = E y B z e x から S x = E 0 2 cos 2 / ( μ 0 c ) = c ε 0 E 0 2 cos 2 = c u S_x=E_0^2\cos^2/(\mu_0c)=c\,\varepsilon_0E_0^2\cos^2=cu S x = E 0 2 cos 2 / ( μ 0 c ) = c ε 0 E 0 2 cos 2 = c u 。応力は E x = B x = 0 E_x=B_x=0 E x = B x = 0 なので σ x x = − ε 0 2 E y 2 − 1 2 μ 0 B z 2 = − u \sigma_{xx}=-\frac{\varepsilon_0}{2}E_y^2-\frac{1}{2\mu_0}B_z^2=-u σ xx = − 2 ε 0 E y 2 − 2 μ 0 1 B z 2 = − u で、T x x = − σ x x = u T^{xx}=-\sigma_{xx}=u T xx = − σ xx = u 。
(2) T x x T^{xx} T xx は「x x x 成分の運動量が x x x の向きに流れる密度」なので、板に毎秒・単位面積あたり流れ込む運動量は ⟨ T x x ⟩ = ⟨ u ⟩ = I / c \langle T^{xx} \rangle =\langle u \rangle =I/c ⟨ T xx ⟩ = ⟨ u ⟩ = I / c である。吸収されるとこの運動量がすべて板に移り、圧力は I / c I/c I / c になる。反射されると、入ってきた運動量 I / c I/c I / c に加えて逆向きに出ていく波の運動量 I / c I/c I / c の分も板が受け持つので、圧力は 2 I / c 2I/c 2 I / c になる。
(3) I / c = 1360.8 / ( 2.998 × 10 8 ) = 4.54 × 10 − 6 P a I/c=1360.8/(2.998\times10^8)=4.54\times10^{-6}\,\mathrm{Pa} I / c = 1360.8/ ( 2.998 × 1 0 8 ) = 4.54 × 1 0 − 6 Pa (反射なら 9.08 × 10 − 6 P a 9.08\times10^{-6}\,\mathrm{Pa} 9.08 × 1 0 − 6 Pa )。1 m 2 1\,\mathrm{m^2} 1 m 2 の板が受ける力は 4.54 × 10 − 6 N 4.54\times10^{-6}\,\mathrm N 4.54 × 1 0 − 6 N で、1日に受け取る運動量は 4.54 × 10 − 6 × 86400 = 0.392 k g m s − 1 4.54\times10^{-6}\times86400=0.392\,\mathrm{kg\,m\,s^{-1}} 4.54 × 1 0 − 6 × 86400 = 0.392 kg m s − 1 である。
例題 2.7.6 作用反作用を満たさない磁気力 ★★
オリジナル 電荷 q 1 q_1 q 1 が原点を + x +x + x の向きに速さ v v v で、電荷 q 2 q_2 q 2 が点 ( 0 , d , 0 ) (0,d,0) ( 0 , d , 0 ) を + y +y + y の向きに速さ u u u で通過する瞬間を考える。v , u ≪ c v,u\ll c v , u ≪ c とし、点電荷 q q q が速度 v \bm{v} v で動くときにつくる磁場を、最低次の B = μ 0 4 π q v × e r r 2 \bm{B}=\frac{\mu_0}{4\pi}\frac{q\,\bm{v}\times\bm{e}_r}{r^2} B = 4 π μ 0 r 2 q v × e r (e r \bm{e}_r e r は電荷から測る点への単位ベクトル)で近似する。
(1) q 2 q_2 q 2 が受ける磁気力と、q 1 q_1 q 1 が受ける磁気力を求めよ。
(2) q 1 = q 2 = 1.0 μ C q_1=q_2=1.0\,\mu\mathrm C q 1 = q 2 = 1.0 μ C , d = 0.10 m d=0.10\,\mathrm m d = 0.10 m , v = u = 1.0 × 10 3 m / s v=u=1.0\times10^3\,\mathrm{m/s} v = u = 1.0 × 1 0 3 m/s で、磁気力と Coulomb 力の大きさを比べよ。
(3) 粒子の運動量の和はこの瞬間にどう変わるか。保存則とどう折り合うかを述べよ。
解答 (1) q 1 q_1 q 1 から q 2 q_2 q 2 への単位ベクトルは e y \bm{e}_y e y で、q 2 q_2 q 2 の位置での磁場は B 1 = μ 0 q 1 v 4 π d 2 e x × e y = μ 0 q 1 v 4 π d 2 e z \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 B 1 = 4 π d 2 μ 0 q 1 v e x × e y = 4 π d 2 μ 0 q 1 v e z 。q 2 q_2 q 2 が受ける力は q 2 u e y × B 1 = μ 0 q 1 q 2 v u 4 π d 2 e x q_2u\,\bm{e}_y\times\bm{B}_1=\frac{\mu_0q_1q_2vu}{4\pi d^2}\bm{e}_x q 2 u e y × B 1 = 4 π d 2 μ 0 q 1 q 2 v u e x である。q 2 q_2 q 2 から q 1 q_1 q 1 への単位ベクトルは − e y -\bm{e}_y − e y で、q 2 q_2 q 2 の速度 u e y u\,\bm{e}_y u e y と平行なので B 2 = 0 \bm{B}_2=\bm{0} B 2 = 0 、q 1 q_1 q 1 が受ける磁気力は 0 0 0 である(図 2.7.3 )。
(2) 磁気力は μ 0 4 π q 1 q 2 v u d 2 = 10 − 7 × 10 − 12 × 10 6 10 − 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 4 π μ 0 d 2 q 1 q 2 v u = 1 0 − 7 × 1 0 − 2 1 0 − 12 × 1 0 6 = 1.0 × 1 0 − 11 N 、Coulomb 力は 1 4 π ε 0 q 1 q 2 d 2 = 8.99 × 10 9 × 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 4 π ε 0 1 d 2 q 1 q 2 = 8.99 × 1 0 9 × 1 0 − 12 /1 0 − 2 = 0.899 N 。比は v u / c 2 = 1.1 × 10 − 11 vu/c^2=1.1\times10^{-11} v u / c 2 = 1.1 × 1 0 − 11 で、磁気力は電気力に比べて ( v / c ) 2 (v/c)^2 ( v / c ) 2 の程度に小さい。
(3) Coulomb 力は(この近似で)大きさが等しく向きが逆なので、粒子の運動量の和の変化は磁気力だけで決まり、μ 0 q 1 q 2 v u 4 π d 2 e x \frac{\mu_0q_1q_2vu}{4\pi d^2}\bm{e}_x 4 π d 2 μ 0 q 1 q 2 v u e x の速さで増える。粒子だけの運動量は保存しない。式 (2.7.26) によれば、その分だけ場の運動量 ∫ ε 0 E × B d 3 x \int\varepsilon_0\bm{E}\times\bm{B}\,\mathrm{d}^3x ∫ ε 0 E × B d 3 x が減り、粒子と場を合わせた運動量は保存する。
図 2.7.3 例題 2.7.6 の配置。q 1 q_1 q 1 の磁場は q 2 q_2 q 2 の位置で紙面の手前向き(⊙ \odot ⊙ )で、q 2 q_2 q 2 は + x +x + x の向きの磁気力を受ける。q 1 q_1 q 1 は q 2 q_2 q 2 の速度の延長線上にあるので q 2 q_2 q 2 の磁場を受けず、磁気力は片方にしか働かない
近似の裏付け
型 B(極限) | 小さいもの v / c → 0 v/c\to0 v / c → 0 (電荷の速さが光速に比べて小さい。磁場は Biot–Savart の最低次) | 証明 要確認 | 行き先 近似の数理 第1章
1 G. Kopp and J. L. Lean (2011), Geophys. Res. Lett. 38, L01706(全太陽放射照度 1360.8 ± 0.5 W m − 2 1360.8\pm0.5\,\mathrm{W\,m^{-2}} 1360.8 ± 0.5 W m − 2 )。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com