グラフェンメタサーフェスによる渦波の生成と操舵の二機能同時実現

概要

原子レベルの厚さを持つ革新的な二次元材料であるグラフェンは、様々な応用において非常に有望な候補として注目を集めています。グラフェンメタサーフェスは、様々な波面の動的制御を可能にし、優れた機能性を実現します。グラフェンメタサーフェスの柔軟性により、多機能デバイスを容易に実装できます。本研究では、渦波の生成と操舵の機能を組み合わせた、新規な多機能グラフェンメタサーフェスの設計を提案します。このメタサーフェスは、グラフェン反射ユニットセルの大規模アレイで構成され、各ユニットセルはそのサイズと外部静電ゲート電圧によって独立に制御されます。グラフェンセルの反射特性を精査することで、多機能性を実現するメタサーフェスを設計しました。シミュレーション結果により、渦波の生成と操舵が可能であることを示しました。この研究は、多機能グラフェンメタサーフェスの設計手法を確立し、グラフェンの調整可能性が再構成可能なグラフェンデバイスの設計・製作への道を開きます。

はじめに

原子レベルの厚さを持つ革新的な二次元材料であるグラフェンは、生物学、光電子工学、テラヘルツ通信などの分野でますます注目されています[1]。テラヘルツ領域では、グラフェンは表面プラズモンポラリトン（SPP）の伝播をサポートするため[2]、従来の貴金属よりも優れた性能を発揮し、テラヘルツ技術における非常に有望な候補となっています。そのため近年、テラヘルツおよび中赤外領域において、変調器[3-6]、検出器[7]、吸収体[8,9]、レーザー[10,11]など、多数のグラフェンベースデバイスが提案されてきました。

電磁波の振る舞いを制御するための再構成可能なメタマテリアルの設計・製作は極めて重要です[12,13]。そのため、様々な周波数帯域で多くの調整機構が実現されてきました[14]：電気的再構成可能メタマテリアル[15]、機械的再構成可能メタマテリアル[16]、非線形材料[17]、液晶[18]、マイクロ流体[19]、半導体構造[20]、グラフェン[21]などです。グラフェンは革新的な材料として、電気/磁気的に制御可能な伝導率を持つことから、微小化可能な制御デバイスの設計・製作に適した傑出した候補です[14,22]。したがって、再構成可能なメタサーフェスの設計に大きな潜在能力があり、その調整可能性を利用した多くの応用が[23]および[24]で提案されています。一般化スネルの法則[25,26]を適用することで、グラフェンメタサーフェスにより異常反射を調整・実現できます[27]。これらの研究は、調整可能なテラヘルツデバイスの設計・製作への道を開いています。

通信において、軌道角運動量（OAM）は無限の状態を提供できるため、チャネル容量の向上に重要です[28,29]。三次元メタマテリアルはOAM波の生成に利用できます[30]。二次元メタマテリアルと見なせるメタサーフェスは、サブ波長厚で優れた性能を発揮します。マイクロ波領域では、メタサーフェスは様々な偏光と利得特性を持つ波を生成するサブ波長デバイスの設計・製作に広く利用されています[31-34]。テラヘルツ領域では、反射型グラフェンメタサーフェスによる調整可能な渦波生成が報告されています[35]。グラフェンメタサーフェスは波面制御の柔軟性を持つ[36]ため、渦波生成と異常反射の機能を組み合わせた実現可能な設計により、渦波の指向性を高精度に調整できることが期待されます。

本研究では、マイクロナノ光学におけるメタサーフェスに関する先行研究[37-41]に基づき、2つのメタサーフェスの機能を組み合わせるメカニズムを研究します。グラフェンセルを解析し、反射係数とその化学ポテンシャルおよびパッチサイズの関係を取得します。グラフェンメタサーフェスを設計するための参照として、完全な360°の反射位相範囲を較正しました。組み合わせメタサーフェスは、反射型グラフェンセルの大規模アレイによって実現されます。シミュレーション結果により、渦波の生成と一定角度での反射による操舵が可能であることを示しました。

手法

