Samouczek dotyczący domeny ekscytonów (część A): uruchom symulację domeny ekscytonów w 3D.

1. Wprowadzenie

Zrozumienie procesów powstawania, transportu i rozpadu ekscytonów ma kluczowe znaczenie dla działania organicznych ogniw słonecznych z heterospójeniem objętościowym (BHJ). W takich układach światło jest absorbowane głównie w fazie donora, co prowadzi do powstawania związanych par elektron-dziura, znanych jako ekscytony. Ekscytony te muszą dyfundować przez materiał donora i dotrzeć do granicy między donorem a akceptorem, zanim ulegną rozpadowi. Na granicy mogą one rozpaść się na swobodne elektrony i dziury, przy czym prawdopodobieństwo tego procesu zależy od lokalnej morfologii, wymiarowości oraz odpowiednich stałych kinetycznych.

Równowaga między dyfuzją ekscytonów, dysocjacją na granicy faz i konkurującymi procesami strat odgrywa zatem decydującą rolę w określaniu efektywnej wydajności fotogeneracji. Jednak w przeszłości unikało się szczegółowego modelowania ekscytonów nie ze względu na trudności koncepcyjne, lecz z powodu niepewności dotyczącej parametrów. Długości dyfuzji ekscytonów, czasy życia, szybkości dysocjacji na granicy faz, kanały radiacyjnej i nieradiacyjnej rekombinacji oraz procesy anihilacji były często słabo określone dla danego systemu materiałowego. W wielu praktycznych procesach modelowania efekty te były zatem uwzględniane w jednym skalarnym współczynniku wydajności fotonów, \(\eta_{\mathrm{photon}}\), który reprezentował sumaryczne straty, bez uwzględniania podstawowych procesów transportu i kinetyki.

Obecnie sytuacja się zmienia. Postępy w metodach eksperymentalnej charakterystyki pozwalają na coraz dokładniejsze pomiary czasu życia ekscytonów, długości dyfuzji i kanałów strat w systemach podobnych do BHJ. Wraz z udoskonalaniem tych parametrów, szczegółowe modelowanie domen ekscytonowych staje się coraz bardziej użyteczne – nie jako narzędzie do przewidywania bezwzględnej wydajności urządzeń, ale jako sposób na zbadanie, w jaki sposób rozmiar domen, ich wymiarowość i kinetyka materiału wpływają łącznie na efektywność generowania ładunków. Model domen ekscytonowych, wykorzystany w tym samouczku, został opracowany i opublikowany w ramach tych najnowszych badań eksperymentalnych i modelowych (patrz np. Nature Materials 21, 55–61 (2022)).

W niniejszym poradniku przyjmujemy zatem idealizowaną, ale w pełni trójwymiarową geometrię elementu podstawowego: obszar donora osadzony w matrycy akceptora, początkowo przedstawiony jako sfera donora wewnątrz sześcianu akceptora. Chociaż geometria jest uproszczona, model uwzględnia w sposób wyraźny generowanie, dyfuzję, dysocjację na granicy faz oraz konkurencyjne procesy strat eksytonów. Dzięki temu stanowi on praktyczną i zrozumiałą z punktu widzenia fizyki podstawę do badania, w jaki sposób mierzone eksperymentalnie parametry oddziałują ze sobą, do sprawdzania ich wewnętrznej spójności oraz do rozwijania intuicji na temat tego, jak morfologia i kinetyka wspólnie wpływają na efektywność fotogeneracji w trzech wymiarach.

2. Równania rządzące.

W tym samouczku transport ekscytonów jest modelowany w sposób bezpośredni poprzez rozwiązywanie trójwymiarowego równania dyfuzji ekscytonów w całym obszarze symulacji. Po włączeniu modelu ekscytonów, absorpcja światła bezpośrednio wpływa na populację ekscytonów, która następnie zmienia się zgodnie z następującym równaniem:

\[ \frac{\partial X}{\partial t} = \nabla \cdot \left( D \nabla X \right) + G_{\mathrm{optical}} - k_{\mathrm{dis}} X - k_{\mathrm{FRET}} X - k_{\mathrm{PL}} X - \alpha X^2. \]

