Metoda FDTD a metoda macierzy przenoszenia (TMM): wybór odpowiedniej metody symulacji optycznej

1. Wprowadzenie

Zarówno metoda macierzowa metoda macierzy przenoszenia (TMM), jak i metoda różnic skończonych w dziedzinie czasu metoda FDTD są powszechnie stosowane do rozwiązywania problemów optycznych. Jednakże, mimo że obie metody opierają się ostatecznie na równaniach Maxwella, rozwiązują one bardzo różne problemy matematyczne i są przeznaczone do analizy różnych klas układów fizycznych.

Jednym z najczęstszych błędów popełnianych przez nowych użytkowników jest założenie, że „najbardziej zaawansowana” metoda symulacji jest zawsze najlepsza. W praktyce tak zwykle nie jest. W fizyce obliczeniowej celem nie jest zazwyczaj rozwiązywanie najbardziej skomplikowanych problemów, lecz raczej rozwiązywanie najprostszych problemów, które pozwalają na dokładne opisanie interesujących nas zjawisk fizycznych.

To ważna koncepcja, ponieważ koszt obliczeniowy rośnie bardzo szybko wraz ze wzrostem liczby wymiarów i złożoności modelu. Starannie dobrany, uproszczony model może często generować wyniki, które są bardziej zrozumiałe z punktu widzenia fizyki, łatwiejsze do interpretacji i wielokrotnie szybsze w obliczeniach.

Metoda macierzy przenoszenia stanowi doskonały przykład tej zasady. Zamiast rozwiązywać pełne równanie elektromagnetyczne w przestrzeni i czasie, metoda ta zakłada, że struktura jest w zasadzie jednowymiarowa i składa się z warstw, zgodnie z opisem w Wyprowadzenie metody macierzy przenoszenia. Przy takim przybliżeniu pole elektromagnetyczne wewnątrz każdej warstwy można zapisać analitycznie jako superpozycję fal rozchodzących się w kierunku do przodu i do tyłu.

Dzięki temu niezwykle skomplikowane problemy związane z interferencją można często rozwiązywać w ciągu milisekund, zachowując bardzo wysoką dokładność numeryczną, a jednocześnie uwzględniając kluczowe aspekty fizyki rezonatorów, charakterystykę fal stojących, widma odbicia i transmisji oraz efekty interferencji wielowarstwowych.

W przeciwieństwie do metody FDTD, która bezpośrednio rozwiązuje równania Maxwella na siatce przestrzennej, śledząc zmiany pól elektromagnetycznych w czasie, metoda FDTD jest znacznie bardziej ogólna i w większym stopniu oddaje fizyczną rzeczywistość, ale jednocześnie wymaga znacznie większych nakładów obliczeniowych. Szczegółowe wyprowadzenie matematyczne znajduje się w Wyprowadzenia i podstawy matematyczne metody FDTD.

2. Podstawowa różnica między metodą TMM a metodą FDTD.

Kluczowa różnica między tymi dwoma metodami polega na tym, że metoda TMM upraszcza model fizyczny przed rozwiązaniem równań, natomiast metoda FDTD próbuje rozwiązać pełne równanie elektromagnetyczne bezpośrednio.

W metodzie macierzy przenoszenia zakłada się, że struktura składa się z płaskich warstw ułożonych jedna na drugiej wzdłuż jednej osi. Ponieważ właściwości materiału zmieniają się tylko wzdłuż jednej osi, równania Maxwella upraszczają się do znacznie prostszego jednowymiarowego problemu falowego.

Oznacza to, że wiele skomplikowanych efektów elektromagnetycznych po prostu nie występuje w modelach matematycznych. Na przykład, standardowy jednowymiarowy model TMM nie uwzględnia naturalnie następujących elementów:

• Dyfrakcja
• Rozprzestrzenianie się w bok
• Zakrzywione fronty fal
• Rozpraszanie na losowych powierzchniach.
• Zakrzywienia światłowodów
• Lokalne defekty występujące w dwóch lub trzech wymiarach.

To nie jest słabość tej metody. Wręcz przeciwnie, to właśnie dlatego metoda ta jest tak skuteczna.

Eliminując zbędne stopnie swobody, metoda TMM koncentruje obliczenia wyłącznie na aspektach fizycznych istotnych dla warstwowych układów optycznych. W przypadku optyki cienkowarstwowej, dielektrycznych luster, rezonatorów Fabry-Pérota, układów OLED, efektów interferencyjnych w ogniwach słonecznych oraz powłok wielowarstwowych, to przybliżenie jest często niezwykle dokładne.

