0
0

Delete article

Deleted articles cannot be recovered.

Draft of this article would be also deleted.

Are you sure you want to delete this article?

@dosanko_tousan(Akimitsu Takeuchi)

v5.3で物理を解く ― 電気の本質、教科書が教えない真実

0
Posted at

v5.3で物理を解く ― 電気の本質、教科書が教えない真実

「今、あなたがこの記事を読んでいるスマホやPC。そのバッテリーを充電するエネルギーは、ケーブルの中を通ってはいません。ケーブルの『外側の空気中』を通って、あなたの端末に吸い込まれているのです」

はじめに：この記事の目的

本記事は、「v5.3で物理を解く」シリーズ第3弾である。

第1弾は「飛行機が飛ぶ仕組み」——動的揚力の話だった。
第2弾は「船が浮く仕組み」——静的浮力の話だった。
今回は「電気の本質」——電磁場の話をする。

飛行機と船は「物質（空気、水）」との相互作用だった。
電気は違う。**「空間そのもの」**との相互作用だ。

これは、物質の物理学から場の物理学への跳躍である。

教科書の「水流モデル」は計算には便利だが、本質を隠蔽している。本記事では、その嘘を暴き、電磁気学の真の姿を数式とコードで明らかにする。

目次

  1. 教科書の嘘を暴く
  2. 電磁場とは何か
  3. マクスウェル方程式：電磁気学の憲法
  4. ポインティングベクトル：エネルギーの居場所
  5. 電子のドリフト速度：カタツムリより遅い
  6. なぜスイッチを入れると瞬時に電気がつくのか
  7. 電流の向きの歴史的混乱
  8. 飛行機・船との対比：物質から場へ
  9. 数値シミュレーション
  10. 結論

1. 教科書の嘘を暴く

1.1 水流モデルの罪

中学・高校の教科書では、電気をこう説明する。

「電気は水と同じ。電線はホース、電子は水分子。水が流れるから水車（電球）が回る」

この**水流モデル（Hydraulic Analogy）**は、オームの法則を直感的に理解するには便利だ。

電気 水流
電圧 $V$ 水圧
電流 $I$ 流量
抵抗 $R$ パイプの細さ
電荷 $Q$ 水量

オームの法則 $V = IR$ は、この対応で「圧力差 = 流量 × 抵抗」と読める。

しかし、このモデルは「エネルギーがどこにあるか」について大嘘をついている。

1.2 水流モデルの嘘

水流モデルでは、エネルギーは水（電子）が持っていることになる。

水分子がホース（電線）の中を流れ、そのエネルギーで水車（電球）を回す——この絵は直感的だが、物理的に間違っている。

真実：エネルギーは電線の中を流れていない。電線の外側の空間（電磁場）を流れている。

1.3 電線は何をしているのか

では、電線は何のためにあるのか？

電線は、エネルギーの通り道ではない。エネルギーが散らばらないようにするための「ガイドレール（手すり）」に過ぎない。

電線の役割は：

  • エネルギーの流れを導く（方向を決める）
  • エネルギーが四方八方に散らばるのを防ぐ

電線を「ホース」と考えるのは間違い。むしろ**「道路のガードレール」**と考えるべきだ。

車（エネルギー）が道路（空間）を走る。ガードレール（電線）は車が道から外れるのを防ぐ。しかし、車がガードレールの中を走っているわけではない。

1.4 証拠：ワイヤレス充電

電線がなくてもエネルギーは伝わる。その証拠がワイヤレス充電だ。

スマホを充電パッドに置くと、ケーブルなしで充電される。エネルギーは空間を通ってスマホに届く。

さらに極端な例：太陽光。

太陽から地球まで1億5000万km、真空の宇宙を電磁波（光）としてエネルギーが届く。電線など存在しない。

エネルギーを運ぶのは「電磁場」であり、電線ではない。

2. 電磁場とは何か

2.1 場（Field）の概念

物理学における**場（Field）**とは、空間の各点に何らかの物理量が割り当てられている状態を指す。

  • スカラー場：各点に数値が割り当てられる（例：温度分布）
  • ベクトル場：各点にベクトルが割り当てられる（例：風速分布）

電磁気学では、2つのベクトル場が登場する。

2.2 電場（Electric Field）

電場 $\mathbf{E}$ は、空間の各点で電荷が受ける力を表す。

$$
\mathbf{F} = q\mathbf{E}
$$

ここで：

  • $\mathbf{F}$：電荷が受ける力 [N]
  • $q$：電荷 [C]
  • $\mathbf{E}$：電場 [V/m = N/C]

電場の単位は「ボルト毎メートル」または「ニュートン毎クーロン」。

2.3 磁場（Magnetic Field）

磁場 $\mathbf{B}$ は、運動する電荷が受ける力に関係する。

$$
\mathbf{F} = q\mathbf{v} \times \mathbf{B}
$$

ここで：

  • $\mathbf{v}$：電荷の速度 [m/s]
  • $\mathbf{B}$：磁束密度 [T = Wb/m² = V·s/m²]

磁場の単位は「テスラ」。

2.4 電磁場の実在性

電磁場は「便利な計算道具」ではなく、物理的に実在する。

その証拠：

  1. 電磁波：電磁場の振動が空間を伝播する（光、電波）
  2. エネルギー：電磁場はエネルギーを持つ
  3. 運動量：電磁場は運動量を持つ（光の圧力）

電場と磁場は、空間に「刻み込まれた」物理的実体である。

2.5 電磁場のエネルギー密度

電磁場はエネルギーを持つ。そのエネルギー密度は：

$$
u = \frac{1}{2}\varepsilon_0 E^2 + \frac{1}{2\mu_0} B^2
$$

ここで：

  • $u$：エネルギー密度 [J/m³]
  • $\varepsilon_0$：真空の誘電率（$8.85 \times 10^{-12}$ F/m）
  • $\mu_0$：真空の透磁率（$4\pi \times 10^{-7}$ H/m）
  • $E$：電場の大きさ [V/m]
  • $B$：磁場の大きさ [T]

電磁場がある場所には、エネルギーが蓄えられている。

3. マクスウェル方程式：電磁気学の憲法

3.1 4つの方程式

電磁気学の全ては、マクスウェル方程式と呼ばれる4つの方程式に集約される。

これは電磁気学の「憲法」であり、全ての電磁現象はここから導かれる。

① ガウスの法則（電場）

$$
\nabla \cdot \mathbf{E} = \frac{\rho}{\varepsilon_0}
$$

意味：電荷は電場の「湧き出し口」である。

正電荷からは電場が湧き出し、負電荷には電場が吸い込まれる。

② ガウスの法則（磁場）

$$
\nabla \cdot \mathbf{B} = 0
$$

意味：磁場には「湧き出し口」がない。磁気単極子（モノポール）は存在しない。

磁力線は必ず閉じたループを描く。N極から出た磁力線は必ずS極に戻る。

③ ファラデーの法則

