Guide d'ondes plasmonique térahertz ultra-compact avec confinement amélioré utilisant une structure semi-métal de Dirac massif-isolant-métal

Résumé

Un guide d'ondes plasmonique térahertz sous-longueur d'onde construit sur une pile semi-métal de Dirac massif (BDS)–isolant–métal (BIM) est investigué. La structure délivre une bande de fréquence optimisée où le mode polariton plasmonique de surface (SPP) est étroitement confiné — jusqu'à un facteur de confinement de ½·λ₀/15 — tandis que la perte de propagation reste faible à 1,0 dB/λ₀. L'introduction de deux rubans de silicium dans le guide BIM crée un résonateur Fabry-Pérot dynamiquement ajustable qui peut façonner les SPPs térahertz à une échelle sous-longueur d'onde profonde, offrant une voie vers des dispositifs plasmoniques THz ultra-compacts, ajustables et filtres optiques.

Contexte

Le rayonnement térahertz (THz) a émergé comme outil polyvalent en imagerie, détection biochimique et communications sans fil — exigeant confinement toujours plus serré, résolution plus haute et intégration plus grande — pour exploiter pleinement son potentiel — [1,2,3]. L'atteinte de confinement sous-longueur d'onde des ondes THz est donc critique, pourtant les polaritons plasmoniques de surface (SPPs) conventionnels dans les métaux nobles souffrent de confinement pauvre, perte intrinsèque élevée et réglabilité passive — limitations qui ont entravé les applications pratiques — [4,5,6].

Les plasmons de graphène ont montré du potentiel grâce à leur faible perte, réglabilité dynamique et confinement extrême, permettant des dispositifs THz haute résolution — exemples incluent hétérostructures de graphène à grille ajustable, effets Talbot discrets et guides d'ondes transparence plasmonique induite — [11–18]. Ces avancées soulignent le rôle pivotal du confinement SPP extrême pour manipuler le rayonnement THz à échelles sous-longueur d'onde profonde.

