磁石をコイルに近づけたり遠ざけたりすると、動かしているあいだだけコイルに電流が流れます。また、コンデンサを充電している回路では、極板のあいだに何も流れていないのに磁場が現れます。この記事は、この二つの不思議な現象を一つの流れで説明し、電気と磁気の法則を 4 本の方程式(マクスウェル方程式)にまとめるところまでを追います。読み終えると、「変化する磁場が電場をつくる」ことと「変化する電場が磁場をつくる」ことの両方がなぜ必要なのかを、自分の言葉で説明できるようになります。
磁束の時間変化は起電力を生みます。これがファラデーの電磁誘導の法則で、積分形は ∮ C E ⋅ d l = − d d t ∫ S B ⋅ d S \displaystyle \oint_C \boldsymbol{E}\cdot d\boldsymbol{l} = -\frac{d}{dt}\int_S \boldsymbol{B}\cdot d\boldsymbol{S} ∮ C E ⋅ d l = − d t d ∫ S B ⋅ d S です。ストークスの定理を使うと微分形 rot E = − ∂ B / ∂ t \operatorname{rot}\boldsymbol{E} = -\partial \boldsymbol{B}/\partial t rot E = − ∂ B / ∂ t になります。
この一本の式によって、静電場で成り立っていた rot E = 0 \operatorname{rot}\boldsymbol{E}=\boldsymbol{0} rot E = 0 は特殊な場合の話に格下げされます。時間変化する場では電位(スカラーポテンシャル)だけでは電場を記述できません。
回路が動く場合、起電力の一部は電場ではなくローレンツ力(磁場の中を動く電荷が受ける力)の v × B \boldsymbol{v}\times\boldsymbol{B} v × B 項から来ます。「磁束の変化率」という一つの量が両者をまとめて与えることが、フラックス則の内容です。
静磁場のアンペールの法則 rot B = μ 0 j \operatorname{rot}\boldsymbol{B}=\mu_0\boldsymbol{j} rot B = μ 0 j は、電荷が溜まる状況では電荷保存則と矛盾します。コンデンサの充電がその典型例です。
矛盾を消す最小の修正が変位電流 ε 0 ∂ E / ∂ t \varepsilon_0 \partial\boldsymbol{E}/\partial t ε 0 ∂ E / ∂ t の追加です。これでマクスウェル方程式が閉じ、そこから真空中の電磁波と c = 1 / ε 0 μ 0 c = 1/\sqrt{\varepsilon_0\mu_0} c = 1/ ε 0 μ 0 が出てきます。
この記事で使う事前知識
電場 : 電荷が周りの空間につくる、他の電荷に力を及ぼす場のこと。→ 電場
電流密度 : 単位面積あたりに流れる電流を、向きも含めて表したベクトル場。→ 電流密度
発散 : ベクトル場の各点での「湧き出しの密度」。→ 発散
回転 : ベクトル場の各点での「渦の強さ」を表すベクトル。→ 回転
ストークスの定理 : 曲面上での回転の面積分が、境界の閉曲線に沿った線積分に等しいという関係。→ 定常電流と静磁場
1820 年にエルステッドが「電流のそばに置いた磁針が振れる」ことを発見しました。電流は磁場をつくる。ではその逆、磁場は電流をつくるのか。この素朴な問いに答えるまでに 11 年かかりました。答えが遅れた理由は単純です。磁石をコイルのそばに置いておくだけでは、何も起こらない からです。
1831 年、ファラデーは磁石を「動かした」瞬間だけ検流計の針が振れることを見つけました。定常状態ではなく、変化 が原因だったのです。この発見は、それまでの電磁気学の構造を根本から変えました。静電場と静磁場は、静電場とガウスの法則 と 定常電流と静磁場 で見たとおり、互いに独立な二つの理論として並立していました。電磁誘導は、その二つを初めて結びつける橋になります。
もう一つ、ファラデーの実験には注目すべき点があります。磁石を止めてコイルを動かしても、コイルを止めて磁石を動かしても、同じ起電力 が観測されます。しかし当時の理論的な説明は二つの場合で全く別物でした。コイルが動く場合は導体中の電荷が受けるローレンツ力で説明され、磁石が動く場合は「渦電場」の発生で説明されます。原因の説明が違うのに結果が同じ。この不自然さこそが、1905 年のアインシュタインの論文「動いている物体の電気力学」の冒頭で問題にされたことでした。この記事の §5 では、その二つの説明が一つの式にまとまる様子を見ます。
さらに、電磁誘導だけでは話は終わりません。ファラデーの法則を加えた 4 本の法則を並べてみると、そのうちの 1 本、静磁場のアンペールの法則が電荷保存則と両立しない ことが分かります。この矛盾を解消するためにマクスウェルが導入したのが変位電流です。§6 でこれを扱います。歴史の順序としては、電磁誘導の発見(1831)から変位電流の導入(1861–1865)まで 30 年です。この記事はその 30 年分を追体験する章になります。
SI 単位系を使います。電場 E \boldsymbol{E} E (V/m)、磁束密度 B \boldsymbol{B} B (T)、電荷密度 ρ \rho ρ (C/m³)、電流密度 j \boldsymbol{j} j (A/m²)、真空の誘電率 ε 0 \varepsilon_0 ε 0 、真空の透磁率 μ 0 \mu_0 μ 0 とします。
前章までに得ている法則は次の 3 本です。
法則 積分形 微分形 ガウスの法則 ∮ ∂ V E ⋅ d S = 1 ε 0 ∫ V ρ d V \displaystyle\oint_{\partial V} \boldsymbol{E}\cdot d\boldsymbol{S} = \frac{1}{\varepsilon_0}\int_V \rho\, dV ∮ ∂ V E ⋅ d S = ε 0 1 ∫ V ρ d V div E = ρ / ε 0 \operatorname{div}\boldsymbol{E} = \rho/\varepsilon_0 div E = ρ / ε 0 磁場のガウスの法則 ∮ ∂ V B ⋅ d S = 0 \displaystyle\oint_{\partial V} \boldsymbol{B}\cdot d\boldsymbol{S} = 0 ∮ ∂ V B ⋅ d S = 0 div B = 0 \operatorname{div}\boldsymbol{B} = 0 div B = 0 アンペールの法則(静磁場) ∮ C B ⋅ d l = μ 0 ∫ S j ⋅ d S \displaystyle\oint_C \boldsymbol{B}\cdot d\boldsymbol{l} = \mu_0 \int_S \boldsymbol{j}\cdot d\boldsymbol{S} ∮ C B ⋅ d l = μ 0 ∫ S j ⋅ d S rot B = μ 0 j \operatorname{rot}\boldsymbol{B} = \mu_0\boldsymbol{j} rot B = μ 0 j
そして、電荷の保存を表す連続の式(定理 2.3)[定常電流と静磁場]
div j + ∂ ρ ∂ t = 0 \operatorname{div}\boldsymbol{j} + \frac{\partial \rho}{\partial t} = 0 div j + ∂ t ∂ ρ = 0
があります。これは実験法則というより、電荷が生成も消滅もしないという要請を局所的に書いたものです。ある領域から電流として流れ出た分だけ、その領域の電荷が減る、と読みます。
積分形と微分形を行き来するために、ベクトル解析の 2 つの定理を使います。V V V を有界な領域、S S S を区分的に滑らかな有向曲面、C = ∂ S C = \partial S C = ∂ S をその境界とします。
∫ V div F d V = ∮ ∂ V F ⋅ d S , ∫ S rot F ⋅ d S = ∮ C F ⋅ d l \int_V \operatorname{div}\boldsymbol{F}\, dV = \oint_{\partial V} \boldsymbol{F}\cdot d\boldsymbol{S}, \qquad
\int_S \operatorname{rot}\boldsymbol{F}\cdot d\boldsymbol{S} = \oint_{C} \boldsymbol{F}\cdot d\boldsymbol{l} ∫ V div F d V = ∮ ∂ V F ⋅ d S , ∫ S rot F ⋅ d S = ∮ C F ⋅ d l
前者がガウスの発散定理(定理 5.3)[静電場とガウスの法則] 、後者がストークスの定理(定理 5.3)[定常電流と静磁場] です。面積分の計算そのものについては 重積分と累次積分 を参照してください。
定義 3.1 (磁束 ) ↗
有向曲面 S S S (法線 n \boldsymbol{n} n )を貫く磁束 Φ B \Phi_B Φ B を
Φ B = ∫ S B ⋅ d S = ∫ S B ⋅ n d A \Phi_B = \int_S \boldsymbol{B}\cdot d\boldsymbol{S} = \int_S \boldsymbol{B}\cdot \boldsymbol{n}\, dA Φ B = ∫ S B ⋅ d S = ∫ S B ⋅ n d A で定義する。単位はウェーバ(Wb = = = T·m²)である。
磁束が曲面 S S S の取り方によらず境界 C C C だけで決まることは、div B = 0 \operatorname{div}\boldsymbol{B}=0 div B = 0 (定理 4.4[定常電流と静磁場] )の帰結です。同じ境界 C C C を持つ 2 つの曲面 S 1 , S 2 S_1, S_2 S 1 , S 2 を貼り合わせると閉曲面ができ、そこを通る正味の磁束が 0 0 0 なので、∫ S 1 B ⋅ d S = ∫ S 2 B ⋅ d S \int_{S_1}\boldsymbol{B}\cdot d\boldsymbol{S} = \int_{S_2}\boldsymbol{B}\cdot d\boldsymbol{S} ∫ S 1 B ⋅ d S = ∫ S 2 B ⋅ d S となります。したがって「閉曲線 C C C を貫く磁束」と言ってよいわけです。この事実は §6 でもう一度、今度は電流について問い直されます。
定義 3.2 (起電力 ) ↗
閉曲線 C C C に沿った起電力 (electromotive force, EMF)E \mathcal{E} E を、単位電荷を C C C に沿って一周させる間に非静電的な力がする仕事として定義する。C C C 上の各点が速度 v \boldsymbol{v} v で動いているとき、電荷 q q q が受ける力は F = q ( E + v × B ) \boldsymbol{F} = q(\boldsymbol{E} + \boldsymbol{v}\times\boldsymbol{B}) F = q ( E + v × B ) なので
E = ∮ C ( E + v × B ) ⋅ d l \mathcal{E} = \oint_C (\boldsymbol{E} + \boldsymbol{v}\times\boldsymbol{B})\cdot d\boldsymbol{l} E = ∮ C ( E + v × B ) ⋅ d l である。単位はボルト(V)。C C C が静止しているときは E = ∮ C E ⋅ d l \mathcal{E} = \oint_C \boldsymbol{E}\cdot d\boldsymbol{l} E = ∮ C E ⋅ d l となる。
名前に「力」と付いていますが次元は電圧です。19 世紀の用語がそのまま残ったもので、力ではありません。
公理 3.3 (ファラデーの電磁誘導の法則 ) ↗
任意の閉曲線 C C C (時間とともに変形・移動してよい)と、C C C を境界とする任意の有向曲面 S ( t ) S(t) S ( t ) について
E = − d Φ B d t , Φ B ( t ) = ∫ S ( t ) B ⋅ d S \mathcal{E} = -\frac{d\Phi_B}{dt}, \qquad \Phi_B(t) = \int_{S(t)} \boldsymbol{B}\cdot d\boldsymbol{S} E = − d t d Φ B , Φ B ( t ) = ∫ S ( t ) B ⋅ d S が成り立つ。ここで C C C の向きと S S S の法線は右ねじの規則で対応させる。
右辺のマイナス符号がレンツの法則 です。誘導起電力は、それが駆動する電流がつくる磁場によって、元の磁束変化を打ち消す向きに生じます。もしこの符号がプラスだったらどうなるか考えてみてください。磁束がわずかに増えると、それを増やす向きの電流が流れ、磁束がさらに増え、電流がさらに増える。エネルギーが無限に湧き出します。マイナス符号はエネルギー保存の要請そのものです。
例 3.4 (一様磁場中で回転するコイル(交流発電機) )
一様な磁束密度 B = B x ^ \boldsymbol{B} = B\hat{\boldsymbol{x}} B = B x ^ の中に、面積 A A A 、巻数 N N N の平面コイルを置き、z z z 軸のまわりに一定の角速度 ω \omega ω で回転させます。時刻 t t t でのコイルの法線を n ( t ) = ( cos ω t , sin ω t , 0 ) \boldsymbol{n}(t) = (\cos\omega t,\ \sin\omega t,\ 0) n ( t ) = ( cos ω t , sin ω t , 0 ) とすると、1 巻きあたりの磁束は
Φ 1 ( t ) = B ⋅ n ( t ) A = B A cos ω t \Phi_1(t) = \boldsymbol{B}\cdot\boldsymbol{n}(t)\, A = BA\cos\omega t Φ 1 ( t ) = B ⋅ n ( t ) A = B A cos ω t です。N N N 巻きのコイルは同じ磁束を N N N 回貫かせるので Φ B = N B A cos ω t \Phi_B = NBA\cos\omega t Φ B = N B A cos ω t 、したがって 公理 3.3 より
E = − d d t ( N B A cos ω t ) = N B A ω sin ω t \mathcal{E} = -\frac{d}{dt}\bigl(NBA\cos\omega t\bigr) = NBA\,\omega \sin\omega t E = − d t d ( N B A cos ω t ) = N B A ω sin ω t となります。振幅 N B A ω NBA\omega N B A ω 、角周波数 ω \omega ω の正弦波、つまり交流です。数値を入れてみます。N = 200 N = 200 N = 200 、A = 2.0 × 10 − 2 A = 2.0\times10^{-2} A = 2.0 × 1 0 − 2 m²、B = 0.50 B = 0.50 B = 0.50 T、f = 50 f = 50 f = 50 Hz(ω = 2 π f = 314 \omega = 2\pi f = 314 ω = 2 π f = 314 rad/s)とすると
N B A ω = 200 × 0.50 × 2.0 × 10 − 2 × 314 = 6.3 × 10 2 V NBA\omega = 200 \times 0.50 \times 2.0\times 10^{-2} \times 314 = 6.3\times 10^{2}\ \text{V} N B A ω = 200 × 0.50 × 2.0 × 1 0 − 2 × 314 = 6.3 × 1 0 2 V で、およそ 630 V の振幅になります。発電所の同期発電機はこの原理そのものです。
定理 4.1 (ファラデーの法則の微分形 ) ✓
場 E , B \boldsymbol{E}, \boldsymbol{B} E , B が C 1 C^1 C 1 級であるとする。静止した 任意の閉曲線 C C C とそれを境界とする任意の有向曲面 S S S について 公理 3.3 が成り立つことは、空間の各点で
rot E = − ∂ B ∂ t \operatorname{rot}\boldsymbol{E} = -\frac{\partial \boldsymbol{B}}{\partial t} rot E = − ∂ t ∂ B が成り立つことと同値である。
次の証明の中心は、「どんな曲面で積分しても 0 0 0 になるなら、各点で 0 0 0 である」という論法です。曲面は自由に選べるので、問題の点のまわりの小さな円板だけを選びます。もしその点で値が 0 0 0 でなければ、円板を十分小さくすれば円板の上でも同じ符号が続き、積分が 0 0 0 にならなくなって矛盾します。円板を使うのは、小さくすれば被積分関数の符号が変わらないと言える形にできるからです。
証明(定理 4.1) (積分形 ⇒ \Rightarrow ⇒ 微分形) C C C と S S S は静止しているので、時間微分と面積分の順序を交換できます。B \boldsymbol{B} B が C 1 C^1 C 1 級で S S S が有界であることから、微分と積分の交換が正当化されます。よって
− d d t ∫ S B ⋅ d S = − ∫ S ∂ B ∂ t ⋅ d S -\frac{d}{dt}\int_S \boldsymbol{B}\cdot d\boldsymbol{S} = -\int_S \frac{\partial\boldsymbol{B}}{\partial t}\cdot d\boldsymbol{S} − d t d ∫ S B ⋅ d S = − ∫ S ∂ t ∂ B ⋅ d S です。一方、C C C が静止しているので 定義 3.2 より E = ∮ C E ⋅ d l \mathcal{E} = \oint_C \boldsymbol{E}\cdot d\boldsymbol{l} E = ∮ C E ⋅ d l であり、ストークスの定理(§2)から
∮ C E ⋅ d l = ∫ S rot E ⋅ d S \oint_C \boldsymbol{E}\cdot d\boldsymbol{l} = \int_S \operatorname{rot}\boldsymbol{E}\cdot d\boldsymbol{S} ∮ C E ⋅ d l = ∫ S rot E ⋅ d S です。公理 3.3 を代入すると、任意の有向曲面 S S S について
∫ S ( rot E + ∂ B ∂ t ) ⋅ d S = 0 \int_S \left(\operatorname{rot}\boldsymbol{E} + \frac{\partial\boldsymbol{B}}{\partial t}\right)\cdot d\boldsymbol{S} = 0 ∫ S ( rot E + ∂ t ∂ B ) ⋅ d S = 0 が成り立ちます。ここで G : = rot E + ∂ B / ∂ t \boldsymbol{G} := \operatorname{rot}\boldsymbol{E} + \partial\boldsymbol{B}/\partial t G := rot E + ∂ B / ∂ t と置き、ある点 P P P で G ( P ) ≠ 0 \boldsymbol{G}(P) \ne \boldsymbol{0} G ( P ) = 0 と仮定して矛盾を導きます。n 0 : = G ( P ) / ∣ G ( P ) ∣ \boldsymbol{n}_0 := \boldsymbol{G}(P)/|\boldsymbol{G}(P)| n 0 := G ( P ) /∣ G ( P ) ∣ とおくと G ( P ) ⋅ n 0 = ∣ G ( P ) ∣ > 0 \boldsymbol{G}(P)\cdot\boldsymbol{n}_0 = |\boldsymbol{G}(P)| > 0 G ( P ) ⋅ n 0 = ∣ G ( P ) ∣ > 0 です。G \boldsymbol{G} G は連続なので、P P P を含む十分小さい半径 a a a の円板 D a D_a D a (法線 n 0 \boldsymbol{n}_0 n 0 、中心 P P P )を取れば、その上で G ⋅ n 0 > ∣ G ( P ) ∣ / 2 > 0 \boldsymbol{G}\cdot\boldsymbol{n}_0 > |\boldsymbol{G}(P)|/2 > 0 G ⋅ n 0 > ∣ G ( P ) ∣/2 > 0 とできます。すると
∫ D a G ⋅ d S > ∣ G ( P ) ∣ 2 π a 2 > 0 \int_{D_a} \boldsymbol{G}\cdot d\boldsymbol{S} > \frac{|\boldsymbol{G}(P)|}{2}\,\pi a^2 > 0 ∫ D a G ⋅ d S > 2 ∣ G ( P ) ∣ π a 2 > 0 となり、上の等式に反します。したがってすべての点で G = 0 \boldsymbol{G} = \boldsymbol{0} G = 0 、すなわち rot E = − ∂ B / ∂ t \operatorname{rot}\boldsymbol{E} = -\partial\boldsymbol{B}/\partial t rot E = − ∂ B / ∂ t です。
(微分形 ⇒ \Rightarrow ⇒ 積分形) 逆向きは上の計算をそのまま逆にたどります。任意の静止した S S S について両辺を面積分し、左辺にストークスの定理、右辺に微分と積分の交換を適用すれば ∮ C E ⋅ d l = − d Φ B / d t \oint_C\boldsymbol{E}\cdot d\boldsymbol{l} = -d\Phi_B/dt ∮ C E ⋅ d l = − d Φ B / d t を得ます。
∎
この「任意の曲面で積分が 0 0 0 なら被積分関数が 0 0 0 」という論法は電磁気学で何度も使うので、上のように一度きちんと書いておく価値があります。使った仮定は G \boldsymbol{G} G の連続性 だけです。連続性がなければ、たとえば一点でだけ値を変えた関数を考えると結論は成り立ちません。
次の注意は、静電場で使った電位が、磁場が時間変化する場では電場の記述に足りなくなる、という話です。電位は場所ごとに決まる一つの値ですが、電場が渦を巻くと、閉じた経路を一周して戻っても元の値に戻らなくなります。
例 4.3 (ソレノイド内外に誘導される電場 )
半径 R R R の無限に長いソレノイドがあり、内部に一様な磁束密度 B ( t ) = B ( t ) z ^ \boldsymbol{B}(t) = B(t)\hat{\boldsymbol{z}} B ( t ) = B ( t ) z ^ 、外部では B = 0 \boldsymbol{B}=\boldsymbol{0} B = 0 とします。軸対称性から、誘導電場は方位角方向の成分 E φ ( r , t ) E_\varphi(r,t) E φ ( r , t ) のみを持つと仮定できます。
半径 r r r の円 C r C_r C r (軸に垂直、中心は軸上、z ^ \hat{\boldsymbol{z}} z ^ を法線とする向き)に 公理 3.3 を適用します。C r C_r C r は静止しているので E = ∮ C r E ⋅ d l = 2 π r E φ \mathcal{E} = \oint_{C_r}\boldsymbol{E}\cdot d\boldsymbol{l} = 2\pi r E_\varphi E = ∮ C r E ⋅ d l = 2 π r E φ です。
r < R r < R r < R のとき、円板を貫く磁束は Φ B = π r 2 B ( t ) \Phi_B = \pi r^2 B(t) Φ B = π r 2 B ( t ) なので
2 π r E φ = − π r 2 d B d t ⟹ E φ = − r 2 d B d t 2\pi r E_\varphi = -\pi r^2 \frac{dB}{dt} \quad\Longrightarrow\quad E_\varphi = -\frac{r}{2}\frac{dB}{dt} 2 π r E φ = − π r 2 d t d B ⟹ E φ = − 2 r d t d B r > R r > R r > R のとき、磁束はソレノイド内部の分しかないので Φ B = π R 2 B ( t ) \Phi_B = \pi R^2 B(t) Φ B = π R 2 B ( t ) であり
2 π r E φ = − π R 2 d B d t ⟹ E φ = − R 2 2 r d B d t 2\pi r E_\varphi = -\pi R^2 \frac{dB}{dt} \quad\Longrightarrow\quad E_\varphi = -\frac{R^2}{2r}\frac{dB}{dt} 2 π r E φ = − π R 2 d t d B ⟹ E φ = − 2 r R 2 d t d B となります。注目すべきは後者です。ソレノイドの外では B = 0 \boldsymbol{B} = \boldsymbol{0} B = 0 なのに、電場は 0 0 0 ではありません。「磁場のあるところに電場が生じる」のではなく、「磁束の変化を囲む 経路に沿って電場が回る」のです。この誘導電場を使って電子を加速する装置がベータトロンです。
内側の円:ソレノイド(半径 R) 点:B は紙面から手前向き(磁束が貫く範囲) 破線の円 C_r:半径 r の経路(r が R より大きい) 矢印:dB/dt が正のときの誘導電場の向き ソレノイドを軸の方向から見た断面。内側の円は磁束が貫く範囲、破線の円 C_r は半径 r(> R)の経路。dB/dt > 0 のとき誘導電場は矢印のように時計回りに回る
定理 4.1 は静止した回路の話でした。回路が動く場合、E \mathcal{E} E には v × B \boldsymbol{v}\times\boldsymbol{B} v × B の寄与も入ります。公理 3.3 はこの場合も含めて成り立つと主張していますが、それが微分形と整合することを確かめる必要があります。鍵は次の補題です。
補題 5.1 (運動する曲面上の磁束の時間微分 ) ✓
C ( t ) C(t) C ( t ) を、各点が速度場 v ( r , t ) \boldsymbol{v}(\boldsymbol{r},t) v ( r , t ) に従って運動する閉曲線とし、S ( t ) S(t) S ( t ) をそれを境界とする有向曲面(同じ速度場で運動する)とする。B \boldsymbol{B} B が C 1 C^1 C 1 級で div B = 0 \operatorname{div}\boldsymbol{B}=0 div B = 0 が成り立つとき
d d t ∫ S ( t ) B ⋅ d S = ∫ S ( t ) ∂ B ∂ t ⋅ d S − ∮ C ( t ) ( v × B ) ⋅ d l \frac{d}{dt}\int_{S(t)} \boldsymbol{B}\cdot d\boldsymbol{S}
= \int_{S(t)} \frac{\partial\boldsymbol{B}}{\partial t}\cdot d\boldsymbol{S} - \oint_{C(t)} (\boldsymbol{v}\times\boldsymbol{B})\cdot d\boldsymbol{l} d t d ∫ S ( t ) B ⋅ d S = ∫ S ( t ) ∂ t ∂ B ⋅ d S − ∮ C ( t ) ( v × B ) ⋅ d l が成り立つ。
この補題は、動く回路を貫く磁束の変化が「磁場そのものの時間変化」と「回路が動いたことによる分」の二つに分けられる、という式です。証明では、時刻 t t t の曲面、時刻 t + d t t+dt t + d t の曲面、そして境界の曲線が d t dt d t のあいだに掃く細い帯の 3 つで一つの閉じた立体を作ります。閉じた面を通る磁束の合計は 0 0 0 なので、3 つの面の磁束がつりあう関係から、帯を通る分が「回路が動いたことによる分」として取り出せます。
証明(補題 5.1) 時刻 t t t から t + d t t+dt t + d t までの間に S S S が掃く領域 V V V を考えます。V V V の境界は 3 つの部分からなります。S ( t ) S(t) S ( t ) (外向き法線は − n -\boldsymbol{n} − n )、S ( t + d t ) S(t+dt) S ( t + d t ) (外向き法線は + n +\boldsymbol{n} + n )、そして境界 C C C が掃く帯 Σ \Sigma Σ です。Σ \Sigma Σ の外向き面積要素は d S Σ = d l × ( v d t ) d\boldsymbol{S}_\Sigma = d\boldsymbol{l}\times(\boldsymbol{v}\,dt) d S Σ = d l × ( v d t ) です(C C C の向きと n \boldsymbol{n} n の右ねじの関係から、この外積が外向きになります)。
時刻 t + d t t+dt t + d t の場 B \boldsymbol{B} B について div B = 0 \operatorname{div}\boldsymbol{B}=0 div B = 0 とガウスの発散定理を V V V に適用すると、閉曲面を通る正味の磁束が 0 0 0 なので
∫ S ( t + d t ) B ( t + d t ) ⋅ d S − ∫ S ( t ) B ( t + d t ) ⋅ d S + ∮ C B ( t + d t ) ⋅ ( d l × v d t ) = 0 \int_{S(t+dt)}\boldsymbol{B}(t+dt)\cdot d\boldsymbol{S} - \int_{S(t)}\boldsymbol{B}(t+dt)\cdot d\boldsymbol{S} + \oint_C \boldsymbol{B}(t+dt)\cdot\bigl(d\boldsymbol{l}\times\boldsymbol{v}\,dt\bigr) = 0 ∫ S ( t + d t ) B ( t + d t ) ⋅ d S − ∫ S ( t ) B ( t + d t ) ⋅ d S + ∮ C B ( t + d t ) ⋅ ( d l × v d t ) = 0 が成り立ちます。第 1 項は Φ B ( t + d t ) \Phi_B(t+dt) Φ B ( t + d t ) です。第 2 項は、B ( t + d t ) = B ( t ) + ( ∂ B / ∂ t ) d t + O ( d t 2 ) \boldsymbol{B}(t+dt) = \boldsymbol{B}(t) + (\partial\boldsymbol{B}/\partial t)\,dt + O(dt^2) B ( t + d t ) = B ( t ) + ( ∂ B / ∂ t ) d t + O ( d t 2 ) を代入して
∫ S ( t ) B ( t + d t ) ⋅ d S = Φ B ( t ) + d t ∫ S ( t ) ∂ B ∂ t ⋅ d S + O ( d t 2 ) \int_{S(t)}\boldsymbol{B}(t+dt)\cdot d\boldsymbol{S} = \Phi_B(t) + dt\int_{S(t)}\frac{\partial\boldsymbol{B}}{\partial t}\cdot d\boldsymbol{S} + O(dt^2) ∫ S ( t ) B ( t + d t ) ⋅ d S = Φ B ( t ) + d t ∫ S ( t ) ∂ t ∂ B ⋅ d S + O ( d t 2 ) です。第 3 項はすでに d t dt d t の 1 次なので B ( t + d t ) \boldsymbol{B}(t+dt) B ( t + d t ) を B ( t ) \boldsymbol{B}(t) B ( t ) で置き換えても誤差は O ( d t 2 ) O(dt^2) O ( d t 2 ) です。スカラー三重積の巡回性 B ⋅ ( d l × v ) = d l ⋅ ( v × B ) \boldsymbol{B}\cdot(d\boldsymbol{l}\times\boldsymbol{v}) = d\boldsymbol{l}\cdot(\boldsymbol{v}\times\boldsymbol{B}) B ⋅ ( d l × v ) = d l ⋅ ( v × B ) を使うと、第 3 項は d t ∮ C ( v × B ) ⋅ d l dt\oint_C(\boldsymbol{v}\times\boldsymbol{B})\cdot d\boldsymbol{l} d t ∮ C ( v × B ) ⋅ d l になります。以上をまとめて d t dt d t で割り、d t → 0 dt\to 0 d t → 0 とすれば主張の式を得ます。
∎
定理 5.2 (フラックス則 ) ✓
div B = 0 \operatorname{div}\boldsymbol{B}=0 div B = 0 と rot E = − ∂ B / ∂ t \operatorname{rot}\boldsymbol{E} = -\partial\boldsymbol{B}/\partial t rot E = − ∂ B / ∂ t が成り立つとする。このとき、速度場 v \boldsymbol{v} v で運動する任意の閉曲線 C ( t ) C(t) C ( t ) について
∮ C ( t ) ( E + v × B ) ⋅ d l = − d d t ∫ S ( t ) B ⋅ d S \oint_{C(t)} (\boldsymbol{E} + \boldsymbol{v}\times\boldsymbol{B})\cdot d\boldsymbol{l} = -\frac{d}{dt}\int_{S(t)}\boldsymbol{B}\cdot d\boldsymbol{S} ∮ C ( t ) ( E + v × B ) ⋅ d l = − d t d ∫ S ( t ) B ⋅ d S が成り立つ。すなわち 公理 3.3 は、静止回路の微分形 定理 4.1 とローレンツ力から導かれる。
証明(定理 5.2) 補題 5.1 より
− d Φ B d t = − ∫ S ( t ) ∂ B ∂ t ⋅ d S + ∮ C ( t ) ( v × B ) ⋅ d l -\frac{d\Phi_B}{dt} = -\int_{S(t)}\frac{\partial\boldsymbol{B}}{\partial t}\cdot d\boldsymbol{S} + \oint_{C(t)}(\boldsymbol{v}\times\boldsymbol{B})\cdot d\boldsymbol{l} − d t d Φ B = − ∫ S ( t ) ∂ t ∂ B ⋅ d S + ∮ C ( t ) ( v × B ) ⋅ d l です。ここで仮定 rot E = − ∂ B / ∂ t \operatorname{rot}\boldsymbol{E} = -\partial\boldsymbol{B}/\partial t rot E = − ∂ B / ∂ t を第 1 項に使うと
− ∫ S ( t ) ∂ B ∂ t ⋅ d S = ∫ S ( t ) rot E ⋅ d S = ∮ C ( t ) E ⋅ d l -\int_{S(t)}\frac{\partial\boldsymbol{B}}{\partial t}\cdot d\boldsymbol{S} = \int_{S(t)}\operatorname{rot}\boldsymbol{E}\cdot d\boldsymbol{S} = \oint_{C(t)}\boldsymbol{E}\cdot d\boldsymbol{l} − ∫ S ( t ) ∂ t ∂ B ⋅ d S = ∫ S ( t ) rot E ⋅ d S = ∮ C ( t ) E ⋅ d l となります(最後の等号はストークスの定理)。この 2 式を足し合わせれば主張が出ます。なお div B = 0 \operatorname{div}\boldsymbol{B}=0 div B = 0 は 補題 5.1 の証明で使われています。
∎
この定理が §1 で述べた「不自然さ」の正体です。磁石を動かすと ∂ B / ∂ t ≠ 0 \partial\boldsymbol{B}/\partial t \ne \boldsymbol{0} ∂ B / ∂ t = 0 で誘導電場が生じ、コイルを動かすと v × B \boldsymbol{v}\times\boldsymbol{B} v × B 項が効きます。物理的な起源は別ですが、両者の和がちょうど磁束の全微分になる ので、結果は同じ式で書けます。この「偶然」は、後に相対論において E \boldsymbol{E} E と B \boldsymbol{B} B が一つのテンソル の成分であることから説明されます。
例 5.3 (レール上を滑る導体棒 )
間隔 L L L の 2 本の平行なレール上を、導体棒が速度 v v v で滑っています。一様な磁束密度 B B B が回路面に垂直(紙面に入る向き)にかかっているとします。
棒の中の自由電子は棒とともに速度 v \boldsymbol{v} v で動くので、ローレンツ力 q v × B q\boldsymbol{v}\times\boldsymbol{B} q v × B を受けます。この力の大きさは単位電荷あたり v B vB v B で、向きは棒に沿った方向です。棒の長さ L L L にわたって積分すると、定義 3.2 より
E = ∫ 0 L v B d l = B L v \mathcal{E} = \int_0^L vB\, dl = BLv E = ∫ 0 L v B d l = B Lv です(レールの部分は動いていないので v × B \boldsymbol{v}\times\boldsymbol{B} v × B の寄与がなく、静電場の寄与も一周すると打ち消し合います)。
一方、フラックスで計算すると、回路が囲む面積は A ( t ) = L x ( t ) A(t) = L\,x(t) A ( t ) = L x ( t ) 、d x / d t = v dx/dt = v d x / d t = v なので Φ B = B L x ( t ) \Phi_B = BLx(t) Φ B = B Lx ( t ) 、したがって E = − d Φ B / d t = − B L v \mathcal{E} = -d\Phi_B/dt = -BLv E = − d Φ B / d t = − B Lv です。符号の違いは向きの取り方の違いで、大きさは一致します。定理 5.2 が保証するとおりです。
回路の全抵抗を R R R とすると電流は I = B L v / R I = BLv/R I = B Lv / R 、棒が受ける磁気力は F = B I L = B 2 L 2 v / R F = BIL = B^2L^2v/R F = B I L = B 2 L 2 v / R で、レンツの法則どおり運動を妨げる向きです。単位時間あたりに外力がする仕事は F v = B 2 L 2 v 2 / R Fv = B^2L^2v^2/R F v = B 2 L 2 v 2 / R で、これは回路で消費されるジュール熱 I 2 R = B 2 L 2 v 2 / R I^2R = B^2L^2v^2/R I 2 R = B 2 L 2 v 2 / R にちょうど等しくなります。エネルギーは帳尻が合っています。
例 6.1 (平行板コンデンサを充電する回路 )
電池と抵抗と平行板コンデンサを直列につないだ回路を考えます。スイッチを入れた直後、導線には電流 I ( t ) I(t) I ( t ) が流れ、コンデンサの極板に電荷が溜まっていきます。ただし極板の間は真空(あるいは絶縁体)なので、そこには電流は流れません。
導線を囲む閉曲線 C C C を取り、静磁場のアンペールの法則(系 6.3[定常電流と静磁場] )
∮ C B ⋅ d l = μ 0 ∫ S j ⋅ d S \oint_C \boldsymbol{B}\cdot d\boldsymbol{l} = \mu_0 \int_S \boldsymbol{j}\cdot d\boldsymbol{S} ∮ C B ⋅ d l = μ 0 ∫ S j ⋅ d S を適用します。C C C を境界とする曲面として、まず導線を垂直に横切る平らな円板 S 1 S_1 S 1 を取ると、右辺は μ 0 I ( t ) \mu_0 I(t) μ 0 I ( t ) です。次に、C C C を境界としたまま袋のようにふくらませて、極板の間 を通る曲面 S 2 S_2 S 2 を取ります。S 2 S_2 S 2 はどの導線も横切らないので、右辺は 0 0 0 です。
同じ左辺に対して、右辺が μ 0 I \mu_0 I μ 0 I と 0 0 0 の二通り出てきました。矛盾です。
同じ境界 C をもつ 2 つの曲面。S₁ は導線を横切り、S₂ は極板の間を通る
命題 6.2 (静磁場のアンペールの法則の限界 ) ✓
B \boldsymbol{B} B が C 2 C^2 C 2 級であるとする。微分形 rot B = μ 0 j \operatorname{rot}\boldsymbol{B} = \mu_0\boldsymbol{j} rot B = μ 0 j が空間の各点で成り立つならば、その点で ∂ ρ / ∂ t = 0 \partial\rho/\partial t = 0 ∂ ρ / ∂ t = 0 でなければならない。したがって、電荷密度が時間変化する状況ではこの法則は成立しえない。
証明(命題 6.2) 任意の C 2 C^2 C 2 級ベクトル場 F \boldsymbol{F} F について div ( rot F ) = 0 \operatorname{div}(\operatorname{rot}\boldsymbol{F}) = 0 div ( rot F ) = 0 が恒等的に成り立ちます。成分で確認すると
div ( rot F ) = ∂ x ( ∂ y F z − ∂ z F y ) + ∂ y ( ∂ z F x − ∂ x F z ) + ∂ z ( ∂ x F y − ∂ y F x ) \operatorname{div}(\operatorname{rot}\boldsymbol{F}) = \partial_x(\partial_y F_z - \partial_z F_y) + \partial_y(\partial_z F_x - \partial_x F_z) + \partial_z(\partial_x F_y - \partial_y F_x) div ( rot F ) = ∂ x ( ∂ y F z − ∂ z F y ) + ∂ y ( ∂ z F x − ∂ x F z ) + ∂ z ( ∂ x F y − ∂ y F x ) で、∂ x ∂ y F z \partial_x\partial_y F_z ∂ x ∂ y F z と − ∂ y ∂ x F z -\partial_y\partial_x F_z − ∂ y ∂ x F z のように 6 項が偏微分の交換可能性によって 2 つずつ打ち消し合い、0 0 0 になります。
そこで rot B = μ 0 j \operatorname{rot}\boldsymbol{B} = \mu_0\boldsymbol{j} rot B = μ 0 j の両辺の発散を取ると 0 = μ 0 div j 0 = \mu_0\operatorname{div}\boldsymbol{j} 0 = μ 0 div j 、すなわち div j = 0 \operatorname{div}\boldsymbol{j} = 0 div j = 0 です。これを連続の式 div j + ∂ ρ / ∂ t = 0 \operatorname{div}\boldsymbol{j} + \partial\rho/\partial t = 0 div j + ∂ ρ / ∂ t = 0 (§2)に代入すると ∂ ρ / ∂ t = 0 \partial\rho/\partial t = 0 ∂ ρ / ∂ t = 0 を得ます。
∎
つまりアンペールの法則は、静磁場という前提を外した瞬間に、電荷保存則と衝突します。例 6.1 の逆理は、この一般的な事実が具体的な回路として現れたものにすぎません。
両辺の div をとる div(rot B) = 0 は恒等式
補正項を足す rot B = μ0 (j + jd)
変位電流が要請される論理
定義 6.3 (変位電流密度 ) ↗
真空中で
j d : = ε 0 ∂ E ∂ t \boldsymbol{j}_{\mathrm{d}} := \varepsilon_0 \frac{\partial \boldsymbol{E}}{\partial t} j d := ε 0 ∂ t ∂ E を変位電流密度 という。次元は電流密度と同じ A/m² である。曲面 S S S を通る変位電流 I d I_{\mathrm{d}} I d は I d = ∫ S j d ⋅ d S I_{\mathrm{d}} = \int_S \boldsymbol{j}_{\mathrm{d}}\cdot d\boldsymbol{S} I d = ∫ S j d ⋅ d S で定義される。
「電流」と呼びますが、電荷が流れているわけではありません。マクスウェルは当時、真空を満たすエーテルの微小な変位(displacement)が生む電流としてこれを解釈しました。エーテルは否定されましたが、名前だけが残っています。実体は「電場の時間変化」です。
命題 6.4 (変位電流による修復 ) ✓
アンペールの法則を
rot B = μ 0 ( j + ε 0 ∂ E ∂ t ) \operatorname{rot}\boldsymbol{B} = \mu_0\left(\boldsymbol{j} + \varepsilon_0\frac{\partial\boldsymbol{E}}{\partial t}\right) rot B = μ 0 ( j + ε 0 ∂ t ∂ E ) と修正すれば、ガウスの法則 div E = ρ / ε 0 \operatorname{div}\boldsymbol{E} = \rho/\varepsilon_0 div E = ρ / ε 0 のもとで、この式は連続の式と矛盾しない。さらに、修正項を j d \boldsymbol{j}_{\mathrm{d}} j d という未知のベクトル場として「任意の電荷分布に対して整合する」ことだけを要求すると、j d \boldsymbol{j}_{\mathrm{d}} j d は ε 0 ∂ E / ∂ t \varepsilon_0\partial\boldsymbol{E}/\partial t ε 0 ∂ E / ∂ t と回転場(発散が 0 0 0 の場)の分だけしか違わない。
証明(命題 6.4) (整合性) 修正した式の両辺の発散を取ります。左辺は 命題 6.2 の証明で示した恒等式により 0 0 0 です。右辺は
μ 0 ( div j + ε 0 ∂ ∂ t div E ) = μ 0 ( div j + ε 0 ∂ ∂ t ρ ε 0 ) = μ 0 ( div j + ∂ ρ ∂ t ) \mu_0\left(\operatorname{div}\boldsymbol{j} + \varepsilon_0 \frac{\partial}{\partial t}\operatorname{div}\boldsymbol{E}\right)
= \mu_0\left(\operatorname{div}\boldsymbol{j} + \varepsilon_0\frac{\partial}{\partial t}\frac{\rho}{\varepsilon_0}\right)
= \mu_0\left(\operatorname{div}\boldsymbol{j} + \frac{\partial\rho}{\partial t}\right) μ 0 ( div j + ε 0 ∂ t ∂ div E ) = μ 0 ( div j + ε 0 ∂ t ∂ ε 0 ρ ) = μ 0 ( div j + ∂ t ∂ ρ ) となり(2 番目の等号でガウスの法則、および時間微分と空間微分の交換可能性を使いました)、連続の式によりこれは 0 0 0 です。両辺とも 0 0 0 なので矛盾はありません。
(ほぼ一意) 修正項を j d \boldsymbol{j}_{\mathrm{d}} j d と書くと、上と同じ計算から整合の条件は div j + div j d = 0 \operatorname{div}\boldsymbol{j} + \operatorname{div}\boldsymbol{j}_{\mathrm{d}} = 0 div j + div j d = 0 、すなわち連続の式と合わせて
div j d = ∂ ρ ∂ t = ε 0 ∂ ∂ t div E = div ( ε 0 ∂ E ∂ t ) \operatorname{div}\boldsymbol{j}_{\mathrm{d}} = \frac{\partial\rho}{\partial t} = \varepsilon_0\frac{\partial}{\partial t}\operatorname{div}\boldsymbol{E} = \operatorname{div}\left(\varepsilon_0\frac{\partial\boldsymbol{E}}{\partial t}\right) div j d = ∂ t ∂ ρ = ε 0 ∂ t ∂ div E = div ( ε 0 ∂ t ∂ E ) です。よって div ( j d − ε 0 ∂ E / ∂ t ) = 0 \operatorname{div}\bigl(\boldsymbol{j}_{\mathrm{d}} - \varepsilon_0\partial\boldsymbol{E}/\partial t\bigr) = 0 div ( j d − ε 0 ∂ E / ∂ t ) = 0 となり、差は発散が 0 0 0 の場に限られます。
∎
例 6.6 (コンデンサの極板間の磁場 )
半径 R R R の円形極板が間隔 d d d (d ≪ R d \ll R d ≪ R )で向かい合っている平行板コンデンサに、一定電流 I I I が流れ込んでいるとします。極板上の電荷を Q ( t ) Q(t) Q ( t ) 、d Q / d t = I dQ/dt = I d Q / d t = I とします。
極板間の電場は(縁の効果を無視すれば)一様で、大きさは面電荷密度 σ = Q / ( π R 2 ) \sigma = Q/(\pi R^2) σ = Q / ( π R 2 ) を使って E = σ / ε 0 = Q / ( ε 0 π R 2 ) E = \sigma/\varepsilon_0 = Q/(\varepsilon_0 \pi R^2) E = σ / ε 0 = Q / ( ε 0 π R 2 ) 、向きは正極板から負極板へ向かう軸方向です。よって変位電流密度の大きさは
j d = ε 0 ∂ E ∂ t = ε 0 ⋅ 1 ε 0 π R 2 d Q d t = I π R 2 j_{\mathrm{d}} = \varepsilon_0\frac{\partial E}{\partial t} = \varepsilon_0 \cdot \frac{1}{\varepsilon_0 \pi R^2}\frac{dQ}{dt} = \frac{I}{\pi R^2} j d = ε 0 ∂ t ∂ E = ε 0 ⋅ ε 0 π R 2 1 d t d Q = π R 2 I です。極板の全断面 π R 2 \pi R^2 π R 2 にわたって積分すると I d = I I_{\mathrm{d}} = I I d = I 。極板間を「流れる」変位電流は、導線を流れる伝導電流にちょうど等しい わけです。これで 例 6.1 の逆理は消えます。S 1 S_1 S 1 では伝導電流 I I I 、S 2 S_2 S 2 では変位電流 I I I で、どちらも同じ答えになります。
極板間、軸から距離 r r r の点での磁束密度を求めます。軸を中心とする半径 r r r の円 C r C_r C r に修正されたアンペールの法則(命題 6.4 )の積分形を適用すると、r < R r < R r < R では
B ( r ) ⋅ 2 π r = μ 0 ⋅ I π R 2 ⋅ π r 2 ⟹ B ( r ) = μ 0 I r 2 π R 2 B(r)\cdot 2\pi r = \mu_0 \cdot \frac{I}{\pi R^2}\cdot \pi r^2
\quad\Longrightarrow\quad
B(r) = \frac{\mu_0 I r}{2\pi R^2} B ( r ) ⋅ 2 π r = μ 0 ⋅ π R 2 I ⋅ π r 2 ⟹ B ( r ) = 2 π R 2 μ 0 I r r > R r > R r > R では変位電流の全量 I I I が囲まれるので
B ( r ) ⋅ 2 π r = μ 0 I ⟹ B ( r ) = μ 0 I 2 π r B(r)\cdot 2\pi r = \mu_0 I \quad\Longrightarrow\quad B(r) = \frac{\mu_0 I}{2\pi r} B ( r ) ⋅ 2 π r = μ 0 I ⟹ B ( r ) = 2 π r μ 0 I となります。前者は一様な電流が流れる太い導線の内部の磁場(例 6.5[定常電流と静磁場] )と、後者は細い導線がつくる磁場(例 3.3[定常電流と静磁場] )と、それぞれ同じ形です。極板間には電荷の流れが一切ないのに、あたかも電流が流れているかのような磁場ができます。
数値を入れます。R = 5.0 R = 5.0 R = 5.0 cm、I = 1.0 I = 1.0 I = 1.0 A、r = R r = R r = R のとき
B = μ 0 I 2 π R = 4 π × 10 − 7 × 1.0 2 π × 0.050 = 2 × 10 − 7 0.050 = 4.0 × 10 − 6 T B = \frac{\mu_0 I}{2\pi R} = \frac{4\pi\times10^{-7}\times 1.0}{2\pi\times 0.050} = \frac{2\times10^{-7}}{0.050} = 4.0\times10^{-6}\ \text{T} B = 2 π R μ 0 I = 2 π × 0.050 4 π × 1 0 − 7 × 1.0 = 0.050 2 × 1 0 − 7 = 4.0 × 1 0 − 6 T で、4.0 μ T 4.0\ \mu\text{T} 4.0 μ T です。地磁気(およそ 50 μ T 50\ \mu\text{T} 50 μ T )の 10 分の 1 程度なので、注意深い実験なら測定できます。実際、この磁場は測定されています。
例 6.7 (変位電流はいつ効くか(数値評価) )
導体の中では、伝導電流と変位電流のどちらが優勢でしょうか。オームの法則 j = σ E \boldsymbol{j} = \sigma\boldsymbol{E} j = σ E (σ \sigma σ は電気伝導率)が成り立ち、場が角周波数 ω \omega ω で振動しているとすると、大きさの比は
∣ ε 0 ∂ E / ∂ t ∣ ∣ σ E ∣ = ε 0 ω σ \frac{|\varepsilon_0 \partial\boldsymbol{E}/\partial t|}{|\sigma\boldsymbol{E}|} = \frac{\varepsilon_0\omega}{\sigma} ∣ σ E ∣ ∣ ε 0 ∂ E / ∂ t ∣ = σ ε 0 ω です。銅の場合 σ = 5.96 × 10 7 \sigma = 5.96\times10^7 σ = 5.96 × 1 0 7 S/m なので、f = 1.0 f = 1.0 f = 1.0 GHz(ω = 6.28 × 10 9 \omega = 6.28\times10^9 ω = 6.28 × 1 0 9 rad/s)でも
ε 0 ω σ = 8.85 × 10 − 12 × 6.28 × 10 9 5.96 × 10 7 = 5.56 × 10 − 2 5.96 × 10 7 ≈ 9.3 × 10 − 10 \frac{\varepsilon_0\omega}{\sigma} = \frac{8.85\times10^{-12}\times 6.28\times10^{9}}{5.96\times10^{7}} = \frac{5.56\times10^{-2}}{5.96\times10^{7}} \approx 9.3\times10^{-10} σ ε 0 ω = 5.96 × 1 0 7 8.85 × 1 0 − 12 × 6.28 × 1 0 9 = 5.96 × 1 0 7 5.56 × 1 0 − 2 ≈ 9.3 × 1 0 − 10 で、変位電流は伝導電流の 10 億分の 1 以下です。金属導体の中では、光の振動数に近づくまで変位電流は無視できます。
逆に、真空や良い絶縁体の中では σ = 0 \sigma = 0 σ = 0 なので、電流はすべて変位電流です。電磁波が真空中を伝わるとき、磁場をつくっているのは変位電流だけです。つまり変位電流項は、実験室の回路ではほとんど見えないのに、光の存在そのものを担っているという、極端な二面性を持っています。
これで 4 本が揃いました。
名称 微分形 積分形 ガウスの法則 div E = ρ ε 0 \operatorname{div}\boldsymbol{E} = \dfrac{\rho}{\varepsilon_0} div E = ε 0 ρ ∮ ∂ V E ⋅ d S = Q in ε 0 \displaystyle\oint_{\partial V}\boldsymbol{E}\cdot d\boldsymbol{S} = \frac{Q_{\text{in}}}{\varepsilon_0} ∮ ∂ V E ⋅ d S = ε 0 Q in 磁場のガウスの法則 div B = 0 \operatorname{div}\boldsymbol{B} = 0 div B = 0 ∮ ∂ V B ⋅ d S = 0 \displaystyle\oint_{\partial V}\boldsymbol{B}\cdot d\boldsymbol{S} = 0 ∮ ∂ V B ⋅ d S = 0 ファラデーの法則 rot E = − ∂ B ∂ t \operatorname{rot}\boldsymbol{E} = -\dfrac{\partial\boldsymbol{B}}{\partial t} rot E = − ∂ t ∂ B ∮ C E ⋅ d l = − d d t ∫ S B ⋅ d S \displaystyle\oint_C \boldsymbol{E}\cdot d\boldsymbol{l} = -\frac{d}{dt}\int_S\boldsymbol{B}\cdot d\boldsymbol{S} ∮ C E ⋅ d l = − d t d ∫ S B ⋅ d S アンペール・マクスウェルの法則 rot B = μ 0 j + ε 0 μ 0 ∂ E ∂ t \operatorname{rot}\boldsymbol{B} = \mu_0\boldsymbol{j} + \varepsilon_0\mu_0\dfrac{\partial\boldsymbol{E}}{\partial t} rot B = μ 0 j + ε 0 μ 0 ∂ t ∂ E ∮ C B ⋅ d l = μ 0 I in + ε 0 μ 0 d d t ∫ S E ⋅ d S \displaystyle\oint_C \boldsymbol{B}\cdot d\boldsymbol{l} = \mu_0 I_{\text{in}} + \varepsilon_0\mu_0\frac{d}{dt}\int_S\boldsymbol{E}\cdot d\boldsymbol{S} ∮ C B ⋅ d l = μ 0 I in + ε 0 μ 0 d t d ∫ S E ⋅ d S
これに、電荷が受ける力を与えるローレンツ力の式 F = q ( E + v × B ) \boldsymbol{F} = q(\boldsymbol{E} + \boldsymbol{v}\times\boldsymbol{B}) F = q ( E + v × B ) を加えると、古典電磁気学は完結します。
系 7.1 (電荷保存則はマクスウェル方程式の帰結 ) ✓
マクスウェル方程式のうちガウスの法則とアンペール・マクスウェルの法則が成り立つならば、連続の式 div j + ∂ ρ / ∂ t = 0 \operatorname{div}\boldsymbol{j} + \partial\rho/\partial t = 0 div j + ∂ ρ / ∂ t = 0 が自動的に成り立つ。
証明(系 7.1) アンペール・マクスウェルの法則の両辺の発散を取ります。左辺は恒等式 div ( rot B ) = 0 \operatorname{div}(\operatorname{rot}\boldsymbol{B}) = 0 div ( rot B ) = 0 (命題 6.2 の証明を参照)により 0 0 0 です。右辺は
μ 0 div j + ε 0 μ 0 ∂ ∂ t div E = μ 0 div j + ε 0 μ 0 ∂ ∂ t ( ρ ε 0 ) = μ 0 ( div j + ∂ ρ ∂ t ) \mu_0\operatorname{div}\boldsymbol{j} + \varepsilon_0\mu_0\frac{\partial}{\partial t}\operatorname{div}\boldsymbol{E}
= \mu_0\operatorname{div}\boldsymbol{j} + \varepsilon_0\mu_0\frac{\partial}{\partial t}\left(\frac{\rho}{\varepsilon_0}\right)
= \mu_0\left(\operatorname{div}\boldsymbol{j} + \frac{\partial\rho}{\partial t}\right) μ 0 div j + ε 0 μ 0 ∂ t ∂ div E = μ 0 div j + ε 0 μ 0 ∂ t ∂ ( ε 0 ρ ) = μ 0 ( div j + ∂ t ∂ ρ ) です(2 番目の等号でガウスの法則を使いました)。μ 0 ≠ 0 \mu_0 \ne 0 μ 0 = 0 なので div j + ∂ ρ / ∂ t = 0 \operatorname{div}\boldsymbol{j} + \partial\rho/\partial t = 0 div j + ∂ ρ / ∂ t = 0 を得ます。
∎
これは注目に値します。§6 では電荷保存則を外から課す制約 として使って変位電流を導きました。ところが完成した方程式系では、電荷保存則は方程式から導かれる定理 になっています。方程式系がそれだけ強く縛られている、ということです。
次章では、真空中(ρ = 0 \rho = 0 ρ = 0 , j = 0 \boldsymbol{j} = \boldsymbol{0} j = 0 )でこの 4 本から E \boldsymbol{E} E と B \boldsymbol{B} B が波動方程式を満たすことを示します。そこで現れる速度 1 / ε 0 μ 0 1/\sqrt{\varepsilon_0\mu_0} 1/ ε 0 μ 0 が光速に一致することが、電磁気学と光学を統一する結論になります。
演習 8.1 標準
半径 R = 5.0 R = 5.0 R = 5.0 cm の円形極板をもつ平行板コンデンサに、一定の電流 I = 2.0 I = 2.0 I = 2.0 A が流れ込んでいます。縁の効果は無視できるとして、(1) 極板間の変位電流密度の大きさ、(2) 軸から r = 2.0 r = 2.0 r = 2.0 cm の点における磁束密度の大きさ、を求めてください。
解答 (1) 例 6.6 と同じ計算です。極板間の電場は一様で E = Q / ( ε 0 π R 2 ) E = Q/(\varepsilon_0\pi R^2) E = Q / ( ε 0 π R 2 ) なので
j d = ε 0 ∂ E ∂ t = I π R 2 = 2.0 π × ( 0.050 ) 2 = 2.0 7.854 × 10 − 3 = 2.5 × 10 2 A/m 2 j_{\mathrm{d}} = \varepsilon_0\frac{\partial E}{\partial t} = \frac{I}{\pi R^2} = \frac{2.0}{\pi\times(0.050)^2} = \frac{2.0}{7.854\times10^{-3}} = 2.5\times10^{2}\ \text{A/m}^2 j d = ε 0 ∂ t ∂ E = π R 2 I = π × ( 0.050 ) 2 2.0 = 7.854 × 1 0 − 3 2.0 = 2.5 × 1 0 2 A/m 2 およそ 255 255 255 A/m² です。
(2) r = 2.0 cm < R = 5.0 cm r = 2.0\ \text{cm} < R = 5.0\ \text{cm} r = 2.0 cm < R = 5.0 cm なので、半径 r r r の円が囲む変位電流は I d ( r ) = j d ⋅ π r 2 = I r 2 / R 2 I_{\mathrm{d}}(r) = j_{\mathrm{d}}\cdot\pi r^2 = I r^2/R^2 I d ( r ) = j d ⋅ π r 2 = I r 2 / R 2 です。アンペール・マクスウェルの法則(命題 6.4 )から
B ⋅ 2 π r = μ 0 I r 2 R 2 ⟹ B = μ 0 I r 2 π R 2 = ( 4 π × 10 − 7 ) ( 2.0 ) ( 0.020 ) 2 π × ( 0.050 ) 2 B\cdot 2\pi r = \mu_0 \frac{I r^2}{R^2}
\quad\Longrightarrow\quad
B = \frac{\mu_0 I r}{2\pi R^2}
= \frac{(4\pi\times10^{-7})(2.0)(0.020)}{2\pi\times(0.050)^2} B ⋅ 2 π r = μ 0 R 2 I r 2 ⟹ B = 2 π R 2 μ 0 I r = 2 π × ( 0.050 ) 2 ( 4 π × 1 0 − 7 ) ( 2.0 ) ( 0.020 ) 分子は 4 π × 10 − 7 × 0.040 = 1.6 π × 10 − 8 4\pi\times10^{-7}\times 0.040 = 1.6\pi\times10^{-8} 4 π × 1 0 − 7 × 0.040 = 1.6 π × 1 0 − 8 、分母は 2 π × 2.5 × 10 − 3 = 5.0 π × 10 − 3 2\pi\times 2.5\times10^{-3} = 5.0\pi\times10^{-3} 2 π × 2.5 × 1 0 − 3 = 5.0 π × 1 0 − 3 なので
B = 1.6 × 10 − 8 5.0 × 10 − 3 = 3.2 × 10 − 6 T B = \frac{1.6\times10^{-8}}{5.0\times10^{-3}} = 3.2\times10^{-6}\ \text{T} B = 5.0 × 1 0 − 3 1.6 × 1 0 − 8 = 3.2 × 1 0 − 6 T すなわち 3.2 μ T 3.2\ \mu\text{T} 3.2 μ T です。
演習 8.2 標準
鉛直に立てた間隔 L L L の 2 本の導体レールの上を、質量 m m m の導体棒が摩擦なく滑り落ちます。水平方向に一様な磁束密度 B B B が回路面に垂直にかかっており、回路の全抵抗は R R R (一定)です。棒の自己インダクタンスは無視できるとします。棒の終端速度 v ∞ v_\infty v ∞ を求め、m = 10 m = 10 m = 10 g, L = 0.20 L = 0.20 L = 0.20 m, B = 0.50 B = 0.50 B = 0.50 T, R = 0.10 Ω R = 0.10\ \Omega R = 0.10 Ω , g = 9.8 g = 9.8 g = 9.8 m/s² のときの値を計算してください。
解答 速度 v v v で落下しているとき、例 5.3 と同じ計算で起電力の大きさは E = B L v \mathcal{E} = BLv E = B Lv 、電流は I = B L v / R I = BLv/R I = B Lv / R です。この電流が流れる棒は磁場から力 F = B I L = B 2 L 2 v / R F = BIL = B^2L^2v/R F = B I L = B 2 L 2 v / R を受け、レンツの法則(公理 3.3 のマイナス符号)より向きは運動を妨げる向き、つまり上向きです。
運動方程式は
m d v d t = m g − B 2 L 2 R v m\frac{dv}{dt} = mg - \frac{B^2L^2}{R}v m d t d v = m g − R B 2 L 2 v です。終端速度は d v / d t = 0 dv/dt = 0 d v / d t = 0 となる速度なので
v ∞ = m g R B 2 L 2 v_\infty = \frac{mgR}{B^2L^2} v ∞ = B 2 L 2 m g R 数値を入れると
v ∞ = ( 0.010 ) ( 9.8 ) ( 0.10 ) ( 0.50 ) 2 ( 0.20 ) 2 = 9.8 × 10 − 3 0.25 × 0.040 = 9.8 × 10 − 3 1.0 × 10 − 2 = 0.98 m/s v_\infty = \frac{(0.010)(9.8)(0.10)}{(0.50)^2(0.20)^2} = \frac{9.8\times10^{-3}}{0.25\times 0.040} = \frac{9.8\times10^{-3}}{1.0\times10^{-2}} = 0.98\ \text{m/s} v ∞ = ( 0.50 ) 2 ( 0.20 ) 2 ( 0.010 ) ( 9.8 ) ( 0.10 ) = 0.25 × 0.040 9.8 × 1 0 − 3 = 1.0 × 1 0 − 2 9.8 × 1 0 − 3 = 0.98 m/s 約 0.98 0.98 0.98 m/s です。なお上の 1 階線形微分方程式を解くと v ( t ) = v ∞ ( 1 − e − t / τ ) v(t) = v_\infty\bigl(1 - e^{-t/\tau}\bigr) v ( t ) = v ∞ ( 1 − e − t / τ ) 、τ = m R / ( B 2 L 2 ) = 0.10 \tau = mR/(B^2L^2) = 0.10 τ = m R / ( B 2 L 2 ) = 0.10 s で、v ∞ v_\infty v ∞ に指数関数的に近づきます。
演習 8.3 標準
半径 R R R の無限に長いソレノイドが単位長さあたり n n n 回巻かれ、電流 I ( t ) = I 0 sin ω t I(t) = I_0\sin\omega t I ( t ) = I 0 sin ω t が流れています。内部の磁束密度は B ( t ) = μ 0 n I ( t ) B(t) = \mu_0 n I(t) B ( t ) = μ 0 n I ( t ) (軸方向、一様)、外部では 0 0 0 です。軸から距離 r > R r > R r > R の点における誘導電場の方位角成分 E φ ( r , t ) E_\varphi(r,t) E φ ( r , t ) を求め、その振幅が r r r とともにどう変化するかを述べてください。
解答 例 4.3 の r > R r > R r > R の場合をそのまま使います。軸を中心とする半径 r r r の円 C r C_r C r に 公理 3.3 を適用すると、C r C_r C r を貫く磁束はソレノイド内部の分だけなので Φ B = π R 2 B ( t ) \Phi_B = \pi R^2 B(t) Φ B = π R 2 B ( t ) であり
2 π r E φ = − d d t ( π R 2 μ 0 n I 0 sin ω t ) = − π R 2 μ 0 n I 0 ω cos ω t 2\pi r E_\varphi = -\frac{d}{dt}\bigl(\pi R^2 \mu_0 n I_0\sin\omega t\bigr) = -\pi R^2\mu_0 n I_0\omega\cos\omega t 2 π r E φ = − d t d ( π R 2 μ 0 n I 0 sin ω t ) = − π R 2 μ 0 n I 0 ω cos ω t したがって
E φ ( r , t ) = − μ 0 n R 2 I 0 ω 2 r cos ω t E_\varphi(r,t) = -\frac{\mu_0 n R^2 I_0 \omega}{2r}\cos\omega t E φ ( r , t ) = − 2 r μ 0 n R 2 I 0 ω cos ω t です。振幅は r r r に反比例して減衰します。B = 0 \boldsymbol{B} = \boldsymbol{0} B = 0 の領域なのに電場が生じており、しかもその大きさは囲まれた磁束の変化率だけで決まっています。これは電場が「その場の磁場」ではなく「囲んだ磁束の変化」に応答することの表れです。
(補足:ここでは準定常近似、すなわち ω r / c ≪ 1 \omega r/c \ll 1 ω r / c ≪ 1 を仮定しています。r r r が大きくなると変位電流の効果と遅延が効き、上の式は補正を受けます。)
演習 8.4 難
アンペール・マクスウェルの法則の積分形に現れる量
∫ S ( j + ε 0 ∂ E ∂ t ) ⋅ d S \int_S \left(\boldsymbol{j} + \varepsilon_0\frac{\partial\boldsymbol{E}}{\partial t}\right)\cdot d\boldsymbol{S} ∫ S ( j + ε 0 ∂ t ∂ E ) ⋅ d S が、境界 C C C を共有する任意の 2 つの曲面 S 1 , S 2 S_1, S_2 S 1 , S 2 について同じ値をとることを、ガウスの法則と連続の式から示してください(S 1 , S 2 S_1, S_2 S 1 , S 2 は互いに交わらず、貼り合わせて閉曲面 ∂ V \partial V ∂ V をつくるとします)。
解答 K : = j + ε 0 ∂ E / ∂ t \boldsymbol{K} := \boldsymbol{j} + \varepsilon_0\partial\boldsymbol{E}/\partial t K := j + ε 0 ∂ E / ∂ t と置きます。まず div K = 0 \operatorname{div}\boldsymbol{K} = 0 div K = 0 を示します。ガウスの法則 div E = ρ / ε 0 \operatorname{div}\boldsymbol{E} = \rho/\varepsilon_0 div E = ρ / ε 0 と、時間微分と空間微分の交換可能性から
div K = div j + ε 0 ∂ ∂ t div E = div j + ∂ ρ ∂ t = 0 \operatorname{div}\boldsymbol{K} = \operatorname{div}\boldsymbol{j} + \varepsilon_0\frac{\partial}{\partial t}\operatorname{div}\boldsymbol{E}
= \operatorname{div}\boldsymbol{j} + \frac{\partial\rho}{\partial t} = 0 div K = div j + ε 0 ∂ t ∂ div E = div j + ∂ t ∂ ρ = 0 (最後の等号は連続の式)です。
次に S 1 S_1 S 1 と S 2 S_2 S 2 を貼り合わせて閉曲面 ∂ V \partial V ∂ V をつくります。両者は同じ境界 C C C を持ち、C C C 上で同じ向きづけを共有しているので、閉曲面 ∂ V \partial V ∂ V の外向き法線で見ると、一方(S 2 S_2 S 2 とします)は元の法線と同じ向き、他方(S 1 S_1 S 1 )は逆向きになります。したがってガウスの発散定理より
∫ S 2 K ⋅ d S − ∫ S 1 K ⋅ d S = ∮ ∂ V K ⋅ d S = ∫ V div K d V = 0 \int_{S_2}\boldsymbol{K}\cdot d\boldsymbol{S} - \int_{S_1}\boldsymbol{K}\cdot d\boldsymbol{S} = \oint_{\partial V}\boldsymbol{K}\cdot d\boldsymbol{S} = \int_V \operatorname{div}\boldsymbol{K}\, dV = 0 ∫ S 2 K ⋅ d S − ∫ S 1 K ⋅ d S = ∮ ∂ V K ⋅ d S = ∫ V div K d V = 0 となり、2 つの面積分は等しくなります。
これが 例 6.1 の逆理が解消される理由の一般形です。逆理が生じたのは K \boldsymbol{K} K ではなく j \boldsymbol{j} j だけを使ったからで、div j ≠ 0 \operatorname{div}\boldsymbol{j}\ne 0 div j = 0 のときには j \boldsymbol{j} j の面積分は曲面の取り方に依存してしまいます。
R. P. Feynman, R. B. Leighton, M. Sands, The Feynman Lectures on Physics , Vol. II — Chapter 17 “The Laws of Induction”, Chapter 18 “The Maxwell Equations”. 全文が 公式サイト で公開されています。フラックス則が使えない例(ファラデーの円盤など)の議論は第 17 章にあります。
D. J. Griffiths, Introduction to Electrodynamics , 4th ed., Cambridge University Press, 2017 — Chapter 7 “Electrodynamics”。電磁誘導と変位電流の標準的な扱い。
J. D. Jackson, Classical Electrodynamics , 3rd ed., Wiley, 1999 — Chapter 5, Chapter 6。より進んだ扱いと、準定常近似の適用範囲。
砂川重信『理論電磁気学 第3版』紀伊國屋書店、1999 — 日本語の標準的な教科書。積分形と微分形の対応を丁寧に扱っています。
J. C. Maxwell, “A Dynamical Theory of the Electromagnetic Field”, Philosophical Transactions of the Royal Society of London 155 (1865), 459–512. 変位電流と電磁波の予言が提示された原論文。
A. Einstein, “Zur Elektrodynamik bewegter Körper”, Annalen der Physik 17 (1905), 891–921. 冒頭で「磁石と導体の相対運動」の非対称性が問題提起されています。
フラックス則は「回路」の概念に依存しています。 定理 5.2 は、速度場 v \boldsymbol{v} v で連続的に運動する閉曲線 C ( t ) C(t) C ( t ) について証明されました。証明で本質的だったのは、時刻 t t t の曲線と時刻 t + d t t+dt t + d t の曲線が、掃かれる帯 Σ \Sigma Σ によって滑らかにつながることです。この前提が崩れる状況では、− d Φ B / d t -d\Phi_B/dt − d Φ B / d t という計算が起電力を与えるとは限りません。
典型例はファラデーの円盤(単極誘導)です。 一様な磁場 B B B の中で、導体の円板を軸のまわりに角速度 ω \omega ω で回転させ、軸と円板の縁にブラシで導線を接触させて回路をつくります。この装置には確かに起電力が生じます。半径 r r r の位置にある電子は速度 v = ω r v = \omega r v = ω r で動き、単位電荷あたり ω r B \omega r B ω r B の径方向の力を受けるので、定義 3.2 に従って軸から縁(半径 a a a )まで積分すると
E = ∫ 0 a ω r B d r = 1 2 ω B a 2 \mathcal{E} = \int_0^a \omega r B\, dr = \frac{1}{2}\omega B a^2 E = ∫ 0 a ω r B d r = 2 1 ω B a 2
となり、これが実測値と一致します。ところがこの回路で「磁束」を計算しようとすると困ります。回路の形(軸 → 円板 → ブラシ → 外部導線 → 軸)は時間とともに変化しておらず、磁場も一定なので、素朴に描いた回路が囲む磁束は変化しません。それでも起電力は出ます。
この矛盾は、円板の中で「どの経路が回路なのか」が一意に決まらないことから来ています。 円板は連続体なので、電流の経路を特定の閉曲線として指定できません。定理 5.2 の前提である「物質とともに動く閉曲線」が定義できないのです。
したがって、より基本的なのは微分形とローレンツ力です。 rot E = − ∂ B / ∂ t \operatorname{rot}\boldsymbol{E} = -\partial\boldsymbol{B}/\partial t rot E = − ∂ B / ∂ t と F = q ( E + v × B ) \boldsymbol{F} = q(\boldsymbol{E}+\boldsymbol{v}\times\boldsymbol{B}) F = q ( E + v × B ) の 2 つは、回路の概念を一切使わずに場と粒子の言葉だけで書かれており、上の円盤にもそのまま適用できます。実際、上で行った計算は v × B \boldsymbol{v}\times\boldsymbol{B} v × B を積分しただけで、磁束は一度も使っていません。フラックス則は多くの場合に便利な導かれた規則 であって、法則の本体ではない、と理解しておいてください。