Метод FDTD проти методу передатної матриці (TMM): вибір оптимального методу оптичного моделювання.

1. Вступ

Обидва методи – метод матриці передачі (TMM) і метод скінченних різниць у часовій області (FDTD) – широко використовуються для розв’язання оптичних задач. Однак, попри те, що в основі обох методів лежать рівняння Максвелла, вони розв’язують різні математичні задачі та призначені для різних класів фізичних систем.

Однією з найпоширеніших помилок, яку роблять нові користувачі, є припущення, що «найсучасніший» метод моделювання завжди є найкращим. На практиці це рідко відповідає дійсності. У обчислювальній фізиці метою зазвичай є не розв’язання найскладнішої можливої задачі, а розв’язання найпростішої задачі, яка здатна точно описувати потрібні фізичні явища.

Це важлива ідея, оскільки обчислювальні витрати надзвичайно швидко зростають зі збільшенням розмірності та складності моделі. Ретельно підібрана спрощена модель часто може давати результати, які є більш зрозумілими з фізичної точки зору, легшими для інтерпретації та значно швидшими в обчисленні.

Метод матриці передачі є чудовим прикладом цього принципу. Замість того, щоб розв’язувати повну електромагнітну задачу в просторі та часі, метод матриці передачі передбачає, що структура є, по суті, одновимірною та складається з шарів, як описано в Виведення методу матриці передачі. За таких наближень електромагнітне поле всередині кожного шару можна аналітично представити як суперпозицію хвиль, що поширюються в прямому та зворотному напрямках.

У результаті надзвичайно складні проблеми, пов’язані з інтерференцією, часто можна вирішити за лічені мілісекунди з дуже високою чисельною точністю, при цьому зберігаючи важливі фізичні характеристики резонатора, поведінку стоячих хвиль, спектри відбиття та пропускання, а також ефекти багатошарової інтерференції.

На відміну від цього, метод FDTD безпосередньо розв’язує рівняння Максвелла на просторовій сітці, відстежуючи зміни електромагнітних полів у часі. Детальний математичний виклад наведено в математичні обґрунтування та теоретичні основи методу FDTD. Це робить метод FDTD значно більш загальним і набагато повнішим з фізичної точки зору, але водночас і значно більш обчислювально затратним.

2. Основна відмінність між методами TMM і FDTD.

Основна концептуальна відмінність між цими двома методами полягає в тому, що метод TMM спрощує фізичні процеси перед розв’язанням рівнянь, тоді як метод FDTD намагається безпосередньо розв’язати повну електромагнітну задачу.

У методі матриці передачі передбачається, що структура складається з плоских шарів, розташованих один над одним в одному напрямку. Оскільки властивості матеріалу змінюються лише вздовж однієї осі, рівняння Максвелла спрощуються до значно простішої одновимірної хвильової задачі.

Це означає, що багато складних електромагнітних ефектів просто не відображаються в математичних моделях. Наприклад, стандартна одномірна матрична теорія передачі (TMM) за своєю суттю не враховує:

• Дифракція
• Бічне поширення
• Викривлені хвильові фронти
• Розсіювання на довільних геометричних поверхнях
• Вигин хвилеводу
• Локалізовані дефекти у двох або трьох вимірах

Це не є недоліком даного методу. Навпаки, саме завдяки цьому метод працює так швидко.

Усуваючи зайві ступені свободи, метод TMM зосереджує обчислювальні зусилля лише на фізичних процесах, які мають значення для шаруватих оптичних систем. Для тонкоплівкової оптики, діелектричних дзеркал, резонаторів Фабрі-Перо, багатошарових OLED-структур, інтерференційних ефектів у сонячних батареях і багатошарових покриттів ця наближена модель часто є надзвичайно точною.

