Guía de onda plasmónica de terahercios ultracompacta con confinamiento mejorado mediante estructura de semimetal de Dirac en masa-aislante-metal

Resumen

Se investiga una guía de onda plasmónica de terahercios sublongitud de onda construida sobre una pila de semimetal de Dirac en masa (BDS)-aislante-metal (BIM). La estructura entrega una banda de frecuencia optimizada donde el modo polaritón de plasmon superficial (SPP) está fuertemente confinado —hasta factor de confinamiento de ½·λ₀/15— mientras la pérdida de propagación permanece baja a 1.0 dB/λ₀. La introducción de dos cintas de silicio en la guía BIM crea un resonador Fabry-Pérot dinámicamente sintonizable que puede moldear SPPs de terahercios en escala profunda sublongitud de onda, ofreciendo vía hacia dispositivos plasmónicos THz ultracompactos, sintonizables y filtros ópticos.

Antecedentes

La radiación terahercios (THz) ha emergido como herramienta versátil en imagen, sensado bioquímico y comunicaciones inalámbricas —demandando confinamiento cada vez más estrecho, mayor resolución e integración mayor — para aprovechar su pleno potencial. Lograr confinamiento sublongitud de onda de ondas THz es crítico, sin embargo, los polaritones de plasmon superficial (SPPs) convencionales en metales nobles sufren confinamiento pobre, pérdida intrínseca alta y sintonización pasiva —limitaciones que han obstaculizado aplicaciones prácticas.

Los plasmons de grafeno han mostrado promesa debido a su baja pérdida, sintonización dinámica y confinamiento extremo, permitiendo dispositivos THz de alta resolución —ejemplos incluyen heteroestructuras de grafeno sintonizables por compuerta, efectos Talbot discretos y guías de onda de transparencia inducida por plasmon. Estos avances subrayan el papel crucial del confinamiento extremo de SPPs para manipular radiación THz a escalas profunda sublongitud de onda.

Los semimetales de Dirac en masa (BDS), a menudo llamados “grafeno 3-D”, exhiben movilidades de portadores ultrahigh (hasta 9×10⁶ cm² V⁻¹ s⁻¹) que superan con creces al mejor grafeno (2×10⁵ cm² V⁻¹ s⁻¹) —ventaja clave para plasmons de baja pérdida. Su respuesta dieléctrica también es sintonizable vía modulación de energía de Fermi, y materiales como Na₃Bi, Cd₃As₂ y cuasicristales AlCuFe son más fáciles de procesar y más estables que grafeno —posicionando BDS como material plasmónico prometedor de próxima generación. Sin embargo, el confinamiento del modo SPP en una interfaz BDS-aislante por sí sola sigue siendo insuficiente para manipulación profunda sublongitud de onda.

En este estudio, presentamos una guía BIM que mejora marcadamente el confinamiento SPP manteniendo baja pérdida y permitiendo sintonización dinámica. La dispersión, características de pérdida y capacidades de filtrado de este modo altamente confinado se investigan sistemáticamente, revelando una ventana de frecuencia optimizada con confinamiento mejorado y pérdida de propagación reducida —comportamiento raro en SPPs basados en metales tradicionales. Además, el modo BIM puede propagarse a través de rendijas más estrechas que λ₀/2000, y al incrustar cintas de silicio como espejos, realizamos un filtro pasa-banda sintonizable cuya resonancia puede desplazarse ajustando la energía de Fermi del BDS.

Teoría y Simulación

La guía BIM consiste en una película monocapa de BDS (espesor 0.2 µm) separada de un sustrato de plata por un espaciador dieléctrico (permitividad εᵣ). En régimen THz, la plata se comporta como conductor eléctrico perfecto. Para excitación TM-polarizada, el modo SPP se propaga junto a x con vector de onda k_SPP y decae junto a y en espacio libre. Aplicando condiciones de frontera se obtiene la relación de dispersión:

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

donde k₀ es el vector de onda en espacio libre y σ es la conductividad compleja de la película BDS (Ecs. 3–4). Resolviendo Ec. (1) se obtiene el índice de refracción efectivo n_eff = k_SPP/k₀ = Re(n_eff) + i Im(n_eff). Re(n_eff) cuantifica confinamiento del modo, mientras Im(n_eff) se relaciona directamente con pérdida de propagación.

Para huecos suficientemente grandes (tanh[…]≈1), Ec. (1) se simplifica a dispersión de monocapa BDS (Ec. 2), destacando el beneficio añadido del sustrato metálico en configuración BIM.

Resultados y Discusión