W tym równaniu \(X(\mathbf{r},t)\) oznacza gęstość ekscytonów (\(\mathrm{m^{-3}}\)), a \(D\) – współczynnik dyfuzji ekscytonów (\(\mathrm{m^2\,s^{-1}}\)). Człon źródłowy \(G_{\mathrm{optical}}\) (\(\mathrm{m^{-3}\,s^{-1}}\)) reprezentuje generację ekscytonów w wyniku absorpcji światła. Pozostałe człony opisują konkurencyjne procesy, które powodują zanik ekscytonów:

• Dysocjacja (\(k_{\mathrm{dis}} X\)): przekształcenie ekscytonów w swobodne elektrony i dziury na granicy między donorami a akceptorami. Jest to proces, który ostatecznie przyczynia się do generowania ładunków.
• Transfer Förstera (\(k_{\mathrm{FRET}} X\)): niefotonowy transfer energii od wzbudzonych cząsteczek do sąsiednich miejsc lub domen, w wyniku którego energia wzbudzenia jest redystrybuowana bez bezpośredniego generowania nośników swobodnych.
• Straty radiacyjne (\(k_{\mathrm{PL}} X\)): rekombinacja ekscytonów poprzez emisję fotonów (fotoluminescencja), stanowiąca wewnętrzną ścieżkę strat.
• Anihilacja (\(\alpha X^2\)): bimolekularna anihilacja ekscytonów zachodząca przy wysokich stężeniach, która staje się istotna przy silnym wzbudzeniu lub w ograniczonych przestrzeniach.

Współczynnik dyfuzji \(D\) jest określany za pomocą długości dyfuzji ekscytonu \(L\) (m) i czasu życia ekscytonu \(\tau\) (s) w następujący sposób:

\[ D = \frac{L^2}{\tau}. \]

W tym samouczku przeprowadzamy symulację, w której uwzględniamy jedynie eksytony: solver oblicza trójwymiarową gęstość eksytonów \(X(\mathbf{r})\) oraz powiązane z nią składniki reakcji, ale nie rozwiązuje równań dotyczących transportu elektronów i dziur. W związku z tym składnik dysocjacji \(k_{\mathrm{dis}} X\) interpretowany jest jako przestrzenna funkcja opisująca szybkość generowania ładunku, wskazująca, gdzie i z jaką efektywnością eksytony przekształcałyby się w swobodne nośniki ładunku na granicach między donorami a akceptorami.

3. Utwórz model symulujący domenę eksytonową.

Otwórz okno nowej symulacji (rycina ??). Można to zrobić, klikając opcję „Nowa symulacja” w głównym oknie. Jeśli klikniesz dwukrotnie opcję „Symulacje ekscytonów”, wyświetli się biblioteka ekscytonów (rycina ??). Kliknij dwukrotnie „Domena ekscytonów”, aby otworzyć przykładowy projekt.

4. Sprawdź geometrię i parametry.

Po otwarciu przykładu w głównym oknie wyświetlana jest prosta, trójwymiarowa scena: kula dawcy umieszczona wewnątrz pudełka akceptora (Rysunek ??). Zielone promienie wskazują na oświetlenie padające z góry. Aby sprawdzić parametry ekscytonów, kliknij opcję „Parametry elektryczne” w lewym panelu (w zakładce „Struktura urządzenia”). Otworzy to edytor parametrów elektrycznych. Jeśli przewiniesz do dołu tego okna, znajdziesz sekcję oznaczoną jako „Ekscytony”, w której parametry specyficzne dla ekscytonów są wymienione oddzielnie dla każdego obiektu w scenie (otaczającej warstwy i umieszczonej w niej kuli), jak pokazano na Rysunkach ?? i ??.

