Symulacja morfologii struktur heterospójeniowych za pomocą dwuwymiarowych przekrojów elektrycznych.

1. Wprowadzenie

Struktury oparte na heterospójeniach (BHJ) stanowią podstawowy element wielu organicznych urządzeń elektronicznych – najbardziej znane to organiczne ogniwa słoneczne, ale także fotodetektory, czujniki i inne mieszane systemy jonowo-elektroniczne, w których istotne są ścieżki przewodzenia i lokalne ograniczenia. W strukturze BHJ dwie fazy przenikają się w skali 10–100 nm, tworząc duży wewnętrzny interfejs i krajobraz transportu, który zasadniczo nie jest jednowymiarowy. Kompromisem jest złożoność: po zdefiniowaniu mikrostruktury nieuchronnie pojawiają się ograniczenia związane z przewodzeniem, ślepe zaułki, zwężenia i znaczne przestrzenne różnice w gęstości nośników i rekombinacji. W większości przypadków, przy codziennym dostosowywaniu parametrów urządzeń, nie zaczęto by od tego. Zazwyczaj pracowano by na poziomie jednowymiarowego urządzenia – szybko, efektywnie i często w pełni odpowiednio – ewentualnie uwzględniając dodatkowe zjawiska, takie jak pułapkowanie i rekombinacja. Na przykład, przewodnik szybkiego uruchomienia OPV demonstruje standardowy jednowymiarowy schemat działania (w tym fizykę pułapek, jeśli to konieczne), który jest odpowiednim narzędziem dla większości rutynowych badań i ekstrakcji parametrów.

Na drugim końcu hierarchii abstrakcji, jeśli Twoim celem jest analiza zagadnień na poziomie całego systemu, takich jak skalowalność, straty szeregowe/równoległe lub zachowanie modułu/panelu, zazwyczaj przechodzisz do obejmujący duży obszar / uwzględniający strukturę obwodu opisanie, a nie analiza nanostruktury. szczegółowy przewodnik dotyczący PM6:Y6 o dużej powierzchni stanowi przykład bardziej zaawansowanego podejścia do modelowania.

Niniejszy poradnik celowo idzie w przeciwną stronę: zachowuje mikrostrukturę (poprzez określoną morfologię BHJ), ale pozostaje praktyczny obliczeniowo, rozwiązując problem elektryczny na dwuwymiarowych przekrojach. Chodzi o to, że nie każde urządzenie powinno być zaprojektowane w ten sposób, ale że posiadanie metody analizy, która uwzględnia mikrostrukturę, jest wartościowe z punktu widzenia naukowego i technologicznego. Pozwala to na sprawdzenie założeń dotyczących efektywnego medium, zidentyfikowanie, gdzie „proste” modele przestają działać, oraz na rozwinięcie intuicji dotyczącej ograniczeń wydajności wynikających z morfologii – szczególnie w przypadkach, gdy dominują procesy przewodzenia i lokalnej absorpcji/rekombinacji.

2. Utwórz model symulujący zmiany morfologiczne i przeanalizuj strukturę BHJ.

W głównym oknie kliknij „Nowa symulacja”. Spowoduje to otwarcie biblioteki typów urządzeń, jak pokazano na ilustracji ??. Kliknij dwukrotnie „Morfologia”, a następnie zapisz projekt w wybranym miejscu na dysku lokalnym. Zapisywanie lokalnie jest ważne, ponieważ symulacje uwzględniające morfologię mogą generować wiele plików pośrednich (siatki, migawki, pliki pamięci podręcznej), a lokalna pamięć jest zazwyczaj znacznie szybsza i bardziej niezawodna niż foldery synchronizowane w sieci lub w chmurze.

Po utworzeniu projektu otwiera się główne okno urządzenia. Powinien być widoczny centralny, czerwony element, reprezentujący jedną fazę struktury wielozłączowej (BHJ), umieszczony wewnątrz otaczającej go aktywnej warstwy, jak pokazano na ilustracji. ??Chociaż pełny widok 3D jest pouczający, wizualizacja morfologii może być zbyt skomplikowana. Aby lepiej zobrazować strukturę wewnętrzną, kliknij prawym przyciskiem myszy obiekt BHJ i wybierz Wyświetl / Obejrzyj → Pokaż przekrój wzdłuż osi Y..

Włączenie opcji wyświetlania przekroju powoduje wygenerowanie obrazu przedstawionego na ??, który pokazuje dwuwymiarowy przekrój przez trójwymiarową strukturę. Dzięki temu znacznie łatwiej można obserwować połączenia faz, powierzchnie styku i ścieżki przepływu w strukturze BHJ.

