Ultra-breitbandige Terahertz-Polarisationskontrolle durch anisotrope dielektrische Metasurfaces

Abstract

Wir demonstrieren Terahertz-dielektrische Metasurfaces mit anisotropen Multipolen im Rahmen des verallgemeinerten Huygens-Prinzips, wobei die Interferenz zwischen diesen Multipolen riesige Phasenschiebungen mit verbreiteter Bandbreite und hohen Transmissionskoeffizienten erreicht. Noch wichtiger: Durch das anisotrope Design werden verschiedene Phasenverzögerungen zwischen π/2 und 3π/2 erzielt, die linear polarisierte Terahertzwellen in rechts-/linkszirkular polarisiertes Licht, elliptisch polarisiertes Licht und kreuzpolarisiertes Licht umwandeln. Sowohl Simulationen als auch Experimente verifizieren die vollständige Terahertz-Polarisationskontrolle mit Elliptizität von 1 bis −1, was den Weg für polarisationsbezogene Anwendungen von Terahertz-Metageräten ebnet.

Einleitung

Polarisation ist ein Schlüsselparameter zur Quantifizierung des Zustands elektromagnetischer Wellen. Besonders Polarisationskontrolle im Terahertzbereich hat großes Forschungsinteresse aufgrund potenzieller Anwendungen in der Terahertztechnologie geweckt. Allerdings ist Terahertzstrahlung aus den meisten Quellen linear polarisiert, was komplexe polarimetrische Terahertzsysteme nicht erfüllen kann. Konventionelle Ansätze zur Manipulation der Terahertzpolarisation nutzen birefringente Materialien, die unter Nachteilen wie voluminöser Größe und schmalbandigem Betrieb leiden – Hindernisse für die Integration in moderne kompakte, breitbandige Terahertz-photonische Systeme.

In jüngerer Zeit ermöglichen Metasurfaces – künstlich designte Antennenarrays – einen effizienten Ansatz zur Formung der Terahertzpolarisation. Verschiedene metallische Metasurfaces wurden vorgeschlagen, basierend auf anisotropen, chiralen oder mehrlagigen Designs. Aktive Medien (Phasenwechselmaterialien, Halbleiter, 2D-Materialien, Flüssigkristalle, Supraleiter) wurden integriert, um Funktionalitäten zu erweitern. Mehrlagige und aktive metallische Metasurfaces können die Leistung steigern, jedoch auf Kosten hoher Verluste und komplexer Fertigung. Kürzlich bieten dielektrische Metasurfaces, bestehend aus dielektrischen Antennen, einen neuen Ansatz zur Kontrolle elektromagnetischer Wellen. Unterstützt durch Interferenz zwischen elektrischen und magnetischen Mie-Resonanzen erreichen dielektrische Metasurfaces 2π-Phasenkontrolle mit hoher Effizienz. Große Anstrengungen zielten auf die Verbesserung dielektrischer Metasurfaces für Terahertz-Polarisationskontrolle ab. Doch die meisten Arbeiten basieren auf elektrischen und magnetischen Dipolresonanzen, was begrenzte Phasenverzögerungsbereiche und im Prinzip Einfrequenzbetrieb zur Folge hat – was die vollständige Manipulation der Terahertzpolarisation mit hoher Leistung behindert.

Hier schlagen wir dielektrische Metasurfaces mit Multipolen vor, die die Phasendispersion stark anheben (Phasenschub bis 4π) und riesige Phasenverzögerung, verbreiterte Bandbreite und hohe Effizienz ermöglichen, was vollständige Terahertz-Polarisationskontrolle erlaubt. Aus elliptischen Siliciumsäulen-Arrays bestehend, unterstützen die Metasurfaces verschiedene elektrische und magnetische Multipole. Durch Überlagerung dieser Multipole lassen sich nahezu perfekte Transmission in verbreiterter Bandbreite und bis zu 4π Phasenschub mittels verallgemeinertem Huygens-Prinzip realisieren. Aufgrund der Anisotropie der Siliciumsäulen kann die Überlagerung der Multipole unabhängig entlang der kurzen und langen Achsen der elliptischen Säulen variiert werden. So wird riesige Phasenverzögerung in einem breiten Band erreichbar, was eine überlegene Leistung im Vergleich zu anderen metallischen und dielektrischen Designs zeigt (siehe Zusatzdatei 1: Abb. S1). Da unsere Designs vollständige Polarisationskontrolle in einem einfachen Designrahmen ermöglichen, können die Meta-Atome künstlich angeordnet werden, um den Polarisationsgrad räumlich zu variieren und komplexe Terahertzstrahlen zu erzeugen, wie elliptizitätsvariante Vektorfelder.