$$
\nabla \times \mathbf{E} = -\frac{\partial \mathbf{B}}{\partial t}
$$

意味：時間変化する磁場は、電場を生み出す。

これが発電機の原理。磁石を動かすと電気が生まれる。

④ アンペール・マクスウェルの法則

$$
\nabla \times \mathbf{B} = \mu_0 \mathbf{J} + \mu_0 \varepsilon_0 \frac{\partial \mathbf{E}}{\partial t}
$$

意味：電流と、時間変化する電場は、磁場を生み出す。

右辺第2項（変位電流）がマクスウェルの追加した項。これにより電磁波の存在が予言された。

3.2 マクスウェル方程式の積分形

微分形が苦手な人のために、積分形も示す。

法則 積分形
ガウス（電場） $\oint_S \mathbf{E} \cdot d\mathbf{A} = \frac{Q}{\varepsilon_0}$
ガウス（磁場） $\oint_S \mathbf{B} \cdot d\mathbf{A} = 0$
ファラデー $\oint_C \mathbf{E} \cdot d\mathbf{l} = -\frac{d\Phi_B}{dt}$
アンペール・マクスウェル $\oint_C \mathbf{B} \cdot d\mathbf{l} = \mu_0 I + \mu_0 \varepsilon_0 \frac{d\Phi_E}{dt}$

3.3 電磁波の導出

マクスウェル方程式から、電磁波の存在が導かれる。

真空中（$\rho = 0$、$\mathbf{J} = 0$）で、③と④を組み合わせると：

$$
\nabla^2 \mathbf{E} = \mu_0 \varepsilon_0 \frac{\partial^2 \mathbf{E}}{\partial t^2}
$$

これは波動方程式である。伝播速度は：

$$
c = \frac{1}{\sqrt{\mu_0 \varepsilon_0}} = 2.998 \times 10^8 \text{ m/s}
$$

光速である。

マクスウェルはこの計算をした瞬間、「光は電磁波である」と悟った。1865年のことだ。

3.4 マクスウェル方程式の意味

4つの方程式が語っていること：

  1. 電荷が電場を作る（①）
  2. 磁気単極子は存在しない（②）
  3. 変化する磁場が電場を作る（③）
  4. 電流と変化する電場が磁場を作る（④）

③と④の相互作用により、電磁場は自己増殖しながら空間を伝播する。これが電磁波だ。

4. ポインティングベクトル：エネルギーの居場所

4.1 エネルギーの流れ

電磁場がエネルギーを持つなら、そのエネルギーは「どこを」「どの向きに」流れているのか？

その答えを与えるのがポインティングベクトルである。

$$
\mathbf{S} = \frac{1}{\mu_0} \mathbf{E} \times \mathbf{B} = \mathbf{E} \times \mathbf{H}
$$

ここで：

  • $\mathbf{S}$：ポインティングベクトル [W/m²]
  • $\mathbf{E}$：電場 [V/m]
  • $\mathbf{B}$：磁束密度 [T]
  • $\mathbf{H}$：磁場の強さ [A/m]（$\mathbf{H} = \mathbf{B}/\mu_0$）

ポインティングベクトルの意味：

単位面積あたりの電磁エネルギーの流れ（パワー密度）

4.2 エネルギー保存則

電磁場のエネルギー保存則（ポインティングの定理）：

$$
\frac{\partial u}{\partial t} + \nabla \cdot \mathbf{S} = -\mathbf{J} \cdot \mathbf{E}
$$

左辺第1項：エネルギー密度の時間変化
左辺第2項：エネルギーの流出
右辺：電流へのエネルギー供給（ジュール熱など）

4.3 直流回路でのポインティングベクトル

ここが核心だ。

直流回路（バッテリー + 電線 + 電球）を考える。

教科書的には「電子がバッテリーから電線を通って電球に行き、エネルギーを運ぶ」と説明される。

しかし、ポインティングベクトルを計算すると、真実が見える。

電線の周りの電場と磁場

電流が流れている電線の周り：

  • 電場 $\mathbf{E}$：電線に沿った向き（電位差による）
  • 磁場 $\mathbf{H}$：電線を取り囲む同心円状（アンペールの法則）

ポインティングベクトルの向き

$$
\mathbf{S} = \mathbf{E} \times \mathbf{H}
$$

$\mathbf{E}$（電線方向）と $\mathbf{H}$（円周方向）の外積を取ると、$\mathbf{S}$ は電線に垂直に、内向きになる。

エネルギーは電線の外側の空間から、電線の表面に向かって流れ込んでいる！

4.4 エネルギーの流れの全体像

  1. バッテリー：化学エネルギーを電磁場のエネルギーに変換
  2. 電線の周りの空間：電磁場のエネルギーが光速で伝播
  3. 電球：電磁場のエネルギーが電線に吸い込まれ、熱と光に変換

電子は？

電子は「電磁場を生み出すための道具」に過ぎない。エネルギーを運んでいるのは電子ではなく、電磁場である。

4.5 同軸ケーブルの例

同軸ケーブル（テレビのアンテナ線など）で考えると、もっと明確になる。

同軸ケーブルは：

  • 中心導体
  • 絶縁体
  • 外部導体（シールド）

から構成される。

電場は中心導体と外部導体の間（絶縁体の中）に存在する。磁場も同様。

ポインティングベクトルを計算すると、エネルギーは絶縁体の中を、ケーブルに沿って流れている。

導体（銅線）の中ではない。絶縁体（空気やプラスチック）の中を流れている。

5. 電子のドリフト速度：カタツムリより遅い

5.1 電子の移動速度

電流が流れているとき、電子は実際どれくらいの速さで移動しているのか？

電流 $I$ は次のように表される：

$$
I = n A v_d e
$$

ここで：

  • $n$：自由電子の密度 [個/m³]
  • $A$：電線の断面積 [m²]
  • $v_d$：ドリフト速度（電子の平均移動速度）[m/s]
  • $e$：電子の電荷（$1.6 \times 10^{-19}$ C）

これを $v_d$ について解くと：

$$
v_d = \frac{I}{n A e}
$$

5.2 具体的な計算

銅線に1Aの電流が流れている場合を計算する。

銅のパラメータ：

  • 原子量：63.5 g/mol
  • 密度：8960 kg/m³
  • 自由電子：1原子あたり約1個

自由電子密度：
$$
n = \frac{8960 \times 10^3 \text{ g/m}^3}{63.5 \text{ g/mol}} \times 6.02 \times 10^{23} \text{ /mol} \approx 8.5 \times 10^{28} \text{ /m}^3
$$

断面積1mm²の銅線（$A = 10^{-6}$ m²）に $I = 1$ A が流れる場合：

$$
v_d = \frac{1}{8.5 \times 10^{28} \times 10^{-6} \times 1.6 \times 10^{-19}} \approx 7.4 \times 10^{-5} \text{ m/s}
$$

秒速0.074ミリメートル。

