A Tight-binding Approach for Computing Subwavelength Guided Modes in Crystals with Line Defects

技術要約：結晶内のサブ波長導波モードを計算するためのタイトバインディング・アプローチ

問題の定式化
本論文は、非コンパクトな線欠陥を持つ高コントラスト周期媒体における、サブ波長導波モードを決定するという計算上の課題に取り組んでいる。物理系は、密度コントラスト $\delta$ が小さい（$\delta \ll 1$）高コントラストパラメータを持つヘルムホルツ方程式によってモデル化される。このようなレジームでは、単位セルのサイズよりも長い波長に対してバンドギャップが存在する。周期的な共振器配列に線欠陥を導入すると、欠陥付近に局在するモードに対応するサブ波長周波数のバンドが生じる。

数学的に、この問題は、非有界領域上の演算子 $A_\delta = -\kappa \nabla \cdot (\rho^{-1} \nabla \cdot)$ に関する固有値問題として定式化される。著者らは、導波モード（横方向のファイバーに対して平方可積分なもの）および周波数 $\omega = O(\sqrt{\delta})$ に焦点を当てている。標準的な数値的手法（有限要素法 (FEM) や境界要素法 (BEM) など）は、このレジームにおいて、材料界面の解像や結果として生じる大規模な離散系が必要となるため、計算コストが高くなる。

手法
著者らは、キャパシタンス行列定式化から導出されたタイトバインディング近似に基づくフレームワークを提案している。その手法は以下のステップで進行する：

1. 層ポテンシャル理論： 連続的なヘルムホルツ問題を、非有界領域上の単層ポテンシャルを用いて再定式化する。これにより、包含物の境界上の密度に関する境界積分方程式の系へと問題が還元される。
2. キャパシタンス行列の定義： 著者らは、完全なキャパシタンス行列 $C$ および準周期キャパシタンス行列 $\hat{C}^\alpha$ （$\alpha$ はブロッホ準運動量）を定義する。これらの行列は、各包含物上でポテンシャルが一定となる調和境界値問題を解くことで得られる。
3. 打ち切りと減衰解析： 本手法の中核は、これらのキャパシタンス行列の非対角要素が指数関数的に減衰することを証明することにある。この特性により、無限の稠密行列を帯行列（有限範囲の行列）として近似することが可能となり、実質的にタイトバインディング・モデルが構築される。
4. スペクトルの近接性： 小さな $\omega$ に対する層ポテンシャル演算子の漸近展開を用いることで、連続的な固有値問題と、キャパシタンス行列および質量行列を含む離散的な一般化固有値問題との間の定量的な関連性を確立する。

主な貢献
本論文には、以下の3つの理論的および計算的な貢献がある：

1. キャパシタンス行列要素の指数減衰：
   著者らは（定理2.1および2.2）、完全および準周期キャパシタンス行列の両方の非対角要素が、対角成分から離れるにつれて指数関数的に減衰することを証明している。具体的には、$|C_{m,n}^{p,q}| \leq C \min\{\rho^{|m-p|}, \rho^{|n-q|}\}$ （ここで $\rho \in (0,1)$）が成り立つ。これは、半空間およびストリップ上に定義された補助的な調和関数を用いた打ち切り法、ならびに最大原理とハルナックの不等式を組み合わせることで確立されている。この減衰は、無限行列を帯形式に打ち切ることの厳密な正当化を与える。

2. サブ波長周波数の漸近展開：
   固定されたブロッホパラメータ $\alpha \neq 0$ に対して、著者らは導波周波数の漸近公式を導出する：
   $$\omega(\alpha) = \mu_1 \sqrt{\delta \lambda_\alpha} + O(\delta)$$
   ここで $\lambda_\alpha$ は準周期キャパシタンス行列 $\hat{C}^\alpha$ の一般化固有値である（定理5.1）。この結果は、連続演算子と離散行列モデルとの間のスペクトル近接性を確立し、$O(\delta)$ の誤差範囲を与える。

3. 帯行列近似に関する定量的誤差制御：
   系2.2は、完全な準周期キャパシタンス行列とその帯行列への打ち切り（$N$ 個の近傍との相互作用を保持する）との差に関する明示的な誤差界を提供している。誤差は $O(\rho^N)$ で減衰するため、打ち切り誤差と漸近誤差のバランスを取るために $N$ を $\delta$ の関数として選択することが可能となる。

結果

• 理論的： 本論文は、タイトバインディング近似が高コントラスト・サブ波長系（線欠陥を含む）に対して有効であることを厳密に証明している。これは、共振器間の相互作用がキャパシタンス行列表現において短距離であることを示しており、最近接または有限範囲モデルの使用を正当化している。
• 数値的： セクション6では、理論的な減衰率を検証するための数値実験を示す。
  • 減衰の検証： $\hat{C}^\alpha$ の非対角要素の計算により、予測された指数減衰が確認された。
  • 欠陥モード： 線欠陥に対する完全なキャパシタンス行列の固有値を、打ち切り（最近接）行列と比較した。結果として、打ち切り行列は、最大相対誤差約2.6%で欠陥固有値を正確に捉えることが示された。
  • 欠陥間の相互作用： 二量体配置（ダイマー）における2つの欠陥の研究では、欠陥間の距離が増加するにつれて固有値の分裂が指数関数的に減少することが示されており、これはキャパシタンス行列の減衰特性と一致している。
  • ダイマー化された界面： 本手法は、Su-Schrieffer-Heeger (SSH) モデルの2次元版に応用され、打ち切り行列のスペクトルにおける界面モードの候補を正常に特定した。

意義と主張
本論文は、サブ波長導波モードを計算するための計算効率が高く正確なアプローチを提供すると主張しており、高コントラスト極限においてコストがかかる直接的な離散化手法（FEM/BEM）に代わる厳密な選択肢を提示している。

• 厳密な正当化： メタマテリアルにおけるタイトバインディング・モデルの従来のヘリスティックな使用とは異なり、本研究はキャパシタンス行列要素の指数減衰を数学的に証明し、それによって打ち切りによって導入される誤差を定量化している。
• 一般化： 本アプローチは、摂動のない結晶に対する従来の減衰解析を、非コンパクトな線欠陥を持つシステムへと一般化している。
• 限界： 著者らは、サブ波長周波数と $\hat{C}^\alpha$ の固有値との対応関係が $\alpha \neq 0$ について証明されていることを明記している。$\alpha = 0$ の場合は、単層ポテンシャル演算子がその極限において逆行列を持たないため、別途の漸近解析が必要となる。
• 今後の課題： 本手法は曲線の線欠陥や高次元システムへ拡張できる可能性があることを示唆しているが、これらは今後の調査事項として残されている。また、ダイマー化された界面モードのトポロジカルな特性付けも、バルク・シンボルに関するさらなる分析を要する別個の問題として指摘されている。

要約すると、本論文は、連続的な高コントラスト・ヘルムホルツ問題と離散的なタイトバインディング・モデルとの間に数学的に健全な架け橋を築き、誤差を制御しながら欠陥バンドを効率的に計算することを可能にしている。