Guia de Ondas Plasmonico Terahertz Ultracompacto com Confinamento Aprimorado usando Estrutura Semicondutor de Dirac Bulk-Isolante-Metal

Resumo

Um guia de ondas plasmonico terahertz subcomprimento de onda construído sobre uma pilha semicondutor de Dirac bulk (BDS)–isolante–metal (BIM) é investigado. A estrutura fornece banda de frequência otimizada onde o modo polariton de plasmon de superfície (SPP) é fortemente confinado — até fator de confinamento de ½·λ₀/15 — enquanto a perda de propagação permanece baixa em 1,0 dB/λ₀. Introduzindo duas fitas de silício no guia BIM cria um ressonador Fabry–Pérot dinamicamente sintonizável que pode moldar SPPs terahertz em escala deep-subcomprimento de onda, oferecendo caminho para dispositivos plasmonicos THz ultracompactos, sintonizáveis e filtros ópticos.

Introdução

Radiação terahertz (THz) emergiu como ferramenta versátil em imagem, sensoriamento bioquímico e comunicações sem fio — demandando confinamento cada vez mais apertado, maior resolução e maior integração — para aproveitar seu potencial total — [1,2,3]. Alcançar confinamento subcomprimento de onda de ondas THz é crítico, mas plasmons de superfície polaritons (SPPs) em metais nobres sofrem de confinamento fraco, perda intrínseca alta e sintonizabilidade passiva — limitações que dificultaram aplicações práticas — [4,5,6].

Plasmons de grafeno mostraram promessa devido a baixa perda, sintonizabilidade dinâmica e confinamento extremo, permitindo dispositivos THz de alta resolução — exemplos incluem heteroestruturas de grafeno sintonizáveis por porta, efeitos Talbot discretos e guias de onda de transparência induzida por plasmon — [11–18]. Estes avanços sublinham papel crucial do confinamento extremo de SPP para manipular radiação THz em escalas deep-subcomprimento de onda.

Semicondutores de Dirac bulk (BDS), frequentemente chamados "grafeno 3-D", exibem mobilidades de portadores ultra-altas (até 9×10⁶ cm² V⁻¹ s⁻¹) que excedem em muito o melhor grafeno (2×10⁵ cm² V⁻¹ s⁻¹) — vantagem chave para plasmons de baixa perda — [19]. Sua resposta dielétrica também é sintonizável via modulação de energia de Fermi, e materiais como Na₃Bi, Cd₃As₂ e quasicristais AlCuFe são mais fáceis de processar e mais estáveis que grafeno — posicionando BDS como material plasmonico promissor de próxima geração — [20,21]. No entanto, confinamento de modo SPP em interface BDS-isolante sozinho ainda é insuficiente para manipulação deep-subcomprimento de onda.

Neste estudo, apresentamos guia de ondas BIM que melhora marcadamente confinamento de SPP mantendo baixa perda e permitindo sintonizabilidade dinâmica. Dispersão, características de perda e capacidades de filtragem deste modo altamente confinado são investigadas sistematicamente, revelando janela de frequência otimizada com confinamento aprimorado e perda de propagação reduzida — comportamento raramente observado em SPPs baseados em metal tradicionais. Além disso, modo BIM pode propagar através de fendas mais estreitas que λ₀/2000, e incorporando fitas de silício como espelhos, realizamos filtro passa-banda sintonizável cuja ressonância pode ser deslocada ajustando energia de Fermi do BDS.

Teoria e Simulação

O guia de ondas BIM consiste em filme monocamada BDS (espessura 0,2 µm) separado de substrato de prata por espaçador dielétrico (permissividade εᵣ). No regime THz, prata comporta-se como condutor elétrico perfeito. Para excitação TM-polarizada, modo SPP propaga ao longo de x com vetor de onda k_SPP e decai ao longo de y no espaço livre. Aplicando condições de contorno obtém-se relação de dispersão:

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

onde k₀ é vetor de onda no espaço livre e σ é condutividade complexa do filme BDS (Eq. 3–4). Resolvendo Eq. (1) obtém-se índice de refração efetivo n_eff = k_SPP/k₀ = Re(n_eff) + i Im(n_eff). Re(n_eff) quantifica confinamento do modo, enquanto Im(n_eff) relaciona-se diretamente com perda de propagação.

Para gaps suficientemente grandes (tanh[…]≈1), Eq. (1) simplifica para dispersão de camada única BDS (Eq. 2), destacando benefício adicional do substrato metálico na configuração BIM.

Esquemático do guia de ondas plasmonico BIM: filme monocamada BDS posicionado a gap g de substrato de prata, separado por espaçador dielétrico com permissividade εᵣ. Modo SPP TM-polarizado propaga ao longo de x e decai ao longo de y. Distribuição de campo Eₓ mostrada pela linha vermelha.

Resultados e Discussão