Design und Simulation

Die von einer dielektrischen Antenne gestreute elektromagnetische Welle kann in Multipole unterschiedlicher Symmetrie zerlegt werden. Werden dielektrische Antennen zu Arrays in Metasurfaces angeordnet, lässt sich das gestreute Feld E als Summe einer symmetrischen Komponente E_s und einer antisymmetrischen Komponente E_as ausdrücken. Daraus leiten sich Transmission und Reflexion der Metasurfaces entlang der Ausbreitungsrichtung ẑ allgemein her (Gleichungen 1–2). Um hohe Transmission und vernachlässigbare Reflexion zu erreichen, müssen E_s(−ẑ) und E_as(−ẑ) in Rückwärtsrichtung gleiche Amplituden, aber entgegengesetzte Phasen besitzen. Insbesondere wenn die dielektrische Antenne zwei Multipole unterstützt – z. B. eine symmetrische Resonanz (elektrischer Dipol) und eine antisymmetrische Resonanz (magnetischer Dipol) – kann ihre Überlagerung die Forderung nach destruktiver Interferenz erfüllen. Diese destruktive Interferenz führt zu nullliger Reflexion, wenn beide Dipolmoden dieselbe Resonanzfrequenz, Amplitude und Gütefaktor besitzen – das Prinzip von Huygens-Metasurfaces. Solche destruktive Interferenz tritt jedoch nur in einem schmalen Band auf, was die Realisierung breitbandiger Bauelemente fundamental beschränkt. Um das Betriebsband zu verbreitern, sollten die gestreuten Felder E_s und E_as Beiträge höherer Multipole enthalten, wobei die resultierende Transmission ein Gleichgewicht der multipolaren Interferenz darstellt. Dieses Szenario ähnelt dem Konzept der verallgemeinerten Kerker-Bedingung mit multipolarer Interferenz.

Um alle Polarisationzustände abzudecken – rechts-/linkszirkulare Polarisation, elliptische Polarisation und lineare Polarisation – sollte die Phasenverzögerung 90° bis 270° abdecken, was Elliptizität von 1 bis −1 entspricht. Wir schlagen daher anisotrope dielektrische Metasurfaces aus elliptischen Siliciumsäulen-Arrays vor (Abb. 1a). Im Terahertzbereich wird intrinsisches Silicium verwendet, um Absorptionsverluste zu eliminieren. Abb. 1a zeigt, wie linear polarisiertes einfallendes Licht in zirkular, elliptisch und kreuzpolarisiertes Licht umgewandelt wird, wenn multipolare Interferenz unterschiedliche Kombinationen bezüglich verschiedener Geometriegrößen beibehält. Abb. 1b zeigt die Einheitszellenparameter: kurze und lange Achsen a und b, Höhe h, Periodizitäten P_x und P_y. Abb. 1c zeigt ein REM-Bild typischer Siliciumsäulen-Arrays.

Abb. 1. a) Schematik der dielektrischen Metasurfaces für vollständige Polarisationskontrolle. b) Einheitszelle. c) REM-Bild typischer Siliciumsäulen-Arrays in Schrägansicht mit Zoom.

