Институт в фотографиях

Научные и прикладные результаты лаборатории 11

2021 г.

Эффективный метод решения нелинейного векторного уравнения Шредингера

Разработан численный метод решения задачи Коши для нелинейного векторного уравнения Шредингера (модель Манакова), учитывающего, наряду с дисперсией и нелинейностью, еще и поляризацию волн. Метод основан на новых векторных алгоритмах решения обратной и прямой Задач Рассеяния (ЗР) для системы Манакова. Новая структура 4-блочных матриц, с векторными матрицами в недиагональных блоках, позволила обобщить известные скалярные алгоритмы. Обратная ЗР сводится к обращению алгоритмом типа Левинсона системы вложенных блочных матриц дискретизованных интегральных уравнений Гельфанда–Левитана–Марченко. Теплицева симметрия системы резко ускоряет расчеты: алгоритм требует всего O(N2) арифметических операций, где N – размер сетки. Обращение шагов алгоритма решения обратной задачи решает прямую ЗР. Численные тесты подтвердили эффективность новых векторных алгоритмов. Примером решения задачи Коши служит расчет столкновения ортогонально поляризованных солитонов Манакова (Рис. 1). Система Манакова возникает также при описании движения ультракоротких поляризованных оптические импульсов в резонансной среде. Новый метод пригоден и для этой важной задачи нелинейной оптики.

Пространственное распределение интенсивности поля в результате столкновения двух поляризационно-ортогональных солитонов

Рис.1. Расчет столкновения двух ортогонально поляризованных солитонов Манакова.

Frumin L.L. Algorithms for solving scattering problems for the Manakov model of nonlinear Schrödinger equations // Journal of Inverse and Ill-posed Problems – 2021, V.29, No.3, pp.369-383. [ DOI ]


 

2019 г.

Теория возмущений в задаче рассеяния света на диэлектрических телах

(ИАиЭ СО РАН совместно с ИФП СО РАН)

Можно построить теорию возмущений, считая поле в свободном пространстве примерно равным невозмущенному полю падающей волны.  Это приближение работает, если показатель преломления объекта близок к показателю среды. Такая ситуация типична, например, для биологических объектов в воде. Аналогичное приближение используется в квантовой механике для построения борновского ряда. Однако квантовомеханический ряд не дает корректного решения уравнений Максвелла для тел с резкой границей, потому что не соблюдаются граничные условия (скачки производной поля).

Построен ряд теории возмущений для задачи рассеяния электромагнитной волны на параллельных цилиндрах. Выведена система интегральных уравнений, точно учитывающая граничные условия на поверхности. Интегральные уравнения обобщены на случай нескольких тел произвольного сечения. Получен борновский ряд для продольной и поперечной поляризаций волны. Установлено, что в случае одного цилиндра ряд сходится к известному аналитическому решению. Сравнение с численным расчетом для двух цилиндров показало, что в третьем порядке погрешность индикатрисы рассеяния меньше 1%, рис.1. Предложенный метод позволяет получать простые аналитические формулы для сечений рассеяния света на наночастицах, которые полезны для оптической диагностики наноструктур, в частности, для анализа зависимости интенсивности рассеяния от параметров волны и диэлектрической частицы.

Диаграмма рассеяния плоской э/м волны на двух цилиндрах

Рис.1. Угловая диаграмма рассеяния плоской э/м волны с λ = 1.5 мкм на двух одинаковых параллельных стеклянных цилиндрах с радиусами 0.1 мкм, зазор между цилиндрами 0.1 мкм, угол плоскости с осями цилиндров составляет 450 от волнового вектора (горизонта), расстояние от точки наблюдения до оси одного из цилиндров 3 мкм: 3-е борновское приближение (черная пунктирная), BEM (синяя сплошная), DDA (красная пунктирная).

  1. Bereza A.S., Nemykin A.V., Perminov S.V., Frumin L.L., Shapiro D.A. Light scattering by dielectric bodies in the Born approximation // Physical Review A. – 2017, V.95, p.063839.
  2. Немыкин А.В., Берёза А.С., Шапиро Д.А. Теория возмущений в задаче рассеяния света на диэлектрическом теле // Письма в ЖЭТФ – 2018, Т.108, №9. стр.613-617.
  3. Bereza A.S., Frumin L.L., Nemykin A.V., Perminov S.V., Shapiro D.A. Perturbation series for the scattering of electromagnetic waves by parallel cylinders // EPL (Europhys. Lett.) – 2019, V.127, p.20002.

 

2016 г.

Поле у входа в субволновую щель

Совместно с Физико-техническим институтом, Брауншвейг (Германия)

Недавно была предсказана светоиндуцированная сила между металлическими телами. Подготовка эксперимента требует расчета оптического поля между двумя проводящими плоскостями, который до сих пор не был выполнен. В данной работе задача решена для идеального металла. Найдена переходная кривая между хорошо известными предельными случаями широкой и узкой щели по сравнению с длиной волны. Показано, что кривая содержит осцилляции, вызванные рождением в щели новых волноводных мод.Для случая субволновой щели найдены простые асимптотические формулы. Теоретическая кривая хорошо согласуется с численным расчетом для золота, см. рис. 1. Исследование плазмонной силы притяжения открывает возможность практических приложений, в частности, в манипуляции микронными металлическими частицами и в технике микропереключателей.

