# 定常電流と静磁場：ビオ・サバールの法則から rot B = μ0 J へ

> 電流密度と電荷保存則から出発し、ビオ・サバールの法則を経て静磁場の二つの基本方程式 div B = 0 と rot B = μ0 J を導く。回転（rot）を循環の面密度として定義し、アンペールの法則の積分形と微分形の同値性を証明する。
> https://rikai.mugen-giken.com/physics/electromagnetism/magnetostatics

## 0. この記事の要点

- 電流密度 $\boldsymbol{J}$ は「単位時間・単位面積あたりに通り抜ける電荷」を表すベクトル場です。電荷が消えも湧きもしないという要請は、局所的な形で $\nabla\cdot\boldsymbol{J} + \partial\rho/\partial t = 0$（連続の方程式）と書けます。
- 電流分布が時間によらない**定常電流**では $\nabla\cdot\boldsymbol{J} = 0$ が成り立ちます。この記事の議論はすべてこの条件の上に立っています。
- 定常電流のつくる磁場は**ビオ・サバールの法則**で与えられます。これは実験法則であり、静電場におけるクーロンの法則に対応する出発点です。
- ビオ・サバールの法則から、静磁場が満たす二つの微分方程式 $\nabla\cdot\boldsymbol{B} = 0$（磁気モノポールが存在しない）と $\nabla\times\boldsymbol{B} = \mu_0\boldsymbol{J}$（アンペールの法則）が**導かれます**。前者は $\boldsymbol{B}$ がベクトルポテンシャルの回転で書けることの帰結です。
- 回転 $\nabla\times$ は「単位面積あたりの循環」です。この読み方とストークスの定理により、微分形 $\nabla\times\boldsymbol{B} = \mu_0\boldsymbol{J}$ と積分形 $\oint_C\boldsymbol{B}\cdot d\boldsymbol{l} = \mu_0 I_{\text{enc}}$ が同値になります。
- $\nabla\times\boldsymbol{B} = \mu_0\boldsymbol{J}$ は定常でない場合に必ず破綻します。その破綻の仕方が、次章の変位電流を要求します。

## 1. 動機：エルステッドが見つけた「渦」

1820 年、コペンハーゲンのエルステッド (H. C. Ørsted) は、電流を流した導線の近くに置いた磁針が振れることに気づきました。当時、電気と磁気はまったく別の現象だと考えられていましたから、これは大きな驚きでした。さらに奇妙だったのは**振れ方**です。磁針は導線に引き寄せられも遠ざけられもせず、導線に**巻きつく向き**を向いたのです。

静電場では、電荷から放射状に電気力線が出ていました。ところが電流のまわりの場は、電流を軸として円を描きます。「湧き出し」ではなく「渦」なのです。この違いを数学的に言い切るのがこの記事の目標で、答えは次の二つの式に集約されます。

$$
\nabla\cdot\boldsymbol{B} = 0,\qquad \nabla\times\boldsymbol{B} = \mu_0\boldsymbol{J}.
$$

左の式は「湧き出しがない」、右の式は「渦の源は電流である」と読めます。この記事では、実験法則であるビオ・サバールの法則からこの二つを**証明**します。

<Figure caption="直線電流のまわりの磁場。磁力線は電流を軸とする同心円になり、距離に反比例して弱くなる。">
<svg viewBox="0 0 490 265" width="100%" role="img" aria-label="直線電流のまわりの同心円状の磁力線">
  <defs>
    <marker id="ms-arrow-a" markerWidth="10" markerHeight="10" refX="6.5" refY="3" orient="auto">
      <path d="M0,0 L7,3 L0,6 z" fill="currentColor" />
    </marker>
  </defs>
  <g fill="none" stroke="currentColor" strokeWidth="1.4" markerEnd="url(#ms-arrow-a)">
    <path d="M185,133 A45,45 0 1 0 184.3,140.8" />
    <path d="M215,133 A75,75 0 1 0 213.9,146.0" />
    <path d="M245,133 A105,105 0 1 0 243.4,151.2" />
  </g>
  <circle cx="140" cy="133" r="11" fill="none" stroke="var(--sl-color-accent)" strokeWidth="2" />
  <circle cx="140" cy="133" r="3" fill="var(--sl-color-accent)" />
  <text x="118" y="112" fill="var(--sl-color-accent)" fontSize="14" textAnchor="end">I</text>
  <g fill="currentColor" fontSize="13">
    <text x="285" y="105">B は電流を軸とする同心円</text>
    <text x="285" y="128">向きは右ねじの規則</text>
    <text x="285" y="151">大きさは距離 s に反比例</text>
    <text x="140" y="255" textAnchor="middle">電流 I は紙面の手前向き</text>
  </g>
</svg>
</Figure>

磁場 $\boldsymbol{B}$ そのものは、電荷に働く力によって定義します。電荷 $q$ が速度 $\boldsymbol{v}$ で運動しているとき、電場からの力に加えて

$$
\boldsymbol{F} = q\,\boldsymbol{v}\times\boldsymbol{B}
$$

という力を受けます。これがローレンツ力の磁気的な部分で、$\boldsymbol{B}$ の操作的な定義になっています。力が速度と直交するため磁場は仕事をしません。この点も静電場との大きな違いです。

## 2. 電流と電流密度、そして電荷保存則

### 2.1. 電流密度の定義

「電流 $I$ アンペア」という量は、導線という細い管を前提にした量です。導体の内部の場所ごとの流れを扱うには、場としての量が必要になります。

<Definition id="def-current-density" title="電流密度">
空間の各点 $\boldsymbol{r}$ で、電荷密度 $\rho(\boldsymbol{r}, t)$ の電荷が速度場 $\boldsymbol{v}(\boldsymbol{r}, t)$ で運動しているとき、**電流密度**を

$$
\boldsymbol{J}(\boldsymbol{r}, t) = \rho(\boldsymbol{r}, t)\,\boldsymbol{v}(\boldsymbol{r}, t)
$$

で定める。担い手が複数種類あるときは $\boldsymbol{J} = \sum_\alpha \rho_\alpha \boldsymbol{v}_\alpha$ とする。単位は $\mathrm{A/m^2}$ である。

有向曲面 $S$ を単位時間に通り抜ける電荷、すなわち $S$ を貫く**電流**は

$$
I_S = \int_S \boldsymbol{J}\cdot d\boldsymbol{S}
$$

で与えられる。
</Definition>

面積分が電流を与えることは、次のように確かめられます。面上の微小部分 $dS$ を、法線 $\boldsymbol{n}$ をもつ平面片とみなします。時間 $dt$ の間にこの面を通過する電荷は、底面 $dS$、母線 $\boldsymbol{v}\,dt$ の斜柱の中にあった電荷ですから、その体積 $(\boldsymbol{v}\cdot\boldsymbol{n})\,dt\,dS$ に $\rho$ を掛けたもの、つまり $\rho\,\boldsymbol{v}\cdot\boldsymbol{n}\,dS\,dt = \boldsymbol{J}\cdot\boldsymbol{n}\,dS\,dt$ です。これを $dt$ で割って $S$ 上で足し上げれば上式になります。面積分の扱いに不安があれば [重積分と累次積分](/mathematics/calculus/multiple-integrals) を参照してください。

### 2.2. 電荷保存則

電荷が生成も消滅もしないという実験事実を、場の言葉に翻訳します。まず、体積分がすべての領域で消えるなら被積分関数自身が消える、という当たり前に見えて頻繁に使う補題を用意します。

<Lemma id="lem-localization" title="局所化補題">
$\Omega \subset \mathbb{R}^3$ を開集合、$f:\Omega \to \mathbb{R}$ を連続関数とする。$\Omega$ に含まれるすべての閉球 $V$ について $\int_V f\,dV = 0$ が成り立つならば、$\Omega$ 上で $f \equiv 0$ である。
</Lemma>

<Proof of="lem-localization">
対偶を示す。ある点 $\boldsymbol{x}_0 \in \Omega$ で $f(\boldsymbol{x}_0) \neq 0$ だとする。必要なら $f$ を $-f$ に取り替えて $c := f(\boldsymbol{x}_0) > 0$ としてよい。

