Ultrakompakter Terahertz-Plasmonenwellenleiter mit verbesserter Einsperrung durch Bulk-Dirac-Halbmetall-Isolator-Metall-Struktur

Zusammenfassung

Ein subwellenlängen Terahertz-Plasmonenwellenleiter, basierend auf einem Bulk-Dirac-Halbmetall (BDS)-Isolator-Metall (BIM) Stapel, wird untersucht. Die Struktur liefert ein optimiertes Frequenzband, in dem der Oberflächenplasmonenpolariton (SPP) Modus stark eingedämmt ist – bis zu einem Einsperrfaktor von ½·λ₀/15 – während der Ausbreitungsverlust niedrig bei 1,0 dB/λ₀ bleibt. Die Einführung von zwei Siliziumbändern in den BIM-Wellenleiter schafft einen dynamisch abstimmbaren Fabry-Pérot-Resonator, der Terahertz-SPPs im tief-subwellenlängigen Maßstab formen kann und einen Weg zu ultrakompakten, abstimmbaren THz-Plasmonengeräten und optischen Filtern bietet.

Hintergrund

Terahertz (THz) Strahlung hat sich als vielseitiges Werkzeug in Bildgebung, biochemischem Sensorik und drahtloser Kommunikation etabliert – was immer engere Einsperrung, höhere Auflösung und größere Integration erfordert, um ihr volles Potenzial auszuschöpfen – [1,2,3]. Die Erreichung subwellenlängiger Einsperrung von THz-Wellen ist daher kritisch, doch konventionelle Oberflächenplasmonenpolaritonen (SPPs) in Edelmetallen leiden unter schlechter Einsperrung, hohem intrinsischem Verlust und passiver Abstimmbarkeit – Einschränkungen, die praktische Anwendungen behindert haben – [4,5,6].

Graphen-Plasmonen haben aufgrund ihres geringen Verlusts, dynamischer Abstimmbarkeit und extremer Einsperrung Versprechen gezeigt, was hochauflösende THz-Geräte ermöglicht – Beispiele sind gate-abstimmbare Graphen-Heterostrukturen, diskrete Talbot-Effekte und plasmon-induzierte Transparenz-Wellenleiter – [11–18]. Diese Fortschritte unterstreichen die entscheidende Rolle extremer SPP-Einsperrung für die Manipulation von THz-Strahlung im tief-subwellenlängigen Maßstab.

Bulk-Dirac-Halbmetalle (BDS), oft als „3-D-Graphen“ bezeichnet, weisen ultrahohe Trägerbeweglichkeiten (bis zu 9×10⁶ cm² V⁻¹ s⁻¹) auf, die das beste Graphen (2×10⁵ cm² V⁻¹ s⁻¹) weit übertreffen – ein Schlüsselvorteil für niedrigverlustige Plasmonen – [19]. Ihre dielektrische Antwort ist auch über Fermi-Energie-Modulation abstimmbar, und Materialien wie Na₃Bi, Cd₃As₂ und AlCuFe-Quasikristalle sind einfacher zu verarbeiten und stabiler als Graphen – was BDS als vielversprechendes nächstgenerationelles Plasmonenmaterial positioniert – [20,21]. Allerdings ist die SPP-Moden-Einsperrung an einer reinen BDS-Isolator-Grenzfläche noch unzureichend für tief-subwellenlängige Manipulation.

In dieser Studie präsentieren wir einen BIM-Wellenleiter, der die SPP-Einsperrung markant verbessert, während er niedrigen Verlust und dynamische Abstimmbarkeit beibehält. Die Dispersion, Verlustcharakteristiken und Filterfähigkeiten dieses hoch eingedämmten Modus werden systematisch untersucht, was ein optimiertes Frequenzfenster mit sowohl verbesserter Einsperrung als auch reduziertem Ausbreitungsverlust offenbart – ein Verhalten, das in traditionellen metallbasierten SPPs selten beobachtet wird. Zudem kann der BIM-Modus durch Schlitze schmaler als λ₀/2000 propagieren, und durch Einbetten von Siliziumbändern als Spiegel realisieren wir einen abstimmbaren Bandpassfilter, dessen Resonanz durch Anpassung der BDS-Fermi-Energie verschoben werden kann.

Theorie und Simulation

Der BIM-Wellenleiter besteht aus einem Monolagen-BDS-Film (Dicke 0,2 µm), der von einem Silbersubstrat durch einen dielektrischen Spacer (Permittivität εᵣ) getrennt ist. Im THz-Regime verhält sich Silber als perfekter elektrischer Leiter. Für TM-polarisierte Anregung propagiert der SPP-Modus entlang x mit Wellenzahl k_SPP und fällt entlang y in den Freiraum ab. Anwendung der Randbedingungen liefert die Dispersionsrelation:

$$-\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),$$

wobei k₀ der Freiraum-Wellenzahlvektor und σ die komplexe Leitfähigkeit des BDS-Films ist (Gleichungen 3–4). Das Lösen von Gleichung (1) gibt den effektiven Brechungsindex n_eff = k_SPP/k₀ = Re(n_eff) + i Im(n_eff). Re(n_eff) quantifiziert die Moden-Einsperrung, während Im(n_eff) direkt mit dem Ausbreitungsverlust zusammenhängt.

Für hinreichend große Lücken (tanh[…]≈1) vereinfacht sich Gleichung (1) zur Einzelschicht-BDS-Dispersion (Gleichung 2), was den zusätzlichen Nutzen des Metallsubstrats in der BIM-Konfiguration hervorhebt.