Metoda FDTD stosuje odmienne podejście. Zamiast upraszczać geometrię analitycznie, przestrzeń dzieli się na numeryczną siatkę, a równania Maxwella rozwiązuje się bezpośrednio w całym obszarze symulacji. Pola elektromagnetyczne rozprzestrzeniają się w strukturze w sposób naturalny, jako funkcja czasu.

Ponieważ zachowywana jest pełna ewolucja pola, metoda FDTD pozwala modelować dyfrakcję, propagację pod kątem, rozpraszanie, zachowanie w czasie przejściowym, interferencję w przestrzeni o wyższej liczbie wymiarów, mody falowodu oraz dowolne nanostruktury fotoniczne.

3. Porównanie wiedzy na temat właściwości fizycznych, którą dostarcza każda z metod.

Metoda macierzy przenoszenia często pozwala uzyskać bardzo precyzyjne zrozumienie fizycznych właściwości rezonansowych układów optycznych. Na przykład, mapa gęstości fotonów w rezonatorze Fabry-Pérota, przedstawiona w transfer matrix Fabry-Pérot cavity tutorial i pokazana w ??, bezpośrednio ujawnia rezonansowe długości fal i strukturę fal stojących wewnątrz rezonatora.

Ponieważ problem ogranicza się do jednej przestrzeni, rezonanse można łatwo interpretować w kategoriach fizycznych. Widocznie uwidaczniają się mody fal stojących, długości fal rezonansowych są wyraźnie określone, a związek między długością rezonatora a warunkiem rezonansu jest od razu widoczny.

W wielu sytuacjach taka prostota jest bardzo pożądana. Celem symulacji jest często nie tylko generowanie obrazów, ale przede wszystkim pogłębianie zrozumienia zjawisk fizycznych. Zbyt skomplikowane symulacje mogą czasami zaciemniać podstawowe zasady fizyki, zamiast je wyjaśniać.

W przeciwieństwie do tego, wielowymiarowa symulacja FDTD przedstawiona w FDTD Fabry-Pérot cavity tutorial (patrz ??) charakteryzuje się znacznie bogatszym zachowaniem przestrzennym. W tym przypadku fala rozchodzi się w wielu wymiarach, rozprzestrzenia się na boki, ulega dyfrakcji i tworzy złożone wzory interferencyjne w całej przestrzeni rezonatora.

To dodatkowe uwzględnienie realizmu jest niezwykle ważne, gdy sama struktura jest rzeczywiście wielowymiarowa. Jednak oznacza to również, że wyniki mogą stać się znacznie trudniejsze do interpretacji. Użytkownik musi teraz odróżniać rezonanse wewnątrz wnęk, efekty dyfrakcji, propagację pod określonym kątem, efekty związane z ograniczonymi rozmiarami, zachowanie w czasie przejściowym oraz artefakty numeryczne.

Innymi słowy, metoda FDTD często pozwala na bardziej szczegółowe modelowanie zjawisk fizycznych, ale nie zawsze prowadzi to do bardziej użytecznych wyników.

4. Koszty obliczeniowe i skalowalność

Jedną z najważniejszych praktycznych różnic między tymi dwiema metodami jest koszt obliczeniowy.

Metoda macierzy przenoszenia jest niezwykle wydajna obliczeniowo, ponieważ rozwiązuje jedynie jednowymiarowy problem warstwowy. Typowe obliczenia dotyczące wielowarstwowych układów optycznych zajmują ułamki sekundy, nawet w przypadku struktur zawierających wiele warstw i szeroki zakres długości fal.

W przeciwieństwie do tego, w symulacjach FDTD złożoność obliczeniowa bardzo szybko rośnie wraz ze wzrostem liczby wymiarów. Jeśli podwoimy rozdzielczość przestrzenną, liczba punktów siatki wzrasta w przybliżeniu zgodnie z następującym wzorem:

• \(N\) w perspektywie 1D
• \(N^2\) w 2D
• \(N^3\) w 3D

Ponadto, krok czasowy musi być odpowiednio zmniejszony, aby spełnić warunki stabilności numerycznej. W rezultacie, całkowity koszt obliczeniowy metody FDTD może być bardzo wysoki w przypadku dużych lub wymagających wysokiej rozdzielczości symulacji.

Dlatego wiele realistycznych symulacji 3D opartych na metodzie FDTD wymaga:

• Akceleracja sprzętowa za pomocą GPU.
• Duża pojemność pamięci.
• Obliczenia równoległe
• Długi czas pracy na jednym ładowaniu/zasilaniu.

Co ważne, nie oznacza to, że metoda FDTD jest „zła”. Oznacza to jedynie, że rozwiązuje ona znacznie bardziej złożony problem matematyczny.

