このコンテンツを表示するには、JoVEへの購読が必要です。 または、無料トライアルをお申し込みください。

方法論記事

テラヘルツ波ビームステアリングに向けたカーボンナノチューブベース再構成可能知能表面のモデリングのためのハイブリッドANN-Z法

53 回視聴

⸱

DOI:

10.3791/70498

⸱

2026年8月28日

 ,  , 
1Technical Engineering College, Al-Bayan University, 2College of Computer Science and Mathematics, Tikrit University, 3Avicenna E-Learning Center, Tikrit University, 4Department of Automation and Artificial Intelligence Engineering, College of Information Engineering at Al-Nahrain University, 5International Applied and Theoretical Research Center (IATRC)

この記事について

サマリー

本プロトコルでは、6G無線通信アプリケーション向けにテラヘルツ帯（0.5～30 THz）で動作する単層カーボンナノチューブベースの光再構成可能知的反射面の正確な電磁界モデリングを実現するため、ハイブリッド人工ニューラルネットワークとZ変換を用いた手法を提示し、180× 反射位相調整能の向上による計算加速 310° およびビームステアリング範囲は ±45°.

要約

単層カーボンナノチューブに基づく光再構成可能インテリジェント表面は、次世代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ユニットセルの設計

  1. SWCNTカイラリティの選択
    1. 量子導電性分析に基づき、SWCNTのカイラリティを(10,5)に選択する。
    2. 以下に示す式を用いて、ナノチューブの直径を算出します。
      六方格子定数、d=acc*√(n²+nm+m²)/π、構造方程式図。
      ここでaはcc炭素-炭素結合長は0.142 nmである。(10,5)のカイラリティにより、直径0.60 nm、バンドギャップ1.762 eVとなり、テラヘルツ動作に最適である。
    3. 以下を用いてカイラル角を算出する 結晶学の角度公式 θ = tan⁻¹(√3 × m/(2n + m)) = 23.4°、数学方程式。
  2. 単位格子の幾何学的定義
    1. 以下のパラメータを用いて、石英基板上に交差型SWCNTナノストリップ共振器を設計してください。基板誘電率 3.8、基板厚 50 nm、単位格子の寸法 1.5 × 1.5 µm2SWCNTストリップ幅50 nm、SWCNTストリップ長70 nm、およびストリップ間隔120 nm。これらの寸法により、0.5〜30 THz帯全域でサブ波長動作が保証されます。
  3. CST Microwave Studio シミュレーション設定
    1. CST Microwave Studioを起動し、「Microwave and RF frequency domain solver」を使用して新規プロジェクトを作成します。
    2. 定義されたパラメータを用いて単位格子の形状を構築します。x方向およびy方向の単位格子の境界条件に周期境界を適用します。平面波励起のため、z方向に沿ってフローケポートを定義します。
    3. 周波数スイープ範囲を0.5 THzから30 THzに設定し、ステップサイズを0.05 THzとして590個の周波数ポイントを生成します。時間領域ソルバーを適応メッシュリファインメントで構成し、収束ターゲットを-40 dBに設定します。
    4. 全波電磁界シミュレーションを実行する。複素S₁反射係数、S₂₁透過係数、反射位相、および吸収スペクトルをCSVファイルとしてエクスポートする。

2. 久保形式を用いた量子伝導度モデリング

  1. 物理定数の初期化
    1. 物理定数を初期化します：電気素量 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Ω。
  2. バンド内導電率の算出
    1. 以下の式を用いて、各周波数点におけるバンド内導電率を算出します。
      静的平衡式、複素導電率方程式、物理学研究、数学的表記
      ここで、τ=0.5psはキャリア緩和時間であり、EはF フェルミエネルギー（0.1eVから0.4eVまで変化させた）。
  3. バンド間導電率の計算
    1. 以下の式を用いて、バンド間導電率を算出します。
      複雑な公式を用いた静的平衡の概念。物理学的解析のための数式を含む。
    2. 全表面電導率をσとして算出する。合計(ω) = σ細胞内(ω) + σ間(ω).
  4. 光学応答計算
    1. 石英（εr=3.8、厚さ=500nm)をφを用いてsub=2πf(nsub)dsub/c（ここでnは）sub=. 比誘電率の平方根、√εr、公式；電磁気学、材料特性。
    2. 導電率から反射係数を算出する方法
      分光分析結果、S1(ω)式、電磁波挙動の特性評価。
    3. 反射率を R(ω)=|S₁₁(ω)| として算出する。2吸収率を A(ω)=1-R(ω)-T(ω) として算出する。

3. CSTデータの生成と前処理

  1. エミュレートされたCSTデータの生成
    1. 量子S行列に現実的な共鳴を加えることで、エミュレートされたCSTフルウェーブシミュレーション結果を生成する11.
    2. 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°).
    3. CSTの数値変動をシミュレートするため、ノイズレベル0.03（SNR=30.5dB）の数値ノイズを加える。
  2. データの前処理
    1. 以下を用いて周波数軸を[0,1]に正規化します
      正規化周波数式 \(z_{f_{norm}}=\frac{f-f_{min}}{f_{max}-f_{min}}\) の数式。
    2. S₁の実部と虚部、およびS₂₁の大きさを抽出します。
    3. ランダムパーミュテーションを用いて、データセットをトレーニング用（80%）、検証用（10%）、およびテスト用（10%）に分割する。入力を平均0、分散1となるように正規化する。

