異方性誘電体メタサーフェスによる超広帯域テラヘルツ偏光制御

要旨

一般化ホイヘンス原理の枠組みで異方性多重極を持つテラヘルツ誘電体メタサーフェスを実証し、これら多重極の干渉により広帯域・高透過率で巨大な位相シフトを実現。さらに異方性設計によりπ/2から3π/2までの多様な位相遅延が得られ、入射直線偏光テラヘルツ波を右/左旋円偏光、楕円偏光、直交偏光へ変換可能。シミュレーションと実験の双方で楕円率1から-1までの完全なテラヘルツ偏光制御を検証し、テラヘルツメタデバイスの偏光関連応用への道を開いた。

はじめに

偏光は電磁波の状態を定量化する主要パラメータの一つである。特にテラヘルツ領域での偏光制御は、テラヘルツ技術への潜在的応用から大きな研究関心を集めている。しかし、ほとんどのテラヘルツ光源から生成される波は直線偏光であり、複雑な偏光テラヘルツシステムの要件を満たせない。従来の偏光制御手法は複屈折材料を用いるが、大型化と狭帯域動作という本質的欠点があり、現代のコンパクト・広帯域テラヘルツフォトニックシステムへの統合を妨げている。

近年、人工的に設計されたアンテナアレイであるメタサーフェスがテラヘルツ波偏光操作の効率的アプローチとして台頭している。様々な金属メタサーフェスが提案されており、異方性メタサーフェス、カイラルメタサーフェス、多層メタサーフェスに基づく設計が一般的。相変化材料、半導体、二次元材料、液晶、超伝導体などの能動媒質を統合し機能拡張も図られている。多層・能動金属メタサーフェスは性能向上が可能だが、高損失と複雑な製造プロセスを伴う。一方、誘電体アンテナで構成される誘電体メタサーフェスは電磁波制御の新たなアプローチを提供し、電気・磁気ミー共鳴の干渉により高効率で2π位相制御が可能。テラヘルツ偏光制御向け誘電体メタサーフェスの性能向上に多くの努力が払われてきたが、従来は電気・磁気双極子共鳴ベースで位相遅延制御範囲が限定的、原理的に単一周波数動作に留まり、高性能な完全偏光操作を阻害していた。

本研究では多重極を持つ誘電体メタサーフェスを提案し、位相分散を大幅に改善して最大4π位相シフト、巨大位相遅延、広帯域、高効率を実現し、完全テラヘルツ偏光制御を可能にした。楕円形シリコン柱アレイで構成され、異なる電気・磁気多重極を支持。これら多重極の重ね合わせにより、一般化ホイヘンス原理を用いて広帯域で近似完全透過と最大4π位相シフトを達成。シリコン柱の異方性により、短軸・長軸方向で多重極の重ね合わせを独立に制御可能。これにより広帯域巨大位相遅延が実現し、他の金属・誘電体設計を凌駕する性能を示す（補足図S1参照）。単純な設計フレームワークで完全偏光制御が可能なため、メタ原子を人工的に配列し偏光度を空間的に変調し、楕円率変動ベクトル場など複雑なテラヘルツビーム生成にも応用できる。

設計とシミュレーション

誘電体アンテナからの散乱電磁波は、異なる対称性を持つ多重極に分解可能。メタサーフェスとして配列すると、散乱場Eは対称成分E_sと反対称成分E_asの和で表現され、波伝播方向ẑに沿った透過率Tと反射率Rは以下で導出される（式1、2）。高透過・無反射を実現するには、後方へのE_s(-ẑ)とE_as(-ẑ)が同振幅・逆位相である必要がある。特に、対称共鳴（電気双極子など）と反対称共鳴（磁気双極子など）の2つの多重極を支持する場合、これらが同一共鳴周波数・振幅・品質因子を持てば破壊的干渉で零反射が実現される（ホイヘンスメタサーフェス）。しかしこの破壊的干渉は狭帯域のみで生じ、広帯域デバイス実現の根本的制約となる。帯域拡大にはE_sとE_asに高次多重極の寄与を含め、これら多重極間の多重極干渉のバランスで結果的透過を制御する必要がある。これは多重極干渉を用いた一般化ケルカー条件の概念に類似する。

全偏光状態（右/左旋円偏光、楕円偏光、直線偏光）を網羅するには位相遅延が90°-270°（楕円率1から-1）をカバーする必要がある。そこで楕円形シリコン柱アレイからなる異方性誘電体メタサーフェスを提案（図1a）。テラヘルツ領域では吸収損失排除のため本質シリコンを採用。図1aに示すように、多重極干渉が幾何パラメータに応じて異なる組み合わせを維持することで、入射直線偏光を円偏光、楕円偏光、直交偏光へ変換可能。図1bはユニットセルパラメータ：楕円柱の短軸・長軸長a、b、高さh、周期P_x、P_y。図1cは典型的シリコン柱アレイのSEM像（傾斜ビュー、拡大インセット付）。

