楕円グラフェン-黒リンナノ構造による高度可変デュアルバンド赤外線メタマテリアル吸収体

概要

周期的な楕円グラフェン-黒リン（BP）ペアに基づく赤外線領域デュアルバンド吸収体を数値提案。グラフェンとBPの組み合わせにより、両共鳴でニアユニティ異方性吸収を実現。各共鳴は幾何学パラメータ調整で独立にチューニング可能。さらに、グラフェンとBPのドーピングレベルでも共鳴特性を効果的に制御可能。電場分布解析でグラフェン-BP楕円内の表面プラズモン共鳴を観測、強力で異方性なプラズモニック応答に寄与。また、入射角ロバスト性と偏光感度も解説。

はじめに

グラフェンはハニカム格子に配列した炭素原子からなる二次元材料[1, 2]。超コンパクトサイズと独自の光-グラフェン相互作用により、近年様々なグラフェンベースフォトニックデバイスが開発[3-6]。最も重要な応用の一つとして、グラフェンベースメタマテリアル吸収体が強力でチューナブルなプラズモニック応答により爆発的関心を集めています[7-10]。しかし、グラフェンのゼロまたは近ゼロバンドギャップにより、高オン/オフ比を必要とする用途が制限されます[11]。代替二次元材料として、黒リン（BP）はプッカード構造を持つ六方最密充填格子に配列したリン原子単層[12]として、近年研究関心が急増。面内異方性、厚み依存的チューナブルバンドギャップ[13]、高キャリア密度と移動度[14]など卓越した光学的・電子的特性を持ちます。過去数年、赤外線領域でBPベースメタマテリアルにおける光-BP相互作用強度向上のため多くの構造が研究されてきました[15-17]。しかし、BPベース吸収体のプラズモン共鳴は柔軟かつ効果的なチューニングが困難で、中程度ドーピングでは吸収率が比較的低い傾向。これは単層BPの共鳴強度が弱く、異方性ポテンシャルを制限するため。グラフェン-BPベースプラズモニック吸収体が、グラフェンとBPのハイブリダイゼーションを利用し強力で異方性なプラズモニック吸収達成のため提案されてきました[18-20]。しかし、従来報告されたグラフェン-BP吸収体は概して複雑な製造技術を必要とし、または単一吸収帯のみで、イメージング、バイオセンシング、通信システムへの更なる応用を阻害。

本研究では、周期的楕円グラフェン-BPペアを用いた異方性デュアルバンド赤外線吸収体を数値提案。製造容易性を備え、幾何学サイズとドーピングレベルによる共鳴の独立チューナビリティを実証。電場分布プロットで物理メカニズムを解明。入射角許容度と偏光感度も解説。

方法

提案吸収体は、図1に示す通りSiO₂層上に堆積した横方向と縦方向の楕円グラフェン-BPペアで構成。単層グラフェンとBPの間に六方窒化ホウ素（hBN）層を絶縁スペーサーとして挿入し、キャリア輸送を防止し高キャリア移動度を保証。SiO₂とhBNパラメータは文献[21]と[22]から取得。シミュレーションはCOMSOL Multiphysicsで実施、周波数領域有限要素法（FEM）ベース。x-とy-方向両方にフロッケ周期境界条件適用。計算領域上部表面に赤外線波励起ポート、底面に完全電気伝導体（PEC）境界条件設定。ユーザー制御メッシュ密度の四面体メッシュを全領域に適用。

シミュレーションではグラフェンとBPを、体積材料の誘電率テンソルではなく表面伝導率を持つ二次元面として扱う。この仮定により超薄膜材料の厚さ定義問題と低計算効率問題を解決[23]。

グラフェン表面伝導率σ(ω)記述に、よく知られたKubo公式を使用[24]：

σ(ω, μ_c, Γ, T) = σ_intra + σ_inter （1）

σ_intra = (je²/πħ²(ω - j2Γ)) × ∫₀^∞ ξ (∂f_d(ξ, μ_c, T)/∂ξ - ∂f_d(-ξ, μ_c, T)/∂ξ) dξ （2）

σ_inter = (-je²(ω - j2Γ)/πħ²) × ∫₀^∞ (f_d(-ξ, μ_c, T) - f_d(ξ, μ_c, T))/((ω - j2Γ)² - 4(ξ/ħ)²) dξ （3）

f_d(ξ, μ_c, T) = (e^((ξ-μ_c)/k_BT) + 1)⁻¹ （4）

式1によりσ(ω)は帯内と帯間成分σ_intraとσ_interで構成。ωは角周波数、μ_cは化学ポテンシャル、Γは散乱率、Tはケルビン温度。ħ、e、ξ、k_Bは換算プランク定数、電子電荷、電子エネルギー、ボルツマン定数。

赤外線領域では入射光子が帯間遷移を励起しにくいため、光-グラフェン相互作用は帯内遷移が支配。特にμ_c ≫ k_BT時、Kubo公式は式5に簡略化：

σ_g = i e² μ_c / (π ħ² (ω + i2Γ)) （5）

したがってグラフェン表面伝導率はω、Γ、μ_cに依存。ここではΓ=0.3 meV、μ_c=0.7 eVと仮定（先行研究[25, 26]による）。

一方、BP表面伝導率σ_jをシンプルな半古典的Drudeモデル[27]で計算：

σ_j = iD / (π (ω + iΓ_BP/ħ)) （6）

D_j = π e² n_s / m_j （7）

n_sはドーピングレベルに関連するキャリア密度。n_s = 1.9×10¹³ cm⁻²、Γ_BP = 10 meVを先行参照[16]から採用。jは関心方向で、σ_xとσ_yはx-とy-方向電子質量で決定。m_xとm_yはさらに：

m_x = ħ² / (2γ²/Δ + η_c) （8）

m_y = ħ² / (2ν_c) （9）

η_c = ħ² / (0.4 m₀) （10）

ν_c = ħ² / (1.4 m₀) （11）

γ = 4a/π （12）

m₀は標準電子質量、ΔとaはBP単層のバンドギャップとスケール長。式10-12を式8、9に代入し、アームチェア（x-）とジグザグ（y-）方向電子質量を取得。それらの差異がBPの異方性表面伝導率に寄与。

結果と考察

提案吸収体の異方性吸収特性を説明するため、単独グラフェン層、単独BP層、グラフェン-BPペアの吸収スペクトルをシミュレーション比較。図2aに示す通り、グラフェンプラズモニック応答は等方的で9.9 µmと15.4 µmに明確な2吸収ピーク、偏光非依存。一方BPプラズモン共鳴は異方的だが、TE（<12.7%）またはTM（<0.7%）入射いずれも強度がかなり弱い。グラフェンとBPの利点を組み合わせ、グラフェン-BPペアは強力かつ異方性プラズモニック応答を示す。TE入射で2吸収ピークが8.8 µmと14.1 µmに位置、吸収率90%以上。TM入射では最大吸収波長が9.5 µmと15.4 µmにシフト。偏光消光比PER = 10×log(R₁/R₀)で定義（R₁、R₀は同波長異偏光反射率、R=1-A、Aは吸収率）。各共鳴の最大PERはλ=9.5 µmで23 dB、λ=14.1 µmで25 dBに達する。したがって提案吸収体は高性能デュアルバンド反射偏光子として利用可能。