メタマテリアルとは何ですか？

波長より十分小さい人工の構造（共振器や金属線など）を周期的に並べて、自然の材料にはない電磁的な性質を持たせたものです。電磁波から見ると、細かい構造は平均化されて一様な材料のように見えるので、実効的な誘電率と透磁率で表せます。両方を負にすると、屈折率が負になります。

負の屈折率はどうやって作りますか？

金属の細い線を並べると、ある周波数（プラズマ周波数）より下で誘電率が負になります。分割リング共振器を並べると、共振周波数のすぐ上で透磁率が負になります。両方が負になる周波数帯で屈折率が負になります。この例では透磁率が10〜11.95 GHzで負になり、誘電率は15 GHzより下で負なので、10〜11.95 GHzが負の屈折率の帯域です。

メタマテリアルの解析で注意することは？

単位セルが波長より十分小さい（10分の1程度以下）ことが、実効的な材料とみなせる条件です。この例のセル3 mmは、10 GHzの波長30 mmの約10分の1で、ぎりぎりです。また共振を使うので損失が大きく、共振のすぐ近くでは屈折率の虚部（損失）が大きくなります。実際の構造は電磁界解析（FDTDや有限要素法）で単位セルを計算し、反射と透過の係数から実効定数を求めます。

メタマテリアル電磁解析

カテゴリ: 電磁界解析・高周波 | 2026-10-01 改稿

$Metamaterial unit cell electromagnetic simulation showing negative refractive index and split-ring resonator field distribution$

メタマテリアルのユニットセル解析：SRR構造の電磁場分布と等価物性値抽出の概念図

理論：細かい構造を平均した材料

概要

🙋

メタマテリアルで屈折率が負になるというのは、どういう仕組みですか？

🎓

屈折率は誘電率と透磁率の積の平方根で決まる。普通の材料ではどちらも正だが、金属の細い線を並べると、電子が自由に動けるので、ある周波数より下では電界と逆向きに分極して誘電率が負になる。一方、切れ目の入った小さなリング（分割リング共振器）を並べると、共振のすぐ上の周波数では磁界と逆向きに磁化して透磁率が負になる。両方が同時に負になる周波数では、平方根の符号を負にとるのが物理的に正しく、屈折率が負になる。構造が波長より十分小さければ、波はそれを一様な材料として感じるので、実効的な誘電率と透磁率で扱える。

実効透磁率と実効誘電率のモデル

$$ \mu_{eff}(\omega) = 1 - \frac{F\,\omega^2}{\omega^2 - \omega_0^2 + i\gamma\omega},\qquad \varepsilon_{eff}(\omega) = 1 - \frac{\omega_p^2}{\omega(\omega + i\gamma_w)},\qquad n = \pm\sqrt{\varepsilon_{eff}\,\mu_{eff}} $$

F は分割リングの充填率（0.3）、f₀ は共振周波数（10 GHz）、γ は損失（0.1 GHz相当）、f_p は金属線配列のプラズマ周波数（15 GHz）。屈折率の符号は、損失による減衰が正（Im n > 0）になるようにとる。

負になる帯域と性能指数

$$ \mu_{eff} < 0:\ \ \omega_0 < \omega < \frac{\omega_0}{\sqrt{1 - F}},\qquad \varepsilon_{eff} < 0:\ \ \omega < \omega_p,\qquad \mathrm{FOM} = \frac{|\mathrm{Re}\,n|}{\mathrm{Im}\,n} $$

Coffee Break よもやま話

ヴェセラゴの予言

1968年、ヴェセラゴは誘電率と透磁率がともに負の材料があれば、光が逆向きに屈折するなど不思議な性質を示すことを理論的に予言した。しかし自然にはそのような材料がなく、長い間忘れられていた。1999年にペンドリーが分割リング共振器で負の透磁率を作れることを示し、2001年にスミスらが金属線と組み合わせて、マイクロ波で負の屈折を実験で確かめた。その後、電磁波を物体のまわりに迂回させるクローキングなどへ研究が広がった。

計算例

例：周波数ごとの実効定数と屈折率（単位セル3 mm）

周波数	Re ε	Re μ	n	性能指数	セル/波長
9.0 GHz	−1.78	+2.28	−0.02 + 2.01i	0.0	0.090
10.2 GHz	−1.16	−6.26	−2.72 + 0.40i	6.8	0.102
10.5 GHz	−1.04	−2.19	−1.51 + 0.12i	12.7	0.105
11.0 GHz	−0.86	−0.72	−0.79 + 0.05i	14.9	0.110
11.5 GHz	−0.70	−0.23	−0.40 + 0.04i	10.0	0.115
12.5 GHz	−0.44	+0.17	−0.01 + 0.27i	0.0	0.125
16.0 GHz	+0.12	+0.51	+0.25 + 0.00i	61	0.160

9 GHzと12.5 GHzでは誘電率だけが負で、屈折率はほぼ純虚数になり、波は伝わらずに減衰する。透磁率が負になる帯域は10〜11.95 GHz（損失なしの場合）。

🙋

共振のすぐ上の10.2 GHzが一番大きな負の屈折率なのに、性能指数は11 GHzのほうが高いんですね。

🎓

共振のすぐ近くでは透磁率が大きく負になるが、同時に損失（虚部）も大きくなる。屈折率の大きさだけ見ると10.2 GHzが魅力的だが、波は短い距離で減衰してしまう。性能指数は屈折率の実部と虚部の比で、何波長分伝わるかの目安になる。この例では11 GHz付近で約15と最も高い。実際の設計では、共振器の損失（金属の抵抗や基板の誘電損）を下げることと、使う周波数を共振から少し離すことのかね合いになる。もう一つの注意は、セルの大きさが波長の約10分の1しかないことで、セルがこれより大きくなると、実効的な材料とみなす前提が崩れ、反射と透過から求めた実効定数が入射角などで変わってしまう。

解析の進め方

目標とする周波数帯から、共振周波数とプラズマ周波数を決める。
単位セルを周期境界条件で電磁界解析する（FDTDや有限要素法）。
反射と透過の係数（Sパラメータ）から実効的な誘電率・透磁率を求める。
性能指数と、セルの大きさと波長の比を確かめる。
有限の大きさの試料全体を解析し、端の影響と入射角の依存を確かめる。

Coffee Break よもやま話

「抽出した屈折率が厚さで変わる」

あるメタマテリアルの解析で、1層の単位セルから求めた実効屈折率と、3層重ねたときに求めた値が大きく違った。調べると、セルの大きさが波長の約5分の1あり、層と層の間の結合が強く、1層の値は実効的な材料の性質を表していなかった。セルを小さくして波長の10分の1程度にすると、層の数を変えても実効定数がほぼ一定になった。実効定数を求めるときは、層の数を変えて値が変わらないことを確かめるとよい。

よくある間違い

間違いと対策

間違い	起きること	対策
セルが大きすぎる	実効定数が定まらない	波長の1/10以下
平方根の符号を誤る	屈折率の符号が逆	Im n > 0 で選ぶ
損失を無視	性能を過大評価	性能指数を確認
1層だけで抽出	結合の影響を見逃す	層数を変えて確認
垂直入射だけで評価	斜め入射で性質が変わる	入射角を変える

🙋

関連する内容も知りたいです。

🎓

関連するページはFDTD法、高周波の有限要素法、空洞共振器、誘電体共振器、電磁波の伝搬が参考になる。