Primero examinamos cómo el hueco BIM g y la energía de Fermi del BDS E_F afectan confinamiento y pérdida. Con E_F=70 meV, los índices de refracción efectivos para varios valores g se grafican en Fig. 2a (Re(n_eff)) y Fig. 2b (Im(n_eff)). Para g≥10 µm, las curvas convergen arriba de 0.05 THz, indicando que el campo SPP ya está confinado dentro de ~10 µm y la contribución de plata es despreciable. Reducir g por debajo de 1 µm aumenta marcadamente confinamiento y reduce pérdida, produciendo una banda de frecuencia optimizada donde el confinamiento alcanza pico mientras la pérdida disminuye —rasgo poco común en SPPs basados en metales convencionales.

Figura 2c,d ilustran dependencia en E_F (g=1 µm). Aumentar E_F reduce tanto confinamiento como pérdida, atribuible a mayor metallicidad y tiempos de relajación de portadores más largos. Por ejemplo, a 2.5 THz, confinamiento de λ₀/15 y pérdida de 1.0 dB/λ₀ se logran con g=10 nm y E_F=70 meV.

Simulaciones de transmisión (Fig. 3) confirman pérdida reducida: a 1.56 THz, la guía BIM exhibe transmisión de 0.97, superando caso de monocapa BDS. Las correspondientes distribuciones de campo magnético (Fig. 3b,c) muestran Re(n_eff)=2.45 para BIM versus 1.002 para monocapa, con campo fuertemente localizado dentro de rendijas <λ₀/2000.

Para aprovechar el modo altamente confinado, integramos dos cintas de silicio (n_Si=3.4) en el espaciador dieléctrico, formando espejos reflectivos (Fig. 4a). La cavidad Fabry-Pérot resultante soporta resonancias de onda estacionaria, produciendo dos picos de transmisión a 1.56 y 2.22 THz con FWHM de 0.12 y 0.09 THz, respectivamente. Mapas de campo magnético (Fig. 4c,e) confirman modos de cavidad de primer y segundo orden. Fuera de resonancia, ondas estacionarias no pueden formarse y la entrada se refleja (Fig. 4d).

La transmisión puede describirse analíticamente por teoría de modos acoplados (CMT):

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

donde κ_w y κ_i son tasas de acoplamiento de guía y pérdida intrínseca. Los factores de calidad Q_t=ω₀/FWHM y Q_oi=−Re(n_eff)/(2Im(n_eff)) permiten extraer κ_w y κ_i, con resultados analíticos coincidiendo con espectros numéricos (Fig. 4b).

Figura 5 explora dependencia de frecuencia resonante en longitud de cavidad L. Primera y segunda resonancias se desplazan al rojo con L creciente, consistente con condición de onda estacionaria 2k_SPP(ω_r)L+θ=2mπ. Similarmente, aumentar g desplaza resonancias al azul (Fig. 6a) debido a Re(n_eff) reducido, mientras aumentar E_F también desplaza picos al azul (Fig. 6b) y estrecha ancho de banda, reflejando menor Im(n_eff). Estos mecanismos de sintonización permiten control dinámico del filtro sin rediseñar estructura.

Conclusiones

La guía BIM entrega modo plasmónico de terahercios ultracompacto con confinamiento λ₀/15 y pérdida 1.0 dB/λ₀, alcanzable en rendijas más estrechas que λ₀/2000. Al incorporar cintas de silicio, realizamos un filtro pasa-banda sintonizable dinámicamente cuya frecuencia resonante se desplaza con la energía de Fermi del BDS, permitiendo conmutación y filtrado THz sin re-optimización estructural.

Métodos

Simulaciones numéricas emplearon solver 2-D de diferencias finitas en dominio temporal (FDTD) con capas perfectamente acopladas en direcciones x e y. Malla BDS se estableció a 1 µm × 0.02 µm para convergencia. Conductividad dependiente de frecuencia de BDS sigue fórmula de Kubo bajo aproximación de fase aleatoria —Ecs. (3) y (4) describen Reσ(Ω) y Imσ(Ω). Parámetros para AlCuFe: t=40, ε_c=3, μ=3×10⁴ cm² V⁻¹ s⁻¹, E_F=70 meV.

No participaron humanos, datos ni animales en esta investigación.

Abreviaturas

BDS

Semimetales de Dirac en masa

BIM

BDS-aislante-metal

CMT

Teoría de modos acoplados

FDTD

Diferencias finitas en dominio temporal

FWHM

Ancho a media altura

SPPs

Polaritones de plasmon superficial