На симуляции метаповерхность сдвинула частоту света и повернула луч

Зеркало отражает свет по закону геометрии: угол падения задаёт угол выхода. А можно ли заставить тот же тонкий отражающий слой создать новые частоты и разложить их по нужным лучам? Авторы препринта arXiv:2608.09905 рассчитали именно такую ITO ENZ-метаповерхность. Две боковые полосы, появившиеся около исходной оптической частоты, в их модели преимущественно уходят в один и тот же дифракционный порядок, хотя под немного разными углами. Вопрос оказывается сложнее привычного расчёта неподвижной решётки: материал сам меняется, пока через него проходит поле.

Периодическая модуляция плазменной частоты ITO работает как насос частотного смесителя. Поле с несущей частотой ω0 раскладывается на Floquet-компоненты ωn = ω0 + nΩ. Модулируемая проводимость или диэлектрическая проницаемость связывает эти компоненты: несущая обменивается энергией с боковыми полосами, а пространственная решётка назначает каждой из них доступные направления выхода.

Для метаповерхности из трёх TiO2-элементов над ENZ-слоем ITO и золотым зеркалом авторы задали модуляцию на 3 ТГц. После численной оптимизации полосы n = −1 и n = +1 показали наибольшую амплитуду отражённого поля в порядке m = +1. Его расчётные углы составили 54.9° и 53.4°. Получается компактная оптическая структура, где частота и направление луча проектируются совместно: временная накачка создаёт канал, литографический рисунок выбирает траекторию.

Обычный частотный расчёт здесь легко даёт неверную картину. ITO около epsilon-near-zero режима, графеновые листы и резонансные метаструктуры резко меняют отклик с частотой. Боковую полосу нельзя заставить распространяться с параметрами несущей. Поэтому метод собирает все удержанные гармоники в одну матрицу harmonic balance, но для каждой подставляет её физическую частоту и собственную дисперсию материала. Внедиагональные блоки этой матрицы описывают частотное смешение, диагональные — распространение каждой отдельной полосы.

Авторы сначала проверили схему на диэлектрическом цилиндре с модулируемым графеновым покрытием. Для параметрической накачки при Ω = 2ω0 численный расчёт рассеяния и поглощения совпал с аналитическим Floquet-решением. Компонента n = −1 в этой конфигурации соответствует сопряжённой отрицательно-частотной части комплексного описания реального поля; отдельного света с физически отрицательной частотой здесь нет.

Это пока вычислительная работа. Образец не изготовили, измерений нет. Ещё существеннее то, что авторы сравнили амплитуды полей по дифракционным каналам, но не посчитали, какая доля оптической мощности действительно попадает в нужные лучи. Без этого нельзя судить о потерях в ITO и зеркале, паразитных порядках и цене временной накачки. Зато сама постановка даёт инженерный язык для устройства, которое одновременно ведёт свет по частоте и по углу, вместо того чтобы рассматривать эти задачи по отдельности.

Первоисточники: arXiv:2608.09905, «Full-Wave Harmonic Balance Framework for Dispersive Time-Varying Photonic Structures», arXiv:2608.09905. Полный текст с параметрами моделей и результатами оптимизации: arXiv:2608.09905.