Na tym etapie przydatny jest następujący model koncepcyjny: otaczający obszar w kształcie prostokąta określa objętość aktywnej warstwy, natomiast osadzony w nim czerwony obiekt definiuje jedną konkretną fazę mieszaniny heterojunkcyjnej (np. polimer). Uzupełniająca faza jest definiowana pośrednio jako pozostała część objętości aktywnej warstwy, która nie jest zajęta przez osadzony obiekt. To rozróżnienie staje się wyraźne podczas analizy właściwości obiektów w edytorze obiektów, co zostanie omówione w kolejnych sekcjach.

3. Budowa urządzenia: warstwa aktywna i osadzona w niej morfologia BHJ.

Otwórz zakładkę „Struktura urządzenia” i kliknij „Edytor warstw”. Urządzenie składa się z trzech warstw: górnej warstwy kontaktowej, pojedynczej warstwy aktywnej i dolnej warstwy kontaktowej (rycina ??). Ta celowo prosta struktura ułatwia zrozumienie, w jaki sposób morfologia wpływa na problem elektryczny, bez wprowadzania dodatkowych interfejsów lub warstw transportowych. Kluczowe jest to, że sama warstwa aktywna jest prostą, prostokątną objętością zdefiniowaną w edytorze warstw. W tym przykładzie przypisano jej materiał molekularny Y6, który pełni funkcję akceptora, a sama w sobie jest tylko obszarem. Cała złożoność strukturalna związana z heterospójeniem objętościowym jest wprowadzana oddzielnie, poprzez umieszczenie drugiego obiektu wewnątrz tej objętości warstwy aktywnej. Osadzony obiekt to złożona siatka CAD zaimportowana z bazy danych Shape, która reprezentuje fazę polimeru PM6. Siatka polimerowa definiuje zamkniętą, trójwymiarową objętość wewnątrz warstwy aktywnej. Uzupełniająca faza (obszar bogaty w akceptor) jest następnie definiowana pośrednio jako część warstwy aktywnej, która nie jest zajęta przez siatkę polimerową. Innymi słowy, heterospójenie objętościowe jest reprezentowane przez wykluczenie objętościowe, a nie przez wyraźne zdefiniowanie dwóch wzajemnie przeplatających się siatek.

To rozróżnienie staje się wyraźne podczas oglądania obiektów w widoku 3D. Istnieją dwa istotne obiekty: otaczający obszar warstwy aktywnej (zwykle wyświetlany jako niebieska siatka) oraz osadzona w nim struktura polimeru (wyświetlana na czerwono). Aby sprawdzić ich właściwości, kliknij prawym przyciskiem myszy niebieski obszar i wybierz opcję „Edytuj obiekt”, jak pokazano na ilustracji ??.

W edytorze obiektów wyraźnie zaznaczona jest różnica. W przypadku obiektu aktywnego, który stanowi element otaczający, kształt obiektu to prosty sześcian (z wyłączoną opcją złożonych siatek). W przypadku obiektu polimerowego kształt obiektu odnosi się do wpisu w bazie danych kształtów, który zawiera siatkę trójkątów definiującą morfologię PM6. Zmiana tego odwołania powoduje zmianę morfologii polimeru, przy czym otaczająca warstwa pozostaje niezmieniona.

Ta konstrukcja – prosta, aktywna warstwa w formie pudełka, w którą wbudowana jest objętość polimeru – stanowi kluczowy element w sposobie, w jaki OghmaNano zarządza morfologią BHJ. Ułatwia ona analizę geometrii, zapewnia, że faza polimerowa tworzy wyraźnie zdefiniowany, zamknięty obszar oraz umożliwia badanie, w jaki sposób procesy perkolacji, połączenia i zwężenia w sieci polimerowej wpływają na transport ładunku po uruchomieniu solvera elektrycznego.

4. Parametry elektryczne: ruchliwość nośników ładunku i pasma pułapek SRH.

Otwórz edytor parametrów elektrycznych w głównym oknie. Dostępne są oddzielne panele parametrów dla dwóch faz aktywnej warstwy: otaczającego obszaru aktywnego (akceptor/matryca) i osadzonej w nim morfologii polimerowej, co przedstawiono na rysunkach ?? i ??. W tym miejscu definiuje się właściwości elektryczne obu wzajemnie przenikających się faz, odróżniając je od siebie.