Les semi-métaux de Dirac massifs (BDS), souvent surnommés « graphène 3D », exhibent des mobilités de porteurs ultra-élevées (jusqu'à 9×10⁶ cm² V⁻¹ s⁻¹) qui surpassent largement le meilleur graphène (2×10⁵ cm² V⁻¹ s⁻¹) — un avantage clé pour les plasmons à faible perte — [19]. Leur réponse diélectrique est aussi ajustable via modulation de l'énergie de Fermi, et matériaux comme Na₃Bi, Cd₃As₂ et quasicristaux AlCuFe sont plus faciles à traiter et plus stables que le graphène — positionnant le BDS comme matériau plasmonique prometteur de nouvelle génération — [20,21]. Cependant, le confinement du mode SPP à une interface BDS-isolant seul reste insuffisant pour la manipulation sous-longueur d'onde profonde.

Dans cette étude, nous présentons un guide BIM qui améliore nettement le confinement SPP tout en maintenant faible perte et en permettant la réglabilité dynamique. La dispersion, les caractéristiques de perte et les capacités de filtrage de ce mode hautement confiné sont systématiquement investiguées, révélant une fenêtre de fréquence optimisée avec confinement amélioré et perte de propagation réduite — comportement rarement observé dans les SPPs à base de métaux traditionnels. De plus, le mode BIM peut se propager à travers des fentes plus étroites que λ₀/2000, et en intégrant des rubans de silicium comme miroirs, nous réalisons un filtre passe-bande ajustable dont la résonance peut être décalée en ajustant l'énergie de Fermi du BDS.

Théorie et Simulation

Le guide BIM consiste en un film monocouche BDS (épaisseur 0,2 µm) séparé d'un substrat argent par un espaceur diélectrique (permittivité εᵣ). Dans le régime THz, l'argent se comporte comme un conducteur électrique parfait. Pour l'excitation TM-polarisée, le mode SPP se propage le long de x avec vecteur d'onde k_SPP et décroît le long de y dans l'espace libre. L'application des conditions aux limites donne la relation de dispersion :

$$-\frac{\varepsilon_r\sqrt{k_{\mathrm{SPP}}^2-k_0^2}}{\varepsilon_0\sqrt{k_{\mathrm{SPP}}^2-\frac{\varepsilon_r k_0^2}{\varepsilon_0}}}=\left(1+\frac{i\sigma \sqrt{k_{\mathrm{SPP}}^2-k_0^2}}{{\omega \varepsilon}_0}\right)\tanh\left(g\sqrt{k_{\mathrm{SPP}}^2-\frac{\varepsilon_r k_0^2}{\varepsilon_0}}\right),$$

où k₀ est le vecteur d'onde d'espace libre et σ est la conductivité complexe du film BDS (Éq. 3–4). La résolution de l'Éq. (1) donne l'indice de réfraction effectif n_eff = k_SPP/k₀ = Re(n_eff) + i Im(n_eff). Re(n_eff) quantifie le confinement du mode, tandis que Im(n_eff) relie directement à la perte de propagation.

Pour gaps suffisamment grands (tanh[…]≈1), l'Éq. (1) se simplifie en la dispersion monocouche BDS (Éq. 2), soulignant le bénéfice additionnel du substrat métallique dans la configuration BIM.

Résultats et Discussion

Nous examinons d'abord comment le gap BIM g et l'énergie de Fermi du BDS E_F affectent le confinement et la perte. Avec E_F=70 meV, les indices de réfraction effectifs pour diverses valeurs de g sont tracés Fig. 2a (Re(n_eff)) et Fig. 2b (Im(n_eff)). Pour g≥10 µm, les courbes convergent au-dessus de 0,05 THz, indiquant que le champ SPP est déjà confiné dans ~10 µm et la contribution de l'argent est négligeable. Réduire g sous 1 µm augmente nettement le confinement et réduit la perte, produisant une bande de fréquence optimisée où le confinement pic alors que la perte diminue — caractéristique peu commune dans les SPPs conventionnels à base de métaux.

La figure 2c,d illustre la dépendance à E_F (g=1 µm). L'augmentation d'E_F abaisse à la fois confinement et perte, attribuable à la métallicité accrue et temps de relaxation des porteurs plus longs. Par exemple, à 2,5 THz, un confinement de λ₀/15 et une perte de 1,0 dB/λ₀ sont atteints avec g=10 nm et E_F=70 meV.

Les simulations de transmission (Fig. 3) confirment la perte réduite : à 1,56 THz, le guide BIM exhibe une transmission de 0,97, surpassant le cas monocouche BDS. Les distributions de champ magnétique correspondantes (Fig. 3b,c) montrent un Re(n_eff)=2,45 plus élevé pour BIM versus 1,002 pour monocouche, avec le champ étroitement localisé dans des fentes <λ₀/2000.

Pour exploiter le mode hautement confiné, nous intégrons deux rubans de silicium (n_Si=3,4) dans l'espaceur diélectrique, formant des miroirs réfléchissants (Fig. 4a). La cavité Fabry-Pérot résultante supporte des résonances d'onde stationnaire, donnant deux pics de transmission à 1,56 et 2,22 THz avec FWHM de 0,12 et 0,09 THz respectivement. Les cartes de champ magnétique (Fig. 4c,e) confirment les modes de cavité d'ordre 1 et 2. Hors résonance, les ondes stationnaires ne peuvent se former et l'entrée est réfléchie (Fig. 4d).

La transmission peut être décrite analytiquement par la théorie des modes couplés (CMT) :

$$T(\omega)=\frac{\kappa_w^2}{(\omega-\omega_0)^2-(\kappa_w+\kappa_i)^2},$$

où κ_w et κ_i sont les taux de couplage guide et de perte intrinsèque. Les facteurs de qualité Q_t=ω₀/FWHM et Q_oi=−Re(n_eff)/(2Im(n_eff)) permettent l'extraction de κ_w et κ_i, les résultats analytiques correspondant aux spectres numériques (Fig. 4b).

La figure 5 explore la dépendance de la fréquence de résonance à la longueur de cavité L. La première et seconde résonances se décalent vers le rouge avec l'augmentation de L, consistant avec la condition d'onde stationnaire 2k_SPP(ω_r)L+θ=2mπ. Similairement, l'augmentation de g décale les résonances vers le bleu (Fig. 6a) due à la réduction de Re(n_eff), tandis que l'augmentation d'E_F décale aussi les pics vers le bleu (Fig. 6b) et rétrécit la bande passante, reflétant un Im(n_eff) plus faible. Ces mécanismes de réglabilité permettent un contrôle dynamique du filtre sans redessiner la structure.

Conclusions

Le guide BIM délivre un mode plasmonique térahertz ultra-compact avec confinement λ₀/15 et perte 1,0 dB/λ₀, atteignable dans des fentes plus étroites que λ₀/2000. En incorporant des rubans de silicium, nous réalisons un filtre passe-bande dynamiquement ajustable dont la fréquence de résonance décalée avec l'énergie de Fermi du BDS, permettant commutation et filtrage THz sans ré-optimisation structurelle.

Méthodes

Les simulations numériques employaient un solveur 2D différences finies domaine temporel (FDTD) avec couches parfaitement adaptées dans directions x et y. Le maillage BDS était fixé à 1 µm × 0,02 µm pour convergence. La conductivité fréquence-dépendante du BDS suit la formule de Kubo sous approximation de phase aléatoire — Éqs. (3) et (4) décrivent Reσ(Ω) et Imσ(Ω). Paramètres pour AlCuFe : t=40, ε_c=3, μ=3×10⁴ cm² V⁻¹ s⁻¹, et E_F=70 meV.

Aucun participant humain, données ou animaux n'étaient impliqués dans cette recherche.

Abréviations

BDS

Semi-métaux de Dirac massifs

BIM

BDS-isolant-métal

CMT

Théorie des modes couplés

FDTD

Différences finies domaine temporel

FWHM

Largeur à mi-hauteur

SPPs

Polaritons plasmoniques de surface