Méthode hybride ANN-Z pour la modélisation de surfaces intelligentes reconfigurables basées sur des nanotubes de carbone pour le pilotage de faisceau térahertz

Résumé

Ce protocole présente une méthode hybride combinant un réseau neuronal artificiel et une transformation en Z pour la modélisation électromagnétique précise de surfaces intelligentes reconfigurables optiques à base de nanotubes de carbone monoparois fonctionnant dans la bande térahertz (0,5-30 THz) pour des applications de communication sans fil 6G, permettant une accélération computationnelle de 180×, une possibilité de réglage de la phase de réflexion supérieure à 310° et une plage de guidage du faisceau de ±45°.

Résumé

Les surfaces intelligentes reconfigurables par voie optique à base de nanotubes de carbone à paroi unique offrent des solutions prometteuses pour le guidage de faisceaux térahertz et la manipulation d'ondes photoniques dans les futurs systèmes sans fil 6G. Toutefois, la modélisation précise de ces structures reste difficile en raison des effets de transport quantique, de l'inductance cinétique et du comportement excitonique multi-résonant sur de larges plages de fréquences. Ce protocole décrit un cadre computationnel hybride intégrant la modélisation de la conductivité quantique par le formalisme de Kubo, un lissage des données fondé sur une régression polynomiale et une analyse de fonction de transfert dans le domaine Z, afin de caractériser précisément les cellules unitaires de surfaces intelligentes reconfigurables optiquement à base de nanotubes de carbone à paroi unique. La méthode commence par la conception de résonateurs nano-bandes croisés en nanotubes de carbone à paroi unique sur un substrat de quartz, de chiralité (10,5) (diamètre 0,60 nm, bande interdite 1,762 eV), suivie d'une simulation électromagnétique en régime d'ondes complètes dans CST Microwave Studio sur la bande 0,5-30 THz. Un modèle de régression polynomiale d'ordre 8 traite les paramètres S extraits afin d'éliminer les fluctuations numériques et de prédire des réponses électromagnétiques lissées. Une fonction de transfert discrète H(z), d'ordre 6 au numérateur et d'ordre 7 au dénominateur, est ensuite ajustée par optimisation aux moindres carrés avec décomposition QR, permettant une analyse de stabilité par pôles et zéros ainsi que la vérification de la passivité. Le protocole intègre en outre un réglage de la conductivité quantique par modulation du potentiel chimique afin d'optimiser le guidage du faisceau. Des résultats représentatifs montrent une plage de réglage de la phase en réflexion supérieure à 310°, un gain d'absorption atteignant 92,3 %, une plage de guidage du faisceau de ±45° avec des niveaux de lobes secondaires inférieurs à -12 dB, ainsi qu'une accélération computationnelle de 180× par rapport aux méthodes conventionnelles d'optimisation par ondes complètes. La régression polynomiale a atteint une erreur quadratique moyenne de 0,0688 avec un coefficient R² de 0,994, tandis que l'ajustement de H(z) a donné une erreur quadratique moyenne de 0,89 dB. L'analyse de stabilité a confirmé que tous les pôles se situent à l'intérieur du cercle unité. Ce protocole fournit une voie efficace et reproductible pour concevoir des métasurfaces photoniques programmables et des systèmes de communication térahertz intelligents pour la 6G et au-delà.

Introduction

L'évolution rapide vers les systèmes sans fil de sixième génération (6G) a accéléré l'exploration des bandes de fréquences térahertz (THz) et optiques afin d'atteindre des débits de données ultra-élevés dépassant 1 Tbps, ainsi que des fonctionnalités telles que la détection intelligente, le formage de faisceau holographique et l'ingénierie adaptative du front d'onde^1,2,3. La bande térahertz (0,1-30 THz) offre une largeur de bande abondante, mais souffre d'une forte atténuation en espace libre (environ 20 dB·km^-1 à 1 THz), d'une absorption atmosphérique moléculaire due à la vapeur d'eau à 557 GHz, 752 GHz, 988 GHz et 1,13 THz, ainsi que d'une sensibilité extrême au brouillage causé par les particules atmosphériques et la pluie⁴.

Les surfaces intelligentes reconfigurables sont apparues comme une technologie permettant de manipuler dynamiquement les fronts d'onde électromagnétiques grâce à une réflexion programmable, une absorption, un contrôle de phase et une focalisation des ondes^5,6,7. Lorsqu'elles sont étendues aux fréquences optiques, ces surfaces offrent un contrôle sans précédent de la propagation de la lumière, permettant des applications telles que le LiDAR, les communications optiques en espace libre, les affichages holographiques et le guidage adaptatif de faisceaux⁸.