W obu fazach stosuje się ten sam podstawowy formalizm transportu: dryfowo-dyfuzyjny, z wykorzystaniem statystyki nośników zgodnie z rozkładem Maxwella-Boltzmanna oraz symetrycznych tensorów ruchliwości. Kluczową różnicą jest celowe wykorzystanie różnicy w ruchliwości. W fazie polimerowej ruchliwość dziur jest wysoka, natomiast ruchliwość elektronów jest ograniczona; w uzupełniającej fazie aktywnej sytuacja jest odwrotna. To zapewnia selektywny transport: dziury przemieszczają się przez sieć polimerową, a elektrony przez matrycę.

Ta prosta asymetria wystarcza, aby wygenerować fizycznie zrozumiałe zachowanie mieszaniny heterojunkcyjnej (BHJ). Już przed uruchomieniem symulacji określa ona, czego można się spodziewać w wynikach przestrzennych: ciągłych ścieżek przewodzenia, ślepych zaułków, w których nośniki uwięzione są geometrycznie, oraz lokalnego gromadzenia ładunku w miejscach zwężenia lub w słabo połączonych obszarach struktury. Oprócz transportu swobodnych nośników, w obu fazach aktywowane jest nieliniowe wychwytywanie zgodnie z mechanizmem Shockleya–Reada–Halla (SRH). Jak pokazano na panelach parametrów, gęstość stanów jest reprezentowana za pomocą wykładniczych rozkładów ogona zarówno dla elektronów, jak i dla dziur, przy identycznej gęstości pułapek i energii aktywacji w każdej fazie. Solver śledzi pięć dyskretnych pasm pułapek w każdym punkcie siatki.

Każdy pas pułapek dynamicznie wymienia nośniki z populacją swobodnych nośników poprzez wyraźne parametry wychwytu i emisji (stałe sprzężenia między nośnikami swobodnymi i uwięzionymi, widoczne w edytorze). W rezultacie, rozwiązywany system obejmuje nie tylko gęstość swobodnych elektronów i dziur, ale także wiele populacji uwięzionych nośników w każdym punkcie przestrzennym. To znacząco zwiększa sztywność numeryczną i koszt obliczeniowy, ale jest niezbędne do modelowania dyspersyjnego transportu, rekombinacji wspomaganej przez pułapki oraz realistycznego zachowania w czasie w nieuporządkowanych organicznych półprzewodnikach. Użyte tutaj wartości numeryczne mają służyć jako stabilny, ilustracyjny przykład, a nie skalibrowany model materiału. Po zrozumieniu przebiegu obliczeń, parametry te można doprecyzować, zweryfikować w oparciu o wyniki eksperymentalne lub uprościć, jeśli celem jest izolacja efektów morfologii i perkolacji, a nie fizyki pułapek.

5. Tworzenie siatki.

Zanim uruchomisz symulację, warto sprawdzić siatkę elektryczną, aby upewnić się, że jest poprawna. W głównym oknie otwórz zakładkę „Elektryczność” (rycina ??) i kliknij „Siatka elektryczna”. Spowoduje to otwarcie edytora siatki, jak pokazano na rysunku ??.

W tym przykładzie siatka ma wymiary 40 × 40 w płaszczyźnie xy, co oznacza, że problem elektryczny jest rozwiązywany w dwóch wymiarach. Oś z została celowo pominięta. W programie OghmaNano możliwe są w pełni trójwymiarowe symulacje elektryczne, ale wymagają one rzeczywistej trójwymiarowej geometrii i wiążą się ze znacznym wzrostem kosztów obliczeniowych i złożoności. Tworzenie i interpretacja takich modeli wykracza poza zakres tego samouczka.

Zamiast tego, należy postrzegać to obliczenie jako dwuwymiarowy przekrój przez urządzenie. Przekrój ten przebiega przez środek struktury i zachowuje pełną, boczną mikrostrukturę BHJ w tej płaszczyźnie. Takie podejście pozwala uwzględnić przewodzenie, wąskie gardła i połączenia fazowe, przy jednoczesnym zachowaniu szybkości, czytelności i przydatności do analiz na poziomie ilustracyjnym.

6. Uruchamianie symulacji.

Kliknij niebieski przycisk „Uruchom symulację” w głównym oknie, aby rozpocząć działanie programu. Podczas symulacji w oknie terminala wyświetlane są wartości punktu pracy i gęstości prądu w każdym kroku (??). W tym przykładzie napięcie przyłożone do górnego kontaktu zmienia się stopniowo od 0,00 V do około 0,54 V, co pozwala na uzyskanie charakterystyki prądowo-napięciowej. Przy niskim napięciu zgłaszana gęstość prądu jest ujemna podczas oświetlenia, a gdy napięcie zbliża się do stanu obwodu otwartego, wypadkowy prąd przechodzi przez zero i zmienia znak, co jest zgodne z zachowaniem ogniwa słonecznego.