Zur Realisierung vollständiger Terahertz-Polarisationskontrolle wurde numerische Simulation zur Optimierung der Dimensionen durchgeführt, die gleichzeitig 90°–270° Phasenverzögerungsvariation mit hoher Transmission und verbreiterter Bandbreite erfüllen. Zwischen 90° und 270° wurde ein Schrittweite von 45° gewählt, um verschiedene Polarisationskontrollschemata zu demonstrieren. Wir benennen die Designs nach ihren Phasenverzögerungen: P90, P135, P180, P225, P270. Die Simulation erfolgte in CST Microwave Studio. Silicium wurde als verlustfreies Dielektrikum mit ε_Si = 11,7 modelliert, der Substrat als verlustfrei mit ε_sub = 4,5. Periodische Randbedingungen entlang x- und y-Achse, offene Randbedingungen vorne und hinten. Abb. 2a zeigt simulierte Transmission und Phasenverzögerungen der fünf Metasurfaces. Detaillierte Parameter finden sich in Zusatzdatei 1: Tabelle S1. Alle Metasurfaces zeigen hohe Transmissionskoeffizienten für sowohl x- als auch y-polarisierten Einfall von 1,2 bis 1,3 THz, während die Phasenverzögerungen 90°, 135°, 180°, 225°, 270° betragen. Gleiche Transmissionskoeffizienten bei 90° Phasenverzögerung bedeuten Umwandlung in linkszirkular polarisiertes (LCP) Licht. Entsprechend decken 135°, 180°, 225°, 270° elliptische, kreuz- und rechtspolarisiertes Licht ab. Damit ist vollständige Polarisationskontrolle der Terahertzwelle in den vorgeschlagenen dielektrischen Metasurfaces erreicht – mit überlegener Leistung gegenüber Meta-Geräten mit begrenzter Bandbreite, niedriger Effizienz und begrenztem Phasenverzögerungsbereich.

Abb. 2. a) Simulierte und b) experimentell gemessene Transmissionskoeffizienten und Phasenverzögerungen der Designs P90, P135, P180, P225, P270.

Ergebnisse und Diskussion

Zur Validierung der Polarisationskontrolle wurden Siliciumsäulen-Arrays gefertigt und in einem Terahertz-Zeitbereichsspektroskop (THz-TDS) charakterisiert. Der Fertigungsprozess ist im Methodenabschnitt beschrieben. Ein dünnes Borosilikatglas (BF33, 300 µm) diente als Substrat. Das REM-Bild eines typischen Samples für das 270°-Design ist in Abb. 1c (schräg mit Zoom) gezeigt. Zur Charakterisierung wurden die elektrischen Felder der transmittierten Terahertzwelle entlang der kurzen und langen Achsen der Siliciumsäulen als E_x und E_y bezeichnet. Ein blankes Glassubstrat diente als Referenz mit E_x(ref) und E_y(ref). Transmissionskoeffizienten: t_x = E_x/E_x(ref), t_y = E_y/E_y(ref). Phasenverzögerung φ = φ_y − φ_x = arg(t_y) − arg(t_x). Details zum Messsystem finden sich im Methodenabschnitt.

Die gemessenen Transmissionskoeffizienten und Phasenverzögerungen sind in Abb. 2b dargestellt. Alle Designs (P90–P270) zeigen hohe Transmissionskoeffizienten in den entworfenen Frequenzbereichen mit Phasenverzögerungen nahe 90°, 135°, 180°, 225°, 270°. Kleine Diskrepanzen zwischen Simulation und Experiment können auf Größenfluktuationen bei der Fertigung zurückzuführen sein. Oberflächenrauheit kann zusätzlichen Verlust und verringerte Transmission verursachen. Substrateffekte (Verluste, Reflexionen) werden in Zusatzdatei 1: Abb. S2 detailliert diskutiert. Dennoch validieren die ähnlichen Variationstrends die Leistung der dielektrischen Metasurfaces für Polarisationskontrolle.

Um die Polarisationskonversionsleistung vollständig zu untersuchen, wurde die Elliptizität der transmittierten Welle berechnet: χ = S₃/S₀, wobei S₀ und S₃ Stokes-Parameter sind, die direkt aus Transmissionskoeffizienten und Phasenverzögerungen berechenbar sind. Abb. 3 zeigt, dass die Simulationsergebnisse eine vollständige Abdeckung der Elliptizität von 1 bis −1 aufweisen. Im Bereich 1,2–1,3 THz zeigen sowohl Simulation als auch Experiment ähnliche Trends. Diskrepanzen um 1,4 THz können zwei Aspekte haben: Erstens wird in der Simulation der Substrat als verlustfrei mit unendlicher Dicke behandelt, im Experiment weist der 300 µm dicke Substrat deutliche Verluste auf, die hohe-Q-Resonanzen dämpfen (z. B. MD bei 1,4 THz) und Transmissionsspektren glätten. Zweitens weichen geometrische Parameter im Experiment von der Simulation ab – z. B. variiert die Säulenbreite allmählich mit der Höhe durch den tiefen reaktiven Ionenätzprozess. Diese Variationen verbreitern die Multipole und erhöhen ihre Überlappung, wodurch einzelne hohe-Q-Resonanzen durch Superposition und Interferenz verschlechtern. Zusammenfassend führen Substrateffekte und Geometrievariationen zu den Diskrepanzen um 1,4 THz. Diese lassen sich durch verlustarme Substrate (Quarz, Polyimid, SU8) mit geringer Dicke und optimierte Fertigung minimieren. Es sei angemerkt, dass die Betriebsfrequenzen generell im Off-Resonanz-Bereich entworfen wurden, weshalb sie schwach von der Verschlechterung hoher-Q-Resonanzen betroffen sind.