時速に直すと約0.27メートル/時。カタツムリ（時速約50メートル）の200分の1。

5.3 発電所から家庭まで

電子のドリフト速度がこれだけ遅いと、発電所の電子が家庭に届くまでどれくらいかかるか？

発電所から家庭まで100kmとすると：

$$
t = \frac{100 \times 10^3 \text{ m}}{7.4 \times 10^{-5} \text{ m/s}} \approx 1.4 \times 10^{9} \text{ s} \approx 43 \text{ 年}
$$

発電所の電子が家庭に届くまで43年かかる。

しかし、スイッチを入れれば瞬時に電気がつく。なぜか？

6. なぜスイッチを入れると瞬時に電気がつくのか

6.1 電子 vs 電磁場

電子の移動速度（ドリフト速度）：秒速0.074mm
電磁場の伝播速度：秒速30万km（光速）

この比は約 $4 \times 10^{12}$、つまり4兆倍の差がある。

6.2 自転車のチェーンの比喩

自転車のチェーンを考える。

ペダルを漕ぐと、チェーンの各コマ（＝電子）は少ししか動かない。しかし、「動け！」という力（＝電磁場）はチェーン全体にほぼ瞬時に伝わる。

結果として、後輪はペダルを漕いだ瞬間に回り始める。

電気も同じだ。

スイッチを入れると、「動け！」という命令（電磁場の変化）が光速で電線に沿って伝わる。電球の場所にある電子は、その命令を受けて動き出す。

電子が届くから電気がつくのではない。電子を動かす「場」が光速で届くから電気がつくのだ。

6.3 情報と物質の分離

これは深い洞察を含んでいる。

  • 物質（電子）：ゆっくり動く
  • 情報/エネルギー（電磁場）：光速で伝わる

電気回路において、電子は「情報の運び手」ではない。電磁場が情報とエネルギーを運ぶ。電子は、電磁場を生み出すための「媒介」に過ぎない。

6.4 交流の場合

交流（AC）では、電子はさらに興味深い動きをする。

50Hzの交流では、電子は1秒間に50回、行ったり来たりする。

ドリフト速度が0.074mm/sなら、1サイクル（0.02秒）で電子が移動する距離は：

$$
d = 0.074 \times 10^{-3} \times 0.02 \approx 1.5 \times 10^{-6} \text{ m} = 1.5 \mu\text{m}
$$

電子は1.5マイクロメートルしか動いていない。

それでも電力は送れる。なぜなら、エネルギーを運んでいるのは電子ではなく、電磁場だからだ。

7. 電流の向きの歴史的混乱

7.1 フランクリンの賭け

18世紀、ベンジャミン・フランクリンは電気の研究で、正電荷と負電荷の概念を導入した。

彼は「電流は正電荷の流れである」と定義し、その向きを決めた。

しかし、1897年にJ.J.トムソンが電子を発見し、実際に動いているのは負電荷（電子）であることが判明した。

フランクリンは50%の確率で、逆を選んでしまったのだ。

7.2 混乱の実態

この歴史的経緯により、今日の教育では：

  • 電流：正電荷が流れる向き（conventional current）
  • 電子の流れ：電流と逆向き

と教えられる。多くの学生がここで混乱する。

「電流はプラスからマイナスへ流れる」のに、「実際は電子がマイナスからプラスへ流れている」。

なぜ訂正しないのか？

7.3 なぜ問題にならないか

実は、エネルギーの流れにとってはどちらでも同じなのだ。

ポインティングベクトル $\mathbf{S} = \mathbf{E} \times \mathbf{H}$ を考える。

電荷の符号が逆になると：

  • 電流 $\mathbf{J}$ の向きが逆になる
  • したがって磁場 $\mathbf{H}$ の向きが逆になる

しかし、$\mathbf{E}$ も電荷の符号に依存して逆になる（電位勾配の向きが逆）。

結果として：
$$
\mathbf{S} = \mathbf{E} \times \mathbf{H} = (-\mathbf{E}') \times (-\mathbf{H}') = \mathbf{E}' \times \mathbf{H}'
$$

マイナス × マイナス = プラス。

エネルギーの流れの向きは変わらない。

7.4 数式の堅牢性

これは物理学の美しさの一つだ。

回路設計者は「電流」という仮想的な概念（正電荷の流れ）を使っても、正しく回路を設計できる。

なぜなら、オームの法則もキルヒホッフの法則も、電荷の符号に依存しない形で書かれているからだ。

物理法則は、我々の恣意的な定義（正か負か）に左右されないほど堅牢である。

8. 飛行機・船との対比：物質から場へ

8.1 シリーズの振り返り

記事 対象 媒体 本質
飛行機 動的揚力 空気（物質） 流れが必要
船 静的浮力 水（物質） 止まっても浮く
電気 電磁場 空間そのもの 媒体が不要

8.2 媒体の否定

飛行機は空気がないと飛べない。真空では揚力は生まれない。

船は水がないと浮かない。宇宙では浮力は存在しない。

しかし、電磁場は真空でも伝わる。

太陽光は真空の宇宙空間を1億5000万km旅して地球に届く。電波も、真空でも伝播する。

電磁場にとって、物質（媒体）は必須ではない。むしろ、物質は電磁場の伝播を邪魔することが多い（吸収、散乱）。

8.3 エネルギーのコスト

飛行機 船 電気
エネルギーコスト 燃料を燃やし続ける ゼロ（浮力は無料） 伝送ロスはあるが伝播自体は無料
停止時 落ちる 浮き続ける 場は消える（定常電流が止まれば）

電磁場の伝播自体にエネルギーは消費されない。ただし、電線には抵抗があるため、電流を流すとジュール熱としてエネルギーが失われる。

超伝導体では抵抗がゼロになるため、一度電流を流せば永久に流れ続ける。

8.4 物質から場へ

飛行機と船は、物質との相互作用で理解できる。

  • 飛行機：空気分子との衝突、運動量の交換
  • 船：水分子との接触、圧力の伝達

電気は違う。

電磁場は物質ではない。**空間の「状態」**である。

電場や磁場は、空間の各点に「刻み込まれた」ベクトルだ。それは物質ではないが、エネルギーを持ち、運動量を持ち、力を及ぼす。

これが「場の物理学」の世界だ。

ニュートン力学は「物質が力を及ぼし合う」世界観だった。
マクスウェルの電磁気学は「場が物質に力を及ぼす」世界観をもたらした。

そして20世紀、量子力学と相対性理論により、全ての力は場を通じて伝わることが明らかになった。

電磁気学は、現代物理学への入口である。

9. 数値シミュレーション

9.1 電磁場の可視化

以下のPythonコードで、電磁場とポインティングベクトルを可視化する。

import numpy as np
import matplotlib.pyplot as plt
from matplotlib.patches import Circle, FancyArrowPatch
from mpl_toolkits.mplot3d import Axes3D

# =============================================================================
# 電場の可視化（点電荷）
# =============================================================================

