結合共鳴グラフェンメタマテリアル吸収体：中赤外領域での調整可能な変調

概要

中赤外領域におけるグラフェン搭載メタマテリアル吸収体を研究しました。クロス形状のスロットによる結合共鳴を通じて、光とグラフェンの相互作用が大幅に強化されています。フェルミレベルの上昇に伴い吸収ピークが著しくブルーシフトし、吸収体に幅広い調整範囲をもたらします。シンプルな回路モデルがこの変調挙動を適切に説明・予測します。本提案は、スイッチング、センシング、変調、生化学検出など多様な分野での応用が期待されます。

背景

プラズモニックメタマテリアル（PM）吸収体は、深いサブ波長スケールの金属ナノ構造を用いて動作します。特定波長で完全吸収を実現・制御でき、発光素子/検出器、センサー、光熱療法、光機械相互作用、ハイパースペクトルイメージングなど幅広い用途につながります[1-7]。PM吸収体は、液晶、半導体、相変化材料などを導入することで電気的・光学的・熱的に光学応答を変調できる、調整可能な機能デバイス設計の有望なプラットフォームでもあります[8-13]。これにより新型の変調器、スイッチ、マルチスペクトル検出器が可能になります。

近年、グラフェンはプラズモニック材料としての高速変調能力と調整可能性から大きな注目を集めています[14-20]。グラフェンの伝導率はフェルミレベル（E_F）に依存し、バイアス電圧でナノ秒オーダーで連続調整可能なため、近赤外・中赤外領域で高い変調率を実現できます[17,19-24]。しかし、単層グラフェンは原子レベルの厚さしかないため、入射光とプラズモン共鳴の相互作用は極めて弱く、中赤外領域ではパウリブロッキングによりさらに弱まります[22]。その結果、波長調整範囲と変調深度は限定的で、共鳴波長の10%未満にとどまります[21,22,25-28]。これは光通信や広帯域スペクトル検出への実用応用において課題となっています。したがって、効率的な電気光学変調にはグラフェン-光相互作用の大幅な強化が必要です。ナノアンテナやスプリットリング共振器などの複雑なナノ構造設計により、グラフェン-光相互作用の強化が理論・実験的に示されてきました[19,21,22,25,27,28]が、これらは構造が複雑、偏光依存、動作周波数範囲が狭い、調整範囲が限定的といった問題がありました。

本研究では、9〜14 μmの変調範囲を持つグラフェン搭載吸収体を提案します。これは生化学センシングやサーモグラフィなどの応用に適しています[5,29-31]。クロス形状スロット内の結合共鳴により電界が4桁強化され、グラフェン-光相互作用が強く増大し、中心波長で最大25%のシフトを実現しました。また、電圧と幾何学的パラメータで制御されるグラフェン誘起変調を説明・予測するシンプルなLC回路モデルを提案します。このような広い調整範囲は多くの応用で有望です。

方法

図1aに示すように、周期Λ=8 μmでパターン化された金属パッチが、誘電体スペーサーを挟んで金属基板上に配置されています。単層グラフェンがパッチとスペーサーの間に挟まれています。基板は十分厚く反射鏡として機能します。スペーサー層の厚さはt_d=520 nm、金属パッチの厚さはt_m=100 nmです。図1bは1単位セルの上面図です。偏光独立性を実現するため、2つのサブユニットが対角対称に配置されています。各正方形パッチにクロス形状のスロットをエッチングし、4つの同一の小片に分割しています。S₁とS₂の小片サイズはそれぞれl₁=1.5 μm、l₂=1.7 μm、両サブユニットのスロット幅はa=20 nmです。金属材料には金（Au）を用い、その光学特性はプラズマ周波数ω_p=1.369×10¹⁶ Hz、衝突周波数τ=1.224×10¹⁴ Hzのドゥルードモデルで記述します[32]。誘電体スペーサーは硫化亜鉛（ZnS）で、中赤外領域では屈折率n=2.2、損失は無視できます[33]。

a 提案するグラフェン搭載メタマテリアルの概略図。各サブユニットのクロス形状スロットにより、偏光依存性なしでグラフェン-光相互作用を大幅に強化。b 1周期内の構造上面図。異なるパッチサイズの2サブユニットが対角配置。

有限差分時間領域法（FDTD；Lumerical FDTD Solutions）を用いて反射率スペクトルと電磁界分布を計算しました。シミュレーションはx、y方向に周期境界条件、z方向に完全整合層条件で実施。単層グラフェンは表面伝導率アプローチにより二次元構造としてモデル化[34]。グラフェン層の表面伝導率σ_g（帯間項σ_interと帯内項σ_intraを含む）はクーボ公式で計算[35]。

