二層ツイスト半円ナノ構造を用いた選択的円偏光対応高性能プラズモニックカイラルメタサーフェス吸収体

概要

本論文では、光学周波数において右回り（RCP）および左回り（LCP）円偏光に対して高い選択的吸収を達成するプラズモニックカイラルメタサーフェス吸収体（CMSA）を提案する。このCMSAは、誘電体基板を四重ツイスト半円金属ナノ構造の二層で挟んだ構造を持つ。提案されたCMSAは強い選択的吸収帯を示し、LCPおよびRCP光の吸収ピークが異なる共鳴周波数で現れ、顕著な円二色性（CD）効果の存在を反映している。LCP光に対する吸収率は93.2%、RCP光に対しては91.6%に達し、最大CD大きさはそれぞれ288.5 THz付近で0.85、404 THz付近で0.91となる。単位セルナノ構造の電場分布により、強いカイラル光学応答のメカニズムを解明した。さらに、提案CMSAの幾何学的パラメータが円偏光選択的吸収特性に及ぼす影響を体系的に調査した。

導入

カイラリティは、鏡像対称性や反転対称性を持たない物体の幾何学的特性として普遍的に存在し、科学技術において常に重要なテーマである[1,2]。カイラル媒質は通常、互いに鏡像関係にあり並進や回転では重ね合わせ不可能な2つのエナンチオマー形態をとり、RCP光とLCP光に対して異なる応答を示す[1]。カイラル媒質に由来する円偏光（CP）の円二色性（CD）は、最も独特なカイラル光学特性の一つである。CD効果とは、カイラル媒質におけるRCP光とLCP光の応答差であり、生物学、医学、化学、偏光関連光電子デバイスへの幅広い応用が期待される[3-5]。しかし、天然材料のCD効果は極めて弱く、実用化を大きく妨げている。メタサーフェスは、単層または数層の平面構造からなるメタマテリアルの一種であり、電磁波や光の任意な操作に大きな可能性を示している[6-10]。特にカイラルメタサーフェス（CMS）は、カイラル光学効果を桁違いに増強できる。

CMSは、負の屈折率や光学活性[11,12]、非対称透過[13,14]、巨大CD効果[15-17]、偏光変換[18,19]、波面操作[20,21]などの異常な電磁特性を示すため、大きな注目を集めてきた。これまでに、スプリットリング、らせんワイヤ、ガンマディオン、L字型など様々なCMS構造が提案され、LCPまたはRCP光に対する高効率カイラル選択的場増強が実現されてきた[22-33]。しかし、従来のCMS設計の多くは透過特性におけるカイラリティ性能に焦点を当てており、工学応用で同等に重要なCP光の反射・吸収特性については十分な検討がなされてこなかった。従来の吸収体研究の大半は直線偏光に適用されており、CP光向けの設計はほとんど研究されていない。実際、CMSを用いてCP光用の新規吸収体を構築することも可能である[25,26,29-33]。例えば、Liら[25]はL字型ワイヤベースの超薄膜吸収体を提案したが、マイクロ波領域でLCP波のみ減衰可能だった。Wangら[29]は、カイラルメタミラーにより赤外領域でLCP光をほぼ全反射し、RCP光を完全吸収することを示した。Tangら[30]は、ŋ字型共振子を用いた吸収体で可視光領域での異なるCP光の選択的吸収を達成した。その後、プラズモニックメタサーフェスを用いた近赤外カイラル吸収体が提案され、LCPまたはRCP光を選択的に吸収することが実証された。しかし、ほとんどのCMSの吸収率は90%未満に留まっている。したがって、高い選択的吸収を持つ効果的なカイラルメタサーフェス吸収体（CMSA）の設計が強く求められている。

本研究では、近赤外・可視領域で動作する二層四重ツイスト半円ナノ構造ベースの高効率CMSAを提案する。このCMSAは、異なる共鳴周波数において異なるヘリシティのCP光に対して90%以上の選択的吸収を達成できる。提案CMSAの強い選択的吸収により、約0.9の高いCD値が実現される。CP光の選択的吸収の物理メカニズムを電場分布により詳細に解析した。さらに、単位セル幾何パラメータが選択的吸収に及ぼす影響を体系的に調べた。本成果は、熱吸収体、光通信デバイス、光検出器、光フィルタ、イメージング、ホログラムなど多くの実用応用における強い吸収とCD効果を持つCMSA設計の指針となると考えられる。

