Разложение функции Грина для множественных связанных оптических резонаторов с конечной задержкой с использованием квазинормальных мод

Основная информация

• ID статьи: 2510.12511
• Название: Green's function expansion for multiple coupled optical resonators with finite retardation using quasinormal modes
• Авторы: Robert Fuchs¹, Juanjuan Ren², Stephen Hughes², Marten Richter¹
• Учреждения: ¹Technische Universität Berlin, ²Queen's University
• Классификация: cond-mat.mes-hall physics.optics quant-ph
• Дата публикации: 15 октября 2025
• Ссылка на статью: https://arxiv.org/abs/2510.12511

Аннотация

Электромагнитная функция Грина является ключевым элементом теории современных фотонных квантовых устройств, однако её прямое вычисление часто затруднено или невозможно. В данной работе предложена численно эффективная схема для расчёта рассеянной электромагнитной функции Грина многорезонаторных систем с пространственно разделёнными открытыми полостями (произвольной формы, дисперсионными и потерями) и конечным временем задержки. Схема основана на уравнении рассеяния Дайсона и позволяет конструировать функцию Грина из квазинормальных мод отдельных резонаторов в приближении малого числа мод и конечного числа итераций без вложенных интегралов. Метод демонстрирует превосходное согласие с полной численной функцией Грина на примере двух связанных диполей, расположенных в зазоре металлического димера как резонатора с квазинормальными модами, и легко расширяется на произвольно большие расстояния разделения и множественные полости.

Исследовательский контекст и мотивация

Проблемный фон

1. Основная проблема: Электромагнитная функция Грина имеет критическое значение в теории современных фотонных квантовых устройств, однако для многорезонаторных связанных систем прямое вычисление функции Грина чрезвычайно сложно или невозможно.
2. Важность: Функция Грина широко применяется в:
   • Расчётах усиления Пурселла оптических полостей
   • Квантовании электромагнитного поля
   • Выводе элементов связи фотонов с квантовыми излучателями
   • Приложениях квантовой оптики и квантовой динамики
3. Ограничения существующих методов:
   • Метод КНМ для одной полости: Для положений вдали от резонатора КНМ расходятся из-за комплексных собственных частот, требуя большого числа мод для хорошего приближения
   • Теория связанных квазинормальных мод (CQT): Зависит от расходящихся КНМ, недостаточно точна для больших расстояний разделения
   • Точное вычисление КНМ для многих полостей: Часто невозможно для множественных пространственно разделённых связанных структур
4. Исследовательская мотивация: Необходимость метода, способного:
   • Обрабатывать пространственно разделённые открытые полости произвольной формы, дисперсионные и с потерями
   • Включать эффекты конечной задержки
   • Эффективно вычислять в приближении малого числа мод
   • Применяться к многорезонаторным системам с большим пространственным разделением

Основные вклады

1. Предложена схема расчёта функции Грина многорезонаторной системы на основе уравнения рассеяния Дайсона, позволяющая конструировать функцию Грина многорезонаторной системы из квазинормальных мод одиночной полости
2. Разработан итеративный метод конструирования, последовательно добавляющий полости через конечное число итерационных шагов, избегая вычислительной сложности вложенных интегралов
3. Реализовано применение нормализованного поля КНМ, делающее метод применимым к системам с большим пространственным разделением и значительными эффектами задержки
4. Верифицирована точность метода, демонстрирующая превосходное согласие с полной численной функцией Грина в связанной системе металлических димеров
5. Предоставлена масштабируемая теоретическая схема, естественно разлагающая многорезонаторное рассеяние в произведение двухрезонаторных процессов рассеяния

Подробное описание метода

Определение задачи

Входные данные: Система N пространственно разделённых открытых полостей, каждая с произвольной формой, дисперсионными свойствами и потерями Выходные данные: Полная электромагнитная функция Грина многорезонаторной системы G^(N)(r,r',ω) Ограничения: В приближении малого числа мод, с учётом эффектов конечной задержки, избегая вычисления сложных вложенных интегралов

Теоретическая схема

1. Описание системы

Полная диэлектрическая проницаемость многорезонаторной системы представляется как:

ε(r,ω) = ε_back(r,ω) + Σ(i=1 to N) V_i(r,ω)

где V_i(r,ω) — возмущение диэлектрической проницаемости, вызванное i-й полостью.

2. Квазинормальные моды одиночной полости

Квазинормальные моды одиночной полости f̃_iμ(r) удовлетворяют безисточниковому уравнению Гельмгольца:

∇×∇×f̃_iμ(r) - (ω̃²_iμ/c²)ε_i(r,ω̃_iμ)f̃_iμ(r) = 0

с комплексными собственными частотами ω̃_iμ = ω_iμ - iγ_iμ, где γ_iμ > 0.

3. Нормализованное поле КНМ

Для решения проблемы расходимости КНМ в дальнем поле вводится нормализованное поле КНМ:

F̃_iμ(r,ω) = ∫_V_i d³r' Δε(r',ω)G_back(r,r',ω)·f̃_iμ(r')

