Na czym polega metoda macierzy przenoszenia?

1. Wprowadzenie

Metoda macierzy przenoszenia (ang. Transfer Matrix Method, TMM) to szybka i wydajna obliczeniowo technika, która służy do modelowania propagacji fal elektromagnetycznych w wielowarstwowych, cienkowarstwowych strukturach optycznych. Jest szeroko stosowana w optyce, fotonice i symulacjach urządzeń półprzewodnikowych do obliczania zależności odbicia, transmisji, absorpcji, interferencji i rozkładu pola elektrycznego od długości fali w układach optycznych składających się z warstwowych materiałów o różnych współczynnikach załamania. Przykłady zastosowań to cienkowarstwowe ogniwa słoneczne, mikrostruktury OLED, lustra dielektryczne, powłoki antyrefleksyjne, rozproszone zwierciadła Bragg'a (DBR), filtry optyczne, powłoki światłowodowe i inne planarne urządzenia fotoniczne.

Metoda TMM opiera się bezpośrednio na równaniach Maxwella i stanowi uproszczone, jednowymiarowe rozwiązanie równania fal elektromagnetycznych dla warstwowych ośrodków. Zamiast rozwiązywać pełne równania pola elektromagnetycznego w wielowymiarowej przestrzeni symulacji, metoda ta przedstawia światło w każdej warstwie jako fale rozchodzące się do przodu i do tyłu, przemieszczające się przez złożony profil współczynnika załamania.

Ponieważ metoda macierzy przenoszenia sprowadza problem optyczny do sekwencji prostych operacji macierzowych, jest ona znacznie szybsza niż pełne metody numeryczne, takie jak... Metoda FDTDPomimo swojej prostej konstrukcji obliczeniowej, metoda TMM precyzyjnie oddaje spójne efekty interferencji cienkich warstw, w tym powstawanie fal stojących, rezonanse Fabry’ego-Perota, właściwości optycznych rezonatorów oraz wzmocnienie pola wewnątrz układów dielektrycznych i wielowarstwowych powłok optycznych. Dzięki temu metoda ta jest jednym z najczęściej stosowanych sposobów modelowania spójnej optyki w warstwowych urządzeniach cienkowarstwowych.

2. Fizyczna zasada działania metody macierzy przenoszenia

Metoda macierzy przenoszenia przedstawia pole optyczne wewnątrz każdej warstwy jako superpozycję fal elektromagnetycznych rozchodzących się w kierunku do przodu i do tyłu:

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

gdzie \(E^{+}\) i \(E^{-}\) oznaczają odpowiednio amplitudy fal rozchodzących się w kierunku do przodu i do tyłu, \(z\) to kierunek propagacji w strukturze wielowarstwowej, a

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

jest to złożony wektor falowy, gdzie \(n\) oznacza współczynnik załamania, \(\kappa\) opisuje absorpcję światła, a \(\lambda\) to długość fali w próżni. Jak pokazano na rysunku ??, fale rozchodzące się w kierunku do przodu i do tyłu oddziałują na sobie na każdej granicy materiału poprzez częściowe odbicie i transmisję.

Przedstawienie falowe uzyskujemy, rozwiązując jednowymiarowe równanie falowe Helmholtza, wyprowadzone z równań Maxwella dla jednorodnego ośrodka:

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

Na każdym z elementów granicznych stosuje się warunki brzegowe elektromagnetyczne Maxwella, aby zapewnić ciągłość stycznych składowych pola elektrycznego i magnetycznego. Warunki te generują współczynniki odbicia i przepuszczalności Fresnela:

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

gdzie \(n_1\) i \(n_2\) oznaczają zespolone współczynniki załamania światła sąsiednich warstw. Współczynniki te określają, w jakim stopniu fale rozchodzące się w kierunku do przodu i do tyłu są częściowo odbijane i przepuszczane na każdej granicy między materiałami.

3. Tworzenie macierzy transformacji

Jak przedstawiono w ??, pole optyczne w każdej warstwie jest reprezentowane przez falę rozchodzącą się w kierunku do przodu \(E^+\) oraz falę rozchodzącą się w kierunku przeciwnym \(E^-\):

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

