可視光域で95%効率の超狭帯域完全吸収体によるプラズモニックセンシング

概要

吸収帯域幅わずか1.82 nm、可視域で吸収効率95%を超える超狭帯域完全吸収体を提案し、数値的に検証しました。この完全超狭帯域吸収は、局在表面プラズモン共鳴（LSPR）によって誘起される結合効果に起因することを実証。構造寸法が光学性能に及ぼす影響も調査し、共鳴ディップの反射率が極めて低い（0.001）最適構造を得ました。この完全吸収体は屈折率センサーとして動作可能で、感度は約425 nm/RIU、フィギュア・オブ・メリット（FOM）は233.5に達し、可視域プラズモニックセンサーの精度を大幅に向上させます。さらに、固定周波数での強度変化検出性能を表すFOM*も計算し、最大1.4 × 10⁵を達成。この高いセンシング性能により、本メタマテリアル構造は生物結合、集積フォトディテクタ、化学応用などに大きな潜在能力を持ちます。

背景

近年、局在表面プラズモン共鳴（LSPR）ベースのプラズモニックメタマテリアルは、その電磁特性と単極共振器[1,2,3,4,5,6,7,8]、光透過増強[9,10,11,12,13]、プラズモニックセンサー[14,15,16,17,18,19,20,21]などの有望な応用により大きく進展しています。メタマテリアル吸収体では、電磁波吸収の増強が有利ですが、他デバイス設計では金属の固有光学損失を慎重に考慮する必要があります。最初の完全メタマテリアル吸収体はLandyによって提案・実証されました[22]。以降、完全メタマテリアル吸収体は急速に発展し[23,24,25,26,27,28,29,30,31]、応用要件に応じて狭帯域吸収体と広帯域吸収体に大別されます。一般に広帯域吸収体はエネルギー収集、狭帯域吸収体はセンサーや単色フォトディテクタに用いられます。

センシング応用では、狭帯域吸収体ベースのプラズモニック屈折率センサーが注目されています。これまでに、赤外線・テラヘルツ領域で動作する様々なタイプのプラズモニック屈折率センサーが報告されています：ハイブリッドマイクロキャビティ[32, 33]、ナノディスク[34]、ネットワーク型メタサーフェス[24]、金属格子[28]など[35,36,37]。可視帯域で動作するセンサーは、赤外線・テラヘルツ・マイクロ波領域のものと比較して、メタマテリアル構造のより小さな周期性を実現でき、化学・生物学など多くの実用応用でのデバイス利用効率を向上できます[38]。しかし、これまで報告された可視域プラズモニック屈折率センサーは概して比較的低いFOMであり、更なる開発・応用を大きく妨げています。理論研究では、2015年にZhouらがサブ波長金属格子構造を用いた可視域屈折率センサーを理論提案し、感度300 nm/RIUを達成しましたがFOMはわずか2でした[28]。Liuらは深サブ波長プラズモニックナノキャビティを用いたマルチスペクトルセンサーを設計し、FOM 58を示しました[34]。Liuらの努力により、ネットワーク型メタサーフェスを持つプラズモニック構造で最小FWHM 3 nm、FOM 68.57の屈折率センサーが得られました[24]。実験研究では、2014年にEmikoとTetsuが単一AuナノスターベースのLSPRセンサーを実験実証し、感度665 nm/RIU、FWHM最大40 nmを達成[39]。Choらは感度378 nm/RIUに達するプラズモニック屈折率センサーを実験実証[40]。理論・実験両面で多くの研究者が可視域屈折率センサーのFOM向上に努めてきましたが、高FOM設計は依然として大きな課題であり、応用を深刻に制限しています。

センサーではFOM向上が非常に意義深いです。例えば生物学分野では、より高いFOMは分子検出におけるより強い性能を意味します。本研究のセンサーFOMは233.5に達し、可視域で報告されたプラズモニック屈折率センサー[24, 28, 34]よりはるかに高い値です。本プラズモニックセンサーは金属-誘電体-金属（MDM）周期構造に基づき、同時に可視域で吸収効率95%以上、FWHMわずか1.82 nmの完全超狭帯域プラズモニック吸収体としても動作します。構造寸法と材料パラメータがメタマテリアルの光学特性に及ぼす影響も調査。さらに、一般的なMDM構造と比較して、三角ナノリボンの使用が吸収性能向上に寄与することを実証。同時に吸収メカニズムも詳細に解析。提案構造の作製を考慮すると、三角ナノリボンは電子ビームリソグラフィ[41]、モールディング[42]、インプリントリソグラフィ[43]など多くの方法で製造可能。本成果がプラズモニックセンサー設計の指針となることが期待されます。