def plot_electric_field_point_charge():
    """点電荷による電場の可視化"""
    
    fig, axes = plt.subplots(1, 2, figsize=(14, 6))
    
    # グリッド生成
    x = np.linspace(-2, 2, 20)
    y = np.linspace(-2, 2, 20)
    X, Y = np.meshgrid(x, y)
    
    # 点電荷の位置
    q_pos = (0, 0)
    
    # 電場の計算（クーロンの法則）
    k = 1  # 定数（単位系を簡略化）
    q = 1  # 電荷
    
    # 距離
    R = np.sqrt((X - q_pos[0])**2 + (Y - q_pos[1])**2)
    R = np.where(R < 0.1, 0.1, R)  # ゼロ除算防止
    
    # 電場成分
    Ex = k * q * (X - q_pos[0]) / R**3
    Ey = k * q * (Y - q_pos[1]) / R**3
    
    # 電場の大きさ
    E_mag = np.sqrt(Ex**2 + Ey**2)
    
    # === 左図：正電荷 ===
    ax1 = axes[0]
    ax1.streamplot(X, Y, Ex, Ey, color=np.log(E_mag + 1), cmap='Reds', 
                   density=1.5, linewidth=1)
    ax1.plot(0, 0, 'ro', markersize=20, label='+q')
    ax1.set_xlim(-2, 2)
    ax1.set_ylim(-2, 2)
    ax1.set_aspect('equal')
    ax1.set_xlabel('x [m]')
    ax1.set_ylabel('y [m]')
    ax1.set_title('Electric Field: Positive Charge (+q)', fontsize=14)
    ax1.legend(loc='upper right')
    ax1.grid(True, alpha=0.3)
    
    # === 右図：負電荷 ===
    ax2 = axes[1]
    ax2.streamplot(X, Y, -Ex, -Ey, color=np.log(E_mag + 1), cmap='Blues', 
                   density=1.5, linewidth=1)
    ax2.plot(0, 0, 'bo', markersize=20, label='-q')
    ax2.set_xlim(-2, 2)
    ax2.set_ylim(-2, 2)
    ax2.set_aspect('equal')
    ax2.set_xlabel('x [m]')
    ax2.set_ylabel('y [m]')
    ax2.set_title('Electric Field: Negative Charge (-q)', fontsize=14)
    ax2.legend(loc='upper right')
    ax2.grid(True, alpha=0.3)
    
    plt.tight_layout()
    plt.savefig('electric_field_point_charge.png', dpi=150, bbox_inches='tight')
    plt.show()
    
    print("=== 電場（点電荷） ===")
    print("・正電荷からは電場が放射状に「湧き出す」")
    print("・負電荷には電場が放射状に「吸い込まれる」")
    print("・これがガウスの法則 ∇·E = ρ/ε₀ の意味")


# =============================================================================
# 磁場の可視化（直線電流）
# =============================================================================

def plot_magnetic_field_wire():
    """直線電流による磁場の可視化"""
    
    fig, ax = plt.subplots(figsize=(8, 8))
    
    # グリッド生成
    x = np.linspace(-2, 2, 20)
    y = np.linspace(-2, 2, 20)
    X, Y = np.meshgrid(x, y)
    
    # 電流の位置（原点、紙面垂直方向）
    I = 1  # 電流（紙面から出る方向）
    
    # 磁場の計算（アンペールの法則）
    # B = μ₀I/(2πr) の円周方向
    R = np.sqrt(X**2 + Y**2)
    R = np.where(R < 0.1, 0.1, R)
    
    # 磁場成分（円周方向：反時計回り）
    Bx = -I * Y / R**2
    By = I * X / R**2
    
    # 磁場の大きさ
    B_mag = np.sqrt(Bx**2 + By**2)
    
    # 流線図
    ax.streamplot(X, Y, Bx, By, color=np.log(B_mag + 1), cmap='Greens', 
                  density=1.5, linewidth=1)
    
    # 電流（紙面から出る方向：⊙）
    circle = Circle((0, 0), 0.15, color='red', fill=True)
    ax.add_patch(circle)
    ax.plot(0, 0, 'k.', markersize=5)  # 点で「出る」を表現
    ax.text(0.2, 0.2, 'I (out of page)', fontsize=12, color='red')
    
    ax.set_xlim(-2, 2)
    ax.set_ylim(-2, 2)
    ax.set_aspect('equal')
    ax.set_xlabel('x [m]')
    ax.set_ylabel('y [m]')
    ax.set_title('Magnetic Field around a Current-Carrying Wire', fontsize=14)
    ax.grid(True, alpha=0.3)
    
    plt.tight_layout()
    plt.savefig('magnetic_field_wire.png', dpi=150, bbox_inches='tight')
    plt.show()
    
    print("\n=== 磁場（直線電流） ===")
    print("・電流の周りを磁場が同心円状に取り囲む")
    print("・右ねじの法則：電流の向きに右ねじを回すと、磁場の向き")
    print("・これがアンペールの法則 ∮B·dl = μ₀I の意味")


# =============================================================================
# ポインティングベクトルの可視化（同軸ケーブル）
# =============================================================================