W miarę jak fala elektromagnetyczna rozchodzi się przez warstwę o grubości \(d\), pole optyczne kumuluje fazę i ulega osłabieniu. Proces rozchodzenia się fali przez warstwę można zapisać w następujący sposób:

\[ \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} \]

Na granicy między dwoma materiałami część pola optycznego jest transmitowana, a część odbijana. Tę zależność można również zapisać w postaci macierzowej:

\[ \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} \]

W związku z tym cała, wielowarstwowa struktura optyczna jest przedstawiana jako sekwencja macierzy propagacji i macierzy interfejsów:

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

gdzie \(M\) oznacza propagację przez warstwę, a \(I\) – transmisję i odbicie na granicy. W związku z tym, odpowiedź optyczną całego urządzenia można obliczyć, stosując sekwencję prostych mnożeń macierzy.

To podejście umożliwia niezwykle efektywne obliczanie procesów odbicia, transmisji, absorpcji, interferencji oraz powstawania fal stojących w złożonych, wielowarstwowych strukturach optycznych, przy jednoczesnym zachowaniu bezpośredniego związku z równaniami Maxwella. Ponieważ metoda macierzy przenoszenia rozwiązuje bezpośrednio równanie fali optycznej, naturalnie uwzględnia efekty interferencji cienkich warstw. W miarę odbijania się światła od wielu powierzchni, odbite fale interferują ze sobą konstruktywnie lub destruktywnie, w zależności od fazy optycznej, która gromadzi się w każdej warstwie \( \phi = \frac{2\pi n d}{\lambda} \), gdzie \(n\) to współczynnik załamania, \(d\) grubość warstwy, a \(\lambda\) długość fali. W rezultacie nawet zmiany grubości rzędu nanometrów mogą znacząco modyfikować właściwości optyczne urządzenia. Ta fizka interferencji reguluje działanie powłok antyrefleksyjnych, rezonatorów Fabry-Perota, rozproszonych zwiercideł Bragg'a (DBR), filtrów optycznych, mikrorezonatorów OLED oraz cienkowarstwowych ogniw słonecznych.

4. Co oblicza metoda macierzowa transferu?

Po skonstruowaniu globalnej macierzy transformacji dla wielowarstwowej struktury można obliczyć pole elektromagnetyczne w dowolnym punkcie wewnątrz urządzenia. Umożliwia to określenie przestrzennego rozkładu pola optycznego, współczynników odbicia, transmisji i absorpcji w całej strukturze.

Ponieważ metoda macierzy przenoszenia rozwiązuje równanie fal optycznych bezpośrednio, w symulacji naturalnie pojawiają się efekty koherentnej interferencji. W miarę jak światło odbija się wielokrotnie między powierzchniami, konstruktywna i destruktywna interferencja tworzy wewnątrz urządzenia wzory fal stojących.

Rozkład gęstości fotonów przedstawiony na rysunku ?? ilustruje, jak intensywność pola optycznego zmienia się w całej warstwowej strukturze. Obszary o wysokiej intensywności pola optycznego odpowiadają zwiększonej gęstości fotonów, co jest spowodowane interferencją optyczną i efektami rezonansowymi.

Lokalny rozkład pola optycznego bezpośrednio określa, w których miejscach wewnątrz urządzenia dochodzi do absorpcji fotonów. Łącząc pole elektromagnetyczne z współczynnikiem absorpcji materiału, można obliczyć profil generacji optycznej:

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

gdzie \(G(z)\) oznacza lokalną szybkość generowania ładunku, \(\alpha(z)\) to współczynnik absorpcji światła, a \(\Phi(z)\) to lokalny strumień fotonów. Przykładowy profil generowania, obliczony za pomocą metody macierzy przenoszenia, przedstawiono na rysunku ??.