手法

図1に提案メタマテリアル構造の1単位セル断面を示します。構造は、誘電体層と基板に挟まれた薄い金層上の2つの金ナノリボンアレイと、それらの間の三角金ナノリボンで構成されます。シミュレーションでは、金の誘電率をドゥルードモデルで特性評価。中間誘電体層と基板の誘電体をそれぞれNaF（n = 1.3）とMgF₂（n = 1.4）に設定。二次元有限差分時間領域（FDTD）法で提案構造の透過と反射を計算し、構造全体の吸収をA = 1 − R − Tで定義。x方向に周期境界条件を設定し、x方向偏光の横磁界（TM）波を構造に垂直入射させました。

提案メタマテリアル構造の1単位セル模式図

周知の通り、等価LC回路モデルは完全吸収体におけるLSPR励起磁気共鳴を定性的に予測するため広く用いられています[44,45,46]。LCモデル議論の便宜上、メタマテリアル吸収体構造の模式図を図2a、等価LCモデルを図2bに示します。ここでは、隣接単位のナノリボン間ギャップ容量をC_g = ε₀t₁/(P − d − 2w)で表し、ε₀は周囲環境の誘電率。ナノリボンと金膜間の容量をC_m = c₁ε₃ε₀(2w + d)/t₃で表し、c₁は金属表面の非均一電荷分布による係数、ε₃は誘電体層の誘電率[44,45,46]。金ナノリボンと金膜の相互インダクタンスはL_m = 0.5μ₀(2w + d)t₃で与えられ、μ₀は周囲環境の透磁率。金ナノリボンと金膜の漂遊電荷寄与を考慮し、運動インダクタンスを\( {L}_e=\left(2 w+ d\right)/\left(\gamma {\varepsilon}_0{t}_1{\omega}_p^2\right) \)で与え、γは金ナノリボンの実効断面積を考慮する係数、ω_pは金のプラズマ周波数[44,45,46]。すると、等価LC回路モデルの総インピーダンスは以下で表せます：

a メタマテリアル吸収体構造の模式図。b 図6a構造の等価LC回路模式図

$$ {Z}_{\mathrm{t} ot}=\frac{i\omega \left({L}_m+{L}_e\right)}{1-{\omega}^2{C}_g\left({L}_m+{L}_e\right)}-\frac{2 i}{\omega {C}_m}+ i\omega \left({L}_m+{L}_e\right) $$ (1)

共鳴波長はZ_tot = 0のときに得られます。

$$ {\lambda}_r=2\pi {c}_0{\left(\frac{C_m{C}_g\left({L}_m+{L}_e\right)}{C_m+{C}_g-\sqrt{C_m^2+{C}_g^2}}\right)}^{\frac{1}{2}} $$ (2)

隣接単位のナノリボン間結合は、大きなギャップ（P − d − 2w）のため非常に弱いです。C_gがC_mの5%未満のとき、C_gの影響は無視できます。この状況では、共鳴波長は以下のように簡略化されます：

$$ {\lambda}_r\approx 2\pi {c}_0\sqrt{\left({L}_m+{L}_e\right){C}_m} $$ (3)

ここで、L_m = 0.5μ₀(2w + d)t₃、\( {L}_e=\left(2 w+ d\right)/\left(\gamma {\varepsilon}_0{t}_1{\omega}_p^2\right) \)、C_m = c₁ε₃ε₀(2w + d)/t₃。LC回路モデルでは、式(3)により構造寸法が共鳴波長に及ぼす影響を定性的に予測できます。共鳴波長λ_rは、誘電体層の誘電率（ε₃）が大きいほどC_mが増加するため増加します。同様に、幅wが大きいとL_m、L_e、C_mの値が大きくなり、共鳴波長のレッドシフトをもたらします。周囲環境の誘電率（ε₀）の増加はL_mC_m値を増加させますが、式(3)の他項L_eC_mはε₀に依存しません。したがって、共鳴波長はε₀の増加とともに増加します。