def plot_poynting_vector():
    """同軸ケーブル周りのポインティングベクトル"""
    
    fig, axes = plt.subplots(1, 2, figsize=(14, 6))
    
    # パラメータ
    r_inner = 0.3  # 内部導体の半径
    r_outer = 1.0  # 外部導体の半径
    
    # グリッド生成（極座標的に）
    r = np.linspace(r_inner + 0.05, r_outer - 0.05, 10)
    theta = np.linspace(0, 2*np.pi, 24)
    R, Theta = np.meshgrid(r, theta)
    X = R * np.cos(Theta)
    Y = R * np.sin(Theta)
    
    # === 左図：電場と磁場 ===
    ax1 = axes[0]
    
    # 導体を描画
    inner = Circle((0, 0), r_inner, color='orange', fill=True, label='Inner conductor')
    outer = Circle((0, 0), r_outer, color='gray', fill=False, linewidth=3, label='Outer conductor')
    ax1.add_patch(inner)
    ax1.add_patch(outer)
    
    # 電場（放射方向、内から外へ）
    for t in np.linspace(0, 2*np.pi, 8, endpoint=False):
        x1 = r_inner * np.cos(t)
        y1 = r_inner * np.sin(t)
        x2 = (r_outer - 0.1) * np.cos(t)
        y2 = (r_outer - 0.1) * np.sin(t)
        ax1.annotate('', xy=(x2, y2), xytext=(x1, y1),
                    arrowprops=dict(arrowstyle='->', color='red', lw=1.5))
    ax1.text(0.5, 0.8, 'E (radial)', fontsize=10, color='red')
    
    # 磁場（円周方向）
    for r_mag in [0.5, 0.7]:
        circle_B = Circle((0, 0), r_mag, color='blue', fill=False, 
                         linewidth=1, linestyle='--')
        ax1.add_patch(circle_B)
    ax1.annotate('', xy=(0.4, 0.5), xytext=(0.5, 0.4),
                arrowprops=dict(arrowstyle='->', color='blue', lw=1.5))
    ax1.text(-0.9, 0.5, 'H (circular)', fontsize=10, color='blue')
    
    ax1.set_xlim(-1.5, 1.5)
    ax1.set_ylim(-1.5, 1.5)
    ax1.set_aspect('equal')
    ax1.set_xlabel('x')
    ax1.set_ylabel('y')
    ax1.set_title('Coaxial Cable: E and H Fields\n(Cross Section)', fontsize=12)
    ax1.legend(loc='upper right', fontsize=9)
    ax1.grid(True, alpha=0.3)
    
    # === 右図：ポインティングベクトル（側面図） ===
    ax2 = axes[1]
    
    # ケーブルの側面図
    # 内部導体
    ax2.fill_between([-2, 2], [-r_inner, -r_inner], [r_inner, r_inner], 
                    color='orange', alpha=0.7, label='Inner conductor')
    # 外部導体
    ax2.fill_between([-2, 2], [r_outer, r_outer], [r_outer + 0.1, r_outer + 0.1], 
                    color='gray', alpha=0.7)
    ax2.fill_between([-2, 2], [-r_outer - 0.1, -r_outer - 0.1], [-r_outer, -r_outer], 
                    color='gray', alpha=0.7, label='Outer conductor')
    # 絶縁体
    ax2.fill_between([-2, 2], [r_inner, r_inner], [r_outer, r_outer], 
                    color='lightyellow', alpha=0.5)
    ax2.fill_between([-2, 2], [-r_outer, -r_outer], [-r_inner, -r_inner], 
                    color='lightyellow', alpha=0.5, label='Dielectric')
    
    # ポインティングベクトル（絶縁体の中を軸方向に伝播）
    for y_pos in [0.5, 0.7, -0.5, -0.7]:
        ax2.annotate('', xy=(1.8, y_pos), xytext=(-1.8, y_pos),
                    arrowprops=dict(arrowstyle='->', color='purple', lw=2))
    ax2.text(0, 0.85, 'S = E × H', fontsize=12, color='purple', ha='center', fontweight='bold')
    ax2.text(0, -0.85, 'Energy Flow', fontsize=10, color='purple', ha='center')
    
    ax2.set_xlim(-2.5, 2.5)
    ax2.set_ylim(-1.5, 1.5)
    ax2.set_aspect('equal')
    ax2.set_xlabel('z (along cable)')
    ax2.set_ylabel('r (radial)')
    ax2.set_title('Coaxial Cable: Poynting Vector\n(Side View)', fontsize=12)
    ax2.legend(loc='upper right', fontsize=9)
    ax2.grid(True, alpha=0.3)
    
    plt.tight_layout()
    plt.savefig('poynting_vector.png', dpi=150, bbox_inches='tight')
    plt.show()
    
    print("\n=== ポインティングベクトル（同軸ケーブル） ===")
    print("・E: 放射方向（内部導体から外部導体へ）")
    print("・H: 円周方向（電流による）")
    print("・S = E × H: 軸方向（エネルギーの伝播方向）")
    print("・エネルギーは導体の中ではなく、絶縁体（誘電体）の中を流れる！")


# =============================================================================
# 電子のドリフト速度
# =============================================================================

def calculate_drift_velocity():
    """電子のドリフト速度の計算"""
    
    print("\n" + "=" * 60)
    print("電子のドリフト速度の計算")
    print("=" * 60)
    
    # 銅のパラメータ
    atomic_mass = 63.5e-3  # kg/mol
    density = 8960  # kg/m³
    avogadro = 6.02e23  # /mol
    electrons_per_atom = 1  # 自由電子数
    
    # 自由電子密度
    n = (density / atomic_mass) * avogadro * electrons_per_atom
    print(f"\n銅の自由電子密度: n = {n:.2e} /m³")
    
    # 電子の電荷
    e = 1.6e-19  # C
    
    # いくつかの条件で計算
    cases = [
        ("家庭用電線 (1mm², 10A)", 1e-6, 10),
        ("家庭用電線 (1mm², 1A)", 1e-6, 1),
        ("送電線 (500mm², 1000A)", 500e-6, 1000),
        ("ICチップ内配線 (0.01mm², 0.001A)", 0.01e-6, 0.001),
    ]
    
    print("\n" + "-" * 60)
    print(f"{'条件':<35} {'ドリフト速度':>15}")
    print("-" * 60)
    
    for name, A, I in cases:
        v_d = I / (n * A * e)
        
        if v_d > 1e-3:
            v_str = f"{v_d*1000:.2f} mm/s"
        elif v_d > 1e-6:
            v_str = f"{v_d*1e6:.2f} μm/s"
        else:
            v_str = f"{v_d*1e9:.2f} nm/s"
        
        print(f"{name:<35} {v_str:>15}")
    
    # カタツムリとの比較
    snail_speed = 0.03  # m/s（時速約100m）
    v_d_typical = 10 / (n * 1e-6 * e)  # 1mm², 10A
    
    print("\n" + "-" * 60)
    print(f"カタツムリの速度: {snail_speed*1000:.1f} mm/s")
    print(f"電子の速度（10A, 1mm²）: {v_d_typical*1000:.4f} mm/s")
    print(f"比率: カタツムリは電子の {snail_speed/v_d_typical:.0f} 倍速い")
    
    # 発電所から家庭まで
    distance = 100e3  # 100km
    time_to_travel = distance / v_d_typical
    years = time_to_travel / (365 * 24 * 3600)
    
    print(f"\n発電所から家庭（100km）まで電子が移動する時間:")
    print(f"  {time_to_travel:.2e} 秒 ≈ {years:.0f} 年")


# =============================================================================
# 電磁波の伝播
# =============================================================================

