携帯電話で通話しているとき、あなたの声は電気信号に変換され、最終的にアンテナから電磁波として空間に飛び出していきます。しかし、なぜ金属の棒に交流電流を流すだけで、目に見えない波が空間を伝わっていくのでしょうか？導線を流れる電流が外部に電磁波を「生み出す」メカニズムは、一見すると魔法のようにも見えます。
この疑問に答えるのがアンテナの放射メカニズムです。鍵となるのは、加速する電荷は必ず電磁波を放射するという物理法則です。この原理を理解すれば、アンテナ設計の根幹が見えてきます。
アンテナの放射メカニズムを理解すると、以下のような幅広い分野に応用できます。
- 無線通信システムの設計: 携帯基地局、Wi-Fiルーター、衛星通信など、あらゆる無線システムでアンテナの放射特性を最適化する基盤となります
- レーダー工学: 電磁波の放射と散乱の理解は、航空管制レーダーや気象レーダーの性能予測に不可欠です
- EMC（電磁両立性）: 意図しない電磁波放射（ノイズ）の原因を理解し、電子機器の設計品質を高めることができます
本記事の内容
- 加速する電荷がなぜ電磁波を放射するか — 直感的理解
- マクスウェル方程式と電磁波放射の関係
- 遅延ポテンシャルの導出と物理的意味
- 微小ダイポール（ヘルツダイポール）の放射界の導出
- 放射パターンの数学的表現
- Pythonによる放射パターンの可視化と考察
前提知識
この記事を読む前に、以下の記事を読んでおくと理解が深まります。
加速する電荷はなぜ電磁波を放射するのか
静止した電荷 — 放射はない
まず、空間に1つの電荷が静止している状況を考えましょう。静止した電荷は、周囲にクーロンの法則に従う静電場 $\bm{E}$ を作ります。この電場は時間的に変化しないので、磁界 $\bm{H}$ は生じません。エネルギーの流れを表すポインティングベクトル $\bm{S} = \bm{E} \times \bm{H}$ もゼロです。つまり、静止した電荷からはエネルギーが外に飛び出していくことはありません。
等速度で動く電荷 — やはり放射はない
次に、電荷が一定の速度 $\bm{v}$ で動いている場合を考えます。等速直線運動する電荷は電流を作り、電場と磁場の両方を生じます。しかし、等速運動する電荷を「一緒に動く座標系」から見れば、それは静止した電荷に過ぎません。特殊相対論の原理により、この座標系で放射がないなら、元の座標系でも放射はありません（ローレンツ変換でポインティングベクトルの遠方成分がゼロのままであることを示せます）。
加速する電荷 — 電磁波が放射される
状況が一変するのは、電荷が加速（速度が変化）するときです。加速する電荷の周囲では、電場の構造が「歪み」を起こします。イメージとしては、電荷が急に動きを変えると、遠方の電場はまだ「古い情報」を持っています。この「新しい位置の情報」と「古い位置の情報」の境界に、電場の急激な変化（キンク）が生じます。このキンクこそが電磁波の正体です。
もう少し具体的に説明しましょう。電荷が加速すると、電場の変化 $\partial \bm{E}/\partial t$ が生じます。マクスウェル方程式のアンペール-マクスウェルの法則によれば、この時間変化する電場は磁場を生み出します。すると今度は、時間変化する磁場がファラデーの法則によって電場を生み出します。この相互作用が自己維持的に伝播していく — これが電磁波です。
つまり、加速する電荷は「電場の変化の種」を空間にまき、その種が電場と磁場の相互誘導によって波として伝播していくのです。定量的には、放射される電力は電荷の加速度の2乗に比例します（ラーモアの公式）。
$$ P = \frac{q^2 a^2}{6\pi \epsilon_0 c^3} $$
ここで $q$ は電荷、$a$ は加速度の大きさ、$\epsilon_0$ は真空の誘電率、$c$ は光速です。加速度がゼロなら放射もゼロ — これが「等速運動では放射しない」ことの定量的な裏付けです。
アンテナにおける加速電荷
アンテナに交流電流 $I(t) = I_0 \cos(\omega t)$ を流すと、導線内の自由電子は正弦的に振動します。振動運動は常に加速・減速を伴うため、電子は絶えず加速状態にあります。特に、電流が最大値と最小値の間を行き来する際に加速度が最大となり、電磁波の放射が生じます。
このように、アンテナの放射メカニズムの本質は「加速する電荷」にあります。交流電流が流れる限り、電荷は常に加速しているため、継続的に電磁波が放射されるのです。
ここまでで、放射の物理的な直感を得ました。次に、この現象を定量的に記述するマクスウェル方程式との関係を見ていきましょう。
マクスウェル方程式と電磁波放射
マクスウェル方程式の4本
電磁波の放射を定量的に扱うためには、電磁気学の基本法則であるマクスウェル方程式に立ち返る必要があります。真空中（自由空間中）のマクスウェル方程式は次の4本です。
$$ \nabla \cdot \bm{E} = \frac{\rho}{\epsilon_0} $$
$$ \nabla \cdot \bm{B} = 0 $$
$$ \nabla \times \bm{E} = -\frac{\partial \bm{B}}{\partial t} $$
$$ \nabla \times \bm{B} = \mu_0 \bm{J} + \mu_0 \epsilon_0 \frac{\partial \bm{E}}{\partial t} $$
ここで $\rho$ は電荷密度、$\bm{J}$ は電流密度、$\mu_0$ は真空の透磁率です。
放射の源泉 — 変位電流の役割
4番目の式（アンペール-マクスウェルの法則）に注目してください。右辺の第2項 $\mu_0 \epsilon_0 \partial \bm{E}/\partial t$ はマクスウェルが理論的に追加した変位電流です。この項がなければ、電場と磁場は独立に存在するだけで、波として伝播することはできません。
変位電流の存在により、次のようなフィードバックループが成立します。
- 加速する電荷（時間変化する電流 $\bm{J}$）が磁場 $\bm{B}$ を作る
- 時間変化する磁場 $\partial \bm{B}/\partial t$ が電場 $\bm{E}$ を誘導する（ファラデーの法則）
- 時間変化する電場 $\partial \bm{E}/\partial t$ が再び磁場を誘導する（変位電流）
- 2と3が繰り返され、電磁波が空間を伝播する
このフィードバックが自立的に維持されるためには、電場と磁場が互いに直交し、かつ伝播方向に垂直でなければなりません。これがまさに横波としての電磁波の性質です。
波動方程式の導出
マクスウェル方程式から電磁波の波動方程式を導出しましょう。源（電荷・電流）がない領域 $(\rho = 0, \bm{J} = 0)$ を考えます。
ファラデーの法則の回転を取ります。
$$ \nabla \times (\nabla \times \bm{E}) = -\frac{\partial}{\partial t}(\nabla \times \bm{B}) $$
左辺にベクトル恒等式 $\nabla \times (\nabla \times \bm{E}) = \nabla(\nabla \cdot \bm{E}) – \nabla^2 \bm{E}$ を適用します。源がない領域では $\nabla \cdot \bm{E} = 0$ なので、左辺は $-\nabla^2 \bm{E}$ となります。
右辺にアンペール-マクスウェルの法則（$\bm{J} = 0$ として）を代入すると、
$$ -\nabla^2 \bm{E} = -\mu_0 \epsilon_0 \frac{\partial^2 \bm{E}}{\partial t^2} $$
整理すると、電場に対する波動方程式が得られます。
$$ \nabla^2 \bm{E} = \mu_0 \epsilon_0 \frac{\partial^2 \bm{E}}{\partial t^2} $$
この方程式の解は、速度 $c = 1/\sqrt{\mu_0 \epsilon_0} \approx 3 \times 10^8$ m/s で伝播する波です。これが光速であり、電磁波が光と同じ性質を持つことを示しています。
同様の波動方程式は磁場 $\bm{B}$ についても成立します。
源がある場合 — ポテンシャルの導入
上の波動方程式は源がない領域で成り立つものでした。アンテナのように電流源がある場合には、源を含む形で電磁場を記述する必要があります。そこで登場するのが電磁ポテンシャルです。
磁場のダイバージェンスがゼロ（$\nabla \cdot \bm{B} = 0$）であることから、磁場はベクトルポテンシャル $\bm{A}$ の回転として表せます。
$$ \bm{B} = \nabla \times \bm{A} $$
これをファラデーの法則に代入すると、
$$ \nabla \times \bm{E} = -\frac{\partial}{\partial t}(\nabla \times \bm{A}) $$
$$ \nabla \times \left(\bm{E} + \frac{\partial \bm{A}}{\partial t}\right) = 0 $$
回転がゼロのベクトル場はスカラーポテンシャル $\phi$ の勾配で表せるので、
$$ \bm{E} + \frac{\partial \bm{A}}{\partial t} = -\nabla \phi $$
すなわち、
$$ \bm{E} = -\nabla \phi – \frac{\partial \bm{A}}{\partial t} $$
これにより、電磁場は2つのポテンシャル $(\phi, \bm{A})$ で完全に記述できます。ポテンシャルを使うメリットは、4本のマクスウェル方程式を2本のポテンシャルに関する方程式にまとめられることです。
ここまでで、マクスウェル方程式が電磁波の存在を予言し、ポテンシャルによって源を含む問題を扱える形にできることがわかりました。次に、このポテンシャルが源からどのように決まるかを示す「遅延ポテンシャル」を導出します。
遅延ポテンシャルの導出
ローレンツゲージ条件
ポテンシャル $(\phi, \bm{A})$ にはゲージ自由度（ある種の任意性）があります。計算を簡潔にするために、ローレンツゲージと呼ばれる条件を課します。
$$ \nabla \cdot \bm{A} + \mu_0 \epsilon_0 \frac{\partial \phi}{\partial t} = 0 $$
この条件のもとで、マクスウェル方程式はポテンシャルに対する2本の独立な波動方程式に帰着します。
ベクトルポテンシャル $\bm{A}$ に対して、
$$ \nabla^2 \bm{A} – \mu_0 \epsilon_0 \frac{\partial^2 \bm{A}}{\partial t^2} = -\mu_0 \bm{J} $$
スカラーポテンシャル $\phi$ に対して、
$$ \nabla^2 \phi – \mu_0 \epsilon_0 \frac{\partial^2 \phi}{\partial t^2} = -\frac{\rho}{\epsilon_0} $$
これらは非同次波動方程式であり、右辺の源（電流密度 $\bm{J}$ と電荷密度 $\rho$）がポテンシャルを生成します。
遅延ポテンシャルの物理的意味
上の波動方程式の解は、静電ポテンシャルや静磁ポテンシャルの公式に「遅延時間」を導入したものになります。直感的に言えば、ある場所の電流が作るポテンシャルは、光速で伝わるために遅れて観測点に届く — これが「遅延（retarded）」の意味です。
源の位置を $\bm{r}’$、観測点を $\bm{r}$ とすると、両者間の距離は $R = |\bm{r} – \bm{r}’|$ です。情報が光速 $c$ で伝わるため、観測点での「現在の」ポテンシャルは、源の「過去の」状態（時刻 $t – R/c$ での状態）によって決まります。
ベクトルポテンシャルの遅延ポテンシャルは、次のように書けます。
$$ \bm{A}(\bm{r}, t) = \frac{\mu_0}{4\pi} \int_V \frac{\bm{J}(\bm{r}’, t – R/c)}{R} \, dV’ $$
スカラーポテンシャルについても同様に、
$$ \phi(\bm{r}, t) = \frac{1}{4\pi \epsilon_0} \int_V \frac{\rho(\bm{r}’, t – R/c)}{R} \, dV’ $$
これらの式は、空間の各点にある源の寄与を、適切な遅延時間を考慮して足し合わせたものです。
時間調和場での遅延ポテンシャル
アンテナ工学では、電流が正弦波的に振動する時間調和場（角周波数 $\omega$）を扱うことが一般的です。電流密度を $\bm{J}(\bm{r}’, t) = \bm{J}(\bm{r}’) e^{j\omega t}$ と書くと、遅延時間 $t – R/c$ による位相遅れは $e^{-jkR}$（$k = \omega/c = 2\pi/\lambda$）として表されます。
$$ \bm{A}(\bm{r}) = \frac{\mu_0}{4\pi} \int_V \frac{\bm{J}(\bm{r}’)}{R} e^{-jkR} \, dV’ $$
この式は、源の電流分布がわかれば、空間の任意の点でのベクトルポテンシャルを計算できることを意味しています。そしてベクトルポテンシャルから $\bm{B} = \nabla \times \bm{A}$ と $\bm{E} = -\nabla \phi – \partial \bm{A}/\partial t$ を使って電磁場を求めることができます。
遅延ポテンシャルは、源と観測点の間の「情報伝達の遅れ」を数学的に表現したものです。電磁波が有限の速度（光速）で伝わるという相対論的な効果を、古典電磁気学の枠組みで正確に取り込んでいます。
ここまでで、電流分布からポテンシャルを経由して電磁場を計算する方法が整いました。いよいよ、最も基本的なアンテナモデルである「微小ダイポール」にこの方法を適用して、放射界を具体的に導出しましょう。
微小ダイポール（ヘルツダイポール）の放射界
モデルの設定
最も基本的なアンテナモデルとして、微小ダイポール（ヘルツダイポール）を考えます。これは、原点に置かれた長さ $dl$（$dl \ll \lambda$）の微小な電流素で、z軸方向に一様な電流 $I_0$ が流れているとします。
このモデルが「基本」と呼ばれるのは、任意の電流分布を持つ複雑なアンテナでも、微小ダイポールの重ね合わせ（積分）として放射界を計算できるからです。ちょうど、複雑な関数をテイラー展開で基本的な項の和として表現するのと同じ発想です。
ベクトルポテンシャルの計算
微小ダイポールの電流密度は、原点に集中したz方向の電流 $I_0 dl$ です。遅延ポテンシャルの公式に代入すると、
$$ \bm{A} = \hat{z} \frac{\mu_0 I_0 dl}{4\pi r} e^{-jkr} $$
ここで $r$ は原点から観測点までの距離です。因子 $1/r$ は球面波としての振幅の減衰を、$e^{-jkr}$ は位相の遅れを表しています。
球面座標への変換
放射界を議論するには球面座標 $(r, \theta, \phi)$ が自然です。z方向の単位ベクトル $\hat{z}$ を球面座標の単位ベクトルで表すと、
$$ \hat{z} = \hat{r}\cos\theta – \hat{\theta}\sin\theta $$
したがって、ベクトルポテンシャルの各成分は次のようになります。
$$ A_r = \frac{\mu_0 I_0 dl}{4\pi r} e^{-jkr} \cos\theta $$
$$ A_\theta = -\frac{\mu_0 I_0 dl}{4\pi r} e^{-jkr} \sin\theta $$
$$ A_\phi = 0 $$
$A_\phi = 0$ はダイポールがz軸方向であるため、問題が $\phi$ に関して対称であることを反映しています。
磁界 $\bm{H}$ の計算
磁界は $\bm{H} = \frac{1}{\mu_0} \nabla \times \bm{A}$ から求めます。軸対称性から $H_r = 0$、$H_\theta = 0$ であり、非ゼロの成分は $H_\phi$ のみです。
球面座標での回転の $\phi$ 成分を計算すると、
$$ H_\phi = \frac{1}{\mu_0} \left[ \frac{1}{r}\frac{\partial(rA_\theta)}{\partial r} – \frac{1}{r}\frac{\partial A_r}{\partial \theta} \right] $$
$rA_\theta$ の $r$ による微分を計算します。$rA_\theta = -\frac{\mu_0 I_0 dl}{4\pi} e^{-jkr} \sin\theta$ なので、
$$ \frac{\partial(rA_\theta)}{\partial r} = \frac{jk\mu_0 I_0 dl}{4\pi} e^{-jkr} \sin\theta $$
$A_r$ の $\theta$ による微分を計算します。
$$ \frac{\partial A_r}{\partial \theta} = -\frac{\mu_0 I_0 dl}{4\pi r} e^{-jkr} \sin\theta $$
これらを代入して整理すると、磁界の完全な表現が得られます。
$$ H_\phi = \frac{I_0 dl}{4\pi} \sin\theta \left( \frac{jk}{r} + \frac{1}{r^2} \right) e^{-jkr} $$
電界 $\bm{E}$ の計算
電界はマクスウェル方程式 $\nabla \times \bm{H} = j\omega \epsilon_0 \bm{E}$（時間調和場）から求めます。
$$ \bm{E} = \frac{1}{j\omega \epsilon_0} \nabla \times \bm{H} $$
$H_\phi$ のみが非ゼロであるため、回転を計算すると $E_r$ と $E_\theta$ の2成分が得られます。
$E_r$ 成分について、$\nabla \times \bm{H}$ の $r$ 成分は $\frac{1}{r\sin\theta}\frac{\partial(\sin\theta \cdot H_\phi)}{\partial \theta}$ で与えられます。計算を実行すると、
$$ E_r = \frac{I_0 dl}{4\pi} \cdot 2\cos\theta \left( \frac{1}{j\omega \epsilon_0 r^3} + \frac{k}{\omega \epsilon_0 r^2} \right) e^{-jkr} $$
$\eta = \sqrt{\mu_0/\epsilon_0}$ （自由空間の波動インピーダンス、約377 $\Omega$）と $k = \omega/c = \omega\sqrt{\mu_0\epsilon_0}$ の関係を使って整理すると、
$$ E_r = \frac{\eta I_0 dl}{4\pi} \cdot 2\cos\theta \left( \frac{1}{(jkr)^2} – \frac{j}{(kr)^3} \right) k \, e^{-jkr} $$
同様に $E_\theta$ 成分も計算すると、
$$ E_\theta = \frac{\eta I_0 dl k}{4\pi} \sin\theta \left( \frac{j}{kr} + \frac{1}{(kr)^2} – \frac{j}{(kr)^3} \right) e^{-jkr} $$
近傍界と遠方界
得られた電磁場の表現には、$1/r$、$1/r^2$、$1/r^3$ の項が含まれています。これらは物理的に異なる性質を持ちます。
近傍界（$kr \ll 1$、すなわち $r \ll \lambda$）: $1/r^3$ の項が支配的です。これは静的な電気双極子の電場と同じ形をしており、誘導界と呼ばれます。エネルギーは源の近傍で蓄積と放出を繰り返し、遠方には伝播しません。
遠方界（$kr \gg 1$、すなわち $r \gg \lambda$）: $1/r$ の項が支配的になります。これが放射界であり、エネルギーが源から離れて永続的に伝播していく成分です。
放射界のみを取り出すと（$1/r^2$ 以降の項を無視）、
$$ E_\theta = j \frac{\eta k I_0 dl}{4\pi} \frac{\sin\theta}{r} e^{-jkr} $$
$$ H_\phi = j \frac{k I_0 dl}{4\pi} \frac{\sin\theta}{r} e^{-jkr} $$
$$ E_r \approx 0 \quad (\text{遠方界では無視}) $$
遠方界ではいくつかの重要な特徴があります。
- $\bm{E}$ と $\bm{H}$ は互いに直交し、かつ伝播方向 $\hat{r}$ にも直交する（横波）
- 電場と磁場の比は波動インピーダンス: $E_\theta / H_\phi = \eta \approx 377 \; \Omega$
- 振幅は $1/r$ に比例して減衰する（球面波の性質）
- 角度依存性は $\sin\theta$ — これが放射パターンを決定する
ここまでで、微小ダイポールの放射界を導出しました。次に、この $\sin\theta$ の依存性をもとに、放射パターンを定量的に議論しましょう。
放射パターンの導出
放射パターンとは
放射パターン（radiation pattern）とは、アンテナからの放射強度の角度分布を表したものです。放射パワー密度（ポインティングベクトルの大きさ）が方向 $(\theta, \phi)$ によってどう変化するかを示します。
直感的には、アンテナを中心に置いた大きな球面上で、各方向に飛んでいく電磁波の「強さ」を地図のように表したものと考えることができます。懐中電灯のように特定の方向に集中するアンテナもあれば、白熱電球のように全方向にほぼ均等に放射するアンテナもあります。
ポインティングベクトルと放射強度
時間平均ポインティングベクトルは、遠方界の電磁場から次のように計算されます。
$$ \bm{S}_{\text{av}} = \frac{1}{2} \text{Re}[\bm{E} \times \bm{H}^*] $$
微小ダイポールの遠方界では $E_\theta$ と $H_\phi$ のみが非ゼロなので、
$$ S_r = \frac{1}{2} \text{Re}[E_\theta H_\phi^*] $$
先に求めた遠方界の式を代入します。
$$ S_r = \frac{1}{2} \frac{\eta k^2 (I_0 dl)^2}{(4\pi)^2} \frac{\sin^2\theta}{r^2} $$
この式から、放射パワー密度は $\sin^2\theta$ に比例することがわかります。
放射強度
放射強度（radiation intensity） $U(\theta, \phi)$ は、ポインティングベクトルに $r^2$ を掛けたもので定義されます。
$$ U(\theta, \phi) = r^2 S_r = \frac{\eta k^2 (I_0 dl)^2}{2(4\pi)^2} \sin^2\theta $$
放射強度は距離 $r$ に依存しない量であり、純粋に方向のみの関数です。これを最大値で正規化したものが正規化放射パターンです。
$$ F(\theta) = \frac{U(\theta)}{U_{\text{max}}} = \sin^2\theta $$
$\sin^2\theta$ は $\theta = 90°$（赤道面、ダイポールに垂直な方向）で最大値1をとり、$\theta = 0°$ と $\theta = 180°$（ダイポールの軸方向）でゼロになります。
放射パターンの特徴
微小ダイポールの放射パターン $\sin^2\theta$ の物理的意味を考えましょう。
なぜダイポール軸方向に放射しないのか？ ダイポールの軸方向（$\theta = 0°$）から見ると、電流素は「点」に見えます。電流の振動方向（z方向）と観測方向（z方向）が平行なので、観測方向に垂直な電場成分（横波条件を満たす成分）が存在しません。電磁波は横波ですから、この方向には放射できないのです。
なぜ赤道面で最大なのか？ $\theta = 90°$ の方向から見ると、電流の振動が最もよく「見えます」。電流の振動方向と観測方向が直交しているため、横波としての電場成分が最大になります。
全放射電力と放射抵抗
全放射電力 $P_{\text{rad}}$ は、ポインティングベクトルを閉曲面で積分して求めます。
$$ P_{\text{rad}} = \int_0^{2\pi} \int_0^{\pi} S_r \cdot r^2 \sin\theta \, d\theta \, d\phi $$
$\sin^2\theta$ を含む積分を実行します。$\phi$ 積分は $2\pi$ を与え、$\theta$ 積分は $\int_0^\pi \sin^3\theta \, d\theta$ です。
$\int_0^\pi \sin^3\theta \, d\theta$ を計算するために、$\sin^3\theta = \sin\theta(1 – \cos^2\theta)$ と変形し、$u = \cos\theta$ と置換します。
$$ \int_0^\pi \sin^3\theta \, d\theta = \int_{-1}^{1} (1-u^2) \, du = \left[u – \frac{u^3}{3}\right]_{-1}^{1} = \frac{4}{3} $$
したがって、
$$ P_{\text{rad}} = \frac{\eta k^2 (I_0 dl)^2}{2(4\pi)^2} \cdot 2\pi \cdot \frac{4}{3} = \frac{\eta k^2 (I_0 dl)^2}{12\pi} $$
$k = 2\pi/\lambda$、$\eta = 120\pi \; \Omega$ を代入すると、
$$ P_{\text{rad}} = \frac{120\pi \cdot (2\pi/\lambda)^2 \cdot (I_0 dl)^2}{12\pi} = 40\pi^2 \left(\frac{dl}{\lambda}\right)^2 I_0^2 $$
放射抵抗 $R_r$ は、アンテナが放射として消費する等価的な抵抗で、$P_{\text{rad}} = \frac{1}{2} R_r I_0^2$ から定義されます。
$$ R_r = \frac{2P_{\text{rad}}}{I_0^2} = 80\pi^2 \left(\frac{dl}{\lambda}\right)^2 $$
例えば $dl = 0.01\lambda$ のとき、$R_r = 80\pi^2 \times (0.01)^2 \approx 0.079 \; \Omega$ と非常に小さな値になります。微小ダイポールの放射抵抗が小さいことは、効率的な放射のためにはある程度の長さのアンテナが必要であることを意味しています。
指向性
指向性 $D$ は、全方向に均等に放射した場合と比べて、最大放射方向にどれだけ集中しているかを示す指標です。
$$ D = \frac{U_{\text{max}}}{P_{\text{rad}} / 4\pi} = \frac{4\pi U_{\text{max}}}{P_{\text{rad}}} $$
微小ダイポールでは $U_{\text{max}}$ は $\theta = 90°$ での値で、
$$ U_{\text{max}} = \frac{\eta k^2 (I_0 dl)^2}{2(4\pi)^2} $$
これを代入すると、
$$ D = \frac{4\pi \cdot \frac{\eta k^2 (I_0 dl)^2}{2(4\pi)^2}}{\frac{\eta k^2 (I_0 dl)^2}{12\pi}} = \frac{12\pi}{2 \cdot 4\pi} = \frac{3}{2} = 1.5 $$
微小ダイポールの指向性は $D = 1.5$（約1.76 dBi）です。これは等方性アンテナの1.5倍の最大放射強度を持つことを意味します。
ここまでで、微小ダイポールの放射パターンの理論的な導出が完了しました。次に、Pythonを使ってこの放射パターンを実際に可視化し、理論で得られた結果を確認しましょう。
Pythonによる放射パターンの可視化
2D放射パターン（極座標プロット）
まず、微小ダイポールの放射パターン $F(\theta) = \sin^2\theta$ を極座標で可視化します。
import numpy as np
import matplotlib.pyplot as plt
# 角度の配列（0からπまで、十分細かく）
theta = np.linspace(0, 2 * np.pi, 361)
# 微小ダイポールの正規化放射パターン F(θ) = sin²θ
# ここでは E面パターン（電場面）を描画
# E面は phi=0 の面（xz平面）
F_theta = np.sin(theta) ** 2
# 極座標プロット
fig, axes = plt.subplots(1, 2, figsize=(14, 6), subplot_kw={"projection": "polar"})
# (a) 線形スケール
axes[0].plot(theta, F_theta, "c-", linewidth=2)
axes[0].set_title("Radiation Pattern (Linear)", fontsize=13, pad=20)
axes[0].set_rticks([0.25, 0.5, 0.75, 1.0])
axes[0].set_rlabel_position(45)
axes[0].grid(True, alpha=0.3)
# (b) dBスケール
F_dB = 10 * np.log10(np.where(F_theta > 1e-10, F_theta, 1e-10))
F_dB_clipped = np.clip(F_dB, -30, 0) # -30 dB以下はクリップ
axes[1].plot(theta, F_dB_clipped + 30, "c-", linewidth=2) # 0〜30の範囲にシフト
axes[1].set_title("Radiation Pattern (dB)", fontsize=13, pad=20)
axes[1].set_rticks([0, 10, 20, 30])
axes[1].set_yticklabels(["-30", "-20", "-10", "0"])
axes[1].set_rlabel_position(45)
axes[1].grid(True, alpha=0.3)
plt.tight_layout()
plt.savefig("dipole_radiation_pattern_2d.png", dpi=150, bbox_inches="tight")
plt.show()
上のグラフから、いくつかの重要な特徴が読み取れます。
-
ドーナツ形状（8の字）: 線形スケールで見ると、放射パターンは典型的な「8の字」型です。赤道方向（$\theta = 90°$ と $270°$）で最大値をとり、ダイポール軸方向（$\theta = 0°$ と $180°$）ではゼロになります。これは $\sin^2\theta$ の理論予測と完全に一致しています。
-
dBスケールでのパターン: dBスケールでは、軸方向のヌル（放射ゼロ）が非常に深い谷として鮮明に表れます。実際のアンテナ測定では、ヌル方向でもわずかな放射がありますが、理想的な微小ダイポールでは完全なヌルです。
-
ビーム幅: 半値角（3 dBビーム幅）は約90度と広く、微小ダイポールの放射は比較的広い範囲に広がっていることがわかります。これは指向性 $D = 1.5$ が比較的低い値であることと整合します。
3D放射パターン
次に、放射パターンの3次元表現を描画します。ドーナツ状の放射パターンを立体的に確認しましょう。
import numpy as np
import matplotlib.pyplot as plt
# 3D放射パターンの計算
theta_3d = np.linspace(0, np.pi, 91)
phi_3d = np.linspace(0, 2 * np.pi, 181)
THETA, PHI = np.meshgrid(theta_3d, phi_3d)
# 放射強度 F(θ) = sin²θ
R = np.sin(THETA) ** 2
# 直交座標に変換
X = R * np.sin(THETA) * np.cos(PHI)
Y = R * np.sin(THETA) * np.sin(PHI)
Z = R * np.cos(THETA)
fig = plt.figure(figsize=(10, 8))
ax = fig.add_subplot(111, projection="3d")
# カラーマップで放射強度を色分け
surf = ax.plot_surface(X, Y, Z, facecolor=plt.cm.cool(R), alpha=0.8,
rstride=1, cstride=1, linewidth=0)
ax.set_xlabel("X", fontsize=12)
ax.set_ylabel("Y", fontsize=12)
ax.set_zlabel("Z (Dipole axis)", fontsize=12)
ax.set_title("3D Radiation Pattern of Hertzian Dipole", fontsize=14)
ax.set_xlim([-1, 1])
ax.set_ylim([-1, 1])
ax.set_zlim([-1, 1])
ax.view_init(elev=20, azim=45)
plt.tight_layout()
plt.savefig("dipole_radiation_pattern_3d.png", dpi=150, bbox_inches="tight")
plt.show()
3D表示により、微小ダイポールの放射パターンがドーナツ（トーラス）形状であることが直感的に把握できます。
-
z軸（ダイポール軸）方向のくぼみ: ダイポールの上下方向には電磁波が放射されないため、z軸方向にパターンがつぶれています。これは横波条件から導かれた $\sin^2\theta$ の帰結です。
-
赤道面（xy平面）の膨らみ: 赤道方向に最も強く放射されています。3D表示で見ると、この膨らみがダイポールの周囲を一周する環状のパターンであることがわかります。
-
$\phi$ 方向の対称性: ダイポールがz軸に沿っているため、z軸周りの回転対称性があり、方位角 $\phi$ には依存しません。これは式中に $\phi$ が含まれないことの視覚的な確認です。
電磁場の各成分の距離依存性
放射界（$1/r$）、誘導界（$1/r^2$）、静電界（$1/r^3$）の距離依存性を比較するグラフを描きましょう。
import numpy as np
import matplotlib.pyplot as plt
# 正規化距離 r/λ
r_norm = np.linspace(0.05, 3.0, 500)
# 各項の相対的な大きさ（kr = 2π(r/λ)）
kr = 2 * np.pi * r_norm
# E_θの各項（大きさ）
radiation = 1 / kr # 放射界 (1/r)
induction = 1 / kr**2 # 誘導界 (1/r²)
electrostatic = 1 / kr**3 # 静電界 (1/r³)
fig, ax = plt.subplots(figsize=(10, 6))
ax.semilogy(r_norm, radiation, "c-", linewidth=2.5, label=r"Radiation field $\propto 1/r$")
ax.semilogy(r_norm, induction, "m--", linewidth=2.5, label=r"Induction field $\propto 1/r^2$")
ax.semilogy(r_norm, electrostatic, "y:", linewidth=2.5, label=r"Electrostatic field $\propto 1/r^3$")
# 遠方界の境界線（r = λ/(2π)、すなわち kr = 1）
ax.axvline(x=1 / (2 * np.pi), color="white", linestyle="-.", alpha=0.5, linewidth=1.5)
ax.text(1 / (2 * np.pi) + 0.02, 5, r"$r = \lambda/2\pi$", fontsize=11, color="white")
ax.set_xlabel(r"Normalized distance $r/\lambda$", fontsize=13)
ax.set_ylabel("Relative field amplitude", fontsize=13)
ax.set_title("Distance Dependence of EM Field Components", fontsize=14)
ax.legend(fontsize=12, loc="upper right")
ax.set_xlim([0.05, 3.0])
ax.set_ylim([1e-3, 1e2])
ax.grid(True, alpha=0.3, which="both")
plt.tight_layout()
plt.savefig("field_distance_dependence.png", dpi=150, bbox_inches="tight")
plt.show()
このグラフから、電磁場の「領域」の性質が明確に読み取れます。
-
近傍界（$r \ll \lambda/2\pi$）: 静電界（$1/r^3$）が圧倒的に大きく、放射界は無視できるほど小さい。この領域でのエネルギーは蓄積・放出を繰り返す反応性のエネルギーです。
-
遷移領域（$r \approx \lambda/2\pi$）: 3つの項がほぼ同じ大きさになる境界です。$r = \lambda/(2\pi)$、すなわち $kr = 1$ の点が遠方界と近傍界の境界の目安となります。
-
遠方界（$r \gg \lambda/2\pi$）: 放射界（$1/r$）が支配的になります。この領域では電磁波としてエネルギーが永続的に伝播していきます。放射界が他の項に対して十分大きくなる $r > 2\lambda$ 程度が、放射パターン測定に適した距離です。
この距離依存性の理解は、アンテナ測定の実務においても重要です。放射パターンを正確に測定するには、近傍界の影響を避けるために十分な距離を確保する必要があります。
半波長ダイポールとの比較
最後に、微小ダイポールと半波長ダイポールの放射パターンを比較しましょう。半波長ダイポールの放射パターンは $\cos(\frac{\pi}{2}\cos\theta)/\sin\theta$ で与えられます。
import numpy as np
import matplotlib.pyplot as plt
theta = np.linspace(0.001, np.pi, 1000)
# 微小ダイポール: F(θ) = sin²θ → E面パターン ∝ sinθ
F_hertz = np.sin(theta)
# 半波長ダイポール: F(θ) = cos(π/2 cosθ) / sinθ
F_half = np.cos(np.pi / 2 * np.cos(theta)) / np.sin(theta)
# 正規化（最大値で割る）
F_hertz_norm = F_hertz / np.max(F_hertz)
F_half_norm = F_half / np.max(F_half)
# 全角度（0〜2π）に拡張して極座標プロット
theta_full = np.concatenate([theta, theta + np.pi])
F_hertz_full = np.concatenate([F_hertz_norm, F_hertz_norm])
F_half_full = np.concatenate([F_half_norm, F_half_norm])
fig, ax = plt.subplots(figsize=(8, 8), subplot_kw={"projection": "polar"})
ax.plot(theta_full, F_hertz_full, "c-", linewidth=2, label="Hertzian dipole")
ax.plot(theta_full, F_half_full, "m--", linewidth=2, label=r"Half-wave dipole ($\lambda/2$)")
ax.set_title("E-plane Pattern Comparison", fontsize=14, pad=20)
ax.legend(loc="upper right", bbox_to_anchor=(1.35, 1.1), fontsize=11)
ax.set_rticks([0.25, 0.5, 0.75, 1.0])
ax.set_rlabel_position(45)
ax.grid(True, alpha=0.3)
plt.tight_layout()
plt.savefig("dipole_comparison.png", dpi=150, bbox_inches="tight")
plt.show()
微小ダイポールと半波長ダイポールの放射パターンを比較すると、以下の特徴がわかります。
-
形状の類似性: 両者とも基本的な「8の字」形状を持ち、軸方向にヌルがあります。これはダイポール型アンテナに共通の特徴です。
-
ビーム幅の違い: 半波長ダイポールの方がわずかにビームが狭く（赤道面に集中し）、指向性が高い（$D = 1.64$ vs $D = 1.50$）ことが視覚的に確認できます。
-
実用上の差: 両者の放射パターンの差は小さく、多くの実用計算では微小ダイポールの $\sin\theta$ パターンで十分な近似が得られます。しかし、利得や放射抵抗の計算では違いが重要になります。
放射メカニズムの全体像
ここまでの議論を整理して、アンテナの放射メカニズムの全体像をまとめましょう。
ステップ1: 加速する電荷
アンテナに交流電流が流れると、導線内の自由電子は正弦的に振動します。振動運動には必ず加速が伴い、加速する電荷はラーモアの公式に従って電磁波を放射します。これが放射の最も根本的な原因です。
ステップ2: 電磁ポテンシャルの生成
加速する電荷（時間変化する電流）は、遅延ポテンシャルを通じて空間にベクトルポテンシャル $\bm{A}$ とスカラーポテンシャル $\phi$ を生成します。遅延時間 $R/c$ は、情報が光速で伝わることの帰結です。
ステップ3: 電磁場の形成
ポテンシャルから $\bm{B} = \nabla \times \bm{A}$ および $\bm{E} = -\nabla\phi – \partial\bm{A}/\partial t$ により電磁場が決まります。遠方では $1/r$ に比例する放射界が支配的になり、これが電磁波として伝播します。
ステップ4: 放射パターンの決定
放射界の角度依存性は、アンテナ上の電流分布によって決まります。微小ダイポールでは $\sin\theta$ の依存性が現れ、ドーナツ型の放射パターンとなります。実際のアンテナは微小ダイポールの重ね合わせで表現でき、電流分布を適切に設計することで所望の放射パターンを実現できます。
なぜ導線だけでは電磁波を放射しないのか
ここで1つの疑問が浮かぶかもしれません。「送電線にも交流電流が流れているのに、なぜ電磁波を大量に放射しないのか？」
答えは、2線式送電線では2本の導線を逆方向に流れる電流が作る電磁場が互いにほぼ打ち消し合うからです。2本の線が近接している（線間距離が波長に比べて十分小さい）限り、遠方では放射がほとんどキャンセルされます。
アンテナはこの打ち消しを意図的に崩した構造です。ダイポールアンテナは2線式伝送路の端を開いて広げた構造であり、電流の方向が揃った部分（放射に寄与する部分）を作り出しています。
この「打ち消し」と「意図的な不均衡」の理解は、アンテナ設計の根底にあるアイデアです。アレイアンテナでは、逆に複数のアンテナ素子からの放射を特定の方向で「強め合い」、他の方向で「打ち消し合い」させることで、高い指向性を実現します。
まとめ
本記事では、アンテナの放射メカニズムについて、加速する電荷という根本原理から出発し、マクスウェル方程式、遅延ポテンシャル、微小ダイポールの放射界の導出まで一貫した流れで解説しました。
- 放射の本質: 加速する電荷が電磁波を放射する。アンテナに交流電流を流すと電子が振動（加速）し、これが放射の源となる
- 遅延ポテンシャル: 電流が作るポテンシャルは光速で伝播するため、遅延時間 $R/c$ が生じる。この遅延効果が波動としての性質を生む
- 近傍界と遠方界: $1/r^3$、$1/r^2$、$1/r$ の各項が支配する領域があり、遠方界（$r \gg \lambda$）で放射界が支配的になる
- 微小ダイポールの放射パターン: $\sin^2\theta$ のドーナツ型パターン。指向性 $D = 1.5$、放射抵抗 $R_r = 80\pi^2(dl/\lambda)^2$
- 放射の制御: 電流分布を設計することで放射パターンを制御できる。これがアンテナ設計の本質
この概念を基盤として、次はアンテナの実用的な評価指標や設計手法を学んでいくことになります。
次のステップとして、以下の記事も参考にしてください。