Sekcja dotycząca ekscytonów definiuje parametry występujące w równaniu dyfuzji i reakcji. Obejmują one średnią drogę swobodną ekscytonu \(L\), która określa współczynnik dyfuzji poprzez \(D = L^{2}/\tau\), oraz czas życia ekscytonu \(\tau\), który kontroluje ogólną skalę czasową transportu w przypadku braku dodatkowych kanałów strat. Radiacyjna rekombinacja jest opisywana przez szybkość fotoluminescencji \(k_{\mathrm{PL}}\), natomiast nieradiacyjne procesy transferu energii są uwzględniane za pomocą współczynnika typu Förstera \(k_{\mathrm{FRET}}\). Przy wyższych gęstościach ekscytonów uwzględnia się bimolekularną anihilację poprzez współczynnik \(\alpha\), który odpowiada członowi \(-\alpha X^{2}\) w równaniu podstawowym. Na koniec, szybkość dysocjacji \(k_{\mathrm{dis}}\) określa, jak efektywnie ekscytony przekształcają się w stany transferu ładunku na granicach donor-akceptor i jest kluczowym parametrem kontrolującym efektywność generowania ładunku w tym modelu.

Kluczową cechą koncepcyjną tego przykładu jest to, że sfera (donor) i warstwa (matryca akceptora) wykorzystują różne mechanizmy związane z ekscytonami. Na powyższych zrzutach ekranu sfera donora zawiera niezerową wartość \(k_{\mathrm{PL}}\) (kanał strat). Otaczający obszar jest skonfigurowany z bardzo dużą wartością \(k_{\mathrm{dis}}\), więc gdy ekscytony docierają do tego obszaru (lub opuszczają go), szybko ulegają rozkładowi. Dzięki temu granica między donorem a akceptorem zachowuje się jak wydajny obszar pochłaniania, tworząc charakterystyczny profil: „wysoka koncentracja w centrum, niska w pobliżu granicy”.

5. Uruchom symulację.

Kliknij niebieski trójkąt „Uruchom symulację”, aby rozpocząć działanie solvera. Wyniki będą wyświetlane w głównym oknie, jak pokazano na rysunku ??. Początkowe wiersze zawierają ogólne informacje o konfiguracji i geometrii symulacji. Kluczowe informacje diagnostyczne pojawiają się, gdy solver eksytonów rozpoczyna działanie i mają postać Exciton f() = … , steps = …. W tym wierszu wyświetlany jest postęp działania trójwymiarowego solvera eksytonów. Wielkość f() to wartość resztowa, która mierzy, jak bardzo aktualne pole gęstości eksytonów odbiega od stanu zgodności; w miarę iteracji solvera wartość ta powinna maleć. Towarzyszący licznik steps wskazuje numer iteracji. W przedstawionym przykładzie wartość resztowa zmniejsza się od początkowej wartości do 2.97 × 104 w kroku 2, po czym solver informuje, że równanie zostało rozwiązane i kończy działanie.

Na typowym, nowoczesnym laptopie wykonanie tego przykładu powinno zająć około 5–10 sekund. Jeśli czas przetwarzania wydłuża się do kilku minut, zazwyczaj oznacza to, że siatka jest zbyt gęsta lub że geometria lub parametry zostały zmodyfikowane w sposób, który znacznie zwiększa liczbę punktów siatki. W takim przypadku warto ponownie sprawdzić ustawienia siatki lub powtórzyć poprzednie kroki konfiguracji przed kontynuowaniem.

6. Wykres przedstawiający wyniki.

Po zakończeniu obliczeń przejdź do zakładki „Wyniki” (rycina ??). Najważniejsze elementy w kontekście tego samouczka to katalog exciton_output/, który zawiera wyniki dotyczące przestrzennego rozkładu ekscytonów, oraz plik exciton_sim_info.json, w którym podsumowano globalne statystyki dotyczące generowania i strat, a który zostanie wykorzystany później w części C. Kliknięcie dwukrotne exciton_output/ wyświetla zawartość przedstawioną na rycinie ??.

Otwórz G.csv, aby zobaczyć szybkość generowania ekscytonów, obliczoną za pomocą modelu optycznego (rycina ??). W oknie wyświetlającym trójwymiarowy wykres możesz obracać sceną za pomocą myszy, a także używać suwaków Z i Y znajdujących się na dole, aby tworzyć przekroje objętości. Jest to często najszybszy sposób na sprawdzenie, czy geometria, kierunek oświetlenia i siatka zachowują się zgodnie z oczekiwaniami.

