Strona główna Przykłady Zrzuty ekranu Instrukcja obsługi Logo Bluesky. YouTube ☰
OghmaNano Platforma do symulacji wielofizycznych procesów zachodzących w optoelektronicznych urządzeniach i systemach fotonicznych. POBIERZ Przewodnik szybkiego uruchomienia

Samouczek dotyczący fotonicznych kryształów: symulacja fotonicznego kryształu i światłowodu.

Scanning electron microscope style image of a two-dimensional photonic crystal consisting of a periodic array of air holes etched into a dielectric slab, illustrating the type of structure commonly used to control and guide light
Obraz fotonicznej struktury krystalicznej uzyskany za pomocą skaningowego mikroskopu elektronowego (SEM). Okresowa sieć otworów tworzy przestrzennie zmienny współczynnik załamania, który można wykorzystać do skupiania, kierowania i manipulowania światłem. Fotoniczne struktury krystaliczne stanowią podstawę wielu zintegrowanych urządzeń fotonicznych, w tym światłowodów, rezonatorów i filtrów optycznych.

1. Wprowadzenie do fotonicznych kryształów i światłowodów

A fotoniczna struktura okresowa jest to struktura optyczna, której współczynnik załamania zmienia się okresowo w skali długości fali światła. Typowym przykładem jest płytka dielektryczna, na której ułożono regularną siatkę otworów wypełnionych powietrzem, jak pokazano na ilustracji. ??Ta okresowa struktura tworzy fotoniczna przerwa pasmowazakres długości fal, które nie mogą rozprzestrzeniać się przez kryształ. W wielu aspektach jest to analogiczne do przerwy energetycznej występującej w półprzewodnikach, z tą różnicą, że dotyczy ona fotonów, a nie elektronów. Poprzez wprowadzenie celowego defektu w strukturze, na przykład poprzez usunięcie rzędu „dziur”, staje się możliwe tworzenie modów optycznych wewnątrz przerwy energetycznej. Światło może być wówczas skupiane i kierowane przez obszar defektu, nawet wokół ostrych zakrętów, przy bardzo małych stratach. światłowody o strukturze fotonicznej stanowią podstawę wielu zintegrowanych urządzeń fotonicznych, w tym optycznych połączeń, rezonatorów, filtrów, czujników oraz elementów selektywnych pod względem długości fali.

W tym samouczku użyjesz narzędzia rozwiązywacza FDTD (metoda różnic skończonych w dziedzinie czasu) firmy OghmaNano do zbudowania i symulacji dwuwymiarowej fotonicznej struktury okresowej. Nauczysz się, jak generować okresową sieć za pomocą edytora kroków, jak wprowadzać światło do struktury za pomocą źródeł światła FDTD, jak wizualizować rozkład pola elektromagnetycznego oraz jak obliczać widmo transmisji w funkcji długości fali za pomocą detektorów FDTD i narzędzi do ekstrakcji widma. Zmieniając geometrię sieci, zobaczysz, jak przesuwa się pasmo fotoniczne, a wprowadzając defekt, stworzysz fotoniczny przewodnik światłowodowy i zaobserwujesz propagację światła w strukturze.

Fotoniczne kryształy są powszechnie stosowane we współczesnej zintegrowanej fotonice, ponieważ umożliwiają precyzyjną kontrolę nad przepływem światła, zajmując przy tym niewielką powierzchnię na układzie scalonym. Zrozumienie, w jaki sposób powstają fotoniczne przerwy pasmowe i w jaki sposób można wykorzystać defekty do manipulowania modami optycznymi, jest zatem ważnym krokiem w kierunku projektowania zaawansowanych urządzeń fotonicznych.

2. Tworzenie nowego modelu symulacji fotonicznej struktury okresowej

Uruchom program OghmaNano z menu Start systemu Windows, a następnie kliknij „Nowa symulacja”, aby otworzyć bibliotekę dostępnych typów urządzeń, którą można zobaczyć w ??.

W nowym oknie symulacji (??) kliknij dwukrotnie przykłady FDTD. Otworzy się lista szablonów optycznych/FDTD, przedstawiona na ilustracji ??, która zawiera rezonator Fabry’ego-Perota, przestrzeń, zintegrowany rezonator pierścieniowy fotoniki, modulator Macha-Zehndera oraz kilka światłowodów. W ramach tego samouczka kliknij dwukrotnie „Kryształ fotoniczny”. Po wyświetleniu odpowiedniego komunikatu zapisz symulację w folderze, do którego masz uprawnienia zapisu.

