体ディラック半金属-絶縁体-金属構造を用いた超コンパクトテラヘルツプラズモニック導波路：閉じ込め強化と動的チューニング

概要

バルクディラック半金属（BDS）-絶縁体-金属（BIM）スタック上に構築されたサブ波長テラヘルツプラズモニック導波路を理論的に検討した。この構造は、表面プラズモンポラリトン（SPP）モードが強く閉じ込められる（閉じ込め因子最大½·λ₀/15）最適化周波数帯を提供し、伝播損失を低く（1.0 dB/λ₀）抑制する。BIM導波路に2本のシリコンリボンを導入することで、動的にチューナブルなファブリ-ペロー共振器が形成され、ディープサブ波長スケールでテラヘルツSPPを制御可能、超コンパクト・チューナブルTHzプラズモニックデバイスと光フィルタへの道を開く。

背景

テラヘルツ（THz）放射は、イメージング、生化学センシング、無線通信で多用途ツールとして台頭し、より強い閉じ込め、高分解能、高集積化が求められている（文献1–3）。THz波のサブ波長閉じ込め達成は極めて重要だが、貴金属における従来の表面プラズモンポラリトン（SPP）は、弱い閉じ込め、高い固有損失、受動的チューナビリティのみという限界で実用応用を妨げてきた（文献4–6）。

グラフェンプラズモンは低損失、動的チューナビリティ、極限的閉じ込めで有望視され、ゲートチューナブルグラフェンヘテロ構造、離散タルボット効果、プラズモン誘起透明性導波路など高分解能THzデバイスを実現（文献11–18）。これら進展は、ディープサブ波長スケールでのTHz放射操作における極限的SPP閉じ込めの重要性を強調する。

バルクディラック半金属（BDS）は「3Dグラフェン」と称され、最高グラフェン（2×10⁵ cm² V⁻¹ s⁻¹）を大きく超える超高キャリア移動度（最大9×10⁶ cm² V⁻¹ s⁻¹）を示し、低損失プラズモンの鍵となる（文献19）。誘電応答もフェルミエネルギー変調でチューナブル、Na₃Bi、Cd₃As₂、AlCuFe準結晶などはグラフェンより加工容易で安定、BDSを次世代プラズモニック材料として位置づける（文献20,21）。しかし、BDS-絶縁体界面単独でのSPPモード閉じ込めは、ディープサブ波長操作にはまだ不十分である。

本研究では、SPP閉じ込めを大幅に改善しつつ低損失と動的チューナビリティを維持するBIM導波路を提示。この高閉じ込めモードの分散、損失特性、フィルタリング能力を体系的に調査し、強化閉じ込めと伝播損失低減を両立する最適化周波数ウィンドウを明らかにした—これは従来金属ベースSPPでは稀な挙動。さらに、BIMモードはλ₀/2000未満の狭スリットを伝播可能、シリコンリボンをミラーとして埋め込むことで、BDSフェルミエネルギー調整で共振周波数シフト可能なチューナブルバンドパスフィルタを実現。

理論とシミュレーション

BIM導波路は、単層BDS膜（厚さ0.2 µm）が誘電体スペーサー（誘電率εᵣ）を介して銀基板からギャップgで分離された構造。THz領域では銀は完全電気伝導体として振る舞う。TM偏光励起で、SPPモードがx方向に波数k_SPPで伝播し、y方向に自由空間へ減衰。境界条件適用で分散関係導出：

$$-\frac{\varepsilon_r\sqrt{k_{\mathrm{SPP}}^2-k_0^2}}{\varepsilon_0\sqrt{k_{\mathrm{SPP}}^2-\frac{\varepsilon_r k_0^2}{\varepsilon_0}}}=\left(1+\frac{i\sigma \sqrt{k_{\mathrm{SPP}}^2-k_0^2}}{{\omega \varepsilon}_0}\right)\tanh\left(g\sqrt{k_{\mathrm{SPP}}^2-\frac{\varepsilon_r k_0^2}{\varepsilon_0}}\right),$$

ここでk₀は自由空間波数、σはBDS膜の複素伝導度（式3–4）。式(1)を解き実効屈折率 n_eff = k_SPP/k₀ = Re(n_eff) + i Im(n_eff)を得る。Re(n_eff)がモード閉じ込め度、Im(n_eff)が伝播損失に直結。

十分大きなギャップ（tanh[…]≈1）では式(1)が単層BDS分散関係（式2）に簡約、BIM構成での金属基板の追加利益を強調。

BIMプラズモニック導波路の模式図：単層BDS膜がギャップgで銀基板から分離、誘電率εᵣの誘電体スペーサー介在。TM偏光SPPモードがx方向伝播、y方向減衰。Eₓ場分布を赤線で表示。