Метод FDTD використовує протилежний підхід. Замість того, щоб аналітично спрощувати геометрію, простір поділяється на чисельну сітку, і рівняння Максвелла розв’язуються безпосередньо в межах усієї області моделювання. Електромагнітні поля поширюються через структуру природним чином як функція часу.

Оскільки зберігається повна інформація про зміни поля в часі, метод FDTD може моделювати дифракцію, поширення під кутом, розсіювання, перехідні процеси, інтерференцію у багатовимірному просторі, хвилеводні режими та довільні нанофотонні структури.

3. Порівняння фізичного розуміння, яке забезпечує кожен метод.

Метод матриці передачі часто дає надзвичайно чітке фізичне уявлення про резонансні оптичні системи. Наприклад, карта щільності фотонів у резонаторі Фабрі-Перо, наведена в навчальний матеріал про резонатор Фабрі-Перо та метод матриці передачі та показана в ??, безпосередньо демонструє резонансні довжини хвиль і структуру стоячих хвиль всередині резонатора.

Оскільки проблема зводиться до однієї розмірності, фізичну інтерпретацію резонансів здійснити легко. Чітко проявляються режими стоячих хвиль, резонансні довжини хвиль чітко визначені, і зв’язок між довжиною резонатора та умовою резонансу одразу стає помітним.

У багатьох випадках така простота є дуже бажаною. Мета моделювання часто полягає не лише в тому, щоб створити зображення, а й у тому, щоб сформувати фізичне розуміння. Надмірно складні моделювання іноді можуть затьмарювати основні фізичні принципи, замість того щоб їх прояснювати.

На відміну від цього, багатовимірна FDTD-симуляція, наведена в FDTD Fabry-Pérot cavity tutorial та показана на ??, містить набагато складнішу просторову поведінку. Тут хвиля поширюється в багатьох вимірах, розтікається вбік, дифрактує та утворює складні інтерференційні картини по всій області резонатора.

Ця додаткова реалістичність має надзвичайно важливе значення, коли сама структура є справді багатовимірною. Однак це також означає, що результати можуть стати значно складнішими для інтерпретації. Користувачу тепер необхідно розрізняти резонанси в порожнинах, ефекти дифракції, кутову поширеність, ефекти, пов’язані з кінцевими розмірами, перехідні явища та числові артефакти.

Іншими словами, метод FDTD часто дозволяє більш детально моделювати фізичні процеси, але не завжди ці моделі є більш корисними.

4. Обчислювальні витрати та масштабованість

Однією з найважливіших практичних відмінностей між цими двома методами є обчислювальні витрати.

Метод матриці передачі є надзвичайно ефективним з обчислювальної точки зору, оскільки він розв’язує лише одновимірну задачу для шаруватих систем. Типові обчислення для багатошарових оптичних систем завершуються за частки секунди, навіть для структур, що містять багато шарів і широкий діапазон довжин хвиль.

На відміну від цього, обчислювальна складність методів FDTD дуже швидко зростає зі збільшенням кількості вимірів. Якщо просторову роздільну здатність збільшити вдвічі, то кількість точок сітки зросте приблизно пропорційно наступному:

• \(N\) за 1 день
• \(N^2\) у двовимірному просторі
• \(N^3\) у 3D.

Крім того, крок часу також має зменшуватися, щоб забезпечити чисельну стабільність. Як наслідок, загальні обчислювальні витрати для методу FDTD можуть значно зрости у випадку великих або високороздільних обчислень.

Саме тому для багатьох реалістичних тривимірних чисельних симуляцій, що базуються на методі FDTD, потрібні:

• Апаратне прискорення за допомогою графічного процесора (GPU).
• Великий обсяг пам’яті
• Паралельні обчислення
• Тривалий час роботи

Важливо зазначити, що це не означає, що метод FDTD є «поганим». Це просто означає, що цей метод використовується для розв’язання значно складнішої математичної задачі.

Якщо структура за своєю суттю є одновимірною, використання методу FDTD може лише незначно підвищити фізичну точність, водночас значно збільшуючи обчислювальні витрати.