グラフェンの伝導率は、バンド間遷移とバンド内遷移で構成されます。バンド内遷移がテラヘルツおよび赤外領域を支配し、バンド間遷移が可視光領域を支配します。テラヘルツおよび赤外領域では、伝導率をドゥルーデモデル[24]でモデル化できます：

$$ \sigma(\omega)=\frac{2e^{2}}{\pi\hbar^{2}}k_{B}T\cdot\ln\left[2\cosh\left(\frac{E_{f}}{2k_{B}T}\right)\right]\frac{i}{\omega+i\tau^{-1}}, $$

ここで、$k_B$はボルツマン定数、$T$は温度、$\tau$は緩和時間、$E_f$はフェルミエネルギーです。

本研究ではデバイスがテラヘルツ領域で動作し、$E_f \gg k_B T$であるため、室温$T=300K$、グラフェンの緩和時間$\tau=1$ psとして、式を以下のように簡略化できます：

$$ \sigma(\omega)=\frac{e^{2}E_{f}}{\pi\hbar^{2}}\frac{i}{\omega+i\tau^{-1}}, $$

フェルミエネルギー$E_f$は外部静電ゲート電圧で制御されます。シミュレーションでは、グラフェンを原子厚のため3Dメタマテリアルブロックではなく、2D表面伝導条件としてモデル化します。

グラフェンメタサーフェスはグラフェンセルの大規模アレイで構成され、表面上で集団的なプラズモニック振る舞いが励起され、特異な電磁特性を実現します。周波数は1.3 THzであり、プラズモニックモードに伴う遅い波の伝播により、共振が非常に小さなサイズ、すなわち$\lambda$/10以下で発生します[23,42]。グラフェンセルのメタサーフェスを設計するため、単一グラフェンセルの反射振る舞いの較正グラフを抽出し、各パラメータの詳細な影響を研究します。

典型的なユニットグラフェンセル（図1参照）は、原子厚のグラフェンパッチが上部に配置された多層構造で構成されます。サイズ$w_x \times w_y$のグラフェンパッチは、一辺$p=14\mu m$の正方形基板スタックの中央上部に配置されます。底部の金属接地層の上に、厚さ25μmの石英基板（$\epsilon_r=3.75, \tan\delta=0.0184$）が配置されます。グラフェンパッチと50nm厚の多結晶シリコン層の間に外部バイアスDC電圧を印加します。その間に10nm厚のAl₂O₃（アルミナ、$\epsilon_r=8.9, \tan\delta=0.01$）スペーサー層を挿入します。外部バイアスDC電圧を0〜14.7Vで制御することで、化学ポテンシャルを0.01〜1.0eVの範囲で調整可能です[23,35]。多結晶シリコン層とアルミナスペーサーは本論文のシミュレーションではモデル化していません。理由は以下の通りです。第一に、これらの層の厚さが石英基板よりはるかに薄いため、反射振る舞いへの影響は無視できることを示す別途の2Dシミュレーション（計算コストが大幅に低い）を実施しました。第二に、有限要素シミュレーションでは、サイズ差が大きい隣接物体を扱う際に極めて多くの要素が必要となり、これら2層をモデル化した3Dシミュレーションは極めて高コストになるためです。

周期境界条件を$x$および$y$方向に割り当て、$x$偏光（電場が$x$方向に偏光）で正入射させることで、$w_x$と化学ポテンシャル$\mu_c$の影響を精査しました。$w_y$は$x$方向の伝導率に無視できる影響しか与えないため、全シミュレーションで4μmに固定しました。$w_x$を0.2〜13.8μm（0.2μmステップ）、$\mu_c$を0.01〜1.00eV（0.01eVステップ）でスイープし、周波数を1.3THzに固定しました。$S_{11}$の位相と振幅を図2にプロットし、これらを較正グラフと呼びます。メタサーフェスの効率を保証するため、反射係数の大きさが0.7以上の領域のみを有効とし、それ以外をブランクとしています。較正グラフでは、360°の完全な位相カバレッジが得られ、グラフェンメタサーフェスの構築に十分です。

