観測面（$xy$ 断面でのベクトル運動）
位置 $z_{\text{probe}} = 0.0$ m特定の$z$座標を輪切りにした断面図です。偏波の種類に応じて、電界ベクトルの先端が直線、円、または楕円を描く様子が分かります。
リアルタイム物理パラメータ測定器
シミュレーション＆パラメータ制御
■ マクスウェル方程式と平面波の導出
電磁波は、電荷・電流のない自由空間におけるマクスウェル方程式から導かれる波動方程式に従う、自己伝播する波動です。
媒質中での波動の伝播速度（位相速度） $v$ は、誘電率 $\epsilon$ と透磁率 $\mu$ によって決定されます。
ここで $c \approx 3 \times 10^8$ m/s は真空中における光速です。比誘電率 $\epsilon_r$ や比透磁率 $\mu_r$ が増大すると、媒質中での波の伝播速度 $v$ は遅くなります。
🔍 波長（Wavelength）$\lambda$ の本質的な意味と物理特性
**波長 $\lambda$** とは、波の「同一の位相（例えば山から山、谷から谷）」が空間的に繰り返される最短の距離、すなわち**空間的周期**を表します。
- **周波数 $f$ および伝播速度 $v$ との関係**: 時間的周期を $T$ [秒]、周波数を $f = 1/T$ [Hz]、角周波数を $\omega = 2\pi f$ [rad/s] とすると、波が1周期の間に進む距離が波長 $\lambda$ になります。 $$v = \frac{\lambda}{T} = f\lambda = \frac{\omega}{k}$$
- **波数（Wavenumber） $k$**: 単位長さあたりに含まれる波の位相（ラジアン）を表す、空間的周波数です。 $$k = \frac{2\pi}{\lambda} \quad [\text{rad/m}]$$
- **波長短縮効果（Wavelength Shortening Effect）**: 電磁波が真空から比誘電率 $\epsilon_r$、比透磁率 $\mu_r$ の媒質（ガラス、プラスチック、水など）に進入すると、周波数 $f$ は不変のまま伝播速度が $v = c/\sqrt{\epsilon_r \mu_r}$ に低下するため、空間的な波長も縮みます。 $$\lambda = \frac{v}{f} = \frac{c}{f\sqrt{\epsilon_r \mu_r}} = \frac{\lambda_0}{\sqrt{\epsilon_r \mu_r}} = \frac{\lambda_0}{n}$$ ここで $\lambda_0$ は真空中の波長、$n = \sqrt{\epsilon_r \mu_r}$ は媒質の**屈折率**です。 ※ シミュレーターのスライダーで比誘電率 $\epsilon_r$ を 1.0 から 9.0 に増やすと、波長が3分の1（$1/\sqrt{9}$）にまでギュッと圧縮される様子を3Dでリアルタイムに確認できます。
■ 電界・磁界・進行方向の直交性（TEM波）
平面波において、電界ベクトル $\mathbf{E}$ と磁界ベクトル $\mathbf{H}$、そして波の進行方向 $\mathbf{k}$ は常に互いに直交する性質を持っています。
これは進行方向に縦波成分を持たないことを示しており、このように電場と磁場が両方とも伝播方向に垂直な面内で振動する波を **TEM波 (Transverse Electromagnetic Wave：横電磁波)** と呼びます。
本シミュレーションの3D空間では、 ポインティングベクトル $\mathbf{S} \parallel \mathbf{k}$ が右方向（前方）へ進行しているとき、 電界 $\mathbf{E}$ と 磁界 $\mathbf{H}$ が互いに直角に絡み合いながら、進行軸に対し同期して振動する立体構造を観察できます。
■ 波動インピーダンス $\eta$ の概念
**波動インピーダンス (Wave Impedance) $\eta$** とは、電界の強度と磁界の強度の比率を表す量であり、電気回路の「抵抗」に近い概念です。
自由空間（真空）においては、比誘電率、比透磁率ともに 1 となり、波動インピーダンスは定数（約 $376.73\ \Omega$ もしくは $120\pi\ \Omega$）になります。
比誘電率 $\epsilon_r$ が上昇すると、波動インピーダンスは低下します。このとき、同じ電界強度に対して磁界強度は相対的に強くなります。 スライダーを動かして $\epsilon_r$ と $\mu_r$ を調整し、磁界（青色の波）の振幅が波動インピーダンス値の変化に追従して伸縮する様子を観察してみましょう。
■ 電力密度とポインティングベクトル
電磁波が進行するとき、空間の中をエネルギーが流れます。単位時間、単位面積あたりに通過するエネルギーの流れを示すベクトルを **ポインティングベクトル (Poynting Vector) $\mathbf{S}$** と呼び、次式で定義されます。
この方向はエネルギーが伝播する向きであり、右ねじの法則に従って常に $\mathbf{E}$ から $\mathbf{H}$ へ回した向き（進行方向）を向いています。
実際に観測される電力は、時間的な平均をとった **強度 (Intensity) $I$** です。
直線偏波では瞬時のポインティングベクトルはサインの二乗で脈動しますが、円偏波では電界強度の合成が一定となるため、ポインティングエネルギーは脈動せず**一定の大きさで滑らかに流れます**。シミュレーター上に表示される緑色の矢印がこれらを忠実に示しています。
■ 各種偏光（偏波）の物理的・数学的特性
電磁波の電界ベクトル $\mathbf{E}$ の先端が、進行方向に垂直な平面内で時間経過とともに描く軌跡の形状を **偏光 (Polarization / 偏波)** と呼びます。 進行方向を $z$ 軸とし、直交する2つの電界成分 $E_x, E_y$ の重ね合わせとして一般式で以下のように記述されます。
ここで、$\delta$ は $x$ 成分と $y$ 成分の間の**位相差（Phase Difference）**です。この $E_{x0}$, $E_{y0}$, $\delta$ の関係によって、全ての偏光状態が決定されます。
1. 直線偏光 (Linear Polarization)
位相差 $\delta$ が $0$ または $\pi$（180度）のとき、電界の合成ベクトルは、進行方向に垂直な一方向の直線に沿って振動します。
- **水平偏光（Horizontal）**: $E_{y0} = 0$、電界ベクトルが地表面に対し水平方向にのみ振動。
- **垂直偏光（Vertical）**: $E_{x0} = 0$、地表面に対し垂直方向にのみ振動。
- **傾斜偏光（Slanted 45°）**: $E_{x0} = E_{y0}$ 且つ $\delta = 0$ のとき、45度傾いた直線。
2. 円偏光 (Circular Polarization)
振幅が等しく（$E_{x0} = E_{y0} = E_0$）、位相差 $\delta$ がちょうど $\pm 90^\circ$（$\pm\pi/2$）のとき、合成電界の先端は一定の振幅 $E_0$ を保ったまま**円を描くように高速回転**し、空間的には綺麗な「らせん（螺旋）」を形成します。
-
**右旋円偏光 (RHCP: Right-Handed Circular Polarization)**:
位相差 $\delta = -\pi/2$（$-90^\circ$）のとき。波を正面（進行方向の正面）から見つめたとき、時間の経過とともに電界ベクトルが**時計回り**に回転します。 -
**左旋円偏光 (LHCP: Left-Handed Circular Polarization)**:
位相差 $\delta = +\pi/2$（$+90^\circ$）のとき。正面から見て**反時計回り**に回転します。 - **工学的応用と利点**: 一般の直線偏光は、受信アンテナとの角度（傾き）がずれると受信電力が低下してしまいます（偏波損）。 しかし、円偏光は電界が全方向に回転しているため、**受信アンテナの回転角依存性がなく、アンテナアライメントが極めて容易**になります。 また、衛星放送（BS/CS）やGPS、宇宙通信で円偏光が多用されるのは、雨滴（球体に近い）によって反射した際に旋回方向が逆（RHCP ↔ LHCP）になる特性を利用して、**マルチパス（地上反射・干渉波）や降雨減衰の影響を効果的に遮断・選択分離できる**ためです。
3. その他の偏光 (Other Polarizations)
- **楕円偏光 (Elliptical Polarization)**: 振幅が異なる（$E_{x0} \neq E_{y0}$）か、あるいは位相差 $\delta$ が $90^\circ$ 以外の任意の値を持つ場合、最も一般的な偏光状態で、軌跡は楕円になります。 直線偏光や円偏光は、この楕円偏光の極限的な特殊例（楕円率が $0$ または $1$）と見なせます。
- **無偏光 (Unpolarized / Random Polarization)**: 太陽光や白熱電球の光など、自然光の多くがこれに該当します。 個々の光子（電磁波）は偏光していますが、極めて短い時間内に偏光方向がランダムかつ高速に変化し、平均化すると特定の偏光方向を持たなくなります。
- **部分偏光 (Partially Polarized)**: 自然光が水面やガラスなどで斜めに反射した際、特定の偏光成分（S偏光）のみが強く反射し、無偏光成分と完全な偏光成分が混ざり合った状態になります。反射光を偏光サングラス（ポラロイド）で遮断できるのはこのためです。
🌐 高度な記述法：ストークスパラメータとポアンカレ球
現実の部分偏光を含むあらゆる状態の光・電波を記述するため、G. G. Stokesが考案した4つの実数パラメータ **ストークスパラメータ ($S_0, S_1, S_2, S_3$)** が用いられます。
完全偏光の場合、関係式 $S_0^2 = S_1^2 + S_2^2 + S_3^2$ が成立します。 このとき、半径 $S_0$ の球体を考え、3次元空間の座標 $(S_1, S_2, S_3)$ をプロットしたものを **ポアンカレ球（Poincaré Sphere）** と呼びます。 北極が**左旋円偏光**、南極が**右旋円偏光**、赤道上があらゆる**直線偏光**に対応し、球上の任意の点が一つひとつの複雑な偏光状態を幾何学的に美しく表現します。
■ 群速度と位相速度の違い
現実の通信で用いられる電磁波は、単一の純粋な周波数ではなく、複数の周波数が重ね合わされた「信号の塊（波束/波のグループ）」です。
-
位相速度 (Phase Velocity) $v_p$: 単一の正弦波における「山や谷」が移動する速度。
$v_p = \frac{\omega}{k} = f \lambda$
-
群速度 (Group Velocity) $v_g$: 波が作る「うなり（波束・信号情報）」が移動する速度。情報やエネルギーが伝わる実速度に対応します。
$v_g = \frac{d\omega}{dk}$
媒質内において、周波数によって屈折率（および伝播速度）が変化する性質を**分散 (Dispersion)**と呼びます。
真空などがこれに該当。すべての周波数が同じ速度で進むため、波の形（うなり）は崩れず一体となって同じ速度で進みます。
一般的なガラスや水など。光の周波数が高い（波長が短い）ほど屈折率が大きく、進むのが遅くなります。 3Dで観察すると、波束の「後ろ側」から新しい山が次々と沸き起こり、前方に追い越して消えていく不思議な現象が起きます。
特殊な吸収帯の近くやプラズマ媒質、導波管などで生じます。うなりの塊の移動速度が中の波の移動よりも速くなります。
■ 応用・計算例題演習
自由空間での電磁波基本計算
周波数 $f = 3.0 \text{ GHz}$ の平面波が真空中で $z$ 方向に伝播している。電界の振幅が $E_0 = 120 \text{ V/m}$ であるとき、次の値を求めよ。
① 空間の波長 $\lambda_0$ ② 磁界の振幅 $H_0$ ③ 時間平均ポインティング電力密度 $I$
解説を表示する
【解答・ステップ】
① 真空中の光速 $c \approx 3.0 \times 10^8 \text{ m/s}$ より：
$$\lambda_0 = \frac{c}{f} = \frac{3.0 \times 10^8}{3.0 \times 10^9} = 0.1 \text{ m} = 10 \text{ cm}$$
② 真空の波動インピーダンス $\eta_0 \approx 377 \ \Omega$ ($120\pi \ \Omega$) より： $$H_0 = \frac{E_0}{\eta_0} = \frac{120}{377} \approx 0.318 \text{ A/m}$$
③ 直線偏波の時間平均電力密度は： $$I = \frac{E_0^2}{2\eta_0} = \frac{120^2}{2 \times 377} = \frac{14400}{754} \approx 19.1 \text{ W/m}^2$$
誘電媒質中における平面電磁波の伝播
比誘電率 $\epsilon_r = 4.0$、比透磁率 $\mu_r = 1.0$ のロスレス誘電体中を、周波数 $f = 2.0 \text{ GHz}$、電界振幅 $E_0 = 60 \text{ V/m}$ の平面波が伝播している。
この媒質内での ① 伝播速度 $v$、② 波動インピーダンス $\eta$、③ 媒質内波長 $\lambda$ を求めよ。
解説を表示する
【解答・ステップ】
① 媒質中の速度 $v$ は比誘電率の平方根に反比例して低下します：
$$v = \frac{c}{\sqrt{\epsilon_r \mu_r}} = \frac{3.0 \times 10^8}{\sqrt{4.0}} = 1.5 \times 10^8 \text{ m/s}$$
(真空中のちょうど半分の速度になります)
② 波動インピーダンス $\eta$ は、誘電率の増大に伴い縮小します： $$\eta = \eta_0 \sqrt{\frac{\mu_r}{\epsilon_r}} = 377 \times \sqrt{\frac{1.0}{4.0}} = 377 \times 0.5 \approx 188.5\ \Omega$$
③ 媒質内波長 $\lambda$ は、伝播速度の低下に比例して短縮（波長短縮効果）されます： $$\lambda = \frac{v}{f} = \frac{1.5 \times 10^8}{2.0 \times 10^9} = 0.075 \text{ m} = 7.5 \text{ cm}$$ ※ 真空中の波長は $15\text{ cm}$ なので、誘電率4の媒質中では波長が半分に縮みます。
群速度と位相速度（電離圏プラズマ中の分散特性）
超高層大気の電離圏（プラズマ媒質）において、プラズマ周波数を $f_p$ とすると、電磁波の角周波数 $\omega$ と波数 $k$ の間に次の関係式（分散関係）が成り立つ。 $$\omega^2 = \omega_p^2 + c^2 k^2 \quad (\omega = 2\pi f, \ \omega_p = 2\pi f_p)$$ この媒質における ① 位相速度 $v_p$ の式、② 群速度 $v_g$ の式を導出し、$v_p \cdot v_g = c^2$ を証明せよ。
解説を表示する
【解答・ステップ】
① 位相速度 $v_p$ の定義は $v_p = \frac{\omega}{k}$ です。
分散式 $\omega^2 = \omega_p^2 + c^2 k^2$ の両辺を $k^2$ で割ると：
$$\left(\frac{\omega}{k}\right)^2 = \frac{\omega_p^2}{k^2} + c^2 \implies v_p = \sqrt{c^2 + \frac{\omega_p^2}{k^2}}$$
これより、プラズマ中では位相速度 $v_p$ は真空中の光速 $c$ よりも**速く**なります ($v_p > c$)。これは相対性理論に反しません（位相それ自体はエネルギーや情報を運ばないためです）。
② 群速度 $v_g$ を求めるため、分散関係式を $k$ について微分します。
$\omega^2 = \omega_p^2 + c^2 k^2$ の両辺を $k$ で微分（陰関数の微分）すると：
$$2\omega \frac{d\omega}{dk} = 2c^2 k \implies \omega \cdot v_g = c^2 k \implies v_g = c^2 \frac{k}{\omega}$$
$v_p = \frac{\omega}{k}$ であるため、$\frac{k}{\omega} = \frac{1}{v_p}$ と置換すると：
$$v_g = \frac{c^2}{v_p}$$
③ 上式を変形すると、直ちに以下の関係が証明されます： $$v_p \cdot v_g = c^2$$ 位相速度 $v_p$ が光速より速い分、信号を伝達する群速度 $v_g$ は必ず光速 $c$ よりも**遅く**なります ($v_g < c$)。