1. Введение

Микроскопические светодиоды (µLED) с размерами около 1 мкм имеют критическое значение для таких приложений следующего поколения, как дисплеи дополненной реальности (AR), где высокая яркость и энергоэффективность являются первостепенными. Ключевой проблемой является достижение высокой эффективности вывода света (LEE), поскольку большая часть генерируемого света оказывается запертой внутри устройства из-за полного внутреннего отражения. Хотя обратное проектирование — вычислительная техника, автоматически оптимизирующая геометрию устройства — обещает многое, оно было вычислительно неразрешимым для µLED из-за необходимости моделирования тысяч пространственно некогерентных источников (например, от спонтанного излучения). Стандартные методы, такие как метод конечных разностей во временной области (FDTD), для этой задачи непозволительно медленны. Данная работа представляет возможности моделирования на основе метода Фурье-мод (FMM), которые преодолевают этот барьер, делая возможным эффективное обратное проектирование µLED, усиленных метаповерхностями.

2. Методология

2.1 Основы метода Фурье-мод (FMM)

Метод FMM, также известный как строгий анализ связанных волн (RCWA), моделирует электромагнитные поля в периодических, стратифицированных средах, разлагая их в усечённый базис Фурье. Уравнения Максвелла решаются в частотной области. Ключевое преимущество заключается в редукции трёхмерной задачи: размерности в плоскости (x,y) обрабатываются через разложение Фурье, в то время как z-размерность (стратификация) рассматривается аналитически. Это приводит к линейной системе, размер которой зависит только от фурье-гармоник в плоскости, а не от объёмной сетки, что даёт относительно компактную систему, решаемую прямыми методами.

2.2 Расширения для моделирования некогерентных источников

Стандартный FMM предполагает периодические источники, что для изолированного µLED в массиве создаёт нефизичную интерференцию. Чтобы смоделировать локализованный некогерентный источник (например, диполь в отдельном µLED), авторы используют векторную формулировку FMM. Это включает представление источника в виде суперпозиции мод Блоха. Затем полный отклик вычисляется суммированием вкладов от всех соответствующих векторов Блоха, эффективно моделируя одиночный излучатель в периодической среде без искусственной связи с его периодическими образами.

2.3 Интегрирование по зоне Бриллюэна

Для точного вычисления отклика локализованного источника выполняется интегрирование по зоне Бриллюэна (BZ) обратной решётки. Эта техника, заимствованная из смежных работ [17–19], производит выборку различных волновых векторов Блоха ($\mathbf{k}$) для построения полного отклика изолированного источника, обеспечивая физически корректные результаты для конфигурации массива µLED.

3. Техническая реализация и FMMAX

Метод реализован в инструменте под названием FMMAX. Ключевые инновации включают улучшенный алгоритм для автоматического вычисления векторных полей внутри слоёв и обработки структур, содержащих металлы, которые традиционно страдают от плохой сходимости в FMM [16]. Реализация позволяет эффективно повторно использовать вычислительно затратные собственные разложения при оптимизации параметров, что является критически важной особенностью для циклов обратного проектирования.

Прирост производительности

> 107 x

Быстрее, чем FDTD на CPU

Улучшение эффективности

~ 2 x

Увеличение LEE в спроектированном устройстве

4. Результаты и производительность

4.1 Сравнение скорости с FDTD

Моделирование на основе FMM даёт результаты, отлично согласующиеся с эталонными расчётами FDTD. Ключевой результат — вычислительная скорость: сообщается, что метод более чем в 107 раз быстрее, чем FDTD на CPU для задачи моделирования µLED. Это колоссальное ускорение превращает обратное проектирование из неразрешимой задачи в высокопрактичную.

4.2 Повышение эффективности вывода света

Используя свою структуру обратного проектирования, авторы оптимизировали метаповерхность, интегрированную поверх µLED. Оптимизированная конструкция удвоила эффективность вывода света (LEE) по сравнению с неоптимизированным базовым устройством. Это демонстрирует мощь метода для обнаружения неочевидных высокопроизводительных наноструктур.