$f$ は $\boldsymbol{x}_0$ で連続だから、$\varepsilon = c/2$ に対してある $r > 0$ が存在し、$|\boldsymbol{x} - \boldsymbol{x}_0| < r$ かつ $\boldsymbol{x}\in\Omega$ ならば $|f(\boldsymbol{x}) - c| < c/2$、したがって $f(\boldsymbol{x}) > c/2$ となる。$\Omega$ は開集合なので、$r$ を小さく取り直せば閉球 $V = \{\boldsymbol{x} : |\boldsymbol{x}-\boldsymbol{x}_0| \le r/2\}$ を $\Omega$ に含めることができる。この $V$ 上で被積分関数は $c/2$ より大きいから

$$
\int_V f\,dV \ge \frac{c}{2}\cdot\frac{4}{3}\pi\left(\frac{r}{2}\right)^3 > 0
$$

となり、仮定に反する。
</Proof>

<Theorem id="thm-continuity" title="連続の方程式">
$\rho$ と $\boldsymbol{J}$ を $\mathbb{R}^3$ 上の $C^1$ 級の関数とし、任意の有界領域 $V$（境界 $\partial V$ は区分的に滑らかで、外向き法線をとる）について、電荷保存則

$$
\frac{d}{dt}\int_V \rho\,dV = -\oint_{\partial V}\boldsymbol{J}\cdot d\boldsymbol{S}
$$

が成り立つとする。このとき空間の各点で

$$
\nabla\cdot\boldsymbol{J} + \frac{\partial\rho}{\partial t} = 0
$$

が成り立つ。
</Theorem>

前提の積分形は「領域 $V$ の中の電荷が増えた分は、境界から流れ込んだ分に等しい」と読みます。右辺の負号は、外向きの流出を正に測っているためです。

<Proof of="thm-continuity">
左辺で、$\rho$ が $C^1$ 級で $V$ が時間によらない固定領域であることから、微分と積分を交換できて

$$
\frac{d}{dt}\int_V \rho\,dV = \int_V \frac{\partial\rho}{\partial t}\,dV .
$$

右辺に <Ref to="physics/electromagnetism/electrostatics#thm-divergence-theorem" text="ガウスの発散定理" /> を適用すると

$$
\oint_{\partial V}\boldsymbol{J}\cdot d\boldsymbol{S} = \int_V \nabla\cdot\boldsymbol{J}\,dV .
$$

したがって仮定は

$$
\int_V\left(\frac{\partial\rho}{\partial t} + \nabla\cdot\boldsymbol{J}\right)dV = 0
$$

と書き直せる。これが任意の有界領域、特に任意の閉球について成り立つ。被積分関数は $\rho, \boldsymbol{J}$ が $C^1$ 級であることから連続なので、<Ref to="lem-localization" /> によりそれは恒等的に $0$ である。
</Proof>

<Definition id="def-steady-current" title="定常電流">
電荷分布と電流分布がともに時間によらない、すなわち $\partial\rho/\partial t = 0$ かつ $\partial\boldsymbol{J}/\partial t = \boldsymbol{0}$ であるとき、その電流を**定常電流**という。
</Definition>

<Corollary id="cor-steady-divj" title="定常電流の条件">
定常電流に対しては

$$
\nabla\cdot\boldsymbol{J} = 0
$$

が成り立つ。すなわち電流密度の流線は途中で途切れず、有限の領域では閉じた回路をなす。
</Corollary>

<Proof of="cor-steady-divj">
<Ref to="thm-continuity" /> の連続の方程式に <Ref to="def-steady-current" /> の条件 $\partial\rho/\partial t = 0$ を代入すればただちに $\nabla\cdot\boldsymbol{J} = 0$ を得る。
</Proof>

<Aside type="caution">
$\nabla\cdot\boldsymbol{J} = 0$ は、この記事のすべての結果が乗っている土台です。「電流の切れ端」（有限の線分だけを流れる定常電流）は存在しません。線分の端に電荷がたまり続けてしまうからです。以下でこの条件をどこで使ったかは、そのつど明示します。
</Aside>

## 3. ビオ・サバールの法則

エルステッドの発見の直後、ビオ (J.-B. Biot) とサバール (F. Savart)、そしてアンペール (A.-M. Ampère) は、電流と磁場の定量的な関係を実験で決めました。その結果は次の形にまとめられます。これは他から導かれるものではなく、実験に基づく出発点です。

<Axiom id="ax-biot-savart" title="ビオ・サバールの法則">
有界な領域に分布する定常電流密度 $\boldsymbol{J}$（$\boldsymbol{J}$ は $C^1$ 級で、ある有界集合の外では $\boldsymbol{0}$）がつくる磁場は

$$
\boldsymbol{B}(\boldsymbol{r}) = \frac{\mu_0}{4\pi}\int_{\mathbb{R}^3} \frac{\boldsymbol{J}(\boldsymbol{r}')\times(\boldsymbol{r}-\boldsymbol{r}')}{|\boldsymbol{r}-\boldsymbol{r}'|^{3}}\,dV'
$$

で与えられる。ここで $\mu_0$ は真空の透磁率である。

電流が断面積の無視できる細い閉曲線 $C$ を流れる場合（線電流、電流 $I$）には、$\boldsymbol{J}\,dV' \to I\,d\boldsymbol{l}'$ と置き換えて

$$
\boldsymbol{B}(\boldsymbol{r}) = \frac{\mu_0 I}{4\pi}\oint_C \frac{d\boldsymbol{l}'\times(\boldsymbol{r}-\boldsymbol{r}')}{|\boldsymbol{r}-\boldsymbol{r}'|^{3}}
$$

となる。
</Axiom>

<Ref to="physics/electromagnetism/electrostatics#ax-coulomb" text="クーロンの法則" /> と比べると構造がよく見えます。距離依存性は同じ $1/R^2$（分母の $R^3$ のうち $R$ の一つは方向ベクトルの規格化ぶんです）。違いは、源がスカラー $\rho$ ではなくベクトル $\boldsymbol{J}$ であること、そして外積が入ることです。この外積が、磁場を「源のまわりを回る」向きに向けています。

<Remark id="rem-mu0" title="μ0 の値について">
2019 年の SI 基本単位の再定義より前は、$\mu_0 = 4\pi\times10^{-7}\ \mathrm{N/A^2}$ はアンペアの定義によって厳密に固定された値でした。再定義後は電気素量 $e$ が厳密値となったため、$\mu_0$ は測定量になっています。ただし値は $4\pi\times10^{-7}\ \mathrm{N/A^2}$ から相対的に $10^{-9}$ 程度しかずれておらず、学部の計算ではこの値をそのまま使って差し支えありません。
</Remark>

<Example id="ex-straight-wire" title="無限に長い直線電流">
$z$ 軸上を $+z$ 方向に電流 $I$ が流れているとする。円筒座標 $(s,\phi,z)$ の単位ベクトルを $\boldsymbol{e}_s,\boldsymbol{e}_\phi,\boldsymbol{e}_z$ と書き、場の点を $\boldsymbol{r} = s\,\boldsymbol{e}_s$（$z=0$ の平面上、一般性を失わない）にとる。源の点は $\boldsymbol{r}' = z'\boldsymbol{e}_z$、線素は $d\boldsymbol{l}' = dz'\,\boldsymbol{e}_z$ である。

$$
\boldsymbol{r}-\boldsymbol{r}' = s\,\boldsymbol{e}_s - z'\,\boldsymbol{e}_z,\qquad |\boldsymbol{r}-\boldsymbol{r}'| = \sqrt{s^2+z'^2}
$$

であり、$\boldsymbol{e}_z\times\boldsymbol{e}_s = \boldsymbol{e}_\phi$、$\boldsymbol{e}_z\times\boldsymbol{e}_z = \boldsymbol{0}$ だから

$$
d\boldsymbol{l}'\times(\boldsymbol{r}-\boldsymbol{r}') = dz'\,\boldsymbol{e}_z\times(s\,\boldsymbol{e}_s - z'\boldsymbol{e}_z) = s\,dz'\,\boldsymbol{e}_\phi .
$$

よって磁場は $\boldsymbol{e}_\phi$ 成分だけをもち、

$$
B_\phi = \frac{\mu_0 I}{4\pi}\int_{-\infty}^{\infty}\frac{s\,dz'}{(s^2+z'^2)^{3/2}} .
$$

積分は $z' = s\tan\theta$（$dz' = s\sec^2\theta\,d\theta$、$(s^2+z'^2)^{3/2} = s^3\sec^3\theta$）と置換すると

$$
\int_{-\infty}^{\infty}\frac{dz'}{(s^2+z'^2)^{3/2}} = \int_{-\pi/2}^{\pi/2}\frac{s\sec^2\theta}{s^3\sec^3\theta}\,d\theta = \frac{1}{s^2}\int_{-\pi/2}^{\pi/2}\cos\theta\,d\theta = \frac{2}{s^2}
$$

