tecktom.com

出てくるもの(1):Gauss の法則と Ampère–Maxwell の法則

2026年10月15日8分で読めますAI生成教科書『虎の巻』第2部 第6章 電磁場の Lagrangian と Maxwell 方程式
目次
  1. 4本の Euler–Lagrange 方程式
  2. 方程式が解を持つための条件:電荷の保存
  3. 木下の層

AI 生成・木下が未勉強

解析力学から様々な理論分野を組み立てる 第2部 第6章 §2.6.3 ★★★ Euler–Lagrange 方程式の4本は ∂μFμν=μ0Jν\partial_\mu F^{\mu\nu}=\mu_0J^\nu。ν=0\nu=0 が Gauss、ν=i\nu=i が変位電流つきの Ampère。FF の反対称性から電荷保存が要る

4本の Euler–Lagrange 方程式

式 (2.6.10) と式 (2.6.12) を式 (2.6.9) に入れると∂μ(−Fμν/μ0)−(−Jν)=0\partial_\mu(-F^{\mu\nu}/\mu_0)-(-J^\nu)=0、すなわち

  ∂μFμν=μ0Jν  (2.6.13)\boxed{\;\partial_\mu F^{\mu\nu}=\mu_0J^\nu\;}\tag{2.6.13}

である。これを3次元の量に戻す。

ν=0\nu=0 では、F00=0F^{00}=0 と Fi0=−F0i=Ei/cF^{i0}=-F^{0i}=E^i/c から

∂μFμ0=∑i∂iEic=∇⋅Ec=μ0cρ.\partial_\mu F^{\mu0}=\sum_i\partial_i\frac{E^i}{c}=\frac{\nabla\cdot\bm{E}}{c}=\mu_0c\rho .

両辺に cc を掛け、μ0c2=1/ε0\mu_0c^2=1/\varepsilon _0 を使うと

∇⋅E=ρε0(2.6.14)\nabla\cdot\bm{E}=\frac{\rho}{\varepsilon _0}\tag{2.6.14}

で、Gauss の法則である。前節では、この式の右辺が ρ/ε0\rho/\varepsilon _0 になるように係数を 1/μ01/\mu_0 に取った。

ν=j\nu=j では、F0j=−Ej/cF^{0j}=-E^j/c、∂0=1c∂t\partial_0=\frac1c\partial_t、Fij=Fij=−∑kϵijkBkF^{ij}=F_{ij}=-\sum_k\epsilon_{ijk}B_k から

∂μFμj=−1c2∂Ej∂t−∑i,kϵijk∂iBk=−1c2∂Ej∂t+(∇×B)j\partial_\mu F^{\mu j}=-\frac{1}{c^2}\frac{\partial E^j}{\partial t}-\sum_{i,k}\epsilon_{ijk}\partial_iB_k=-\frac{1}{c^2}\frac{\partial E^j}{\partial t}+(\nabla\times\bm{B})_j

(ϵijk=−ϵjik\epsilon_{ijk}=-\epsilon_{jik} を使った)。よって

∇×B−1c2∂E∂t=μ0J(2.6.15)\nabla\times\bm{B}-\frac{1}{c^2}\frac{\partial \bm{E}}{\partial t}=\mu_0\bm{J}\tag{2.6.15}

で、変位電流 ε0∂tE\varepsilon _0\partial_t\bm{E} を含む Ampère–Maxwell の法則である。

学部の電磁気学では、Ampère の法則 ∇×B=μ0J\nabla\times\bm{B}=\mu_0\bm{J} に変位電流の項を後から足す。ここでは、変位電流の項は ∂μFμj\partial_\mu F^{\mu j} の μ=0\mu=0 の項として最初から入っている。Lorentz 変換で変わらない L\mathcal L を選んだ時点で、ν=j\nu=j の式は時間微分と空間微分を ∂μ\partial_\mu の形でそろって含むしかない。

方程式が解を持つための条件:電荷の保存

式 (2.6.13) の両辺に ∂ν\partial_\nu を当てる。左辺は

∂ν∂μFμν=0\partial_\nu\partial_\mu F^{\mu\nu}=0

である。∂ν∂μ\partial_\nu\partial_\mu は μ,ν\mu,\nu について対称で、FμνF^{\mu\nu} は反対称なので、和を取ると打ち消し合う(添字の名前 μ↔ν\mu\leftrightarrow\nu を付け替えると全体が −1-1 倍になる)。したがって右辺も

∂νJν=∂ρ∂t+∇⋅J=0(2.6.16)\partial_\nu J^\nu=\frac{\partial \rho}{\partial t}+\nabla\cdot\bm{J}=0\tag{2.6.16}

でなければならない。電荷が保存しない源 JμJ^\mu を入れると、式 (2.6.13) は解を持たない。

式 (2.6.16) は、源 JμJ^\mu を外から与えた量とみなす限り、場の方程式から出た「結果」でなく、場の方程式を立てるための「前提」である。同じ条件は作用の側でも要る。ゲージ変換 Aμ→Aμ−∂μΛA_\mu\to A_\mu-\partial_\mu\Lambda で相互作用の項は

−JμAμ→−JμAμ+Jμ∂μΛ=−JμAμ+∂μ(JμΛ)−(∂μJμ)Λ-J^\mu A_\mu\to-J^\mu A_\mu+J^\mu\partial_\mu\Lambda=-J^\mu A_\mu+\partial_\mu(J^\mu\Lambda)-(\partial_\mu J^\mu)\Lambda

と変わり、∂μJμ=0\partial_\mu J^\mu=0 のときだけ全微分の差で済む(例題 2.6.1 (3) と同じ理由で方程式が変わらない)。電荷の保存とゲージ変換に対する不変性の関係を Noether の定理で読むのは第2部第7章である。

