Prédiction de l'impédance de surface de structures hybrides méta‑surface‑graphène dans le régime térahertz

Résumé

Le contrôle précis de l'impédance de surface dans les méta‑surfaces hybrides à base de graphène est essentiel pour l'optoélectronique térahertz (THz) de prochaine génération. Dans ce travail, nous dérivons des expressions analytiques en forme fermée pour l'impédance de méta‑surfaces à patches carrés et les combinons avec l'impédance de feuille intrinsèque du graphène pour modéliser l'impédance hybride globale. Des simulations à ondes complètes confirment les prédictions analytiques sur une large gamme de fréquences. Nous étudions comment les dimensions des patches et le potentiel chimique du graphène influencent l'impédance composite, fournissant des directives de conception pour les absorbeurs et autres dispositifs THz exploitant l'impédance de surface ingénierée.

Introduction

Les méta‑surfaces d'impédance artificielles ont récemment attiré une attention significative en raison de leur capacité unique à manipuler les ondes électromagnétiques, permettant des applications telles que l'holographie, l'imagerie haute résolution, le camouflage tapis et l'absorption large bande. Cependant, la nature intrinsèquement dispersive des méta‑surfaces conventionnelles limite leur performance à un fonctionnement bande étroite. Le graphène, avec sa conductivité ajustable par tension, offre une plateforme dynamique pour surmonter cette limitation, permettant un contrôle en temps réel des caractéristiques plasmoniques et de méta‑surface à travers les spectres THz et optiques. Par conséquent, les structures hybrides méta‑surface‑graphène ont émergé comme candidats prometteurs pour des composants THz compacts et accordables.

Des modèles analytiques pour l'impédance de surface équivalente des méta‑surfaces et feuilles de graphène ont été développés dans la littérature, pourtant un modèle unifié capturant avec précision l'interaction entre une méta‑surface à patches carrés et une feuille de graphène reste insaisissable. Dans cette étude, nous comblons cette lacune en dérivant un cadre analytique complet et en le validant contre des simulations à ondes complètes.

Méthodes

Impédances pour patches carrés et feuilles de graphène

L'absorbeur canonique méta‑surface‑graphène consiste en une couche structurée mince au‑dessus d'un plan de masse métallique, séparée par un espaceur diélectrique d'épaisseur h. La cellule unité est un patch carré de côté D avec une fente étroite de largeur w (Fig. 1a). Sous la condition D ≪ λ, un circuit équivalent ligne de transmission peut être construit (Fig. 1b). L'impédance d'entrée de la cellule unité s'exprime comme

\(\displaystyle \frac{1}{Z_{\text{in}}}=\frac{1}{Z_{1}}+\frac{1}{Z_{\text{mg}}}=\frac{1}{jZ_{h}\tan(k_{zh}h)}+\frac{1}{Z_{\text{mg}}}\) (1)

où Z_h et k_zh sont l'impédance du substrat et la constante de propagation, respectivement. L'absorptivité globale à incidence normale découle du coefficient de réflexion S₁₁ :

\(\displaystyle A(\omega)=1-|S_{11}|^{2}=1-\left|\frac{Z_{\text{in}}-120\pi}{Z_{\text{in}}+120\pi}\right|^{2}\) (2)

L'Équation (1) permet l'extraction de l'impédance de surface hybride Z_mg à partir des S-paramètres simulés, révélant l'interaction entre géométrie du patch et conductivité du graphène.

Impédance des patches carrés

Pour un réseau dense électriquement de patches parfaitement conducteurs, l'impédance de grille est donnée par

\(\displaystyle Z_{m}=\frac{D}{w}Z_{s}-j\frac{\eta_{\text{eff}}}{2\alpha}\) (3)

avec l'impédance d'onde effective η_eff = √(μ₀/ε₀ε_eff) et le facteur géométrique