a 完全偏光制御を実現する誘電体メタサーフェスの模式図。b ユニットセル。c 典型的シリコン柱アレイのSEM像（傾斜ビュー、拡大インセット）

提案メタサーフェスで完全テラヘルツ偏光制御を実現するため、90°-270°位相遅延変化を高透過・広帯域で同時に満たす寸法を数値シミュレーションで最適化。90°-270°の間で45°ステップで5設計（P90、P135、P180、P225、P270）を定義。CST Microwave Studioでシミュレーション実施。シリコンはテラヘルツ領域で損失なし誘電体（ε_Si=11.7）、基板も損失なし（ε_sub=4.5）とモデル化。x、y軸に周期境界条件、z方向にテラヘルツ波入射、背面に開放境界条件。図2aに5設計のシミュレーション透過率と位相遅延を示す。詳細パラメータは補足表S1。全設計で1.2-1.3 THz帯域でx・y偏光入射両方に高透過係数を示し、位相遅延がそれぞれ90°、135°、180°、225°、270°と変化。等しい透過係数で90°位相遅延は左旋円偏光（LCP）変換を、135°、180°、225°、270°は楕円、直交、右旋円偏光（RCP）を実現。提案誘電体メタサーフェスで完全テラヘルツ偏光制御が達成され、帯域制限・低効率・位相遅延範囲限定の既存メタデバイスを凌駕する（文献18、29）。

a シミュレーション、b 実験測定によるP90、P135、P180、P225、P270設計の透過係数と位相遅延

結果と考察

偏光変換性能検証のため、シリコン柱アレイを製造しテラヘルツ時間領域分光（THz-TDS）で評価。製造プロセスは方法論セクション参照。基板には薄いホウケイ酸ガラス（BF33、厚さ300 μm）を使用。270°位相遅延設計の典型サンプルSEM像を図1cに傾斜ビュー・拡大インセットで示す。メタサーフェス性能評価では、透過テラヘルツ波の楕円柱短軸・長軸方向電場をE_x、E_y、裸基板参照での透過電場をE_x(ref)、E_y(ref)とし、透過係数をt_x=E_x/E_x(ref)、t_y=E_y/E_y(ref)で定義。直交偏光間位相遅延をφ=φ_y-φ_x=arg(t_y)-arg(t_x)で算出。測定システム詳細は方法論セクション参照。

測定透過係数と位相遅延を図2bに示す。設計周波数帯域でP90、P135、P180、P225、P270全ケースで高透過係数を取得、対応位相遅延はそれぞれ90°、135°、180°、225°、270°に近い値。シミュレーションと測定にわずかな乖離が見られるが、製造時のサイズ変動や表面粗さによる損失増大、基板効果（損失・反射）が原因と考えられる（補足図S2参照）。それでも変化傾向が類似しており、誘電体メタサーフェスの偏光制御性能を検証している。

偏光変換性能を完全に調査するため、透過波の楕円率を算出（式3）。ストークスパラメータS₀、S₃は透過係数と位相遅延から直接計算可能。図3に示すように、シミュレーションでは楕円率1から-1まで完全にカバー。一般に1.2-1.3 THz付近の偏光変換性能はシミュレーション・実験とも類似傾向。1.4 THz付近で乖離が生じるが、これはシミュレーションでは無損失無限厚基板を仮定するのに対し実験では300 μm厚・損失あり基板を用いること、深部反応性イオンエッチングによる幾何パラメータ変動（柱幅が高さ方向に変化）が多重極を広げ重なりを増やし高Q共鳴を劣化させること、の2要因が複合的に作用するため。これらは低損失基板（石英、ポリイミド、SU8等）の薄膜化と製造プロセス最適化でさらに低減可能。なお動作周波数は一般にオフ共鳴で設計されているため、高Q共鳴劣化の影響を弱く抑えられる。

a シミュレーション、b 実験測定による各誘電体メタサーフェスの楕円率