5. Анализ сходимости

В статье рассматриваются исторические проблемы FMM, такие как медленная сходимость в металлических структурах и для локализованных источников. Показано, что их векторная формулировка и техники интегрирования по BZ значительно улучшают скорость сходимости, делая FMM надёжным и точным для геометрии µLED, которая включает полупроводниковые слои и потенциально металлические контакты или зеркала.

6. Демонстрация обратного проектирования

Продемонстрировано основное применение: автоматизированное обратное проектирование метаповерхности для повышения LEE. Пространство проектирования, вероятно, включало такие параметры, как форма, размер и расположение мета-атомов. Цикл оптимизации, ставший возможным благодаря быстрому моделированию, успешно исследовал это высокоразмерное пространство, чтобы найти структуру, максимизирующую долю света, выходящего из устройства.

7. Ключевая идея и аналитическая перспектива

Ключевая идея:

Прорыв статьи заключается не в новом алгоритме как таковом, а в стратегическом возрождении и улучшении существующего (FMM) для проблемы (обратное проектирование с некогерентными источниками), где сообщество уперлось в вычислительную стену. В то время как другие исследовали новые методы факторизации [13,14] или следовые формулировки [15] для снижения затрат, эта работа доказывает, что с правильными численными корректировками — векторные поля, интегрирование по BZ — «стандартный» метод может быть не просто достаточным, но и феноменально эффективным. Это классический случай, когда инженерная изобретательность превосходит погоню за чисто теоретической новизной.

Логическая последовательность:

Аргументация убедительна: 1) µLED нуждаются в обратном проектировании для эффективности, 2) некогерентные источники делают его слишком медленным, 3) FMM имеет врождённые скоростные преимущества для стратифицированных задач, 4) но у него есть известные недостатки для металлов и локализованных источников, 5) вот наши исправления, 6) теперь он в 10^7 раз быстрее и работает, 7) смотрите, мы спроектировали лучшее устройство. Последовательность от идентификации проблемы через техническое решение к осязаемому результату безупречна.

Сильные стороны и недостатки:

Сильные стороны: Ускорение в 10^7 раз — это нокаутирующий удар. Демонстрация реального устройства с удвоенной производительностью переводит работу из теории в практическую плоскость. Фокус на исправлении исторических слабостей FMM показывает глубокое техническое понимание.
Недостатки и вопросы: В статье мало деталей о самом алгоритме обратного проектирования (например, какой метод сопряжённых переменных, оптимизатор?). Утверждение о «сопоставимой точности» с FDTD требует проверки — по каким метрикам? Диаграммы направленности в дальней зоне? Интенсивности в ближней зоне? Производительность FMMAX на чрезвычайно сложных, нестратифицированных 3D-геометриях остаётся неподтверждённой. Как и во многих работах по обратному проектированию в фотонике, технологичность и устойчивость (например, к погрешностям изготовления) спроектированной метаповерхности не обсуждаются, что является критическим пробелом, отмеченным в обзорах области, таких как работы Molesky et al. (Nature Photonics, 2018).

Практические выводы:

Для компаний AR/VR: Этот инструмент может радикально ускорить цикл НИОКР для дисплеев на µLED. Инвестиции в такую технологию моделирования или лицензирование её — это высокоэффективный шаг.
Для исследователей: Урок ясен — пересмотрите «решённые» численные методы с современной точки зрения и конкретными ограничениями задачи; огромные выгоды могут скрываться на виду. Следующий шаг — интегрировать этот решатель с надёжными структурами обратного проектирования, учитывающими технологические ограничения, такие как минимальный размер элементов, как это исследуется в работах типа «Inverse design in nanophotonics» Jiang and Fan (Nature Reviews Materials, 2020).
Для разработчиков инструментов: FMMAX представляет собой эталон. Задача состоит в том, чтобы распространить его принципы на ещё более широкие классы устройств, возможно, интегрируя суррогатные модели машинного обучения для самых затратных шагов, чтобы добиться дальнейшего ускорения.

