文献あり

電磁波の存在証明と光速の導出

11218

はじめに

今回の電磁波が存在することの証明は、私が今まで電磁気学を学んできた中で、一番感動した証明です。なんと、電磁気学の実験結果をまとめた綺麗な4本の式(マクスウェル方程式)から、電磁波が存在することが示せるんです！これは、本当に圧巻なのでみなさんにも味わって頂きたいと思い、この記事を作成しています笑。

前提知識

マクスウェル方程式

最初にマクスウェル方程式を確認していきましょう。

マクスウェル方程式は以下の4つの式のことを指す。
$\begin{array}{r} {\begin{cases} div E = 0 \\ div B = 0 \\ rot E + \frac{\partial B}{\partial t} = 0 \\ rot B - ε_{0} μ_{0} \frac{\partial E}{\partial t} = 0 \end{cases} \end{array}$
※ただし、電流、電荷が存在しない真空を仮定している。

これらの法則は今回既知としますが、上から順番にガウスの法則、単極磁荷の非存在性、ファラデーの法則、アンペール・マクスウェルの法則を表す式です。これらのイメージについては、また別の記事で解説しようと思います。

波動方程式

波動方程式は以下のように表される。
$\begin{array}{r} \frac{1}{v^{2}} \frac{\partial^{2} u}{\partial t^{2}} = \frac{\partial^{2} u}{\partial x^{2}} + \frac{\partial^{2} u}{\partial y^{2}} + \frac{\partial^{2} u}{\partial z^{2}} = Δ u \end{array}$

波動方程式は、このような式を満たす $u$ や $v$ が存在した時に、 $u$ が波の関数、 $v$ が波の速度とするような波が発生している考えることができるという式です。

試しに $u (x, t) = \sin (x + v t)$ (1次元の波)は波動方程式を満たす。実際、
$\frac{\partial^{2} u}{\partial t^{2}} = - v^{2} \sin (x + v t), \frac{\partial^{2} u}{\partial x^{2}} = - \sin (x + v t)$
ですから、波動方程式の左辺と右辺を比べれば明らかです。

「波以外の関数も波動関数を満たしてしまうのでは？」と思う方向けにもう少し深めておきましょう。1次元に限った場合は波動関数一般解は
$u = f (x - v t) + g (x + v t)$
と表されます。これは、順方向に進む波と逆方向に進む波の重ね合わせであると考えられるので、波動関数は妥当であると捉えることができます。

ベクトル解析

電磁波の存在を証明するためには、ベクトル解析の知識が必要不可欠であるから、今回使うものを挙げておく。

ベクトル解析の演算子に対して、以下の公式が成り立つ。
$rot rot A = grad div A - Δ A$

基本的にただ展開するだけになりますが、少々お付き合いください。この式は直感的にはイメージしづらいですね。証明を確認していきましょう。

全成分同じ計算となるので、両辺の $x$ 成分のみ比較します。また以下では、 $A = (A_{x}, A_{y}, A_{z})$ とする。
$\begin{array}{rcl} rot A & = & (\frac{\partial A_{z}}{\partial y} - \frac{\partial A_{y}}{\partial z}, \frac{\partial A_{x}}{\partial z} - \frac{\partial A_{z}}{\partial x}, \frac{\partial A_{y}}{\partial x} - \frac{\partial A_{x}}{\partial y}) \\ div A & = & \frac{\partial A_{x}}{\partial x} + \frac{\partial A_{y}}{\partial y} + \frac{\partial A_{z}}{\partial z} \end{array}$
より、左辺と右辺の $x$ 成分はそれぞれ、
$\begin{array}{rcl} (左辺の x 成分) & = & \frac{\partial}{\partial y} (\frac{\partial A_{y}}{\partial x} - \frac{\partial A_{x}}{\partial y}) - \frac{\partial}{\partial z} (\frac{\partial A_{x}}{\partial z} - \frac{\partial A_{z}}{\partial x}) \\ = & \frac{\partial^{2} A_{z}}{\partial z \partial x} + \frac{\partial^{2} A_{y}}{\partial y \partial x} - \frac{\partial^{2} A_{x}}{\partial y^{2}} - \frac{\partial^{2} A_{x}}{\partial z^{2}} \\ (右辺の x 成分) & = & \frac{\partial}{\partial x} (\frac{\partial A_{x}}{\partial x} + \frac{\partial A_{y}}{\partial y} + \frac{\partial A_{z}}{\partial z}) - (\frac{\partial^{2} A_{x}}{\partial x^{2}} + \frac{\partial^{2} A_{y}}{\partial y^{2}} + \frac{\partial^{2} A_{z}}{\partial z^{2}}) \\ = & \frac{\partial^{2} A_{z}}{\partial z \partial x} + \frac{\partial^{2} A_{y}}{\partial y \partial x} - \frac{\partial^{2} A_{x}}{\partial y^{2}} - \frac{\partial^{2} A_{x}}{\partial z^{2}} \end{array}$
となり一致する。

電磁波の存在性

いよいよ、前提知識が揃ったので電磁波の存在性を証明していきたいと思います。

電場

まずは、電場が波であることを示していきます。ファラデーの法則において両辺の $rot$ をとります。
$rot rot E + \frac{\partial}{\partial t} rot B = 0$
ここでベクトル解析の補題とアンペール・マクスウェルの法則の時間微分を用いることにより、
$grad div E - Δ E + ε_{0} μ_{0} \frac{\partial^{2} E}{\partial t^{2}} = 0$
となりますが、 $div E = 0$ ですから、 $c = \frac{1}{\sqrt{ε_{0} μ_{0}}}$ とすれば、
$\frac{1}{c^{2}} \frac{\partial^{2} E}{\partial t^{2}} = Δ E$
が導かれます。

磁場

次に、磁場が波であることを示していきます。アンペール・マクスウェルの法則において両辺の $rot$ をとります。
$rot rot B - ε_{0} μ_{0} \frac{\partial}{\partial t} rot E = 0$
こちらも同様にベクトル解析の補題とファラデーの法則の時間微分を用いることにより、

$grad div B - Δ B + ε_{0} μ_{0} \frac{\partial^{2} B}{\partial t^{2}} = 0$
となりますが、 $div B = 0$ より、
$Δ B = ε_{0} μ_{0} \frac{\partial^{2} B}{\partial t^{2}} = \frac{1}{c^{2}} \frac{\partial^{2} B}{\partial t^{2}}$
が導かれます。

電磁波と光速

電場と磁場に関して導かれた
$\begin{array}{r} {\begin{cases} \frac{1}{c^{2}} \frac{\partial^{2} E}{\partial t^{2}} = Δ E \\ \frac{1}{c^{2}} \frac{\partial^{2} B}{\partial t^{2}} = Δ B \end{cases} \end{array}$
をよく見てください。これは波動方程式の形になっています。つまり、これが電磁波の正体です。

上述した波動方程式と比較すると $v = c$ ということになりますが、 $v$ はいくつでしょうか。真空の誘電率 $ε_{0}$ と真空の透磁率は $μ_{0}$ はそれぞれ、 $ε_{0} ≃ 8.854 \times 10^{- 12}, μ_{0} ≃ 1.257 \times 10^{- 6}$ です(単位省略)。さて早速代入していきましょう。
$c = \frac{1}{\sqrt{ε_{0} μ_{0}}} ≃ 2.998 \times 10^{8}$
みなさんいかがですか。最高ですよね、光速とものすごい近い値が導かれました。