結果と考察

以下の構造寸法で議論を開始します。構造はx方向に格子周期P = 580 nmを持ちます。正方形金ナノリボンと三角金ナノリボンの高さはそれぞれt₁ = 45 nmとt₂ = 30 nm。誘電体層、金膜、基板の厚さはそれぞれt₃ = 10 nm、t₄ = 25 nm、t₅ = 165 nm。三角金ナノリボンと正方形金ナノリボンの幅はそれぞれd = 75 nmとw₁ = w₂ = w = 142 nm。図3aに設計構造のシミュレート吸収、反射、透過スペクトルを示します。図3aの通り、吸収効率は95%に達し、751.225 nmで構造の反射率ディップが0.001未満となります。FWHMは1.82 nmで、可視域で報告された従来の狭帯域吸収体[24, 28, 34, 39]よりはるかに狭いです。

a 提案構造の吸収、反射、透過スペクトル。b 共鳴ピークでのMDM構造の電界E分布。c 共鳴ピークでのMDM構造の磁界H分布。d MDM構造と純金属格子構造の反射と吸収スペクトル。e 共鳴ピークでの純金属格子構造の電界E分布。f 共鳴ピークでの純金属格子構造の磁界H分布

吸収ピークの物理メカニズムを解明するため、共鳴ピークでの電界Eと磁界H分布を計算し、図3b, cに示します。図3bの通り、ギャップ内の電界振幅は入射光の35倍に達することが明確です。したがって、提案構造は完全吸収だけでなく、ナノスリットでの電界増強も実現し、これはバイオセンシング応用で重要な現象です。図3cの通り、磁界の大部分は2つの金ナノリボン間の空間に集中し、一部は誘電体層に浸透しており、LSPRに起因する結合効果を示しています。さらに、誘電体層と金膜が超狭FWHMと高吸収性能に及ぼす影響を理解するため、同じ寸法パラメータを持つMDM構造と純金属格子構造の吸収と反射スペクトルを分析・比較しました（図3d）。明らかに、MDM構造はより狭いFWHMとより低い共鳴ディップ反射率を持ちます。金属格子構造の電界と磁界をシミュレートし、図3e, fに示します。明らかに、MDM構造の磁界分布と比較して、金属格子構造の磁界は三角金ナノリボン表面のみに局在し、金属を貫通する磁界がないため、MDM構造と金属格子構造の吸収比較結果を説明できます。さらに、構造内の結合挙動により、図3bの通り、2つの金ナノリボンと薄い金膜間の電界強度は入射波の約40倍で、文献[25]報告値よりはるかに大きいです。

図4aに入射光偏光配置が提案メタマテリアル吸収体の吸収スペクトルに及ぼす影響を示します。TM配置では鋭い吸収ピークを示しますが、TE配置では現れません。明らかに、非対称吸収体構造のため、TE配置入射光ではLSPRを励起できません。また、実際のシステムでは薄い金膜の表面散乱と粒界効果により、薄膜金の減衰定数はバルク金より高くなる可能性があります。薄膜金の減衰定数の影響を考慮するため、図4bに金膜の減衰定数がバルク金の3倍と5倍の場合の吸収スペクトルを示します。明らかに、異なる振幅とFWHMの吸収ピークが観測されます。結果は、金属の材料損失増大が提案狭帯域吸収体の吸収特性向上に不利であることを示し、先行研究[17]と整合的です。

a 提案構造のTEとTM偏光配置での吸収スペクトル。b 金膜減衰定数依存の計算吸収スペクトル

メタマテリアル吸収体の特性は、構造の幾何形状と寸法に強く影響されることが知られています。まず、三角金ナノリボンが設計構造の反射スペクトルに及ぼす効果を調査。三角金ナノリボンを除去、または正方形と半楕円金ナノリボンに変更し（図5c–e）、他パラメータを固定してシミュレーション。これら3構造の反射スペクトルを分析・比較し、元構造と比較しました（図5f–h）。元構造が最も狭いFWHMと最も低い反射率ディップを達成できることが容易に観察されます。これら結果をより理解するため、図5i–lにこれら4構造の共鳴ピークでの磁界（H）分布をプロットし、色で磁界強度を表示。元構造の磁界強度が他3構造より明らかに強いです。これは元構造でLSPRがより効率的に励起されることを意味し、より狭いFWHMとより低い反射率ディップをもたらします。