8. Технические детали и математическая формулировка

Суть FMM включает разложение периодической диэлектрической проницаемости $\epsilon(x,y)$ и электромагнитных полей в ряды Фурье:

$$ \epsilon(x,y) = \sum_{m,n} \tilde{\epsilon}_{mn} e^{j(mG_x x + nG_y y)} $$ $$ \mathbf{E}(x,y,z) = \sum_{m,n} \tilde{\mathbf{E}}_{mn}(z) e^{j[(k_x+mG_x)x + (k_y+nG_y)y]} $$ где $G_x, G_y$ — векторы обратной решётки, а $\mathbf{k}=(k_x, k_y)$ — волновой вектор Блоха. Подстановка в уравнения Максвелла приводит к системе связанных обыкновенных дифференциальных уравнений по $z$ для фурье-амплитуд $\tilde{\mathbf{E}}_{mn}(z)$, которая решается нахождением собственных мод в каждом слое и сопоставлением граничных условий.

Мощность для некогерентного источника вычисляется интегрированием по положениям источников и векторам Блоха: $$ P_{\text{ext}} \propto \int_{\text{BZ}} d\mathbf{k} \sum_{\text{sources}} |\mathbf{E}_{\text{far}}(\mathbf{k}, \mathbf{r}_s)|^2 $$ где некогерентность учитывается суммой интенсивностей (а не полей).

9. Экспериментальные результаты и описание графиков

Рисунок (концептуальное описание): В статье, вероятно, содержится ключевой рисунок, сравнивающий LEE базового и спроектированного обратным методом µLED. По оси абсцисс может быть представлена длина волны (например, 450–650 нм для синего/зелёного/красного светодиода), а по оси ординат — LEE (0–100%). Можно ожидать увидеть две кривые: 1) более низкую, плоскую кривую для неоптимизированного планарного или простого µLED и 2) значительно более высокую кривую для устройства, усиленного метаповерхностью, возможно, с резонансными пиками, где метаповерхность особенно эффективно выводит свет. Второй график может показывать сходимость метода FMM в зависимости от числа фурье-гармоник, демонстрируя быструю сходимость к стабильному значению LEE с улучшенной формулировкой, в отличие от более медленной или нестабильной сходимости для классического подхода FMM.

10. Аналитическая структура: рабочий процесс обратного проектирования

Пример: Проектирование метаповерхности для синего µLED

  1. Определение задачи: Цель: Максимизировать LEE на длине волны 450 нм для µLED с заданной структурой эпитаксиальных слоёв (например, на основе GaN). Ограничения: Период метаповерхности фиксирован шагом пикселя (например, 1 мкм), высота мета-атома ограничена возможностями изготовления.
  2. Параметризация: Определите элементарную ячейку метаповерхности. Простая параметризация может представлять собой прямоугольный наностолбик с переменными: ширина $w_x$, ширина $w_y$, угол поворота $\theta$ и материал (например, TiO$_2$).
  3. Моделирование: Для заданного набора параметров $(w_x, w_y, \theta)$ используйте FMMAX для вычисления LEE. Это включает решение для полей от ансамбля некогерентных диполей, размещённых в активной области квантовой ямы, и интегрирование вектора Пойнтинга, направленного вверх.
  4. Цикл оптимизации: Используйте градиентный оптимизатор (например, метод сопряжённых переменных) или алгоритм глобального поиска (например, байесовскую оптимизацию) для изменения $(w_x, w_y, \theta)$ и максимизации LEE. Ускорение FMMAX в 10^7 раз позволяет этому циклу выполняться за часы вместо лет.
  5. Верификация и вывод: Оптимизатор сходится к оптимальной форме столбика. Финальный шаг — полное верификационное моделирование и генерация файлов для изготовления (GDSII).

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

  • Полноцветные дисплеи на µLED: Одновременное обратное проектирование метаповерхностей для красных, зелёных и синих субпикселей для баланса эффективности и чистоты цвета.
  • Формирование луча: Расширение целевой функции за пределы общей LEE для включения управления профилем луча в дальней зоне (например, коллимации для проекционных приложений), аналогично целям в проектировании макроскопических светодиодов.
  • Интеграция с активной настройкой: Проектирование метаповерхностей, совместимых с жидкими кристаллами или материалами с фазовым переходом для динамически настраиваемых µLED после изготовления.
  • Совместное проектирование с тепловым менеджментом: Обратное проектирование, учитывающее как фотонные характеристики, так и тепловое рассеяние, поскольку падение эффективности при высоких токах является серьёзной проблемой для µLED.
  • Совместное проектирование алгоритмов и аппаратного обеспечения: Реализация основного решателя FMMAX на GPU или специализированных AI-ускорителях для достижения дальнейшего ускорения, приближаясь к исследованию проектов в реальном времени.
  • Более широкий круг фотонных задач: Применение улучшенной структуры FMM к другим задачам с некогерентными источниками, таким как оптимизация светоизлучающих электрохимических ячеек (LEC), улавливание света в солнечных элементах или инфракрасные излучатели для сенсоров.