式(1)および(2)より、グラフェンの光学定数はE_Fに応じて変化し、この変化が吸収周波数の調整可能性をもたらします。ナノ構造の結合共鳴によりこの調整範囲が大幅に拡大し、デバイスに必要な印加電圧を大幅に低減できます。

結果と考察

図2aは、正入射時のx偏光波（φ=0）における吸収スペクトルを示します。フェルミレベルE_F=0 eVでは、波長λ=12.4 μmと13.3 μmに2つの吸収ピークが観測され、12.1〜13.5 μmの入射光がほぼ完全に吸収されます。E_Fが増加すると共鳴は短波長側へシフトします。E_F=0.2 eVではピークが11.8 μmと12.46 μmに移動し、それぞれ4.8%と6%の相対シフトです。一方、ピーク2の吸収率は高E_Fで低下しますが、これはメタマテリアルと空気のインピーダンスミスマッチに起因します[28]。興味深いことに、E_Fがさらに増加するとピーク2の方がピーク1より速くブルーシフトします。この挙動は後で回路モデルで説明します。

正入射時の吸収スペクトル。a 異なるE_Fでのスペクトル（φ=0）、E_F増加に伴うピークの大きなブルーシフト。b 異なるφでのスペクトル（E_F=0.2 eV）、偏光独立性を示す。偏光角φは図1aで定義。

変調度をパラメータM=Δλ/λ₀で定量化します。λ₀はE_F=0 eVでの共鳴波長、ΔλはE_F変化による波長シフトです。図2aより、E_F=0.6 eVでピーク1のM₁=20.1%、ピーク2のM₂=25.5%となり、先行研究[19,21,22,25-28]よりはるかに広い変調範囲です。このような低E_Fでの大変調は多くの応用で極めて望ましい特性です。別計算により、スペーサー厚を薄くすると吸収ピークがブルーシフトすることも確認しました（補足資料1）。したがって、厚みを最適化して変調の出発点を設定可能です。また、図2bに示すように提案メタマテリアルの光学応答は偏光独立で、偏光角φを0〜90°まで変化させても吸収スペクトルは不変です。これは設計の対称性によるものです。

完全吸収のメカニズムは共鳴時の電界分布で明確に説明できます。図1に示すよく知られた金属-絶縁体-金属（MIM）構造[3,32,36-38]により、局在表面プラズモンポラリトン（SPP）が励起され、各パッチにコンパクトな磁気共鳴が形成されます。図3a、bはE_F=0.2 eVでの共鳴波長λ₁=11.8 μm、λ₂=12.46 μmにおけるグラフェン層の規格化磁界|H|²分布です。SPPが強く局在しているため、2サブユニットは独立して動作します。しかし、各サブユニット内部の狭い分割スロットにより、4つの小片の共鳴は互いに結合しています。この結合によりスロット内の電界が著しく増大します（図3c、d）。ここでは入射光がx偏光であるため、y方向スロットのE界のみが顕著です。共鳴結合により強化されたE界強度は入射光E_incの4桁以上大きくなります。対照的に、先行研究で変調に用いられた最も強化された界はパッチ端部にあります。図3e、fは図3c、dの白線に沿ったスロットと端部の強化比較を示します。

E_F=0.2 eV、正入射x偏光でのグラフェン層電界分布。結合共鳴によるスロット内の巨大な強化を示す。a,b 規格化磁界|H|²（λ₁=11.8 μm、λ₂=12.46 μm）。c,d 対応する|E/E_inc|²分布。e,f 白破線に沿った|E/E_inc|²。スロット内とパッチ端部の強度の鮮明な対比が、先行研究より遥かに広い調整範囲を示唆。

このような電界分布は、本提案で変調が非常に大きい理由をよく説明しています。摂動理論に基づき、グラフェン誘起共鳴シフトはΔω=−iσ_g∫_S|E_s|²dS/W₀で評価されます[22]。ここで|E_s|²はグラフェン層の電界強度、W₀は蓄積エネルギー、Sはグラフェンで覆われた面積です。共鳴のスペクトルシフト（Re(Δω)）はσ_gの虚部で決まり、中赤外領域では実部よりはるかに大きいです[22,28]。図3c-fに明確に示されるように、狭いスロット内の電界強化は端部の10倍以上です。したがって、積分値は主にスロット内の巨大なE界強化に由来し、金属端部でのE界強化のみを持つ先行ケース[21,22,25,27,28]より遥かに大きなピークシフトをもたらします。