Primeiro examinamos como gap BIM g e energia de Fermi E_F do BDS afetam confinamento e perda. Com E_F=70 meV, índices de refração efetivos para vários valores de g são plotados na Fig. 2a (Re(n_eff)) e Fig. 2b (Im(n_eff)). Para g≥10 µm, curvas convergem acima de 0,05 THz, indicando que campo SPP já está confinado dentro de ~10 µm e contribuição da prata é insignificante. Reduzir g abaixo de 1 µm aumenta marcadamente confinamento e reduz perda, gerando banda de frequência otimizada onde confinamento atinge pico enquanto perda diminui — característica incomum em SPPs baseados em metal convencionais.

Figura 2c,d ilustram dependência em E_F (g=1 µm). Aumentar E_F reduz tanto confinamento quanto perda, atribuível a metalicidade aprimorada e tempos de relaxamento de portadores mais longos. Por exemplo, a 2,5 THz, confinamento de λ₀/15 e perda de 1,0 dB/λ₀ são alcançados com g=10 nm e E_F=70 meV.

Re(n_eff) e Im(n_eff) versus (a) gap g e (b) energia de Fermi E_F para E_F=70 meV (g variável) e g=1 µm (E_F variável).

Simulações de transmissão (Fig. 3) confirmam perda reduzida: a 1,56 THz, guia BIM exibe transmissão de 0,97, superando caso BDS monocamada. Distribuições de campo magnético correspondentes (Fig. 3b,c) mostram Re(n_eff)=2,45 para BIM versus 1,002 para monocamada, com campo firmemente localizado dentro de fendas <λ₀/2000.

Espectros de transmissão (a) e distribuições de campo magnético (b,c) para guias BIM (vermelho) e monocamada (azul) a E_F=70 meV, g=50 µm e 1,56 THz.

Para aproveitar modo altamente confinado, integramos duas fitas de silício (n_Si=3,4) no espaçador dielétrico, formando espelhos reflexivos (Fig. 4a). Cavidade Fabry–Pérot resultante suporta ressonâncias de onda estacionária, gerando dois picos de transmissão a 1,56 e 2,22 THz com FWHM de 0,12 e 0,09 THz, respectivamente. Mapas de campo magnético (Fig. 4c,e) confirmam modos de cavidade de primeira e segunda ordem. Fora de ressonância, ondas estacionárias não se formam, e entrada é refletida (Fig. 4d).

Transmissão pode ser descrita analiticamente por teoria de modos acoplados (CMT):

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

onde κ_w e κ_i são taxas de acoplamento do guia e perda intrínseca. Fatores de qualidade Q_t=ω₀/FWHM e Q_oi=−Re(n_eff)/(2Im(n_eff)) permitem extração de κ_w e κ_i, com resultados analíticos correspondendo a espectros numéricos (Fig. 4b).

(a) Guia BIM com fitas de silício (largura d, separação L). (b) Espectros de transmissão numérico (azul) e ajustado CMT (vermelho) para g=1 µm, d=5 µm, L=120 µm. (c–e) Distribuições de campo magnético a 1,56, 1,90 e 2,22 THz.

Figura 5 explora dependência de frequência ressonante em comprimento de cavidade L. Primeira e segunda ressonâncias deslocam para vermelho com aumento de L, consistente com condição de onda estacionária 2k_SPP(ω_r)L+θ=2mπ. Similarmente, aumentar g desloca ressonâncias para azul (Fig. 6a) devido a Re(n_eff) reduzido, enquanto aumentar E_F também desloca picos para azul (Fig. 6b) e estreita largura de banda, refletindo menor Im(n_eff). Esses mecanismos de sintonizabilidade permitem controle dinâmico do filtro sem redesenhar estrutura.

(a) Espectros de transmissão para L variável. (b) Frequências ressonantes dos modos 1 e 2 versus L (g=1 µm, d=5 µm, E_F=70 meV).

Espectros de transmissão versus (a) g e (b) E_F (outros parâmetros fixos). (c) Re(n_eff) vs. E_F e g. (d) Im(n_eff) vs. E_F e g.

Conclusões

Guia BIM entrega modo plasmonico terahertz ultracompacto com confinamento λ₀/15 e perda 1,0 dB/λ₀, alcançável em fendas mais estreitas que λ₀/2000. Incorporando fitas de silício, realizamos filtro passa-banda dinamicamente sintonizável cuja frequência ressonante desloca com energia de Fermi do BDS, permitindo comutação e filtragem THz sem re-otimização estrutural.

Métodos

Simulações numéricas empregaram solucionador 2-D de diferenças finitas no domínio do tempo (FDTD) com camadas perfeitamente casadas nas direções x e y. Malha BDS definida para 1 µm × 0,02 µm para convergência. Condutividade dependente de frequência de BDS segue fórmula de Kubo sob aproximação de fase aleatória – Eqs. (3) e (4) descrevem Reσ(Ω) e Imσ(Ω). Parâmetros para AlCuFe: t=40, ε_c=3, μ=3×10⁴ cm² V⁻¹ s⁻¹, e E_F=70 meV.

Nenhum participante humano, dados ou animais envolvidos nesta pesquisa.

Abreviações

BDS:

Semicondutores de Dirac bulk

BIM:

BDS-isolante-metal

CMT:

Teoria de modos acoplados

FDTD:

Diferenças finitas no domínio do tempo

FWHM:

Largura total a meia altura

SPPs:

Polaritons de plasmon de superfície