12. Список литературы

  1. Z. Zhou et al., "Augmented reality and virtual reality displays: emerging technologies and future perspectives," Light: Science & Applications, 2021.
  2. H. S. Chen et al., "Micro-LED technology for next-generation displays," Journal of the Society for Information Display, 2020.
  3. J. A. Fan et al., "Inverse design of nanophotonic structures," Nature Photonics, 2010.
  4. S. Molesky et al., "Inverse design in nanophotonics," Nature Photonics, 2018.
  5. J. Jiang and J. A. Fan, "Global optimization of dielectric metasurfaces using a physics-driven neural network," Nano Letters, 2019.
  6. K. J. Vahala, "Optical microcavities," Nature, 2003.
  7. M. L. Brongersma et al., "Plasmonics for improved photovoltaic devices," Nature Materials, 2010.
  8. P. Bermel et al., "Design and global optimization of high-efficiency thermophotovoltaic systems," Optics Express, 2010.
  9. J. D. Joannopoulos et al., "Photonic Crystals: Molding the Flow of Light," Princeton University Press, 2008.
  10. L. Li, "New formulation of the Fourier modal method for crossed surface-relief gratings," Journal of the Optical Society of America A, 1997.
  11. A. David et al., "Simulation of light extraction from light-emitting diodes with patterned surfaces by the finite-difference time-domain method," IEEE Transactions on Electron Devices, 2010.
  12. T.-Y. Huang et al., "Design and simulation of GaN-based micro-LEDs with vertical sidewalls," IEEE Photonics Technology Letters, 2016.
  13. R. Pestourie et al., "A computational framework for infinite-dimensional inverse design using factorization," arXiv preprint, 2022.
  14. O. D. Miller et al., "Photonic design: From fundamental solar cell physics to computational inverse design," IEEE Journal of Photovoltaics, 2012.
  15. H. Chung and O. D. Miller, "Tunable metasurfaces via subwavelength phase shifters with uniform amplitude," Scientific Reports, 2020.
  16. G. Granet and B. Guizal, "Efficient implementation of the coupled-wave method for metallic lamellar gratings in TM polarization," Journal of the Optical Society of America A, 1996.
  17. P. Lalanne and G. M. Morris, "Highly improved convergence of the coupled-wave method for TM polarization," Journal of the Optical Society of America A, 1996.
  18. E. Silberstein et al., "Use of grating theories in integrated optics," Journal of the Optical Society of America A, 2001.
  19. M. G. Moharam et al., "Stable implementation of the rigorous coupled-wave analysis for surface-relief gratings: enhanced transmittance matrix approach," Journal of the Optical Society of America A, 1995.