Full-Wave Harmonic Balance Framework for Dispersive Time-Varying Photonic Structures

技術要約：分散性を持つ時変フォトニック構造のためのフルウェーブ・ハーモニックバランス・フレームワーク

問題提起
時変フォトニック構造は、時間並進対称性を破ることにより、周波数変換、非可逆性、動的な波面制御といったユニークな機能を提供します。しかし、関与する材料が強い分散性を持つ場合（例：透明導電酸化物、グラフェン、イプシロン近傍ゼロ（ENZ）材料）、これらのシステムを正確にモデリングすることは依然として大きな課題です。このようなプラットフォームでは、時間変調を単なる周波数に依存しないパラメータの瞬時的な摂動として扱うことはできません。代わりに、因果律によって材料応答は励起の履歴に依存することになります。つまり、生成された各フローケ調和波は、変調を通じて他のすべての調和波と結合しつつも、それぞれ独自の材料応答を受けることになります。既存の計算手法にはトレードオフが存在します。時間領域手法（FDTDなど）は定常状態に達するまでに長いシミュレーションを必要とし、因果的かつ分散的な更新の扱いに苦慮します。一方で、特殊な解析的手法や結合波法は、複雑な3D構造、不均質媒体、あるいは任意の変調スペクトルを扱うための幾何学的な柔軟性に欠けることが多くあります。

手法
著者らは、時間進行（time-marching）を回避し、周波数領域で定常状態の電磁応答を直接解く、フルウェーブ・ハーモニックバランス・フレームワークを提案しています。この手法の核となるのは、時変の問題を、拡張されたフローケ空間内における結合された周波数領域の部分問題のセットとして再構成することです。

手法の主要要素は以下の通りです：

• 二次ソース表現： 時間変調成分は、時間依存する境界条件としてではなく、誘導された二次ソースとして表現されます。変調されたバルク媒体には体積分極密度（$P$）が、変調された導電シートには表面電流密度（$J_s$）が用いられます。
• フローケ展開： すべての電場および誘導ソースは、切り捨てられたフローケ級数（$F(r, t) = \sum F_n(r)e^{j\omega_n t}$）へと展開されます（ここで $\omega_n = \omega_0 + n\Omega$）。
• 結合システム定式化： 問題は、グローバルな線形システム $(L - SC)X = Y_{ext}$ として解かれます。
  • $L$ は、材料分散を考慮するために、各調和波の物理的周波数（$\omega_n$）において評価された、静的な背景となるブロック周波数領域のマックスウェル演算子を表します。
  • $S$ は、ソース注入演算子を表します。
  • $C$ は、材料パラメータの時間フーリエ展開から導出されるソース生成演算子です。これらの演算子は、調和空間におけるオフダイアゴナル（非対角）結合を導入し、周波数変換メカニズムを表現します。
• 実装： このフレームワークは有限要素法（FEM）を用いて実装されています。標準的なFEMの未知数ベクトルは、ハーモニック領域のベクトルへと持ち上げられ、これにより各空間的な自由度が複数の調和係数を保持できるようになります。これにより、フルウェーブ・モデリングの幾何学的柔軟性を維持しながら、調和間の結合を捉えることが可能になります。

主な貢献

1. 統一されたフレームワーク： 本手法は、変調された体積媒体と変調された導電シートの両方を扱う単一の代数構造を提供し、任意の時間変調スペクトル（単一トーン、マルチトーン、または非正弦波）をサポートします。
2. 分散の取り扱い： 弱い分散を仮定したりエンベロープ近似を用いたりする手法とは異なり、本フレームワークは、生成された各調和波の特定の周波数において材料分散および放射演算子を評価します。
3. 調和分解された観測量： 解は、すべての保持されたフローケ調和波に対する定常振幅、位相、および電場プロファイルを直接導出するため、時間領域データの後処理を行うことなく、調和固有の散乱、吸収、パワーフロー、および近接場分布の計算を可能にします。

結果
本フレームワークは、2つの代表的な例を通じて検証および実証されました。

1. パラメトリック励起された巻きグラフェンシリンダー：

   • セットアップ： 誘電体シリンダーにグラフェンシートを巻き付け、時間変調されたフェルミ準位（$\Omega = 2\omega_0$）の下で動作する、縮退パラメトリック増幅領域の系。
   • 検証： 調和分解された散乱および吸収に関する数値結果は、円筒形状に対して導出された解析的なフローケ散乱解と極めて良好に一致しました。
   • 知見： シミュレーションは、キャリアからダウンコンバートされた調和波（$n=-1$）へのエネルギー再分配を正常に捉え、共鳴表面プラズモンシステムにおける強い分散と急速な変調をモデル化する本フレームワークの能力を裏付けました。

2. 光時空間メタ表面（ITOベース）：

   • セットアップ： 金の背面に配置された、時間変調されたインジウム錫酸化物（ITO）ENZ層を備えた二進 $TiO_2$ 格子からなる反射型メタ表面。ITOのプラズマ周波数が変調され、誘電率の超高速変化を駆動します。
   • 適用： フレームワークは、周波数変換ビーム偏向のためのメタ表面の設計と最適化に使用されました。目標は、入射光波をサイドバンドへと変換し、特定の空間回折次数（$m=1$）へと偏向させることです。
   • 知見： 最適化された設計は、生成されたサイドバンド（$\omega_{\pm 1}$）を第1空間回折次数へと正常に偏向させました。結果は、ENZ分散、時間依存Drude結合、および空間回折の相互作用を実証し、複雑で現実的な光学デバイスを設計するための本手法の有用性を裏付けました。

意義
本論文は、このフレームワークが分散性を持つ能動的フォトニック構造をモデリングおよび設計するための一般的なプラットフォームを確立すると主張しています。ハーモニックバランス法による物理的な透明性と、フルウェーブFEMによる幾何学的な柔軟性を組み合わせることで、本手法は時変システムにおける既存のソルバーの限界を克服します。これにより、研究者は変調がどのように周波数シフトされたチャネル間で電磁エネルギーを再分配するかを直接定量化することが可能となり、グラフェン、透明導電酸化物、および共鳴メタ表面に基づく次世代の時空間フォトニックデバイスの設計に不可欠な観測量（調和分解されたパワー、吸収、および回折）を提供します。本研究は、このフレームワークを単なる解析ツールとしてだけでなく、周波数変換、パラメトリック増幅、および動的な波面制御を必要とするデバイスの逆設計のためのツールとして位置付けています。