Abb. 3. a) Simulierte und b) experimentell gemessene Elliptizität der verschiedenen dielektrischen Metasurfaces.

Zur Veranschaulichung der multipolaren Interferenz wurden Streuquerschnitte (SCSs) verschiedener Multipole durch sphärische Multipolzerlegung bezüglich zweier orthogonaler Polarisationrichtungen berechnet. Details zur Multipolzerlegung finden sich im Methodenabschnitt. Abb. 4 zeigt die berechneten SCSs für die verschiedenen Metasurfaces unter x- und y-polarisiertem Einfall. Für P90 dominiert unter x-Polarisation der magnetische Dipol (MD) bei 1,4 THz, unter y-Polarisation bei 1,18 THz. Im höheren Frequenzbereich (~1,42 THz) zeigen elektrischer Dipol (ED), elektrisches Quadrupol (EQ) und magnetisches Quadrupol (MQ) unter y-Polarisation offensichtliche Beiträge. Im Überlappungsbereich 1,2–1,3 THz sichern Off-Resonanz-Bedingungen hohe Transmissionskoeffizienten, während Interferenz zwischen verschiedenen Multipolen unterschiedliche Phasendispersionskurven für die beiden orthogonalen Polarisationen anhebt. Mit einer geeigneten Balance zwischen den Multipolen lässt sich eine bestimmte Phasenverzögerung mit hohen Transmissionskoeffizienten und verbreiterten Bandbreiten erreichen – in unserem Fall 90°. Für P135, P180, P225 zeigen ED, MD, EQ, MQ ähnliche Trends wie P90 mit subtilen Änderungen der Resonanzfrequenzen und Modenüberlappung, was die Funktionalität multipolaren Interferenz für Polarisationskontrolle klar demonstriert. Im Gegensatz dazu erfordert P270 (270° Phasenverzögerung) riesige Phasendispersion mit hoher Transmission in breitem Band, was über Off-Resonanz-Bedingungen schwer realisierbar ist. Hier designen wir eine In-Resonanz-Bedingung: Unter x-Polarisation tragen ED, MD, MQ zwischen 1,2–1,3 THz signifikant zu den SCSs bei; unter y-Polarisation dominiert MD bei 1,39 THz. So führen multipolare Interferenzeffekte zu hoher Transmission in breitem Band mit 270° Phasenverzögerung.

Abb. 4. Multipolzerlegung der SCSs für ED, MD, EQ, MQ unter a) x- und b) y-polarisiertem Einfall.

Verglichen mit anderen existierenden Designs ermöglicht unser Vorschlag eine Einzellagen-Plattform für vollständige Terahertz-Polarisationskontrolle. Noch wichtiger: Die Phasenverzögerung unseres Designs kann von 90° bis 270° variiert werden, was verschiedene Polarisationzustände – Zirkularpolarisation, elliptische Polarisation, kreuzlineare Polarisation – abdeckt, was in anderen Designs schwer zu erreichen ist (Tabelle 1). Gleichzeitig übertreffen Bandbreite und Effizienz unseres Designs andere Einzellagen-Designs. Zwar zeigen Mehrlagendesigns bessere Leistung, doch erfordern sie komplexe Design- und Fertigungsprozesse, was ihre Anwendung in kompakten Terahertz-optischen Systemen einschränkt. Zudem realisieren unsere Designs verschiedene Polarisationskonversionen, während die meisten Mehrlagendesigns begrenzte Phasenverzögerungen mit单一 Polarisationskonversionsfunktion erreichen.

Fazit

