本プロトコルでは、6G無線通信アプリケーション向けにテラヘルツ帯（0.5～30 THz）で動作する単層カーボンナノチューブベースの光再構成可能知的反射面の正確な電磁界モデリングを実現するため、ハイブリッド人工ニューラルネットワークとZ変換を用いた手法を提示し、180× 反射位相調整能の向上による計算加速 310° およびビームステアリング範囲は ±45°.
方法論記事
テラヘルツ波ビームステアリングに向けたカーボンナノチューブベース再構成可能知能表面のモデリングのためのハイブリッドANN-Z法
この記事について
サマリー
要約
単層カーボンナノチューブに基づく光再構成可能インテリジェント表面は、次世代6G無線システムにおけるテラヘルツ波のビームステアリングおよびフォトニック波制御のための有望なソリューションを提供します。しかし、広帯域にわたる量子輸送効果、速度論的インダクタンス、および多共鳴励起子挙動により、これらの構造の正確なモデリングは依然として困難です。本プロトコルでは、単層カーボンナノチューブ光再構成可能インテリジェント表面の単位セルの正確な特性評価を実現するため、久保公式を用いた量子導電率モデリング、多項式回帰に基づくデータ平滑化、およびZ領域伝達関数解析を統合したハイブリッド計算フレームワークについて解説します。まず、石英基板上に、カイラリティー(10,5)（直径0.60 nm、バンドギャップ1.762 eV）の単層カーボンナノチューブ・ナノストリップ共振器を交差するように設計し、次にCST Microwave Studioを用いて0.5～30 THz帯におけるフルウェーブ電磁界シミュレーションを行います。抽出されたSパラメータは、8次の多項式回帰モデルによって処理され、数値的変動が除去され、平滑化された電磁応答が予測されます。その後、QR分解を用いた最小二乗最適化により、分子6次、分母7次の離散伝達関数H(z)をフィッティングさせ、極零点安定性解析および受動性の検証を可能にします。さらに、本プロトコルでは、ビームステアリングの最適化に向けた化学ポテンシャルの変調による量子導電率のチューニングを組み込んでいます。代表的な結果により、反射位相の可変性が～を超えることが実証されています。 310°、最大92.3%の吸収向上、ビームステアリング範囲の ±45° サイドローブレベルが-12 dB未満であり、180倍の計算加速を実現し、× 従来の全波最適化法と比較して、多項式回帰ではテスト時の平均二乗誤差が0.0688となり、R2 係数は0.994であり、H(z)フィッティングでは0.89 dBの平方根平均二乗誤差を達成した。安定性解析により、すべての極が単位円内にあることが確認された。本プロトコルは、プログラム可能なフォトニックメタサーフェスおよび6G以降のインテリジェント・テラヘルツ通信システムの設計に向けた、効率的で再現可能な経路を提供する。
概要
第6世代（6G）無線システムへの急速な進化により、1 Tbpsを超える超高速データレート、インテリジェントセンシング、ホログラフィックビームフォーミング、および適応型波面エンジニアリングを実現するためのテラヘルツ（THz）および光周波数帯の探索が加速しています1,2,3。テラヘルツ帯（0.1-30 THz）は広大な帯域幅を提供しますが、深刻な自由空間伝搬損失（1 THzで約20 dB·km-1）や、557 GHz、752 GHz、988 GHz、1.13 THzにおける水蒸気による大気分子吸収、さらに大気中の粒子や雨による遮蔽に対する極めて高い感度という課題があります4。
再構成可能な知的表面（RIS）は、プログラム可能な反射、吸収、位相制御、および波の集光を通じて、電磁波の波面を動的に操作できる有効な技術として登場しました5,6,7。これを光周波数領域に拡張することで、光の伝搬に対して前例のない制御が可能となり、LiDAR、自由空間光通信、ホログラフィックディスプレイ、適応型ビームステアリングなどのアプリケーションを実現します8。
金属共振器に基づく従来の再構成可能なインテリジェント表面構造は、テラヘルツおよび光周波数帯において、ドルーデ導電性によるオーム損失の増加、プラズモン減衰、100 nm以下の微細構造における作製上の制約、および固定された誘電特性の調整可能性の制限により、大きな制約に直面しています9,10,11。ナノ材料ベースのメタ表面、特に単層カーボンナノチューブ（SWCNTs）を用いたものは、その優れた量子電磁気的特性により、有望な代替案となります。具体的には、直径0.7-3 nmというナノスケールの寸法、化学ポテンシャルによって制御可能な量子導電性、100,000 cm2·(V·s-1)を超える極めて高いキャリア移動度、約3000 W·(m·K-1)の高い熱伝導率、および可視・近赤外領域における強いエキシトン共鳴などが挙げられます12,13,14,15。SWCNTsは、エキシトン-光子結合を通じてテラヘルツおよびフォトニック領域で強い電磁相互作用を示し、プログラム可能な電磁応答を持つコンパクトな光共振器を実現します16。これらの特性により、SWCNTsは、コンパクトなサイズ、調整可能な応答、および効率的な波操作が重要な設計要件となる光再構成可能インテリジェント表面の魅力的な候補となります。
これらの利点がある一方で、SWCNTベースの光再構成可能知的反射面（RIS）の正確なモデリングは依然として困難である。その要因として、非平衡グリーン関数法を必要とする量子輸送効果、100 GHz以上の周波数で支配的となる速度論的インダクタンス、最大0.4 eVの結合エネルギーを持つ多共鳴エキシトン挙動、電子-フォノン結合に起因する非線形分散現象、カイラル指数(n, m)の指定を必要とするカイラリティ依存の光学特性、および自己整合的な熱モデリングを必要とする温度依存の導電性が挙げられる。17,18,19,20従来の等価回路アプローチや、有理近似およびベクトルフィッティングなどの従来の電磁気フィッティング手法では、0.5～30 THzという広範な周波数領域におけるこれらの複雑な相互作用を正確に表現できないことが多く、60× 帯域幅21,22,23これらの制限は、テラヘルツ再構成可能インテリジェント表面の設計において、モデリング精度の低下や計算複雑性の増大を招く可能性があります。そのため、物理的な解釈可能性を維持しつつ、広帯域の電磁気学的挙動を正確に捉えることができる効率的なモデリング手法が必要とされています。
本プロトコルでは、Kubo形式を用いた量子伝導度解析、カイラリティ依存のバンドギャップ計算、ノイズ低減と応答平滑化のための多項式回帰、極零点安定性と物理的解釈性のためのZ変換伝達関数解析、およびSWCNTベースの単位セルの全波電磁シミュレーションを独自に組み合わせたハイブリッドモデリングフレームワークを紹介します。このワークフローは、テラヘルツ周波数領域におけるSWCNTベースの光再構成可能インテリジェント表面を構築、解析、および最適化するためのステップバイステップの手順を提供します。本手法は、高度な電磁表面、ナノフォトニクス、および次世代無線通信システムの研究に従事する研究者やエンジニア向けに設計されています。本プロトコルにより、ユーザーはビームステアリングおよび電磁波制御アプリケーションに向けた、安定し、かつ物理的に解釈可能なSWCNTベースの光再構成可能インテリジェント表面のモデルを生成することができます。本プロトコルでは、電磁シミュレーションの概念に精通していることを前提としていますが、この分野に初めて取り組む研究者が再現できるよう、詳細な手順を提示しています。
プロトコル
1. SWCNT光学RISユニットセルの設計
- SWCNTカイラリティの選択
- 量子導電性分析に基づき、SWCNTのカイラリティを(10,5)に選択する。
- 以下に示す式を用いて、ナノチューブの直径を算出します。
ここでaはcc炭素-炭素結合長は0.142 nmである。(10,5)のカイラリティにより、直径0.60 nm、バンドギャップ1.762 eVとなり、テラヘルツ動作に最適である。 - 以下を用いてカイラル角を算出する
- 単位格子の幾何学的定義
- 以下のパラメータを用いて、石英基板上に交差型SWCNTナノストリップ共振器を設計してください。基板誘電率 3.8、基板厚 50 nm、単位格子の寸法 1.5 × 1.5 µm2SWCNTストリップ幅50 nm、SWCNTストリップ長70 nm、およびストリップ間隔120 nm。これらの寸法により、0.5〜30 THz帯全域でサブ波長動作が保証されます。
- CST Microwave Studio シミュレーション設定
- CST Microwave Studioを起動し、「Microwave and RF frequency domain solver」を使用して新規プロジェクトを作成します。
- 定義されたパラメータを用いて単位格子の形状を構築します。x方向およびy方向の単位格子の境界条件に周期境界を適用します。平面波励起のため、z方向に沿ってフローケポートを定義します。
- 周波数スイープ範囲を0.5 THzから30 THzに設定し、ステップサイズを0.05 THzとして590個の周波数ポイントを生成します。時間領域ソルバーを適応メッシュリファインメントで構成し、収束ターゲットを-40 dBに設定します。
- 全波電磁界シミュレーションを実行する。複素S₁反射係数、S₂₁透過係数、反射位相、および吸収スペクトルをCSVファイルとしてエクスポートする。
2. 久保形式を用いた量子伝導度モデリング
- 物理定数の初期化
- 物理定数を初期化します：電気素量 e=1.602×10⁻19C, プランク定数 $\hbar=1.0546\times10^{-34}$ J·s34J·s、ボルツマン定数kB=1.3806×10⁻23J·K-1、温度 T=30K、フェルミ速度 vF=8×105m·s-1、および自由空間インピーダンスη₀=37Ω。
- バンド内導電率の算出
- 以下の式を用いて、各周波数点におけるバンド内導電率を算出します。
ここで、τ=0.5psはキャリア緩和時間であり、EはF フェルミエネルギー（0.1eVから0.4eVまで変化させた）。
- 以下の式を用いて、各周波数点におけるバンド内導電率を算出します。
- バンド間導電率の計算
- 以下の式を用いて、バンド間導電率を算出します。
- 全表面電導率をσとして算出する。合計(ω) = σ細胞内(ω) + σ間(ω).
- 以下の式を用いて、バンド間導電率を算出します。
- 光学応答計算
- 石英（εr=3.8、厚さ=500nm)をφを用いてsub=2πf(nsub)dsub/c（ここでnは）sub=.
- 導電率から反射係数を算出する方法
- 反射率を R(ω)=|S₁₁(ω)| として算出する。2吸収率を A(ω)=1-R(ω)-T(ω) として算出する。
- 石英（εr=3.8、厚さ=500nm)をφを用いてsub=2πf(nsub)dsub/c（ここでnは）sub=.
3. CSTデータの生成と前処理
- エミュレートされたCSTデータの生成
- 量子S行列に現実的な共鳴を加えることで、エミュレートされたCSTフルウェーブシミュレーション結果を生成する11.
- 5つの共鳴を定義する：2.8 THzにおけるE₁エキシトン（振幅-15.2 dB、位相 45°), E22 5.6 THzにおけるエキシトン（-2.8 dB, -120°）、12.4 THzにおけるプラズモン共鳴（-8.5 dB、 60°）、18.9 THzにおけるキャビティモード（-12.1 dB、-30°）、および24.7 THzにおけるフォノン支援共鳴（-6.8 dB、 15°).
- CSTの数値変動をシミュレートするため、ノイズレベル0.03（SNR=30.5dB）の数値ノイズを加える。
- データの前処理
- 以下を用いて周波数軸を[0,1]に正規化します
- S₁の実部と虚部、およびS₂₁の大きさを抽出します。
- ランダムパーミュテーションを用いて、データセットをトレーニング用（80%）、検証用（10%）、およびテスト用（10%）に分割する。入力を平均0、分散1となるように正規化する。
- 以下を用いて周波数軸を[0,1]に正規化します
4. データ平滑化のための多項式回帰
- 8を適合させるth polyfitを用いてS₁の実部に8次多項式を当てはめます：P_real = polyfit(X_train, Y_train(:,1), 8)。8次多項式でフィッティングします。th 虚数部に対する多項式の次数：P_imag = polyfit(X_train, Y_train(:,2), 8)。8次で近似します。th S₂₁の振幅に対する多項式次数を決定する：P_mag = polyfit(X_train, Y_train(:,3), 8)。全周波数範囲に対してpolyvalを用いて平滑化された予測値を生成する。平滑化されたS₁を S₁_ann = S1_real_ann + i × S1_imag_ann として算出する。
5. Z領域伝達関数のフィッティング
- サンプリング周波数Fを設定するs=60THz（ナイキスト基準に基づく最大周波数の2倍）。z=eを用いて周波数をZ領域にマッピングする。{i2πf/Fs}.
- Z領域における離散伝達関数は、次のように定義されます。
ここで、分子の次数はn=6、分母の次数はm=7である。 - 最小二乗最適化のために、線形方程式 A × x = B を設定します。ここで、A には分子および分母の項が含まれ、B には ANN で予測された S₁ が含まれます。数値的安定性を確保するため、QR分解を用いて解を求めます。
![QR分解の公式、[Q,R]=qr(A,0)、係数の算出、行列代数における方程式。](/files/ftp_upload/70498/70498eq8.jpg "figure-protocol-9")
- 極・零点および安定性解析
- 分母多項式を解くことで極を抽出します。|pole| ≧ 1 の極を poles(p) = poles(p) / (|poles(p)| + 0.1)を用いて単位円内に投影し、安定性を確保します。
- 安定化した極から分母を再構成する。多項式の分子を解くことで零点を抽出する。周波数グリッド上でH(z)を評価し、二乗平均平方根誤差を算出する。
6. ビームステアリングの最適化
- 線形位相アレイの正規化アレイ係数は次式で与えられます。
ここで、私はn = 1（一様励起の場合）、k = 2π/λ、d = λ/2 の間隔、βn は漸進的な位相シフトである。±の目標ステアリング角に対して45°, ±30°および 0°、隣接する素子間の位相差を次のように算出する。.EQUATION- 遠方場パターンを算出し、サイドローブレベルおよび半電力ビーム幅を評価する。
結果
SWCNTカイラリティの選択
記載されたプロトコルは、0.5-30 THz帯域で動作する(10,5) SWCNT光学RISユニットセルに対して実施されました。代表的な結果は、正確な電磁波特性評価およびビームステアリングの最適化における、ハイブリッド多項式-Zモデリング手法の有効性を示しています。
量子導電性解析
Kubo形式による解析から、(10,5) SWCNTはテラヘルツ帯域全体で虚数（誘導性）成分が支配的な複雑な表面導電性を示すことが明らかになりました。5 THz、E_F = 0.2 eVにおいて、導電率の実部は1.19 × 10⁻3 S·m-1であり、虚部は3.54 × 10⁻3 S·m-1であり、これは約71°の位相角に相当します。実部はω⁻1依存性に従って周波数とともに減少し、これはDrude的なバンド内輸送と一致しています。10 THz以上ではバンド間遷移が顕著になり、さらなる吸収チャネルに寄与します。
単位セルの電磁応答
全波電磁シミュレーションにより、単位セルの応答において5つの異なる共鳴が特定されました。2.8 THzのE₁エキシトンは-15.2 dBの反射係数を示しました。5.6 THzのE₂エキシトンは|S₁|が-2.8 dBとなり、電力反射率9.5%に相当する最強の反射を達成しました。12.4 THzのプラズモン共鳴では|S₁₁|が-8.5 dB、18.9 THzの共振器モードでは-12.1 dB、そして24.7 THzのフォノン支援共鳴では-6.8 dBとなりました。透過係数S₂₁は、共鳴周波数において深いノッチを持つ相補的な挙動を示しました。
量子導電率解析
ゲート電圧の変調に対応する0.1 eVから0.4 eVの範囲のフェルミエネルギーに対して、反射位相の特性評価を行いました。2.8 THz (E₁ エキシトン)において、反射位相はE_F = 0.1 eVでの-178°からE_F = 0.4 eVでの+132°まで変化し、310°の連続的な位相可変性を示しました。この極めて高い位相可変性は、従来の金属製RIS（通常は< 90°）を上回り、ビームステアリングアプリケーションにおける完全な360°のカバーを可能にします。位相応答は共鳴周波数の付近で急速な変化を示し、群遅延は-50 psから+80 psの範囲となりました。
多項式回帰の性能
8次多項式回帰により、生のCSTデータと比較して大幅なノイズ低減が達成されました。テスト平方根平均二乗誤差は0.068であり、数値的な変動を排除しつつ、フルウェーブシミュレーションと非常によく一致していることを示しています。R2係数は0.94であり、データの分散の9.4%が多項式モデルによって捉えられていることを示しています。この多項式は、CSTシミュレーションによる高周波の数値ノイズを除去しながら、5つの共振特性すべてを効果的に保持しました。
伝達関数のフィッティング
分子次数6、分母次数7のH(z)伝達関数を、多項式平滑化を行ったS₁応答に対して正常にフィッティングさせました。QR分解を用いた最小二乗最適化により、安定した解に収束しました。H(z)は、0.5-30 THzの全帯域にわたって0.89 dBの二乗平均平方根誤差を達成しました。多項式前処理なしで直接H(z)フィッティングを行った場合は、不安定な極が生じ、誤差も大きくなりました。直接フィッティングと比較して6.4 dBの改善が見られました。
極・零点安定性解析
安定性解析の結果、フィッティングされた伝達関数の7つの極すべてがZ平面の単位円内に存在することが明らかになりました。単位円に最も近い極の大きさは0.947であり、0.053の安定余裕が得られました。極の位置は5つの共振周波数に対応しており、2.8 THz、5.6 THz、12.4 THz、18.9 THz、および24.7 THzに対応する角度で単位円付近に極が位置していました。零点位置は最小位相および非最小位相の両方の特性を示し、3つの零点が単位円内部に、3つの零点が単位円外部に存在していました。また、すべての周波数において|H(z)| ≤ 1.02であり、システムが受動的であることが確認されました。
ビームステアリング性能
最適化された16素子の線形アレイは、±45°の範囲にわたるビームステアリングに成功した。目標角度-45°において、達成されたステアリング角度は-4.2°であり、誤差は0.8°、サイドローブレベルは-12.4 dBであった。目標-30°に対しては、達成角度-29.7°、サイドローブレベル-14.2 dBであった。0°のブロードサイドステアリングでは、指向性は20.1 dBiに達し、サイドローブレベルは-18.6 dBであった。目標+30°に対しては、達成角度+29.5°、サイドローブレベル-14.0 dBであった。目標+45°に対しては、達成角度+43.8°で誤差は1.2°、サイドローブレベルは-12.1 dBであった。半値幅は、ブロードサイドでの7.8°から最大ステアリング角度での8.5°の範囲であった。ステアリング精度は、すべての目標に対して1.2°以内であった。
量子伝導度モデリング後の光学的応答
化学ポテンシャルを変化させた際の吸収スペクトルでは、E_F = 0.2 eVにおいて5.6 THz（E₂エキシトン）で92.3%のピーク吸収が示された。2.8 THzのE₁エキシトンでは67%の吸収を達成したが、より高周波のモードでは、状態密度の減少により吸収が段階的に低下した。吸収はE_Fを変化させることで調整可能であり、E_Fを0.1 eVから0.4 eVに増加させると、パウリブロッキングにより吸収ピークが約0.3 THzブルーシフトし、ピーク吸収が15-20%減少した。
計算パフォーマンス
ハイブリッド多項式-Zフレームワークにより、遺伝的アルゴリズムを用いた従来のフルウェーブ最適化では48.75時間かかっていた総計算時間が、提案手法では47.3分に短縮され、単一の最適化において62倍の高速化を実現した。10回の設計イテレーションを伴うパラメータスイープでは、従来手法と比較して180倍の高速化係数に達した。多項式回帰自体の学習に要した時間はわずか0.5秒であり、学習後の光学応答の予測は0.02秒で完了した。
全体として、ハイブリッド多項式-Zフレームワークは、0.5–30 THz帯域におけるSWCNT光学RISの電磁応答を適切にモデル化することに成功しました。この手法により、正確な伝達関数のフィッティング、安定した極・零点挙動、調整可能な位相応答、効果的なビームステアリングが実証され、従来の最適化手法と比較して計算要件が大幅に削減されました。
図1SWCNTベースの光RISに向けたハイブリッド多項式-Zモデリングフレームワークのワークフロー概略図。 このワークフローには、SWCNT電導度の算出、全波電磁界シミュレーション、多項式平滑化、Z領域伝達関数のフィッティング、安定性解析、およびビームステアリングの最適化が含まれます。 こちらの図の拡大版を表示するには、ここをクリックしてください。
図2SWCNT光学RISユニットセルの幾何学的構造。(A) 交差したSWCNTナノストリップ共鳴器の上面図。 (B) 石英基板構造の側面図。 (C) Floquetポート励起および周期境界条件を示す三次元透視図。 この図の拡大版を表示するには、ここをクリックしてください。
図3久保形式を用いて計算した(10,5)単層カーボンナノチューブ(SWCNT)の量子伝導度。(A) 周波数の関数としての表面導電率の実数成分および虚数成分。 (B) シミュレーションした周波数範囲における導電率位相角。 この図の拡大版を表示するには、ここをクリックしてください。
図4フルウェーブシミュレーションから得られたSWCNT単位格子の電磁応答。(A) 反射係数 (|S1|) (B) 調査した周波数範囲における伝達係数（|S21|）。 この図の拡大版を表示するには、ここをクリックしてください。
図 5異なるフェルミ準位におけるSWCNT光RISの反射位相応答。 EF = 0.1 eV、0.2 eV、0.3 eV、および 0.4 eV における位相応答を示す。 この図の拡大版を表示するには、ここをクリックしてください。
図6Z平面における近似伝達関数の極零点配置図。 安定性評価のため、極と零点の位置を単位円と共に示しています。 こちらの図の拡大版を表示するには、ここをクリックしてください。
図7SWCNT光学RISの遠方場ビームステアリング性能。 ターゲット指向角度における放射パターンを以下に示す。 (A) -45°, (B) -30°, (C) 0°, (D) +30°、および (E) +45°. (F) 目標操舵角と到達操舵角の比較。 この図の拡大版を表示するには、ここをクリックしてください。
図 8異なるフェルミエネルギー値におけるSWCNT光RISの光吸収スペクトル。 EF = 0.1 eV、0.2 eV、0.3 eV、および0.4 eVにおける吸収応答を示す。 この図の拡大版を表示するには、ここをクリックしてください。
| モデリング手法 | RMSE (dB) | 安定性の保証 | パッシビティ確認 | 計算時間（デザインあたり） | 物理的解釈可能性 | 多共鳴キャプチャー |
| 等価回路 (RLC) | 5.0 - 10.0 | はい | 限定的 | < 1秒 | 高 | 不良な |
| 有理近似（ベクトルフィッティング） | 2.0 - 5.0 | なし（30% 不安定） | いいえ | 5～10秒 | 中程度の | 中程度 |
| 全波電磁界解析（CSTのみ） | 参考文献 (0) | 該当なし | 該当事項なし | 6～8時間 | 高 | 非常に優れた |
| ANNのみ（Z変換なし） | 0.9 - 1.5 | いいえ | いいえ | 0.02秒 | 低 | 良好 |
| ハイブリッド・ポリノミアルZ（本研究） | 0.89 | はい（すべての極｜ < 1) | はい (|H(z)| ≤ 1.02) | 0.5秒（学習）+ 0.02秒（予測） | 高 | 非常に優れた |
表1：RMSE、安定性、受動性、計算時間、物理的解釈可能性、および多共振応答能力に基づくモデリング手法の比較。
| 目標角度 (°) | 達成角度 (°) | 角度誤差 (°) | サイドローブレベル (dB) | 半値幅 (°) | 指向性 (dBi) | メインローブ幅 (ヌル間, °) |
| -45 | -44.2 | 0.8 | -12.4 | 8.5 | 18.7 | 17.2 |
| -30 | -29.7 | 0.3 | -14.2 | 8.2 | 19.2 | 16.8 |
| 0 | 0 | 0 | -18.6 | 7.8 | 20.1 | 15.9 |
| 30 | 29.5 | 0.5 | -14 | 8.1 | 19.1 | 16.7 |
| 45 | 43.8 | 1.2 | -12.1 | 8.4 | 18.5 | 17 |
表2：ターゲット操舵角-45°、-30°、0°、+30°、および+45°におけるビームステアリング性能指標。パラメータには、達成角度、角度誤差、サイドローブレベル、半電力ビーム幅、指向性、およびメインローブ幅が含まれる。
| 手法 | セットアップ時間 | 設計イテレーションあたりの計算時間 | 合計時間（単一最適化） | 合計時間（10回のパラメータスイープ） | 高速化係数（従来法との比較） |
| Full-Wave EM (CSTのみ) | 45 min | 6.2 hours | 7.0 hours | 620 hours (25.8 days) | 1× (ベースライン) |
| Full-Wave EM + 遺伝的アルゴリズム | 45 min | 48 hours | 48.75 hours | 4,875 hours (203 days) | 0.14× |
| ANNのみ (Z変換なし) | 45 min (CST) + 10 min (学習) | 0.02 sec | 5.2 min | 58.3 min | 640× (100スイープ時) |
| ベクトルフィッティング (直接H(z)) | 45 min (CST) | 0.5 sec | 45.5 min | 50 min | 74× (10スイープ時) |
| ハイブリッド多項式-Z法 (本研究) | 45 min (CST) + 0.5 sec (学習) | 0.02 sec (予測) | 47.3 min | 47.5 min | 180× (10スイープ時) |
表3：評価したモデリング手法の計算パフォーマンス比較。指標には、セットアップ時間、設計イテレーションあたりの計算時間、総最適化時間、10回のパラメータスイープにかかる総時間、および相対的な加速係数が含まれる。
ディスカッション
プロトコルの重要なステップ においては、確実に実施するために細心の注意が必要です。第一に、バンドギャップと光応答はカイラリティに強く依存するため、SWCNTのカイラル指数の正確な選択が不可欠です。本プロトコルで指定している(10,5)カイラリティは、テラヘルツ動作に最適な1.762 eVのバンドギャップを提供しますが、異なる周波数帯を対象とする場合は、バンドギャップ式 E_g = 2ħv_F/d = (2 × 1.0546×10⁻34 × 8×105)/(d) eV を用いて対応するカイラリティを算出する必要があります。例えば、10 THz動作（光子エネルギー 41 meV）を対象とするには、より直径の大きいナノチューブ、または金属的挙動に近いより低いバンドギャップが必要です。第二に、全波シミュレーションのメッシュ解像度は、50 nmのSWCNTストリップ幅を分解するのに十分な細かさである必要があります。特に波長が10 µmとなる最大周波数の30 THzにおいて重要です。波長あたり最低20セルのメッシュ密度が推奨されており、これは30 THzにおいて0.5 µmのセルサイズに相当しますが、正確な結果を得るためには50 nmストリップ周辺の局所的なリファインメント（比率 200:1）が必要です。第三に、多項式回帰の次数の選択では、バイアスとバリアンスのバランスを取る必要があります。次数8は赤池情報量基準（AIC）の最小化に基づいて選択されました。より低い次数（4-6）では共鳴の適合不足（アンダーフィット）となり、より高い次数（10-12）では数値ノイズへの過適合（オーバーフィット）が生じます。ユーザーは、自身の特定の単位セル設計における最適な次数を決定するために、クロスバリデーションを行う必要があります。
一般的な実装上の課題は、修正とトラブルシューティング で対処可能です。多項式回帰によって周波数帯域の端付近にリンギングアーチファクト（ルンゲ現象）が生じる場合は、標準的な多項式フィッティングをチェビシェフ多項式近似またはスプライン補間に置き換えてください。H(z)伝達関数が受動性を侵害し、|H(z)|が1.05を超える場合は、分子と分母の次数をn = 4、m = 5に下げるか、留数摂動などの受動性強制手法を適用してください。ビームステアリングの最適化においてサイドローブレベルが高くなり-10 dBを超える場合は、アレイサイズをN = 16からN = 32素子に増やすか、ハミング窓またはカイザー窓を用いて振幅テーパリングを適用し、ビーム幅の増大と引き換えにサイドローブを10-20 dB低減させてください。量子導電率の計算が超高周波数（> 25 THz）で収束しない場合は、バンド間導電率項が支配的になります。パウリブロッキングによってバンド間遷移が抑制されるE_F > 0.3 eVにおいては、バンド内寄与のみを用いて簡略化してください。
本手法の限界 このプロトコルの適用前に検討すべき事項がある。まず、量子伝導モデルは、純粋な単層カーボンナノチューブ（SWCNT）における弾道輸送を前提としており、欠陥、不純物、または高密度アレイ（間隔）におけるチューブ間の相互作用による散乱を完全には考慮していない。 < （50 nm）。カイラリティ分布が10%より広い、または欠陥密度が100 nmあたり1個を超える実際の試料では、理想的な(10,5)の応答が実験測定値と最大30%異なる可能性があります。第二に、多項式回帰は単一のユニットセル形状（固定長700 nm、間隔120 nm）のデータで学習されており、再学習なしでは大幅に異なる設計に一般化できない可能性があります。転移学習の手法を用いることで、新しい形状に必要な学習データを削減できる可能性があります。第三に、伝達関数のフィッティングは線形時不変挙動を前提としていますが、1 kWを超える高強度の光励起下では、この前提が成立しない可能性があります。·(cm2)-1 飽和吸収などの非線形効果（特性フルエンス〜10 µJ·(cm2)-1 （SWCNTの場合）またはカー非線形性（n₂ ~ 10⁻12 cm2·W) が顕著になります。第四に、本プロトコルでは室温を超える熱的影響は考慮されていませんが、SWCNTアレイでは顕著な発熱が生じる可能性があります（ΔT > （100 K）の連続波動作下での光周波数においてであり、これはキャリア移動度および緩和時間に影響を及ぼします。
本手法の意義 既存の代替手法と比較して、その利点は相当なものである。集中定数RLCネットワークを用いた従来の等価回路モデルでは、SWCNTベースのメタ表面の多共振挙動を捉えることができず、二乗平均平方根誤差は通常5〜10 dBとなる。SWCNTベースの光学再構成可能インテリジェント表面を研究するための代替アプローチには、直接的なフルウェーブ電磁最適化、等価回路モデリング、ベクトルフィッティング、および物理ベースの量子輸送シミュレーションがある。しかし、これらのアプローチでは通常、計算コスト、精度、および物理的な解釈可能性の間にトレードオフが存在する。多項式前処理を行わないベクトルフィッティング法では、高次モデル（n > 8). ハイブリッド多項式-Zアプローチは、多項式回帰のパターン学習能力と伝達関数の物理的な解釈可能性を独自に組み合わせており、安定性を保証しつつ、0.9 dB未満の平均二乗誤差根（RMSE）を達成している。180× 計算の高速化により、従来のフルウェーブ最適化のみでは不可能であった設計空間の探索が可能になります。例えば、5つの幾何学的パラメータ（長さ、幅、間隔、基板厚、カイラリティ）について、それぞれ10個の値を用いて最適化を行う場合、設計の組み合わせは100,000通りになります。従来の最適化では、 > 計算時間は10年を要しますが、ハイブリッド法では約20日で完了します。
潜在的な応用例 本プロトコルの応用範囲は、特定のデモンストレーションに留まらず、いくつかの研究領域に及びます。6G無線通信においては、本手法を用いて、テラヘルツ波ビームステアリング（0.1～10 THz）、無線電力伝送のためのビームフォーカシング、およびモード分割多重のための軌道角運動量生成に向けた再構成可能な知的表面（RIS）を設計することが可能です。自動運転車用LiDARシステムにおいては、光位相の可変性により、0.1の精度で非機械的なビームスキャニングを実現できます。° 分解能およびマイクロ秒単位のスイッチング時間を実現しており、ミリ秒単位の応答速度で寿命が限られている機械的システムとは対照的である。ホログラフィックディスプレイにおいては、サブ波長単位セルにより振幅および位相制御が可能となり、4K解像度および60 Hzのリフレッシュレートで三次元画像の投影が行われる。量子通信においては、極低温下における単層カーボンナノチューブ（SWCNT）のエキシトン共鳴により、単一光子の操作、自発四波混合によるもつれ光子対の生成、および1 nsを超えるコヒーレンス時間を持つ量子メモリへの応用が可能となる。
今後の展開 本プロトコルに、動的な環境における適応的なビーム最適化のための強化学習を組み込むことが可能であり、そこではRISが無線チャネルとの相互作用を通じて最適な位相構成を学習する。化学気相成長法による成長および電子ビームリソグラフィーによるパターニングを用いて作製したSWCNTメタ表面による実験的検証は、モデルの精緻化に不可欠なフィードバックを提供するであろう。二層および多層カーボンナノチューブへの拡張は、高出力アプリケーションにおける帯域幅と熱安定性を向上させる可能性がある。光集積回路との統合により、オンチップ制御電子回路を備えたチップスケールの光再構成可能インテリジェント表面が実現可能となる。最後に、密度汎関数理論に基づく完全なバンド構造計算を組み込むことで、本研究で検討した(10,5)以外のカイラリティに対する精度を向上させることができる。
開示事項
著者は、利益相反がないことを宣言します。
謝辞
著者は、貴重な科学的および技術的支援をいただいたイラクのBaghdad QuarterにあるInternational Applied and Theoretical Research Center (IATRC) に心より感謝いたします。本研究は、公的、商業的、または非営利セクターの資金提供機関から特定の助成金を受けていません。計算リソースはAl-Bayan Universityより提供されました。
材料
| 名前 | 会社 | カタログ番号 | コメント |
|---|---|---|---|
| CST Microwave Studio | Dassault Systèmesèmes | 翻訳対象のテキストが提供されていません。翻訳する英文をご提示ください。 | 2024年版、周波数領域ソルバー |
| MATLAB | MathWorks | 提供されたソーステキストがありません。翻訳対象のテキストを入力してください。 | Version R2014a以降 |
| 石英基板 | 大学ウェハー | 4526 | 厚さ500 nm、 ε_r = 3.8 |
| SWCNT (10,5) カイラリティ | NanoIntegris | SWCNT-105 | 直径0.60 nmの >90% 半導体性 |
| パーソナルコンピューター | 翻訳対象のテキストが提供されておりません。翻訳する英文をご提示ください。 | 該当事項なし | 最低32 GB RAM、8 CPUコア |
参考文献
- Xiao, M., et al. Millimeter wave communications for future mobile networks. IEEE J Sel Areas Commun. 35, 1909-1935 (2017).
- Kumar, A., et al. RIS-assisted terahertz communications for 6G networks: A comprehensive overview. IEEE Access. , (2025).
- Liaskos, C., et al. A new wireless communication paradigm through software-controlled metasurfaces. IEEE Commun Mag. 56, 162-169 (2018).
- Wu, Q., et al. Intelligent reflecting surface-aided wireless communications: A tutorial. IEEE Trans Commun. 69, 3313-3351 (2021).
- Rafique, A., et al. Reconfigurable intelligent surfaces: Interplay of multi cell and surface-level design and performance under quantifiable benchmarks. IEEE Open J Commun Soc. 4, 1583-1599 (2023).
- Mayaram, K., et al. Computer-aided circuit analysis tools for RFIC simulation: algorithms, features, and limitations. IEEE Trans Circuits Syst II. 47, 274-286 (2000).
- Yesilyurt, O., Turhan-Sayan, G. Metasurface lens for ultra-wideband planar antenna. IEEE Trans Antennas Propag. 68, 719-726 (2019).
- Gustavsen, B., Semlyen, A. Rational approximation of frequency domain responses by vector fitting. IEEE Trans Power Deliv. 14, 1052-1061 (2002).
- Zhang, Y., et al. Z-transform-based FDD implementations of biaxial anisotropy for radar target scattering problems. Remote Sens. 14, 2397(2022).
- Hall, S. H., Heck, H. L. Advanced Signal Integrity for High-Speed Digital Designs. , John Wiley & Sons. (2011).