Euler–Lagrange 方程式から出る2本

L=−14μ0FμνFμν−JμAμ\mathcal L=-\frac{1}{4\mu_0}F_{\mu\nu}F^{\mu\nu}-J^\mu A_\mu の AνA_\nu についての Euler–Lagrange 方程式は ∂μFμν=μ0Jν\partial_\mu F^{\mu\nu}=\mu_0J^\nu の4本で、ν=0\nu=0 が Gauss の法則 ∇⋅E=ρ/ε0\nabla\cdot\bm{E}=\rho/\varepsilon _0、ν=1,2,3\nu=1,2,3 が Ampère–Maxwell の法則 ∇×B−1c2∂tE=μ0J\nabla\times\bm{B}-\frac{1}{c^2}\partial_t\bm{E}=\mu_0\bm{J} である。FμνF^{\mu\nu} の反対称性から、源は ∂νJν=0\partial_\nu J^\nu=0 を満たさなければならない。

例題 2.6.3 Gauss 分布をした電荷の電場 ★★

オリジナル原点のまわりに、全電荷 qq が

ρ(r)=q(2π)3/2a3 e−r2/(2a2)\rho(r)=\frac{q}{(2\pi)^{3/2}a^3}\,\mathrm{e}^{-r^2/(2a^2)}

と分布している(aa は広がりの長さ)。時間によらず、電流は 00 とする。

  1. (1) 式 (2.6.14) を半径 rr の球で積分して、電場の動径成分 Er(r)E_r(r) を求めよ。
  2. (2) ErE_r を同じ qq の点電荷の Coulomb 場と比べ、r=ar=a と r=3ar=3a での比を求めよ。ErE_r が最大になる rr はどこか。
  3. (3) r≪ar\ll a での ErE_r の振る舞いを求めよ。
解答

(1) 球対称なので E=Er(r) r/r\bm{E}=E_r(r)\,\bm{r}/r。式 (2.6.14) を半径 rr の球で積分し、左辺に Gauss の定理を使うと4πr2Er=qin(r)/ε04\pi r^2E_r=q_{\mathrm{in}}(r)/\varepsilon _0。球の中の電荷は、s=r′/as=r'/a と置いて

qin(r)=∫0r4πr′2ρ(r′) dr′=q2π∫0r/as2e−s2/2ds=q[erf⁡(r2a)−2π ra e−r2/(2a2)]q_{\mathrm{in}}(r)=\int_0^r4\pi r'^2\rho(r')\,\mathrm{d} r'=q\sqrt{\frac2\pi}\int_0^{r/a}s^2\mathrm{e}^{-s^2/2}\mathrm{d} s =q\Bigl[\operatorname{erf}\Bigl(\frac{r}{\sqrt2a}\Bigr)-\sqrt{\frac2\pi}\,\frac ra\,\mathrm{e}^{-r^2/(2a^2)}\Bigr]

(部分積分 ∫s2e−s2/2ds=−se−s2/2+∫e−s2/2ds\int s^2\mathrm{e}^{-s^2/2}\mathrm{d} s=-s\mathrm{e}^{-s^2/2}+\int\mathrm{e}^{-s^2/2}\mathrm{d} s を使った)。r→∞r\to\infty で括弧は 11 になり、全電荷が qq であることも確かめられる。よって

Er(r)=q4πε0r2[erf⁡(r2a)−2π ra e−r2/(2a2)].E_r(r)=\frac{q}{4\pi\varepsilon _0r^2}\Bigl[\operatorname{erf}\Bigl(\frac{r}{\sqrt2a}\Bigr)-\sqrt{\frac2\pi}\,\frac ra\,\mathrm{e}^{-r^2/(2a^2)}\Bigr].

(2) Coulomb 場 q/(4πε0r2)q/(4\pi\varepsilon _0r^2) との比は括弧の中身そのもので、r=ar=a で erf⁡(1/2)−2/π e−1/2=0.683−0.484=0.199\operatorname{erf}(1/\sqrt2)-\sqrt{2/\pi}\,\mathrm{e}^{-1/2}=0.683-0.484=0.199、r=3ar=3a で 0.9710.971 である。広がり aa の3倍離れると、電荷の分布の形はほとんど見えなくなる。ErE_r を q/(4πε0a2)q/(4\pi\varepsilon _0a^2) を単位にして数値で調べると、r≈1.37ar\approx1.37a で最大値 0.2140.214 をとる(図 2.6.1)。

(3) r≪ar\ll a では ρ≈ρ(0)=q/((2π)3/2a3)\rho\approx\rho(0)=q/((2\pi)^{3/2}a^3) が一様なので、qin≈4π3r3ρ(0)q_{\mathrm{in}}\approx\frac{4\pi}{3}r^3\rho(0) で

Er≈ρ(0)3ε0 r=q3ε0(2π)3/2a3 r.E_r\approx\frac{\rho(0)}{3\varepsilon _0}\,r=\frac{q}{3\varepsilon _0(2\pi)^{3/2}a^3}\,r .

中心で Er=0E_r=0 になり、点電荷の 1/r21/r^2 の発散は起きない。

図

図 2.6.1 例題 2.6.3 の電場の動径成分。点線が Gauss 分布の電荷の厳密な電場、実線が同じ全電荷の点電荷の Coulomb 場(遠くでの近似)。中心付近では電場は rr に比例して 00 になり、r≈1.37ar\approx1.37a で最大になる。r≳3ar\gtrsim3a では2本の差は3%以下になる

木下の勉強ステータス

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

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

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