Les structures conventionnelles de surfaces intelligentes reconfigurables basées sur des résonateurs métalliques rencontrent des limitations importantes aux fréquences térahertz et optiques en raison des pertes ohmiques accrues dues à la conductivité de Drude, de l'amortissement plasmonique, des contraintes de fabrication liées à des dimensions inférieures à 100 nm, et de la faible possibilité de réglage des propriétés diélectriques fixes^9,10,11. Les métasurfaces à base de nanomatériaux, en particulier celles utilisant des nanotubes de carbone monoparois (SWCNT), offrent des alternatives prometteuses grâce à leurs propriétés quantiques et électromagnétiques exceptionnelles : dimensions nanométriques avec des diamètres de 0,7 à 3 nm, conductivité quantique ajustable contrôlée par le potentiel chimique, mobilité des porteurs extraordinaire dépassant 100 000 cm²·(V·s^-1), conductivité thermique élevée d’environ 3000 W·(m·K^-1), et résonances excitoniques fortes dans les régimes visible et proche infrarouge^12,13,14,15. Les SWCNT présentent des interactions électromagnétiques fortes dans les régimes térahertz et photoniques grâce au couplage exciton-photon, permettant la réalisation de résonateurs optiques compacts dotés de réponses électromagnétiques programmables¹⁶. Ces propriétés rendent les SWCNT particulièrement attractifs pour des surfaces intelligentes reconfigurables en optique, où la compacité, la réponse ajustable et une manipulation efficace des ondes constituent des critères de conception essentiels.

Malgré ces avantages, la modélisation précise des surfaces intelligentes reconfigurables optiques basées sur des nanotubes de carbone (SWCNT) reste difficile en raison des effets de transport quantique nécessitant des approches par fonction de Green hors équilibre, de l'inductance cinétique prédominante à des fréquences supérieures à 100 GHz, du comportement excitonique multi-résonant avec des énergies de liaison allant jusqu'à 0,4 eV, des phénomènes de dispersion non linéaire découlant du couplage électron-phonon, des propriétés optiques dépendant de la chiralité et exigeant la spécification des indices chiraux (n, m), ainsi que de la conductivité dépendant de la température, nécessitant une modélisation thermique auto-cohérente^17,18,19,20. Les approches traditionnelles par circuits équivalents et les méthodes conventionnelles d'ajustement électromagnétique, telles que l'approximation rationnelle et l'ajustement vectoriel, échouent souvent à représenter avec précision ces interactions complexes sur de larges plages de fréquences couvrant 0,5-30 THz, une bande passante de 60×^21,22,23. Ces limitations peuvent réduire la précision de la modélisation et augmenter la complexité computationnelle lors de la conception de surfaces intelligentes reconfigurables en térahertz, soulignant ainsi le besoin d'approches de modélisation efficaces capables de préserver l'interprétabilité physique tout en capturant avec exactitude le comportement électromagnétique large bande.

Ce protocole présente un cadre de modélisation hybride qui combine de manière unique l'analyse de la conductivité quantique à l'aide de la formalisation de Kubo, le calcul de la bande interdite dépendant de la chiralité, une régression polynomiale pour la réduction du bruit et la lissage de la réponse, une analyse de fonction de transfert par transformation en Z afin d'évaluer la stabilité des pôles et des zéros ainsi que l'interprétabilité physique, et une simulation électromagnétique en ondes complètes de cellules unitaires à base de nanotubes de carbone monoparois (SWCNT). Le flux de travail fournit une procédure étape par étape pour concevoir, analyser et optimiser des surfaces intelligentes reconfigurables optiquement à base de SWCNT dans la gamme de fréquences térahertz. La méthode s'adresse aux chercheurs et ingénieurs travaillant sur les surfaces électromagnétiques avancées, la nanophotonique et les systèmes de communication sans fil de nouvelle génération. Ce protocole permet aux utilisateurs de générer des modèles stables et physiquement interprétables de surfaces intelligentes reconfigurables optiquement à base de SWCNT, destinés à des applications de guidage de faisceau et de contrôle des ondes électromagnétiques. Le protocole suppose une connaissance des concepts de simulation électromagnétique, mais fournit des étapes détaillées permettant sa reproduction par des chercheurs novices dans le domaine.

Protocole

1. Conception de la cellule unitaire optique RIS en nanotubes de carbone SWCNT