位相図は位相を精密に制御するため十分に滑らかである必要があります。図3に示すように、$w_x$と$\mu_c$の7つの組み合わせを選択し、0°から360°までの完全な位相カバレッジを実現しました。

結果と考察

様々な機能を実現するため、2つのメタサーフェスの機能を組み合わせたり、既存のメタサーフェスに新機能を追加したりすることは非常に有用です。この方法論は、新しいメタサーフェスを設計するための汎用的な方法を提供します。本論文では、渦波生成と異常反射による波の偏向の機能を組み合わせます。

以下に、2つのメタサーフェスMS₁とMS₂を1つの多機能メタサーフェスMSₜに組み合わせる一般化された方法論を提案します。組み合わせを実現するため、一般化反射の法則[25]から始めます。図4に示すように、自由空間波長$\lambda$の平面波が入射角$\theta_i$で入射する場合、一般化反射の法則は以下の式で記述されます：

$$ \sin\theta_r - \sin\theta_i = \frac{\lambda}{2\pi n_i} \frac{d\phi}{dx}, \quad (1) $$

ここで、$\theta_r$は反射角、$n_i$は上部空間の屈折率、$\phi(x)$は界面に沿った位相不連続性を表します。

波が正入射し、上部空間が自由空間（$n_i=1$）である簡略化されたケースを考えます（図5）。最初の2つのメタサーフェスMS₁とMS₂について、式(1)はさらに簡略化されます：

$$ \frac{d\phi_m}{dx} = \frac{2\pi}{\lambda} \sin\theta_{rm}(x) \quad m=1,2. \quad (2) $$

MSₜの$\phi_t$を得るため、界面に沿ったセグメント$D_x$を選び、以下の問題を考えます：$x \in D_x$において$-\pi/2 < \theta_{r1}(x) + \theta_{r2}(x) < \pi/2$が成り立つとき、全ての$x \in D_x$に対して以下を満たす$\phi_t$を見つけます：

$$ \frac{d\phi_t}{dx} = \frac{2\pi}{\lambda} \sin\theta_{rt}, \quad \theta_{rt}(x) = \theta_{r1}(x) + \theta_{r2}(x). \quad (3) $$

式(2)と(3)から、以下が導かれます：

$$ \phi_t(x) = \cos\theta_{r2}\phi_1(x) + \cos\theta_{r1}\phi_2(x) - \int_{D_x} \left( \sin\theta_{r2} \frac{d\theta_{r2}}{dx} \phi_1 + \sin\theta_{r1} \frac{d\theta_{r1}}{dx} \phi_2 \right) dx. \quad (5) $$

ここで、積分項は$\theta_{ri}(x)$の変化の寄与を計算し、主に数値的に計算できます。式(5)は、2つのメタサーフェスの機能を組み合わせる上で重要な役割を果たします。

さらに、操舵角が一定であれば、式(6)の積分項は消えます。式(5)は以下のように大幅に簡略化されます：

$$ \phi_t(x) = \cos\theta_{r2}\phi_1(x) + \cos\theta_{r1}\phi_2(x) + C. \quad (6) $$

これがメタサーフェスを組み合わせる支配方程式であり、位相分布を計算して渦波生成と異常反射を組み合わせることができます。

本論文では、MS₁を渦波生成メタサーフェス、MS₂を波操舵メタサーフェスとします。

[35]で示されているように、モード$l$の渦波は、位相シフトが連続的に増加する$N$セクタのプレートによって生成できます。第$n$セクタの位相シフト$\phi_n$は$\phi_n = \phi_0 + 2\pi n l / N$で計算され、$\phi_0$は初期セクタの位相シフトです。渦波を生成するには、$-N/2 < l < N/2$を満たす必要があります。したがって、$N=4$でモード$l=0, \pm1$を生成するのに十分です。

$l=1$の渦波を生成するため、プレートを図6aのように4セクタに分割します。位相条件$\phi_1(x,y)$は、セクタを通じて反時計回りに90°ずつ減少する区分的定数関数です：

