屈折率の第一原理計算とFDTDの連携によるメタマテリアルの磁気的な光学応答解析#

Advance/PHASE Logo

メタマテリアル（Metamaterials）は、光の波長よりも小さいナノ構造を人工的に配列することで、負の屈折率や完全レンズなど、自然界の物質にはない特異な光学特性を実現する技術です。特に、本来は非磁性である誘電体材料を用いて、光周波数帯域で人工的に「磁気的な応答（透磁率制御）」を引き出す技術は、メタマテリアル設計の核心となります。これらの特性は主に構造の形状に依存しますが、設計の基礎となるのは構成材料の正確な光学定数（誘電率・透磁率）です。本事例では、第一原理計算ソフトウェアAdvance/PHASEを用いてシリコン（Si）の精密な屈折率（誘電関数）を算出し、そのデータを電磁界解析（FDTD法）に適用することで、シリコンナノディスクアレイにおける磁気双極子共鳴（Magnetic Dipole Resonance）を再現した事例を紹介します。

Keywords: 第一原理計算（DFT）, FDTD法, MEEP, メタマテリアル, 磁気双極子共鳴, 人工磁性, シリコンナノディスク, Yee格子

解析モデルと計算手法#

1. 解析モデル：シリコンナノディスクアレイ#

本解析で対象とするモデルを図1に示します。シリコンナノディスク（半径 $R$ ）が周期 $P$ で一次元配列された構造（メタサーフェスの一部）を想定しています。左側から平面波が入射し、ディスク内部で励起される電磁場モードを解析します。第一原理計算（DFT）で得られたシリコンの物性値は、このナノディスクの材料定義として反映されています。算出されたDFTデータを、PythonスクリプトでPyMEEPを介してFDTDシミュレータMEEP [1] に入力して、シリコンナノディスクアレイの光学応答計算を実行します。

シリコンナノディスクアレイの解析モデル模式図
図1. シリコンナノディスクアレイの解析モデル。中央が解析領域（ユニットセル）であり、Y方向の周期境界条件により無限アレイ構造（上下の円）を再現しています。

2. FDTD法の理論基礎と離散化#

FDTD法（Finite-Difference Time-Domain method：有限差分時間領域法）は、マクスウェルの方程式を空間および時間領域で直接解く手法です [1]。電場 $\mathbf{E}$ と磁場 $\mathbf{H}$ を空間的・時間的に半ステップずらして配置するYee格子（Yee Lattice）[2] を採用することで、以下のマクスウェル方程式（回転）を差分方程式として効率的に計算します。 $\frac{\partial \mathbf{H}}{\partial t} = -\frac{1}{\mu} \nabla \times \mathbf{E}$

$\frac{\partial \mathbf{E}}{\partial t} = \frac{1}{\varepsilon} \nabla \times \mathbf{H} - \frac{\sigma}{\varepsilon}\mathbf{E}$

具体的には、例えば2次元TE（Transverse Electric）偏光（電場 $E_x, E_y$ 、磁場 $H_z$ ）の場合、磁場 $H_z$ の時間更新式は以下のように離散化されます。

$H_z^{n+\frac{1}{2}}(i, j) = H_z^{n-\frac{1}{2}}(i, j) - \frac{\Delta t}{\mu \Delta x} \left[ E_y^n(i+1, j) - E_y^n(i, j) - (E_x^n(i, j+1) - E_x^n(i, j)) \right]$

ここで、 $(i, j)$ は空間グリッドのインデックス、 $n$ は時間ステップ、 $\Delta x$ は空間刻み幅（解像度）、 $\Delta t$ は時間刻み幅です。この式により、ある瞬間の電場分布から次の瞬間の磁場が計算され、さらにその磁場から次の電場が計算される「リープフロッグ（蛙飛び）法」によって、波の伝搬がシミュレーションされます。

3. 解析パラメータと数理モデルの対応#

表1に示す解析パラメータは、上記FDTDの計算式において以下のように反映されます。

解像度 (Resolution): 計算における空間刻み幅 $\Delta x$ を決定します。Resolution = 120 pixels/µm の場合、 $\Delta x \approx 8.3$ nm となり、微細なナノディスク形状（半径65 nm）を十分な精度（直径に対し約16グリッド）で表現します。
材料データ (DFT): 更新式中の誘電率 $\varepsilon$ として使用されます。Advance/PHASEで算出された複素屈折率 $n, k$ から $\varepsilon = (n-ik)^2$ として取り込まれ、媒質中での光の速度と減衰を決定します。
周期境界条件 (Y方向): グリッド端 $j=N_y$ での微分計算において、反対側の端 $j=0$ の値を参照することで、無限に続くアレイ構造を数学的に再現します。

表1. 解析パラメータ

カテゴリー	パラメーター	設定値と物理的意味
材料データ (DFT算出)	対象材料	シリコン (Si)
	ターゲット波長	650 nm (屈折率のDFT計算値を入力)