1. Sélection de la chiralité des SWCNT
   1. Sélectionner la chiralité (10,5) des SWCNT en se basant sur l'analyse de la conductivité quantique.
   2. Calculer le diamètre du nanotube à l'aide de la formule suivante
      
      où a_cc=0,142 nm est la longueur de la liaison carbone-carbone. La chiralité (10,5) donne un diamètre de 0,60 nm et un gap de bande de 1,762 eV, optimal pour le fonctionnement en térahertz.
   3. Calculer l'angle chiral à l'aide de
2. Définition de la géométrie de la cellule unitaire
   1. Concevoir des résonateurs nano-bandes en croix à base de SWCNT sur un substrat de quartz avec les paramètres suivants : permittivité du substrat de 3,8, épaisseur du substrat de 500 nm, dimensions de la cellule unitaire de 1,5 × 1,5 µm², largeur de la bande de SWCNT de 50 nm, longueur de la bande de SWCNT de 700 nm et espacement entre les bandes de 120 nm. Ces dimensions assurent un fonctionnement sous-longueur d'onde dans la bande 0,5-30 THz.
3. Configuration de la simulation dans CST Microwave Studio
   1. Lancer CST Microwave Studio et créer un nouveau projet en utilisant le solveur en domaine fréquentiel Micro-ondes et RF.
   2. Construire la géométrie de la cellule unitaire en utilisant les paramètres définis. Appliquer des conditions aux limites de cellule unitaire dans les directions x et y avec des frontières périodiques. Définir des ports de Floquet selon la direction z pour l'excitation par onde plane.
   3. Définir la plage de balayage en fréquence de 0,5 THz à 30 THz avec un pas de 0,05 THz, générant 590 points de fréquence. Configurer le solveur temporel avec un maillage adaptatif et fixer l'objectif de convergence à -40 dB.
   4. Exécuter la simulation électromagnétique en régime d'ondes complètes. Exporter le coefficient de réflexion complexe S₁₁, le coefficient de transmission S₂₁, la phase de réflexion et les spectres d'absorption sous forme de fichiers CSV.

2. Modélisation de la conductivité quantique à l'aide du formalisme de Kubo

1. Initialisation des constantes physiques
   1. Initialiser les constantes physiques : charge élémentaire e=1,602×10⁻¹⁹C, constante de Planck réduite ħ=1,0546×10⁻³⁴J·s, constante de Boltzmann k_B=1,3806×10⁻²³J·K^-1, température T=300K, vitesse de Fermi v_F=8×10⁵m·s^-1, et impédance du vide η₀=377Ω.
2. Calcul de la conductivité intra-bande
   1. Calculer la conductivité intra-bande pour chaque point de fréquence à l’aide de :
      
      où τ=0,5ps est le temps de relaxation des porteurs et E_F est l’énergie de Fermi (variant de 0,1eV à 0,4eV).
3. Calcul de la conductivité inter-bande
   1. Calculer la conductivité inter-bande à l’aide de :
      
   2. Calculer la conductivité surfacique totale selon σ_total(ω) = σ_intra(ω) + σ_inter(ω).
4. Calcul de la réponse optique
   1. Calculer le déphasage du substrat pour le quartz (ε_r=3,8, épaisseur=500nm) à l’aide de φ_sub=2πf(n_sub)d_sub/c, où n_sub=.
   2. Calculer le coefficient de réflexion à partir de la conductivité à l’aide de
      
   3. Calculer la réflectance selon R(ω)=|S₁₁(ω)|². Calculer l’absorption selon A(ω)=1-R(ω)-T(ω).

3. Génération et prétraitement des données CST

1. Génération de données CST simulées
   1. Générer des résultats de simulation en onde complète CST simulés en ajoutant des résonances réalistes au S₁₁ quantique.
   2. Définir cinq résonances : exciton E₁₁ à 2,8 THz (amplitude -15,2 dB, phase 45°), exciton E₂₂ à 5,6 THz (-22,8 dB, -120°), résonance plasmonique à 12,4 THz (-8,5 dB, 60°), mode de cavité à 18,9 THz (-12,1 dB, -30°) et résonance assistée par phonon à 24,7 THz (-6,8 dB, 15°).
   3. Ajouter un bruit numérique avec un niveau de bruit de 0,03 (RSB = 30,5 dB) pour simuler les fluctuations numériques de CST.
2. Prétraitement des données
   1. Normaliser l'axe de fréquence à [0,1] en utilisant
      
   2. Extraire les parties réelle et imaginaire de S₁₁ ainsi que l'amplitude de S₂₁.
   3. Diviser l'ensemble de données en ensemble d'apprentissage (80 %), de validation (10 %) et de test (10 %) à l'aide d'une permutation aléatoire. Normaliser les entrées à moyenne nulle et variance unitaire.