単位セル設計

図1に提案CMSAの模式図を示す。これはツイスト半円ナノ構造の周期アレイで構成される。誘電体基板の両面に配置された四重ツイスト半円ナノ構造は、それぞれ隣接するものと90°回転して配置され、底面側の半円ナノ構造も上面側に対して90°回転している（図1b）。先行研究[32]と同様に、上面の4つの半円ナノ構造は銅シリンダーで底面のものと接続され、その半径は半円ワイヤ幅と同一で、導電的結合を増強する。ツイスト半円ナノ構造は結合共振子系と見なせ、2つのツイスト接続半円間の電気的・磁気的誘導結合から強いカイラル応答が生じる[34,35]。上面と下面に鏡対称性を持つこの単純なツイスト半円ナノ構造により、提案CMSAのカイラリティが増強される。

提案CMSAの全単位セルは、光伝播方向に対して四重回転（C₄）対称性を示す。中間層の誘電体基板は損失のないMgF₂（相対誘電率1.9）で作られる。カイラル金属ナノ構造には金を用い、その材料パラメータはDrudeモデル[36]で記述される。最適化された単位セル構造パラメータは以下の通り: pₓ = pᵧ = 600 nm, r = 70 nm, w = 40 nm, tₛ = 120 nm, tₘ = 30 nm。単位セルはx軸およびy軸方向に600 nmの周期で配置され、入射光周波数が500 THzに達しても回折を回避する。提案CMSAの効率を検証するため、CST Microwave Studioの周波数領域ソルバーを用いた有限要素法（FEM）ベースの全波高周波電磁シミュレーションを実施した。単位セルナノ構造、適切な境界条件、メッシュサイズ、周波数範囲を適切に設定した後、周波数領域シミュレーションを実行した。

結果と考察

法線入射LCPおよびRCP光に対する共偏光透過係数（t₊₊(ω), t₋₋(ω)）と反射係数（r₊₊(ω), r₋₋(ω)）を図2に示す。2つのカイラルプラズモニック共振モード（モード1とモード2）が、それぞれ周波数 f₁ = 288.5 THz と f₂ = 404 THz で明確に観測される。RCP光の共偏光反射係数 r₊₊(ω) と LCP光の r₋₋(ω) は等しく、関心周波数範囲全体で0.4未満である。さらに、共鳴では両者の大きさが約0.15まで低下し、CMSAと自由空間のインピーダンス整合がRCPおよびLCP光の両方で達成されていることを示す。一方、共偏光透過係数 t₊₊(ω) （RCP）と t₋₋(ω) （LCP）は共鳴で大きく異なり、提案CMSAのカイラル性を反映している。低周波数 f₁ 付近では、RCP光の t₊₊(ω) 大きさは約0.93と高いが、LCP光の t₋₋(ω) は約0.075に過ぎない。高周波数 f₂ 付近では、RCP光の t₊₊(ω) が最小値0.018まで減少する一方、LCP光の t₋₋(ω) は最大値約0.92に達する。つまり、低周波数では入射RCP光のみがCMSAを透過しLCP光は遮断され、高周波数 f₂ では逆にLCP光のみが透過しRCP光が遮断される。このカイラル選択現象は、RCP光とLCP光の異なる吸収をもたらし、共鳴における高効率選択的吸収と巨大CD効果の存在を意味する。