Następnie otwórz exciton.csv (rycina ??), która przedstawia stabilną gęstość ekscytonów. Powinieneś zaobserwować charakterystyczny wzór dyfuzji do obszaru na granicy faz: wysoką gęstość ekscytonów wewnątrz sfery donora oraz jej zmniejszenie w kierunku granicy między donorem a akceptorem, gdzie dysocjacja jest silna. Użyj suwaków do regulacji przekrojów, aby potwierdzić, że to zmniejszenie gęstości odpowiada geometrii granicy faz, a nie jest artefaktem wynikającym z kąta obserwacji.

Na koniec otwórz Gn.csv (rycina ??), która przedstawia rozkład w przestrzeni szybkości generowania elektronów wynikający z dysocjacji ekscytonów. W tym przykładzie powinieneś zobaczyć wyraźny pierścień (lub warstwę) generacji zlokalizowanej w pobliżu granicy między donorem a akceptorem. Użyj suwaków Z i Y, aby przesuwać płaszczyznę przekroju przez objętość i zbadać, jak ten profil dysocjacji zmienia się w zależności od położenia.

7. Jakie pliki znajdują się w folderze exciton_output/?

W katalogu exciton_output/ znajdują się dwa różne typy plików. Pierwsza grupa zawiera pola obliczone, generowane przez moduł rozwiązujący równania eksytonów, takie jak gęstość eksytonów i współczynniki generacji wynikające z dysocjacji. Druga grupa zawiera mapy parametrów: trójwymiarowe kopie parametrów eksytonów zdefiniowanych w interfejsie graficznym. Są one zapisywane głównie w celu zapewnienia spójności i ułatwienia debugowania, co pozwala sprawdzić, czy zamierzone parametry materiału zostały poprawnie zastosowane w przestrzeni. Tabele 1 i 2 podsumowują te dwie kategorie. Wszystkie wartości podawane są w jednostkach SI. (Niektóre pliki mogą nie być zapisywane, jeśli odpowiadające im podmodele są wyłączone).

Pliki zawierające mapy parametrów, wymienione w tabeli 2, nie są wynikami działania solvera w tradycyjnym sensie. Stanowią one jedynie zapis parametrów ekscytonów zdefiniowanych w sekcji „Ekscytony” interfejsu graficznego, zapisanych na tej samej trójwymiarowej siatce, na której obliczane są pola rozwiązania. Ich celem jest zapewnienie przejrzystego sposobu sprawdzania, czy parametry są przypisane do właściwych obiektów i obszarów domeny, szczególnie podczas pracy ze złożonymi geometriami.

8.1 Interpretacja exciton_sim_info.json (współczynnik efektywności dysocjacji)

Chociaż trójwymiarowe wykresy są niezwykle przydatne do wizualizacji miejsc, w których powstają, zanikają i rozpadają się ekscytony, same w sobie nie dostarczają zwięzłych danych liczbowych. Plik exciton_sim_info.json wypełnia tę lukę, prezentując zintegrowane w przestrzeni wartości generacji i zaniku dla każdego obiektu w symulacji (otaczającej warstwy, osadzonej sfery oraz wszelkich dodatkowych kształtów, które można dodać).

Aby wyświetlić te informacje, otwórz główną zakładkę „Wyniki” i kliknij dwukrotnie exciton_sim_info.json. Plik zawiera zagnieżdżoną strukturę JSON, która dla każdego obiektu zawiera informacje o całkowitej szybkości generowania ekscytonów oraz o całkowitych wartościach związanych z każdym konkurencyjnym mechanizmem strat lub dysocjacji. Te dane pozwalają określić, w jaki sposób zmiany w parametrach materiału, geometrii lub morfologii przekładają się na zmiany w efektywności generowania ładunków.

Główny element segment0 oznacza po prostu początek raportu. W jego obrębie pole name identyfikuje obiekt nadrzędny (w tym przypadku warstwę), natomiast zagnieżdżony blok shape0 odpowiada osadzonemu obiektowi (sferze). Każdy blok zawiera informacje o całkowitym wkładzie danego obiektu w procesy generowania, rozpadu i utraty.