5. У яких випадках варто використовувати метод матриці передачі?

Загалом, метод матриці передачі слід використовувати насамперед, коли структура є приблизно плоскою та шаруватою.

Метод передачі матриц (TMM) особливо добре підходить для багатошарових оптичних моделювань, таких як оптичні фільтри, антиблікові покриття та резонатори Фабрі-Перо. Загалом, TMM особливо добре підходить для:

• Оптика тонких плівок
• Резонатори Фабрі-Перо
• Оптичні покриття
• Діелектричні дзеркала
• Дифузійні дзеркала Брегга (ДДБ)
• Оптичне моделювання OLED-дисплеїв
• Обчислення взаємного впливу сонячних батарей
• Обчислення поглинання в багатошарових структурах.
• Швидке послідовне змінювання параметрів

TMM також є ідеальним рішенням у таких випадках:

• Вам потрібні надзвичайно швидкі обчислювальні моделі.
• Вам потрібно оптимізувати багато параметрів.
• Вам потрібні чіткі резонансні спектри.
• Вам потрібне просте фізичне пояснення.
• За своєю суттю, ця структура є одновимірною.

У багатьох практичних задачах оптичної інженерії метод передачі матрицею (TMM) може давати результати, які майже не відрізняються від результатів значно складніших обчислювальних методів.

6. В яких випадках варто використовувати метод FDTD?

Метод FDTD набуває важливого значення, коли припущення, на яких базується метод TMM, перестають бути обґрунтованими.

Зокрема, метод FDTD слід використовувати в таких випадках:

• Геометрія є справді двовимірною або тривимірною.
• Дифракція має важливе значення.
• Хвиля поширюється під великим кутом.
• Конструкція містить вигини або криві лінії.
• Цей пристрій містить ізольовані наноструктури.
• Розсіювання має значення.
• Важливо враховувати динаміку змін у перехідних полях.
• Вам потрібні повні розподіли електромагнітного поля.

Наприклад:

• Photonic crystal waveguides
• Bragg gratings
• Mach-Zehnder interferometers
• Ring resonators
• Tilted facet waveguides
• Tapered waveguides and mode expansion
• Розсіювання наночастинок
• Локалізовані дефекти
• Довільні нанофотонні структури

У цих системах спрощення, які використовує метод передатних матриць, більше не відповідають фізичній реальності, тому виникає необхідність у додаткових обчисленнях за допомогою методу скінченних різниць у часовій області.

7. Підсумок: вирішення найпростішої, але реальної проблеми.

Одне з найважливіших вмінь у обчислювальній фізиці – це вміння визначати оптимальний рівень складності моделі.

Більш складні моделювання не обов’язково є кращими. У багатьох випадках простіша модель може забезпечити:

• Більш чітке фізичне розуміння.
• Прискорення процесу ітерацій
• Простіше налагодження.
• Покращена оптимізація параметрів
• Більш надійна інтерпретація

Метод матриці передачі є ефективним саме тому, що він усуває зайву складність із багатошарових оптичних задач. Він використовує математичне спрощення таким чином, що дозволяє надзвичайно ефективно розв’язувати основні фізичні задачі, пов’язані з інтерференцією.

Метод FDTD стає цінним, коли фізичні явища, які раніше не враховувалися, набувають важливого значення. У цей момент додаткові обчислювальні витрати виправдані, оскільки простіші наближення вже не є достатніми.

На практиці досвідчені дослідники часто використовують обидва методи одночасно: TMM для швидкого дослідження та оптимізації, а потім FDTD для детальної багатовимірної перевірки та аналізу поля.

Отже, мета полягає не в тому, щоб обрати «найсучаснішу» методику моделювання, а скоріше в тому, щоб обрати метод, який найкраще відповідає фізичним принципам, що лежать в основі розв’язуваної задачі.