\(\displaystyle \alpha=\frac{k_{\text{eff}}D}{\pi}\ln\left[\frac{1}{\sin(\pi w/2D)}\right]\) (4)

où ε_eff ≈ (ε_r+1)/2 et k_eff = k₀√ε_eff. Ces expressions valent quand la période du patch est bien plus petite que la longueur d'onde de fonctionnement.

Pour valider le modèle analytique, nous avons effectué des simulations à ondes complètes dans Ansoft HFSS avec une feuille résistive (Z_s = 35 Ω/sq), épaisseur d'espaceur 20 µm, et permittivité relative ε_r = 3,2(1−j0,045). L'impédance de patch extraite de l'impédance d'entrée simulée correspondait étroitement à la prédiction analytique, confirmant l'exactitude du modèle.

Impédance des feuilles de graphène

Le graphène est traité comme une surface conductrice infiniment mince. Sa conductivité de surface est décrite par la formule de Kubo :

\(\displaystyle \sigma_{g}=\frac{j e^{2}k_{B}T}{\pi\hbar^{2}(\omega+j/\tau)}\left[\frac{\mu_{c}}{k_{B}T}+2\ln\left(e^{-\mu_{c}/k_{B}T}+1\right)\right]+\frac{j e^{2}}{4\pi\hbar}\ln\left[\frac{2|\mu_{c}|-\hbar(\omega+j/\tau)}{2|\mu_{c}|+\hbar(\omega+j/\tau)}\right]\) (5)

En supposant température ambiante (T = 300 K) et temps de relaxation τ = 0,1 ps, l'impédance de feuille est Z_g = 1/σ_g = R_g + jX_g. La Figure 4 montre que la résistance et la réactance diminuent monotoni­quement avec l'augmentation du potentiel chimique μ_c, tout en restant presque plates sur la bande 0,2–6 THz pour μ_c fixé.

Résultats et discussion

Contrairement à l'hypothèse conventionnelle que l'impédance hybride est simplement la combinaison en parallèle des impédances patch et graphène, nos simulations révèlent une interaction plus complexe. La Figure 5 compare les valeurs analytiques et simulées des parties réelle et imaginaire de Z_mg pour divers potentiels chimiques (w = 19 µm). Bien que la tendance analytique capture la dépendrance fréquentielle globale, les valeurs absolues diffèrent significativement, surtout à bas μ_c où l'impédance du graphène domine.

Pour élucider davantage la dépendance à la taille, nous avons varié la largeur du patch de 17 à 19,5 µm et extrait Z_mg pour μ_c fixé = 0,4 eV (Fig. 7). La partie réelle de l'impédance diminue avec l'augmentation de w jusqu'à 0,31 THz, puis augmente, reflétant l'équilibre entre couplage capacitif et résistance de feuille du graphène. La partie imaginaire suit une tendance similaire, confirmant l'interaction entre géométrie et réponse matérielle.

Les cartes de densité de courant à 3 THz (Fig. 8) illustrent que les patches plus grands supportent des courants de surface plus faibles, cohérent avec la relation inverse entre longueur de patch et magnitude de courant donnée par l'Éq. (11). Les valeurs de courant intégrées depuis HFSS corroborent cette tendance, renforçant l'interprétation physique des variations d'impédance.

Conclusions

Nous avons développé et validé un cadre analytique complet pour prédire l'impédance de surface des hybrides méta‑surface‑graphène dans le régime THz. Les formules dérivées pour l'impédance de patch carré, combinées avec l'impédance de feuille de graphène basée sur Kubo, permettent une conception rapide d'absorbeurs, antennes et autres dispositifs sans simulations à ondes complètes exhaustives. L'étude clarifie aussi comment la taille du patch et le potentiel chimique du graphène influencent conjointement l'impédance effective, offrant des leviers de conception clairs pour composants THz accordables.

Abréviations

HFSS :

Simulation structure haute fréquence

TEM :

Électromagnétique transverse

THz :

Térahertz