誘電体メタサーフェス内の多重極干渉を説明するため、球面多重極分解により2直交偏光入射での各多重極散乱断面積（SCS）を算出（文献19、24）。詳細は方法論セクション参照。図4に各設計のx・y偏光入射時のED、MD、EQ、MQ共鳴SCSを示す。P90ではx偏光入射でMD共鳴が1.4 THzに寄与、y偏光では1.18 THzに主寄与。1.42 THz付近高周波域ではy偏光入射でED、EQ、MQ成分が顕著に寄与。x・y偏光SCS比較で、1.2-1.3 THz重なり領域ではオフ共鳴条件が高透過を確保し、異なる多重極間干渉が2直交偏光で異なる位相分散曲線を持ち上げる。多重極間の適切なバランスにより、所望位相遅延を高透過・広帯域で実現可能（本ケースでは90°）。P135、P180、P225でもED、MD、EQ、MQ寄与がP90と類似傾向で共鳴周波数・モード重なりに微妙な変化あり、多重極干渉の偏光制御機能を明確に示す。一方P270では270°位相遅延に巨大位相分散と広帯域高透過が必要で、オフ共鳴条件では困難。そこでP270はイン共鳴条件で設計：x偏光入射でED、MD、MQ共鳴が1.2-1.3 THzに顕著寄与、y偏光ではMD共鳴が1.39 THzで支配。これにより多重極干渉効果で広帯域高透過と270°位相遅延を両立。

ED、MD、EQ、MQ共鳴のSCS多重極分解：a x偏光、b y偏光入射

既存設計と比較し、本提案は単層プラットフォームで完全テラヘルツ偏光制御を実現。位相遅延を90°-270°まで変更可能で、円偏光、楕円偏光、直交直線偏光など多様な偏光状態をカバー。これは他設計では困難（表1）。帯域幅と効率も既存単層設計を上回る。多層設計は性能優位だが、複雑な設計・製造がコンパクトテラヘルツ光学システムへの応用を制限。また本設計は多様な偏光変換を実現するが、ほとんどの多層設計は限定位相遅延・単一偏光変換機能に留まる。

結論

全誘電体メタサーフェスにより、広帯域・高効率で完全テラヘルツ偏光制御を提案・実験実証した。楕円形シリコン柱アレイで構成され、x・y軸で等しい高透過係数を維持しながら位相遅延を90°-270°まで45°ステップで連続調整可能。対応楕円率は1から-1まで変化し、LCP光、楕円偏光、直交偏光、RCP光を完全網羅。多重極分解結果が偏光制御への各多重極寄与を検証。多重極干渉支援誘電体メタサーフェスは、高性能テラヘルツ機能偏光制御デバイス実装へのエキゾチックな戦略を約束する。

方法論

誘電体メタサーフェス製造は標準フォトリソグラフィと深部反応性イオンエッチングで実施。まず抵抗率5000 Ω·cm超の本質シリコンウェハ（厚さ500 μm）を陽極接合でガラスウェハ（BF33、300 μm）に接合。シリコンウェハを180 μmまで薄膜化。アセトンと脱イオン水で30分洗浄。フォトレジストAZ4620をスピンコート、100°Cで10分ソフトベーク。露光時間40秒で楕円アレイをパターニング（MA6）、3分現像。110°Cで5分ハードベーク。深部反応性イオンエッチングで56分シリコンエッチング。最後にアセトン、イソプロパノール、脱イオン水でレジスト除去。

THz-TDSで特性評価。自作スピントロニックテラヘルツエミッタ（800 nm、100 fsパルスレーザー、80 MHz繰り返し）でテラヘルツ波生成。4つのオフアキシ放物面鏡でコリメート・集光。ビーム径約3 mmでサンプル位置に集光。テラヘルツ偏光完全特性評価のため、サンプル前後に2つのテラヘルツ偏光子を配置。検出は電気光学サンプリング法、1 mm厚ZnTe（110）電気光学結晶使用。プローブレーザーはテラヘルツ生成と同一レーザーシステム、プローブパワー20 mW。室温・窒素ガス環境（テラヘルツ領域水吸収除去）で測定。

多重極分解は自社開発Matlabコードで実施。数値シミュレーション結果から楕円シリコン柱内電場分布を抽出。電流密度を導出し、各電流多重極モーメントを分解（式4）。1次・2次モーメント（双極子・四重極子）を計算、高次モーメントは寄与無視。多重極係数a_E(l)、a_M(l)を取得し、散乱断面積を算出（式5）。

データと材料の入手可能性

本研究で使用・分析したデータセットは、責任著者への合理的な要請により入手可能。

略語

SEM:

走査型電子顕微鏡

LCP:

左旋円偏光

RCP:

右旋円偏光

SCS:

散乱断面積

MD:

磁気双極子

ED:

電気双極子

EQ:

電気四重極子

MQ:

磁気四重極子

THz-TDS:

テラヘルツ時間領域分光