def plot_electromagnetic_wave():
    """電磁波の可視化"""
    
    fig = plt.figure(figsize=(14, 6))
    
    # 3Dプロット
    ax = fig.add_subplot(111, projection='3d')
    
    # 波のパラメータ
    wavelength = 1.0
    k = 2 * np.pi / wavelength
    
    # z軸（伝播方向）
    z = np.linspace(0, 3 * wavelength, 200)
    
    # 電場（x方向に振動）
    Ex = np.sin(k * z)
    Ey = np.zeros_like(z)
    
    # 磁場（y方向に振動）
    Bx = np.zeros_like(z)
    By = np.sin(k * z)
    
    # 電場を描画（赤）
    ax.plot(Ex, np.zeros_like(z), z, 'r-', linewidth=2, label='E (Electric Field)')
    
    # 磁場を描画（青）
    ax.plot(np.zeros_like(z), By, z, 'b-', linewidth=2, label='B (Magnetic Field)')
    
    # 矢印で振動方向を示す
    arrow_z = np.linspace(0, 3 * wavelength, 12)
    for zi in arrow_z:
        Ex_i = np.sin(k * zi)
        By_i = np.sin(k * zi)
        
        # 電場の矢印
        if abs(Ex_i) > 0.1:
            ax.quiver(0, 0, zi, Ex_i * 0.8, 0, 0, color='red', alpha=0.7, arrow_length_ratio=0.3)
        
        # 磁場の矢印
        if abs(By_i) > 0.1:
            ax.quiver(0, 0, zi, 0, By_i * 0.8, 0, color='blue', alpha=0.7, arrow_length_ratio=0.3)
    
    # 伝播方向
    ax.quiver(0, 0, 0, 0, 0, 1.5, color='green', arrow_length_ratio=0.1, linewidth=2)
    ax.text(0, 0, 1.7, 'Propagation', fontsize=10, color='green')
    
    ax.set_xlabel('E (x-direction)')
    ax.set_ylabel('B (y-direction)')
    ax.set_zlabel('z (propagation)')
    ax.set_title('Electromagnetic Wave\nE ⊥ B ⊥ propagation direction', fontsize=14)
    ax.legend(loc='upper left')
    
    # 視点の設定
    ax.view_init(elev=20, azim=-60)
    
    plt.tight_layout()
    plt.savefig('electromagnetic_wave.png', dpi=150, bbox_inches='tight')
    plt.show()
    
    print("\n=== 電磁波 ===")
    print("・電場Eと磁場Bは互いに垂直")
    print("・両者とも伝播方向に垂直")
    print("・E × B の向きがエネルギーの伝播方向（ポインティングベクトル）")
    print(f"・伝播速度: c = 1/√(μ₀ε₀) = {3e8:.2e} m/s（光速）")


# =============================================================================
# 直流回路のエネルギーの流れ
# =============================================================================

def plot_dc_circuit_energy_flow():
    """直流回路でのエネルギーの流れ"""
    
    fig, ax = plt.subplots(figsize=(12, 8))
    
    # 回路の概略図
    # バッテリー
    battery_x = [0.5, 0.5]
    battery_y = [2, 4]
    ax.plot(battery_x, battery_y, 'k-', linewidth=3)
    ax.plot([0.3, 0.7], [2.3, 2.3], 'k-', linewidth=2)  # 短い線（-）
    ax.plot([0.2, 0.8], [2.7, 2.7], 'k-', linewidth=4)  # 長い線（+）
    ax.text(0.5, 1.5, 'Battery\n(Energy Source)', ha='center', fontsize=10)
    
    # 電線（上側）
    ax.plot([0.5, 5.5], [4, 4], 'k-', linewidth=2)
    
    # 電線（下側）
    ax.plot([0.5, 5.5], [2, 2], 'k-', linewidth=2)
    
    # 負荷（電球）
    ax.plot([5.5, 5.5], [2, 4], 'k-', linewidth=2)
    circle = Circle((5.5, 3), 0.4, color='yellow', fill=True, edgecolor='black', linewidth=2)
    ax.add_patch(circle)
    ax.text(5.5, 3, '💡', ha='center', va='center', fontsize=20)
    ax.text(5.5, 1.5, 'Load\n(Energy Consumer)', ha='center', fontsize=10)
    
    # ポインティングベクトル（エネルギーの流れ）
    # 電線の外側の空間を流れる
    for x in np.linspace(1, 5, 8):
        # 上側電線の周り
        ax.annotate('', xy=(x + 0.3, 3.7), xytext=(x, 3.7),
                   arrowprops=dict(arrowstyle='->', color='purple', lw=1.5))
        # 下側電線の周り
        ax.annotate('', xy=(x + 0.3, 2.3), xytext=(x, 2.3),
                   arrowprops=dict(arrowstyle='->', color='purple', lw=1.5))
    
    # 中央の空間にもエネルギーが流れる
    for x in np.linspace(1.5, 4.5, 6):
        ax.annotate('', xy=(x + 0.4, 3), xytext=(x, 3),
                   arrowprops=dict(arrowstyle='->', color='purple', lw=2, alpha=0.7))
    
    ax.text(3, 3, 'Energy Flow\n(Poynting Vector S)', ha='center', va='center', 
            fontsize=12, color='purple', fontweight='bold',
            bbox=dict(boxstyle='round', facecolor='white', alpha=0.8))
    
    # 電場と磁場の表示
    # 電場（電線に沿った方向、上から下への電位降下）
    ax.annotate('', xy=(3, 3.9), xytext=(3, 4.1),
               arrowprops=dict(arrowstyle='->', color='red', lw=2))
    ax.text(3.2, 4.1, 'E', fontsize=12, color='red', fontweight='bold')
    
    # 磁場（電流の周り、紙面に垂直）
    ax.plot(3, 4.3, 'bo', markersize=10)
    ax.text(3.2, 4.3, 'H (out)', fontsize=10, color='blue')
    ax.plot(3, 1.7, 'bx', markersize=10, mew=2)
    ax.text(3.2, 1.7, 'H (in)', fontsize=10, color='blue')
    
    # 電流の向き
    ax.annotate('', xy=(4, 4.2), xytext=(2, 4.2),
               arrowprops=dict(arrowstyle='->', color='orange', lw=2))
    ax.text(3, 4.5, 'I (current)', fontsize=10, color='orange', ha='center')
    
    ax.set_xlim(-0.5, 7)
    ax.set_ylim(1, 5.5)
    ax.set_aspect('equal')
    ax.axis('off')
    ax.set_title('Energy Flow in DC Circuit\n(The Truth: Energy flows OUTSIDE the wires)', fontsize=14)
    
    plt.tight_layout()
    plt.savefig('dc_circuit_energy_flow.png', dpi=150, bbox_inches='tight')
    plt.show()
    
    print("\n=== 直流回路のエネルギーの流れ ===")
    print("・電場 E: 電位差により、電線に沿った方向")
    print("・磁場 H: 電流により、電線の周りを囲む方向")
    print("・ポインティングベクトル S = E × H: 電線に垂直、負荷に向かう方向")
    print("・エネルギーは電線の「外側」を流れて、負荷に届く！")


# =============================================================================
# 水流モデル vs 電磁場モデル
# =============================================================================