a–e 異なるナノ構造を持つ提案メタマテリアルの1単位セル模式図。f–h 異なる構造の反射スペクトル。i–l 対応構造の共鳴ピークでの磁界H分布

図5より、三角ナノリボン使用の元構造の光学性能が他構造より優れています。三角ナノリボンが光学性能に及ぼす影響をさらに深く理解するため、元構造の三角ナノリボンと同角度θを持つ台形ナノリボンを含む修正構造（図6a）について詳細計算・解析を行いました。まず、図6b, cに角度θを固定し、台形ナノリボン高さhを変えた場合の修正構造光学性能を調査。明らかに、高さhが10 nm以上では光学性能はほぼ不変で、構造の光学性能が製造上ロバストであることを示します。高さhが5 nm未満では反射率ディップが増加し、これは高さhが小さすぎるとLSPR励起の実効面積が低下するためと説明できます。図6d, eに高さh = 15 nm固定で角度θを変えた場合の光学性能を調査。35°〜68°の広い角度範囲で光学性能はほとんど変化しません。しかし、角度θが30°未満では反射率ディップが明確に増加し、これは小さすぎる角度θがLSPR励起効率を低下させるためと理解できます。このように台形ナノリボンと正方形ナノリボン間の角度パラメータが光学性能に及ぼす影響を詳細に分析することで、元構造の完全吸収性能は三角ナノリボンと正方形ナノリボンのコーナーでのLSPR励起に起因することがわかり、これは図5iの磁界結果とよく一致します。同時に、構造は高さhと角度θの広い範囲で良好な光学性能を維持でき、製造ロバスト性の大幅な緩和を示唆し、実験的観点からナノ構造をより現実的なものにします。最後に、実際のナノ構造製造プロセスを考慮し、図6fに金/誘電体表面粗さと全鋭角のパッシベーション処理を施した構造形状を示します。修正構造と元構造の光学性能比較を計算し、図6gに示します。明らかに、製造公差がナノ構造性能に及ぼす影響は非常に小さく、製造におけるロバストな光学性能を示しています。

a 元構造の三角ナノリボンと同角度θを持つ台形ナノリボンを含む修正構造。b, c 角度θ固定で異なる高さhのナノ構造間反射スペクトル比較。d, e 高さh = 15 nm固定で異なる角度θのナノ構造間反射スペクトル比較。f 金/誘電体表面粗さと全鋭角パッシベーション処理を施した修正構造。g l = 3 nm設定時の修正構造と元構造の反射スペクトル比較

さらに、FDTD法を用いて構造寸法と材料パラメータが、設計構造の反射率ディップ、FWHM、共鳴波長に及ぼす影響を調査。誘電体層の屈折率、金ナノリボン幅w、三角金ナノリボン幅d、金ナノリボン厚さt₁など複数パラメータを検討。図7に誘電体層の屈折率がメタマテリアル構造の反射スペクトルに及ぼす効果を示します。図7aの通り、n_dielectric増加とともに共鳴波長が明確にレッドシフトし、LC回路モデルの予測と整合的。図7bの通り、n_dielectric増加に伴い反射率ディップは最初減少後増加、一方FWHMは狭くなります。反射率ディップとFWHMは、ナノリボンと金膜間の結合強度に強く依存し、ナノリボンと金膜間の誘電体スペーサー材料によって異なる光学性能が生じます。反射率ディップは誘電体層の屈折率が約1.3で最小値を取ります。同時にFWHMは約1.85 nmで、可視域報告狭帯域吸収体[24, 28, 34, 39]よりはるかに狭いです。

a 誘電体層屈折率関数としての反射スペクトル。b 誘電体層屈折率関数としての共鳴ディップ反射率とFWHM

図8に金ナノリボン幅wがメタマテリアル構造の反射スペクトルに及ぼす影響を示します。図8aの通り、金ナノリボン幅wが140から177 nmに変化すると、共鳴波長がブルーシフトし、等価LC回路モデルの結果とよく一致。図8bでは、w増加に伴いFWHMは狭くなり反射率ディップは増加。反射率ディップ増加は、w増加による入射光反射のための実効金属面積増加に起因する可能性。反射率ディップ最小値とFWHM最小値を同時に得ることはできません。しかし、本設計ではwの広い範囲（140〜162 nm）で両値がわずかにしか変化せず、実用応用に有利です。

