Circuit models for quasi-static bands in metallic photonic crystals

技術要約：金属フォトニック結晶における準静的バンドの回路モデル

問題提起
通常の誘電体フォトニック結晶では、2つの横波バンドがブリルアンゾーン中心（$\Gamma$）からゼロ周波数から線形に現れる。しかし、金属フォトニック結晶（MPC）は、より複雑な低周波スペクトルを示す。孤立した金属包含物は誘電体のように振る舞うが、連結された金属ネットワークは、プラズマ的なカットオフ以下ですべての伝搬モードを除去するか、あるいは追加の準静的バンドや縦波プラズマモード、および有限のブロッホ運動量におけるゼロ周波数の接点（touching sets）をサポートする。既存の理論は、特定の幾何学的形状や材料パラメータに基づく非局所的な構成応答や伝送線路記述を導出しているが、どの準静的バンドが幾何学的な詳細や材料パラメータではなく、導体の配線（接続性）のみによって決定されるかを判断するための一般的な規則を欠いている。したがって、$\Gamma$における準静的な零空間（null spaces）の完全なカウント、有限運動量における接点の基準、および任意の相互接続パターンに対して線形バンドと$\Gamma$二次的な遮蔽バンドを区別する方法を提供する、一般的な理論が必要とされている。

手法
著者らは、細長い無損失金属フォトニック結晶に対し、準静的なワイヤーモデルを用いた配線ベースの回路記述を定式化している。このアプローチは以下の要素を組み合わせている：

1. 回路ラプラシアン定式化： システムは、セル内の結合を記述する容量（$M$）と、ユニットセル間で導体ポテンシャルの置換をエンコードする誘導リンクを備えたLCネットワークとしてモデル化される。誘導ブロッホ・ラプラシアン（$K_L(k)$）が定義され、その零空間がゼロ周波数の解を決定する。
2. 帰還並進部分群（Return-Translation Subgroups）： 周期的な導体グラフに対して、著者らは導体の「並進軌道」を特定する。各軌道に対して、格子ベクトルが物理的な導体を自身へと写す「帰還並進部分群（$H_c$）」を定義する。
3. 置換定理： 動的な並進軌道は、すべての帰還並進にわたって蓄積されるブロッホ位相が1であるとき、回路剛性零ベクトル（circuit stiffness null vector）に寄与する。この条件は、正定値の誘導パラメータとは独立して、ゼロ周波数の接点集合とその多重度を固定する。
4. $\Gamma$における静電カウント： ブリルアンゾーン中心において、物理的なモード数は、ゲージ方向を除去し、一様な平均場セクターを考慮することによって導出される。このカウントは、静電的な導体ポテンシャル形式のランクと、導体の帰還によって短絡されない一様電場方向の次元に依存する。
5. マクスウェル–テレグラフ固有値問題： 1つの導体につき1つの電流チャネルを持つ三周期構造に対し、著者らはマクロな電磁場と導体ポテンシャルおよび電流を結合させる。彼らは、線形バンドと$\Gamma$二次的な遮蔽バンドを区別するための、方向分解されたゼロネットセル電流の基準を導出している。

主要な貢献

• 一般的な配線ルール： 本論文は、ゼロ周波数の接点集合とその多重度が、特定のインダクタンス／キャパシタンスの値やワイヤーの形状（接続性を維持する場合）に依存せず、導体の接続性と帰還並進部分群によって完全に決定されることを確立している。
• 有限運動量における接点： 著者らは、ゼロ周波数の接点が有限のブロッキ運動量において発生することは、導体の接続性によって選択される直接的なブロッホ解であることを示している。これらは非スノモルフィック対称性（スクリュー軸やグライド軸など）やより小さな原始並進を必要とせず、帰還並進の条件から自然に生じる。
• $\Gamma$点におけるモードカウント： 一様電場由来のバンドと遮蔽された相対ポテンシャル由来のバンドを区別する、$\Gamma$における物理的な静的零方向を計算するための公式（式15）が導出されている。
• 分散次数の選択： 遮蔽されたバンドが線形か$\Gamma$二次的かを決定するための基準が提供されている。これは、適合する主要電流が導体方向行列の核（$\ker S$）に含まれるかどうかに依存する。もしそのような電流が存在する場合（ゼロネットセル電流）、バンドは線形となり、そうでなければ$\Gamma$二次的となる。
• 統合： このフレームワークは、ワイヤー束、交差メッシュ、およびインターレースド・メッシュに関する先行研究の結果を、単一の回路論的記述の下に統合している。

結果

• ねじれたワイヤー束： 巡回置換 $\sigma = (123456)$ を持つ6本のワイヤー束について、本モデルは全ブリルアンゾーンにわたって6つのゼロ周波数運動量（$k_z h = 0, \pm \pi/3, \pm 2\pi/3, \pi$）を予測する。フルウェーブシミュレーションはこの接点を裏付けており、モデルは幾何学的な変形（不規則なワイヤー）が置換を維持する限り、ゼロ周波数の運動量は変化せず、有限周波数の分散のみが変化することを示している。
• ねじれたロッド： ねじれた一対のロッドを持つ立方晶セルにおいて、帰還並進周期が倍増し、$k_x = \pi/a$ および $k_y = \pi/a$ において追加のゼロ周波数接点が生じる。これらは回路計算とフルウェーブシミュレーションの両方によって確認されている。
• 非連結の交差ロッド： 2本の非連結な直交ロッドについて、モデルは1つの光子的な線形バンドと1つの遮蔽された$\Gamma$二次的バンドを予測する。後者は、ゾーンコーナー付近で双曲型の等周波数等高線（$|k_x k_y| = \text{const}$）を生じさせる。
• 3本の非連結直交ロッド： モデルは、3本の非連結直交ロッド（$N=3$, ランク $W=3$）を持つ構造が、2つの遮蔽された$\Gamma$二次的バンドを持ち、プラズマカットオフ以下に一様電場由来のバンドを持たないことを予測する。フルウェーブシミュレーションは、$\Gamma$で傾きが消失する2つの低周波バンドの存在を支持している。
• 検証： 回路計算と選択されたフルウェーブシミュレーション（COMSOL）は、予測された接点集合および分散次数と一致している。モデルは低周波バンドを正確に再現し、接続性に基づいて線形挙動と二次的挙動を区別している。

意義
本論文は、金属の配線（wiring）が、準静的なフォトニックバンドのための離散的な設計変数を提供することを主張している。設計を離散的な制約（配線／接続性）と連続的な係数（幾何学／材料パラメータ）に分離することで、この理論は、導体のトポロジカルな接続のみに基づいてバンド構造を予測することを可能にする。このアプローチは、対称性と結合をチューニングすることでフォトニック・ハミルトニアンを合成する既存の戦略を補完するものである。著者らは、結果として得られる接点シートや遮蔽バンドは、共鳴による平坦化ではなく、電流の連続性と静電遮蔽から生じるものであることを強調している。この理論は、連続的な最適化の前に、まず離散的な配線設計から始めるという、金属ネットワークにおける準静的バンドエンジニアリングの出発点を提供する。配線制約と動的係数の分離は中心的な結果であり、特定の幾何学的対称性を必要とせずに、任意の相互接続パターンに対するモード数と分散次数を予測することを可能にしている。