図2bは、入射LCPおよびRCP光に対する吸収率スペクトル（A₊(ω), A₋(ω)）を示す。LCPおよびRCP光の吸収率はそれぞれ約93.2%と91.6%の最大値に達し、一方でRCPおよびLCP光の吸収率は上記2共鳴付近でそれぞれ約8.7%と4.8%まで低下する。したがって、提案CMSAは低周波数 f₁ でLCP光を強く吸収しRCP光を高透過させ、高周波数 f₂ では完全に逆の挙動を示す。つまり、特定のヘリシティを持つ2つのCP光に対して異なる共鳴で選択的吸収を示し、もう一方を反射する。また、単一サイズのカイラルナノ構造のみで2つの強い吸収周波数帯を持つことは、異なる幾何サイズに依存して特定のCP光に適応する単一孤立吸収帯を持つ従来のカイラル吸収体[25,26,29,31-33]に比べて明らかに優れている。したがって、設計されたカイラルナノ構造は、低周波数では完全なLCP光吸収体、高周波数では完全なRCP光吸収体として機能する。なお、提案CMSAの選択的吸収性能は入射角の増大（斜め入射）とともに悪化する（高次多極プラズモン共鳴のため、データ未提示）。さらに、CP光に対する高いカイラル選択的吸収は、提案CMSAにおいて巨大CD効果をもたらすと推測される。

LCP光とRCP光の吸収または透過の差は、CDパラメータ Δ で特徴付けられる。図3aにCMSAのCDスペクトルを示す。CDパラメータの主ピークは、2つの共鳴周波数でそれぞれ約0.85と0.91であり、報告されているカイラルナノ構造[17-42]よりはるかに大きい。巨大CD効果はCMSAの強いカイラリティによるもので、透明な円偏光子として応用可能である。さらに、CMSAを円偏光子として用いた場合のCP純度を調べるため、楕円率角 η と偏光方位回転角 θ を図3bに示す。η の値は低周波数で約40.4°、高周波数で約-43.9°、θ の値は両共鳴で約0°である。これは、2つの共鳴周波数でCMSAスラブを透過した光が顕著なRCPおよびLCP特性を示すことを意味する。このCMSAベースの円偏光子は、単位セルの高いC₄対称性により、任意の偏光入力に対して有効である。したがって、設計されたカイラルナノ構造により均質な円偏光子が実現可能と考えられる。

CMSAの選択的吸収と巨大CD効果を完全に理解するため、CP光の透過・反射係数から標準的な抽出手順[43,44]により屈折率 Re(n), Re(n₋), Re(n₊) およびカイラルパラメータ Re(κ) を抽出し、図4a,bに示す。設計CMSAには強いカイラリティに関連する2つの共鳴が現れる。低周波数共鳴は約288.5 THz、高周波数共鳴は404 THz であり、選択的吸収とCDピークの特性周波数と一致する。図4aに示すように、Re(n) は最大絶対値 -2.3 と -1.1 で負となり、Re(κ) は上記2共鳴周波数付近で最大絶対値 6.4 と -5.1 に達する。カイラルパラメータ κ も、n± = n ± κ の関係により、RCP/LCP光の屈折率を共鳴で負に押し上げることに寄与する。したがって図4bに示すように、LCP光の Re(n₋) と RCP光の Re(n₊) はそれぞれ286.2-291 THz と 400.2-404 THz で負となり、最大負値はそれぞれ -8.6 と -6.3 に達する。これは、提案CMSAの高い選択的吸収と巨大CD効果が、LCP光とRCP光の負屈折特性に関連していることを示唆する。

提案CMSAの選択的吸収と巨大CD効果の起源をさらに解明するため、288.5 THz と 404 THz におけるRCP光およびLCP光駆動時の単位セル電場（E_z）分布を調べた。半円ナノ構造の直径が異なるCP光の入射波長に比べて小さいため、表面プラズモン共鳴の励起により振動双極子場が生成されることが知られている[45-48]。RCPまたはLCP光が半円ナノ構造に照射されると、提案CMSAにおいて選択的吸収と巨大CD効果が生じ、各層の電場および磁場成分分布が異なることが合理的に推測される[48-53]。