上記の電界分布と考察に基づき、調整挙動を研究するLC回路モデルを提案します。図4aに示すように、L_iとC_i（i=1,2）は図1bのパッチS_iのインダクタンスとキャパシタンスです。スロット幅aが十分大きくグラフェン層がない場合、スロットとグラフェンの効果は無視できます。このときL_iとC_iは、吸収スペクトルから得られる共鳴波長とのフィッティングにより個別に決定できます[37,39,40]。結果はS₁でL₁=0.07 pH、C₁=350 aF、S₂でL₂=0.075 pH、C₂=380 aFです。各サブユニット内のスロット誘起結合効果はシャントキャパシタンスC_cで記述され、スロット幅aの増加とともに減少します。本ケースではa=20 nmでC_c=290 aF、以降aが10 nm増えるごとに200 aF、180 aF、135 aFになります。共鳴波長は回路インピーダンスをゼロとすることで得られます：λ_i⁰=2πc₀√(L_iC_i⁰)。ここでc₀は真空中の光速、「i」はサブユニットS_iを指し、C_i⁰=C_i+C_cです。

a 個別パッチ（L_i、C_i）、スロット（C_c）、グラフェン（L_g）の寄与を含むLC回路モデル。b LCモデル計算とFDTDシミュレーションの共鳴比較。c,d E_F=0.4 eVでの単一パッチ共鳴シフト。幾何学パラメータc スロット幅（l=1.5 μm）、d パッチサイズ（a=20 nm）変更時。

二次元グラフェン層は基本的にインダクタとして振る舞います。図3より、グラフェン層の主寄与は電界が強化されるスロット位置から来ます。スロット幅は動作波長およびグラフェンプラズモン波長よりはるかに小さいため、準静近似が有効です。スロット間の電圧Vと電流IはV=aE、I=2l_it_g(σ_g−iωε₀)Eで評価できます（Eはグラフェン層の電界）。これによりインダクタンスL_g=−1/ω Im(V/I)[41]を導入し、グラフェン層の寄与を記述すると、

L_g=a/(2l_iω²ε₀|Re(ε_g)|t_g) （i=1,2） ……（3）

このインダクタは図4aに示す並列素子として機能します。結果として、1パッチの全インダクタンスは1/L'_i=1/L_i+1/L_gで得られ、グラフェン層ありの各サブユニットの最終共鳴波長は

λ'_i=2πc₀√(L'_iC_i⁰) （i=1,2） ……（4）

各サブユニットは独立に動作するため、メタマテリアルの全インピーダンスは2サブユニットのインピーダンス並列接続から得られます。

このLCモデルは、E_F増加に伴う共鳴のブルーシフトを予測します。式(1)、(2)より、高E_Fではグラフェンの|Re(ε_g)|が大きくなり、式(3)でL_gが小さくなります。インダクタが並列接続されているため、全インダクタンスL'_iが小さくなり、式(4)で共鳴波長が短くなります。計算結果を図4bにまとめ、FDTDシミュレーションで得られた共鳴波長と良好な一致を示します。わずかな偏差は、LCモデルが各パッチ端部の弱い界の寄与を無視しているためです（図3c-f）。LCモデルは幾何学パラメータがブルーシフトにどう影響するかも示します。式(4)を微分すると∂λ'_i/∂L'_i∝1/√L'_iとなり、√L'_iの小さい値がこのブルーシフトの感度を高めることがわかります。インダクタは並列接続でL_iは固定されているため、全インダクタンスL'_iの小さい値はグラフェンインダクタンスL_gの小さい値を意味します。式(3)より、調整範囲を広げるにはスロット幅aを小さく、パッチサイズlを大きくすべきです。図4cは、S₁内のスロット幅を50→20 nmに減らすと、E_F=0.4 eVでの共鳴ブルーシフトが約6%→15%に増加することを示します。一方、スロット幅をa=20 nmで固定し、パッチサイズを1.5→1.8 μmに変えると、共鳴シフトが15%→22%に増加します（図4d）。FDTDシミュレーションとの良好な一致は、このシンプルな回路モデルが関連メタマテリアルデバイスの研究に効率的な手法であることを実証しています。

結論

結論として、偏光独立、広帯域、大変調範囲のメタマテリアル吸収体を設計しました。両共鳴とも、E_Fを0→0.6 eVに増加させると中心波長の20.1%と25.5%まで調整可能です。この大変調は、各金属パッチのクロス形状スロット内の結合共鳴により、グラフェン-光相互作用が著しく強化されることに由来します。この効果は、LCモデルにおけるグラフェン導入インダクタで適切に記述されます。このシンプルなモデルは、異なる幾何学パラメータでの変調挙動を予測し、結果はFDTDシミュレーションとよく一致します。本提案は光通信、センシング、サーモグラフィなどの潜在的応用に有益です。

略語

E_F：

フェルミレベル

FDTD：

有限差分時間領域法

MIM：

金属-絶縁体-金属

PM：

プラズモニックメタマテリアル

ZnS：

硫化亜鉛