Te obliczenia są powszechnie stosowane w cienkowarstwowych ogniwach słonecznych, diodach OLED, powłokach optycznych, zwierciadłach dielektrycznych, czujnikach optycznych i rezonansowych mikrorezonatorach. Ponieważ metoda macierzy przenoszenia jest stosunkowo prosta obliczeniowo, zazwyczaj można w ciągu ułamka sekundy obliczyć pełne widma optyczne i rozkłady pól wewnętrznych.

5. Zalety i ograniczenia metody macierzy przenoszenia.

Metoda macierzy przenoszenia (TMM) jest powszechnie stosowana, ponieważ łączy w sobie wysoką dokładność fizyczną z bardzo niskimi wymaganiami obliczeniowymi. W przypadku warstwowych struktur optycznych, TMM pozwala na odtworzenie tych samych zjawisk interferencji cienkich warstw, co znacznie bardziej złożone symulacje falowe, przy jednoczesnym zużyciu o wiele mniejszej ilości pamięci i czasu obliczeń.

Do najważniejszych zalet metody TMM należą:

• Bardzo szybkie symulacje: idealne do przeprowadzania analiz parametrów i optymalizacji urządzeń.
• Precyzyjne modelowanie interferencji cienkich warstw: dokładnie oddaje spójne zjawiska interferencyjne i efekty fal stojących.
• Bezpośrednie wyjścia optyczne: łatwo można uzyskać informacje o natężeniu pola elektrycznego, gęstości fotonów oraz profilach absorpcji.
• Doskonałe do tworzenia warstwowych struktur planarnych: szczególnie przydatne w produkcji wyświetlaczy OLED, ogniw słonecznych, powłok optycznych i zwierciadeł Braggowskich.
• Łatwe połączenie z modelami elektrycznymi: profile generacji optycznej można bezpośrednio wprowadzać do symulacji opartych na równaniach dryfu i dyfuzji.

Jednak metoda TMM ma również istotne ograniczenia. W założeniach metody przyjmuje się, że strukturę optyczną można traktować jako w przybliżeniu jednowymiarową, przy czym właściwości materiału zmieniają się głównie wzdłuż kierunku wzrostu. W konsekwencji, dokładne modelowanie dowolnych trójwymiarowych geometrii rozpraszania jest utrudnione.

• Przede wszystkim jednowymiarowy: ograniczony do struktur warstwowych.
• Ograniczenia w modelowaniu rozpraszania: rozpraszanie dyfuzjne i losowa nierówność powierzchni nie są uwzględniane w sposób ścisły.
• Nie stosujemy tu dowolnych kształtów geometrycznych: siatki dyfrakcyjne, nanocząstki, powierzchnie o fakturze oraz kryształy fotoniczne wymagają zazwyczaj rozwiązań opartych na pełnej analizie falowej.
• Założenia dotyczące koherencji: standardowy model TMM zakłada koherentną propagację wewnątrz każdej warstwy.

W przypadku złożonych problemów optycznych o wyraźnej trójwymiarowości stosuje się metody numeryczne, takie jak: Metoda FDTDW takich przypadkach zazwyczaj bardziej odpowiednie są metody elementów skończonych (MES) lub rygorystyczna analiza sprzężonych fal (RCWA). Dlatego metoda TMM najlepiej sprawdza się w przypadku płaskich, wielowarstwowych struktur, w których istotne są efekty interferencji, ale struktura pozostaje w przybliżeniu jednowymiarowa.

8. Podsumowanie

Metoda macierzy przenoszenia jest jedną z najważniejszych metod numerycznych w optyce cienkich warstw. Umożliwia ona szybkie i precyzyjne modelowanie odbicia, przepuszczalności, interferencji, absorpcji oraz generowania światła wewnątrz struktur wielowarstwowych.

Chociaż metoda ta jest stosowana głównie w przypadku struktur warstwowych, pozostaje jednym z najbardziej wydajnych obliczeniowo sposobów symulacji płaskich urządzeń optoelektronicznych.

W praktycznych procesach symulacji urządzeń, metoda TMM jest często stosowana w połączeniu z programami do rozwiązywania równań transportu elektrycznego, aby badać wzajemne oddziaływanie między optyką a fizyką półprzewodników.