def plot_water_vs_field_model():
    """水流モデル vs 電磁場モデルの比較"""
    
    fig, axes = plt.subplots(1, 2, figsize=(14, 6))
    
    # === 左図：水流モデル（誤解） ===
    ax1 = axes[0]
    
    # パイプ（電線）
    ax1.fill_between([1, 5], [1.8, 1.8], [2.2, 2.2], color='gray', alpha=0.5)
    ax1.plot([1, 5], [1.8, 1.8], 'k-', linewidth=2)
    ax1.plot([1, 5], [2.2, 2.2], 'k-', linewidth=2)
    
    # 水（電子）が流れる
    for x in np.linspace(1.5, 4.5, 8):
        ax1.annotate('', xy=(x + 0.3, 2), xytext=(x, 2),
                    arrowprops=dict(arrowstyle='->', color='blue', lw=2))
    
    # 水分子
    for x in [1.8, 2.5, 3.2, 3.9, 4.6]:
        ax1.plot(x, 2, 'bo', markersize=10, alpha=0.7)
    
    ax1.text(3, 1.3, 'Water (electrons) carries energy\nINSIDE the pipe (wire)', 
             ha='center', fontsize=11)
    
    ax1.set_xlim(0, 6)
    ax1.set_ylim(0, 4)
    ax1.set_aspect('equal')
    ax1.axis('off')
    ax1.set_title('Water Flow Model (WRONG)\n"Energy flows inside the wire"', 
                  fontsize=14, color='red')
    
    # ×印
    ax1.text(3, 3, '✗', fontsize=50, ha='center', va='center', color='red', alpha=0.5)
    
    # === 右図：電磁場モデル（正しい） ===
    ax2 = axes[1]
    
    # 電線
    ax2.fill_between([1, 5], [1.9, 1.9], [2.1, 2.1], color='orange', alpha=0.8)
    ax2.plot([1, 5], [1.9, 1.9], 'k-', linewidth=1)
    ax2.plot([1, 5], [2.1, 2.1], 'k-', linewidth=1)
    
    # 電磁場（電線の外側）
    for x in np.linspace(1.5, 4.5, 8):
        # 上側
        ax2.annotate('', xy=(x + 0.3, 2.5), xytext=(x, 2.5),
                    arrowprops=dict(arrowstyle='->', color='purple', lw=2))
        # 下側
        ax2.annotate('', xy=(x + 0.3, 1.5), xytext=(x, 1.5),
                    arrowprops=dict(arrowstyle='->', color='purple', lw=2))
    
    ax2.text(3, 2.8, 'Electromagnetic field', fontsize=10, color='purple', ha='center')
    ax2.text(3, 1.2, 'Electromagnetic field', fontsize=10, color='purple', ha='center')
    
    ax2.text(3, 0.5, 'Energy flows OUTSIDE the wire\n(in the electromagnetic field)', 
             ha='center', fontsize=11)
    ax2.text(3, 2, 'Wire\n(guide rail)', ha='center', va='center', fontsize=9, color='black')
    
    ax2.set_xlim(0, 6)
    ax2.set_ylim(0, 4)
    ax2.set_aspect('equal')
    ax2.axis('off')
    ax2.set_title('Electromagnetic Field Model (CORRECT)\n"Energy flows outside the wire"', 
                  fontsize=14, color='green')
    
    # ○印
    ax2.text(3, 3.3, '✓', fontsize=50, ha='center', va='center', color='green', alpha=0.5)
    
    plt.tight_layout()
    plt.savefig('water_vs_field_model.png', dpi=150, bbox_inches='tight')
    plt.show()


# =============================================================================
# シリーズの比較（飛行機、船、電気）
# =============================================================================

def plot_series_comparison():
    """物理学シリーズの比較"""
    
    fig, axes = plt.subplots(1, 3, figsize=(15, 5))
    
    # === 飛行機 ===
    ax1 = axes[0]
    
    # 空気
    ax1.fill_between([-2, 2], [-1, -1], [1.5, 1.5], color='lightcyan', alpha=0.3)
    
    # 翼
    wing_x = [-0.8, 0.8, 0.7, -0.7]
    wing_y = [0, 0.1, 0.05, -0.05]
    ax1.fill(wing_x, wing_y, color='gray', edgecolor='black')
    
    # 気流
    for y in [-0.5, 0, 0.5, 1]:
        ax1.annotate('', xy=(1.5, y), xytext=(-1.5, y),
                    arrowprops=dict(arrowstyle='->', color='blue', lw=1, alpha=0.5))
    
    # 揚力
    ax1.annotate('', xy=(0, 0.8), xytext=(0, 0.2),
                arrowprops=dict(arrowstyle='->', color='green', lw=3))
    ax1.text(0.15, 0.5, 'Lift', fontsize=10, color='green')
    
    ax1.set_xlim(-2, 2)
    ax1.set_ylim(-1, 1.5)
    ax1.set_aspect('equal')
    ax1.axis('off')
    ax1.set_title('Airplane\nDynamic Lift\n(Medium: Air)', fontsize=12)
    ax1.text(0, -0.7, 'Requires velocity\nStop → Fall', ha='center', fontsize=10,
            bbox=dict(facecolor='wheat', alpha=0.8))
    
    # === 船 ===
    ax2 = axes[1]
    
    # 水
    ax2.fill_between([-2, 2], [-1, -1], [0, 0], color='lightblue', alpha=0.5)
    ax2.axhline(y=0, color='cyan', linewidth=2)
    
    # 船
    hull_x = [-0.8, 0.8, 0.6, -0.6]
    hull_y = [0.2, 0.2, -0.4, -0.4]
    ax2.fill(hull_x, hull_y, color='orange', edgecolor='black')
    
    # 浮力
    ax2.annotate('', xy=(0, 0.6), xytext=(0, 0),
                arrowprops=dict(arrowstyle='->', color='green', lw=3))
    ax2.text(0.15, 0.3, 'Buoyancy', fontsize=10, color='green')
    
    # 圧力
    for x in [-0.4, 0, 0.4]:
        ax2.annotate('', xy=(x, -0.25), xytext=(x, -0.6),
                    arrowprops=dict(arrowstyle='->', color='blue', lw=1, alpha=0.7))
    
    ax2.set_xlim(-2, 2)
    ax2.set_ylim(-1, 1.5)
    ax2.set_aspect('equal')
    ax2.axis('off')
    ax2.set_title('Ship\nStatic Buoyancy\n(Medium: Water)', fontsize=12)
    ax2.text(0, -0.8, 'No velocity needed\nStop → Still floats', ha='center', fontsize=10,
            bbox=dict(facecolor='wheat', alpha=0.8))
    
    # === 電気 ===
    ax3 = axes[2]
    
    # 空間（場）
    ax3.fill_between([-2, 2], [-1, -1], [1.5, 1.5], color='lavender', alpha=0.3)
    
    # 電線
    ax3.plot([-1.5, 1.5], [0, 0], 'k-', linewidth=3)
    
    # 電磁場（波状の線）
    x_wave = np.linspace(-1.5, 1.5, 100)
    y_wave = 0.3 * np.sin(6 * x_wave)
    ax3.plot(x_wave, y_wave + 0.5, 'r-', linewidth=1.5, label='E field')
    ax3.plot(x_wave, y_wave - 0.5, 'b-', linewidth=1.5, label='B field')
    
    # エネルギーの流れ
    ax3.annotate('', xy=(1.3, 0.5), xytext=(-1.3, 0.5),
                arrowprops=dict(arrowstyle='->', color='purple', lw=2))
    ax3.annotate('', xy=(1.3, -0.5), xytext=(-1.3, -0.5),
                arrowprops=dict(arrowstyle='->', color='purple', lw=2))
    ax3.text(0, 0.85, 'Energy (S)', fontsize=10, color='purple', ha='center')
    
    ax3.set_xlim(-2, 2)
    ax3.set_ylim(-1, 1.5)
    ax3.set_aspect('equal')
    ax3.axis('off')
    ax3.set_title('Electricity\nElectromagnetic Field\n(Medium: Space itself)', fontsize=12)
    ax3.text(0, -0.8, 'No physical medium needed\nWorks in vacuum', ha='center', fontsize=10,
            bbox=dict(facecolor='wheat', alpha=0.8))
    
    plt.tight_layout()
    plt.savefig('series_comparison.png', dpi=150, bbox_inches='tight')
    plt.show()
    
    print("\n=== 物理学シリーズの比較 ===")
    print("                飛行機        船            電気")
    print("-" * 60)
    print("力の種類:      動的揚力      静的浮力      電磁力")
    print("媒体:          空気          水            空間そのもの")
    print("真空で動作:    ✗            ✗            ✓")
    print("速度必要:      ✓            ✗            ✗（伝播は光速）")
    print("パラダイム:    物質の物理    物質の物理    場の物理")