Zusammenfassend haben wir vollständige Terahertz-Polarisationskontrolle mit verbreiterter Bandbreite und hoher Effizienz mittels alldielektrischer Metasurfaces vorgeschlagen und experimentell demonstriert. Aus elliptischen Siliciumsäulen-Arrays bestehend, erreichen die Metasurfaces gleiche und hohe Transmissionskoeffizienten entlang x- und y-Achse, während ihre Phasenverzögerung kontinuierlich von 90° bis 270° in 45°-Schritten einstellbar ist. Die entsprechende Elliptizität ändert sich von 1 bis −1, was eine vollständige Abdeckung verschiedener polarisierten Lichtzustände – LCP, elliptisch, kreuzpolarisiert, RCP – anzeigt. Darüber hinaus verifizieren Multipolzerlegungsergebnisse die unterschiedlichen Beiträge der Multipole zur Polarisationskontrolle. Solche multipolaren-interferenz-unterstützten dielektrischen Metasurfaces versprechen eine exotische Strategie für die Implementierung hochleistungsfähiger Terahertz-funktionaler Polarisationskontrollbauelemente.

Methoden

Die Fertigung der dielektrischen Metasurfaces umfasst Standard-Photolithografie und tiefes reaktives Ionenätzen. Zunächst wurde ein intrinsischer Siliciumwafer (500 µm) durch anodisches Bonden auf einen Glaswafer (BF33, 300 µm) gebondet. Der Siliciumwafer hatte einen Widerstand >5.000 Ω·cm, um Absorptionsverluste im Terahertzbereich zu eliminieren. Der Wafer wurde auf 180 µm gedünnt, dann mit Aceton und deionisiertem Wasser 30 min gereinigt. Fotolack AZ4620 wurde aufgesponnen, weichgebacken (100 °C, 10 min). Nach Belichtung (MA6, 40 s) und Entwicklung (3 min) folgte Hartbacken (110 °C, 5 min). Dann Silicium-Ätzen durch tiefes reaktives Ionenätzen (56 min). Schließlich Fotolackentfernung mit Aceton, Isopropanol und deionisiertem Wasser.

Die Charakterisierung erfolgte im THz-TDS. Terahertzwelle wurde von einem hausgemachten spintronischen Terahertz-Emitter erzeugt, gepumpt von einem 100 fs Pulslaser bei 800 nm (80 MHz). Die Welle wurde durch vier achsenparabolische Spiegel kollimiert und fokussiert. Das Sample wurde im Fokus (Strahl~3 mm) positioniert. Zwei Terahertz-Polarisatoren vor und nach dem Sample kontrollierten die Polarisation. Detektion erfolgte via elektrooptischem Sampling mit 1 mm ZnTe (110) Kristall. Die Probensonde stammt vom gleichen Lasersystem (20 mW). Messung bei Raumtemperatur in Stickstoffatmosphäre zur Entfernung von Wasserabsorption.

Die Multipolzerlegung erfolgte über hausinternen Matlab-Code. Zunächst wurde die elektrische Feldverteilung E_inter(r̂) im elliptischen Siliciumsäulen aus numerischen Simulationsergebnissen extrahiert. Dann wurde die Stromdichte J(r̂) = −iω[ε(r̂)−ε₀]E_inter(r̂) hergeleitet. Verschiedene Strommultipolmomente wurden zerlegt (Gleichung 4). Erste- und zweitordnungs Momente (Dipol und Quadrupol) wurden berechnet; höhere Ordnungen sind vernachlässigbar. Daraus wurden Multipolkoeffizienten a_E(l), a_M(l) und Streuquerschnitte C_s (Gleichung 5) berechnet.

Verfügbarkeit von Daten und Materialien

Die während der Studie verwendeten und/oder analysierten Datensätze sind auf vernünftige Anfrage beim korrespondierenden Autor erhältlich.

Abkürzungen

SEM:

Rasterelektronenmikroskop

LCP:

Linkszirkular polarisiert

RCP:

Rechtszirkular polarisiert

SCS:

Streuquerschnitt

MD:

Magnetischer Dipol

ED:

Elektrischer Dipol

EQ:

Elektrisches Quadrupol

MQ:

Magnetisches Quadrupol

THz-TDS:

Terahertz-Zeitbereichsspektroskopie