Jeśli struktura jest zasadniczo jednowymiarowa, zastosowanie metody FDTD może przynieść bardzo niewielkie dodatkowe korzyści w zakresie dokładności fizycznej, jednocześnie znacznie zwiększając koszty obliczeniowe.

5. Kiedy należy stosować metodę macierzy przenoszenia?

Ogólnie rzecz biorąc, metoda macierzy przenoszenia powinna być pierwszą metodą, którą należy rozważyć, gdy struktura jest w przybliżeniu płaska i składa się z warstw.

Metoda TMM szczególnie dobrze nadaje się do symulacji warstwowych układów optycznych, takich jak filtry optyczne, powłoki antyrefleksyjne i rezonatory Fabry-Pérota. Ogólnie rzecz biorąc, metoda TMM jest szczególnie przydatna do:

• Optyka cienkich warstw
• Rezonatory Fabry-Pérota
• Powłoki optyczne
• Lustra dielektryczne
• Rozproszone lustra Bragga (DBR)
• Modelowanie optyczne wyświetlaczy OLED.
• Obliczenia dotyczące wzajemnych zakłóceń w działaniu ogniw słonecznych.
• Obliczenia absorpcji w warstwowych układach.
• Szybkie zmiany parametrów

TMM jest również idealny w następujących sytuacjach:

• Potrzebujesz bardzo szybkich symulacji.
• Należy zoptymalizować wiele parametrów.
• Chcesz uzyskać wyraźne, dobrze zdefiniowane widma.
• Chcesz prostego, fizycznego wyjaśnienia.
• Struktura ta jest zasadniczo jednowymiarowa.

W wielu praktycznych problemach z zakresu inżynierii optycznej metoda TMM może dawać wyniki, które są niemal nieodróżnialne od wyników uzyskiwanych za pomocą znacznie bardziej złożonych obliczeniowo metod.

6. Kiedy należy stosować metodę FDTD?

Metoda FDTD staje się istotna, gdy założenia, na których opiera się metoda TMM, przestają być prawdziwe.

Szczególnie warto stosować metodę FDTD, gdy:

• Geometria jest rzeczywiście dwuwymiarowa lub trójwymiarowa.
• Dyfrakcja jest ważna.
• Fala rozchodzi się pod dużymi kątami.
• Konstrukcja zawiera zagięcia lub krzywizny.
• Urządzenie zawiera izolowane nanostruktury.
• Rozpraszanie jest ważne.
• Zmiany w polu przejściowym mają znaczenie.
• Potrzebujesz pełnych danych dotyczących rozkładu pola elektromagnetycznego.

Przykłady obejmują:

• Photonic crystal waveguides
• Bragg gratings
• Mach-Zehnder interferometers
• Ring resonators
• Tilted facet waveguides
• Tapered waveguides and mode expansion
• Rozpraszanie nanocząstek
• Lokalne defekty
• Dowolne nanostruktury fotoniczne

W tego typu systemach uproszczenia stosowane przez metodę TMM przestają być fizycznie uzasadnione, dlatego konieczne staje się zastosowanie bardziej złożonych obliczeń w metodzie FDTD.

7. Podsumowanie: rozwiązanie najprostszego, poprawnego zadania.

Jedną z najważniejszych umiejętności w fizyce obliczeniowej jest nauczenie się, jak dobierać odpowiedni poziom złożoności modelu.

Bardziej złożone symulacje niekoniecznie oznaczają lepsze symulacje. W wielu przypadkach prostszy model może zapewnić:

• Lepsze zrozumienie zasad fizycznych.
• Szybsze iteracje/powtarzanie cykli.
• Łatwiejsze debugowanie.
• Lepsza optymalizacja parametrów
• Bardziej wiarygodna interpretacja.

Metoda macierzy przenoszenia jest skuteczna właśnie dlatego, że eliminuje zbędną złożoność w przypadku problemów optycznych związanych z warstwowymi układami. Wykorzystuje ona inteligentne uproszczenia matematyczne, co pozwala na niezwykle efektywne rozwiązanie podstawowych zagadnień związanych z interferencją fal.

Metoda FDTD staje się użyteczna, gdy pomijane dotąd aspekty fizyczne zaczynają odgrywać istotną rolę. Wówczas dodatkowe koszty obliczeniowe są uzasadnione, ponieważ prostsze przybliżenia przestają być wystarczające.

W praktyce doświadczeni badacze często stosują obie metody łącznie: TMM do szybkiego badania i optymalizacji, a następnie FDTD do szczegółowej, wielowymiarowej weryfikacji i analizy pola.

Celem nie jest zatem wybór „najnowocześniejszej” metody symulacji, lecz raczej wybór metody, która najlepiej odpowiada zasadom fizyki, na których opiera się rozwiązywane zagadnienie.