💡 Wskazówka: Aby uzyskać najlepszą wydajność, zapisz dane na lokalnym dysku, na przykład C:\. Symulacje przechowywane w folderach sieciowych, na dyskach USB lub w chmurze (np. OneDrive) mogą działać wolniej ze względu na dużą liczbę operacji odczytu/zapisu.

OghmaNano new simulation window listing device categories such as organic solar cells, OLEDs, OFETs, perovskites, ray tracing, and FDTD examples
Nowe okno symulacji zawiera bibliotekę typów urządzeń i przykładowych projektów. Kliknij dwukrotnie ikonę, aby otworzyć wstępnie skonfigurowaną symulację. W przypadku kryształów fotonicznych kliknij dwukrotnie ikonę „Przykłady FDTD”, aby wyświetlić dostępne szablony optyczne.
OghmaNano FDTD examples list showing Fabry-Perot cavity, free space, integrated photonics ring resonator, photonic crystal, photonic-xtal trumpet, silicon Mach-Zehnder modulator, and tapered 975 nm waveguide
W folderze „Przykłady FDTD” znajdują się gotowe do uruchomienia modele optycznych struktur. Każdy z nich otwiera wstępnie skonfigurowaną symulację FDTD, którą można modyfikować. Kliknij dwukrotnie „Fotoniczna struktura krystaliczna”, aby załadować model użyty w tym samouczku.

3. Uruchom symulację.

Po zapisaniu otworzy się główne okno (patrz ??). Urządzenie składa się z szeregu dielektrycznych filarów (czerwone rurki), umieszczonych pomiędzy dwiema cienkimi płaszczyznami detektorów: fioletowym detektorem wejściowym, znajdującym się blisko źródła, i zielonym detektorem wyjściowym, umieszczonym po przeciwnej stronie. Użyj przycisków xy/yz/xz, aby zmienić perspektywę. Kliknij „Uruchom symulację” (niebieska ikona odtwarzania) lub naciśnij F9. Na wolniejszych komputerach obliczenia FDTD mogą trwać nieco dłużej.

Po zakończeniu obliczeń otwórz zakładkę „Wyniki” (??), aby przejrzeć pliki, które zostały utworzone przez program. Kliknięcie dwukrotne dowolnego pliku spowoduje jego otwarcie w odpowiednim programie.

Main OghmaNano interface showing the Run Simulation button and a 3D view of a photonic crystal made from a square array of red dielectric pillars, with purple and green detector planes.
Główny OghmaNano interfejs po załadowaniu fotonicznej struktury krystalicznej. Na trójwymiarowym panelu wyświetlana jest siatka dielektrycznych filarów wraz z płaszczyznami wejściowego (fioletowego) i wyjściowego (zielonego) detektora. Naciśnij Uruchom symulację. przycisk (lub) F9) aby rozpocząć obliczenia metodą FDTD.
OghmaNano Output tab while the FDTD simulation is running, showing progress messages, solver initialization and time-step information.
Podczas działania symulacji w zakładce „Wyniki” wyświetlane są informacje o postępie obliczeń wykonywanych przez solver FDTD, w tym dane dotyczące inicjalizacji urządzenia, konfiguracji detektora oraz aktualnego kroku czasowego symulacji.
OghmaNano Output tab displaying the simulation working directory with result files such as detector0, detector1, snapshots, power.dat, and device.csv.
Zakładka „Wyniki” wyświetla katalog roboczy dla symulacji. Najważniejsze elementy to folder snapshots (zawierający zdjęcia pola w różnych momentach czasu) oraz foldery detector0 i detector1, w których przechowywane są widma zarejestrowane przez detektory wejściowe i wyjściowe. Kliknięcie dwukrotne pliku lub folderu powoduje jego otwarcie.

4. Obserwuj, jak pole się rozprzestrzenia.

Kliknij dwukrotnie folder snapshots, aby otworzyć program do wyświetlania danych (??). Wykres przedstawia plik Ey.csv – jeden z elementów pola elektrycznego – w postaci kolorowej siatki, gdzie kolor i wysokość reprezentują natężenie pola w V/m.

Przesuwaj suwak znajdujący się pod wykresem, aby obserwować zmiany w czasie. Zobaczysz, jak emitowany impuls rozchodzi się od źródła po lewej stronie, dociera do fotonicznej struktury krystalicznej i częściowo się w niej rozprasza oraz przechodzi przez nią. Obserwacja zmian pola jest najbardziej intuicyjnym sposobem na potwierdzenie, że symulacja działa zgodnie z oczekiwaniami, zanim zaczniesz ufać wynikom.