4. Régression polynomiale pour la lissage des données

1. Adapter un 8^th ajuster un polynôme d'ordre 8 à la partie réelle de S₁₁ à l'aide de polyfit : P_real = polyfit(X_train, Y_train(:,1), 8). Ajuster un polynôme d'ordre 8^th polynôme d'ordre à la partie imaginaire : P_imag = polyfit(X_train, Y_train(:,2), 8). Ajuster un polynôme d'ordre 8^th ordre du polynôme pour l'amplitude de S₂₁ : P_mag = polyfit(X_train, Y_train(:,3), 8). Générer des prédictions lissées à l'aide de polyval sur toute la plage de fréquences. Calculer S₁₁ lissé comme S₁₁_ann = S11_real_ann + i × S11_imag_ann.

5. Ajustement de la fonction de transfert dans le domaine Z

1. Définir la fréquence d'échantillonnage F_s=60THz (2×fréquence maximale selon le critère de Nyquist). Associer les fréquences au domaine Z à l'aide de z=e^{i2πf/Fs}.
2. Exprimer la fonction de transfert discrète dans le domaine Z comme suit :
   
   où l'ordre du numérateur est n=6 et l'ordre du dénominateur est m=7.
3. Pour l'optimisation des moindres carrés, établir les équations linéaires A × x = B, où A contient les termes du numérateur et du dénominateur, et B contient les valeurs de S₁₁ prédites par le réseau de neurones artificiels (ANN). Résoudre le système à l'aide d'une décomposition QR pour assurer la stabilité numérique :
   
4. Analyse des pôles et zéros, et stabilité
   1. Extraire les pôles en résolvant le polynôme du dénominateur. Assurer la stabilité en projetant les pôles pour lesquels |pôle| ≥ 1 à l'intérieur du cercle unité à l'aide de la relation : pôles(p) = pôles(p) / (|pôles(p)| + 0,1).
   2. Reconstruire le dénominateur à partir des pôles stabilisés. Extraire les zéros en résolvant le polynôme du numérateur. Évaluer H(z) sur une grille de fréquences et calculer l'erreur quadratique moyenne.

6. Optimisation du pilotage du faisceau

1. Le facteur de réseau normalisé pour un réseau linéaire à balayage électronique est donné par :
   ​
   où I_n = 1 pour une excitation uniforme, k = 2π/λ, d = λ/2 d'espacement, β_n est le déphasage progressif. Pour des angles de pointage cibles de ±45°, ±30° et 0°, calculez le déphasage entre les éléments adjacents selon .EQUATION
   1. Calculez les diagrammes en champ lointain et évaluez les niveaux des lobes secondaires ainsi que la largeur du faisceau à mi-puissance.

Résultats

Sélection de la chiralité des SWCNT
Le protocole décrit a été mis en œuvre pour une cellule unitaire optique RIS à base de SWCNT (10,5) fonctionnant dans la bande de 0,5 à 30 THz. Les résultats représentatifs démontrent l'efficacité de l'approche hybride de modélisation polynomiale-Z pour une caractérisation électromagnétique précise et une optimisation du guidage du faisceau.

Analyse de la conductivité quantique
Le formalisme de Kubo a révélé que le S...

Discussion

Les étapes critiques du protocole nécessitent une attention particulière pour assurer une mise en œuvre réussie. Premièrement, la sélection précise des indices chiraux des nanotubes de carbone à paroi unique (SWCNT) est essentielle, car la bande interdite et la réponse optique dépendent fortement de la chiralité. La chiralité (10,5) indiquée dans ce protocole offre une bande interdite optimale de 1,762 eV pour le fonctionnement en térahertz, mais les utilisateurs visant des bandes de fréquence différentes doivent ca...

Matériaux

Liste des matériaux utilisés dans cet article

Nom	Entreprise	Numéro de catalogue	Commentaires
CST Microwave Studio	Dassault Systèmes	N/A	Version 2024, solveur en domaine fréquentiel
MATLAB	MathWorks	N/A	Version R2014a ou ultérieure
Substrat de quartz	University Wafer	4526	épaisseur de 500 nm, ε_r = 3.8
Chiralité SWCNT (10,5)	NanoIntegris	SWCNT-105	diamètre de 0,60 nm, >90 % semiconducteur
Ordinateur personnel	N/A	N/A	32 Go de RAM, 8 cœurs de processeur au minimum