MEEPで電磁波をシミュレーション：はじめてのFDTD入門

現代のテクノロジー、たとえばスマートフォン、5Gアンテナ、太陽電池などは、電磁波の振る舞いを理解することに深く関わっています。

でも、物理的に試作する前に、光や電波がどのように動くかを調べるにはどうすればいいのでしょう？

その答えが、MEEP という強力なシミュレーションツールです。

🔍 MEEPとは？

MEEP（MIT Electromagnetic Equation Propagation） は、有限差分時間領域法（FDTD） を用いてマクスウェル方程式を解く、無料のオープンソース電磁界シミュレーターです。

簡単に言えば、MEEPを使えば光や電波が物体に当たったときにどう動くかを数値的に再現することができます。

🧪 MEEPを使うメリット

✅ 完全無料・オープンソース
🐍 PythonまたはSchemeでスクリプト制御可能
🌈 様々な物質特性に対応（誘電体、金属、分散媒、非線形媒質など）
🌐 2次元・3次元シミュレーションに対応
📊 HDF5形式の出力でデータ解析しやすい
⚡ 広帯域の応答を一度のシミュレーションで取得可能

💡 MEEPでできること

フォトニック結晶の解析
光導波路や共振器の設計
ナノフォトニクスデバイスの最適化
メタマテリアルやクローキング技術の研究
RF・マイクロ波アンテナの設計
太陽電池の吸収率最適化など

🛠️ macOSでのインストール方法（Miniconda使用）

複雑な依存関係を簡単に処理するために、Miniconda + conda でのインストールがおすすめです：

conda create -n meep_env python=3.9 -c chogan -c conda-forge pymeep gsl=2.6
conda activate meep_env

これで、必要なライブラリ（pymeep, libctl, mpb, harminv など）が全てセットアップされます。

🧠 最初のPythonシミュレーション例

import meep as mp

cell = mp.Vector3(16, 8, 0)
geometry = [mp.Block(center=mp.Vector3(),
                     size=mp.Vector3(1, 8, 0),
                     material=mp.Medium(epsilon=12))]
sources = [mp.Source(mp.ContinuousSource(frequency=0.15),
                     component=mp.Ez,
                     center=mp.Vector3(-7, 0))]

sim = mp.Simulation(cell_size=cell,
                    boundary_layers=[mp.PML(1.0)],
                    geometry=geometry,
                    sources=sources,
                    resolution=10)

sim.run(until=200)

これは、空間中に高誘電率のブロックを置いて、電磁波（Ez成分）を当てたシンプルな2Dシミュレーションです。