Мощность поля в щели

Рис.1. Теоретическая переходная кривая для идеального металла с учетом 100 пространственных мод (сплошная кривая) и численный расчет программой JCMsuite для золотой пленки толщиной 100 нм на длине волны 1.5 мкм (точки). Совпадение кривых свидетельствует о применимости модели идеального металла.

David Shapiro, Daniel Nies, Oleg Belai, Matthias Wurm and Vladimir Nesterov. Optical field and attractive force at the subwavelength slit // Optics Express – 2016, V.24, No.14, pp.15977–15982.


 

2015 г.

Исследование новых форматов передачи данных с фазовой модуляцией

Предложен итерационный метод подбора оптимального входного алфавита и расчета пропускной способности канала оптической связи. Продемонстрировано преимущество нового метода по сравнению с классическим алгоритмом Аримото – Блэхута.  Установлено, что новый метод особенно эффективен для каналов с дискретным входным алфавитом, неизвестным априори. На рис. 1 изображены вероятности, рассчитанные новым и традиционным методами при разных размерах интервала d (неопределенности положения символа). Из рисунка видно, что новый метод значительно быстрее сходится и одновременно обеспечивает существенно более высокую точность.

Вероятность символа входного алфавита

Рис.1. Сравнение нового метода с алгоритмом Аримото – Блэхута. Вероятность символа входного алфавита x=0.95±d в зависимости от десятичного логарифма числа итераций для разных интервалов:  d=0.3 (точки), 0.2 (штрихи), 0.1  (сплошная линия). Верхняя кривая (штрих-пунктир) – расчет итерационным методом с интервалом d=0.001.

Shapiro E.G., Shapiro D.A., and Turitsyn S.K. Method for computing the optimal signal distribution and channel capacity // Optics Express – 2015, V.23, No.12., pp.15119–15133.


 

 

2014 г.

Узкая оптическая фильтрация в фазово-модулированных линиях связи

Выполнено сравнение двух форматов передачи оптических битов со сдвигом фазы π и π/2 , соответствующих амплитудной и фазовой модуляции. Найдены выражения для дисперсии амплитуды импульсов, усредненной по возможным реализациям последовательности. Аналитически показано, что в присутствии узкого оптического фильтра π/2-схема обеспечивает меньшую дисперсию и более высокое качество сигнала. Результат объясняет эксперимент и подтверждается численным расчетом. Проанализированы причины искажения сигнала: неидеальный прямоугольный оптический фильтр на приемной стороне, ненайквистовская форма сигнала и сдвиг точки детектирования от центра частотной полосы. Найдены дисперсия и диаметр облака на констелляционной диаграмме. Показано, что при отношении сигнал/шум выше 11 дБ рассмотренные источники становятся главными.

pic1

pic2

Рис.1. Профиль мощности битовой последовательности 000 001 010 100 011 101 110 111 после фильтра (сплошная линия) и усредненный профиль (пунктир). Интервал фазовой модуляции π (а), π/2 (б). Из рисунка видно, что дисперсия в случае (б) заметно меньше.

  1. Shapiro E.G., Shapiro D.A.  Narrow optical filtering in phase modulated communication links // Electronics Letters – 2013, V.49, No.4, pp.279–281.
  2. Shapiro E.G., Shapiro D.A.  Interview // Electronics Letters – 2013, V.49, No.4, pp.279–281.
  3. Shapiro E.G., Shapiro D.A. Imperfect narrow filtering in optical links with phase modulation // Optics Communications – 2014, V.320, pp.27–32.

 

2013 г.

Поверхностный плазмон в периодической структуре нанопроволок

Обнаружены необычные плазмоны, возбуждаемые неоднородной волной в узкой щели в периодическом массиве нанопроволок, лежащем на диэлектрической подложке. Неоднородная волна появляется, когда исходная плоская волна падает на границу подложки под углом полного внутреннего отражения. Обнаружено, что плазмонный резонанс имеет резкую зависимость от угла падения. Это свойство делает резонанс перспективным для приложений в «умных» перестраиваемых устройствах плазмоники.

На рис. 1 показано распределение амплитуды магнитного поля в окрестности двух соседних золотых цилиндров диаметром 100 нм, расположенным над стеклом. Электрическое поле максимально в центре узкой щели размером 5 нм, где магнитное имеет наибольший градиент. На рис.2 показан коэффициент усиления локального поля в зависимости от угла падения для золота на разных длинах волн. Видно резкое изменение поля вблизи угла полного внутреннего отражения 42° для излучения с длиной волны 700–900 нм

cylinder field

Рис. 1. Распределение ближнего поля в окрестности цилиндров. Цветом показана амплитуда магнитного поля

filed amplification

Рис. 2. Коэффициент усиления локального поля в щели при λ = 886 (сплошная линия), 729 (частые штрихи), 653 (точки), 564 (редкие штрихи), 539 нм (штрих-пунктир).

  1. Frumin L.L., Nemykin A.V., Perminov S.V., Shapiro D.A. Plasmons excited by an evanescent wave in a periodic array of nanowires // J. Opt. – 2013, V.15, p.085002, 7 pages.
  2. Frumin L.L., Perminov S.V., Shapiro D.A. Plasmons excited by an evanescent wave // JOSA B – 2013, V.30, No.8, pp.2048–2052.