OghmaNano simulation snapshot viewer showing the Ey component of the electric field as a coloured wireframe surface, with a time slider beneath the plot.
Program do przeglądania obrazów wyświetla Ey.csv w formie siatki linii.
Early-time snapshot of the electric field interacting with the photonic crystal.
Wczesna faza eksperymentu, na której widać, jak impuls światła wchodzi do fotonicznej struktury krystalicznej.
Later-time snapshot showing the electric field after propagating through the photonic crystal.
Późniejsze zdjęcie przedstawia pole po przejściu przez kryształ.

5. Detektory i widmo transmisji

Teraz przyjrzyjmy się danym zarejestrowanym przez oba detektory. W zakładce „Wyniki” otwórz foldery detector0 i detector1; w każdym z nich znajduje się plik o nazwie lam_E.csv, który zawiera informacje o natężeniu pola |E|2 w funkcji długości fali, rejestrowane przez dany detektor. Oba detektory to fioletowa i zielona siatka, którą można zobaczyć w widoku urządzenia.

Otwarcie lam_E.csv w detector0 generuje widmo wejściowe (??): jest to szerokopasmowy impuls wprowadzany do struktury, przy czym większość jego energii skupia się w okolicach 900 nm. Otwarcie lam_E.csv w detector1 generuje widmo wyjściowe (??), które przedstawiono na osi logarytmicznej. Zwróć uwagę, jak silne tłumienie światła występuje dla większości długości fal – fotoniczna struktura je blokuje.

Porównywanie dwóch surowych widm gołym okiem jest niewygodne, dlatego folder detektora zawiera również lam_norm.csv, czyli stosunek |EOutput|2/|EInput|2. Jest to widmo transmisji fotonicznej struktury krystalicznej. Kliknij go dwukrotnie, aby uzyskać ??.

Wyraźny szczyt nieco powyżej 1000 nm odpowiada zakresowi długości fal, które fotoniczna struktura przepuszcza; silnie tłumione obszary po obu stronach to fotoniczna przerwa energetyczna, w której propagacja jest niemożliwa. Innymi słowy, ta pojedyncza krzywa informuje nas zarówno o tym, jakie długości fal przepuszcza kryształ, jak i o położeniu jego fotonicznej przerwy energetycznej.

Input detector spectrum lam_E.csv showing |E| squared versus wavelength on a linear axis, with a broad peak near 900 nm.
Detektor wejściowy (detector0/lam_E.csv): szerokopasmowy impuls padający na kryształ, o maksymalnej amplitudzie w okolicach 900 nm. Służy on jako punkt odniesienia, względem którego porównywany jest sygnał wyjściowy.
Output detector spectrum lam_E.csv showing |E| squared versus wavelength on a logarithmic axis, strongly attenuated across most of the band.
Detektor sygnału wyjściowego (detector1/lam_E.csv) przedstawiony w skali logarytmicznej. Większość długości fal jest tłumiona o kilka rzędów wielkości, ponieważ mieści się ona w zakresie pasma fotonicznego.
Transmission spectrum, the ratio of output to input field intensity, showing a single sharp peak just above 1000 nm.
Widmo transmisji |EOutput|2/|EInput|2 z lam_norm.csv. Maksimum nieco powyżej 1000 nm wskazuje długości fal, które są transmitowane przez fotoniczną strukturę; tłumione długości fal znajdują się wewnątrz pasma zabronionego fotonicznego.

6. Dostosowanie promienia filara.

Położenie przerw energetycznych jest określone przez geometrię sieci krystalicznej. Aby to zobaczyć, kliknij prawym przyciskiem myszy na strukturze fotonicznej w widoku urządzenia i wybierz z menu opcję „Edytor siatki” (??). Edytor siatki (??) umożliwia zmianę kształtu, który jest powtarzany w celu utworzenia struktury – w tym przypadku jest to rura o promieniu 100 nm, wysokości 400 nm oraz określonej liczbie pierścieni i segmentów, które kontrolują gęstość siatki.

Right-click context menu on the photonic crystal with the Mesh editor option highlighted.
Kliknięcie prawym przyciskiem myszy na strukturze fotonicznej powoduje wyświetlenie menu kontekstowego. Wybierz opcję „Edytor siatki”, aby zmienić geometrię powtarzającego się obiektu.
Mesh editor showing a Tube primitive with radius 100 nm, height 400 nm, 10 rings and 32 segments.
Edytor siatki, w którym zaznaczono prymitywny element „rura”. Pole „Promień” (tutaj 100 nm) określa grubość każdego słupka w siatce.
Device view of the photonic crystal after increasing the pillar radius to 150 nm, showing visibly fatter pillars.
Ta sama fotoniczna struktura krystaliczna, po zwiększeniu promienia filarów do 150 nm, charakteryzuje się wyraźnie grubszymi filarami, dzięki czemu w strukturze kratowej mieści się więcej materiału dielektrycznego.

