Ультракомпактный терагерцовый плазмонный волновод с усиленным замыканием на основе объемного дираковского полуметалла–изолятель–металл

Аннотация

Исследуется субволновой терагерцовый плазмонный волновод, построенный на стеке объемный дираковский полуметалл (BDS)–изолятель–металл (BIM). Структура предоставляет оптимизированный частотный диапазон, где режим поверхностного плазмонно-поляритона (SPP) тесно зажат — до фактора замыкания ½·λ₀/15 — в то время как потери распространения остаются низкими, 1,0 дБ/λ₀. Введение двух кремниевых лент в волновод BIM создает динамически настраиваемый Фабри-Перо резонатор, способный формировать терагерцовые SPP в глубоко субволновом масштабе, открывая путь к ультракомпактным, настраиваемым THz плазмонным устройствам и оптическим фильтрам.

Введение

Терагерцовое (THz) излучение стало универсальным инструментом в визуализации, биохимическом сенсиринге и беспроводных коммуникациях — требуя всё более тесного замыкания, более высокого разрешения и большей интеграции — для реализации его полного потенциала. Достижение субволнового замыкания THz волн критично, однако обычные SPP в благородных металлах страдают от слабого замыкания, высоких внутренних потерь и пассивной настраиваемости — ограничений, препятствующих практическим приложениям.

Графеновые плазмоны показали обещание благодаря низким потерям, динамической настраиваемости и экстремальному замыканию, позволяя высокоразрешающие THz устройства — примеры включают затворно-настраиваемые графеновые гетероструктуры, дискретные эффекты Тальбота и волноводы с плазмонно-индуцированной прозрачностью. Эти достижения подчеркивают ключевую роль экстремального SPP замыкания для манипуляции THz излучением в глубоко субволновых масштабах.

Объемные дираковские полуметаллы (BDS), часто называемые «3D графеном», демонстрируют ультравысокую подвижность носителей (до 9×10⁶ см²/В·с), далеко превосходящую лучший графен (2×10⁵ см²/В·с) — ключевое преимущество для низкопотерь плазмоны. Их диэлектрический отклик также настраиваем через модуляцию ферми-энергии, а материалы такие как Na₃Bi, Cd₃As₂ и квазикристаллы AlCuFe легче в обработке и стабильнее графена — позиционируя BDS как перспективный плазмонный материал следующего поколения. Однако замыкание SPP режима на интерфейсе BDS-изолятель в одиночку все еще недостаточно для глубоко субволновой манипуляции.

В этом исследовании представлен BIM волновод, который заметно улучшает SPP замыкание, сохраняя низкие потери и обеспечивая динамическую настраиваемость. Дисперсия, потерь и фильтрующие возможности этого высоко замыкающего режима систематически исследуются, раскрывая оптимизированное частотное окно с одновременно усиленным замыканием и сниженными потерями распространения — поведение, редко наблюдаемое в традиционных металлических SPP. Более того, BIM режим может распространяться через щели уже λ₀/2000, и встраивая кремниевые ленты как зеркала, мы реализуем настраиваемый полосовой фильтр, чей резонанс можно сдвигать настройкой ферми-энергии BDS.

Теория и моделирование

BIM волновод состоит из монослойной пленки BDS (толщина 0,2 мкм), отделенной серебряной подложкой диэлектрическим зазором (проницаемость εᵣ). В THz диапазоне серебро ведет себя как идеальный электрический проводник. Для TM-поляризованного возбуждения SPP режим распространяется вдоль x с волновым вектором k_SPP и затухает вдоль y в свободное пространство. Применение граничных условий дает дисперсионное соотношение:

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

где k₀ — волновой вектор в свободном пространстве, σ — комплексная проводимость пленки BDS (уравнения 3–4). Решение уравнения (1) дает эффективный показатель преломления n_eff = k_SPP/k₀ = Re(n_eff) + i Im(n_eff). Re(n_eff) kvantфицирует замыкание режима, а Im(n_eff) напрямую связано с потерями распространения.

Для достаточно больших зазоров (tanh[…]≈1) уравнение (1) упрощается до дисперсии однослойного BDS (уравнение 2), подчеркивая добавленную выгоду металлической подложки в BIM конфигурации.