Po zakończeniu symulacji otwórz zakładkę „Wyniki” (??), aby sprawdzić wygenerowane pliki. Podsumowanie w oknie terminala zawiera informacje o całkowitym czasie trwania symulacji, liczbie rozwiązanych równań oraz podstawowe statystyki dotyczące operacji wejścia/wyjścia na dysku. W symulacjach uwzględniających morfologię warto monitorować generowane pliki: zapisywanie dużej liczby migawek może znacząco wydłużyć czas trwania symulacji, nawet jeśli sam algorytm jest szybki.

Aby sprawdzić poprawność geometrii, kliknij dwukrotnie device.csv, aby wyświetlić siatkę trójkątów reprezentującą geometrię urządzenia, która jest wykorzystywana wewnętrznie przez solver (??). Dzięki temu można potwierdzić, że morfologia BHJ została zaimportowana i podzielona na elementy, zgodnie z oczekiwaniami.

Aby zobaczyć parametry elektryczne, otwórz jv_contact_0.csv lub jv_contact_1.csv i wyświetl wykres charakterystyki prądowo-napięciowej (??). W przypadku tej symulacji dwuwymiarowej należy użyć danych dotyczących charakterystyki prądowo-napięciowej dla każdego kontaktu; jv.csv można pominąć.

7.0 Analiza migawek z symulacji.

Po zakończeniu symulacji, najbardziej przydatne wyniki można znaleźć w katalogu snapshots/. Pliki zawierające migawki zapisują wewnętrzne zmienne stanu systemu dryfu–dyfuzji–SRH, przedstawione jako funkcja przyłożonego napięcia. Zawierają one przestrzennie rozdzieloną informację o rozwiązaniu równań transportu dla każdego kroku napięcia, a nie tylko podsumowania uzyskane w wyniku dalszej obróbki. Po otwarciu pliku migawki dostępne są dwie opcje sterowania: suwak napięcia, który pozwala na zmianę wartości napięcia w symulacji, oraz suwak pozycji przekroju, który umożliwia wybranie dwuwymiarowego przekroju przez strukturę. W trakcie tego samouczka problem elektryczny jest rozwiązywany dla dwuwymiarowego przekroju, a wszystkie dane z migawek odnoszą się do tego przekroju.

Na rysunkach ??–?? przedstawiono przykładowe dwuwymiarowe obrazy uzyskane przy stałym napięciu. Dane te należy interpretować łącznie. Układ poziomów energetycznych w paśmie przewodnictwa (LUMO), Ec.csv (??), określa warunki energetyczne dla elektronów. W symulacji uwzględniającej morfologię, ten układ zmienia się w płaszczyźnie ze względu na zależność energetyczną od fazy oraz elektrostatyczną odpowiedź na lokalne nagromadzenie ładunku. Gęstość swobodnych elektronów, Q_nfree.csv (??), wskazuje obszary, w których elektrony mogą efektywnie przemieszczać się przez sieć akceptorów/matrycy. Znaczne zmiany przestrzenne wskazują na ścieżki przewodzenia, zwężenia i morfologicznie uwarunkowane przeszkody.

Gęstość elektronów uwięzionych w pułapkach, Q_ntrap.csv (??), wskazuje, gdzie gromadzą się nośniki w stanach pułapkowych. Porównanie gęstości nośników swobodnych i uwięzionych jest często najszybszym sposobem identyfikacji obszarów o słabej efektywności ekstrakcji lub silnym lokalnym uwięzieniu, co wiąże się z morfologią warstwy BHJ. Parametr związany z mobilnością, mu_p_ft.csv (??), pokazuje efektywną mobilność, uśrednioną dla nośników swobodnych i uwięzionych. Zróżnicowanie przestrzenne odzwierciedla zarówno wewnętrzny kontrast fazowy, jak i lokalne uwięzienie, a często wykazuje silną korelację z połączeniami w warstwie bazowej.

Chociaż dwuwymiarowe widoki często stanowią najczytelniejszy sposób analizy przestrzennych zmian, czasami łatwiej jest interpretować te same dane w przestrzeni trójwymiarowej. Włączenie trybu 3D w oknie podglądu powoduje, że dwuwymiarowy przekrój jest ponownie wyświetlany w geometrii urządzenia. Na rysunkach ??–?? przedstawiono przykłady takiej projekcji. Uśrednione wartości ruchliwości elektronów i dziur (mu_n_ft.csv, mu_p_ft.csv), quasi-poziom Fermiego dla elektronów swobodnych (Fn.csv) oraz gęstość prądu elektronów w kierunku przewodzenia (Jn.csv) są szczególnie przydatne do wizualizacji ścieżek przewodzenia i przepływu prądu przez heterospójnik.

