Що таке метод матриці передачі?

1. Вступ

Метод передавальної матриці (МПМ) — це швидкий і обчислювально ефективний метод, який використовується для моделювання поширення електромагнітних хвиль у багатошарових тонкоплівкових оптичних структурах. Він широко застосовується в оптиці, фотоніці та моделюванні напівпровідникових пристроїв для обчислення залежності від довжини хвилі відбиття, пропускання, поглинання, інтерференції та розподілу електричного поля в оптичних структурах, що складаються з шаруватих матеріалів з різними показниками заломлення. До сфер застосування належать тонкоплівкові сонячні батареї, мікрорезонатори OLED, діелектричні дзеркала, антиблікові покриття, розподілені дзеркала Бреджа (DBR), оптичні фільтри, покриття для хвилеводів та інші плоскі фотонні пристрої.

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

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

2. Фізична основа методу матриці передачі

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

\[ E(z)=E^{+}e^{-jkz}+E^{-}e^{jkz} \]

де \(E^{+}\) і \(E^{-}\) – відповідно амплітуди прямо- і зворотньопоширюваних хвиль, \(z\) – напрямок поширення хвилі в багатошаровій структурі, і

\[ k=\frac{2\pi \tilde{n}}{\lambda}, \qquad \tilde{n}=n+j\kappa \]

це комплексний оптичний хвильовий вектор, де \(n\) – показник заломлення, \(\kappa\) описує оптичне поглинання, а \(\lambda\) – довжина хвилі у вакуумі. Як показано на ??, хвилі, що поширюються вперед і назад, взаємодіють на кожній межі між матеріалами за допомогою часткового відбиття та проходження.

Це зображення хвилі отримано шляхом розв’язання одномірного хвильового рівняння Гельмгольца, яке виведено з рівнянь Максвелла для однорідного середовища:

\[ \frac{d^2E(z)}{dz^2}+k^2E(z)=0 \]

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

\[ r=\frac{n_1-n_2}{n_1+n_2}, \qquad t=\frac{2n_1}{n_1+n_2} \]

де \(n_1\) і \(n_2\) — комплексні показники заломлення сусідніх шарів. Ці коефіцієнти визначають, як хвилі, що поширюються в прямому та зворотному напрямках, частково відбиваються та проходять крізь кожен матеріал на межі між ними.

3. Побудова матриці перетворення

Як показано на ??, оптичне поле всередині кожного шару представлене хвилею, що поширюється вперед \(E^+\), і хвилею, що поширюється назад \(E^-\):

\[ \mathbf{E}= \begin{bmatrix} E^{+} \\ E^{-} \end{bmatrix} \]

У міру поширення електромагнітної хвилі крізь шар товщиною \(d\), оптичне поле накопичує фазу та зазнає загасання. Це поширення крізь шар можна описати так:

\[ \begin{bmatrix} E^{+}(d) \\ E^{-}(d) \end{bmatrix} = \begin{bmatrix} e^{-jkd} & 0 \\ 0 & e^{jkd} \end{bmatrix} \begin{bmatrix} E^{+}(0) \\ E^{-}(0) \end{bmatrix} \]

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

\[ \begin{bmatrix} E^{+}_2 \\ E^{-}_2 \end{bmatrix} = \frac{1}{t} \begin{bmatrix} 1 & r \\ r & 1 \end{bmatrix} \begin{bmatrix} E^{+}_1 \\ E^{-}_1 \end{bmatrix} \]

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

\[ \mathbf{E}_{out} = M_N I_N \cdots M_2 I_2 M_1 I_1 \mathbf{E}_{in} \]

де \(M\) позначає поширення через шар, а \(I\) – передачу та відбиття на межі розділу. Таким чином, оптичну характеристику всього пристрою можна обчислити, використовуючи послідовність простих матричних множень.

Завдяки цьому підходу можна надзвичайно ефективно розраховувати процеси відбиття, передавання, поглинання, інтерференції та формування стоячих хвиль у складних багатошарових оптичних структурах, при цьому зберігаючи прямий зв’язок із рівняннями Максвелла. Оскільки метод матриці передавання безпосередньо розв’язує оптичне хвильове рівняння, він природним чином враховує ефекти інтерференції тонких плівок. Коли світло відбивається від кількох меж поділу, відбиті хвилі інтерферують конструктивно або деструктивно, залежно від оптичної фази, накопиченої в кожному шарі, \( \phi = \frac{2\pi n d}{\lambda} \), де \(n\) – показник заломлення, \(d\) – товщина шару, а \(\lambda\) – довжина хвилі. Як наслідок, навіть незначні зміни товщини нанометрового масштабу можуть суттєво змінювати оптичні характеристики пристрою. Ця фізика інтерференції визначає роботу антиблікових покриттів, резонаторів Фабрі–Перо, розподілених дзеркал Брэгга, оптичних фільтрів, мікрорезонаторів OLED і тонкоплівкових сонячних елементів.

4. Що саме обчислює метод матриці передачі?

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

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

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

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

\[ G(z)=\alpha(z)\Phi(z) \]

де \(G(z)\) – локальна швидкість генерації заряду, \(\alpha(z)\) – коефіцієнт оптичного поглинання, а \(\Phi(z)\) – локальний потік фотонів. Приклад профілю генерації, розрахованого за допомогою методу матриці передачі, показано на ??.

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

5. Переваги та обмеження методу матриці передачі

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

До основних переваг TMM належать:

• Надзвичайно швидке моделювання: ідеально підходить для проведення серії розрахунків із різними параметрами та оптимізації пристроїв.
• Точний розрахунок інтерференції тонких плівок: забезпечує ретельне моделювання когерентної інтерференції та ефектів стоячих хвиль.
• Безпосередньо отримувані оптичні вихідні дані: електричне поле, щільність фотонів і профілі поглинання легко визначаються.
• Чудово підходить для створення багатошарових плоских структур: особливо корисний для виробництва OLED-дисплеїв, сонячних батарей, оптичних покриттів і дзеркальних покриттів.
• Просте поєднання з електричними моделями: профілі оптичного генерування можна безпосередньо використовувати в симуляціях, що базуються на моделях дрейфу та дифузії.

Однак, метод TMM також має певні обмеження. За своєю суттю, цей метод передбачає, що оптичну структуру можна розглядати як приблизно одновимірну, причому властивості матеріалу змінюються переважно вздовж напрямку росту. Як наслідок, складно точно змоделювати довільні тривимірні схеми розсіювання.

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

Для розв’язання складних тривимірних оптичних задач застосовуються чисельні методи, такі як: Метод скінченних різниць у часовій областіу таких випадках зазвичай більш доцільно використовувати методи скінченних елементів (МСЕ) або строгий метод аналізу взаємодії хвиль (МСВХ). Отже, метод передачі матрицею (МПМ) найкраще підходить для плоских багатошарових пристроїв, де важливі ефекти інтерференції, але структура залишається приблизно одновимірною.

8. Підсумок.

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

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

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