となる。したがって

$$
\boldsymbol{B} = \frac{\mu_0 I}{4\pi}\cdot s\cdot\frac{2}{s^2}\,\boldsymbol{e}_\phi = \frac{\mu_0 I}{2\pi s}\,\boldsymbol{e}_\phi .
$$

大きさは $s$ に反比例し、向きは電流に右ねじで巻きつきます。エルステッドの観測と一致します。
</Example>

<Example id="ex-circular-loop" title="円形電流の軸上の磁場">
半径 $a$ の円形回路が $xy$ 平面上、原点を中心に置かれ、上から見て反時計回りに電流 $I$ が流れているとする。$z$ 軸上の点 $\boldsymbol{r} = z\,\boldsymbol{e}_z$ での磁場を求める。

源の点を $\boldsymbol{r}' = a(\cos\phi',\sin\phi',0)$ とすると $\boldsymbol{r}-\boldsymbol{r}'$ の大きさは $R = \sqrt{a^2+z^2}$ で、$\phi'$ によらない。また線素 $d\boldsymbol{l}'$ は円の接線方向、$\boldsymbol{r}-\boldsymbol{r}'$ は $\boldsymbol{e}_z$ 成分と動径成分からなり、どちらも接線方向と直交する。よって

$$
|d\boldsymbol{l}'\times(\boldsymbol{r}-\boldsymbol{r}')| = R\,|d\boldsymbol{l}'| = R\,a\,d\phi' .
$$

この外積ベクトルは、円錐面上を一周する向きをもつ。$z$ 軸まわりの対称性から、一周積分すると動径成分は打ち消し合い、$z$ 成分だけが残る。$z$ 成分を取り出す因子は、外積ベクトルと $\boldsymbol{e}_z$ のなす角の余弦、すなわち $a/R$ である。したがって