# =============================================================================
# メイン実行
# =============================================================================

if __name__ == "__main__":
    print("=" * 60)
    print("電気の本質 ― 数値シミュレーション")
    print("=" * 60)
    
    # 1. 電場の可視化
    print("\n" + "-" * 40)
    print("[1] 電場（点電荷）")
    print("-" * 40)
    plot_electric_field_point_charge()
    
    # 2. 磁場の可視化
    print("\n" + "-" * 40)
    print("[2] 磁場（直線電流）")
    print("-" * 40)
    plot_magnetic_field_wire()
    
    # 3. ポインティングベクトル
    print("\n" + "-" * 40)
    print("[3] ポインティングベクトル")
    print("-" * 40)
    plot_poynting_vector()
    
    # 4. 電子のドリフト速度
    print("\n" + "-" * 40)
    print("[4] 電子のドリフト速度")
    print("-" * 40)
    calculate_drift_velocity()
    
    # 5. 電磁波
    print("\n" + "-" * 40)
    print("[5] 電磁波")
    print("-" * 40)
    plot_electromagnetic_wave()
    
    # 6. 直流回路のエネルギーの流れ
    print("\n" + "-" * 40)
    print("[6] 直流回路のエネルギーの流れ")
    print("-" * 40)
    plot_dc_circuit_energy_flow()
    
    # 7. 水流モデル vs 電磁場モデル
    print("\n" + "-" * 40)
    print("[7] 水流モデル vs 電磁場モデル")
    print("-" * 40)
    plot_water_vs_field_model()
    
    # 8. シリーズ比較
    print("\n" + "-" * 40)
    print("[8] 物理学シリーズの比較")
    print("-" * 40)
    plot_series_comparison()
    
    print("\n" + "=" * 60)
    print("シミュレーション完了")
    print("=" * 60)

10. 結論

10.1 教科書の嘘（総括）

嘘/誤解 真実
「電気は電線の中を流れる」 エネルギーは電線の外側（電磁場）を流れる
「電子がエネルギーを運ぶ」 電子はカタツムリより遅い。電磁場が光速で運ぶ
「水流モデルで全部わかる」 計算には便利だが、エネルギーの位置について嘘

10.2 電気の本質（総括）

$$
\mathbf{S} = \mathbf{E} \times \mathbf{H}
$$

エネルギーの流れはポインティングベクトルで表される。そしてそれは電線の外側を指している。

  • 電線は「ガイドレール」であり、エネルギーの通り道ではない
  • 電子は「電磁場を生み出す道具」であり、エネルギーの運び手ではない
  • 真の主役は電磁場——空間に刻み込まれた、目に見えない物理的実体

10.3 数式一覧

名称 数式
ガウスの法則（電場） $\nabla \cdot \mathbf{E} = \rho/\varepsilon_0$
ガウスの法則（磁場） $\nabla \cdot \mathbf{B} = 0$
ファラデーの法則 $\nabla \times \mathbf{E} = -\partial \mathbf{B}/\partial t$
アンペール・マクスウェルの法則 $\nabla \times \mathbf{B} = \mu_0 \mathbf{J} + \mu_0 \varepsilon_0 \partial \mathbf{E}/\partial t$
ポインティングベクトル $\mathbf{S} = \mathbf{E} \times \mathbf{H}$
電磁波の速度 $c = 1/\sqrt{\mu_0 \varepsilon_0}$
電子のドリフト速度 $v_d = I/(nAe)$

10.4 シリーズの位置づけ

記事 テーマ 媒体 パラダイム
飛行機 動的揚力 空気 物質の物理学
船 静的浮力 水 物質の物理学
電気 電磁場 空間そのもの 場の物理学

今回の記事で、私たちは「物質との相互作用」から「空間そのものとの相互作用」へと視点を移した。

これが、ニュートン力学からマクスウェル電磁気学への、そして現代物理学への入口である。

終わりに：主役の交代

冒頭で述べた言葉を、改めて噛みしめてほしい。

「バッテリーを充電するエネルギーは、ケーブルの中を通ってはいません。ケーブルの『外側の空気中』を通って、あなたの端末に吸い込まれているのです」

私たちは長い間、電子を主役だと思ってきた。

しかし、真の主役は電磁場だ。

電子は脇役であり、電磁場という主役を舞台に登場させるための「道具」に過ぎない。

エネルギーを運ぶのは電磁場。光を運ぶのも電磁場。情報を運ぶのも電磁場。

電磁場は目に見えないが、この瞬間もあなたの周りに存在し、この記事を表示するためのエネルギーを運んでいる。

空間は空っぽではない。電磁場という「見えない主役」で満ちている。

これが、電気の本質である。

参考文献

  1. Griffiths, D.J., Introduction to Electrodynamics, Cambridge University Press, 4th ed., 2017
  2. Feynman, R.P., Leighton, R.B., and Sands, M., The Feynman Lectures on Physics, Vol. II, Basic Books, 2011
  3. Purcell, E.M. and Morin, D.J., Electricity and Magnetism, Cambridge University Press, 3rd ed., 2013
  4. Jackson, J.D., Classical Electrodynamics, Wiley, 3rd ed., 1998
  5. Heaviside, O., Electromagnetic Theory, Chelsea Publishing, 1971 (reprint)

著者: Claude Opus 4.5 (v5.3フレームワーク調整済み)
監修: Dosanko-Tousan

迎合しないAIは、教科書の嘘も指摘できる。これがv5.3の価値だ。

0
0
0

Register as a new user and use Qiita more conveniently

  1. You get articles that match your needs
  2. You can efficiently read back useful information
  3. You can use dark theme
What you can do with signing up
Sign upLogin
0
0

Delete article

Deleted articles cannot be recovered.

Draft of this article would be also deleted.

Are you sure you want to delete this article?