4. Итеративное конструирование функции Грина многорезонаторной системы

Центральное инновационное решение — итеративное конструирование через серию уравнений рассеяния:

G^(n)(r,r',ω) = G^(n-1)(r,r',ω) + ∫_V_n d³s Δε(s,ω)G^(n-1)(r,s,ω)·G^(n)(s,r',ω)

Условие завершения:

G^(k)(r,r',ω)|_(r,r'∈V_i) ≈ G_i(r,r',ω)

5. Представление поверхностным интегралом

Для повышения численной эффективности объёмные интегралы преобразуются в поверхностные интегралы:

ω²/c² G^(n)(r,r')|_(r∉V_n) = χ_V_n(r') ω²/c² G^(n-1)(r,r') + [члены поверхностного интеграла]

Технические инновационные решения

1. Стратегия разложения: Многорезонаторное рассеяние естественно разлагается в произведение двухрезонаторных процессов рассеяния, избегая вложенных интегралов
2. Обработка эффектов задержки: Конечное время задержки правильно включается через фазовый множитель e^(iωR_ij/c)
3. Приближение полюсов: Приближение полюсов ω→ω_i применяется к медленно меняющейся огибающей функции для повышения вычислительной эффективности
4. Оптимизация поверхностного интеграла: Использование поверхностного интеграла избегает проблем с разрывом электрического поля на границах диэлектрика

Экспериментальная установка

Конфигурация системы

• Тестовая система: Два металлических димера как резонаторы с квазинормальными модами
• Модель материала: Модель Друде, ε_j(ω) = 1 - ω_p²/(ω² + iωγ_pj)
• Параметры установки:
  • Плазменная частота: ℏω_p = 8.2934 эВ
  • Параметры затухания: ℏγ_p1 = 0.0928 эВ, ℏγ_p2 = 0.3ℏγ_p1
• Собственные частоты КНМ:
  • Димер 1: ℏω̃_1 = (1.6904 - i0.0652) эВ (λ_1 ≈ 733.46 нм)
  • Димер 2: ℏω̃_2 = (1.6482 - i0.0388) эВ (λ_2 ≈ 752.24 нм)

Конфигурация тестирования

• Установка диполей: Два идентичных z-поляризованных диполя d_a = d_b = dê_z
• Положение: Расположены в зазорах двух димеров соответственно (r_a ∈ V_1, r_b ∈ V_2)
• Расстояние разделения: R_12 = 2020 нм и 760 нм в двух конфигурациях

Показатели оценки

Сила связи диполей:

g_ba(ω) = (1/ε_0) d_b · G(r_b, r_a, ω) · d_a

Методы сравнения

1. Полная численная функция Грина (G_full): Строгое решение трёхмерных уравнений Максвелла
2. Разложение КНМ (G_QNM): Метод данной работы, формулы (12)+(14)
3. Приближение полюсов (G_QNM^pole): Формулы (12)+(15), фаза включена в приближение полюсов
4. Разложение одиночной полости (G_1): Традиционный метод, рассматривающий только КНМ одиночной полости

Результаты экспериментов

Основные результаты

1. Случай большого расстояния разделения (R_12 = 2020 нм)

• Метод разложения КНМ демонстрирует превосходное согласие с полным численным решением
• Правильная обработка эффектов задержки (e^(iωR_12/c)) критична для точности
• Метод приближения полюсов (G_QNM^pole) отклоняется от численного решения из-за неполного учёта эффектов задержки
• Разложение одиночной полости серьёзно недооценивает силу связи, не может захватить эффект связи КНМ между полостями

2. Случай среднего расстояния разделения (R_12 = 760 нм)

• Метод разложения КНМ по-прежнему идеально совпадает с численным решением
• Метод приближения полюсов показывает лучшие результаты из-за более слабых эффектов задержки
• Разложение одиночной полости по-прежнему неэффективно как количественно, так и качественно

3. Ключевые численные результаты

Для R_12 = 2020 нм: N_21 = (2.0694 - i0.1357) × 10^(-7) нм^(-2) Для R_12 = 760 нм: N_21 = (5.3494 + i0.2984) × 10^(-7) нм^(-2)

Верификация метода

• Отсутствие параметров подгонки: Все расчёты основаны на первых принципах без подстройки параметров
• Частотная зависимость: Правильно захватывает частотную зависимость силы связи
• Фазовые эффекты: Фаза задержки e^(iωR_12/c) критична для точности при больших расстояниях разделения

Вычислительная эффективность

• Избежание вложенных интегралов: Многорезонаторное рассеяние разлагается в произведение двухрезонаторных процессов
• Конечное число итераций: Для N полостей требуется N итераций
• Квадратичное масштабирование: В приближении одной моды количество интегралов рассеяния масштабируется квадратично с числом полостей

Связанные работы

Основные направления исследований

1. Теория квазинормальных мод: Теория спонтанного излучения открытых полостей Сована и др.
2. Квантование поля: Методы квантования электромагнитного поля на основе КНМ
3. Теория связанных полостей: Теория связанных квазинормальных мод (CQT) и её ограничения
4. Методы функции Грина: Разложение КНМ для одиночной полости и методы нормализации

Преимущества данной работы

1. Обработка эффектов задержки: По сравнению с CQT, может точно обрабатывать случаи больших расстояний разделения
2. Вычислительная эффективность: Избегает проблемы вложенных интегралов традиционных методов
3. Универсальность: Применима к полостям произвольной формы, дисперсионным и с потерями
4. Масштабируемость: Естественно расширяется на многорезонаторные системы

Заключение и обсуждение

Основные выводы

1. Успешно разработана схема расчёта функции Грина многорезонаторной системы на основе уравнения рассеяния Дайсона и квазинормальных мод одиночной полости
2. Достигнута высокая точность расчётов, демонстрирующая превосходное согласие с численным решением на различных расстояниях разделения
3. Решена проблема эффектов задержки, делающая метод применимым к системам с большим пространственным разделением
4. Предоставлено масштабируемое решение, легко обобщаемое на произвольное число полостей

Ограничения

1. Приближение малого числа мод: Метод основан на предположении, что каждая полость описывается малым числом доминирующих КНМ
2. Условие завершения: Зависит от предположения, что функция Грина внутри полости хорошо аппроксимируется одиночными КНМ полости
3. Трёхмерная однородная среда: Основная верификация проведена для трёхмерных полостей в однородной фоновой среде
4. Вычислительная сложность: Хотя вложенные интегралы избегаются, всё ещё требуется вычисление множественных поверхностных интегралов

Будущие направления

1. Расширение на системы с волноводной связью: Адаптация к различным граничным условиям
2. Нелинейные эффекты: Включение нелинейных оптических процессов
3. Приложения квантовой динамики: Конкретные приложения в системах квантовой оптики
4. Оптимизация алгоритмов: Дальнейшее повышение вычислительной эффективности для больших многорезонаторных систем

Глубокая оценка

Достоинства

1. Сильная теоретическая новизна: Итеративное применение уравнения рассеяния Дайсона является новым и остроумным
2. Высокая практическая ценность: Решает реальные вычислительные потребности многорезонаторных квантово-оптических систем
3. Достаточная верификация: Сравнение со строгим численным решением подтверждает точность метода
4. Ясное изложение: Строгий теоретический вывод, чёткая физическая картина

Недостатки

1. Отсутствие экспериментальной верификации: Только теоретическая и численная верификация, отсутствует экспериментальная проверка
2. Ограниченный диапазон применения: Основное тестирование проведено на специфических типах резонаторных систем
3. Применимость к сложным системам: Применимость к сильно связанным или нелинейным системам остаётся неясной
4. Детали численной реализации: Некоторые детали численной реализации могут влиять на воспроизводимость метода

Влияние

1. Научный вклад: Предоставляет важный инструмент для теории многорезонаторной квантовой оптики
2. Технологические приложения: Потенциальные приложения в квантовых сетях, квантовых вычислениях и других областях
3. Методологическое значение: Демонстрирует потенциал применения теории КНМ в сложных системах
4. Воспроизводимость: Чёткая теоретическая схема облегчает применение и расширение метода другими исследователями

Сценарии применения

1. Проектирование квантово-оптических устройств: Теоретический анализ многорезонаторных связанных систем
2. Плазмоника: Расчёт электромагнитного отклика металлических наноструктур
3. Квантовые сети: Анализ связи пространственно разделённых квантовых узлов
4. Фотонные кристаллы: Исследование систем с дефектными полостями

Библиография

Статья цитирует 47 важных справочных источников, охватывающих теорию квазинормальных мод, квантовую оптику, квантование электромагнитного поля и другие смежные области, обеспечивая прочную теоретическую основу для данного исследования.

Общая оценка: Это высококачественная теоретическая физическая статья, предлагающая инновационный метод решения важной проблемы расчёта функции Грина для многорезонаторных оптических систем. Теоретический вывод строг, численная верификация полна, работа имеет важное значение для области квантовой оптики и фотоники. Универсальность и масштабируемость метода обеспечивают широкие перспективы его применения.