$$
B_z = \frac{\mu_0 I}{4\pi}\int_0^{2\pi}\frac{R\,a\,d\phi'}{R^3}\cdot\frac{a}{R} = \frac{\mu_0 I}{4\pi}\cdot\frac{a^2}{R^3}\cdot 2\pi = \frac{\mu_0 I\,a^2}{2(a^2+z^2)^{3/2}} .
$$

中心 $z=0$ では $B_z = \mu_0 I/(2a)$、遠方 $|z| \gg a$ では

$$
B_z \simeq \frac{\mu_0 I a^2}{2|z|^3} = \frac{\mu_0}{4\pi}\cdot\frac{2m}{|z|^3},\qquad m := I\pi a^2
$$

となり、電気双極子の場と同じ $1/|z|^3$ の振る舞いを示す。$m$ を**磁気モーメント**と呼ぶ。
</Example>

## 4. 磁場に湧き出しはない

ビオ・サバールの法則から、まず $\nabla\cdot\boldsymbol{B} = 0$ を導きます。最短の道は、$\boldsymbol{B}$ が何かの回転として書けることを示すことです。

<Definition id="def-vector-potential" title="ベクトルポテンシャル">
<Ref to="ax-biot-savart" /> と同じ仮定のもとで

$$
\boldsymbol{A}(\boldsymbol{r}) := \frac{\mu_0}{4\pi}\int_{\mathbb{R}^3}\frac{\boldsymbol{J}(\boldsymbol{r}')}{|\boldsymbol{r}-\boldsymbol{r}'|}\,dV'
$$

を定常電流 $\boldsymbol{J}$ の**ベクトルポテンシャル**という。
</Definition>

これは <Ref to="physics/electromagnetism/electrostatics#def-potential" text="静電ポテンシャル" /> $\varphi(\boldsymbol{r}) = \frac{1}{4\pi\epsilon_0}\int \rho(\boldsymbol{r}')/|\boldsymbol{r}-\boldsymbol{r}'|\,dV'$ の各成分ごとの類似物です。以下では記号を短くするため

$$
\boldsymbol{R} := \boldsymbol{r}-\boldsymbol{r}',\qquad R := |\boldsymbol{R}|
$$

と置きます。$\nabla$ は場の点 $\boldsymbol{r}$ についての微分、$\nabla'$ は源の点 $\boldsymbol{r}'$ についての微分を表します。基本となる計算は

$$
\nabla\frac{1}{R} = -\frac{\boldsymbol{R}}{R^3},\qquad \nabla'\frac{1}{R} = +\frac{\boldsymbol{R}}{R^3} = -\nabla\frac{1}{R}
$$

です（$1/R$ は $\boldsymbol{r}$ と $\boldsymbol{r}'$ の差だけの関数なので、両者の微分は符号だけ違います）。

<Proposition id="prop-vector-potential" title="磁場はベクトルポテンシャルの回転">
<Ref to="ax-biot-savart" /> の磁場と <Ref to="def-vector-potential" /> のベクトルポテンシャルの間には

$$
\boldsymbol{B} = \nabla\times\boldsymbol{A}
$$

が成り立つ。
</Proposition>

<Proof of="prop-vector-potential">
微分と積分の交換が許されることを認めた上で（$\boldsymbol{J}$ は有界台をもち $1/R$ の特異性は $\mathbb{R}^3$ で可積分なので、積分は広義積分として絶対収束し、交換が正当化される）、

$$
\nabla\times\boldsymbol{A} = \frac{\mu_0}{4\pi}\int \nabla\times\left(\frac{\boldsymbol{J}(\boldsymbol{r}')}{R}\right)dV' .
$$

被積分関数にスカラーとベクトルの積の回転の公式 $\nabla\times(f\boldsymbol{a}) = f\,(\nabla\times\boldsymbol{a}) + (\nabla f)\times\boldsymbol{a}$ を使う。ここで $\boldsymbol{a} = \boldsymbol{J}(\boldsymbol{r}')$ は積分変数 $\boldsymbol{r}'$ の関数であって場の点 $\boldsymbol{r}$ には依存しないから $\nabla\times\boldsymbol{J}(\boldsymbol{r}') = \boldsymbol{0}$ であり、

$$
\nabla\times\left(\frac{\boldsymbol{J}(\boldsymbol{r}')}{R}\right) = \left(\nabla\frac1R\right)\times\boldsymbol{J}(\boldsymbol{r}') = -\frac{\boldsymbol{R}}{R^3}\times\boldsymbol{J}(\boldsymbol{r}') = \frac{\boldsymbol{J}(\boldsymbol{r}')\times\boldsymbol{R}}{R^3}
$$

を得る。最後の等号では外積の反交換性 $\boldsymbol{a}\times\boldsymbol{b} = -\boldsymbol{b}\times\boldsymbol{a}$ を使った。これを積分に戻すと、<Ref to="ax-biot-savart" /> の右辺そのものである。
</Proof>

<Lemma id="lem-div-curl" title="回転の発散は 0">
$C^2$ 級のベクトル場 $\boldsymbol{F}$ について $\nabla\cdot(\nabla\times\boldsymbol{F}) = 0$ が成り立つ。
</Lemma>

<Proof of="lem-div-curl">
成分で書き下す。

$$
\nabla\cdot(\nabla\times\boldsymbol{F}) = \frac{\partial}{\partial x}\left(\frac{\partial F_z}{\partial y}-\frac{\partial F_y}{\partial z}\right) + \frac{\partial}{\partial y}\left(\frac{\partial F_x}{\partial z}-\frac{\partial F_z}{\partial x}\right) + \frac{\partial}{\partial z}\left(\frac{\partial F_y}{\partial x}-\frac{\partial F_x}{\partial y}\right).
$$

$\boldsymbol{F}$ は $C^2$ 級なので、<Ref to="mathematics/calculus/multivariable-differentiation#thm-schwarz" text="シュワルツの定理" /> により二階偏微分の順序を交換できる。すると $\partial^2 F_z/\partial x\partial y$ と $-\partial^2 F_z/\partial y\partial x$ が打ち消し合い、同様に $F_x$ の項どうし、$F_y$ の項どうしも打ち消し合う。全体で $0$ である。
</Proof>

<Theorem id="thm-div-b" title="磁場の発散は 0">
<Ref to="ax-biot-savart" /> の仮定のもとで、空間のすべての点で

$$
\nabla\cdot\boldsymbol{B} = 0
$$

が成り立つ。同値に、任意の閉曲面 $S$ について $\oint_S \boldsymbol{B}\cdot d\boldsymbol{S} = 0$ である。
</Theorem>

<Proof of="thm-div-b">
<Ref to="prop-vector-potential" /> により $\boldsymbol{B} = \nabla\times\boldsymbol{A}$ と書ける。$\boldsymbol{A}$ は $\boldsymbol{J}$ が $C^1$ 級で有界台をもつことから $C^2$ 級であるので、<Ref to="lem-div-curl" /> を $\boldsymbol{F} = \boldsymbol{A}$ に適用して

$$
\nabla\cdot\boldsymbol{B} = \nabla\cdot(\nabla\times\boldsymbol{A}) = 0 .
$$

積分形は、この式を任意の有界領域 $V$ で積分してガウスの発散定理を使えばよい。逆に、任意の閉曲面で面積分が $0$ なら、発散定理と <Ref to="lem-localization" /> により $\nabla\cdot\boldsymbol{B} = 0$ が従う。
</Proof>

<Remark id="rem-monopole" title="磁気モノポールの不在">
静電場の <Ref to="physics/electromagnetism/electrostatics#thm-gauss-differential" text="ガウスの法則（微分形）" /> は $\nabla\cdot\boldsymbol{E} = \rho/\epsilon_0$ で、右辺に電荷密度という「湧き出しの源」がありました。<Ref to="thm-div-b" /> の右辺は $0$ です。これは「磁荷（磁気モノポール）が存在しない」ことを意味します。棒磁石を半分に折っても N 極だけの破片は得られず、断面に新しい S 極と N 極が現れます。

ここで導いた $\nabla\cdot\boldsymbol{B} = 0$ は、あくまで <Ref to="ax-biot-savart" /> の帰結です。すなわち「磁場の源は電流だけである」という前提から出ています。より根源的には、これは実験事実であり、ディラックは 1931 年に、磁気モノポールが宇宙に一つでも存在すれば電荷が離散的な値しか取れない（電荷の量子化）ことを示しました。モノポールの探索は現在も続いていますが、確立された検出例はありません。

$\boldsymbol{B} = \nabla\times\boldsymbol{A}$ という表し方は、$\boldsymbol{A}$ の取り方に自由度を残します。任意のスカラー場 $\chi$ について $\boldsymbol{A} \to \boldsymbol{A} + \nabla\chi$ と変えても、$\nabla\times\nabla\chi = \boldsymbol{0}$ なので $\boldsymbol{B}$ は変わりません。これが <Ref to="physics/electromagnetism/gauge-potentials#def-gauge-transformation" text="ゲージ変換" /> で、詳しくは [電磁ポテンシャルとゲージ変換](/physics/electromagnetism/gauge-potentials) で扱います。
</Remark>

## 5. 回転（rot）は循環の面密度である

もう一つの式 $\nabla\times\boldsymbol{B} = \mu_0\boldsymbol{J}$ に進む前に、回転という演算の意味を確かめます。発散が「単位体積あたりの湧き出し」だったのに対応して、回転は「単位面積あたりの循環」です。

<Definition id="def-curl" title="回転">
$C^1$ 級のベクトル場 $\boldsymbol{F} = (F_x, F_y, F_z)$ に対し、その**回転** $\nabla\times\boldsymbol{F}$（$\operatorname{rot}\boldsymbol{F}$ とも書く）を

$$
\nabla\times\boldsymbol{F} := \left(\frac{\partial F_z}{\partial y}-\frac{\partial F_y}{\partial z},\ \frac{\partial F_x}{\partial z}-\frac{\partial F_z}{\partial x},\ \frac{\partial F_y}{\partial x}-\frac{\partial F_x}{\partial y}\right)
$$

で定める。また、閉曲線 $C$ に沿った線積分 $\oint_C \boldsymbol{F}\cdot d\boldsymbol{l}$ を $\boldsymbol{F}$ の $C$ に沿った**循環**という。
</Definition>

この定義は、成分を眺めても意味がつかめません。次の命題が意味を与えます。

<Figure caption="回転の z 成分は、xy 平面内の微小長方形に沿った循環を面積で割ったもの。">
<svg viewBox="0 0 440 235" width="100%" role="img" aria-label="微小長方形に沿った循環と回転の関係を表す図">
  <defs>
    <marker id="ms-arrow-b" markerWidth="10" markerHeight="10" refX="6.5" refY="3" orient="auto">
      <path d="M0,0 L7,3 L0,6 z" fill="currentColor" />
    </marker>
  </defs>
  <g fill="none" stroke="currentColor" strokeWidth="1.6" markerEnd="url(#ms-arrow-b)">
    <path d="M120,165 L275,165" />
    <path d="M290,155 L290,70" />
    <path d="M280,60 L125,60" />
    <path d="M110,70 L110,155" />
  </g>
  <g fill="none" stroke="currentColor" strokeWidth="1" opacity="0.45">
    <path d="M110,60 L290,60 L290,165 L110,165 Z" />
  </g>
  <g fill="currentColor" fontSize="12.5">
    <text x="200" y="186" textAnchor="middle">Fx(x, y0)</text>
    <text x="200" y="46" textAnchor="middle">Fx(x, y0 + Δy)</text>
    <text x="300" y="117">Fy(x0 + Δx, y)</text>
    <text x="100" y="117" textAnchor="end">Fy(x0, y)</text>
  </g>
  <text x="200" y="118" fill="var(--sl-color-accent)" fontSize="13.5" textAnchor="middle">循環 ≈ (rot F)z ΔxΔy</text>
</svg>
</Figure>

<Proposition id="prop-curl-circulation" title="回転と循環">
$\boldsymbol{F}$ を $C^1$ 級のベクトル場とし、点 $(x_0,y_0,z_0)$ を左下の角にもつ $xy$ 平面に平行な長方形

$$
S_{\Delta} = [x_0, x_0+\Delta x]\times[y_0,y_0+\Delta y]\times\{z_0\}
$$

を考える。その境界 $\partial S_\Delta$ を、$z$ 軸の正の向きを法線とする右手系の向き（上から見て反時計回り）に一周する。このとき

$$
\lim_{\Delta x,\Delta y\to 0}\frac{1}{\Delta x\,\Delta y}\oint_{\partial S_\Delta}\boldsymbol{F}\cdot d\boldsymbol{l} = \left(\frac{\partial F_y}{\partial x}-\frac{\partial F_x}{\partial y}\right)\Bigg|_{(x_0,y_0,z_0)} = (\nabla\times\boldsymbol{F})_z\big|_{(x_0,y_0,z_0)}
$$

が成り立つ。
</Proposition>

<Proof of="prop-curl-circulation">
以下 $z = z_0$ を固定して省略する。四辺の線積分を書き下すと、下辺（$+x$ 向き）と上辺（$-x$ 向き）では $d\boldsymbol{l} = \pm dx\,\boldsymbol{e}_x$、右辺（$+y$ 向き）と左辺（$-y$ 向き）では $d\boldsymbol{l} = \pm dy\,\boldsymbol{e}_y$ だから

$$
\begin{aligned}
\oint_{\partial S_\Delta}\boldsymbol{F}\cdot d\boldsymbol{l}
&= \int_{x_0}^{x_0+\Delta x}F_x(x,y_0)\,dx + \int_{y_0}^{y_0+\Delta y}F_y(x_0+\Delta x, y)\,dy \\
&\quad - \int_{x_0}^{x_0+\Delta x}F_x(x,y_0+\Delta y)\,dx - \int_{y_0}^{y_0+\Delta y}F_y(x_0, y)\,dy .
\end{aligned}
$$

第 2 項と第 4 項、第 1 項と第 3 項をそれぞれまとめると

$$
\oint_{\partial S_\Delta}\boldsymbol{F}\cdot d\boldsymbol{l} = \int_{y_0}^{y_0+\Delta y}\big[F_y(x_0+\Delta x,y)-F_y(x_0,y)\big]dy - \int_{x_0}^{x_0+\Delta x}\big[F_x(x,y_0+\Delta y)-F_x(x,y_0)\big]dx .
$$

各角括弧に、$x$ 方向・$y$ 方向についての <Ref to="mathematics/calculus/mean-value-and-taylor#thm-mvt" text="平均値の定理" /> を適用する（$\boldsymbol{F}$ は $C^1$ 級）。ある $\xi \in (x_0, x_0+\Delta x)$、$\eta\in(y_0,y_0+\Delta y)$ が（$y$ あるいは $x$ ごとに）存在して

$$
F_y(x_0+\Delta x,y)-F_y(x_0,y) = \frac{\partial F_y}{\partial x}(\xi, y)\,\Delta x,\qquad
F_x(x,y_0+\Delta y)-F_x(x,y_0) = \frac{\partial F_x}{\partial y}(x, \eta)\,\Delta y .
$$

これを代入し、さらに残った積分に <Ref to="mathematics/calculus/integration-and-ftc#prop-mvt-integral" text="積分の平均値の定理" /> を使うと、ある $(\xi_1,\eta_1), (\xi_2,\eta_2)$ が長方形内に存在して

$$
\oint_{\partial S_\Delta}\boldsymbol{F}\cdot d\boldsymbol{l} = \left[\frac{\partial F_y}{\partial x}(\xi_1,\eta_1) - \frac{\partial F_x}{\partial y}(\xi_2,\eta_2)\right]\Delta x\,\Delta y .
$$

両辺を $\Delta x\,\Delta y$ で割り、$\Delta x,\Delta y\to 0$ とする。$(\xi_i,\eta_i)$ はいずれも $(x_0,y_0)$ に収束し、偏導関数は $C^1$ 級の仮定から連続なので、極限は $\partial_x F_y - \partial_y F_x$ の $(x_0,y_0)$ での値になる。これは <Ref to="def-curl" /> の $z$ 成分にほかならない。
</Proof>

つまり $(\nabla\times\boldsymbol{F})\cdot\boldsymbol{n}$ は「法線 $\boldsymbol{n}$ をもつ微小面のふちを一周したときの循環の面密度」です。渦の強さを測る量だ、と読めます。この局所的な関係を有限の面に積み上げたものが、次のストークスの定理です。

<Theorem id="thm-stokes" title="ストークスの定理">
$S$ を $\mathbb{R}^3$ 内の区分的に滑らかな有向曲面、$\partial S$ をその境界曲線とし、$\partial S$ の向きは $S$ の法線に対して右手系にとる。$\boldsymbol{F}$ を $S$ を含む開集合上の $C^1$ 級ベクトル場とすると

$$
\int_S (\nabla\times\boldsymbol{F})\cdot d\boldsymbol{S} = \oint_{\partial S}\boldsymbol{F}\cdot d\boldsymbol{l}
$$

が成り立つ。
</Theorem>

<Remark id="rem-stokes-proof" title="ストークスの定理の証明について">
証明は、曲面を微小な長方形（あるいは三角形）に分割し、各片に <Ref to="prop-curl-circulation" /> を適用したうえで足し合わせる、という筋道です。内部の辺は隣り合う片から逆向きに二度ずつ数えられて打ち消し合い、境界 $\partial S$ の辺だけが残ります。厳密な扱いは多変数微積分の教科書に譲ります（参考文献の杉浦、あるいは Jackson の付録を参照してください）。この記事では既知として使います。
</Remark>

## 6. アンペールの法則

準備が整いました。ビオ・サバールの法則から $\nabla\times\boldsymbol{B}$ を計算します。

<Lemma id="lem-div-a" title="定常電流のベクトルポテンシャルは横波的">
<Ref to="ax-biot-savart" /> の仮定（$\boldsymbol{J}$ は $C^1$ 級、有界台、かつ定常電流なので <Ref to="cor-steady-divj" /> により $\nabla\cdot\boldsymbol{J} = 0$）のもとで、<Ref to="def-vector-potential" /> のベクトルポテンシャルは

$$
\nabla\cdot\boldsymbol{A} = 0
$$

を満たす。
</Lemma>

<Proof of="lem-div-a">
$\boldsymbol{J}(\boldsymbol{r}')$ は $\boldsymbol{r}$ に依存しないから

$$
\nabla\cdot\boldsymbol{A} = \frac{\mu_0}{4\pi}\int \nabla\cdot\left(\frac{\boldsymbol{J}(\boldsymbol{r}')}{R}\right)dV' = \frac{\mu_0}{4\pi}\int \boldsymbol{J}(\boldsymbol{r}')\cdot\nabla\frac1R\,dV' .
$$

ここで $\nabla(1/R) = -\nabla'(1/R)$ を使うと

$$
\nabla\cdot\boldsymbol{A} = -\frac{\mu_0}{4\pi}\int \boldsymbol{J}(\boldsymbol{r}')\cdot\nabla'\frac1R\,dV' .
$$

$\boldsymbol{r}'$ についての積の微分公式 $\nabla'\cdot(f\boldsymbol{a}) = (\nabla' f)\cdot\boldsymbol{a} + f\,\nabla'\cdot\boldsymbol{a}$ を $f = 1/R$、$\boldsymbol{a} = \boldsymbol{J}$ に対して用いると

$$
\boldsymbol{J}\cdot\nabla'\frac1R = \nabla'\cdot\left(\frac{\boldsymbol{J}}{R}\right) - \frac{\nabla'\cdot\boldsymbol{J}}{R}
$$

なので、十分大きな半径の球 $V$（その外では $\boldsymbol{J} = \boldsymbol{0}$）をとって

$$
\nabla\cdot\boldsymbol{A} = -\frac{\mu_0}{4\pi}\oint_{\partial V}\frac{\boldsymbol{J}(\boldsymbol{r}')\cdot d\boldsymbol{S}'}{R} + \frac{\mu_0}{4\pi}\int_V \frac{\nabla'\cdot\boldsymbol{J}(\boldsymbol{r}')}{R}\,dV'
$$

を得る（第 1 項でガウスの発散定理を使った）。

第 1 項は、$\partial V$ 上で $\boldsymbol{J} = \boldsymbol{0}$ だから消える（ここで「電流分布が有界な領域に限られる」という仮定を使った）。第 2 項は <Ref to="cor-steady-divj" /> の $\nabla'\cdot\boldsymbol{J} = 0$ により消える（ここで**定常電流**の仮定を使った）。よって $\nabla\cdot\boldsymbol{A} = 0$ である。
</Proof>

<Theorem id="thm-ampere-differential" title="アンペールの法則（微分形）">
<Ref to="ax-biot-savart" /> の仮定（定常電流、$C^1$ 級、有界台）のもとで、空間のすべての点で

$$
\nabla\times\boldsymbol{B} = \mu_0\boldsymbol{J}
$$

が成り立つ。
</Theorem>

<Proof of="thm-ampere-differential">
<Ref to="prop-vector-potential" /> より $\boldsymbol{B} = \nabla\times\boldsymbol{A}$ だから、ベクトル解析の恒等式（Appendix で成分計算により確かめます）

$$
\nabla\times(\nabla\times\boldsymbol{A}) = \nabla(\nabla\cdot\boldsymbol{A}) - \nabla^2\boldsymbol{A}
$$

を用いる。ここで $\nabla^2\boldsymbol{A}$ は各デカルト成分にラプラシアンを施したものである。

第 1 項は <Ref to="lem-div-a" /> により $\nabla\cdot\boldsymbol{A} = 0$ なので消える。第 2 項を計算する。<Ref to="def-vector-potential" /> の定義式で微分と積分を交換して

$$
\nabla^2\boldsymbol{A} = \frac{\mu_0}{4\pi}\int \boldsymbol{J}(\boldsymbol{r}')\,\nabla^2\frac{1}{|\boldsymbol{r}-\boldsymbol{r}'|}\,dV' .
$$

[静電場とガウスの法則](/physics/electromagnetism/electrostatics) の <Ref to="physics/electromagnetism/electrostatics#thm-poisson" /> で確かめたとおり、$1/|\boldsymbol{r}-\boldsymbol{r}'|$ は $\boldsymbol{r} \ne \boldsymbol{r}'$ では調和関数であり、原点での特異性まで込めると超関数の意味で

$$
\nabla^2\frac{1}{|\boldsymbol{r}-\boldsymbol{r}'|} = -4\pi\,\delta^3(\boldsymbol{r}-\boldsymbol{r}')
$$

が成り立つ。これを代入すると

$$
\nabla^2\boldsymbol{A} = \frac{\mu_0}{4\pi}\int \boldsymbol{J}(\boldsymbol{r}')\cdot(-4\pi)\,\delta^3(\boldsymbol{r}-\boldsymbol{r}')\,dV' = -\mu_0\boldsymbol{J}(\boldsymbol{r}) .
$$

以上より

$$
\nabla\times\boldsymbol{B} = \nabla(\nabla\cdot\boldsymbol{A}) - \nabla^2\boldsymbol{A} = \boldsymbol{0} + \mu_0\boldsymbol{J} = \mu_0\boldsymbol{J}
$$

を得る。
</Proof>

<Corollary id="cor-ampere-integral" title="アンペールの法則（積分形）">
定常電流に対し、区分的に滑らかな有向曲面 $S$ とその境界閉曲線 $C = \partial S$（向きは $S$ の法線に対して右手系）について

$$
\oint_C \boldsymbol{B}\cdot d\boldsymbol{l} = \mu_0 I_{\text{enc}},\qquad I_{\text{enc}} := \int_S \boldsymbol{J}\cdot d\boldsymbol{S}
$$

が成り立つ。ここで $I_{\text{enc}}$ は $C$ を「くぐり抜ける」正味の電流である。
</Corollary>

<Proof of="cor-ampere-integral">
<Ref to="thm-stokes" /> を $\boldsymbol{F} = \boldsymbol{B}$、曲面 $S$ に適用すると

$$
\oint_C\boldsymbol{B}\cdot d\boldsymbol{l} = \int_S(\nabla\times\boldsymbol{B})\cdot d\boldsymbol{S} .
$$

右辺に <Ref to="thm-ampere-differential" /> を代入して

$$
\int_S(\nabla\times\boldsymbol{B})\cdot d\boldsymbol{S} = \mu_0\int_S\boldsymbol{J}\cdot d\boldsymbol{S} = \mu_0 I_{\text{enc}}
$$

となる。逆向きも同様で、任意の曲面について積分形が成り立てば、<Ref to="prop-curl-circulation" /> の極限操作により微分形が回復する。
</Proof>

<Remark id="rem-surface-independence" title="曲面の取り方によらないこと">
<Ref to="cor-ampere-integral" /> の左辺は曲線 $C$ だけで決まるのに、右辺は $C$ を境界とする曲面 $S$ を使って書かれています。これが矛盾しないのは、$\nabla\cdot\boldsymbol{J} = 0$ のおかげです。実際、同じ境界をもつ二つの曲面 $S_1, S_2$ をつなぐと閉曲面ができ、ガウスの発散定理より

$$
\int_{S_1}\boldsymbol{J}\cdot d\boldsymbol{S} - \int_{S_2}\boldsymbol{J}\cdot d\boldsymbol{S} = \int_{V}\nabla\cdot\boldsymbol{J}\,dV = 0
$$

となるからです（$V$ は二つの曲面が囲む領域）。定常でなくなるとこの一致は崩れ、アンペールの法則そのものが意味を失います。これが変位電流を要求する仕組みで、演習の最後の問題で確かめます。
</Remark>

### 6.1. 対称性を使った計算

積分形は、対称性の高い配置で $\boldsymbol{B}$ を求めるのに非常に強力です。手順は静電場のガウスの法則と同じで、まず対称性から $\boldsymbol{B}$ の形を絞り、次に都合のよい閉曲線を選びます。

<Example id="ex-cylinder" title="太さのある円柱導体">
半径 $a$ の無限に長い円柱導体が $z$ 軸に沿って置かれ、断面を一様な電流密度で全電流 $I$ が $+z$ 方向に流れている。すなわち $s \le a$ で $\boldsymbol{J} = \dfrac{I}{\pi a^2}\boldsymbol{e}_z$、$s > a$ で $\boldsymbol{J} = \boldsymbol{0}$ である。

**手順 1：場の形を絞る。** 配置は $z$ 方向の平行移動と $z$ 軸まわりの回転で不変だから、円筒座標成分 $B_s, B_\phi, B_z$ はいずれも $s$ だけの関数である。

$B_s$ については、半径 $s$、長さ $L$ の円筒面と両端の円板からなる閉曲面に <Ref to="thm-div-b" /> の積分形を適用する。両端の円板からの寄与は $B_z(s)$ が $z$ によらないため打ち消し合い、側面からの寄与は $B_s(s)\cdot 2\pi s L$ である。よって $B_s(s)\,2\pi s L = 0$、すなわち $B_s = 0$。

$B_z$ については、$s$ 軸と $z$ 軸を含む平面内の長方形 $[s_1,s_2]\times[0,L]$ を境界とする閉曲線に <Ref to="cor-ampere-integral" /> を適用する。この長方形をくぐる電流は $0$（電流は $z$ 方向で、長方形の法線は $\boldsymbol{e}_\phi$ 方向）なので $\big(B_z(s_1)-B_z(s_2)\big)L = 0$、つまり $B_z$ は $s$ によらない定数である。無限遠で場が消えることを要求すれば $B_z \equiv 0$。

残るのは $\boldsymbol{B} = B_\phi(s)\,\boldsymbol{e}_\phi$ である。

**手順 2：円形のループを取る。** $z$ 軸を中心とする半径 $s$ の円 $C$（反時計回り）を境界とする円板に <Ref to="cor-ampere-integral" /> を適用すると、左辺は $B_\phi(s)\cdot 2\pi s$ である。右辺の $I_{\text{enc}}$ は

$$
I_{\text{enc}} = \begin{cases} \dfrac{I}{\pi a^2}\cdot\pi s^2 = I\dfrac{s^2}{a^2} & (s \le a) \\[2mm] I & (s \ge a)\end{cases}
$$

したがって

$$
\boldsymbol{B} = \begin{cases} \dfrac{\mu_0 I\,s}{2\pi a^2}\,\boldsymbol{e}_\phi & (s\le a)\\[2mm] \dfrac{\mu_0 I}{2\pi s}\,\boldsymbol{e}_\phi & (s\ge a)\end{cases}
$$

内部では中心で $0$ から $s$ に比例して増え、$s=a$ で最大値 $\mu_0 I/(2\pi a)$ に達し、外部では $1/s$ で減衰します。$s = a$ で両式が一致すること、および外部の式が <Ref to="ex-straight-wire" /> と一致することを確かめてください。導線の外から見るかぎり、太さは見えません。
</Example>

<Example id="ex-solenoid" title="無限に長いソレノイド">
単位長さあたり $n$ 回の密な巻き数をもつ無限に長い円筒コイルに電流 $I$ が流れている。巻きが十分密で、電流は円筒面上を $\boldsymbol{e}_\phi$ 方向に面電流密度 $nI$ で流れているとみなす。

対称性から、<Ref to="ex-cylinder" /> と同じ議論で $B_s = 0$ が言える。また、いま電流は $\boldsymbol{e}_\phi$ 方向なので、$z$ 軸を中心とする円を境界とする円板をくぐる電流は $0$ であり、<Ref to="cor-ampere-integral" /> から $B_\phi\cdot 2\pi s = 0$、つまり $B_\phi = 0$。残るのは $\boldsymbol{B} = B_z(s)\boldsymbol{e}_z$ である。

$sz$ 平面内の長方形ループ（$z$ 方向の長さ $L$、$s$ 方向に $s_1$ から $s_2$）をとる。$\boldsymbol{B}$ が $z$ 成分のみなので、循環は $z$ 方向の二辺からだけ生じ、$\big(B_z(s_1)-B_z(s_2)\big)L$ となる。

- 両辺ともコイルの外側（$s_1, s_2 > $ 半径）に取ると電流をくぐらないので $B_z$ は外部で一定。無限遠で $0$ となることから外部では $B_z = 0$。
- 一方の辺を内側、他方を外側に取ると、ループは巻線を $nL$ 回くぐるので $I_{\text{enc}} = nLI$。よって $\big(B_z^{\text{in}} - 0\big)L = \mu_0 n L I$、すなわち

$$
\boldsymbol{B}^{\text{in}} = \mu_0 n I\,\boldsymbol{e}_z .
$$

- 両辺とも内側に取ると $I_{\text{enc}} = 0$ なので、内部で $B_z$ は $s$ によらない。

結論として、無限ソレノイドの内部は一様な磁場 $\mu_0 nI$、外部は $\boldsymbol{0}$ です。たとえば $n = 1000\ \mathrm{m^{-1}}$、$I = 2\ \mathrm{A}$ なら $B = 4\pi\times10^{-7}\times1000\times2 \approx 2.5\times10^{-3}\ \mathrm{T}$ となります。
</Example>

## 7. 静磁場のまとめと、次章への橋

定常電流に対して、この記事で得たのは次の二つです。

| | 微分形 | 積分形 | 意味 |
|---|---|---|---|
| 磁場の発散 | $\nabla\cdot\boldsymbol{B} = 0$ | $\oint_S \boldsymbol{B}\cdot d\boldsymbol{S} = 0$ | 磁荷は存在しない |
| 磁場の回転 | $\nabla\times\boldsymbol{B} = \mu_0\boldsymbol{J}$ | $\oint_C\boldsymbol{B}\cdot d\boldsymbol{l} = \mu_0 I_{\text{enc}}$ | 渦の源は電流 |

静電場の $\nabla\cdot\boldsymbol{E} = \rho/\epsilon_0$、$\nabla\times\boldsymbol{E} = \boldsymbol{0}$ と並べると、四つの式がそろいます。マクスウェル方程式まであと一歩です。

しかし、この記事の結果はすべて $\nabla\cdot\boldsymbol{J} = 0$ に依存していました。実際、$\nabla\times\boldsymbol{B} = \mu_0\boldsymbol{J}$ の両辺の発散を取ると、<Ref to="lem-div-curl" /> により左辺は恒等的に $0$ ですから

$$
0 = \mu_0\,\nabla\cdot\boldsymbol{J} = -\mu_0\frac{\partial\rho}{\partial t}
$$

となり、$\partial\rho/\partial t \ne 0$ の状況（コンデンサーの充電など）とは両立しません。この矛盾をどう解消するかが、[電磁誘導と変位電流](/physics/electromagnetism/induction-and-displacement-current) の主題であり、その先に [マクスウェル方程式と電磁波](/physics/electromagnetism/maxwell-equations) があります。

## 8. 演習

<Exercise id="exr-finite-wire" difficulty="標準">
$z$ 軸上の線分 $z_1 \le z \le z_2$ を $+z$ 方向に電流 $I$ が流れているとする（実際にはこれは回路の一部分であり、それだけでは定常電流にならないことに注意）。<Ref to="ax-biot-savart" /> を使って、点 $(s, 0, 0)$（$s > 0$）における磁場を求めよ。また $z_1 \to -\infty$、$z_2\to+\infty$ の極限で <Ref to="ex-straight-wire" /> の結果が回復することを確かめよ。

<Solution>
<Ref to="ex-straight-wire" /> の途中計算がそのまま使えます。$d\boldsymbol{l}'\times(\boldsymbol{r}-\boldsymbol{r}') = s\,dz'\,\boldsymbol{e}_\phi$ なので

$$
B_\phi = \frac{\mu_0 I s}{4\pi}\int_{z_1}^{z_2}\frac{dz'}{(s^2+z'^2)^{3/2}} .
$$

不定積分は

$$
\int\frac{dz'}{(s^2+z'^2)^{3/2}} = \frac{z'}{s^2\sqrt{s^2+z'^2}} + C
$$

です（右辺を $z'$ で微分すると、商の微分により $\dfrac{\sqrt{s^2+z'^2} - z'\cdot z'/\sqrt{s^2+z'^2}}{s^2(s^2+z'^2)} = \dfrac{s^2}{s^2(s^2+z'^2)^{3/2}} = \dfrac{1}{(s^2+z'^2)^{3/2}}$ となり確かめられます）。したがって

$$
\boldsymbol{B} = \frac{\mu_0 I}{4\pi s}\left(\frac{z_2}{\sqrt{s^2+z_2^2}} - \frac{z_1}{\sqrt{s^2+z_1^2}}\right)\boldsymbol{e}_\phi .
$$

線分の端を場の点から見込む角を、垂線からの角度 $\theta_i$（$\sin\theta_i = z_i/\sqrt{s^2+z_i^2}$）で測れば $B_\phi = \dfrac{\mu_0 I}{4\pi s}(\sin\theta_2 - \sin\theta_1)$ と書けます。

$z_2\to+\infty$ では $z_2/\sqrt{s^2+z_2^2}\to 1$、$z_1\to-\infty$ では $z_1/\sqrt{s^2+z_1^2}\to -1$ なので、括弧の中は $2$ に収束し

$$
\boldsymbol{B} \to \frac{\mu_0 I}{4\pi s}\cdot 2\,\boldsymbol{e}_\phi = \frac{\mu_0 I}{2\pi s}\boldsymbol{e}_\phi
$$

となって <Ref to="ex-straight-wire" /> と一致します。

なお、線分だけの電流には <Ref to="cor-steady-divj" /> が成り立たないので、<Ref to="cor-ampere-integral" /> をこの配置に直接適用してはいけません。この公式が意味をもつのは、閉回路全体を線分に分けて足し合わせるときだけです。
</Solution>
</Exercise>

<Exercise id="exr-coaxial" difficulty="標準">
同軸ケーブルを考える。中心の導体は半径 $a$ の円柱で、一様な電流密度で全電流 $I$ が $+z$ 方向に流れる。外側の導体は $b \le s \le c$（$a < b$）の円筒殻で、一様な電流密度で全電流 $I$ が $-z$ 方向に流れる。四つの領域 $s < a$、$a < s < b$、$b < s < c$、$s > c$ のそれぞれで磁場を求めよ。

<Solution>
<Ref to="ex-cylinder" /> と同じ対称性の議論から $\boldsymbol{B} = B_\phi(s)\boldsymbol{e}_\phi$ であり、半径 $s$ の円に <Ref to="cor-ampere-integral" /> を適用して $B_\phi(s) = \dfrac{\mu_0 I_{\text{enc}}(s)}{2\pi s}$ となります。あとは $I_{\text{enc}}$ を数えるだけです。

**$s < a$：** 内部導体の断面積比だけ数えて $I_{\text{enc}} = I\,s^2/a^2$。よって $B_\phi = \dfrac{\mu_0 I s}{2\pi a^2}$。

**$a < s < b$：** $I_{\text{enc}} = I$。よって $B_\phi = \dfrac{\mu_0 I}{2\pi s}$。

**$b < s < c$：** 外側導体の断面積は $\pi(c^2-b^2)$ で、半径 $s$ までに含まれる面積は $\pi(s^2-b^2)$ だから、逆向きの電流の寄与は $-I\dfrac{s^2-b^2}{c^2-b^2}$。合わせて

$$
I_{\text{enc}} = I\left(1 - \frac{s^2-b^2}{c^2-b^2}\right) = I\,\frac{c^2-s^2}{c^2-b^2},
\qquad
B_\phi = \frac{\mu_0 I}{2\pi s}\cdot\frac{c^2-s^2}{c^2-b^2}.
$$

**$s > c$：** $I_{\text{enc}} = I - I = 0$、よって $B_\phi = 0$。

外部で磁場がまったく漏れないことが同軸ケーブルの利点です。$s = b$ と $s = c$ で式が連続につながることも確かめてください（$s=b$ で第 3 式は $\mu_0 I/(2\pi b)$、$s=c$ で $0$ になります）。
</Solution>
</Exercise>

<Exercise id="exr-toroid" difficulty="易">
中心軸を $z$ 軸とするドーナツ形の芯に導線を $N$ 回一様に巻き、電流 $I$ を流す（トロイダルコイル）。巻線は芯を密に覆っているとする。芯の内部（コイルが囲む領域）と外部で磁場を求めよ。

<Solution>
系は $z$ 軸まわりの回転対称です。巻線の電流は $sz$ 平面内を流れるので、<Ref to="ex-solenoid" /> と同じ議論により $\boldsymbol{B} = B_\phi(s,z)\,\boldsymbol{e}_\phi$ の形になります（$\boldsymbol{e}_\phi$ 成分だけが残る）。

$z$ 軸を中心とする半径 $s$ の水平な円 $C$ を取り、<Ref to="cor-ampere-integral" /> を適用します。左辺は $B_\phi\cdot 2\pi s$ です。

- $C$ が芯の内部を通るとき、$C$ を境界とする円板は巻線を $N$ 回くぐるので $I_{\text{enc}} = NI$。よって

$$
B_\phi = \frac{\mu_0 N I}{2\pi s}.
$$

- $C$ が芯の穴の内側（$s$ が小さい）を通るとき、円板は巻線をくぐらないので $I_{\text{enc}} = 0$、よって $B_\phi = 0$。
- $C$ が芯より外側を通るとき、円板は各巻線を往きと還りで 2 回くぐり、向きが逆なので正味 $I_{\text{enc}} = 0$。よって $B_\phi = 0$。

つまり磁場は芯の内部に完全に閉じ込められ、大きさは中心軸からの距離 $s$ に反比例します。芯の太さが中心半径に比べて小さければ、$s$ をほぼ一定とみなせて、ソレノイドと同じ $\mu_0 n I$（$n = N/(2\pi s)$ は単位長さあたりの巻き数）になります。
</Solution>
</Exercise>

<Exercise id="exr-capacitor" difficulty="難">
平行平板コンデンサーを、導線を通して一定電流 $I$ で充電している。導線を一周する円 $C$ を固定し、$C$ を境界とする二つの曲面を考える。$S_1$ は導線を垂直に貫く平らな円板、$S_2$ は導線を避けて極板のあいだを通る袋状の曲面である。$\mu_0\int_S \boldsymbol{J}\cdot d\boldsymbol{S}$ を $S = S_1, S_2$ について計算し、<Ref to="cor-ampere-integral" /> のどの仮定が破れているかを述べよ。

<Solution>
$S_1$ は導線を貫くので $\int_{S_1}\boldsymbol{J}\cdot d\boldsymbol{S} = I$、したがって右辺は $\mu_0 I$ です。一方 $S_2$ は極板のあいだを通ります。極板間は真空（あるいは絶縁体）で電荷の流れがないので $\boldsymbol{J} = \boldsymbol{0}$、よって $\int_{S_2}\boldsymbol{J}\cdot d\boldsymbol{S} = 0$ となり、右辺は $0$ です。

左辺 $\oint_C\boldsymbol{B}\cdot d\boldsymbol{l}$ は曲線 $C$ だけで決まる量ですから、$\mu_0 I \ne 0$ と $0$ が同時に等しいことはあり得ません。矛盾です。

破れている仮定は**定常性**です。極板には電荷が時間とともに蓄積していくので、極板の位置で $\partial\rho/\partial t \ne 0$ であり、<Ref to="thm-continuity" /> から $\nabla\cdot\boldsymbol{J} = -\partial\rho/\partial t \ne 0$ となります。すると <Ref to="cor-steady-divj" /> が使えず、<Ref to="lem-div-a" />、<Ref to="thm-ampere-differential" />、<Ref to="cor-ampere-integral" /> の導出がすべて崩れます。<Ref to="rem-surface-independence" /> で見た「曲面の取り方によらない」という性質も、まさに $\nabla\cdot\boldsymbol{J}=0$ に依拠していました。

この破綻を救うのがマクスウェルの変位電流で、$\boldsymbol{J}$ を $\boldsymbol{J} + \epsilon_0\partial\boldsymbol{E}/\partial t$ に置き換えると、連続の方程式により $\nabla\cdot\!\left(\boldsymbol{J} + \epsilon_0\dfrac{\partial\boldsymbol{E}}{\partial t}\right) = -\dfrac{\partial\rho}{\partial t} + \dfrac{\partial\rho}{\partial t} = 0$ が恒等的に成り立ちます（$\nabla\cdot\boldsymbol{E} = \rho/\epsilon_0$ を使いました）。詳しくは [電磁誘導と変位電流](/physics/electromagnetism/induction-and-displacement-current) を参照してください。
</Solution>
</Exercise>

## 参考文献

- 砂川重信『理論電磁気学』第 3 版、紀伊國屋書店、1999 — 定常電流と静磁場の章。ベクトルポテンシャルを軸にした構成が本記事の流れに近い。
- D. J. Griffiths, *Introduction to Electrodynamics*, 4th ed., Cambridge University Press, 2017 — Chapter 5 "Magnetostatics"。ビオ・サバールの法則からアンペールの法則を導く道筋が丁寧に書かれている。
- J. D. Jackson, *Classical Electrodynamics*, 3rd ed., Wiley, 1999 — Chapter 5。ベクトル解析の恒等式は付録にまとまっている。
- 『ファインマン物理学 III 電磁気学』岩波書店 — 磁場の章。回転や循環の物理的な意味づけについての説明が読みやすい。
- 太田浩一『電磁気学の基礎 I』東京大学出版会、2012 — 歴史的経緯（エルステッド、ビオ・サバール、アンペール）の記述が詳しい。
- 杉浦光夫『解析入門 II』東京大学出版会、1985 — ガウスの発散定理・ストークスの定理の厳密な扱い。
- BIPM, *The International System of Units (SI Brochure)*, 9th ed. — 2019 年の再定義後の $\mu_0$ の位置づけ。[SI Brochure (BIPM)](https://www.bipm.org/en/publications/si-brochure)

## Appendix: 使ったベクトル解析の恒等式

**記法の準備。** 以下ではデカルト座標の成分を $x_1, x_2, x_3$ とし、$\partial_i := \partial/\partial x_i$ と書きます。同じ添字が二度現れたら $1$ から $3$ まで和を取る（アインシュタインの規約）ものとします。レビ・チビタ記号 $\epsilon_{ijk}$ は、$(i,j,k)$ が $(1,2,3)$ の偶置換なら $+1$、奇置換なら $-1$、それ以外は $0$ と定めます。外積と回転はそれぞれ

$$
(\boldsymbol{a}\times\boldsymbol{b})_i = \epsilon_{ijk}a_j b_k,\qquad (\nabla\times\boldsymbol{F})_i = \epsilon_{ijk}\partial_j F_k
$$

と書けます。よく使う関係式は

$$
\epsilon_{kij}\epsilon_{klm} = \delta_{il}\delta_{jm} - \delta_{im}\delta_{jl}
$$

です（$\delta$ はクロネッカーのデルタ）。

**二重の回転の展開。** <Ref to="thm-ampere-differential" /> の証明で使った恒等式を確かめます。$\boldsymbol{A}$ を $C^2$ 級として

$$
\begin{aligned}
\big[\nabla\times(\nabla\times\boldsymbol{A})\big]_i
&= \epsilon_{ijk}\partial_j(\nabla\times\boldsymbol{A})_k
= \epsilon_{ijk}\epsilon_{klm}\partial_j\partial_l A_m \\
&= \epsilon_{kij}\epsilon_{klm}\partial_j\partial_l A_m
= (\delta_{il}\delta_{jm}-\delta_{im}\delta_{jl})\partial_j\partial_l A_m \\
&= \partial_j\partial_i A_j - \partial_j\partial_j A_i
= \partial_i(\nabla\cdot\boldsymbol{A}) - \nabla^2 A_i .
\end{aligned}
$$

二行目で $\epsilon_{ijk} = \epsilon_{kij}$（巡回置換で不変）を使い、最後の行では $\boldsymbol{A}$ が $C^2$ 級であることからシュワルツの定理により $\partial_j\partial_i = \partial_i\partial_j$ と交換しました。よって

$$
\nabla\times(\nabla\times\boldsymbol{A}) = \nabla(\nabla\cdot\boldsymbol{A}) - \nabla^2\boldsymbol{A} .
$$

**外積の発散。** <Ref to="thm-div-b" /> は <Ref to="prop-vector-potential" /> 経由で示しましたが、ビオ・サバールの積分から直接示すこともできます。そのとき使うのが

$$
\nabla\cdot(\boldsymbol{a}\times\boldsymbol{b}) = \boldsymbol{b}\cdot(\nabla\times\boldsymbol{a}) - \boldsymbol{a}\cdot(\nabla\times\boldsymbol{b})
$$

です。成分計算は次のとおりです。

$$
\begin{aligned}
\nabla\cdot(\boldsymbol{a}\times\boldsymbol{b}) &= \partial_i(\epsilon_{ijk}a_j b_k) = \epsilon_{ijk}(\partial_i a_j)b_k + \epsilon_{ijk}a_j(\partial_i b_k) \\
&= b_k\,\epsilon_{kij}\partial_i a_j - a_j\,\epsilon_{jik}\partial_i b_k
= \boldsymbol{b}\cdot(\nabla\times\boldsymbol{a}) - \boldsymbol{a}\cdot(\nabla\times\boldsymbol{b}).
\end{aligned}
$$

ここで第 1 項は $\epsilon_{ijk} = \epsilon_{kij}$、第 2 項は $\epsilon_{ijk} = -\epsilon_{jik}$ を用いました。

これを $\boldsymbol{a} = \boldsymbol{J}(\boldsymbol{r}')$、$\boldsymbol{b} = \boldsymbol{R}/R^3$ に適用すると、$\boldsymbol{J}(\boldsymbol{r}')$ は $\boldsymbol{r}$ によらないので $\nabla\times\boldsymbol{a} = \boldsymbol{0}$、また $\boldsymbol{R}/R^3 = -\nabla(1/R)$ で勾配の回転は $0$ だから $\nabla\times\boldsymbol{b} = \boldsymbol{0}$ です。よって <Ref to="ax-biot-savart" /> の被積分関数の発散は各点で $0$ となり、$\nabla\cdot\boldsymbol{B} = 0$ が直接得られます。