Результаты и обсуждение

Сначала изучаем влияние зазора BIM g и ферми-энергии BDS E_F на замыкание и потери. При E_F=70 мэВ эффективные показатели преломления для различных g показаны на рис. 2a (Re(n_eff)) и 2b (Im(n_eff)). Для g≥10 мкм кривые сходятся выше 0,05 THz, указывая, что SPP поле уже зажато в ~10 мкм и вклад серебра пренебрежим. Уменьшение g ниже 1 мкм резко увеличивает замыкание и снижает потери, давая оптимизированный частотный диапазон, где замыкание пикует, пока потери уменьшаются — особенность, нетипичная для обычных металлических SPP.

Рис. 2c,d иллюстрируют зависимость от E_F (g=1 мкм). Рост E_F снижает и замыкание, и потери, благодаря усиленной металличности и более длинным временам релаксации носителей. Например, при 2,5 THz достигается замыкание λ₀/15 и потери 1,0 дБ/λ₀ при g=10 нм и E_F=70 мэВ.

Симуляции прохождения (рис. 3) подтверждают сниженные потери: при 1,56 THz BIM волновод демонстрирует прохождение 0,97, превосходя однослойный BDS. Соответствующие распределения магнитного поля (рис. 3b,c) показывают более высокий Re(n_eff)=2,45 для BIM против 1,002 для однослойного, с полем тесно локализованным в щелях <λ₀/2000.

Для использования высоко замыкающего режима интегрируем две кремниевые ленты (n_Si=3,4) в диэлектрический зазор, формируя отражающие зеркала (рис. 4a). Полученная Фабри-Перо полость поддерживает стоячие волны резонансы, давая два пика прохождения при 1,56 и 2,22 THz с FWHM 0,12 и 0,09 THz соответственно. Карты магнитного поля (рис. 4c,e) подтверждают первый и второй порядок полостных мод. Вне резонанса стоячие волны не формируются, и вход отражается (рис. 4d).

Прохождение можно аналитически описать теорией связанных мод (CMT):

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

где κ_w и κ_i — скорости связывания волновода и внутренние потери. Качественные факторы Q_t=ω₀/FWHM и Q_oi=−Re(n_eff)/(2Im(n_eff)) позволяют извлечь κ_w и κ_i, аналитические результаты совпадают с числовыми спектрами (рис. 4b).

Рис. 5 исследует зависимость резонансной частоты от длины полости L. Первый и второй резонансы красное смещаются с увеличением L, согласуясь с условием стоячей волны 2k_SPP(ω_r)L+θ=2mπ. Аналогично, увеличение g синее смещает резонансы (рис. 6a) из-за снижения Re(n_eff), тогда как увеличение E_F также синее смещает пики (рис. 6b) и сужает полосу, отражая более низкие Im(n_eff). Эти механизмы настраиваемости позволяют динамическое управление фильтром без перепроектирования структуры.

Заключение

BIM волновод доставляет ультракомпактный терагерцовый плазмонный режим с замыканием λ₀/15 и потерями 1,0 дБ/λ₀, достижимым в щелях уже λ₀/2000. Включая кремниевые ленты, мы реализуем динамически настраиваемый полосовой фильтр, чья резонансная частота сдвигается с ферми-энергией BDS, позволяя THz коммутацию и фильтрацию без структурной ре-оптимизации.

Методы

Численные симуляции использовали 2D FDTD решатель с идеально согласованными слоями в направлениях x и y. Сетка BDS установлена 1 мкм × 0,02 мкм для сходимости. Частотно-зависимая проводимость BDS следует Кубо формуле под случайной фазовой аппроксимацией — уравнения (3) и (4) описывают Reσ(Ω) и Imσ(Ω). Параметры для AlCuFe: t=40, ε_c=3, μ=3×10⁴ см²/В·с, E_F=70 мэВ.

В исследовании не участвовали люди, данные или животные.

Аббревиатуры

BDS

Объемные дираковские полуметаллы

BIM

BDS-изолятель-металл

CMT

Теория связанных мод

FDTD

Метод конечных разностей во временной области

FWHM

Полная ширина на полувысоте

SPPs

Поверхностные плазмонно-поляритоны