$$ \phi_1(x,y) = \begin{cases} 0^\circ & x\geq 0, y\geq 0 \\ -90^\circ & x<0, y\geq 0 \\ -180^\circ & x<0, y<0 \\ -270^\circ & x\geq 0, y<0 \end{cases} \quad (7) $$

$x$偏光波が上から正入射すると、$l=1$の渦波が反射されます。波は垂直に反射されるため、偏向角は0°、すなわち$\theta_{r1}(x)=0^\circ$です。

偏向角$\theta_r$を持つ異常反射を生成するには、式(1)を適用します。図4のように、自由空間で正入射（$\theta_i=0^\circ, n_i=1$）の場合、式(1)は以下のように簡略化されます：

$$ \phi_2(x) = \frac{2\pi \sin\theta_r}{\lambda} x + C. $$

本研究では偏向角を$\theta_r=30^\circ$と設定します。上式より、ユニットセルの周期が14μmであることを考慮すると、隣接パッチ間の位相シフト差は10.9°と計算されます。位相分布を図6bに示します。

MS₁とMS₂を組み合わせるため、$\theta_{r1}(x)=0^\circ$と$\theta_{r2}(x)=30^\circ$を式(6)に代入し、MSₜの設計式を得ます：

$$ \phi_t(x) = \frac{\sqrt{3}}{2} \phi_1(x) + \phi_2(x) + C. $$

この式から位相分布を計算し、図6cに示します。図3に従い、各セルの化学ポテンシャル$\mu_c$とパッチサイズ$w_x$を選択することで、32×32のグラフェンメタサーフェスを構成します。図1aはメタサーフェス上のグラフェンセルの配置を示します。各セクタは16×16のサブドメインで、垂直方向に16列からなります。各列は16個の同一グラフェンパッチで構成され、特定の$w_x$と$\mu_c$の組み合わせが割り当てられます。

プレートは上から入射する$x$偏光波で励起されます。入射波の電場は正規化されています（$\vec{E}_{inc} = \vec{x}$）。シミュレーションは商用有限要素ソルバーCOMSOL Multiphysics 5.2で実施しました。グラフェンは原子厚ですが、基板の厚さはマイクロメートルオーダーです。グラフェンパッチに3Dメッシングを適用すると計算コストが膨大になるため、グラフェンパッチの厚さを無視し、COMSOL Multiphysicsで遷移境界条件として等価な2D表面伝導条件を適用しました。プレート上には32×32のパッチがあり、4セクタに分割され、各セクタに16×16のパッチが配置され、それぞれサイズと化学ポテンシャルで独立制御されます。シミュレーションは710万自由度を消費し、40×2.1GHzスレッド、256GBメモリのサーバで実行されました。

図7bは、入射波で正規化された反射波の電場大きさを示します。グラフェンメタサーフェスは、$l=1$の渦波を生成し、$x$軸方向に30°偏向させます。

結論

要約すると、多機能グラフェンメタサーフェスの設計原理を研究しました。2つのメタサーフェスを組み合わせる方法論を提案しました。例として、渦波生成と波操舵の機能を組み合わせたグラフェンメタサーフェスを設計しました。グラフェンは二次元原子厚材料であり、外部ゲート電圧を印加することで位相条件を動的に調整できます。そのパラメータを精査し、単一グラフェンセルの反射振る舞いを較正して、360°の位相シフトカバレッジを取得しました。32×32ユニットセルからなるグラフェンメタサーフェスを設計し、異常反射と渦波THz波生成を同時に実現しました。シミュレーション結果により、$l=1$の渦波が生成され、操舵されることを示しました。グラフェンはテラヘルツ領域でSPPサポート、高効率、調整可能性などの特異な振る舞いを示すため、テラヘルツ技術における有望な候補です。本研究は、グラフェンによって実装される異なるメタサーフェスの機能を組み合わせるアプローチを調査し、テラヘルツ領域における動的に制御可能な多機能メタサーフェスへの道を開きました。

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

本研究で生成および/または分析されたデータセットは、対応著者から合理的な要請に基づき入手可能です。

略語

OAM:

軌道角運動量

SPP:

表面プラズモンポラリトン