Tabela 2 zawiera podsumowanie pól, które znajdą Państwo w tym pliku, oraz wyjaśnienia dotyczące sposobu ich interpretacji.

Wszystkie trzy elementy – trójwymiarowe wykresy i zestawienia danych liczbowych zawarte w exciton_sim_info.json – razem tworzą spójny obraz symulacji: wykresy pokazują, gdzie zachodzą poszczególne procesy, natomiast plik JSON określa, w jakim stopniu każdy z tych procesów się do nich przyczynia. Dzięki temu stanowi to użyteczne narzędzie do systematycznej analizy wpływu zmian parametrów lub geometrii na dysocjację ekscytonów i efektywną generację ładunków.

9. Montaż.

Analizując wszystkie wyniki w przestrzeni trójwymiarowej, można dostrzec spójną zależność fizyczną. Absorpcja światła generuje ekscytony głównie wewnątrz obszaru donora. Następnie ekscytony te rozprzestrzeniają się w donorze, ich liczba maleje w pobliżu granicy między donorem a akceptorem, a po dotarciu do otaczającego obszaru akceptora efektywnie ulegają rozszczepieniu. Diagramy przestrzenne pokazują, gdzie zachodzą te procesy, a podsumowujący raport exciton_sim_info.json określa, w jakim stopniu każdy z nich się do nich przyczynia.

W tym przykładzie raport podsumowujący wyraźnie oddziela podstawowe zjawiska fizyczne w zależności od obszaru przestrzennego, co umożliwia przejrzystą analizę procesów generowania, strat i dysocjacji.

• Sfera (donor): Większość absorpcji światła – a zatem i większość generacji ekscytonów – zachodzi wewnątrz sfery donora, co odzwierciedla wysoka wartość tot_G w bloku shape0/sphere. Wewnątrz donora dominującym konkurującym mechanizmem strat w tym przykładzie jest radiacyjna rekombinacja, której wielkość określana jest przez tot_Rk_pl. Stosunek tot_Rk_pl do tot_G stanowi bezpośrednią miarę tego, jak silnie radiacyjne straty tłumią populację ekscytonów, zanim dotrą one do granicy między donorem a akceptorem.
• Warstwa (matryca akceptora): W tym przykładzie otaczająca warstwa słabo absorbuje, dlatego jej wkład w generowanie ekscytonów (tot_G) jest niewielki. Jej główną rolą jest natomiast konwersja ekscytonów: warstwa jest skonfigurowana tak, aby posiadała silny kanał dysocjacji, co sprawia, że tot_Rk_dis jest dominującym składnikiem w bloku layer. To bezpośrednio odpowiada obserwowanej na wykresach przestrzennych dysocjacji na granicy faz, gdzie Gn.csv przedstawia pierścieniową lub muszelkową strefę generowania elektronów, skoncentrowaną na granicy między donorem a akceptorem.

Kluczowy wskaźnik dis_eff przedstawia ogólną efektywność dysocjacji w systemie (w tym przypadku w przybliżeniu \(89\%\)). Fizycznie oznacza to, że większość ekscitonów generowanych w donorze dyfunduje do obszaru styku i ulega dysocjacji, przy czym stosunkowo niewielka ich część ulega utracie w wyniku radiacyjnej rekombinacji. Inne mechanizmy, takie jak transfer Förstera lub anihilacja ekscitonów, są nieaktywne w tym konkretnym przypadku.

Innymi słowy, symulacja odbywa się w warunkach, w których dysocjacja na granicy faz przeważa nad zanikiem ekscytonu, co prowadzi do wysokiej efektywności konwersji fotonów w nośniki ładunku. Połączenie trójwymiarowej wizualizacji i ilościowych danych pozwala uzyskać zwięzły, ale jednocześnie pełny obraz tego, w jaki sposób geometria i parametry kinetyczne wpływają na efektywność generowania ładunku w modelu domen ekscytonowych.