# 電磁波と光学：媒質中の波、フレネルの式、導波管、干渉と回折

> 媒質中の電磁波の速さ c/n と平面波の関係を導き、界面の境界条件から反射・屈折の法則とフレネルの式を証明する。矩形導波管のモードと遮断周波数、ヘルムホルツ・キルヒホッフの積分公式から単スリットと二重スリットの回折像までを示す。
> https://rikai.mugen-giken.com/physics/electromagnetism/waves-and-optics

## 0. この記事の要点

ガラスに入った光はなぜ曲がるのか、水面の照り返しはなぜ偏光サングラスで消えるのか、光ファイバーはなぜ曲げても光が漏れないのか。この記事は、こうした身近な現象の土台にある反射・屈折・全反射を、マクスウェル方程式と物質中の境界条件から説明します（回折では、これに開口についての近似を加えます）。読み終えると、反射と屈折の法則、金属管の中を電磁波が通る条件、スリットを通った光がつくる縞の形を、式から導けるようになります。

- 一様で線形・等方な媒質の中でも、電磁場は波動方程式を満たします。速さは真空の $c$ ではなく $v=1/\sqrt{\varepsilon\mu}=c/n$ で、屈折率は $n=\sqrt{\varepsilon_r\mu_r}$ です。導電性がある場合は波数が複素数になり、波は深さ $\delta=\sqrt{2/(\mu\sigma\omega)}$ 程度で減衰します（表皮効果）。
- 2 つの透明な媒質の界面での反射と屈折は、[物質中の電磁場](/physics/electromagnetism/fields-in-matter)で証明した境界条件から導けます。入射波・反射波・透過波の界面方向の位相がそろう条件から**反射の法則**と**スネルの法則**が、振幅の接続条件から**フレネルの式**が出ます。
- フレネルの式の帰結として、$p$ 偏光の反射がちょうど消える**ブリュースター角**（$\tan\theta_B=n_2/n_1$）と、屈折率の大きい側から入るときの**全反射**が説明されます。全反射では反射率が 1 でも、界面の向こう側に減衰する波（エバネッセント波）がしみ出します。
- 完全導体の壁で囲まれた導波管では、横方向の固有値問題が 1 つ解けるたびに 1 つの**モード**が決まり、各モードに**遮断周波数**があります。[位相速度](/physics/mechanics/oscillations-and-waves#def-dispersion)（波の山の動く速さ）は $c$ より速く、[群速度](/physics/mechanics/oscillations-and-waves#def-dispersion)（波のかたまりの動く速さ）とエネルギー速度は $c$ より遅く、積は $c^2$ です。
- 干渉と回折は同じ現象です。ヘルムホルツ方程式（波動方程式を単色の波に限った式）の解に対する**積分公式**（グリーンの恒等式から証明できる）に開口の仮定を入れ、遠方で近似すると、単スリットは $\mathrm{sinc}^2$ の、二重スリットは $\cos^2$ の縞の強度分布を与えます。

<Aside type="tip" title="この記事で使う事前知識">

- **平面波**: 波面が平らで、位相が $\boldsymbol{k}\cdot\boldsymbol{x}-\omega t$ で決まる波。→ [平面波](/physics/electromagnetism/maxwell-equations#def-plane-wave)
- **分散関係**: 角周波数 $\omega$ と波数 $k$ を結ぶ関係式で、位相速度・群速度はここから決まる。→ [分散関係](/physics/mechanics/oscillations-and-waves#def-dispersion)
- **境界条件**: 界面で $\boldsymbol{E}$ と $\boldsymbol{H}$ の接線成分がつながるという条件。→ [物質中の電磁場](/physics/electromagnetism/fields-in-matter)
- **ポインティング・ベクトル**: 電磁場がエネルギーを運ぶ向きと量を表すベクトル。→ [マクスウェル方程式と電磁波](/physics/electromagnetism/maxwell-equations)
- **複素数の指数関数**: $e^{i\theta}=\cos\theta+i\sin\theta$ として、振動を複素数で表す道具。

</Aside>

---

## 1. 動機：光は電磁波である、の先へ

[マクスウェル方程式と電磁波](/physics/electromagnetism/maxwell-equations) で、真空中の電磁場が速さ $c$ で伝わる横波であることを示しました。ところが、私たちが日常で出会う光は、ガラスの中で曲がり、水面で反射し、スリットを抜けて縞をつくり、光ファイバーの中を何キロも進みます。これらはすべて、物質が入ったときのマクスウェル方程式の帰結です。

[物質中の電磁場](/physics/electromagnetism/fields-in-matter) で、物質を $\varepsilon,\mu$ という 2 つの定数にまとめた方程式と、界面で場がつながる条件を手に入れました。この記事ではそれを 3 つの舞台に適用します。

1. 無限に広い一様な媒質（§3）と、2 つの媒質の界面（§4）。
2. 金属の壁に囲まれた管（§5）。
3. 障害物のある空間（§6）。

どれも「方程式を解き、境界で場をつなぐ」という同じ手順で解けます。違いは、境界の形と、使う境界条件だけです。

---

## 2. 準備：複素表示と時間平均

単色の波、つまり時間依存性が $\cos(\omega t+\text{位相})$ の形の場だけを考えます。実の場を複素関数の実部として $\boldsymbol{E}(\boldsymbol{x},t)=\mathrm{Re}\bigl[\boldsymbol{E}(\boldsymbol{x})e^{-i\omega t}\bigr]$ と書きます。マクスウェル方程式は線形なので、複素場がそれを満たすなら実部も満たします。時間微分は $\partial_t\to-i\omega$ に置き換わります。

2 つの量の積（たとえばポインティング・ベクトル）は非線形なので、実部をとってから掛けなければなりません。1 周期の時間平均は

$$
\bigl\langle\mathrm{Re}(ae^{-i\omega t})\,\mathrm{Re}(be^{-i\omega t})\bigr\rangle=\frac12\mathrm{Re}(a\bar b)
$$

です。実際、$\mathrm{Re}(ae^{-i\omega t})=\frac12(ae^{-i\omega t}+\bar ae^{i\omega t})$ なので積は $\frac14(a b e^{-2i\omega t}+a\bar b+\bar ab+\bar a\bar be^{2i\omega t})$、$e^{\pm2i\omega t}$ の項は 1 周期で平均すると消え、残るのは $\frac14(a\bar b+\bar ab)=\frac12\mathrm{Re}(a\bar b)$ です。

媒質は**一様**（$\varepsilon,\mu$ が場所によらない）、**線形・等方**（$\boldsymbol{D}=\varepsilon\boldsymbol{E}$、$\boldsymbol{B}=\mu\boldsymbol{H}$）とし、与える周波数 $\omega$ の範囲では $\varepsilon,\mu$ は定数（分散を無視）とします。自由電荷・自由電流はないとします（導電性は §3.3 で扱います）。

---

## 3. 媒質中の平面波

### 3.1. 波動方程式と屈折率

<Theorem id="thm-medium-wave" title="媒質中の波動方程式">
一様な線形等方媒質で $\rho_f=0$、$\boldsymbol{J}_f=\boldsymbol{0}$ のとき、$C^2$ 級の場 $\boldsymbol{E},\boldsymbol{B}$ は
$$
\nabla^2\boldsymbol{E}=\varepsilon\mu\,\frac{\partial^2\boldsymbol{E}}{\partial t^2},\qquad\nabla^2\boldsymbol{B}=\varepsilon\mu\,\frac{\partial^2\boldsymbol{B}}{\partial t^2}
$$
を満たす。位相速度は $v=1/\sqrt{\varepsilon\mu}=c/n$、$n=\sqrt{\varepsilon_r\mu_r}$（屈折率）である。
</Theorem>

証明の方針は、真空で光速を出した計算の $\varepsilon_0\mu_0$ を $\varepsilon\mu$ に置き換えるだけです。

<Proof of="thm-medium-wave">
<Ref to="physics/electromagnetism/fields-in-matter#thm-macroscopic" /> の方程式は、$\boldsymbol{D}=\varepsilon\boldsymbol{E}$、$\boldsymbol{H}=\boldsymbol{B}/\mu$、$\rho_f=0$、$\boldsymbol{J}_f=\boldsymbol{0}$ のもとで $\nabla\cdot\boldsymbol{E}=0$、$\nabla\cdot\boldsymbol{B}=0$、$\nabla\times\boldsymbol{E}=-\partial_t\boldsymbol{B}$、$\nabla\times\boldsymbol{B}=\varepsilon\mu\,\partial_t\boldsymbol{E}$ となり、真空のマクスウェル方程式の $\varepsilon_0\mu_0$ を $\varepsilon\mu$ に替えただけの形です。そこで <Ref to="physics/electromagnetism/maxwell-equations#thm-wave-equation" /> と同じ計算を行います。<Ref to="physics/electromagnetism/maxwell-equations#lem-curl-curl" /> と $\nabla\cdot\boldsymbol{E}=0$ より $\nabla\times(\nabla\times\boldsymbol{E})=-\nabla^2\boldsymbol{E}$ で、一方 $\nabla\times(\nabla\times\boldsymbol{E})=-\partial_t(\nabla\times\boldsymbol{B})=-\varepsilon\mu\,\partial_t^2\boldsymbol{E}$ です。$\boldsymbol{B}$ も同様です。波動方程式 $\nabla^2f=v^{-2}\partial_t^2f$ の速さは $v=1/\sqrt{\varepsilon\mu}$ で、$\varepsilon_0\mu_0=1/c^2$ より $v=c/\sqrt{\varepsilon_r\mu_r}$ です。
</Proof>

ガラスでは $n\approx1.5$ で、光速は $2\times10^8\,\mathrm{m/s}$ ほどです。水の静電的な $\varepsilon_r\approx80$ から $\sqrt{80}\approx9$ と予想される屈折率が、可視光では実際には $1.33$ なのは、水分子の電気的応答が光の周波数（$10^{15}\,\mathrm{Hz}$）では追従できず、$\varepsilon$ が周波数に依存するからです。これを**分散**といいます。この記事では、考える周波数での $\varepsilon(\omega)$ の値をそのまま $\varepsilon$ として使います。

<Proposition id="prop-medium-plane-wave" title="媒質中の平面波">
実ベクトル $\boldsymbol{k}$ と複素ベクトル $\boldsymbol{E}_0$ について、$\boldsymbol{E}=\boldsymbol{E}_0e^{i(\boldsymbol{k}\cdot\boldsymbol{x}-\omega t)}$ が <Ref to="thm-medium-wave" /> の前提の方程式系の解になる（$\boldsymbol{B}$ を適当に選べる）ための必要十分条件は、
$$
\boldsymbol{k}\cdot\boldsymbol{E}_0=0,\qquad|\boldsymbol{k}|=\frac{n\omega}{c}
$$
である。このとき $\boldsymbol{B}=\boldsymbol{k}\times\boldsymbol{E}/\omega$、$\boldsymbol{H}=\boldsymbol{B}/\mu$ で、$\boldsymbol{E},\boldsymbol{H},\boldsymbol{k}$ は互いに直交し、$|\boldsymbol{H}|=|\boldsymbol{E}|/Z$（$Z=\sqrt{\mu/\varepsilon}$ は媒質の波動インピーダンス）である。時間平均のポインティング・ベクトルは
$$
\langle\boldsymbol{S}\rangle=\frac12\mathrm{Re}(\boldsymbol{E}\times\bar{\boldsymbol{H}})=\frac{|\boldsymbol{E}_0|^2}{2Z}\,\hat{\boldsymbol{k}} ,
$$
非磁性媒質（$\mu=\mu_0$）では強度は $I=|\langle\boldsymbol{S}\rangle|=\frac12n\varepsilon_0c\,|\boldsymbol{E}_0|^2$ である。
</Proposition>

<Proof of="prop-medium-plane-wave">
$\nabla\to i\boldsymbol{k}$、$\partial_t\to-i\omega$ と置きます。$\nabla\cdot\boldsymbol{E}=0$ は $\boldsymbol{k}\cdot\boldsymbol{E}_0=0$、ファラデーの法則 $\nabla\times\boldsymbol{E}=-\partial_t\boldsymbol{B}$ は $i\boldsymbol{k}\times\boldsymbol{E}=i\omega\boldsymbol{B}$ すなわち $\boldsymbol{B}=\boldsymbol{k}\times\boldsymbol{E}/\omega$（これが $\nabla\cdot\boldsymbol{B}=0$ も満たす）、アンペール・マクスウェルの法則 $\nabla\times\boldsymbol{B}=\varepsilon\mu\,\partial_t\boldsymbol{E}$ は $i\boldsymbol{k}\times\boldsymbol{B}=-i\omega\varepsilon\mu\boldsymbol{E}$ となります。左辺に $\boldsymbol{B}$ を入れ、$\boldsymbol{k}\cdot\boldsymbol{E}=0$ のときの $\boldsymbol{k}\times(\boldsymbol{k}\times\boldsymbol{E})=-|\boldsymbol{k}|^2\boldsymbol{E}$ を使うと $-|\boldsymbol{k}|^2\boldsymbol{E}/\omega=-\omega\varepsilon\mu\boldsymbol{E}$、つまり $|\boldsymbol{k}|^2=\varepsilon\mu\omega^2=n^2\omega^2/c^2$ です。逆向きも同じ式変形をたどれます。

直交性は $\boldsymbol{B}=\boldsymbol{k}\times\boldsymbol{E}/\omega$ から、大きさは $|\boldsymbol{H}|=|\boldsymbol{k}||\boldsymbol{E}|/(\omega\mu)=\sqrt{\varepsilon\mu}|\boldsymbol{E}|/\mu=|\boldsymbol{E}|/Z$ です。時間平均は §2 の公式から $\frac12\mathrm{Re}(\boldsymbol{E}\times\bar{\boldsymbol{H}})$ で、$\boldsymbol{E}\times(\hat{\boldsymbol{k}}\times\bar{\boldsymbol{E}})=\hat{\boldsymbol{k}}|\boldsymbol{E}|^2$（$\hat{\boldsymbol{k}}\cdot\boldsymbol{E}=0$ より）から $\frac{|\boldsymbol{E}_0|^2}{2Z}\hat{\boldsymbol{k}}$ を得ます。非磁性では $1/Z=\sqrt{\varepsilon/\mu_0}=n\sqrt{\varepsilon_0/\mu_0}=n\varepsilon_0c$ です。
</Proof>

波動インピーダンス $Z$ とは、波の中で電場と磁場の大きさがどんな比になっているかを表す、媒質ごとの定数です。

### 3.2. 導電性のある媒質

金属のように自由電荷が動ける媒質では、オームの法則 $\boldsymbol{J}_f=\sigma\boldsymbol{E}$（$\sigma$ は電気伝導度）を構成方程式に加えます。電荷保存則と <Ref to="physics/electromagnetism/fields-in-matter#thm-gauss-displacement" /> から $\partial_t\rho_f=-\sigma\rho_f/\varepsilon$ となるので、自由電荷は時間 $\varepsilon/\sigma$（銅で $10^{-19}\,\mathrm{s}$ ほど）で消え、$\rho_f=0$ としてよいことが分かります。

<Proposition id="prop-conducting-wave" title="導電性媒質中の平面波と表皮厚さ">
$\boldsymbol{J}_f=\sigma\boldsymbol{E}$、$\rho_f=0$ の一様な媒質では、$\boldsymbol{E}=\boldsymbol{E}_0e^{i(kz-\omega t)}$（$\boldsymbol{E}_0\perp\boldsymbol{e}_z$）が解になる条件は、複素波数 $k$ が
$$
k^2=\varepsilon\mu\omega^2\Bigl(1+\frac{i\sigma}{\varepsilon\omega}\Bigr)
$$
を満たすことである。良導体（$\sigma\gg\varepsilon\omega$）では $k\approx(1+i)/\delta$、$\delta=\sqrt{2/(\mu\sigma\omega)}$ で、振幅は $e^{-z/\delta}$ で減衰する。
</Proposition>

<Proof of="prop-conducting-wave">
$\nabla\times\boldsymbol{B}=\mu\boldsymbol{J}_f+\varepsilon\mu\partial_t\boldsymbol{E}=\mu\sigma\boldsymbol{E}+\varepsilon\mu\partial_t\boldsymbol{E}$ の回転をとり、<Ref to="thm-medium-wave" /> の証明と同じ手順で $\nabla^2\boldsymbol{E}=\varepsilon\mu\,\partial_t^2\boldsymbol{E}+\mu\sigma\,\partial_t\boldsymbol{E}$ を得ます。$\nabla\to ik\boldsymbol{e}_z$、$\partial_t\to-i\omega$ で $-k^2=-\varepsilon\mu\omega^2-i\mu\sigma\omega$、すなわち $k^2=\varepsilon\mu\omega^2+i\mu\sigma\omega$ です。$\sigma\gg\varepsilon\omega$ では $k^2\approx i\mu\sigma\omega$ で、$i=e^{i\pi/2}$ の平方根 $e^{i\pi/4}=(1+i)/\sqrt2$ から $k\approx(1+i)\sqrt{\mu\sigma\omega/2}=(1+i)/\delta$ です。$e^{ikz}=e^{iz/\delta}e^{-z/\delta}$ より主張を得ます。
</Proof>

銅（$\sigma=5.8\times10^7\,\mathrm{S/m}$、$\mu\approx\mu_0$）の表皮厚さは、$60\,\mathrm{Hz}$ で約 $8.5\,\mathrm{mm}$、$1\,\mathrm{GHz}$ で約 $2.1\,\mu\mathrm{m}$ です。高周波では電流が導体の表面の薄い層にしか流れません。

---

## 4. 反射と屈折

### 4.1. 位相がそろう条件

以下、2 つの媒質はどちらも非磁性（$\mu=\mu_0$）で透明（$n_1,n_2$ は実数）とします。界面を $z=0$ とし、$z<0$ に媒質 1（屈折率 $n_1$）、$z>0$ に媒質 2（$n_2$）があります。媒質 1 から平面波が入射します。

$$
\text{入射:}\ \boldsymbol{E}_ie^{i(\boldsymbol{k}_i\cdot\boldsymbol{x}-\omega t)},\quad\text{反射:}\ \boldsymbol{E}_re^{i(\boldsymbol{k}_r\cdot\boldsymbol{x}-\omega t)},\quad\text{透過:}\ \boldsymbol{E}_te^{i(\boldsymbol{k}_t\cdot\boldsymbol{x}-\omega t)}.
$$

媒質が線形で時間に依らないので、3 つの波の角周波数は同じ $\omega$ です。<Ref to="prop-medium-plane-wave" /> により $|\boldsymbol{k}_i|=|\boldsymbol{k}_r|=n_1\omega/c$、$|\boldsymbol{k}_t|=n_2\omega/c$ です。界面での接続条件は $\boldsymbol{E}_\parallel$、$\boldsymbol{H}_\parallel$ の連続（<Ref to="physics/electromagnetism/fields-in-matter#thm-boundary" />、自由電荷・自由電流なし）で、$z=0$ 上の任意の点 $\boldsymbol{\xi}=(x,y)$ と時刻で成り立つ必要があります。ここで次の補題を使います。補題の中身は「異なる波数の波は互いに打ち消し合えない」というだけのことで、証明は線形代数の道具でそれを確かめる技術的なものです。結論だけ受け取って先へ進んでも構いません。

<Lemma id="lem-phase-matching" title="指数関数の一次独立性">
$q_1,\dots,q_N\in\mathbb{R}^2$ が互いに異なるとき、複素数 $c_1,\dots,c_N$ について、すべての $\boldsymbol{\xi}\in\mathbb{R}^2$ で $\sum_jc_je^{i\boldsymbol{q}_j\cdot\boldsymbol{\xi}}=0$ ならば $c_1=\dots=c_N=0$ である。
</Lemma>

<Proof of="lem-phase-matching">
$\boldsymbol{q}_j-\boldsymbol{q}_l$ ($j\ne l$) のすべてと直交しない方向 $\boldsymbol{v}$ を取れます（有限個の直線以外の方向）。すると $\lambda_j=\boldsymbol{q}_j\cdot\boldsymbol{v}$ は互いに異なります。$\boldsymbol{\xi}=s\boldsymbol{v}$ として $s$ の関数 $\sum_jc_je^{i\lambda_js}=0$ を $s=0$ で $m=0,\dots,N-1$ 回微分すると $\sum_jc_j\lambda_j^m=0$ で、係数行列は互いに異なる $\lambda_j$ のヴァンデルモンド行列です。その行列式は $\prod_{j<l}(\lambda_l-\lambda_j)\ne0$ なので $c_j=0$ です。
</Proof>

接続条件は、たとえば $\boldsymbol{E}_{i\parallel}e^{i\boldsymbol{k}_i^\parallel\cdot\boldsymbol{\xi}}+\boldsymbol{E}_{r\parallel}e^{i\boldsymbol{k}_r^\parallel\cdot\boldsymbol{\xi}}=\boldsymbol{E}_{t\parallel}e^{i\boldsymbol{k}_t^\parallel\cdot\boldsymbol{\xi}}$ の形で（$\boldsymbol{H}$ も同様）、$\boldsymbol{k}^\parallel$ は波数ベクトルの界面に平行な成分です。<Ref to="lem-phase-matching" /> により、接線波数が異なる指数関数の係数は個別に $0$ でなければなりません。入射波は横波で振幅が $0$ でなく、$\boldsymbol{E}_{i\parallel}$ と $\boldsymbol{H}_{i\parallel}$ のうち少なくとも一方は $0$ ではない（両方とも法線方向なら、$\boldsymbol{E}\perp\boldsymbol{H}$ に反する）ので、入射波の接線波数は反射波か透過波のどちらかと一致します。一致しない波は、$\boldsymbol{E}_\parallel$ と $\boldsymbol{H}_\parallel$ の両方の係数が $0$、つまり存在しません。したがって、**存在する波はすべて同じ接線波数 $\boldsymbol{k}^\parallel$ をもちます**。

<Theorem id="thm-snell" title="反射の法則とスネルの法則">
入射角 $\theta_i$、反射角 $\theta_r$、屈折角 $\theta_t$ を法線（$z$ 軸）と波数ベクトルのなす角とする。反射波と透過波がともに存在するとき、波数ベクトルはすべて入射面（$\boldsymbol{k}_i$ と法線がはる平面）内にあり、
$$
\theta_r=\theta_i,\qquad n_1\sin\theta_i=n_2\sin\theta_t .
$$
</Theorem>

<Proof of="thm-snell">
同じ接線波数 $\boldsymbol{k}^\parallel$ をもち、反射波と透過波の接線波数は入射波と一致するので、$\boldsymbol{k}_r,\boldsymbol{k}_t$ も $\boldsymbol{k}_i$ と同じ入射面内にあります（残る成分は法線成分 $k_z$ だけで、符号は進む向きで決まる）。$x$ 軸を入射面内の界面方向に取ると、$|\boldsymbol{k}^\parallel|=\frac{n_1\omega}{c}\sin\theta_i$（入射波）$=\frac{n_1\omega}{c}\sin\theta_r$（反射波）$=\frac{n_2\omega}{c}\sin\theta_t$（透過波）です。$0\le\theta_r\le\pi/2$ としたとき、最初の等号から $\theta_r=\theta_i$、2 つ目の等号から $n_1\sin\theta_i=n_2\sin\theta_t$ です。
</Proof>

<Figure caption="反射と屈折。3 つの波の接線方向の波数（x 成分）が等しいことが、反射の法則とスネルの法則を与える。">
<svg viewBox="0 0 520 330" width="100%" role="img" aria-label="界面 z=0 で入射光が反射光と屈折光に分かれる様子。入射角 θi、反射角 θr、屈折角 θt が描かれている。">
  <defs>
    <marker id="wo-arrow" viewBox="0 0 10 10" refX="9" refY="5" markerWidth="7" markerHeight="7" orient="auto">
      <path d="M 0 0 L 10 5 L 0 10 z" fill="currentColor" />
    </marker>
  </defs>
  <line x1="20" y1="165" x2="500" y2="165" stroke="currentColor" stroke-width="2" />
  <line x1="260" y1="40" x2="260" y2="290" stroke="currentColor" stroke-width="1" stroke-dasharray="4 4" />
  <g stroke="var(--sl-color-accent)" stroke-width="2" fill="none" marker-end="url(#wo-arrow)">
    <line x1="140" y1="50" x2="256" y2="161" />
    <line x1="264" y1="161" x2="380" y2="50" />
    <line x1="264" y1="169" x2="340" y2="280" />
  </g>
  <g fill="currentColor" font-size="15">
    <text x="60" y="140">媒質 1（n₁）</text>
    <text x="60" y="200">媒質 2（n₂）</text>
    <text x="200" y="100">θᵢ</text>
    <text x="290" y="100">θᵣ</text>
    <text x="280" y="225">θₜ</text>
    <text x="110" y="40">入射</text>
    <text x="390" y="45">反射</text>
    <text x="345" y="290">透過</text>
  </g>
</svg>
</Figure>

### 4.2. フレネルの式

振幅の関係は、接続条件を直接解いて求めます。方針は単純で、界面で連続な量が 2 つ（電場の接線成分と磁場の接線成分）あるので、未知数の反射波と透過波の振幅 2 つについての連立一次方程式を解くだけです。後の証明は成分と符号の計算が続きますが、各偏光で書き下す式は 2 本だけです。入射面（$xz$ 面）に対して偏光を 2 つに分けます。

- **$s$ 偏光**（TE）: 電場が入射面に垂直（$y$ 方向）。
- **$p$ 偏光**（TM）: 磁場が入射面に垂直（$y$ 方向）で、電場は入射面内にある。

一般の偏光は両者の重ね合わせで、境界条件は線形なので 2 つは独立に扱えます。

<Theorem id="thm-fresnel" title="フレネルの式">
$a=n_1\cos\theta_i$、$b=n_2\cos\theta_t$ とおく。

**$s$ 偏光**: 反射と透過の電場振幅比は
$$
r_s=\frac{E_r}{E_i}=\frac{n_1\cos\theta_i-n_2\cos\theta_t}{n_1\cos\theta_i+n_2\cos\theta_t},\qquad t_s=\frac{E_t}{E_i}=\frac{2n_1\cos\theta_i}{n_1\cos\theta_i+n_2\cos\theta_t}.
$$

**$p$ 偏光**: 磁場 $H_y$ の反射振幅比と、電場の透過振幅比は
$$
r_p=\frac{H_r}{H_i}=\frac{n_2\cos\theta_i-n_1\cos\theta_t}{n_2\cos\theta_i+n_1\cos\theta_t},\qquad t_p=\frac{E_t}{E_i}=\frac{2n_1\cos\theta_i}{n_2\cos\theta_i+n_1\cos\theta_t}.
$$
</Theorem>

<Proof of="thm-fresnel">
$\boldsymbol{k}=(k_x,0,k_z)=\frac{n\omega}{c}(\sin\theta,0,\pm\cos\theta)$ とします。<Ref to="prop-medium-plane-wave" /> から $\boldsymbol{H}=\dfrac{n}{\mu_0c}\hat{\boldsymbol{k}}\times\boldsymbol{E}$ で、逆に $\boldsymbol{E}=-\dfrac{\mu_0c}{n}\hat{\boldsymbol{k}}\times\boldsymbol{H}$ です。

**$s$ 偏光。** $\boldsymbol{E}=E\boldsymbol{e}_y$。$\hat{\boldsymbol{k}}\times\boldsymbol{e}_y=(\mp\cos\theta,0,\sin\theta)$ （複号は $k_z$ の符号）より、$H_x=\mp\dfrac{n\cos\theta}{\mu_0c}E$ です。入射波（$k_z>0$）では $H_{ix}=-\frac{n_1\cos\theta_i}{\mu_0c}E_i$、反射波（$k_z<0$、$\theta_r=\theta_i$）では $H_{rx}=+\frac{n_1\cos\theta_i}{\mu_0c}E_r$、透過波では $H_{tx}=-\frac{n_2\cos\theta_t}{\mu_0c}E_t$ です。$z=0$ で接線成分 $E_y$ と $H_x$ が連続なので
$$
E_i+E_r=E_t,\qquad-n_1\cos\theta_i(E_i-E_r)=-n_2\cos\theta_tE_t .
$$
すなわち $a(E_i-E_r)=bE_t$ で、$E_t=E_i+E_r$ を代入すると $a(E_i-E_r)=b(E_i+E_r)$、$E_r=\dfrac{a-b}{a+b}E_i$、$E_t=\dfrac{2a}{a+b}E_i$ です。

**$p$ 偏光。** $\boldsymbol{H}=H\boldsymbol{e}_y$。$E_x=-\dfrac{\mu_0c}{n}(\hat{\boldsymbol{k}}\times\boldsymbol{H})_x=-\dfrac{\mu_0c}{n}H\cdot(\mp\cos\theta)$ より、$E_{ix}=\frac{\mu_0c}{n_1}\cos\theta_i\,H_i$、$E_{rx}=-\frac{\mu_0c}{n_1}\cos\theta_i\,H_r$、$E_{tx}=\frac{\mu_0c}{n_2}\cos\theta_t\,H_t$ です。接線成分 $H_y$ と $E_x$ の連続は
$$
H_i+H_r=H_t,\qquad\frac{\cos\theta_i}{n_1}(H_i-H_r)=\frac{\cos\theta_t}{n_2}H_t .
$$
$H_t=H_i+H_r$ を代入して $H_r\Bigl(\dfrac{\cos\theta_i}{n_1}+\dfrac{\cos\theta_t}{n_2}\Bigr)=H_i\Bigl(\dfrac{\cos\theta_i}{n_1}-\dfrac{\cos\theta_t}{n_2}\Bigr)$、よって $r_p=\dfrac{n_2\cos\theta_i-n_1\cos\theta_t}{n_2\cos\theta_i+n_1\cos\theta_t}$ です。透過波の磁場は $H_t=(1+r_p)H_i=\dfrac{2n_2\cos\theta_i}{n_2\cos\theta_i+n_1\cos\theta_t}H_i$ で、電場の大きさは $|E|=\frac{\mu_0c}{n}|H|$ なので $E_t/E_i=(n_1/n_2)(H_t/H_i)$ より $t_p$ を得ます。

法線成分（$D_z$ と $B_z$）の連続は、これらの解が自動的に満たします。$s$ 偏光では $D_z=0$、$B_z=k_xE_y/\omega$ で、$k_x$ が 3 つの波で等しいことと $E_y$ の連続から従います。$p$ 偏光では $B_z=0$、$D_z=\varepsilon E_z=-(n\sin\theta/c)H_y$ で、$n\sin\theta$ が両側で等しい（スネルの法則）ことと $H_y$ の連続から従います。
</Proof>

<Corollary id="cor-brewster" title="ブリュースター角">
$n_1\ne n_2$ のとき、$\tan\theta_B=n_2/n_1$ となる入射角 $\theta_B$ で $r_p=0$ である。このとき $\theta_B+\theta_t=\pi/2$ で、反射光は $s$ 偏光だけになる。
</Corollary>

<Proof of="cor-brewster">
$\tan\theta_i=n_2/n_1$ のとき $\sin\theta_i=n_2/\sqrt{n_1^2+n_2^2}$、$\cos\theta_i=n_1/\sqrt{n_1^2+n_2^2}$ です。スネルの法則から $\sin\theta_t=\dfrac{n_1}{n_2}\sin\theta_i=\dfrac{n_1}{\sqrt{n_1^2+n_2^2}}$、したがって $\cos\theta_t=\dfrac{n_2}{\sqrt{n_1^2+n_2^2}}=\sin\theta_i$（$\theta_t<\pi/2$ のとき）で $\theta_t=\pi/2-\theta_i$ です。このとき $n_2\cos\theta_i=\dfrac{n_1n_2}{\sqrt{n_1^2+n_2^2}}=n_1\cos\theta_t$ で、$r_p$ の分子は $0$ です。
</Proof>

空気からガラス（$n=1.5$）への入射では $\theta_B=\arctan1.5\approx56.3^\circ$ で、水面や路面の反射光（ギラつき）を偏光サングラスが切るのは、反射光がこの付近で $s$ 偏光（水平偏光）に偏っているためです。

<Proposition id="prop-energy-conservation" title="反射率と透過率">
反射率 $R$ を反射波の強度の入射波の強度に対する比、透過率 $T$ を透過波が界面の法線方向に運ぶ時間平均エネルギーの流れの、入射波の同じ量に対する比とする。$s$、$p$ 偏光のいずれでも
$$
R=|r|^2,\qquad T=\frac{n_2\cos\theta_t}{n_1\cos\theta_i}\,|t|^2,\qquad R+T=1 .
$$
</Proposition>

<Proof of="prop-energy-conservation">
<Ref to="prop-medium-plane-wave" /> から、時間平均のエネルギーの流れの大きさは $\frac12n\varepsilon_0c|E|^2$ で、$z$ 方向成分はそれに $\cos\theta$ を掛けたものです。反射波は同じ媒質中を同じ角度で進むので $R=|E_r/E_i|^2=|r|^2$ ($p$ 偏光でも $|E_r/E_i|=|H_r/H_i|$)、透過波は $T=\dfrac{n_2\cos\theta_t|E_t|^2}{n_1\cos\theta_i|E_i|^2}$ です。$s$ 偏光では $a=n_1\cos\theta_i$、$b=n_2\cos\theta_t$ として $T=\dfrac ba\cdot\dfrac{4a^2}{(a+b)^2}=\dfrac{4ab}{(a+b)^2}$、$R=\dfrac{(a-b)^2}{(a+b)^2}$ で、$R+T=\dfrac{(a-b)^2+4ab}{(a+b)^2}=1$ です。$p$ 偏光では $a'=n_2\cos\theta_i$、$b'=n_1\cos\theta_t$ とおくと $T=\dfrac{n_2\cos\theta_t}{n_1\cos\theta_i}\cdot\dfrac{4n_1^2\cos^2\theta_i}{(a'+b')^2}=\dfrac{4a'b'}{(a'+b')^2}$ で、同様に $R+T=1$ です。
</Proof>

垂直入射 $\theta_i=\theta_t=0$ では、どちらの偏光でも $R=\bigl(\dfrac{n_1-n_2}{n_1+n_2}\bigr)^2$ です。空気（$1$）とガラス（$1.5$）の界面では $R=(0.5/2.5)^2=0.04$ で、片面で 4% が反射されます。

### 4.3. 全反射

<Proposition id="prop-tir" title="全反射とエバネッセント波">
$n_1>n_2$ で $n_1\sin\theta_i>n_2$（臨界角 $\theta_c=\arcsin(n_2/n_1)$ を超える）とする。このとき実数の屈折角 $\theta_t$ は存在せず、透過側の場は界面に沿って進みながら $z$ 方向に指数的に減衰する波 $e^{ik_xx}e^{-\kappa z}$、$\kappa=\dfrac\omega c\sqrt{n_1^2\sin^2\theta_i-n_2^2}$ になり、$|r_s|=|r_p|=1$、$T=0$ である。
</Proposition>

<Proof of="prop-tir">
透過波の接線波数は入射波と同じ $k_x=\frac{n_1\omega}{c}\sin\theta_i$ で、$|\boldsymbol{k}_t|^2=k_x^2+k_z^2=\frac{n_2^2\omega^2}{c^2}$ より $k_z^2=\frac{\omega^2}{c^2}(n_2^2-n_1^2\sin^2\theta_i)<0$ です。$k_z=i\kappa$（$\kappa>0$、$z\to\infty$ で有界な枝）となり、場は $e^{ik_xx-\kappa z}$ に比例します。このとき $\cos\theta_t=k_z/|\boldsymbol{k}_t|=i\kappa c/(n_2\omega)$ は純虚数で、$s$ 偏光では $r_s=\dfrac{a-ib'}{a+ib'}$（$a=n_1\cos\theta_i$、$b'=\kappa c/\omega$ はともに実数）、よって $|r_s|=1$ です。$p$ 偏光も同様に分子分母が互いに複素共役となり $|r_p|=1$ です。エバネッセント波のポインティング・ベクトルの $z$ 成分の時間平均は、$\boldsymbol{E}$ と $\boldsymbol{H}$ の位相差が $\pi/2$ になって $0$ です（$R=1$ とも整合します）。
</Proof>

ガラス（$n=1.5$）から空気への臨界角は $\arcsin(1/1.5)\approx41.8^\circ$ です。波長 $600\,\mathrm{nm}$ の光が $45^\circ$ で入射するとき、$\kappa=\frac{2\pi}{600\,\mathrm{nm}}\sqrt{1.5^2\sin^245^\circ-1}\approx3.7\times10^{-3}\,\mathrm{nm}^{-1}$ で、しみ出す距離 $1/\kappa$ は約 $270\,\mathrm{nm}$ です。光ファイバーが曲げても光を漏らさないのは、この全反射によります。

<Figure caption="全反射のときの透過側の場の振幅。界面（z=0）から離れるにつれて e^{-κz} で減衰し、距離 1/κ で 1/e になる。">
<svg viewBox="0 0 520 230" width="100%" role="img" aria-label="界面 z=0 から右へ向かって指数関数的に減衰する振幅の曲線。z=1/κ の位置に縦の目印がある。">
  <line x1="100" y1="20" x2="100" y2="190" stroke="currentColor" stroke-width="2" />
  <line x1="100" y1="160" x2="480" y2="160" stroke="currentColor" stroke-width="1" />
  <path d="M 100 60 L 180 123 L 260 146 L 340 155 L 420 158" stroke="var(--sl-color-accent)" stroke-width="2" fill="none" stroke-linejoin="round" />
  <line x1="180" y1="123" x2="180" y2="160" stroke="currentColor" stroke-width="1" stroke-dasharray="4 4" />
  <g fill="currentColor" font-size="15">
    <text x="20" y="45">媒質 1</text>
    <text x="120" y="40">媒質 2（透過側）</text>
    <text x="108" y="210">z=0</text>
    <text x="165" y="182">1/κ</text>
    <text x="270" y="120">振幅 ∝ e^(−κz)</text>
    <text x="470" y="182">z</text>
  </g>
</svg>
</Figure>

---

## 5. 導波管

金属の管の中を電磁波が伝わる場合を考えます。電子レンジやレーダーで使われる、マイクロ波を運ぶ管の話です。管は壁が完全導体（$\sigma\to\infty$）で、内部は真空で、断面は軸 $z$ に沿って一定とします。

### 5.1. 完全導体の壁での境界条件

<Lemma id="lem-pec-bc" title="完全導体の境界条件">
時間変化する場に対して、完全導体の内部では $\boldsymbol{E}=\boldsymbol{0}$、$\boldsymbol{B}=\boldsymbol{0}$ である。壁の外向き（導体から真空へ向かう）単位法線を $\boldsymbol{n}$ とすると、真空側の壁の表面で
$$
\boldsymbol{n}\times\boldsymbol{E}=\boldsymbol{0},\qquad\boldsymbol{n}\cdot\boldsymbol{B}=0
$$
が成り立つ。
</Lemma>

<Proof of="lem-pec-bc">
導体内の $\boldsymbol{E}$ が $0$ でないと、有限の電場が無限の伝導度 $\sigma$ により無限大の電流を流すので、$\boldsymbol{E}=\boldsymbol{0}$ です。ファラデーの法則から $\partial_t\boldsymbol{B}=-\nabla\times\boldsymbol{E}=\boldsymbol{0}$ で、振動する場には時間によらない成分がないので $\boldsymbol{B}=\boldsymbol{0}$ です。<Ref to="physics/electromagnetism/fields-in-matter#thm-boundary" /> で、導体側の値が $\boldsymbol{E}_2=\boldsymbol{B}_2=\boldsymbol{0}$ なので、真空側の接線成分 $\boldsymbol{E}_\parallel$ と法線成分 $B_n$ は $0$ です。
</Proof>

### 5.2. 横断面の方程式

管の内部では $\rho=0$、$\boldsymbol{J}=\boldsymbol{0}$ で、時間依存性 $e^{-i\omega t}$ の単色の場は $\nabla\times\boldsymbol{E}=i\omega\boldsymbol{B}$、$\nabla\times\boldsymbol{B}=-(i\omega/c^2)\boldsymbol{E}$ を満たします。$z$ 依存性が $e^{ikz}$ の波を考え、場を $\boldsymbol{E}=\boldsymbol{E}_t+E_z\boldsymbol{e}_z$ のように断面内の成分（添字 $t$）と $z$ 成分に分け、断面内の勾配を $\nabla_t=(\partial_x,\partial_y)$ とします。次の命題のねらいは、進行方向の成分 $E_z,B_z$ さえ分かれば、断面内の成分はすべてそこから計算で決まる、ということです。そのため以降は $E_z,B_z$ だけを求めれば足ります。証明の外積の計算は、その代数的な確認です。

<Proposition id="prop-guide-fields" title="縦成分から横成分が決まる">
$\gamma^2:=\omega^2/c^2-k^2\ne0$ のとき、
$$
\boldsymbol{E}_t=\frac{i}{\gamma^2}\Bigl[k\nabla_tE_z-\omega\,\boldsymbol{e}_z\times\nabla_tB_z\Bigr],\qquad\boldsymbol{B}_t=\frac{i}{\gamma^2}\Bigl[k\nabla_tB_z+\frac{\omega}{c^2}\,\boldsymbol{e}_z\times\nabla_tE_z\Bigr]
$$
であり、$f=E_z$ と $f=B_z$ は断面上で $(\nabla_t^2+\gamma^2)f=0$ を満たす。
</Proposition>

<Proof of="prop-guide-fields">
2 つの回転の式の断面内成分をとります。$\nabla=\nabla_t+ik\boldsymbol{e}_z$ で、$\nabla\times(E_z\boldsymbol{e}_z)=\nabla_tE_z\times\boldsymbol{e}_z=-\boldsymbol{e}_z\times\nabla_tE_z$ なので
$$
\text{(a)}\ \ ik\,\boldsymbol{e}_z\times\boldsymbol{E}_t-\boldsymbol{e}_z\times\nabla_tE_z=i\omega\boldsymbol{B}_t,\qquad\text{(b)}\ \ ik\,\boldsymbol{e}_z\times\boldsymbol{B}_t-\boldsymbol{e}_z\times\nabla_tB_z=-\frac{i\omega}{c^2}\boldsymbol{E}_t
$$
です（$\nabla_t\times\boldsymbol{E}_t$ は $z$ 方向なので断面内成分には現れません）。(a) に左から $\boldsymbol{e}_z\times$ を掛け、断面内のベクトル $\boldsymbol{V}$ について $\boldsymbol{e}_z\times(\boldsymbol{e}_z\times\boldsymbol{V})=-\boldsymbol{V}$ であることを使うと、$-ik\boldsymbol{E}_t+\nabla_tE_z=i\omega\,\boldsymbol{e}_z\times\boldsymbol{B}_t$ です。これを (b) に代入して
$$
\frac k\omega\bigl(-ik\boldsymbol{E}_t+\nabla_tE_z\bigr)-\boldsymbol{e}_z\times\nabla_tB_z=-\frac{i\omega}{c^2}\boldsymbol{E}_t
$$
すなわち $\dfrac{i\gamma^2}{\omega}\boldsymbol{E}_t=\boldsymbol{e}_z\times\nabla_tB_z-\dfrac k\omega\nabla_tE_z$ となり、$\boldsymbol{E}_t$ の式を得ます。(a),(b) は置き換え $\boldsymbol{E}\to c\boldsymbol{B}$、$c\boldsymbol{B}\to-\boldsymbol{E}$ で互いに移り合うので、$\boldsymbol{B}_t$ の式は同じ置き換えで得られます。

$z$ 成分については、真空中の各成分が波動方程式 $\nabla^2f=c^{-2}\partial_t^2f$ を満たし（<Ref to="physics/electromagnetism/maxwell-equations#thm-wave-equation" />）、$f\propto e^{i(kz-\omega t)}$ なので $(\nabla_t^2-k^2+\omega^2/c^2)f=0$ です。
</Proof>

### 5.3. 矩形導波管

断面が $0<x<a$、$0<y<b$ の矩形の管を考えます。TM と TE は、磁場（TM）または電場（TE）の進行方向成分が $0$ のモードという意味で、添字の $m,n$ は $x,y$ 方向に断面内で立つ定在波の半波長の数を表します。以下の証明の大半は、この形の場が壁の条件を満たすことの確認です。

<Theorem id="thm-rect-modes" title="矩形導波管のモード">
$\gamma_{mn}^2=\pi^2\bigl(m^2/a^2+n^2/b^2\bigr)$ とし、$k^2=\omega^2/c^2-\gamma_{mn}^2$ とする。次の場は <Ref to="lem-pec-bc" /> の壁の境界条件を満たす真空のマクスウェル方程式の解である。

- **TM$_{mn}$ モード**（$m,n\ge1$）: $B_z=0$、$E_z=E_0\sin\dfrac{m\pi x}{a}\sin\dfrac{n\pi y}{b}\,e^{i(kz-\omega t)}$。
- **TE$_{mn}$ モード**（$m,n\ge0$、$(m,n)\ne(0,0)$）: $E_z=0$、$B_z=B_0\cos\dfrac{m\pi x}{a}\cos\dfrac{n\pi y}{b}\,e^{i(kz-\omega t)}$。

横成分は <Ref to="prop-guide-fields" /> で与えられる。
</Theorem>

<Proof of="thm-rect-modes">
$\psi=\sin(m\pi x/a)\sin(n\pi y/b)$ も $\psi=\cos(m\pi x/a)\cos(n\pi y/b)$ も $\nabla_t^2\psi=-\gamma_{mn}^2\psi$ を満たすので、<Ref to="prop-guide-fields" /> の最後の条件は $\gamma=\gamma_{mn}$ のとき満たされます。$\gamma_{mn}>0$ なので横成分の式も使えます。

壁での条件を確かめます。壁の接線方向（断面内）を $\hat{\boldsymbol{s}}$、法線を $\hat{\boldsymbol{n}}$ とします。**TM**: $E_z=0$ が 4 つの壁すべてで成り立ち（$x=0,a$ では $\sin=0$、$y=0,b$ では $\sin=0$）、断面内の接線成分は $\hat{\boldsymbol{s}}\cdot\boldsymbol{E}_t=(ik/\gamma^2)\partial_sE_z=0$ です。$B_z=0$ なので $\boldsymbol{B}_t=\frac{i\omega}{c^2\gamma^2}\boldsymbol{e}_z\times\nabla_tE_z$ の法線成分は $\hat{\boldsymbol{n}}\cdot(\boldsymbol{e}_z\times\nabla E_z)=\nabla E_z\cdot(\hat{\boldsymbol{n}}\times\boldsymbol{e}_z)\propto\partial_sE_z=0$ です。**TE**: $E_z=0$ は自明、$\hat{\boldsymbol{s}}\cdot\boldsymbol{E}_t=-\frac{i\omega}{\gamma^2}\hat{\boldsymbol{s}}\cdot(\boldsymbol{e}_z\times\nabla_tB_z)=\mp\frac{i\omega}{\gamma^2}\partial_nB_z$ で、$\cos$ 関数の壁での法線微分（$x=0,a$ での $\partial_x\cos(m\pi x/a)\propto\sin=0$ など）は $0$ なので成り立ちます。$B_n=\frac{ik}{\gamma^2}\partial_nB_z=0$ も同じ理由で成り立ちます。
</Proof>

<Remark id="rem-modes-complete">
矩形管の場がこれらのモードの重ね合わせですべて尽くされること（固有関数の完全性）は、本記事では認めます。断面上の関数 $\sin\frac{m\pi x}{a}\sin\frac{n\pi y}{b}$ や $\cos\frac{m\pi x}{a}\cos\frac{n\pi y}{b}$ の系が $L^2$ で完全であることは、偶・奇拡張を周期関数と見たフーリエ級数の完全性（<Ref to="mathematics/fourier-analysis/fourier-series#thm-parseval-torus" text="Parseval の等式" /> の $L^2$ 完全性、多次元は <Ref to="mathematics/fourier-analysis/fourier-series#prop-torus-uniqueness" />）から従います。また、E,B の縦成分がどちらも $0$ でない場は、TM と TE の重ね合わせとして得られます。
</Remark>

<Proposition id="prop-no-tem" title="中空の導波管は TEM 波を運ばない">
単連結な断面をもつ中空の完全導体管の中に、$E_z=B_z=0$ で $k=\omega/c\ne 0$ となる（$\gamma=0$ の）非自明な場は存在しない。
</Proposition>

<Proof of="prop-no-tem">
$E_z=B_z=0$ のとき、(a) より $\boldsymbol{B}_t=\frac1c\boldsymbol{e}_z\times\boldsymbol{E}_t$ です。ファラデーの法則の $z$ 成分は $(\nabla_t\times\boldsymbol{E}_t)_z=i\omega B_z=0$、$\nabla\cdot\boldsymbol{E}=0$ は $\nabla_t\cdot\boldsymbol{E}_t=-ikE_z=0$ です。単連結な断面では回転のないベクトル場はスカラー関数の勾配なので（<Ref to="physics/electromagnetism/gauge-potentials#lem-poincare-curl" /> の 2 次元版）、$\boldsymbol{E}_t=-\nabla_t\phi$、$\nabla_t^2\phi=0$ です。壁で $\hat{\boldsymbol{s}}\cdot\boldsymbol{E}_t=0$ ですから $\phi$ は連結な境界上で一定で、<Ref to="physics/electromagnetism/boundary-value-problems#thm-uniqueness" /> (a) により断面全体で $\phi$ は定数、$\boldsymbol{E}_t=\boldsymbol{0}$ です。
</Proof>

TEM 波とは、電場も磁場も進行方向の成分をもたない波のことで、単連結とは断面に穴がないということです。同軸ケーブルのように導体が 2 つ以上あれば断面が単連結でなくなり、境界上の定数が導体ごとに異なってよいので TEM 波が可能になります。中空の管には遮断周波数があり、それより低い周波数の波は通りません。

### 5.4. 遮断周波数と速度

分散関係は $k^2=\omega^2/c^2-\gamma_{mn}^2$ で、**遮断周波数** $\omega_{mn}=c\gamma_{mn}$ を境に性質が変わります。

- $\omega>\omega_{mn}$: $k$ は実数で、$e^{i(kz-\omega t)}$ は $z$ 方向に伝わる波です。
- $\omega<\omega_{mn}$: $k=i\kappa$ ($\kappa=\sqrt{\gamma_{mn}^2-\omega^2/c^2}$) は純虚数で、場は $e^{-\kappa z}$ で減衰し、波は伝わりません。

最も低い遮断周波数をもつモードが基本モードで、$a>b$ の矩形管では TE$_{10}$ です（$\omega_{10}=\pi c/a$）。

<Proposition id="prop-te10" title="TE$_{10}$ モードの速度">
TE$_{10}$ モードの場は
$$
E_y=E_0\sin\frac{\pi x}{a}\,e^{i(kz-\omega t)},\qquad B_x=-\frac k\omega E_y,\qquad B_z=-\frac{i\pi}{a\omega}E_0\cos\frac{\pi x}{a}\,e^{i(kz-\omega t)}
$$
（他の成分は $0$）である。位相速度 $v_p=\omega/k$ と群速度 $v_g=d\omega/dk$ は
$$
v_p=\frac{c}{\sqrt{1-\omega_{10}^2/\omega^2}}>c,\qquad v_g=c\sqrt{1-\frac{\omega_{10}^2}{\omega^2}}<c,\qquad v_pv_g=c^2
$$
で、時間平均のエネルギーの輸送速度（単位長さあたりのエネルギーに対する、断面を通る電力の比）は $v_g$ に等しい。
</Proposition>

<Proof of="prop-te10">
<Ref to="prop-guide-fields" /> に $B_z=B_0\cos(\pi x/a)$、$E_z=0$、$\gamma=\pi/a$ を代入すると、$E_y=\dfrac{i\omega}{\gamma^2}\cdot\dfrac{\pi}{a}B_0\sin\dfrac{\pi x}{a}=\dfrac{i\omega a}{\pi}B_0\sin\dfrac{\pi x}a$、よって $B_0=-\dfrac{i\pi}{a\omega}E_0$ とおけば $E_y$ の式を得ます。$B_x=\dfrac{ik}{\gamma^2}\partial_xB_z=-\dfrac k\omega E_0\sin\dfrac{\pi x}{a}$ です。

分散関係 $\omega^2=c^2k^2+\omega_{10}^2$ から、$v_p=\omega/k=\omega/\sqrt{\omega^2/c^2-\omega_{10}^2/c^2}=c/\sqrt{1-\omega_{10}^2/\omega^2}$、両辺を微分して $\omega\,d\omega=c^2k\,dk$ だから $v_g=c^2k/\omega=c^2/v_p$ です。

エネルギー輸送を計算します。$z$ 方向の時間平均ポインティング・ベクトルは $\langle S_z\rangle=\frac{1}{2\mu_0}\mathrm{Re}(E_x\bar B_y-E_y\bar B_x)=\frac{1}{2\mu_0}\frac k\omega|E_y|^2$ で、断面積分した電力は $P=\dfrac{k}{2\mu_0\omega}|E_0|^2\dfrac{ab}{2}$ です。単位長さあたりのエネルギーは、時間平均の密度 $\frac14\varepsilon_0|\boldsymbol{E}|^2+\frac{1}{4\mu_0}|\boldsymbol{B}|^2$ の断面積分で、
$$
U=\frac{ab}{8}|E_0|^2\Bigl[\varepsilon_0+\frac{1}{\mu_0\omega^2}\Bigl(k^2+\frac{\pi^2}{a^2}\Bigr)\Bigr]=\frac{ab}{8}|E_0|^2\bigl[\varepsilon_0+\varepsilon_0\bigr]=\frac{ab\,\varepsilon_0|E_0|^2}{4}
$$
です（$\sin^2$ と $\cos^2$ の $x$ での積分が $a/2$、$y$ では $b$、$k^2+\pi^2/a^2=\omega^2/c^2$、$1/(\mu_0c^2)=\varepsilon_0$）。比は $P/U=\dfrac{k}{\mu_0\varepsilon_0\omega}=\dfrac{c^2k}{\omega}=v_g$ です。
</Proof>

位相速度が $c$ を超えるのは、位相の模様が壁で反射を繰り返す波の重ね合わせとして生じる見かけの速度だからで、情報もエネルギーも運びません。実際にエネルギーが進む速さは $v_g<c$ です。

<Example id="ex-wr90" title="X バンドの標準導波管 WR-90">
$a=22.86\,\mathrm{mm}$、$b=10.16\,\mathrm{mm}$ の矩形管では、遮断周波数 $f_{mn}=\dfrac c2\sqrt{m^2/a^2+n^2/b^2}$ は TE$_{10}$ が $6.56\,\mathrm{GHz}$、TE$_{20}$ が $13.1\,\mathrm{GHz}$、TE$_{01}$ が $14.8\,\mathrm{GHz}$、TE$_{11}$/TM$_{11}$ が $16.1\,\mathrm{GHz}$ です。したがって $6.56$–$13.1\,\mathrm{GHz}$ の間では TE$_{10}$ だけが伝わり、単一モードの運用ができます。$10\,\mathrm{GHz}$ では $\sqrt{1-(6.56/10)^2}=0.755$ で、$v_g=0.755c$、$v_p=1.32c$ です。
</Example>

---

## 6. 干渉と回折

### 6.1. ヘルムホルツ・キルヒホッフの積分公式

単色の波の場の 1 成分 $u(\boldsymbol{x})e^{-i\omega t}$ は、真空中で $(\nabla^2+k^2)u=0$（$k=\omega/c$）を満たします（<Ref to="physics/electromagnetism/maxwell-equations#thm-wave-equation" />）。この方程式の解を、領域の境界上の値から復元する公式があります。

<Theorem id="thm-kirchhoff" title="ヘルムホルツ・キルヒホッフの積分公式">
$V$ を境界が区分的に滑らかな有界領域とし、$u\in C^2(\overline V)$ が $V$ で $(\nabla^2+k^2)u=0$ を満たすとする。$P\in V$ とし、$R=|\boldsymbol{x}-P|$、$G=e^{ikR}/R$、$\boldsymbol{n}$ を外向き単位法線とすると
$$
u(P)=\frac{1}{4\pi}\oint_{\partial V}\bigl(G\,\partial_nu-u\,\partial_nG\bigr)\,dS .
$$
</Theorem>

<Proof of="thm-kirchhoff">
$V_\epsilon=V\setminus\overline{B_\epsilon(P)}$ に <Ref to="physics/electromagnetism/boundary-value-problems#prop-green" text="グリーンの第二恒等式" /> を $\varphi=u$、$\psi=G$ で使います。$R>0$ で $(\nabla^2+k^2)G=0$ （$G$ は $R$ のみの関数で、$\nabla^2G=R^{-1}(RG)''=-k^2G$）なので、$\int_{V_\epsilon}(u\nabla^2G-G\nabla^2u)\,dV=\int_{V_\epsilon}(-k^2uG+k^2Gu)\,dV=0$ です。したがって
$$
0=\oint_{\partial V}\bigl(u\,\partial_nG-G\,\partial_nu\bigr)dS+\oint_{|\boldsymbol{x}-P|=\epsilon}\bigl(u\,\partial_\nu G-G\,\partial_\nu u\bigr)dS ,
$$
$\nu$ は $V_\epsilon$ からみた外向き、つまり $P$ に向かう向きで、$\partial_\nu=-\partial_R$ です。小球面上で $G=e^{ik\epsilon}/\epsilon$、$\partial_\nu G=(\epsilon^{-2}-ik\epsilon^{-1})e^{ik\epsilon}$ ですから、面積 $4\pi\epsilon^2$ と $u$ の連続性から $\oint u\,\partial_\nu G\,dS\to4\pi u(P)$、$\bigl|\oint G\,\partial_\nu u\,dS\bigr|\le4\pi\epsilon^2\cdot\epsilon^{-1}\sup|\nabla u|\to0$ です。$\epsilon\to0$ として整理すれば公式を得ます。
</Proof>

この公式は、穴のあいたスクリーンの向こう側の場を、開口上の値から求めるのに使えます。ただし開口上の本当の場の値は分からないので、次の 3 つの仮定で「開口では入射波そのまま、影では 0」と見なして計算を閉じるのがねらいです。**キルヒホッフの近似**として、(i) 開口上では $u$ と $\partial_nu$ が障害物がないときの入射波の値に等しい、(ii) スクリーンの影の面では $u=\partial_nu=0$、(iii) 無限遠の半球面からの寄与は $0$、の 3 つを仮定します（厳密には成り立たない近似です。開口の縁の近くで場が乱れるためで、開口の大きさが波長より十分大きいときに良い近似になります）。スクリーンを $z=0$ の平面にし、$z>0$ 側に観測点 $P$、$z<0$ 側から $+z$ 方向に垂直入射する平面波 $u_{\mathrm{inc}}=u_0e^{ikz}$ を考えると、開口上で $\partial_nu=-\partial_zu=-iku_0$（$V$ は $z>0$ の領域で外向き法線は $-\boldsymbol{e}_z$）、$\partial_nG=-\partial_zG\approx ik\cos\chi\,G$（$\chi$ は法線と開口から $P$ への方向のなす角、$kR\gg1$ で $1/R$ の項を落とした）なので、

$$
u(P)\approx-\frac{iku_0}{4\pi}\int_A\frac{e^{ikR}}{R}\bigl(1+\cos\chi\bigr)dS .
$$

### 6.2. フラウンホーファー回折

遠方で見る回折を**フラウンホーファー回折**といい、ここでは光路の差を一次の項だけで近似して積分を簡単にします。観測点が開口から十分遠く、方向 $\hat{\boldsymbol{r}}$ に距離 $r$ にあるとします。開口上の点 $\boldsymbol{x}'$ からの距離は
$$
R=|r\hat{\boldsymbol{r}}-\boldsymbol{x}'|=r-\hat{\boldsymbol{r}}\cdot\boldsymbol{x}'+\frac{|\boldsymbol{x}'|^2-(\hat{\boldsymbol{r}}\cdot\boldsymbol{x}')^2}{2r}+\cdots
$$
と展開できます。開口の大きさを $D$ として $kD^2/(2r)\ll1$（**フラウンホーファー条件**）なら第 3 項以降の位相への寄与は無視でき、小角度では $1+\cos\chi\approx2$、$1/R\approx1/r$ としてよいので、

<Proposition id="prop-fraunhofer" title="フラウンホーファー回折">
キルヒホッフの近似とフラウンホーファー条件のもとで
$$
u(P)\approx-\frac{iu_0}{\lambda}\frac{e^{ikr}}{r}\int_Ae^{-ik\,\hat{\boldsymbol{r}}\cdot\boldsymbol{x}'}\,dS' ,\qquad\lambda=\frac{2\pi}{k}
$$
である。すなわち回折像の振幅は、開口の形の（空間）フーリエ変換で与えられる（<Ref to="mathematics/fourier-analysis/fourier-transform#def-fourier-transform" text="フーリエ変換の定義" /> を参照）。
</Proposition>

<Proof of="prop-fraunhofer">
上の近似を $u(P)\approx-\frac{iku_0}{4\pi}\int_A\frac{e^{ikR}}{R}(1+\cos\chi)\,dS$ に代入します。$e^{ikR}\approx e^{ikr}e^{-ik\hat{\boldsymbol{r}}\cdot\boldsymbol{x}'}$、$(1+\cos\chi)/R\approx2/r$ で、係数は $-\frac{ik}{4\pi}\cdot2=-\frac{ik}{2\pi}=-\frac i\lambda$ です。
</Proof>

<Example id="ex-single-slit" title="単スリットと二重スリット">
$x$ 方向の幅が $a$、$y$ 方向に十分長いスリットを考え、観測方向が $x$ 軸となす角を使って $\hat{\boldsymbol{r}}\cdot\boldsymbol{x}'=x'\sin\theta$ とします。<Ref to="prop-fraunhofer" /> の積分は
$$
\int_{-a/2}^{a/2}e^{-ikx'\sin\theta}dx'=\frac{2\sin(ka\sin\theta/2)}{k\sin\theta}=a\,\frac{\sin\beta}{\beta},\qquad\beta=\frac{\pi a\sin\theta}{\lambda}
$$
で、強度は $I(\theta)=I(0)\Bigl(\dfrac{\sin\beta}{\beta}\Bigr)^2$ です。暗線は $a\sin\theta=m\lambda$（$m=\pm1,\pm2,\dots$）に現れ、中央の明部の角幅は $2\lambda/a$ です。たとえば $a=0.05\,\mathrm{mm}$、$\lambda=500\,\mathrm{nm}$、スクリーンまで $1\,\mathrm{m}$ なら、中央の明部の幅は $2\lambda L/a=2\,\mathrm{cm}$ で、フラウンホーファー条件 $a^2/(\lambda L)\approx0.005\ll1$ も満たされています。スリットを狭くするほど像が広がります。

中心間隔 $d$ の同じ 2 本のスリットでは、積分が 2 つの区間 $[\pm d/2-a/2,\ \pm d/2+a/2]$ に分かれ、それぞれ $a\,\mathrm{sinc}\,\beta$ に位相 $e^{\mp ikd\sin\theta/2}$ が付くので、和は $2a\,\frac{\sin\beta}{\beta}\cos\bigl(\frac{\pi d\sin\theta}{\lambda}\bigr)$ になり、
$$
I(\theta)=I(0)\cos^2\Bigl(\frac{\pi d\sin\theta}{\lambda}\Bigr)\Bigl(\frac{\sin\beta}{\beta}\Bigr)^2
$$
です。明線は $d\sin\theta=m\lambda$、スクリーン上の縞の間隔は小角度で $\lambda L/d$ です（$\lambda=550\,\mathrm{nm}$、$d=0.5\,\mathrm{mm}$、$L=1\,\mathrm{m}$ で $1.1\,\mathrm{mm}$）。$\cos^2$ が**干渉**縞、$\mathrm{sinc}^2$ が各スリットの**回折**による包絡線です。
</Example>

前半ではスリットの幅 $a$ を広げ、後半では間隔 $d$ を広げます。破線の包絡線（回折）は $a$ だけで決まり、その内側の細かい縞（干渉）は $d$ だけで決まることを見てください。

<Figure caption="二重スリットの強度 I(θ)。a を広げると包絡線が狭まり、d を広げると縞が細かくなる（λ = 500 nm）。">
<img src="/figures/double-slit-intensity.gif" alt="横軸 sinθ、縦軸強度の二重スリットの強度分布。前半はスリット幅 a を 0.05 mm から 0.20 mm に広げると破線の包絡線が狭くなり、後半は間隔 d を 0.10 mm から 0.50 mm に広げると包絡線は変わらず内側の縞が細かくなる。" loading="lazy" />
</Figure>

### 6.3. 薄膜の干渉

2 つの波が干渉するとき、複素振幅を $A_1=\sqrt{I_1}$、$A_2=\sqrt{I_2}e^{i\delta}$ とすると、合成強度は $|A_1+A_2|^2=I_1+I_2+2\sqrt{I_1I_2}\cos\delta$ です（$\delta$ は位相差）。

<Example id="ex-ar-coating" title="反射防止膜の条件">
屈折率 $n_2$ のガラスの上に、屈折率 $n_c$（$1<n_c<n_2$）、厚さ $d$ の薄膜をつけ、空気（$n_1=1$）側から垂直入射させます。膜の表面での反射の振幅比は <Ref to="thm-fresnel" /> の垂直入射の場合で $r_1=\dfrac{n_1-n_c}{n_1+n_c}$、ガラスとの境界では $r_2=\dfrac{n_c-n_2}{n_c+n_2}$ で、$n_1<n_c<n_2$ なので両方とも負（同じ符号）です。膜の中を往復した波は光路 $2n_cd$ ぶん余分に進むので位相が $2kn_cd=\dfrac{4\pi n_cd}{\lambda}$ だけ遅れます。2 つの反射波が逆位相で打ち消し合うのは、遅れが $\pi$ のとき、すなわち
$$
d=\frac{\lambda}{4n_c}
$$
（**4 分の 1 波長膜**）のときで、振幅も等しい（$r_1=r_2$）なら完全に消えます。$r_1=r_2$ を整理すると $(n_c-n_1)(n_2+n_c)=(n_2-n_c)(n_c+n_1)$ から $n_c^2=n_1n_2$ です。膜内の多重反射まで含めた反射振幅は $\dfrac{r_1+r_2e^{i\phi}}{1+r_1r_2e^{i\phi}}$（$\phi$ は往復の位相遅れ）で、$\phi=\pi$ では $\dfrac{r_1-r_2}{1-r_1r_2}$ となり、$r_1=r_2$ で厳密に $0$ になります。以下の数値は 2 つの反射波だけを足した近似です。MgF$_2$（$n_c=1.38$）膜をガラス（$n_2=1.52$）につけた場合は、理想値 $\sqrt{1.52}=1.23$ より少し大きく、残りの反射は約 $1.2\%$（無反射の $4.3\%$ から減る）です（演習）。
</Example>

---

## 7. 演習

<Exercise id="exr-brewster-water" difficulty="易">
空気（$n_1=1$）から水（$n_2=1.33$）への入射で、ブリュースター角と、そのときの屈折角を求めよ。また水から空気へ向かうときの臨界角を求めよ。
<Solution>
<Ref to="cor-brewster" /> より $\theta_B=\arctan1.33\approx53.1^\circ$。屈折角は $\theta_t=90^\circ-\theta_B\approx36.9^\circ$（$\sin\theta_t=\sin\theta_B/1.33=0.7996/1.33\approx0.601$ と整合）。臨界角は $\arcsin(1/1.33)\approx48.8^\circ$ です。
</Solution>
</Exercise>

<Exercise id="exr-skin-depth-sea" difficulty="標準">
海水（$\sigma\approx4\,\mathrm{S/m}$、$\varepsilon\approx80\varepsilon_0$、$\mu\approx\mu_0$）の中を伝わる $1\,\mathrm{kHz}$ の電磁波の表皮厚さを求め、この周波数が良導体の近似の範囲にあることを確かめよ。潜水艦に向けた通信に超長波（VLF）以下の低い周波数が使われる理由を述べよ。
<Solution>
$\omega=2\pi\times10^3$ で $\sigma/(\varepsilon\omega)=4/(80\times8.85\times10^{-12}\times6283)\approx9\times10^5\gg1$ なので <Ref to="prop-conducting-wave" /> の近似が使えます。$\delta=\sqrt{2/(\mu_0\sigma\omega)}=\sqrt{2/(1.257\times10^{-6}\times4\times6283)}\approx8.0\,\mathrm{m}$ です。$\delta\propto\omega^{-1/2}$ なので周波数が高いと数 cm〜数 m で減衰し、潜水艦まで届きません。水中に届かせるには数十 Hz〜数 kHz の低い周波数が必要です（そのかわり情報量は小さい）。
</Solution>
</Exercise>

<Exercise id="exr-ar-numbers" difficulty="標準">
波長 $550\,\mathrm{nm}$ の光に対して、MgF$_2$（$n_c=1.38$）の反射防止膜を $n_2=1.52$ のガラスにつけるとき、膜厚と、多重反射を無視したときの反射率を求めよ。膜がない場合と比較せよ。
<Solution>
$d=\lambda/(4n_c)=550/(4\times1.38)\approx99.6\,\mathrm{nm}$。$r_1=\dfrac{1-1.38}{2.38}=-0.1597$、$r_2=\dfrac{1.38-1.52}{2.90}=-0.0483$。膜内を往復した波は位相が $\pi$ 遅れるので合成の反射振幅は $r_1+r_2e^{i\pi}=-0.1597+0.0483=-0.1114$ で、反射率は $0.1114^2\approx1.2\%$ です。膜がなければ $((1.52-1)/(2.52))^2=4.3\%$ です。完全な無反射にならないのは $r_1\ne r_2$（$n_c\ne\sqrt{n_2}$）のためです。
</Solution>
</Exercise>

<Exercise id="exr-guide-cutoff" difficulty="難">
導波管の壁の間隔が $a$ の一次元的な場合として、2 枚の平行な完全導体板 $x=0$、$x=a$ の間を $z$ 方向に伝わる TE モード（$E_y$ のみ）を考える。$E_y=\sin(m\pi x/a)e^{i(kz-\omega t)}$ が解になる条件を求め、波が 2 つの平面波の重ね合わせであることを示して、遮断周波数の意味を説明せよ。
<Solution>
波動方程式から $-(m\pi/a)^2-k^2+\omega^2/c^2=0$ で $k^2=\omega^2/c^2-(m\pi/a)^2$、壁（$x=0,a$）で $E_y=0$ は満たされています。$\sin(m\pi x/a)=\dfrac{1}{2i}\bigl(e^{im\pi x/a}-e^{-im\pi x/a}\bigr)$ なので、場は波数ベクトル $(\pm m\pi/a,0,k)$ をもつ 2 つの平面波の差で、板の間を斜めに反射しながら進む波です。波数ベクトルの大きさは $\omega/c$ で、$z$ 成分 $k=\sqrt{\omega^2/c^2-(m\pi/a)^2}$ が実数になるには、$x$ 成分 $m\pi/a\le\omega/c$ が必要です。$\omega=cm\pi/a$ は波が板に垂直に往復するだけで $z$ 方向に進まない境目で、これが遮断周波数です。
</Solution>
</Exercise>

---

## 参考文献

- J. D. Jackson『電磁気学（原書第 3 版）』（Classical Electrodynamics, 3rd ed., Wiley, 1999）— 第 7 章（平面波、反射と屈折、分散）、第 8 章（導波管）、第 10 章（回折）。
- M. Born, E. Wolf, *Principles of Optics*, 7th ed., Cambridge University Press, 1999 — 第 1 章（電磁理論）、第 8 章（回折の基礎理論）。
- E. Hecht, *Optics*, 5th ed., Pearson, 2017 — 反射・屈折、干渉、回折の章。
- D. J. Griffiths, *Introduction to Electrodynamics*, 4th ed., Cambridge University Press, 2017 — 第 9 章（電磁波）。
- L. D. Landau, E. M. Lifshitz, *Electrodynamics of Continuous Media*, 2nd ed., Pergamon Press, 1984 — 電磁波の伝播の章。