FDTD構造設定	ナノディスク半径 ( $R$ )	65 nm
	アレイ周期 ( $P$ )	400 nm (Y方向周期境界)
	偏光	TE偏光 ( $H_z$ 成分入射により磁気モードを励起)
	解像度 (Resolution)	120 pixels/µm ( $\Delta x \approx 8.3$ nm)

解析結果と考察#

1. 第一原理計算による誘電関数の算出#

図2に、Advance/PHASEで計算されたシリコンの誘電関数（実部 $\varepsilon_1$ 、虚部 $\varepsilon_2$ ）を示します。標準的なDFT計算ではバンドギャップが実験値（約1.12 eV）よりも小さく見積もられる傾向がありますが、Scissor Operator法を適用することで吸収端のエネルギー位置を補正し、実験値とよく一致するスペクトルを得ました。この結果から、波長 650 nm（約1.9 eV）において、シリコンは高い屈折率（ $n \approx 3.86$ ）を持ちながら、消衰係数 $k$ は十分に小さい（ほぼゼロ）であることが確認されました。これは、光を熱に変えずに散乱・共鳴させるメタマテリアル材料として理想的な特性です。

DFT計算によるSiの誘電関数
図2. 第一原理計算で得られたシリコンの誘電関数スペクトル

2. FDTDによる磁気双極子共鳴の可視化#

DFT由来の材料データを適用し、高解像度FDTDシミュレーションを実施した結果を図3に示します。

FDTDシミュレーション結果：磁場分布 FDTDシミュレーション結果：電場ベクトル
図3. (上) 磁場分布 ( $H_z$ )。ディスク中心（赤色）に強い磁場が局在している。(下) 電場ベクトル分布。磁場を取り囲むような変位電流のループが確認できます。

磁場分布 ( $H_z$ ) の考察：

図3（上）を見ると、ナノディスクの中心に非常に強い磁場（赤色の領域）が局在していることが分かります。通常、シリコンなどの誘電体（透磁率 $\mu=1$ ）は磁気的な性質を持ちませんが、ナノ構造化することによって光の磁場成分と強く相互作用し、あたかも「人工的な磁石（磁気双極子）」のように振る舞っています [3]。

電場ベクトル ( $E$ ) の考察：

図3（下）の流線図では、中心の磁場を取り囲むように、電場ベクトルが綺麗な渦（ループ）を描いている様子が確認できます。これは電磁気学におけるアンペールの法則 $\nabla \times \mathbf{H} = \varepsilon \frac{\partial \mathbf{E}}{\partial t}$ に対応する現象であり、誘電体内部に誘起された「環状の変位電流」が、中心の磁場双極子モーメントを生成していることを物理的に裏付けています。この現象はMie共振（Mie Resonance）の一種であり、誘電体メタマテリアルの基本原理となる重要なモードです。

大規模解析への拡張と展望#

本事例では、FDTD (MEEP) との連携をPythonスクリプトによる自動化で行い、物理現象の基礎的な理解やプロトタイピング（Prototyping）に最適です。しかし、3次元構造やより大規模な計算においては、計算負荷が大幅に増大し、難易度も上がります。産業用途での複雑な3次元デバイス設計には、GUI環境やCADインターフェースが充実した商用FDTDソフトウェア（Advance/ParallelWave等）は、本事例と同様のデータフローで連携可能であり、より高度な産業課題の解決に寄与します。

まとめ#

本事例では、Advance/PHASEによる第一原理計算とFDTD解析をシームレスに連携させることで、デバイス動作波長（650 nm）におけるシリコンの高精度な複素屈折率および低損失性（ $k \approx 0$ ）を確認するとともに、ナノディスクアレイ内部の電磁場分布を詳細に解析しました。特に、高解像度シミュレーションによって「磁気双極子共鳴」に伴う渦構造を鮮明に可視化し、人工磁性の発現機構を明らかにすることに成功しました。この「DFT × FDTD」の連携解析フローは、実験データが存在しない新材料を用いたメタマテリアル探索や、ナノスケールでの物性変化を考慮した高精度なデバイス設計において、強力なツールとなります。

本解析の詳細や、研究への適用可能性に関するご相談はこちら

お問い合わせ

参考文献#

A. F. Oskooi, D. Roundy, M. Ibanescu, P. Bermel, J. D. Joannopoulos, and S. G. Johnson, "Meep: A flexible free-software package for electromagnetic simulations by the FDTD method". Comput. Phys. Commun. 181, 687 (2010).
K. S. Yee, "Numerical solution of initial boundary value problems involving Maxwell's equations in isotropic media", IEEE Trans. Antennas Propag. 14, 302 (1966).
A. I. Kuznetsov, A. E. Miroshnichenko, M. L. Brongersma, Y. S. Kivshar, and B. Luk’yanchuk, "Optically resonant dielectric nanostructures", Science 354, aag2472 (2016).