均質化（Homogenization）アルゴリズム テンソル

メタマテリアルの設計において、この「均質化（Homogenization）」は、微細な人工構造（メタ原子）の電磁気的応答を、回路設計や電磁界シミュレータで扱いやすい「1つのマクロな物質特性（実効テンソル誘電率・透磁率）」へと翻訳するための極めて重要な数理的・物理的プロセスです。

スプリットリング共振器（SRR）やワイヤなどの複雑な金属形状をそのままマクロな空間に敷き詰めてシミュレーションしようとすると、メッシュ数が爆発し、計算が実質不可能になります。そこで、構造を「均質な（中身の詰まった）1つの異方性材料」として置き換える必要があり、そのために均質化アルゴリズムが使われます。

VNAの測定データ（複素 S パラメータ）やフルウェーブシミュレーションから、実効的なテンソル（Effective Tensor）を導き出す代表的なアルゴリズムとそのメカニズムを解説します。

1. 均質化の基本思想：有効メディア近似（EMA）

均質化が成り立つ大前提は、「構造の周期（p）が、伝搬する電磁波の波長（λ）よりも十分に小さいこと（p ＜ λ）」です。

電磁波から見れば、微細な金属の切れ込みやワイヤの1つ1つを個別に認識することはできず、それらが一体となった「1つの連続的なメディア（有効メディア）」として感じられます。均質化アルゴリズムは、この「電磁波が感じる見かけの応答」を、マクスウェル方程式に適合するテンソル形式（ $\bar{\bar{\epsilon}}_{eff}, \bar{\bar{\mu}}_{eff}$ ）に落とし込む作業です。

2. 代表的な均質化アルゴリズムとアプローチ

① Sパラメータ反転法（Nicholson-Ross-Weir: NRW法の拡張）

VNAの測定や、単位セル（Unit Cell）のシミュレーションから得られる複素 $S$ パラメータ（$S_{11}, S_{21}$）から直接逆算する、最も実用的な手法です。

• 等方性の基本NRW法: 透過・反射特性から、インピーダンス $Z$ と伝搬定数 $k$ を求め、 $\epsilon = k/Z$ 、 $\mu = kZ$ を導出します。

• テンソルへの拡張（Anisotropic NRW）: 材料を傾けた斜め入射（TE/TMモード）や、偏波を回転させた複数の $S$ パラメータデータセット（ $S_{vv}, S_{hh}, S_{vh}$ ）をインプットします。それぞれの偏波方向における位相変化と減衰量から、方向ごとの有効インピーダンス行列 $\bar{\bar{Z}}_{eff}$ と伝搬行列 $\bar{\bar{n}}_{eff}$ を算出し、最終的に3×3の複素テンソル成分を代数的に分離・抽出します。

② フィールド平均化法（Field Averaging Method）

単位セル（1周期分）の内部における電磁界（$\boldsymbol{E}, \boldsymbol{H}, \boldsymbol{D}, \boldsymbol{B}$）の空間マクロ平均をとる手法です。フルウェーブシミュレータ（HFSSやCSTなど）の内部データを均質化する際によく用いられます。

周期構造の内部では、金属のエッジなどで電磁界が局所的に激しく不均一（微視的電磁界）になります。これを空間的に積分（平均化）して、マクロな量を定義します。

この方法は、構造が波長に対してやや大きくなり、NRW法がマルチモードの発生によって破綻し始める境界領域でも、物理的に妥当なテンソル値を抽出しやすいというメリットがあります。

③ 漸近展開均質化法（Asymptotic Homhomogenization）

周期構造の微細スケール（$l$）とマクロスケール（$L$）の比 $\delta = l/L \to 0$ という極限を利用した、高度な数学的アプローチです。

マクスウェル方程式の電磁界関数を、微視的座標と巨視的座標の2つの変数に分離し、 $\delta$ のべき級数として展開（摂動法）します。構造の静電幾何学的な形状（境界条件）を解くことで、周波数に依存しない物質固有の「有効テンソルカーネル」を厳密に導出できます。

3. 均質化における「物理的障壁」と注意点

メタマテリアルを均質化してテンソルに落とし込む際、通常の自然界の材料測定にはない特有の現象（非局所性）に直面します。

• 空間分散（Spatial Dispersion）と空間非局所性:

  構造のサイズが波長に対して完全に無視できない大きさ（例えば $\lambda/10$ 程度）になってくると、「ある点の分極が、その点の電界だけでなく、隣のセルの電界からも影響を受ける」という現象が起きます。数理的には、テンソル成分が周波数だけでなく波数ベクトル（波の進む方向 $\boldsymbol{k}$）の関数になってしまい、均質化されたテンソルが「入射角によって変化してしまう」という問題が生じます。

• ビアン異方性（Bianisotropy）の発生:

  スプリットリング共振器（SRR）のように構造が上下非対称な場合、「電界をかけると（誘電分極だけでなく）磁気分極も発生する」「磁界をかけると電気分極が発生する」という、電気と磁気の交差結合が起こります。この場合、誘電率・透磁率テンソルだけでは記述できず、相互作用を表す別のテンソル（結合テンソル $\bar{\bar{\xi}}, \bar{\bar{\zeta}}$）を含めた、より高階な構成方程式が必要になります。

まとめ

「均質化アルゴリズムによって実効的なテンソルへと落とし込む」という作業は、複雑なマイクロ・ナノ構造が持つ電磁気的なポテンシャルを、「3×3の複素行列」というエレガントなマクロ指標に凝縮するプロセスです。

VNAで取得した精密な複素 S パラメータ波形は、これらのアルゴリズム（反転解法や最適化フィッティング）を経ることで、初めて「このメタマテリアル基板は、x 軸方向の ε_r = -2.1 + j0.05、厚み方向の ε_r = 3.2 + j0.01$ の異方性物質として振る舞っている」という、設計者が直接扱える数値へと昇華されます。