4. データ平滑化のための多項式回帰

  1. 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領域伝達関数のフィッティング

  1. サンプリング周波数Fを設定するs=60THz（ナイキスト基準に基づく最大周波数の2倍）。z=eを用いて周波数をZ領域にマッピングする。{i2πf/Fs}.
  2. Z領域における離散伝達関数は、次のように定義されます。
    デジタル信号処理で用いられる、多項式で構成される離散時間伝達関数式 H(z)。
    ここで、分子の次数はn=6、分母の次数はm=7である。
  3. 最小二乗最適化のために、線形方程式 A × x = B を設定します。ここで、A には分子および分母の項が含まれ、B には ANN で予測された S₁ が含まれます。数値的安定性を確保するため、QR分解を用いて解を求めます。
    QR分解の公式、[Q,R]=qr(A,0)、係数の算出、行列代数における方程式。
  4. 極・零点および安定性解析
    1. 分母多項式を解くことで極を抽出します。|pole| ≧ 1 の極を poles(p) = poles(p) / (|poles(p)| + 0.1)を用いて単位円内に投影し、安定性を確保します。
    2. 安定化した極から分母を再構成する。多項式の分子を解くことで零点を抽出する。周波数グリッド上でH(z)を評価し、二乗平均平方根誤差を算出する。

6. ビームステアリングの最適化

  1. 線形位相アレイの正規化アレイ係数は次式で与えられます。
    ​アンテナ設計におけるアレイファクタ式；AF(θ)公式；工学における数学方程式。
    ここで、私はn = 1（一様励起の場合）、k = 2π/λ、d = λ/2 の間隔、βn は漸進的な位相シフトである。±の目標ステアリング角に対して45°, ±30°および 0°、隣接する素子間の位相差を次のように算出する。.EQUATION
    1. 遠方場パターンを算出し、サイドローブレベルおよび半電力ビーム幅を評価する。

アクセスが制限されています。このコンテンツを表示するにはログインするか、トライアルを開始してください。

結果

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以上ではバンド間遷移が顕著になり、さらなる...

アクセスが制限されています。このコンテンツを表示するにはログインするか、トライアルを開始してください。

ディスカッション

プロトコルの重要なステップ においては、確実に実施するために細心の注意が必要です。第一に、バンドギャップと光応答はカイラリティに強く依存するため、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のセルサイズに相当しま...

アクセスが制限されています。このコンテンツを表示するにはログインするか、トライアルを開始してください。

開示事項

著者は、利益相反がないことを宣言します。

謝辞

著者は、貴重な科学的および技術的支援をいただいたイラクのBaghdad QuarterにあるInternational Applied and Theoretical Research Center (IATRC) に心より感謝いたします。本研究は、公的、商業的、または非営利セクターの資金提供機関から特定の助成金を受けていません。計算リソースはAl-Bayan Universityより提供されました。

アクセスが制限されています。このコンテンツを表示するにはログインするか、トライアルを開始してください。

材料

この記事で使用された材料の一覧
名前会社カタログ番号コメント
CST Microwave StudioDassault Systèmesèmes翻訳対象のテキストが提供されていません。翻訳する英文をご提示ください。2024年版、周波数領域ソルバー
MATLABMathWorks提供されたソーステキストがありません。翻訳対象のテキストを入力してください。Version R2014a以降
石英基板大学ウェハー4526厚さ500 nm、 ε_r = 3.8
SWCNT (10,5) カイラリティNanoIntegrisSWCNT-105直径0.60 nmの >90% 半導体性
パーソナルコンピューター翻訳対象のテキストが提供されておりません。翻訳する英文をご提示ください。該当事項なし最低32 GB RAM、8 CPUコア

参考文献

  1. Xiao, M., et al. Millimeter wave communications for future mobile networks. IEEE J Sel Areas Commun. 35, 1909-1935 (2017).
  2. Kumar, A., et al. RIS-assisted terahertz communications for 6G networks: A comprehensive overview. IEEE Access. , (2025).
  3. Liaskos, C., et al. A new wireless communication paradigm through software-controlled metasurfaces. IEEE Commun Mag. 56, 162-169 (2018).
  4. Wu, Q., et al. Intelligent reflecting surface-aided wireless communications: A tutorial. IEEE Trans Commun. 69, 3313-3351 (2021).
  5. 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).
  6. 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).
  7. Yesilyurt, O., Turhan-Sayan, G. Metasurface lens for ultra-wideband planar antenna. IEEE Trans Antennas Propag. 68, 719-726 (2019).
  8. Gustavsen, B., Semlyen, A. Rational approximation of frequency domain responses by vector fitting. IEEE Trans Power Deliv. 14, 1052-1061 (2002).
  9. Zhang, Y., et al. Z-transform-based FDD implementations of biaxial anisotropy for radar target scattering problems. Remote Sens. 14, 2397(2022).
  10. Hall, S. H., Heck, H. L. Advanced Signal Integrity for High-Speed Digital Designs. , John Wiley & Sons. (2011).

アクセスが制限されています。このコンテンツを表示するにはログインするか、トライアルを開始してください。

再版と許可

タグ

カーボンナノチューブ表面量子導電率モデリング久保形式多項式回帰Z領域解析電磁界シミュレーションビームステアリング最適化フォトニックメタ表面
お問い合わせ
図書館に推薦する

プライバシー

利用規約

ポリシー

著作権 © 2026 MyJoVE Corporation. 無断転載を禁じます。