図5に、異なる共鳴周波数におけるRCP光およびLCP光駆動時の提案CMSA単位セルの電場（E_z）分布を示す。上面および下面の半円ナノ構造上の電場分布の詳細プロットは、各表面プラズモンモードの性質を明確に示す[54]。上面および下面の半円ナノ構造上の赤と青の領域は、RCPおよびLCP光励起下での正負の電荷蓄積を表す。正負電荷は各半円ナノ構造の隅に主に分離蓄積し、電気双極子振動のように振る舞う。設計された半円ナノ構造では電気双極子パワーが磁気双極子よりもはるかに強く、電気双極子振動が支配的であることがわかる。共鳴で生成される選択的吸収と巨大CD効果は、LCP光とRCP光励起下での明確な双極子パワー差に起因する。ここでは、各層の4つの半円ナノ構造の電荷振動を1つの双極子振動とみなす等価電気双極子モーメントによる簡略化手法を用いる[48-50]。Born-Kuhn理論[47,48]によれば、同じ電場方向を持つ2つの双極子からハイブリダイズされたモードを結合モード、90°または直交方向を持つ2つの双極子からハイブリダイズされたモードを反結合モードと呼ぶ。図5a1,b1に示すように、共鳴周波数 f₁ = 288.5 THz でRCP光照射下では、上面と下面の電気双極子場が直交方向を示し反結合モードを形成し、Born-KuhnモデルによればRCP光の高透過が生じる。一方、図5c1,d1のLCP光照射下では、CMSAの電場分布が上面と下面の同方向の等価電気双極子モーメントからなる結合モードのハイブリッドと見なせ、LCP光の高吸収をもたらす。したがって、結合モードと反結合モードは、低周波数におけるLCP光とRCP光照射下での異なる共鳴エネルギーと異なる透過・吸収を引き起こす（図2参照）。図5a2,b2,c2,d2に示すように、共鳴周波数 f₂ = 404 THz では、RCP光照射下で上面と下面の電気双極子場が同方向（結合モード）、LCP光照射下で直交方向（反結合モード）を示し、それぞれRCP光の高吸収とLCP光の高透過をもたらす。したがって、2つの異なる周波数での選択的吸収とCD効果は、主に上面と下面の電気双極子モーメントのハイブリッド結合によって誘起される結合モードと反結合モードに起因する。

以下では、単位セルの幾何パラメータが提案CMSAの吸収特性に及ぼす影響を調べる。図6に、単位セルの異なる幾何パラメータ（r, w, tₘ, tₛ）におけるLCP光およびRCP光の吸収スペクトルと共鳴周波数（f₋, f₊）を示す。設計ナノ構造では、明らかに多因子競合的で複雑な、パラメータ依存的選択的吸収特性のスペクトル変化が観測される。本研究では、制御群の幾何パラメータを r = 70 nm, w = 40 nm, tₘ = 30 nm, tₛ = 120 nm とし、一度に1パラメータずつ変更した。

まず、他のパラメータを固定して半円ナノ構造の半径 r を変化（65, 70, 75, 80 nm）させた場合（図6a1,a2）。r が増加すると、LCP光の f₋ と RCP光の f₊ がともに徐々に減少する。これは等価 LC 共振回路理論[55,56]で解釈できる。共鳴周波数は f∓ = 1/(2π√(LC)) で表され、等価容量 C とインダクタンス L は主にCMSAの幾何パラメータで決まる。r が増加すると L が増加し、共鳴周波数が減少する。また図6a1に示すように、r 増加に伴いLCP光の吸収率は徐々に減少するが、RCP光の吸収率はほぼ変化しない。次に、w を30-45 nm（ステップ5 nm）で変化させた場合（図6b1,b2）。w 増加に伴い、LCP光とRCP光の共鳴周波数 (f₋, f₊) はともに徐々に増加する。これは主に C の減少による。LCP光の吸収率は最初増加してからわずかに減少し、RCP光の吸収率は徐々に減少する（図6b2）。次に、tₘ を20-50 nm（ステップ10 nm）で変化させた場合（図6c1,c2）。w の場合と同様に、tₘ 増加に伴いLCP光の f₋ が顕著に増加し、RCP光の f₊ がわずかに増加する。tₘ 増加で L が減少するため共鳴周波数が増加する。また、LCP光とRCP光の吸収率はともに最初増加してから減少する（図6c2）。最後に、tₛ を110, 120, 130, 140 nm で変化させた場合（図6d1,d2）。tₛ 増加に伴いLCP光の吸収率は徐々に増加し、RCP光の吸収率はわずかに減少する（図6d1）。また、LCP光とRCP光の共鳴周波数 (f₋, f₊) はともに徐々に減少する（図6d2）。この場合、tₛ 増加で C が増加し、共鳴周波数が減少する。以上より、設計カイラルナノ構造単位セルの幾何パラメータに対して、LCP光とRCP光の共鳴周波数 (f₋, f₊) と吸収レベルは敏感であると結論付けられる。したがって、提案CMSAの選択的吸収特性は構造パラメータを変化させることで動的に調整可能である。