Ergebnisse und Diskussion

Zunächst untersuchen wir, wie der BIM-Abstand g und die BDS-Fermi-Energie E_F die Moden-Einsperrung und den Verlust beeinflussen. Mit E_F=70 meV sind die effektiven Brechungsindizes für verschiedene g-Werte in Abb. 2a (Re(n_eff)) und Abb. 2b (Im(n_eff)) dargestellt. Für g≥10 µm konvergieren die Kurven oberhalb von 0,05 THz, was anzeigt, dass das SPP-Feld bereits innerhalb von ~10 µm eingedämmt ist und der Silberbeitrag vernachlässigbar ist. Die Reduktion von g unter 1 µm erhöht die Einsperrung markant und reduziert den Verlust, was ein optimiertes Frequenzband ergibt, in dem die Einsperrung ihren Peak erreicht, während der Verlust abnimmt – ein Merkmal, das in konventionellen metallbasierten SPPs ungewöhnlich ist.

Abbildung 2c,d illustrieren die Abhängigkeit von E_F (g=1 µm). Erhöhung von E_F senkt sowohl Einsperrung als auch Verlust, was auf erhöhte Metallizität und längere Trägerrelaxationszeiten zurückzuführen ist. Beispielsweise werden bei 2,5 THz eine Einsperrung von λ₀/15 und ein Verlust von 1,0 dB/λ₀ mit g=10 nm und E_F=70 meV erreicht.

Transmissionssimulationen (Abb. 3) bestätigen den reduzierten Verlust: bei 1,56 THz zeigt der BIM-Wellenleiter eine Transmission von 0,97, überlegen dem Monolagen-BDS-Fall. Die entsprechenden Magnetfeldverteilungen (Abb. 3b,c) zeigen ein höheres Re(n_eff)=2,45 für BIM versus 1,002 für Monolage, mit dem Feld eng in Schlitzen <λ₀/2000 lokalisiert.

Um den hoch eingedämmten Modus zu nutzen, integrieren wir zwei Siliziumbänder (n_Si=3,4) in den dielektrischen Spacer, die reflektierende Spiegel bilden (Abb. 4a). Der resultierende Fabry-Pérot-Resonator unterstützt Stehwellenresonanzen, was zwei Transmissionspeaks bei 1,56 und 2,22 THz mit FWHM von 0,12 und 0,09 THz ergibt. Magnetfeldkarten (Abb. 4c,e) bestätigen erste und zweite Kavitätenmoden. Außerhalb der Resonanz können keine Stehwellen gebildet werden, und die Eingabe wird reflektiert (Abb. 4d).

Die Transmission kann analytisch durch gekoppelte-Moden-Theorie (CMT) beschrieben werden:

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

wobei κ_w und κ_i die Wellenleiter-Kopplungs- und intrinsischen Verlustraten sind. Die Qualitätsfaktoren Q_t=ω₀/FWHM und Q_oi=−Re(n_eff)/(2Im(n_eff)) ermöglichen die Extraktion von κ_w und κ_i, wobei die analytischen Ergebnisse die numerischen Spektren matchen (Abb. 4b).

Abbildung 5 untersucht die Abhängigkeit der Resonanzfrequenz von der Kavitätenlänge L. Die erste und zweite Resonanz rotschiebt mit zunehmender L, konsistent mit der Stehwellenbedingung 2k_SPP(ω_r)L+θ=2mπ. Ähnlich blauverschiebt eine Erhöhung von g die Resonanzen (Abb. 6a) aufgrund reduzierten Re(n_eff), während eine Erhöhung von E_F die Peaks ebenfalls blauverschiebt (Abb. 6b) und die Bandbreite verengt, was niedrigeres Im(n_eff) reflektiert. Diese Abstimmbarkeitsmechanismen erlauben dynamische Kontrolle des Filters ohne Neudesign der Struktur.

Schlussfolgerungen

Der BIM-Wellenleiter liefert einen ultrakompakten Terahertz-Plasmonenmodus mit λ₀/15 Einsperrung und 1,0 dB/λ₀ Verlust, erreichbar in Schlitzen schmaler als λ₀/2000. Durch Einbau von Siliziumbändern realisieren wir einen dynamisch abstimmbaren Bandpassfilter, dessen Resonanzfrequenz sich mit der BDS-Fermi-Energie verschiebt, was THz-Schalten und -Filtern ohne strukturelle Neoptimierung ermöglicht.

Methoden

Numerische Simulationen nutzten einen 2-D finiten-Differenzen-Zeitbereich (FDTD) Löser mit perfekt angepassten Schichten in x- und y-Richtung. Das BDS-Gitter wurde auf 1 µm × 0,02 µm für Konvergenz gesetzt. Die frequenzabhängige Leitfähigkeit von BDS folgt der Kubo-Formel unter der Zufallsphasenapproximation – Gleichungen (3) und (4) beschreiben Reσ(Ω) und Imσ(Ω). Parameter für AlCuFe waren: t=40, ε_c=3, μ=3×10⁴ cm² V⁻¹ s⁻¹, und E_F=70 meV.

Keine menschlichen Teilnehmer, Daten oder Tiere waren an dieser Forschung beteiligt.

Abkürzungen

BDS:

Bulk-Dirac-Halbmetalle

BIM:

BDS-Isolator-Metall

CMT:

Gekoppelte-Moden-Theorie

FDTD:

Finite-Differenzen-Zeitbereich

FWHM:

Volle Breite bei halbem Maximum

SPPs:

Oberflächenplasmonenpolaritonen