Te przykładowe wyniki stanowią główne źródło wartości naukowej symulacji, które pozwalają na analizę morfologii: umożliwiają one bezpośrednie powiązanie obserwowalnych parametrów na poziomie urządzenia z lokalną energią, procesami wychwytywania i połączeniami wewnątrz heterospójnika.

8. Zadania, które warto wypróbować.

Najszybszym sposobem na rozwijanie intuicji jest wprowadzanie zmian w modelu i obserwowanie, co się dzieje. Poniższe zadania zostały zaprojektowane tak, aby można je było szybko wykonać: zmień siatkę morfologiczną, uruchom model ponownie i porównaj krzywą prądu i napięcia (a jeśli masz czas, porównaj również kilka migawek przy tym samym napięciu polaryzacyjnym).

Zadanie 1 — zamień siatkę morfologiczną.

W głównym widoku urządzenia kliknij prawym przyciskiem myszy czerwony obiekt BHJ (osadzoną siatkę polimerową) i wybierz opcję „Edytor siatki”. W edytorze upewnij się, że aktualnie wybrana jest opcja „Kształt z bazy danych”. Zmień wybrany kształt na „morfologia/1”, „morfologia/2” i „morfologia/3” (rysunki ??, ??, ??). Po każdej zmianie uruchom ponownie symulację i porównaj krzywą JV.

Chodzi o to, jak połączenia i wąskie gardła wpływają na proces ekstrakcji i rekombinacji: niektóre struktury tworzą bardziej przejrzyste ścieżki przewodzenia (wyższy prąd), podczas gdy inne zmuszają nośniki do przemieszczania się przez wąskie przejścia i ślepe zaułki (większa akumulacja, większe zatrzymywanie i zazwyczaj gorsze parametry ogniwa słonecznego).

Jako celowe, choć „niepoprawne, ale pouczające” ćwiczenie, w bazie danych znajdziesz również kształt o nazwie „czajnik”. Zastąp siatkę BHJ kształtem czajnika, uruchom symulację ponownie i wygeneruj krzywą JV dla urządzenia z wkładką w kształcie czajnika. To nie jest fizycznie poprawna siatka BHJ, ale pokazuje, że proces obliczeniowy jest rzeczywiście oparty na geometrii: dowolną zamkniętą siatkę CAD można wprowadzić i poddać analizie elektrycznej.

Zadanie 2 — zmiana siły pułapkowania (pasma pułapek SRH).

W sekcji „Parametry elektryczne” dostosuj ustawienia pułapek SRH, które nie znajdują się w stanie równowagi, i uruchom symulację ponownie. Wykonuj zmiany pojedynczo:

• Zwiększ lub zmniejsz gęstość pułapek o 1–2 rzędy wielkości i porównaj Q_nfree z Q_ntrap.

Celem tego zadania jest sprawdzenie, w jaki sposób transport ograniczony przez rozpraszanie lub pułapkowanie konkuruje z morfologią: pułapkowanie ma tendencję do wzmacniania „niekorzystnych” obszarów (ślepych zaułków, zwężeń), ponieważ ładunek nie może się z nich szybko rozładować.

Zadanie 3 — zmień wartości stałych szybkości rekombinacji.

W sekcji „Parametry elektryczne” zmodyfikuj ustawienia rekombinacji i uruchom symulację ponownie. Dwa proste parametry do zmiany:

• Zmień wartość stałej szybkości rekombinacji w procesie swobodnej rekombinacji (od n_free do p_free), zwiększając ją lub zmniejszając o 1–2 rzędy wielkości. Obserwuj, jak zmienia się krzywa charakterystyki prądowo-napięciowej (JV) oraz przestrzenny rozkład ładunku.
• Jeśli w ramach projektu analizowane są parametry związane z wychwytywaniem/łączeniem sygnałów SRH (stałe sprzężenia wolnych nośników z nośnikami uwięzionymi), należy odpowiednio zwiększyć lub zmniejszyć te parametry i porównać wartości Q_nfree i Q_ntrap przy tym samym napięciu polaryzacyjnym.

Przydatnym nawykiem podczas wykonywania zadań 2 i 3 jest utrzymywanie stałej morfologii, zmienianie jednego parametru, a następnie za każdym razem sprawdzanie tych samych dwóch aspektów: (i) jak zmienia się krzywa zależności prądu od napięcia (JV), oraz (ii) czy ścieżki przepływu prądu (np. Jn i pola quasi-poziomu Fermiego) stają się bardziej gładkie, czy też bardziej „rozdrobnione” pod wpływem mikrostruktury.