Zmień promień filarów na 150 nm. W widoku urządzenia filary staną się wyraźnie grubsze (??). Uruchom ponownie symulację i ponownie otwórz lam_norm.csv, aby wygenerować zaktualizowane spektrum transmisji (??).

Następnie ustaw promień na 120 nm i powtórz symulację. Otrzymane spektrum transmisji (??) wykazuje teraz maksimum w pobliżu 1200 nm, co znajduje się pomiędzy początkową wartością 100 nm a wartością 150 nm. Widoczna jest wyraźna tendencja: zwiększenie promienia filara przesuwa pasma fotoniczne w kierunku dłuższych długości fal. Poprzez zmianę jednego parametru geometrycznego można regulować położenie przerwy pasmowej i długości fal transmitowanych przez kryształ.

Transmission spectrum for 150 nm pillars, with the main transmission feature shifted to roughly 1400 to 1500 nm.
Przesunięcie widma dla struktur o wysokości 150 nm. Maksimum transmisji przesunęło się do zakresu około 1400–1500 nm, co stanowi znaczące przesunięcie w kierunku dłuższych fal w porównaniu ze strukturą o wysokości 100 nm.
Transmission spectrum for 120 nm pillars, with the main peak near 1200 nm.
Przesunięcie widma dla struktur o wysokości 120 nm. Maksimum w pobliżu 1200 nm znajduje się pomiędzy wynikami uzyskanymi dla struktur o wysokości 100 nm i 150 nm, co pokazuje płynną zmianę pasm fotonicznych wraz ze wzrostem promienia struktury.

7. Tworzenie fotonicznego kryształowego światłowodu.

Do tej pory wygenerowano dla Ciebie strukturę siatki. Aby zobaczyć, jak to zostało zrobione, i móc ją modyfikować, kliknij prawym przyciskiem myszy na strukturze fotonicznej i wybierz. Edytuj obiekt (??). Edytor obiektów (??) zbiera wszystkie informacje na temat danego obiektu: jego położenie, kształt rury, materiał, z którego jest wykonany (w tym przypadku krzem), oraz – co ważne – Kroki / Etapy / Instrukcje wiersz — zasada, która definiuje kształt i tworzy z niego siatkę. Kliknij przycisk. … naciśnij przycisk w tym wierszu, aby otworzyć edytor krokówgdzie można tworzyć i modyfikować kryształowe struktury, powtarzalne wzory oraz precyzyjnie zdefiniowane kształty geometryczne.

Right-click context menu on the photonic crystal with the Edit object option highlighted.
Kliknij prawym przyciskiem myszy na strukturze fotonicznej i wybierz opcję „Edytuj obiekt”, aby otworzyć pełny edytor obiektów.
Object editor showing object type, offset, the Steps row labelled custom crystal lattice 10x1x10, colour, tube shape, and silicon optical material.
Edytor obiektów. W wierszu „Kroki” widnieje napis „Niestandardowa sieć krystaliczna 10x1x10” – to receptura używana do powtarzania rury i tworzenia kryształu. Kliknij przycisk z trzema kropkami po prawej stronie wiersza, aby go edytować.
Crystal editor showing lattice vectors a1, a2, a3, an origin, and a number of objects of 10 by 1 by 10.
Zakładka „Kryształ” w edytorze etapów. Wektory sieci a1, a2 i a3, wraz z liczbą obiektów, definiują trójkątną sieć, która jest wykorzystywana do generowania fotonicznej struktury krystalicznej.

Kliknięcie przycisku „…” w wierszu „Kroki” otwiera edytor kroków. Na karcie „Kryształ” (??) sieć krystaliczna jest definiowana przez trzy wektory i liczbę. Tutaj a1 = (400 nm, 0, 0) określa odstęp 400 nm wzdłuż osi x, a2 = (0, 1 µm, 0) to pojedyncza warstwa wzdłuż osi y, a a3 = (200 nm, 0, 346 nm) przesuwa każdy wiersz o połowę odstępu – przesunięcie o 200 nm przy odstępie 346 nm = 400 nm × √3/2 – co daje znane, trójkątne (heksagonalne) ułożenie zwartych warstw. Kryształ jest powtarzany 10 × 1 × 10 razy.

