解析力学から様々な理論分野を組み立てる 第2部 第6章 §2.6.1 ★★★ 場の変数は E , B \bm{E},\bm{B} E , B でなく A μ A_\mu A μ 。Lagrangian 密度は − F 2 / ( 4 μ 0 ) − J ⋅ A -F^2/(4\mu_0)-J\cdot A − F 2 / ( 4 μ 0 ) − J ⋅ A で、3次元では ε 0 E 2 / 2 − B 2 / ( 2 μ 0 ) − ρ ϕ + J ⋅ A \varepsilon _0E^2/2-B^2/(2\mu_0)-\rho\phi+\bm{J}\cdot\bm{A} ε 0 E 2 /2 − B 2 / ( 2 μ 0 ) − ρϕ + J ⋅ A
学部の電磁気学では、Maxwell 方程式の4本を実験法則(Coulomb の法則・電磁誘導・Ampère の法則)から1本ずつ組み立て、最後に変位電流を足して完成させる。この章では順序を逆にする。最初に Lagrangian 密度を1つ選び、第2部第5章 の場の Euler–Lagrange 方程式を当てて、4本がどこから出るかを確かめる。
選ぶ場は4元ポテンシャル A μ = ( ϕ / c , − A ) A_\mu=(\phi/c,-\bm{A}) A μ = ( ϕ / c , − A ) の4成分で、選ぶ Lagrangian 密度は
L = − 1 4 μ 0 F μ ν F μ ν − J μ A μ , F μ ν = ∂ μ A ν − ∂ ν A μ \mathcal L=-\frac{1}{4\mu_0}F_{\mu\nu}F^{\mu\nu}-J^\mu A_\mu,\qquad F_{\mu\nu}=\partial_\mu A_\nu-\partial_\nu A_\mu L = − 4 μ 0 1 F μν F μν − J μ A μ , F μν = ∂ μ A ν − ∂ ν A μ
である。当てる道具は場の Euler–Lagrange 方程式と、粒子の側では第2部第2章 の荷電粒子の Lagrangian である。出てくるものは3つある。1つ目は、Euler–Lagrange 方程式から Gauss の法則と Ampère–Maxwell の法則(変位電流つき)が出ること。2つ目は、残りの2本(∇ ⋅ B = 0 \nabla\cdot\bm{B}=0 ∇ ⋅ B = 0 と Faraday の法則)が方程式ではなく、E \bm{E} E と B \bm{B} B をポテンシャルから作ったことの帰結(恒等式)であること。3つ目は、粒子と場を1つの作用にまとめると Lorentz 力と Maxwell 方程式が同時に出て、真空では速さ c = 1 / ε 0 μ 0 c=1/\sqrt{\varepsilon _0\mu_0} c = 1/ ε 0 μ 0 の波が出ること。最後に、学部の授業で習う4本の形と、この章の形を式で並べる。
計量は η μ ν = d i a g ( + 1 , − 1 , − 1 , − 1 ) \eta_{\mu\nu}=\mathrm{diag}(+1,-1,-1,-1) η μν = diag ( + 1 , − 1 , − 1 , − 1 ) 、座標は 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 , ∇ ) とし、上下に1回ずつ現れるギリシャ添字について和を取る(この本の流儀。第2部第3章 )。3次元のラテン添字 i , j , k = 1 , 2 , 3 i,j,k=1,2,3 i , j , k = 1 , 2 , 3 には和の規約を使わず、和を取るときは ∑ \sum ∑ を書くかベクトル記法で書く。
場の変数を何に取るか
第2部第5章 では、時空の各点に値を持つ場 φ a ( x ) \varphi_a(x) φ a ( x ) を一般化座標とみなし、作用を
S [ φ ] = ∫ L d t d 3 x = 1 c ∫ L ( φ a , ∂ μ φ a ) d 4 x (2.6.1) S[\varphi]=\int\mathcal L\,\mathrm{d} t\,\mathrm{d}^3x=\frac1c\int\mathcal L\bigl(\varphi_a,\partial_\mu\varphi_a\bigr)\,\mathrm{d}^4x\tag{2.6.1} S [ φ ] = ∫ L d t d 3 x = c 1 ∫ L ( φ a , ∂ μ φ a ) d 4 x ( 2.6.1 )
と書いた。電磁場に当てるには、まず場の変数 φ a \varphi_a φ a を何に取るかを決める必要がある。
素直な候補は E \bm{E} E と B \bm{B} B の6成分である。しかしこの選び方では、第2部第2章 で荷電粒子の Lagrangian に現れた結合− q ϕ + q r ˙ ⋅ A -q\phi+q\dot{\bm{r}}\cdot\bm{A} − q ϕ + q r ˙ ⋅ A を書くことができない。粒子と場を1つの作用で扱うには、結合に現れる量を場の変数にする必要がある。そこで、粒子の場合(第2部第2章 )にすでに Lagrangian に現れていた量、すなわちスカラーポテンシャル ϕ \phi ϕ とベクトルポテンシャル A \bm{A} A を場の変数に取る。E \bm{E} E と B \bm{B} B はそこから
E = − ∇ ϕ − ∂ A ∂ t , B = ∇ × A (2.6.2) \bm{E}=-\nabla\phi-\frac{\partial \bm{A}}{\partial t},\qquad \bm{B}=\nabla\times\bm{A}\tag{2.6.2} E = − ∇ ϕ − ∂ t ∂ A , B = ∇ × A ( 2.6.2 )
で作る。4成分をまとめた A μ = ( ϕ / c , A ) A^\mu=(\phi/c,\bm{A}) A μ = ( ϕ / c , A ) が4元ポテンシャルで、添字を下げると A μ = ( ϕ / c , − A ) A_\mu=(\phi/c,-\bm{A}) A μ = ( ϕ / c , − A ) になる。
A μ A_\mu A μ の微分を反対称に組んだ
F μ ν ≔ ∂ μ A ν − ∂ ν A μ (2.6.3) F_{\mu\nu}\coloneqq\partial_\mu A_\nu-\partial_\nu A_\mu\tag{2.6.3} F μν : = ∂ μ A ν − ∂ ν A μ ( 2.6.3 )
を電磁場テンソルと呼ぶ。成分を書くと、F 0 i = 1 c ∂ t A i − ∂ i A 0 F_{0i}=\frac1c\partial_tA_i-\partial_iA_0 F 0 i = c 1 ∂ t A i − ∂ i A 0 で A i = − A i A_i=-A^i A i = − A i 、A 0 = ϕ / c A_0=\phi/c A 0 = ϕ / c だからF 0 i = − 1 c ( ∂ t A i + ∂ i ϕ ) = E i / c F_{0i}=-\frac1c(\partial_tA^i+\partial_i\phi)=E^i/c F 0 i = − c 1 ( ∂ t A i + ∂ i ϕ ) = E i / c となる。空間の成分は F 12 = ∂ 1 A 2 − ∂ 2 A 1 = − ( ∂ x A y − ∂ y A x ) = − B 3 F_{12}=\partial_1A_2-\partial_2A_1=-(\partial_xA_y-\partial_yA_x)=-B^3 F 12 = ∂ 1 A 2 − ∂ 2 A 1 = − ( ∂ x A y − ∂ y A x ) = − B 3 で、F 23 = − B 1 F_{23}=-B^1 F 23 = − B 1 、F 31 = − B 2 F_{31}=-B^2 F 31 = − B 2 も同じである。まとめると
F μ ν = ( 0 E x / c E y / c E z / c − E x / c 0 − B z B y − E y / c B z 0 − B x − E z / c − B y B x 0 ) , F 0 i = − E i / c , F i j = F i j (2.6.4) F_{\mu\nu}=\begin{pmatrix}0&E_x/c&E_y/c&E_z/c\\-E_x/c&0&-B_z&B_y\\-E_y/c&B_z&0&-B_x\\-E_z/c&-B_y&B_x&0\end{pmatrix},
\qquad F^{0i}=-E^i/c,\quad F^{ij}=F_{ij}\tag{2.6.4} F μν = 0 − E x / c − E y / c − E z / c E x / c 0 B z − B y E y / c − B z 0 B x E z / c B y − B x 0 , F 0 i = − E i / c , F ij = F ij ( 2.6.4 )
である(行が μ \mu μ 、列が ν \nu ν )。添字を2つとも上げると、時間の添字を1つだけ含む成分の符号が変わり、空間の添字2つの成分は変わらない。
Lagrangian 密度の形を決める条件
L \mathcal L L の形は、次の3つの条件でほぼ決まる。
1つ目は、作用が Lorentz 変換で変わらないこと(第2部第3章 で粒子の作用を固有時に取ったのと同じ要請)。d 4 x \mathrm{d}^4x d 4 x は Lorentz 変換で変わらないので、L \mathcal L L がスカラーであればよい。
2つ目は、場の方程式が A μ A_\mu A μ について線形になること。電磁場は重ね合わせの原理に従うので、L \mathcal L L は A μ A_\mu A μ とその微分の2次式に限る。さらに方程式が時間の2階微分までになるように、微分は1階までにする。
3つ目は、ゲージ変換 A μ → A μ − ∂ μ Λ A_\mu\to A_\mu-\partial_\mu\Lambda A μ → A μ − ∂ μ Λ (A → A + ∇ Λ \bm{A}\to\bm{A}+\nabla\Lambda A → A + ∇Λ 、ϕ → ϕ − ∂ t Λ \phi\to\phi-\partial_t\Lambda ϕ → ϕ − ∂ t Λ )でE \bm{E} E と B \bm{B} B が変わらないので、場だけの部分は F μ ν F_{\mu\nu} F μν を通してだけ A μ A_\mu A μ によること。F μ ν F_{\mu\nu} F μν がゲージ変換で変わらないことは、∂ μ ∂ ν Λ = ∂ ν ∂ μ Λ \partial_\mu\partial_\nu\Lambda=\partial_\nu\partial_\mu\Lambda ∂ μ ∂ ν Λ = ∂ ν ∂ μ Λ から1行で分かる。この条件がなぜ要るのか、という問いは第2部第8章 で扱う。
F μ ν F_{\mu\nu} F μν の2次でスカラーになる組は F μ ν F μ ν F_{\mu\nu}F^{\mu\nu} F μν F μν と ϵ μ ν ρ σ F μ ν F ρ σ \epsilon^{\mu\nu\rho\sigma}F_{\mu\nu}F_{\rho\sigma} ϵ μν ρ σ F μν F ρ σ の2つしかない。後者は全微分になって場の方程式に何も足さないこと(例題 2.6.1 )が分かるので、場だけの部分は F μ ν F μ ν F_{\mu\nu}F^{\mu\nu} F μν F μν の定数倍に決まる。成分で計算すると、式 (2.6.4) から
F μ ν F μ ν = 2 ∑ i F 0 i F 0 i + ∑ i , j F i j F i j = − 2 E 2 c 2 + 2 B 2 (2.6.5) F_{\mu\nu}F^{\mu\nu}=2\sum_iF_{0i}F^{0i}+\sum_{i,j}F_{ij}F^{ij}=-\frac{2E^2}{c^2}+2B^2\tag{2.6.5} F μν F μν = 2 i ∑ F 0 i F 0 i + i , j ∑ F ij F ij = − c 2 2 E 2 + 2 B 2 ( 2.6.5 )
である。前の係数は − 1 4 μ 0 -\frac{1}{4\mu_0} − 4 μ 0 1 に取る。符号は、E 2 E^2 E 2 の項(∂ t A \partial_t\bm{A} ∂ t A の2乗を含む項)が正になるように選ぶ。粒子の L = T − V L=T-V L = T − V で速度の2乗の項が正であるのと同じで、この符号を逆にするとエネルギーが下に有界でなくなる。大きさ 1 / μ 0 1/\mu_0 1/ μ 0 は、次の節で Gauss の法則が ∇ ⋅ E = ρ / ε 0 \nabla\cdot\bm{E}=\rho/\varepsilon _0 ∇ ⋅ E = ρ / ε 0 (Coulomb の法則)の形になるように決まる。
源との結合は、荷電粒子の作用 − q ∫ A μ d x μ -q\int A_\mu\,\mathrm{d} x^\mu − q ∫ A μ d x μ (第2部第2章 ・第2部第3章 )を連続的な電荷に広げた − J μ A μ -J^\mu A_\mu − J μ A μ に取る。ここで J μ = ( c ρ , J ) J^\mu=(c\rho,\bm{J}) J μ = ( c ρ , J ) は4元電流である。以上をまとめて、この章で選ぶ Lagrangian 密度は
L = − 1 4 μ 0 F μ ν F μ ν − J μ A μ (2.6.6) \boxed{\;\mathcal L=-\frac{1}{4\mu_0}F_{\mu\nu}F^{\mu\nu}-J^\mu A_\mu\;}\tag{2.6.6} L = − 4 μ 0 1 F μν F μν − J μ A μ ( 2.6.6 )
である。式 (2.6.5) と J μ A μ = c ρ ⋅ ϕ / c − J ⋅ A J^\mu A_\mu=c\rho\cdot\phi/c-\bm{J}\cdot\bm{A} J μ A μ = c ρ ⋅ ϕ / c − J ⋅ A 、ε 0 = 1 / ( μ 0 c 2 ) \varepsilon _0=1/(\mu_0c^2) ε 0 = 1/ ( μ 0 c 2 ) を使って3次元の量で書き直すと
L = ε 0 2 E 2 − 1 2 μ 0 B 2 − ρ ϕ + J ⋅ A = ε 0 2 ( ∇ ϕ + ∂ A ∂ t ) 2 − 1 2 μ 0 ( ∇ × A ) 2 − ρ ϕ + J ⋅ A (2.6.7) \mathcal L=\frac{\varepsilon _0}{2}E^2-\frac{1}{2\mu_0}B^2-\rho\phi+\bm{J}\cdot\bm{A}
=\frac{\varepsilon _0}{2}\Bigl(\nabla\phi+\frac{\partial \bm{A}}{\partial t}\Bigr)^2-\frac{1}{2\mu_0}(\nabla\times\bm{A})^2-\rho\phi+\bm{J}\cdot\bm{A}\tag{2.6.7} L = 2 ε 0 E 2 − 2 μ 0 1 B 2 − ρϕ + J ⋅ A = 2 ε 0 ( ∇ ϕ + ∂ t ∂ A ) 2 − 2 μ 0 1 ( ∇ × A ) 2 − ρϕ + J ⋅ A ( 2.6.7 )
となる。場の変数 ϕ , A \phi,\bm{A} ϕ , A で書くと、∂ t A \partial_t\bm{A} ∂ t A の2乗が第1項に入り、空間微分の2乗が第2項に入る。第1項は粒子の運動エネルギー、第2項はばねの位置エネルギーに当たり、第2部第4章 で原子の列から作った弦の場L = ρ 2 ( ∂ t u ) 2 − Y 2 ( ∂ x u ) 2 \mathcal L=\frac{\rho}{2}(\partial_tu)^2-\frac{Y}{2}(\partial_xu)^2 L = 2 ρ ( ∂ t u ) 2 − 2 Y ( ∂ x u ) 2 (ρ \rho ρ は線密度、Y Y Y は弾性率で、この節の電荷密度 ρ \rho ρ とは別の量)と同じ「時間微分の2乗 − - − 空間微分の2乗」の形をしている。
なぜ E , B \bm{E},\bm{B} E , B でなく ϕ , A \phi,\bm{A} ϕ , A なのか
作用の停留は「一般化座標をずらしても作用が1次で変わらない」という条件で、出てくる方程式の本数は一般化座標の数(場なら場の成分の数)と同じである。A μ A_\mu A μ は4成分なので、Euler–Lagrange 方程式は4本しか出ない。Maxwell 方程式は成分で数えると 1 + 3 + 1 + 3 = 8 1+3+1+3=8 1 + 3 + 1 + 3 = 8 本ある。残りの4本は、E , B \bm{E},\bm{B} E , B を式 (2.6.2) で作った時点で自動的に満たされる(2.6.4 節)。E , B \bm{E},\bm{B} E , B を変数に取ると、この4本の恒等式を方程式として別に課す必要が生じ、作用1つから全部を出すことができなくなる。
例題 2.6.1 E ⋅ B \bm{E}\cdot\bm{B} E ⋅ B を Lagrangian 密度に足しても方程式は変わらない ★★
オリジナル もう1つのスカラー ϵ μ ν ρ σ F μ ν F ρ σ \epsilon^{\mu\nu\rho\sigma}F_{\mu\nu}F_{\rho\sigma} ϵ μν ρ σ F μν F ρ σ (ϵ 0123 = + 1 \epsilon^{0123}=+1 ϵ 0123 = + 1 )について考える。
(1) ϵ μ ν ρ σ F μ ν F ρ σ = − 8 E ⋅ B / c \epsilon^{\mu\nu\rho\sigma}F_{\mu\nu}F_{\rho\sigma}=-8\,\bm{E}\cdot\bm{B}/c ϵ μν ρ σ F μν F ρ σ = − 8 E ⋅ B / c を示せ。
(2) 式 (2.6.2) を使ってE ⋅ B = − 1 2 ∂ ∂ t ( A ⋅ B ) − ∇ ⋅ ( ϕ B + 1 2 A × ∂ t A ) \bm{E}\cdot\bm{B}=-\frac12\frac{\partial }{\partial t}(\bm{A}\cdot\bm{B})-\nabla\cdot\bigl(\phi\bm{B}+\frac12\bm{A}\times\partial_t\bm{A}\bigr) E ⋅ B = − 2 1 ∂ t ∂ ( A ⋅ B ) − ∇ ⋅ ( ϕ B + 2 1 A × ∂ t A ) を示せ。
(3) 定数 θ \theta θ について L \mathcal L L に θ E ⋅ B \theta\,\bm{E}\cdot\bm{B} θ E ⋅ B を足しても、場の Euler–Lagrange 方程式が変わらない理由を述べよ。
解答 (1) ϵ μ ν ρ σ \epsilon^{\mu\nu\rho\sigma} ϵ μν ρ σ は添字がすべて異なるときだけ 0 0 0 でないので、( μ , ν ) (\mu,\nu) ( μ , ν ) と ( ρ , σ ) (\rho,\sigma) ( ρ , σ ) のどちらか一方に必ず 0 0 0 が入る。0 0 0 が ( μ , ν ) (\mu,\nu) ( μ , ν ) に入る場合、μ = 0 \mu=0 μ = 0 と ν = 0 \nu=0 ν = 0 の2通りが反対称性で同じ値になり、0 0 0 が ( ρ , σ ) (\rho,\sigma) ( ρ , σ ) に入る場合も同じ値なので、全体は4 ∑ i , j , k ϵ 0 i j k F 0 i F j k 4\sum_{i,j,k}\epsilon^{0ijk}F_{0i}F_{jk} 4 ∑ i , j , k ϵ 0 ij k F 0 i F j k になる。ϵ 0 i j k = ϵ i j k \epsilon^{0ijk}=\epsilon_{ijk} ϵ 0 ij k = ϵ ij k (3次元の Levi-Civita、ϵ 123 = + 1 \epsilon_{123}=+1 ϵ 123 = + 1 )で、式 (2.6.4) から F j k = − ∑ l ϵ j k l B l F_{jk}=-\sum_l\epsilon_{jkl}B_l F j k = − ∑ l ϵ j k l B l なので、∑ j , k ϵ i j k ϵ j k l = 2 δ i l \sum_{j,k}\epsilon_{ijk}\epsilon_{jkl}=2\delta_{il} ∑ j , k ϵ ij k ϵ j k l = 2 δ i l を使うと
4 ∑ i F 0 i ∑ j , k ϵ i j k F j k = 4 ∑ i E i c ( − 2 B i ) = − 8 c E ⋅ B . 4\sum_{i}F_{0i}\sum_{j,k}\epsilon_{ijk}F_{jk}=4\sum_i\frac{E_i}{c}(-2B_i)=-\frac{8}{c}\bm{E}\cdot\bm{B} . 4 i ∑ F 0 i j , k ∑ ϵ ij k F j k = 4 i ∑ c E i ( − 2 B i ) = − c 8 E ⋅ B .
(2) − ∇ ϕ ⋅ B = − ∇ ⋅ ( ϕ B ) + ϕ ∇ ⋅ B -\nabla\phi\cdot\bm{B}=-\nabla\cdot(\phi\bm{B})+\phi\nabla\cdot\bm{B} − ∇ ϕ ⋅ B = − ∇ ⋅ ( ϕ B ) + ϕ ∇ ⋅ B で、∇ ⋅ B = ∇ ⋅ ( ∇ × A ) = 0 \nabla\cdot\bm{B}=\nabla\cdot(\nabla\times\bm{A})=0 ∇ ⋅ B = ∇ ⋅ ( ∇ × A ) = 0 だから第2項は消える。残りの − ∂ t A ⋅ B -\partial_t\bm{A}\cdot\bm{B} − ∂ t A ⋅ B には、2つの恒等式
∂ ∂ t ( A ⋅ ∇ × A ) = ∂ t A ⋅ ∇ × A + A ⋅ ∇ × ∂ t A , ∇ ⋅ ( A × ∂ t A ) = ∂ t A ⋅ ∇ × A − A ⋅ ∇ × ∂ t A \frac{\partial }{\partial t}\bigl(\bm{A}\cdot\nabla\times\bm{A}\bigr)=\partial_t\bm{A}\cdot\nabla\times\bm{A}+\bm{A}\cdot\nabla\times\partial_t\bm{A},\qquad
\nabla\cdot(\bm{A}\times\partial_t\bm{A})=\partial_t\bm{A}\cdot\nabla\times\bm{A}-\bm{A}\cdot\nabla\times\partial_t\bm{A} ∂ t ∂ ( A ⋅ ∇ × A ) = ∂ t A ⋅ ∇ × A + A ⋅ ∇ × ∂ t A , ∇ ⋅ ( A × ∂ t A ) = ∂ t A ⋅ ∇ × A − A ⋅ ∇ × ∂ t A
を足して 2 ∂ t A ⋅ B = ∂ t ( A ⋅ B ) + ∇ ⋅ ( A × ∂ t A ) 2\,\partial_t\bm{A}\cdot\bm{B}=\partial_t(\bm{A}\cdot\bm{B})+\nabla\cdot(\bm{A}\times\partial_t\bm{A}) 2 ∂ t A ⋅ B = ∂ t ( A ⋅ B ) + ∇ ⋅ ( A × ∂ t A ) を得る。合わせると示す式になる。
(3) (2) により E ⋅ B = ∂ μ K μ \bm{E}\cdot\bm{B}=\partial_\mu K^\mu E ⋅ B = ∂ μ K μ の形(4元の発散)に書ける。発散を時空で積分すると、Gauss の定理で時空の領域の境界の上の積分になる。場の変分は境界で 0 0 0 に取るので、作用の変分 δ S \delta S δ S は θ E ⋅ B \theta\,\bm{E}\cdot\bm{B} θ E ⋅ B の項から何も受け取らない。これは第1部第4章 で、L L L に時間の全微分を足しても Euler–Lagrange 方程式が変わらなかったことの場の版である。直接計算しても、ϕ \phi ϕ と A x , A y , A z A_x,A_y,A_z A x , A y , A z の4本の Euler–Lagrange 方程式で θ E ⋅ B \theta\,\bm{E}\cdot\bm{B} θ E ⋅ B からの寄与は恒等的に 0 0 0 になる(検算スクリプトで確かめた)。したがって、場の方程式を決めるスカラーは F μ ν F μ ν F_{\mu\nu}F^{\mu\nu} F μν F μν だけになる。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com