結論

結論として、二層四重ツイスト半円ナノ構造ベースのCMSAを提案し、近赤外および可視領域においてRCP光とLCP光に対するほぼ完全なカイラル選択的吸収と巨大CD効果を達成した。シミュレーション結果では、RCP光とLCP光の両方に対するカイラル選択的吸収率が90%以上、CD大きさが0.91に達することを示した。抽出された有効電磁パラメータによると、低周波数吸収とCD効果はLCP光の負屈折特性に関連し、高周波数ではRCP光の負屈折特性に関連する。電場分布から、CMSAのカイラル選択的吸収特性と巨大CD効果は、主に上面と下面の電気双極子モーメントのハイブリッド結合によって誘起される結合モードと反結合モードに起因することが示された。さらに、CMSAの共鳴周波数とカイラル選択的吸収レベルは、単位セルの幾何パラメータ変更により調整可能である。したがって、CMSAの設計は光フィルタ、カイラルイメージング、円偏光子、検出、光通信などの将来応用に有望であると結論付けられる。

数値計算手法

FEMシミュレーション: 有限要素法（FEM）ベースの全波電磁シミュレーションを実施した。シミュレーションでは、金の電気的特性をDrudeモデル[36]で記述した:

ε_Au = 1 - ω_p² / [ω(ω + iγ)] (1)

ここで ω_p = 1.37×10¹⁶ rad/s はプラズマ周波数、γ = 8.04×10¹³ rad/s は光学周波数領域における金の衝突周波数である。シミュレーションでは、x軸およびy軸方向に単位セル境界条件を適用し、2つのCP固有光を直接励起源とした。広帯域CP光を励起源として、設計カイラルナノ構造の単位セルを -z から +z 方向に法線入射させた。その後、LCP光とRCP光の反射係数と透過係数を取得した。一般に、LCP/RCP光の吸収率 A₋(ω)/A₊(ω) は[17,32]以下のように表される: A₋(ω) = 1 - R₋₋(ω) - T₋₋(ω) = 1 - |r₋₋(ω)|² - |t₋₋(ω)|², A₊(ω) = 1 - R₊₊(ω) - T₊₊(ω) = 1 - |r₊₊(ω)|² - |t₊₊(ω)|²。共偏光透過係数 t₋₋(ω) はLCP光、t₊₊(ω) はRCP光用、r₋₋(ω) と r₊₊(ω) はそれぞれ共偏光反射係数である。設計カイラルナノ構造単位セルの高い C₄ 対称性により、クロス偏光透過係数（t₊₋(ω), t₋₊(ω)）と反射係数（r₊₋(ω), r₋₊(ω)）は0.01未満で無視できる。また、CD効果は2つのCP光の選択的吸収により誘起され、Δ = |t₊₊(ω)| - |t₋₋(ω)| [14,29] で表される。楕円率と光学活性は、設計カイラルナノ構造のカイラリティを評価する重要なパラメータである。楕円率はカイラルナノ構造の透過光の偏光状態を特徴づけ、楕円率角 η = arctan[(|t₊₊(ω)| - |t₋₋(ω)|) / (|t₊₊(ω)| + |t₋₋(ω)|)] で記述される。一方、光学活性は透過直線偏光の偏光面回転特性を表し、偏光方位回転角 θ = [arg(t₊₊(ω)) - arg(t₋₋(ω))]/2 で記述される。