結果と考察

まずBIMギャップgとBDSフェルミエネルギーE_Fがモード閉じ込めと損失に与える影響を検証。E_F=70 meVで、様々なgでの実効屈折率を図2a（Re(n_eff)）と図2b（Im(n_eff)）にプロット。g≥10 µmで曲線が0.05 THz以上で収束、SPP場が~10 µm以内に閉じ込められ銀寄与が無視可能。gを1 µm未満に縮小すると閉じ込め顕著増大、損失減少、閉じ込めピークと損失減少が共存する最適化周波数帯が出現—従来金属ベースSPPでは稀な特徴。

図2c,dでE_F依存性（g=1 µm）を示す。E_F増加で閉じ込めと損失とも低下、金属性増強とキャリア緩和時間延長に起因。例：2.5 THzでg=10 nm、E_F=70 meVにてλ₀/15閉じ込めと1.0 dB/λ₀損失達成。

Re(n_eff)とIm(n_eff)の(g)依存性（E_F=70 meV）と(E_F)依存性（g=1 µm）

透過シミュレーション（図3）で損失低減確認：1.56 THzでBIM導波路透過率0.97、単層BDSを上回る。対応する磁場分布（図3b,c）でBIMのRe(n_eff)=2.45、単層1.002、場がλ₀/2000未満スリット内に強く局在化。

E_F=70 meV、g=50 µm、1.56 THzでのBIM（赤）と単層（青）導波路の透過スペクトル（a）と磁場分布（b,c）

高閉じ込めモード活用のため、誘電体スペーサーに2本シリコンリボン（n_Si=3.4）を統合、反射ミラー形成（図4a）。結果ファブリ-ペロー共振器が定常波共振支持、1.56と2.22 THzに2透過ピーク、FWHM 0.12と0.09 THz。磁場マップ（図4c,e）で第1、第2次共振モード確認。共振外では定常波形成不能、入力反射（図4d）。

透過は結合モード理論（CMT）で解析的記述可能：

$$T(\omega)=\frac{\kappa_w^2}{(\omega-\omega_0)^2-(\kappa_w+\kappa_i)^2},$$

κ_wとκ_iは導波路結合と固有損失率。品質因子Q_t=ω₀/FWHMとQ_oi=−Re(n_eff)/(2Im(n_eff))でκ_w、κ_i抽出、解析結果が数値スペクトルと一致（図4b）。

(a) シリコンリボン（幅d、間隔L）入りBIM導波路. (b) g=1 µm、d=5 µm、L=120 µmでの数値（青）とCMTフィット（赤）透過スペクトル. (c–e) 1.56、1.90、2.22 THzでの磁場分布

図5で共振周波数の共振器長L依存性を探索。第1、第2共振ともL増加で赤方シフト、定常波条件2k_SPP(ω_r)L+θ=2mπと整合。同様にg増加で青方シフト（図6a）、Re(n_eff)減少起因。E_F増加でも青方シフト（図6b）、帯域幅狭小化、Im(n_eff)低下反映。これらチューナビリティ機構で構造再設計なしにフィルタ動的制御可能。

(a) 変化するLでの透過スペクトル. (b) 第1、第2モード共振周波数のL依存性（g=1 µm、d=5 µm、E_F=70 meV）

透過スペクトルの(g)依存性(a)と(E_F)依存性(b)（他パラメータ固定）。(c) Re(n_eff) vs E_Fとg。 (d) Im(n_eff) vs E_Fとg。

結論

BIM導波路はλ₀/15閉じ込めと1.0 dB/λ₀損失の超コンパクトTHzプラズモニックモードを実現、λ₀/2000未満スリット伝播可能。シリコンリボン導入で動的チューナブルバンドパスフィルタを構築、共振周波数がBDSフェルミエネルギーでシフト、構造再最適化なしでTHzスイッチングとフィルタリングを可能にする。

実験方法

数値シミュレーションに2D有限差分時間領域（FDTD）ソルバー使用、xとy方向に完全整合層。BDSメッシュ1 µm × 0.02 µmで収束確保。BDS周波数依存伝導度はランダム位相近似下のKubo式従う—式(3)(4)でReσ(Ω)とImσ(Ω)記述。AlCuFeパラメータ：t=40、ε_c=3、μ=3×10⁴ cm² V⁻¹ s⁻¹、E_F=70 meV。

本研究にヒト参加者、データ、動物は関与しない。

略語

BDS

バルクディラック半金属

BIM

BDS-絶縁体-金属

CMT

結合モード理論

FDTD

有限差分時間領域法

FWHM

半値全幅

SPPs

表面プラズモンポラリトン