a 金ナノリボン幅w関数としての反射スペクトル。b 金ナノリボン幅w関数としての反射率ディップとFWHM

さらに、図9aの通り、三角金ナノリボン幅dが55〜75 nmの間では反射率ディップが低い値を維持しますが、dが76 nmを超えると明確に増加。これはナノリボン間距離が大きすぎるとLSPR励起効率が低下し、入射光吸収効率が減少するためと説明できます。FWHMはd増加とともに狭くなり、dの最適サイズは約75 nm。図9bより、金ナノリボン厚さt₁が35〜50 nmで変化すると反射率ディップは低い値を維持しFWHMは狭くなります。しかし、t₁が50〜60 nmに増加すると反射率ディップが明確に増加。これはナノリボンが厚すぎると入射光反射が増加するためと理解できます。図9cより、三角金高さt₂が約30 nmのとき共鳴ディップ最小値を取得。この構造では、三角金高さが15〜40 nmの範囲で反射率ディップが0.025以下であり、優れたロバスト性性能によりメタマテリアル構造設計に有益です。

a 三角金ナノリボン幅d関数としての反射率ディップとFWHM。b 金ナノリボン厚さt₁関数としての反射率ディップとFWHM。c 三角金高さt₂関数としての共鳴ディップ反射率とFWHM

メタマテリアル構造の共鳴波長が環境媒質の屈折率に強く依存することは周知の通りで、これはセンシング応用で広く利用されています。図10aに環境屈折率増加に伴い共鳴波長が明確にレッドシフトすることを示し、これはLCモデルの予測と一致し、同時に反射率ディップが極めて低い値を維持します。RIが1.07から1.12に増加すると、共鳴波長は733.828から755.097 nmへシフト。計算波長感度（S）は約425 nm/RIU、FWHMは1.82 nmまで狭くなります。したがって、FOMは233.5に達します。現時点では、このFOMは可視域で報告された従来のプラズモニック屈折率センサー[24, 28, 34, 39]よりはるかに高いです。提案プラズモニック屈折率センサーは良好な直線性を示します（図10b）。

a 種々の環境屈折率でのプラズモニック屈折率センサー反射スペクトル。b 周囲屈折率に対する共鳴波長シフト

実用応用では、周囲媒質屈折率変化に対する固定波長での相対強度変化検出が必要なことが多く、対応するフィギュア・オブ・メリットはFOM* = max |(dI/dn)/I|で定義されます[17]。図11aの通り、FOM*はw減少とともに明確に変化し、w約358 nmで最大FOM* 1.4 × 10⁵に達します。図11bでは、d減少とともにFOM*が増加し、d = 75 nmで最大。図11cでは、金ナノリボン厚さt₁ = 35 nmでFOM*最大。さらに図11dでは、三角金高さt₂約30 nmでFOM*最大。構造寸法変化に伴うFOMとFOM*の特性を数値的に調査し、高性能プラズモニックセンサー設計への指針を提供します。

a–d 金ナノリボン幅w、三角金ナノリボン幅d、金ナノリボン厚さt₁、三角金高さt₂関数としてのFOMとFOM*

結論

まとめると、可視域で吸収95%に達するほぼ完全な超狭帯域吸収体を提案し、数値的に実証しました。二次元FDTDを用いて、構造形状と寸法がメタマテリアル構造の光学特性に及ぼす影響を詳細に分析。最適化構造寸法では、可視域垂直入射で反射率ディップ0.001、FWHM 1.82 nmを示します。さらにセンシング能力も実証。感度約425 nm/RIU、FOM 233.5を達成。これは可視域で報告された従来センサー[24, 28, 34, 39]よりはるかに優れています。この高いセンシング性能により、本メタマテリアル構造は生物結合、集積フォトディテクタ、化学応用などへの応用が期待されます。

略語

FDTD:

有限差分時間領域法

FOM:

フィギュア・オブ・メリット

FWHM:

半値全幅

LSPR:

局在表面プラズモン共鳴

MDM:

金属-誘電体-金属

S:

感度

TM:

横磁界