Lua script editor generating the photonic crystal lattice with clear_points and add_point inside two nested for loops, with conditions that skip selected pillars.
Zakładka skryptu Lua. Zamiast standardowej siatki, położenia są generowane w kodzie, dzięki czemu możliwe jest dowolne rozmieszczenie. Dwa kluczowe polecenia to clear_points() i add_point(), a wszystkie współrzędne są podawane w metrach.
Device view of a photonic crystal with a U-shaped channel of missing pillars forming a photonic crystal waveguide.
Fotoniczny przewodnik falowodowy. Linia usuniętych filarów tworzy kanał w kształcie litery U, stanowiący defekt. Światło o długości fali znajdującej się w zakresie pasma fotonicznego jest uwięzione w tym defekcie i kierowane wzdłuż zakrętu.
Transmission spectrum of the photonic crystal waveguide showing several sharp guided-mode peaks.
Przesyłanie sygnału |EOutput|2/|EInput|2 w przewodniku falowym. Liczne, wyraźne piki odpowiadają kierowanym modom w zakrzywionym kanale w kształcie litery U – są to długości fal, które są kierowane wokół zakrętu, podczas gdy reszta spektrum jest blokowana.

Regularna sieć to najprostszy przypadek, ale nie jest to jedyna możliwość. Kliknij przycisk skryptu Lua, a edytor przełączy się w tryb skryptów (??). W tym trybie położenie filarów jest generowane programowo, a nie na podstawie zdefiniowanej wcześniej struktury krystalicznej. Wystarczą tylko dwa polecenia: clear_points() usuwa wszystkie istniejące położenia filarów, a add_point(x, y, z) dodaje nowy filar w określonych współrzędnych. W połączeniu z pętlami i instrukcjami warunkowymi, te polecenia umożliwiają generowanie niemal dowolnej geometrii fotonicznej struktury krystalicznej.

Przedstawiony w ?? przykładowy skrypt tworzy siatkę o wymiarach 13 × 10, z odstępem między elementami wynoszącym 400 nm. Niewielka liczba instrukcji if ustawia wartość skip=true dla wybranych filarów, uniemożliwiając ich dodanie do struktury. W rezultacie powstaje linia brakujących filarów, tworząca kanał defektu w strukturze kryształu. W tym przykładzie defekt ma kształt litery U, tworząc fotoniczny przewodnik falowy, jak pokazano na ilustracji ??.

Uruchom ponownie symulację i otwórz ponownie lam_norm.csv. Widmo transmisji światłowodu (??) zawiera teraz kilka wąskich pików transmisji, a nie pojedynczą szeroką wiązkę. Piki te odpowiadają falowodom, które są generowane w kanale defektu wewnątrz przerw fotonicznych. Światło o tych długościach fal jest skupione wewnątrz światłowodu i skutecznie kierowane wokół zakrętu, natomiast długości fal spoza zakresu falowodów są blokowane przez otaczający kryształ.

Gratulacje! Przeprowadziłeś symulację metodą FDTD dla fotonicznej struktury okresowej, zmierzyłeś jej widmo transmisji, zbadałeś, jak geometria wpływa na położenie pasm fotonicznych, oraz stworzyłeś funkcjonalny fotoniczny przewodnik falowodowy, wprowadzając defekt do sieci.

8. Wyniki symulacji FDTD.

Każde uruchomienie symulacji FDTD generuje zestaw danych, które przedstawiają różne aspekty zachowania optycznego – od chwilowych obrazów pola elektromagnetycznego, przez spektrum zarejestrowane przez każdy detektor, po geometrię samego urządzenia. Pliki te są zazwyczaj prostymi plikami CSV, które można otworzyć bezpośrednio w wbudowanych przeglądarkach programu OghmaNano lub przetwarzać w inny sposób (np. poprzez tworzenie wykresów w programach Excel lub Python). Najważniejsze dane uzyskane w ramach tego badania fotonicznej struktury krystalicznej zostały podsumowane w tabeli 1 poniżej.

Table 1: Files produced by the photonic crystal FDTD simulation
Nazwa pliku Opis
detektor0/lam_E.csvWidmo wejściowe |E|2 w funkcji długości fali (fioletowy detektor).
detektor1/lam_E.csvWidmo wyjściowe |E|² w funkcji długości fali (zielony detektor).
detektorN/lam_norm.csvWspółczynnik transmisji |Eout|²/|Ein|² w funkcji długości fali.
zdjęcia/Przykładowe dane z terenu w określonych momentach czasowych (np. Ey.csv); zob. ??.
plik device.csvModel 3D / model geometryczny
escape.csvEnergia opuszczająca obszar symulacji.
power.datZależność natężenia pola od czasu.
plik konwersji.datZbieżność metody FDTD.
plik gui_config.jsonZapisane ustawienia widoku/układu; więcej informacji w ??.