Symulacja ogniwa słonecznego o strukturze heterozłącza AlGaAs/GaAs: struktura urządzenia, dopasowanie pasm energetycznych, charakterystyka prądowo-napięciowa (J–V), wydajność kwantowa (EQE) i zależność napięcia w obwodzie otwartym od natężenia promieniowania słonecznego (Suns–Voc).Symulacja ogniwa słonecznego o strukturze heterozłącza AlGaAs/GaAs: struktura urządzenia, dopasowanie pasm energetycznych, charakterystyka prądowo-napięciowa (J–V), wydajność kwantowa (EQE) i zależność napięcia w obwodzie otwartym od natężenia promieniowania słonecznego (Suns–Voc).
1. Wprowadzenie: dlaczego warto stosować heterozłącze AlGaAs/GaAs?1. Wprowadzenie: dlaczego warto stosować heterozłącze AlGaAs/GaAs?

Arszenek galu (GaAs) to półprzewodnik z grupy III–V o bezpośredniej przerwie energetycznej, powszechnie stosowany w wysokowydajnych urządzeniach fotowoltaicznych. Jego stosunkowo silna absorpcja powyżej krawędzi pasma pozwala na efektywne pochłanianie znacznej części padającego światła słonecznego w odpowiednich zakresach długości fal, nawet przy użyciu cienkiej warstwy absorbującej. Jednak konwencjonalne homozłącze GaAs może powodować niepożądane absorpcje optyczne i rekombinację nośników mniejszościowych w pobliżu oświetlanej powierzchni. Umieszczenie na oświetlanej stronie półprzewodnika o szerszej przerwie energetycznej może zmniejszyć te straty, jeśli odpowiednio dobierze się jego właściwości optyczne i elektroniczne.Arszenek galu (GaAs) to półprzewodnik z grupy III–V o bezpośredniej przerwie energetycznej, powszechnie stosowany w wysokowydajnych urządzeniach fotowoltaicznych. Jego stosunkowo silna absorpcja powyżej krawędzi pasma pozwala na efektywne pochłanianie znacznej części padającego światła słonecznego w odpowiednich zakresach długości fal, nawet przy użyciu cienkiej warstwy absorbującej. Jednak konwencjonalne homozłącze GaAs może powodować niepożądane absorpcje optyczne i rekombinację nośników mniejszościowych w pobliżu oświetlanej powierzchni. Umieszczenie na oświetlanej stronie półprzewodnika o szerszej przerwie energetycznej może zmniejszyć te straty, jeśli odpowiednio dobierze się jego właściwości optyczne i elektroniczne.
W tym samouczku modelowana jest heterospójka typu p – Al0,3Ga0,7As / typu n – GaAs. Dwie górne warstwy AlGaAs tworzą obszar o większej przerwie energetycznej, działający jako okno/emiter typu p, natomiast znajdująca się pod spodem warstwa GaAs typu n stanowi główny obszar absorpcji światła. Różnica w przerwie energetycznej i powinowactwie elektronowym na granicy tych warstw powoduje wystąpienie nieciągłości w pasmach przewodnictwa i walencyjnych, oprócz elektrostatycznego zakrzywienia wywołanego zjonizowanymi domieszkami. W rezultacie, analizując tę heterospójkę, należy uwzględnić zarówno różnice w pasmach energetycznych materiałów, jak i quasi-poziomy Fermiego nośników; traktowanie jej jako diody PN wykonanej wyłącznie z GaAs pomijałoby istotne aspekty fizyczne.W tym samouczku modelowana jest heterospójka typu p – Al0,3Ga0,7As / typu n – GaAs. Dwie górne warstwy AlGaAs tworzą obszar o większej przerwie energetycznej, działający jako okno/emiter typu p, natomiast znajdująca się pod spodem warstwa GaAs typu n stanowi główny obszar absorpcji światła. Różnica w przerwie energetycznej i powinowactwie elektronowym na granicy tych warstw powoduje wystąpienie nieciągłości w pasmach przewodnictwa i walencyjnych, oprócz elektrostatycznego zakrzywienia wywołanego zjonizowanymi domieszkami. W rezultacie, analizując tę heterospójkę, należy uwzględnić zarówno różnice w pasmach energetycznych materiałów, jak i quasi-poziomy Fermiego nośników; traktowanie jej jako diody PN wykonanej wyłącznie z GaAs pomijałoby istotne aspekty fizyczne.
Wykorzystamy jednowymiarowy solver równań dyfuzji i dryfu opracowany przez OghmaNano, połączony z obliczeniami dotyczącymi generacji optycznej. Dostarczony model zawiera dane z pomiarów charakterystyki prądowo-napięciowej (J–V) przy oświetleniu, opcje dotyczące obliczeń zależnych od wydajności kwantowej (EQE) i natężenia oraz zapisane profile wewnętrzne. Proces pracy został celowo oparty na wcześniej skonfigurowanym przykładzie, dlatego celem jest zrozumienie, uruchomienie i analiza działania urządzenia, a nie ponowne wprowadzanie wszystkich parametrów od podstaw.Wykorzystamy jednowymiarowy solver równań dyfuzji i dryfu opracowany przez OghmaNano, połączony z obliczeniami dotyczącymi generacji optycznej. Dostarczony model zawiera dane z pomiarów charakterystyki prądowo-napięciowej (J–V) przy oświetleniu, opcje dotyczące obliczeń zależnych od wydajności kwantowej (EQE) i natężenia oraz zapisane profile wewnętrzne. Proces pracy został celowo oparty na wcześniej skonfigurowanym przykładzie, dlatego celem jest zrozumienie, uruchomienie i analiza działania urządzenia, a nie ponowne wprowadzanie wszystkich parametrów od podstaw.
W urządzeniu fotowoltaicznym pochłonięte fotony generują pary elektron-dziura; selektywne kontakty pobierają nośniki ładunku, a wynikające z tego wewnętrzne rozszczepienie quasi-poziomu Fermiego powoduje powstanie napięcia w obwodzie otwartym. Przydatny model musi łączyć absorpcję światła i generowanie ładunku z elektrostatyką złącza, pobieraniem nośników i rekombinacją. W poniższych sekcjach omówiono bezpośrednio każdy z tych związków.W urządzeniu fotowoltaicznym pochłonięte fotony generują pary elektron-dziura; selektywne kontakty pobierają nośniki ładunku, a wynikające z tego wewnętrzne rozszczepienie quasi-poziomu Fermiego powoduje powstanie napięcia w obwodzie otwartym. Przydatny model musi łączyć absorpcję światła i generowanie ładunku z elektrostatyką złącza, pobieraniem nośników i rekombinacją. W poniższych sekcjach omówiono bezpośrednio każdy z tych związków.
2. Tworzenie symulacji na podstawie demonstracyjnych przykładów wykorzystujących GaAs.2. Tworzenie symulacji na podstawie demonstracyjnych przykładów wykorzystujących GaAs.
Otwórz program OghmaNano i wybierz opcję „Nowa symulacja”. W oknie przeglądarki symulacji otwórz sekcję „Przykłady GaAs”, a następnie wybierz „Heterostrukturalna ogniwo słoneczne AlGaAs/GaAs”. Okno dialogowe wyboru i kategoria przykładów są wyświetlane w ?? i ??. Wybrany przykład otworzy się w głównym oknie struktury urządzenia (??).Otwórz program OghmaNano i wybierz opcję „Nowa symulacja”. W oknie przeglądarki symulacji otwórz sekcję „Przykłady GaAs”, a następnie wybierz „Heterostrukturalna ogniwo słoneczne AlGaAs/GaAs”. Okno dialogowe wyboru i kategoria przykładów są wyświetlane w ?? i ??. Wybrany przykład otworzy się w głównym oknie struktury urządzenia (??).


Główne okno zawiera trójwymiarową, graficzną reprezentację urządzenia. Opisana tutaj symulacja fizyczna wykorzystuje zmienność przestrzenną wzdłuż osi y, prostopadłej do warstw. Zobrazowana warstwa stanowi zatem sposób wizualizacji struktury i połączeń, a nie dowód na pełne, trójwymiarowe rozwiązanie problemu elektrycznego. Aby sprawdzić lub zmodyfikować warstwy, wybierz zakładkę „Struktura urządzenia”; zakładki „Wyniki” i „Złącza” będą używane po uruchomieniu symulacji.Główne okno zawiera trójwymiarową, graficzną reprezentację urządzenia. Opisana tutaj symulacja fizyczna wykorzystuje zmienność przestrzenną wzdłuż osi y, prostopadłej do warstw. Zobrazowana warstwa stanowi zatem sposób wizualizacji struktury i połączeń, a nie dowód na pełne, trójwymiarowe rozwiązanie problemu elektrycznego. Aby sprawdzić lub zmodyfikować warstwy, wybierz zakładkę „Struktura urządzenia”; zakładki „Wyniki” i „Złącza” będą używane po uruchomieniu symulacji.
Przed rozpoczęciem eksperymentów warto zachować kopię zapasową przykładu. Należy zachować niezmienioną wersję początkową, aby można było porównać zmiany dotyczące domieszek, składu stopu lub parametrów interfejsu z działającą wersją referencyjną.Przed rozpoczęciem eksperymentów warto zachować kopię zapasową przykładu. Należy zachować niezmienioną wersję początkową, aby można było porównać zmiany dotyczące domieszek, składu stopu lub parametrów interfejsu z działającą wersją referencyjną.
3. Kontrola struktury warstw epitaksjalnych, domieszkowania i składu stopu.3. Kontrola struktury warstw epitaksjalnych, domieszkowania i składu stopu.
Otwórz edytor warstw, aby zbadać cztery epitaksjalne obszary (??). Poniżej przedstawiono dokładną geometrię początkową i domieszkowanie użyte w tym przykładzie. Wszystkie stężenia wyrażone są jako gęstości objętościowe w jednostkach m−3, zgodnie z danymi zapisanymi w OghmaNano; aby przeliczyć na cm−3, należy podzielić przez 106.Otwórz edytor warstw, aby zbadać cztery epitaksjalne obszary (??). Poniżej przedstawiono dokładną geometrię początkową i domieszkowanie użyte w tym przykładzie. Wszystkie stężenia wyrażone są jako gęstości objętościowe w jednostkach m−3, zgodnie z danymi zapisanymi w OghmaNano; aby przeliczyć na cm−3, należy podzielić przez 106.
| Warstwa, od góry do dołu.Warstwa, od góry do dołu. | GrubośćGrubość | Pozycja od góryPozycja od góry | Dominujący rodzaj domieszkowaniaDominujący rodzaj domieszkowania | Rola wymagająca obecności fizycznej.Rola wymagająca obecności fizycznej. |
|---|---|---|---|---|
| p+ Al0,3Ga0,7Asp+ Al0,3Ga0,7As | 0,10 µm0,10 µm | 0–0,10 µm0–0,10 µm | NA = 2 × 1024 m−3 (2 × 1018 cm−3)NA = 2 × 1024 m−3 (2 × 1018 cm−3) | Warstwa kontaktowa/przezroczysta o wysokim stopniu domieszkowania, znajdująca się na przedniej powierzchni.Warstwa kontaktowa/przezroczysta o wysokim stopniu domieszkowania, znajdująca się na przedniej powierzchni. |
| p-typ Al0,3Ga0,7Asp-typ Al0,3Ga0,7As | 0,30 µm0,30 µm | 0,10–0,40 µm0,10–0,40 µm | NA = 3 × 1023 m−3 (3 × 1017 cm−3)NA = 3 × 1023 m−3 (3 × 1017 cm−3) | emitator typu p o szerokiej przerwie energetycznejemitator typu p o szerokiej przerwie energetycznej |
| absorber z arsenku galu (GaAs)absorber z arsenku galu (GaAs) | 2,80 µm2,80 µm | 0,40–3,20 µm0,40–3,20 µm | ND = 3 × 1022 m−3 (3 × 1016 cm−3)ND = 3 × 1022 m−3 (3 × 1016 cm−3) | Główny obszar pochłaniania fotonów i zbierania ładunków.Główny obszar pochłaniania fotonów i zbierania ładunków. |
| Warstwa n+ z arsenku galu, pełniąca funkcję tylnego kontaktu.Warstwa n+ z arsenku galu, pełniąca funkcję tylnego kontaktu. | 0,20 µm0,20 µm | 3,20–3,40 µm3,20–3,40 µm | ND = 1 × 1024 m−3 (1 × 1018 cm−3)ND = 1 × 1024 m−3 (1 × 1018 cm−3) | Warstwa tylnego kontaktu selektywna dla elektronówWarstwa tylnego kontaktu selektywna dla elektronów |
Heterointerfejs znajduje się w punkcie y = 0,40 µm, gdzie warstwa emitująca typu p, wykonana z AlGaAs, styka się z warstwą absorbującą typu n, wykonaną z GaAs. Jednocześnie następuje zmiana polaryzacji domieszkowania oraz zmiana składu materiału półprzewodnikowego. Natomiast interfejsy o grubości 0,10 µm wewnątrz AlGaAs i 3,20 µm wewnątrz GaAs dotyczą głównie zmian w koncentracji domieszek; granice te nie wprowadzają nowego systemu materiałowego.Heterointerfejs znajduje się w punkcie y = 0,40 µm, gdzie warstwa emitująca typu p, wykonana z AlGaAs, styka się z warstwą absorbującą typu n, wykonaną z GaAs. Jednocześnie następuje zmiana polaryzacji domieszkowania oraz zmiana składu materiału półprzewodnikowego. Natomiast interfejsy o grubości 0,10 µm wewnątrz AlGaAs i 3,20 µm wewnątrz GaAs dotyczą głównie zmian w koncentracji domieszek; granice te nie wprowadzają nowego systemu materiałowego.
Otwarty doping / jony z elementów sterujących elektrycznie (??). Na wykresie powinna być widoczna wysoka koncentracja akceptorów w pobliżu oświetlanego kontaktu, niższa koncentracja akceptorów w pozostałej części emitera AlGaAs, donorowo domieszkowany absorber GaAs oraz silniej domieszkowany tylny kontakt GaAs. Zmiana znaku w wypadkowej koncentracji domieszek wskazuje na obecność złącza PN. Doping sam w sobie nie jest parametrem absorpcji optycznej; wpływa na równanie Poissona poprzez zjonizowany ładunek przestrzenny i wpływa na przewodność, rekombinację oraz zbieranie nośników.Otwarty doping / jony z elementów sterujących elektrycznie (??). Na wykresie powinna być widoczna wysoka koncentracja akceptorów w pobliżu oświetlanego kontaktu, niższa koncentracja akceptorów w pozostałej części emitera AlGaAs, donorowo domieszkowany absorber GaAs oraz silniej domieszkowany tylny kontakt GaAs. Zmiana znaku w wypadkowej koncentracji domieszek wskazuje na obecność złącza PN. Doping sam w sobie nie jest parametrem absorpcji optycznej; wpływa na równanie Poissona poprzez zjonizowany ładunek przestrzenny i wpływa na przewodność, rekombinację oraz zbieranie nośników.



Następnie należy sprawdzić skład stopu (??). Obie warstwy AlGaAs charakteryzują się g_x0 = g_x1 = 0.30: stałą zawartością aluminium na całej ich grubości. Obie warstwy GaAs mają g_x0 = g_x1 = 0 w tym konkretnym przykładzie. W dostarczonej strukturze nie występuje stopniowa zmiana składu. Obszar AlGaAs o liniowo zmieniającym się składzie wymagałby różnych wartości składu na końcach i modelu materiałowego Lua, który byłby w stanie obliczyć lokalny skład w każdym punkcie siatki.Następnie należy sprawdzić skład stopu (??). Obie warstwy AlGaAs charakteryzują się g_x0 = g_x1 = 0.30: stałą zawartością aluminium na całej ich grubości. Obie warstwy GaAs mają g_x0 = g_x1 = 0 w tym konkretnym przykładzie. W dostarczonej strukturze nie występuje stopniowa zmiana składu. Obszar AlGaAs o liniowo zmieniającym się składzie wymagałby różnych wartości składu na końcach i modelu materiałowego Lua, który byłby w stanie obliczyć lokalny skład w każdym punkcie siatki.
Wpływ zastosowania materiału o szerokiej przerwie energetycznej zależy od długości fali. Fotony o energii niższej niż próg absorpcji AlGaAs mogą łatwiej przenikać przez górną warstwę półprzewodnika, natomiast fotony o energii wyższej niż ten próg mogą być również absorbowane w warstwie emisyjnej. Wielkość tego efektu zależy od rzeczywistych stałych optycznych, grubości warstwy i właściwości wewnętrznego zbierania nośników; sama zawartość Al nie zapewnia ilościowego zmniejszenia strat na powierzchni.Wpływ zastosowania materiału o szerokiej przerwie energetycznej zależy od długości fali. Fotony o energii niższej niż próg absorpcji AlGaAs mogą łatwiej przenikać przez górną warstwę półprzewodnika, natomiast fotony o energii wyższej niż ten próg mogą być również absorbowane w warstwie emisyjnej. Wielkość tego efektu zależy od rzeczywistych stałych optycznych, grubości warstwy i właściwości wewnętrznego zbierania nośników; sama zawartość Al nie zapewnia ilościowego zmniejszenia strat na powierzchni.
4. Modele materiałów elektrycznych i przesunięcia poziomów energii w heterospójeniach.4. Modele materiałów elektrycznych i przesunięcia poziomów energii w heterospójeniach.
Otwórz edytor parametrów elektrycznych i sprawdź skrypty materiałów Lua powiązane z obszarami (?? i ??). W tym przykładzie parametry materiału są obliczane na podstawie zarejestrowanych skryptów, a nie poprzez ustawianie każdej wartości parametru oddzielnie na stałą wartość.Otwórz edytor parametrów elektrycznych i sprawdź skrypty materiałów Lua powiązane z obszarami (?? i ??). W tym przykładzie parametry materiału są obliczane na podstawie zarejestrowanych skryptów, a nie poprzez ustawianie każdej wartości parametru oddzielnie na stałą wartość.


Ścieżki Lua dla AlGaAs i GaAs są odrębne: Inorganic/AlGaAs/AlGaAs.lua i Inorganic/GaAs/GaAs.lua. Dla modelu AlGaAs dostępna jest wartość składu warstwy wierzchniej wynosząca 0,30, co umożliwia obliczanie zależnych od składu parametrów elektrycznych. Zestawy danych optycznych są konfigurowane oddzielnie: Inorganic/AlGaAs/Adachi dla dwóch warstw AlGaAs i Inorganic/GaAs/Adachi dla warstw GaAs. Nie należy mylić parametrów elektrycznych pasma przewodnictwa z zestawami danych optycznych dotyczącymi współczynnika załamania i absorpcji.Ścieżki Lua dla AlGaAs i GaAs są odrębne: Inorganic/AlGaAs/AlGaAs.lua i Inorganic/GaAs/GaAs.lua. Dla modelu AlGaAs dostępna jest wartość składu warstwy wierzchniej wynosząca 0,30, co umożliwia obliczanie zależnych od składu parametrów elektrycznych. Zestawy danych optycznych są konfigurowane oddzielnie: Inorganic/AlGaAs/Adachi dla dwóch warstw AlGaAs i Inorganic/GaAs/Adachi dla warstw GaAs. Nie należy mylić parametrów elektrycznych pasma przewodnictwa z zestawami danych optycznych dotyczącymi współczynnika załamania i absorpcji.
Na granicy między różnymi materiałami, lokalna energia krawędzi pasma przewodnictwa jest powiązana z powinowactwem elektronowym za pomocą konwencji wyrażonej wzorem \(E_C=-\chi-q\phi\), a energia krawędzi pasma walencyjnego – za pomocą wzoru \(E_V=E_C-E_g\), pod warunkiem, że po obu stronach stosowana jest ta sama energia odniesienia. Powstałe nieciągłości są zwyczajowo zapisywane jako \(\Delta E_C\) i \(\Delta E_V\). Zależą one zarówno od zależnej od składu przerwy energetycznej, jak i od modelu powinowactwa elektronowego. Większa przerwa energetyczna AlGaAs sama w sobie nie determinuje, w jaki sposób całkowita nieciągłość jest rozdzielana między pasmo przewodnictwa i pasmo walencyjne.Na granicy między różnymi materiałami, lokalna energia krawędzi pasma przewodnictwa jest powiązana z powinowactwem elektronowym za pomocą konwencji wyrażonej wzorem \(E_C=-\chi-q\phi\), a energia krawędzi pasma walencyjnego – za pomocą wzoru \(E_V=E_C-E_g\), pod warunkiem, że po obu stronach stosowana jest ta sama energia odniesienia. Powstałe nieciągłości są zwyczajowo zapisywane jako \(\Delta E_C\) i \(\Delta E_V\). Zależą one zarówno od zależnej od składu przerwy energetycznej, jak i od modelu powinowactwa elektronowego. Większa przerwa energetyczna AlGaAs sama w sobie nie determinuje, w jaki sposób całkowita nieciągłość jest rozdzielana między pasmo przewodnictwa i pasmo walencyjne.
Z punktu widzenia fizyki, odpowiednio duży „wzrost” w paśmie przewodnictwa może utrudniać ekstrakcję elektronów, natomiast odpowiednia przerwa energetyczna w paśmie walencyjnym może ograniczać niepożądany ruch dziur w przeciwnym kierunku. Efekty te konkurują z dryfem, dyfuzją, domieszkowaniem, aktywacją termiczną i rekombinacją. Dlatego symulowana struktura pasmowa dostarcza więcej informacji niż samo podanie różnicy szerokości pasma.Z punktu widzenia fizyki, odpowiednio duży „wzrost” w paśmie przewodnictwa może utrudniać ekstrakcję elektronów, natomiast odpowiednia przerwa energetyczna w paśmie walencyjnym może ograniczać niepożądany ruch dziur w przeciwnym kierunku. Efekty te konkurują z dryfem, dyfuzją, domieszkowaniem, aktywacją termiczną i rekombinacją. Dlatego symulowana struktura pasmowa dostarcza więcej informacji niż samo podanie różnicy szerokości pasma.
4.1. Parametry transportu i rekombinacji4.1. Parametry transportu i rekombinacji
OghmaNano rozwiązuje sprzężone równania Poissona oraz równania ciągłości dla elektronów i dziur. W przybliżeniu dryfowo-dyfuzyjnym gęstości prądu elektronów i dziur zależą od ruchliwości nośników, gradientów gęstości nośników oraz pola elektrostatycznego. Przy użyciu standardowej konwencji znaków, można je wyrazić w następujący sposób:OghmaNano rozwiązuje sprzężone równania Poissona oraz równania ciągłości dla elektronów i dziur. W przybliżeniu dryfowo-dyfuzyjnym gęstości prądu elektronów i dziur zależą od ruchliwości nośników, gradientów gęstości nośników oraz pola elektrostatycznego. Przy użyciu standardowej konwencji znaków, można je wyrazić w następujący sposób:
\[J_n=q\mu_n n\mathcal{E}+qD_n\frac{dn}{dy},\qquad J_p=q\mu_p p\mathcal{E}-qD_p\frac{dp}{dy}.\]Pole elektryczne i koncentracja nośników ładunku są obliczane w sposób samo-spójny. Podczas oświetlania, równania ciągłości muszą dodatkowo uwzględniać fotogenerację \(G\) i rekombinację \(R\). Na przykład, przy użyciu standardowej konwencji prądowej w jednym wymiarze, \(dJ_n/dy=q(R-G)\) i \(dJ_p/dy=q(G-R)\) w stanie ustalonym. W konsekwencji, ich suma powinna pozostawać stała w przestrzeni w idealnym rozwiązaniu stanu ustalonego, bez dodatkowego rozproszonego członu związanego z magazynowaniem ładunku.Pole elektryczne i koncentracja nośników ładunku są obliczane w sposób samo-spójny. Podczas oświetlania, równania ciągłości muszą dodatkowo uwzględniać fotogenerację \(G\) i rekombinację \(R\). Na przykład, przy użyciu standardowej konwencji prądowej w jednym wymiarze, \(dJ_n/dy=q(R-G)\) i \(dJ_p/dy=q(G-R)\) w stanie ustalonym. W konsekwencji, ich suma powinna pozostawać stała w przestrzeni w idealnym rozwiązaniu stanu ustalonego, bez dodatkowego rozproszonego członu związanego z magazynowaniem ładunku.
Format JSON umożliwia modelowanie rekombinacji SRH w stanie ustalonym (ss_srh_enabled = "True") i określa przybliżoną energię pułapki w środkowej części pasma energetycznego, gęstość pułapek \(N_t=5\times10^{19}\,\mathrm{m}^{-3}\) oraz efektywny przekrój poprzeczny wychwytywania \(\sigma_n=\sigma_p=10^{-22}\,\mathrm{m}^2\) w każdej warstwie. Te pola opisują skonfigurowany model defektów; należy sprawdzić, czy dane pole jest pobierane bezpośrednio z pliku JSON, czy też jest przeliczane przez wybrany skrypt materiałowy, i porównać to z wartościami efektywnymi raportowanymi przez OghmaNano. Szybkość standardowego procesu SRH na jednym poziomie wynosi:Format JSON umożliwia modelowanie rekombinacji SRH w stanie ustalonym (ss_srh_enabled = "True") i określa przybliżoną energię pułapki w środkowej części pasma energetycznego, gęstość pułapek \(N_t=5\times10^{19}\,\mathrm{m}^{-3}\) oraz efektywny przekrój poprzeczny wychwytywania \(\sigma_n=\sigma_p=10^{-22}\,\mathrm{m}^2\) w każdej warstwie. Te pola opisują skonfigurowany model defektów; należy sprawdzić, czy dane pole jest pobierane bezpośrednio z pliku JSON, czy też jest przeliczane przez wybrany skrypt materiałowy, i porównać to z wartościami efektywnymi raportowanymi przez OghmaNano. Szybkość standardowego procesu SRH na jednym poziomie wynosi:
W półprzewodniku o bezpośredniej przerwie energetycznej również ważna może być radiacyjna rekombinacja swobodnych elektronów ze swobodnymi dziurami, którą zazwyczaj przedstawia się jako \(R_{\mathrm{rad}}=B(np-n_i^2)\). Program może również dostarczać współczynniki Augera i inne parametry za pośrednictwem modeli materiałowych. Nie należy wnioskować o dominującym mechanizmie jedynie na podstawie obecności jego ustawień: należy wykorzystać wygenerowane pliki rekombinacji, rzeczywiste włączone przełączniki oraz wartości modelu materiałowego w działającym przykładzie. W przypadku tego urządzenia odpowiednią interpretację uzyskuje się poprzez porównanie wartości dla tego samego punktu napięcia i przy użyciu spójnej skali pionowej.W półprzewodniku o bezpośredniej przerwie energetycznej również ważna może być radiacyjna rekombinacja swobodnych elektronów ze swobodnymi dziurami, którą zazwyczaj przedstawia się jako \(R_{\mathrm{rad}}=B(np-n_i^2)\). Program może również dostarczać współczynniki Augera i inne parametry za pośrednictwem modeli materiałowych. Nie należy wnioskować o dominującym mechanizmie jedynie na podstawie obecności jego ustawień: należy wykorzystać wygenerowane pliki rekombinacji, rzeczywiste włączone przełączniki oraz wartości modelu materiałowego w działającym przykładzie. W przypadku tego urządzenia odpowiednią interpretację uzyskuje się poprzez porównanie wartości dla tego samego punktu napięcia i przy użyciu spójnej skali pionowej.
Widok edytora elektrycznego można wykorzystać do analizy właściwości związanych z mobilnością i innymi właściwościami lokalnymi (??). Podczas eksperymentowania z różnymi składami, należy wyłączyć możliwość ręcznej zmiany wartości w bazie danych dotyczących szerokości pasma i powinowactwa elektronowego, chyba że istnieje konkretny powód, aby to zrobić: ręczne ustawienie wartości \(E_g\) lub \(\chi\) może przypadkowo usunąć zamierzoną zależność od składu.Widok edytora elektrycznego można wykorzystać do analizy właściwości związanych z mobilnością i innymi właściwościami lokalnymi (??). Podczas eksperymentowania z różnymi składami, należy wyłączyć możliwość ręcznej zmiany wartości w bazie danych dotyczących szerokości pasma i powinowactwa elektronowego, chyba że istnieje konkretny powód, aby to zrobić: ręczne ustawienie wartości \(E_g\) lub \(\chi\) może przypadkowo usunąć zamierzoną zależność od składu.

5. Elektryczne styki i siatka elementów użyta w symulacji.5. Elektryczne styki i siatka elementów użyta w symulacji.
Model posiada dwa zewnętrzne zaciski elektryczne. Górna anoda typu p, oznaczona jako y0, jest zaciskiem aktywnym, natomiast dolna katoda typu n, oznaczona jako y1, jest zaciskiem uziemionym. Oba zaciski wykorzystują granicę o wysokiej przewodności dla nośników większościowych oraz granicę blokującą dla nośników mniejszościowych. Jest to istotny wybór w procesie modelowania: struktura nie jest reprezentowana przez ogólne, symetryczne zaciski, które losowo pobierają elektrony i dziury. Otwórz edytor zacisków, aby zobaczyć rzędy zacisków (??).Model posiada dwa zewnętrzne zaciski elektryczne. Górna anoda typu p, oznaczona jako y0, jest zaciskiem aktywnym, natomiast dolna katoda typu n, oznaczona jako y1, jest zaciskiem uziemionym. Oba zaciski wykorzystują granicę o wysokiej przewodności dla nośników większościowych oraz granicę blokującą dla nośników mniejszościowych. Jest to istotny wybór w procesie modelowania: struktura nie jest reprezentowana przez ogólne, symetryczne zaciski, które losowo pobierają elektrony i dziury. Otwórz edytor zacisków, aby zobaczyć rzędy zacisków (??).
Na anodzie dominującym nośnikiem ładunku jest dziura, a nośnikiem mniejszościowym jest elektron. Na katodzie role się odwracają. Zastosowane napięcie zmienia więc potencjał na zacisku p w odniesieniu do uziemionego zacisku n. Model kontaktu różni się od silnie domieszkowanych warstw epitaksjalnych, które stanowią część geometrii półprzewodnika. Określone warunki brzegowe definiują, w jaki sposób nośniki opuszczają obszar obliczeniowy.Na anodzie dominującym nośnikiem ładunku jest dziura, a nośnikiem mniejszościowym jest elektron. Na katodzie role się odwracają. Zastosowane napięcie zmienia więc potencjał na zacisku p w odniesieniu do uziemionego zacisku n. Model kontaktu różni się od silnie domieszkowanych warstw epitaksjalnych, które stanowią część geometrii półprzewodnika. Określone warunki brzegowe definiują, w jaki sposób nośniki opuszczają obszar obliczeniowy.
Moduł obliczeniowy dla pól elektrycznych jest skonfigurowany do pracy z jednowymiarową siatką wzdłuż osi y. Jego cztery segmenty siatki śledzą cztery wartości grubości warstw epitaksjalnych, przy czym liczba punktów w każdym segmencie wynosi odpowiednio 10, 17, 156 i 11. Oznacza to łącznie około 194 określone punkty, chociaż rzeczywista liczba unikalnych węzłów, które są wykorzystywane, może zależeć od sposobu, w jaki moduł obliczeniowy obsługuje wspólne granice. Oba kierunki siatki elektrycznej są wyłączone. Ustawienia siatki są widoczne w ??.Moduł obliczeniowy dla pól elektrycznych jest skonfigurowany do pracy z jednowymiarową siatką wzdłuż osi y. Jego cztery segmenty siatki śledzą cztery wartości grubości warstw epitaksjalnych, przy czym liczba punktów w każdym segmencie wynosi odpowiednio 10, 17, 156 i 11. Oznacza to łącznie około 194 określone punkty, chociaż rzeczywista liczba unikalnych węzłów, które są wykorzystywane, może zależeć od sposobu, w jaki moduł obliczeniowy obsługuje wspólne granice. Oba kierunki siatki elektrycznej są wyłączone. Ustawienia siatki są widoczne w ??.


Gęstość siatki ma szczególne znaczenie, zwłaszcza w pobliżu heterointerfejsu, przy wartości 0,40 µm, gdzie elektrostatyczne odkształcenia i przesunięcia pasm mogą zmieniać się na krótkim odcinku. Jeśli obliczone przerwanie pasma lub gęstość nośników tworzą ostry pik, sprawdzenie zbieżności siatki jest bardziej przekonujące niż po prostu przyjęcie domyślnych ustawień. Powtórz obliczenia, używając bardziej szczegółowej siatki dla emitera/złącza, i porównaj wyniki J–V, a także odpowiednie profile przestrzenne. Jeśli te wyniki ulegają znaczącym zmianom, oznacza to, że podstawowa siatka nie była jeszcze wystarczająco gęsta.Gęstość siatki ma szczególne znaczenie, zwłaszcza w pobliżu heterointerfejsu, przy wartości 0,40 µm, gdzie elektrostatyczne odkształcenia i przesunięcia pasm mogą zmieniać się na krótkim odcinku. Jeśli obliczone przerwanie pasma lub gęstość nośników tworzą ostry pik, sprawdzenie zbieżności siatki jest bardziej przekonujące niż po prostu przyjęcie domyślnych ustawień. Powtórz obliczenia, używając bardziej szczegółowej siatki dla emitera/złącza, i porównaj wyniki J–V, a także odpowiednie profile przestrzenne. Jeśli te wyniki ulegają znaczącym zmianom, oznacza to, że podstawowa siatka nie była jeszcze wystarczająco gęsta.
6. Uruchamianie symulacji J–V z podświetleniem.6. Uruchamianie symulacji J–V z podświetleniem.
Domyślny, wybrany tryb w dostarczonej konfiguracji to J–V. Zakres zmian napięcia jest ustawiony od –0,10 V do +1,00 V, ze skokiem 0,02 V. Do oświetlenia wykorzystywany jest widm AM1.5G, a skala intensywności źródła światła bazowego odpowiada natężeniu promieniowania słonecznego (Psun = 1). Generowanie sygnału optycznego jest sprzężone z symulacją elektryczną za pomocą skonfigurowanego modelu generowania macierzy przenoszenia. Ustawiona temperatura wynosi 300 K.Domyślny, wybrany tryb w dostarczonej konfiguracji to J–V. Zakres zmian napięcia jest ustawiony od –0,10 V do +1,00 V, ze skokiem 0,02 V. Do oświetlenia wykorzystywany jest widm AM1.5G, a skala intensywności źródła światła bazowego odpowiada natężeniu promieniowania słonecznego (Psun = 1). Generowanie sygnału optycznego jest sprzężone z symulacją elektryczną za pomocą skonfigurowanego modelu generowania macierzy przenoszenia. Ustawiona temperatura wynosi 300 K.
Przed uruchomieniem sprawdź zakładkę „Wyniki” (??). Kliknij „Uruchom symulację”. W zakładce „Terminal” wyświetlane są informacje o postępie i wyniki obliczeń; na zrzutach ekranu z terminala (?? i ??) można zobaczyć wczesny i późniejszy etap obliczeń. Zmieniający się wyświetlacz w terminalu świadczy o tym, że program działa, ale sam w sobie nie jest dowodem na to, że wszystkie punkty zostały obliczone. Sprawdź, czy obliczenia zostały zakończone, czy wystąpiły błędy oraz czy wygenerowano pliki wynikowe.Przed uruchomieniem sprawdź zakładkę „Wyniki” (??). Kliknij „Uruchom symulację”. W zakładce „Terminal” wyświetlane są informacje o postępie i wyniki obliczeń; na zrzutach ekranu z terminala (?? i ??) można zobaczyć wczesny i późniejszy etap obliczeń. Zmieniający się wyświetlacz w terminalu świadczy o tym, że program działa, ale sam w sobie nie jest dowodem na to, że wszystkie punkty zostały obliczone. Sprawdź, czy obliczenia zostały zakończone, czy wystąpiły błędy oraz czy wygenerowano pliki wynikowe.




Po zakończeniu symulacji wróć do sekcji „Wyniki” i otwórz jv.csv (??). Oświetlona krzywa J–V różni się koncepcyjnie od krzywej diody w stanie przewodzenia w ciemności. Proste wyrażenie pierwszego rzędu to \(J(V)=J_{\mathrm{dark}}(V)-J_{\mathrm{ph}}(V)\), przy czym przyjęto konwencję znaków, zgodnie z którą wyekstrahowany prąd fotowoltaiczny jest ujemny. Wykres przedstawiony przez symulator odzwierciedla zdefiniowaną przez użytkownika konwencję znaków dla danego kontaktu, dlatego fizycznie istotne wartości należy identyfikować na podstawie punktów przecięcia krzywej z osiami napięcia i prądu, zamiast zakładać, że każdy wykres fotowoltaiczny jest rysowany przy użyciu tych samych znaków.Po zakończeniu symulacji wróć do sekcji „Wyniki” i otwórz jv.csv (??). Oświetlona krzywa J–V różni się koncepcyjnie od krzywej diody w stanie przewodzenia w ciemności. Proste wyrażenie pierwszego rzędu to \(J(V)=J_{\mathrm{dark}}(V)-J_{\mathrm{ph}}(V)\), przy czym przyjęto konwencję znaków, zgodnie z którą wyekstrahowany prąd fotowoltaiczny jest ujemny. Wykres przedstawiony przez symulator odzwierciedla zdefiniowaną przez użytkownika konwencję znaków dla danego kontaktu, dlatego fizycznie istotne wartości należy identyfikować na podstawie punktów przecięcia krzywej z osiami napięcia i prądu, zamiast zakładać, że każdy wykres fotowoltaiczny jest rysowany przy użyciu tych samych znaków.
Gęstość prądu zwarciowego, \(J_{SC}\), to gęstość prądu w punkcie \(V=0\), natomiast napięcie w obwodzie otwartym, \(V_{OC}\), to napięcie, przy którym gęstość prądu na zaciskach wynosi zero. Maksymalna moc uzyskiwana jest przy określonym napięciu, które maksymalizuje wartość mocy dostarczanej do zewnętrznego obciążenia, a współczynnik wypełnienia jest definiowany przez:Gęstość prądu zwarciowego, \(J_{SC}\), to gęstość prądu w punkcie \(V=0\), natomiast napięcie w obwodzie otwartym, \(V_{OC}\), to napięcie, przy którym gęstość prądu na zaciskach wynosi zero. Maksymalna moc uzyskiwana jest przy określonym napięciu, które maksymalizuje wartość mocy dostarczanej do zewnętrznego obciążenia, a współczynnik wypełnienia jest definiowany przez:
\[FF=\frac{|J_{\mathrm{mpp}}V_{\mathrm{mpp}}|}{|J_{SC}|V_{OC}}.\]Przy stałej, określonej wartości gęstości mocy optycznej padającej na powierzchnię \(P_{\mathrm{in}}\), sprawność konwersji zmienia się zgodnie z \(\eta=|J_{\mathrm{mpp}}V_{\mathrm{mpp}}|/P_{\mathrm{in}}\). Nie należy podawać wartości sprawności jedynie na podstawie kształtu krzywej; konieczne jest również uwzględnienie konwencji dotyczącej mocy oświetlenia oraz wiarygodne określenie punktu pracy o maksymalnej mocy.Przy stałej, określonej wartości gęstości mocy optycznej padającej na powierzchnię \(P_{\mathrm{in}}\), sprawność konwersji zmienia się zgodnie z \(\eta=|J_{\mathrm{mpp}}V_{\mathrm{mpp}}|/P_{\mathrm{in}}\). Nie należy podawać wartości sprawności jedynie na podstawie kształtu krzywej; konieczne jest również uwzględnienie konwencji dotyczącej mocy oświetlenia oraz wiarygodne określenie punktu pracy o maksymalnej mocy.



Krzywa J–V stanowi najprostsze podsumowanie parametrów, ale nie wyjaśnia, dlaczego dochodzi do strat prądu. Spadek prądu fotokatalitycznego może być spowodowany niewystarczającą absorpcją, słabym zbieraniem nośników ładunku lub zwiększoną rekombinacją; niższe napięcie w obwodzie otwartym może wynikać z silniejszej rekombinacji lub zmian w elektrostatyce złącza. Z tego powodu w kolejnych sekcjach wykorzystywane są dane przestrzenne jako materiał dowodowy, a nie tylko wykres J–V.Krzywa J–V stanowi najprostsze podsumowanie parametrów, ale nie wyjaśnia, dlaczego dochodzi do strat prądu. Spadek prądu fotokatalitycznego może być spowodowany niewystarczającą absorpcją, słabym zbieraniem nośników ładunku lub zwiększoną rekombinacją; niższe napięcie w obwodzie otwartym może wynikać z silniejszej rekombinacji lub zmian w elektrostatyce złącza. Z tego powodu w kolejnych sekcjach wykorzystywane są dane przestrzenne jako materiał dowodowy, a nie tylko wykres J–V.
device.csv Wizualizacja (??) stanowi przydatne, niezależne narzędzie do sprawdzenia, czy symulowane obszary fizyczne i ich granice odpowiadają zamierzonej geometrii. Panel z informacjami o symulacji (??) jest przydatny podczas porównywania różnych wariantów symulacji i sprawdzania ustawień związanych z danym wynikiem.device.csv Wizualizacja (??) stanowi przydatne, niezależne narzędzie do sprawdzenia, czy symulowane obszary fizyczne i ich granice odpowiadają zamierzonej geometrii. Panel z informacjami o symulacji (??) jest przydatny podczas porównywania różnych wariantów symulacji i sprawdzania ustawień związanych z danym wynikiem.
7. Schematy pasm energetycznych, gęstość nośników ładunku i transport w strukturze heterospójenia.7. Schematy pasm energetycznych, gęstość nośników ładunku i transport w strukturze heterospójenia.
Otwórz zapisane migawki w oknie wyświetlającym wyniki i przesuń suwak migawki przez punkty odchylenia. Cztery zrzuty ekranu dla każdej zmiennej przedstawione poniżej pokazują różne zapisane stany z tego samego obliczenia; ich skale numeryczne i etykiety odchylenia należy odczytywać bezpośrednio z odpowiedniego okna wyświetlania. Porównując dwa stany, najpierw sprawdź, czy mają one takie same zakresy osi i która zmienna jest porównywana.Otwórz zapisane migawki w oknie wyświetlającym wyniki i przesuń suwak migawki przez punkty odchylenia. Cztery zrzuty ekranu dla każdej zmiennej przedstawione poniżej pokazują różne zapisane stany z tego samego obliczenia; ich skale numeryczne i etykiety odchylenia należy odczytywać bezpośrednio z odpowiedniego okna wyświetlania. Porównując dwa stany, najpierw sprawdź, czy mają one takie same zakresy osi i która zmienna jest porównywana.
7.1. Pasma przewodnictwa i walencyjne, poziom Fermiego7.1. Pasma przewodnictwa i walencyjne, poziom Fermiego
Schemat pasm energetycznych jest najłatwiejszym sposobem na zidentyfikowanie heterospójenia AlGaAs/GaAs (??–??). Granica między materiałami znajduje się w odległości 0,40 µm. Nagła zmiana lokalnych krawędzi pasma energetycznego w tym miejscu odzwierciedla przesunięcie materiału, natomiast płynna zmiana krzywizny w innych miejscach odzwierciedla właściwości elektrostatyczne. To rozróżnienie jest kluczowe: heterospójenia nie można interpretować wyłącznie w kategoriach ciągłego zakrzywienia pasma energetycznego GaAs.Schemat pasm energetycznych jest najłatwiejszym sposobem na zidentyfikowanie heterospójenia AlGaAs/GaAs (??–??). Granica między materiałami znajduje się w odległości 0,40 µm. Nagła zmiana lokalnych krawędzi pasma energetycznego w tym miejscu odzwierciedla przesunięcie materiału, natomiast płynna zmiana krzywizny w innych miejscach odzwierciedla właściwości elektrostatyczne. To rozróżnienie jest kluczowe: heterospójenia nie można interpretować wyłącznie w kategoriach ciągłego zakrzywienia pasma energetycznego GaAs.
W stanie równowagi termicznej potencjały elektrochemiczne elektronów i dziur są identyczne, co prowadzi do powstania wspólnego poziomu Fermiego. W stanie innym niż równowaga, poziomy quasi-Fermiego elektronów i dziur mogą się rozdzielać. Różnica między nimi jest miarą lokalnej nierównowagi nośników i jest związana z lokalnym iloczynem koncentracji elektronów i dziur. W obwodzie otwartym, napięcie na zaciskach odzwierciedla różnicę między potencjałami elektrochemicznymi ustalonymi na selektywnych kontaktach; nie należy utożsamiać pojedynczej lokalnej różnicy poziomów quasi-Fermiego w dowolnej głębokości z mierzonym napięciem w obwodzie otwartym.W stanie równowagi termicznej potencjały elektrochemiczne elektronów i dziur są identyczne, co prowadzi do powstania wspólnego poziomu Fermiego. W stanie innym niż równowaga, poziomy quasi-Fermiego elektronów i dziur mogą się rozdzielać. Różnica między nimi jest miarą lokalnej nierównowagi nośników i jest związana z lokalnym iloczynem koncentracji elektronów i dziur. W obwodzie otwartym, napięcie na zaciskach odzwierciedla różnicę między potencjałami elektrochemicznymi ustalonymi na selektywnych kontaktach; nie należy utożsamiać pojedynczej lokalnej różnicy poziomów quasi-Fermiego w dowolnej głębokości z mierzonym napięciem w obwodzie otwartym.




Wraz ze zmianą przyłożonego napięcia, potencjał elektrostatyczny i rozkład nośników ładunku zmieniają się jednocześnie. W warunkach przewodzenia, bariera potencjału w obszarze złącza zmniejsza się w porównaniu ze stanem równowagi; podczas pracy w trybie fotowoltaicznym, fotonośniki tworzą populacje w stanie innym niż równowaga, co prowadzi do powstania fotonapięcia. Każde widoczne skokowe zmiany w wykresie quasi-poziomu Fermiego należy analizować w odniesieniu do przepływu prądu i rozdzielczości siatki, zanim zostaną one zinterpretowane jako nagła zmiana potencjału elektrochemicznego.Wraz ze zmianą przyłożonego napięcia, potencjał elektrostatyczny i rozkład nośników ładunku zmieniają się jednocześnie. W warunkach przewodzenia, bariera potencjału w obszarze złącza zmniejsza się w porównaniu ze stanem równowagi; podczas pracy w trybie fotowoltaicznym, fotonośniki tworzą populacje w stanie innym niż równowaga, co prowadzi do powstania fotonapięcia. Każde widoczne skokowe zmiany w wykresie quasi-poziomu Fermiego należy analizować w odniesieniu do przepływu prądu i rozdzielczości siatki, zanim zostaną one zinterpretowane jako nagła zmiana potencjału elektrochemicznego.
7.2. Koncentracja elektronów i dziur.7.2. Koncentracja elektronów i dziur.
Wykresy gęstości nośników (??–??) pozwalają odróżnić obszary, w których przeważają nośniki większościowe, od obszarów, w których znajdują się nośniki mniejszościowe, wprowadzone lub generowane przez światło. W obszarze emitera typu p nośnikami większościowymi są dziury, natomiast w obszarze absorbora typu n nośnikami większościowymi są elektrony. Profile rozkładu nośników mniejszościowych zawierają ważne informacje na temat procesów zbierania i rekombinacji, szczególnie w obszarze GaAs złącza.Wykresy gęstości nośników (??–??) pozwalają odróżnić obszary, w których przeważają nośniki większościowe, od obszarów, w których znajdują się nośniki mniejszościowe, wprowadzone lub generowane przez światło. W obszarze emitera typu p nośnikami większościowymi są dziury, natomiast w obszarze absorbora typu n nośnikami większościowymi są elektrony. Profile rozkładu nośników mniejszościowych zawierają ważne informacje na temat procesów zbierania i rekombinacji, szczególnie w obszarze GaAs złącza.
Stężenie elektronów i dziur jest wrażliwe na efektywną gęstość stanów, przesunięcia pasm energetycznych, potencjał elektrostatyczny oraz równowagę procesów generacji i rekombinacji. Zmiany stężenia nie mogą zatem być przypisywane wyłącznie profilowi domieszkowania. Gwałtowna zmiana w pobliżu heterointerfejsu może wynikać z połączonego wpływu lokalnej konfiguracji pasm energetycznych oraz transportu przez dryf i dyfuzję.Stężenie elektronów i dziur jest wrażliwe na efektywną gęstość stanów, przesunięcia pasm energetycznych, potencjał elektrostatyczny oraz równowagę procesów generacji i rekombinacji. Zmiany stężenia nie mogą zatem być przypisywane wyłącznie profilowi domieszkowania. Gwałtowna zmiana w pobliżu heterointerfejsu może wynikać z połączonego wpływu lokalnej konfiguracji pasm energetycznych oraz transportu przez dryf i dyfuzję.




7.3. Gęstości prądu elektronów i dziur.7.3. Gęstości prądu elektronów i dziur.
Wykresy prądów elektronów i dziur (??–??) pokazują, w jaki sposób całkowity prąd na zaciskach jest rozdzielany między dwa rodzaje nośników ładunku w zależności od głębokości. Fotogeneracja powoduje powstawanie zarówno elektronów, jak i dziur, ale sposób, w jaki każdy z tych rodzajów nośników przyczynia się do prądu, zmienia się w całym urządzeniu w zależności od lokalnej generacji, rekombinacji, ruchliwości i pola elektrycznego.Wykresy prądów elektronów i dziur (??–??) pokazują, w jaki sposób całkowity prąd na zaciskach jest rozdzielany między dwa rodzaje nośników ładunku w zależności od głębokości. Fotogeneracja powoduje powstawanie zarówno elektronów, jak i dziur, ale sposób, w jaki każdy z tych rodzajów nośników przyczynia się do prądu, zmienia się w całym urządzeniu w zależności od lokalnej generacji, rekombinacji, ruchliwości i pola elektrycznego.
W statycznym, jednowymiarowym modelu obliczeniowym suma wartościowych prądów elektronów i dziur powinna być stała w każdym punkcie, z wyjątkiem miejsc, w których występuje wyraźnie zdefiniowane, rozproszone źródło prądu lub miejsce strat. Jednak poszczególne składowe nie muszą być stałe: w miarę generowania i rekombinacji nośników ładunku, składowe prądu elektronów i dziur wymieniają się udziałami, zachowując przy tym całkowitą wartość prądu. Jest to skuteczny sposób weryfikacji spójności uzyskanych wyników.W statycznym, jednowymiarowym modelu obliczeniowym suma wartościowych prądów elektronów i dziur powinna być stała w każdym punkcie, z wyjątkiem miejsc, w których występuje wyraźnie zdefiniowane, rozproszone źródło prądu lub miejsce strat. Jednak poszczególne składowe nie muszą być stałe: w miarę generowania i rekombinacji nośników ładunku, składowe prądu elektronów i dziur wymieniają się udziałami, zachowując przy tym całkowitą wartość prądu. Jest to skuteczny sposób weryfikacji spójności uzyskanych wyników.




Niektóre wykresy mogą wykazywać wyraźne zmiany w pobliżu niskich wartości prądu. Należy unikać diagnozowania fizycznej bariery transportowej na podstawie pojedynczego, zakłóconego sygnału: należy jednocześnie analizować krawędzie pasm, gęstość nośników i składowe prądu, a także sprawdzić, czy dany element zmienia się wraz z udoskonaleniem siatki. Oddzielny wykres mobilności (??) stanowi wygodny sposób na obserwację modelowanych lokalnych współczynników transportu i sprawdzenie, czy zmieniają się one w zależności od materiału lub składu.Niektóre wykresy mogą wykazywać wyraźne zmiany w pobliżu niskich wartości prądu. Należy unikać diagnozowania fizycznej bariery transportowej na podstawie pojedynczego, zakłóconego sygnału: należy jednocześnie analizować krawędzie pasm, gęstość nośników i składowe prądu, a także sprawdzić, czy dany element zmienia się wraz z udoskonaleniem siatki. Oddzielny wykres mobilności (??) stanowi wygodny sposób na obserwację modelowanych lokalnych współczynników transportu i sprawdzenie, czy zmieniają się one w zależności od materiału lub składu.

8. Rekombinacja, szybkość przenoszenia ładunku i fotogeneracja.8. Rekombinacja, szybkość przenoszenia ładunku i fotogeneracja.
Prąd generowany przez ogniwo słoneczne jest wynikiem rywalizacji między procesem generowania ładunków a procesem ich utraty. Fotony pochłaniane przez urządzenie tworzą nośniki ładunku, ale tylko niewielka ich część dociera do elektrod. Lokalna równowaga między elektronami i dziurami oddaje tę zależność. OghmaNano umożliwia uzyskanie różnych wartości prądu wyjściowego, które można porównywać przestrzennie, a nie tylko na podstawie prądu mieranego wewnątrz urządzenia.Prąd generowany przez ogniwo słoneczne jest wynikiem rywalizacji między procesem generowania ładunków a procesem ich utraty. Fotony pochłaniane przez urządzenie tworzą nośniki ładunku, ale tylko niewielka ich część dociera do elektrod. Lokalna równowaga między elektronami i dziurami oddaje tę zależność. OghmaNano umożliwia uzyskanie różnych wartości prądu wyjściowego, które można porównywać przestrzennie, a nie tylko na podstawie prądu mieranego wewnątrz urządzenia.
8.1. Rekombinacja a zastosowane napięcie polaryzacyjne.8.1. Rekombinacja a zastosowane napięcie polaryzacyjne.
Otwórz profile rekombinacji (??–??). Współczynnik rekombinacji to wielkość określająca szybkość procesu, zwykle wyrażana w liczbie nośników na metr sześcienny na sekundę; nie należy go porównywać numerycznie z gęstością prądu bez uwzględnienia integracji wzdłuż odległości i pomnożenia przez elementarny ładunek. W przypadku urządzenia jednowymiarowego, przykładowa gęstość prądu rekombinacyjnego wynosi \(J_{\mathrm{rec}}=q\int R(y)\,dy\) i zależy od wybranej definicji wypadkowej szybkości oraz granic.Otwórz profile rekombinacji (??–??). Współczynnik rekombinacji to wielkość określająca szybkość procesu, zwykle wyrażana w liczbie nośników na metr sześcienny na sekundę; nie należy go porównywać numerycznie z gęstością prądu bez uwzględnienia integracji wzdłuż odległości i pomnożenia przez elementarny ładunek. W przypadku urządzenia jednowymiarowego, przykładowa gęstość prądu rekombinacyjnego wynosi \(J_{\mathrm{rec}}=q\int R(y)\,dy\) i zależy od wybranej definicji wypadkowej szybkości oraz granic.
Porównaj nie tylko maksymalne wartości, ale także szerokość i położenie tych wartości. Wąski, wysoki pik na wykresie może w mniejszym stopniu przyczyniać się do całkowitej straty niż umiarkowana wartość rozciągająca się na całej powierzchni absorbentu. Podczas porównywania mechanizmów należy stosować to samo napięcie i obszar. W przypadku struktury heterospójeniowej, akumulacja nośników w pobliżu przesunięcia pasma może lokalnie zwiększać rekombinację, ale kluczowe jest ustalenie, czy taka akumulacja jest rzeczywiście widoczna na wykresach gęstości i pasm.Porównaj nie tylko maksymalne wartości, ale także szerokość i położenie tych wartości. Wąski, wysoki pik na wykresie może w mniejszym stopniu przyczyniać się do całkowitej straty niż umiarkowana wartość rozciągająca się na całej powierzchni absorbentu. Podczas porównywania mechanizmów należy stosować to samo napięcie i obszar. W przypadku struktury heterospójeniowej, akumulacja nośników w pobliżu przesunięcia pasma może lokalnie zwiększać rekombinację, ale kluczowe jest ustalenie, czy taka akumulacja jest rzeczywiście widoczna na wykresach gęstości i pasm.




8.2. Częstotliwość rekombinacji elektronów i dziur.8.2. Częstotliwość rekombinacji elektronów i dziur.
Powiązane wartości szybkości rekombinacji elektronów i dziur (??–??) ułatwiają porównanie szybkości, które są wyświetlane w oknie podglądu. Ważne jest, aby odróżnić całkowitą rekombinację od poszczególnych procesów wychwytu, emisji, generacji lub przejścia nośników swobodnych. Nie każda przedstawiona na wykresie wartość szybkości jest koniecznie niezależnym mechanizmem strat; sumowanie wyświetlanych krzywych bez sprawdzania ich definicji może prowadzić do podwójnego liczenia ścieżek transferu.Powiązane wartości szybkości rekombinacji elektronów i dziur (??–??) ułatwiają porównanie szybkości, które są wyświetlane w oknie podglądu. Ważne jest, aby odróżnić całkowitą rekombinację od poszczególnych procesów wychwytu, emisji, generacji lub przejścia nośników swobodnych. Nie każda przedstawiona na wykresie wartość szybkości jest koniecznie niezależnym mechanizmem strat; sumowanie wyświetlanych krzywych bez sprawdzania ich definicji może prowadzić do podwójnego liczenia ścieżek transferu.
W przypadku kanału rekombinacji radiacyjnej, znane prawo \(R_{\mathrm{rad}}=B(np-n_i^2)\) ilustruje, dlaczego rekombinacja nasila się wraz ze wzrostem koncentracji nośników. Rekombinacja SRH zależy natomiast od zajętości pułapek i kinetyki wychwytywania. Względna waga tych procesów zależy od efektywnych parametrów materiału w danym procesie. Niniejszy poradnik nie zakłada, że rekombinacja radiacyjna dominuje we wszystkich warstwach lub punktach pracy tylko dlatego, że GaAs jest półprzewodnikiem o bezpośredniej przerwie energetycznej.W przypadku kanału rekombinacji radiacyjnej, znane prawo \(R_{\mathrm{rad}}=B(np-n_i^2)\) ilustruje, dlaczego rekombinacja nasila się wraz ze wzrostem koncentracji nośników. Rekombinacja SRH zależy natomiast od zajętości pułapek i kinetyki wychwytywania. Względna waga tych procesów zależy od efektywnych parametrów materiału w danym procesie. Niniejszy poradnik nie zakłada, że rekombinacja radiacyjna dominuje we wszystkich warstwach lub punktach pracy tylko dlatego, że GaAs jest półprzewodnikiem o bezpośredniej przerwie energetycznej.




8.3. Profil generacji8.3. Profil generacji
Współczynnik generowania nośników optycznych zasadniczo różni się od współczynnika rekombinacji: jest on źródłem, a nie przyczyną strat. W przypadku zależnego od długości fali strumienia fotonów \(\Phi(\lambda,y)\) i współczynnika absorpcji \(\alpha(\lambda,y)\), prosty lokalny związek wyraża się wzorem \(G(y)=\int \alpha(\lambda,y)\Phi(\lambda,y)\,d\lambda\), przy czym dokładne obliczenia zależą od interferencji, odbicia oraz definicji pola wybranego modelu optycznego.Współczynnik generowania nośników optycznych zasadniczo różni się od współczynnika rekombinacji: jest on źródłem, a nie przyczyną strat. W przypadku zależnego od długości fali strumienia fotonów \(\Phi(\lambda,y)\) i współczynnika absorpcji \(\alpha(\lambda,y)\), prosty lokalny związek wyraża się wzorem \(G(y)=\int \alpha(\lambda,y)\Phi(\lambda,y)\,d\lambda\), przy czym dokładne obliczenia zależą od interferencji, odbicia oraz definicji pola wybranego modelu optycznego.
Zapisane profile generacji (?? i ??) pozwalają zobaczyć, gdzie powstają nośniki ładunku generowane przez światło. Silniejsza generacja w pobliżu oświetlonej strony może wynikać z wykładniczej absorpcji; interferencja i zależna od długości fali absorpcja mogą prowadzić do bardziej złożonej zależności od głębokości. Nośnik ładunku generowany optycznie nie jest automatycznie nośnikiem, który został zebrany, dlatego generację należy analizować łącznie z wykresami przedstawiającymi prąd i rekombinację.Zapisane profile generacji (?? i ??) pozwalają zobaczyć, gdzie powstają nośniki ładunku generowane przez światło. Silniejsza generacja w pobliżu oświetlonej strony może wynikać z wykładniczej absorpcji; interferencja i zależna od długości fali absorpcja mogą prowadzić do bardziej złożonej zależności od głębokości. Nośnik ładunku generowany optycznie nie jest automatycznie nośnikiem, który został zebrany, dlatego generację należy analizować łącznie z wykresami przedstawiającymi prąd i rekombinację.


9. Zewnętrzna wydajność kwantowa: zależność od długości fali podczas zbierania nośników ładunku.9. Zewnętrzna wydajność kwantowa: zależność od długości fali podczas zbierania nośników ładunku.
Zewnętrzna wydajność kwantowa (EQE) określa liczbę elektronów zebranych przez zewnętrzny obwód na każdy padający na urządzenie foton o określonej długości fali, przy czym uwzględniany jest znak i sposób normalizacji stosowany w odniesieniu do sygnału wyjściowego. W idealnym modelu, w którym na każdy foton zbierany jest jeden elektron (\(EQE(\lambda)=J_{\mathrm{photo}}(\lambda)/(q\Phi_{\mathrm{incident}}(\lambda))\)), prąd spektralny i strumień fotonów są wyrażone w spójny sposób.Zewnętrzna wydajność kwantowa (EQE) określa liczbę elektronów zebranych przez zewnętrzny obwód na każdy padający na urządzenie foton o określonej długości fali, przy czym uwzględniany jest znak i sposób normalizacji stosowany w odniesieniu do sygnału wyjściowego. W idealnym modelu, w którym na każdy foton zbierany jest jeden elektron (\(EQE(\lambda)=J_{\mathrm{photo}}(\lambda)/(q\Phi_{\mathrm{incident}}(\lambda))\)), prąd spektralny i strumień fotonów są wyrażone w spójny sposób.
Aby uruchomić obliczenia EQE, wybierz opcję symulacji EQE w programie OghmaNano i uruchom ją, a następnie sprawdź zapisany wynik w sekcji Output (??). Dostarczony plik symulacji zawiera również skonfigurowaną długość fali testowej EQE, wynoszącą 710 nm; nie należy mylić tego ustawienia z pełnym, zmierzonym lub obliczonym szerokopasmowym widmem EQE, chyba że faktycznie przeprowadzono skanowanie długości fali, a wynik zawiera to widmo.Aby uruchomić obliczenia EQE, wybierz opcję symulacji EQE w programie OghmaNano i uruchom ją, a następnie sprawdź zapisany wynik w sekcji Output (??). Dostarczony plik symulacji zawiera również skonfigurowaną długość fali testowej EQE, wynoszącą 710 nm; nie należy mylić tego ustawienia z pełnym, zmierzonym lub obliczonym szerokopasmowym widmem EQE, chyba że faktycznie przeprowadzono skanowanie długości fali, a wynik zawiera to widmo.

W przypadku arsenku galu (GaAs) energie fotonów bliskie i wyższe od wartości przerw energetycznych odpowiadają zakresowi bliskiej podczerwieni/widzialnego światła, w którym materiał wykazuje właściwości fotowoltaiczne. Emitter wykonany z AlGaAs, charakteryzujący się większą przerwą energetyczną, może być bardziej przezroczysty w części zakresu długości fal, w którym arsenek galu absorbuje światło, co zwiększa ilość światła dostępnego dla absorbentu. Przy bardzo krótkich długościach fal absorpcja w przedniej warstwie AlGaAs oraz rekombinacja nośników generowanych blisko powierzchni mogą zmniejszyć efektywność zbierania światła. Przy dłuższych długościach fal, w pobliżu granicy absorpcji arsenku galu, niepełna absorpcja w absorberze o określonej grubości (2,8 µm) może stać się istotna. Konkretne wnioski dotyczące widma powinny być wyciągane na podstawie obliczonej krzywej EQE (External Quantum Efficiency) i weryfikowane za pomocą danych optycznych, a nie przypuszczane na podstawie samej wartości przerwy energetycznej.W przypadku arsenku galu (GaAs) energie fotonów bliskie i wyższe od wartości przerw energetycznych odpowiadają zakresowi bliskiej podczerwieni/widzialnego światła, w którym materiał wykazuje właściwości fotowoltaiczne. Emitter wykonany z AlGaAs, charakteryzujący się większą przerwą energetyczną, może być bardziej przezroczysty w części zakresu długości fal, w którym arsenek galu absorbuje światło, co zwiększa ilość światła dostępnego dla absorbentu. Przy bardzo krótkich długościach fal absorpcja w przedniej warstwie AlGaAs oraz rekombinacja nośników generowanych blisko powierzchni mogą zmniejszyć efektywność zbierania światła. Przy dłuższych długościach fal, w pobliżu granicy absorpcji arsenku galu, niepełna absorpcja w absorberze o określonej grubości (2,8 µm) może stać się istotna. Konkretne wnioski dotyczące widma powinny być wyciągane na podstawie obliczonej krzywej EQE (External Quantum Efficiency) i weryfikowane za pomocą danych optycznych, a nie przypuszczane na podstawie samej wartości przerwy energetycznej.
Współczynnik kwantowej efektywności zewnętrznej (EQE) nie jest równoważny współczynnikowi absorpcji optycznej: EQE mierzy liczbę zaabsorbowanych fotonów, natomiast efektywność kwantowa mierzy liczbę zebranych nośników. Porównanie tych dwóch wartości pozwala ocenić straty związane z procesem zbierania ładunków. Jeśli absorpcja jest wysoka, ale EQE jest niska dla określonej długości fali, należy zbadać profil generacji w przestrzeni oraz właściwości warstwy, w której fotony są absorbowane, w zakresie rekombinacji i transportu ładunków.Współczynnik kwantowej efektywności zewnętrznej (EQE) nie jest równoważny współczynnikowi absorpcji optycznej: EQE mierzy liczbę zaabsorbowanych fotonów, natomiast efektywność kwantowa mierzy liczbę zebranych nośników. Porównanie tych dwóch wartości pozwala ocenić straty związane z procesem zbierania ładunków. Jeśli absorpcja jest wysoka, ale EQE jest niska dla określonej długości fali, należy zbadać profil generacji w przestrzeni oraz właściwości warstwy, w której fotony są absorbowane, w zakresie rekombinacji i transportu ładunków.
10. Obliczenia dotyczące natężenia oświetlenia dla modeli Suns–Voc i Suns–Jsc.10. Obliczenia dotyczące natężenia oświetlenia dla modeli Suns–Voc i Suns–Jsc.
10.1. Suns–Voc10.1. Suns–Voc
Suns–Voc bada zależność napięcia w obwodzie otwartym od natężenia oświetlenia. Obwód otwarty definiuje się jako zerowy wypadkowy prąd zewnętrzny, dlatego nośniki generowane przez światło tworzą napięcie zamiast dostarczać energię do obciążenia. Dołączony przykład zawiera konfigurację Suns–Voc, w której natężenie zmienia się od 0,01 do 1000 jednostek natężenia słonecznego, przy czym zmiany następują w krokach o współczynniku 1,2. Uzyskaj dostęp do tego trybu symulacji i sprawdź wybór wyników (??) oraz wykres (??).Suns–Voc bada zależność napięcia w obwodzie otwartym od natężenia oświetlenia. Obwód otwarty definiuje się jako zerowy wypadkowy prąd zewnętrzny, dlatego nośniki generowane przez światło tworzą napięcie zamiast dostarczać energię do obciążenia. Dołączony przykład zawiera konfigurację Suns–Voc, w której natężenie zmienia się od 0,01 do 1000 jednostek natężenia słonecznego, przy czym zmiany następują w krokach o współczynniku 1,2. Uzyskaj dostęp do tego trybu symulacji i sprawdź wybór wyników (??) oraz wykres (??).
W uproszczonym modelu diody, w którym fizyka rekombinacji jest w przybliżeniu niezależna od natężenia oświetlenia, zależność \(V_{OC}\approx(n_{\mathrm{id}}k_BT/q)\ln(J_{\mathrm{ph}}/J_0+1)\) przewiduje w przybliżeniu logarytmiczny wzrost napięcia w obwodzie otwartym wraz ze wzrostem natężenia światła. Tutaj \(n_{\mathrm{id}}\) to efektywny współczynnik idealności. W urządzeniu z heterospójeniem lokalny rozkład nośników, przesunięcia pasm, kanały rekombinacji i zachowanie przy wysokim wstrzykiwaniu mogą zmieniać się wraz z natężeniem, dlatego efektywna wartość nachylenia nie musi być stała w całym zakresie zmian.W uproszczonym modelu diody, w którym fizyka rekombinacji jest w przybliżeniu niezależna od natężenia oświetlenia, zależność \(V_{OC}\approx(n_{\mathrm{id}}k_BT/q)\ln(J_{\mathrm{ph}}/J_0+1)\) przewiduje w przybliżeniu logarytmiczny wzrost napięcia w obwodzie otwartym wraz ze wzrostem natężenia światła. Tutaj \(n_{\mathrm{id}}\) to efektywny współczynnik idealności. W urządzeniu z heterospójeniem lokalny rozkład nośników, przesunięcia pasm, kanały rekombinacji i zachowanie przy wysokim wstrzykiwaniu mogą zmieniać się wraz z natężeniem, dlatego efektywna wartość nachylenia nie musi być stała w całym zakresie zmian.


Metoda Suns–Voc jest przydatna do badania charakterystyki diody przy zerowym prądzie. Sama w sobie nie uwzględnia jednak wszystkich strat napięcia spowodowanych przepływem prądu przez rezystywny kontakt lub warstwę transportową w warunkach generowania energii. Dlatego porównanie zachowania metody Suns–Voc z charakterystyką J–V uzyskaną przy oświetleniu może pomóc w rozróżnieniu ograniczeń związanych z rekombinacją od strat transportowych zależnych od prądu, chociaż interpretacja takiej różnicy wymaga spójności w zakresie oświetlenia, modeli rekombinacji i definicji napięcia.Metoda Suns–Voc jest przydatna do badania charakterystyki diody przy zerowym prądzie. Sama w sobie nie uwzględnia jednak wszystkich strat napięcia spowodowanych przepływem prądu przez rezystywny kontakt lub warstwę transportową w warunkach generowania energii. Dlatego porównanie zachowania metody Suns–Voc z charakterystyką J–V uzyskaną przy oświetleniu może pomóc w rozróżnieniu ograniczeń związanych z rekombinacją od strat transportowych zależnych od prądu, chociaż interpretacja takiej różnicy wymaga spójności w zakresie oświetlenia, modeli rekombinacji i definicji napięcia.
10.2. Suns – Jsc10.2. Suns – Jsc
Suns–Jsc analizuje zależność między prądem zwarciowym a natężeniem oświetlenia. W przykładzie ustawiono zakres od 0,01 jednostek oświetlenia słonecznego do 10 jednostek, przy współczynniku intensywności wynoszącym 1,05 oraz zapisano początkowy przyrost sundp = 0.01. Użyj selektora symulacji Suns–Jsc i sprawdź wynik w ?? i ??.Suns–Jsc analizuje zależność między prądem zwarciowym a natężeniem oświetlenia. W przykładzie ustawiono zakres od 0,01 jednostek oświetlenia słonecznego do 10 jednostek, przy współczynniku intensywności wynoszącym 1,05 oraz zapisano początkowy przyrost sundp = 0.01. Użyj selektora symulacji Suns–Jsc i sprawdź wynik w ?? i ??.
W najprostszym modelu o stałej efektywności zbierania ładunków, prąd zwarciowy zmienia się w przybliżeniu liniowo wraz ze strumieniem padających fotonów, czyli \(|J_{SC}|\propto P_{\mathrm{sun}}\). Systematyczne odchylenie od liniowości może wskazywać na zależność zbierania, rekombinacji lub transportu ładunków od natężenia. Należy jednak wziąć pod uwagę również sposób, w jaki kod zmienia skalę optycznego oświetlenia, oraz numeryczną zbieżność dla każdego punktu.W najprostszym modelu o stałej efektywności zbierania ładunków, prąd zwarciowy zmienia się w przybliżeniu liniowo wraz ze strumieniem padających fotonów, czyli \(|J_{SC}|\propto P_{\mathrm{sun}}\). Systematyczne odchylenie od liniowości może wskazywać na zależność zbierania, rekombinacji lub transportu ładunków od natężenia. Należy jednak wziąć pod uwagę również sposób, w jaki kod zmienia skalę optycznego oświetlenia, oraz numeryczną zbieżność dla każdego punktu.


Oba pomiary przy różnych natężeniach oświetlenia pozwalają odpowiedzieć na różne pytania: pomiar Suns–Voc bada napięcie generowane bez zewnętrznego prądu, natomiast pomiar Suns–Jsc bada prąd generowany przy napięciu bliskim zeru. Należy je interpretować w połączeniu ze standardową krzywą J–V dla jednego źródła światła, a nie traktować jako jej zamiennik.Oba pomiary przy różnych natężeniach oświetlenia pozwalają odpowiedzieć na różne pytania: pomiar Suns–Voc bada napięcie generowane bez zewnętrznego prądu, natomiast pomiar Suns–Jsc bada prąd generowany przy napięciu bliskim zeru. Należy je interpretować w połączeniu ze standardową krzywą J–V dla jednego źródła światła, a nie traktować jako jej zamiennik.
11. Modelowanie optyczne i przestrzenny rozkład generowania fotonów11. Modelowanie optyczne i przestrzenny rozkład generowania fotonów
Rozwiązanie numeryczne dla problemów elektrycznych wymaga określenia przestrzennej szybkości fotogeneracji, która jest dostarczana przez obliczenia optyczne. W przedstawionym przykładzie źródło światła jest skonfigurowane do generowania promieniowania odpowiadającego spektrum słonecznemu AM1.5G, padającego od góry (strona y0), a model generacji opiera się na metodzie macierzy transferu. Spektrum jest dyskretyzowane w zakresie 320–1100 nm, przy czym zdefiniowano 201 punktów długości fali. Oddzielna przestrzennie dyskretyzowana siatka optyczna obejmuje całą powierzchnię urządzenia o wymiarach 3,4 µm i składa się z 200 punktów. Są to dyskretyzacje optyczne, odmienne od dyskretyzacji elektrycznej, o której mowa wcześniej.Rozwiązanie numeryczne dla problemów elektrycznych wymaga określenia przestrzennej szybkości fotogeneracji, która jest dostarczana przez obliczenia optyczne. W przedstawionym przykładzie źródło światła jest skonfigurowane do generowania promieniowania odpowiadającego spektrum słonecznemu AM1.5G, padającego od góry (strona y0), a model generacji opiera się na metodzie macierzy transferu. Spektrum jest dyskretyzowane w zakresie 320–1100 nm, przy czym zdefiniowano 201 punktów długości fali. Oddzielna przestrzennie dyskretyzowana siatka optyczna obejmuje całą powierzchnię urządzenia o wymiarach 3,4 µm i składa się z 200 punktów. Są to dyskretyzacje optyczne, odmienne od dyskretyzacji elektrycznej, o której mowa wcześniej.
Otwórz interfejs optyczny i wybierz dostępne zakładki wyjściowe (??–??). Zrzuty ekranu ilustrują różne widoki graficzne, w tym wykresy, które dla określonych wybranych parametrów mogą zawierać niewiele danych lub być puste. Pusty panel nie powinien być automatycznie interpretowany jako zerowa absorpcja optyczna: najpierw sprawdź, czy odpowiednie obliczenia optyczne zostały wykonane, czy istnieje odpowiedni zestaw danych wyjściowych oraz czy jego zakresy długości fal lub przestrzenne są objęte wybranymi osiami.Otwórz interfejs optyczny i wybierz dostępne zakładki wyjściowe (??–??). Zrzuty ekranu ilustrują różne widoki graficzne, w tym wykresy, które dla określonych wybranych parametrów mogą zawierać niewiele danych lub być puste. Pusty panel nie powinien być automatycznie interpretowany jako zerowa absorpcja optyczna: najpierw sprawdź, czy odpowiednie obliczenia optyczne zostały wykonane, czy istnieje odpowiedni zestaw danych wyjściowych oraz czy jego zakresy długości fal lub przestrzenne są objęte wybranymi osiami.




Metoda macierzy przenoszenia pozwala ocenić, w jaki sposób fale optyczne rozchodzące się w kierunku do przodu i do tyłu oddziałują z warstwami o złożonych współczynnikach załamania \(\tilde n(\lambda)=n(\lambda)+ik(\lambda)\). Rozwiązanie pola uwzględnia odbicie i interferencję; lokalna absorpcja jest następnie przekształcana w szybkość generacji przy użyciu energii fotonów. W przypadku absorbentu GaAs, zależny od długości fali współczynnik absorpcji jest kluczowym parametrem wejściowym, ponieważ określa, w których miejscach w warstwie o grubości 2,8 µm absorbowane są fotony o różnych energiach.Metoda macierzy przenoszenia pozwala ocenić, w jaki sposób fale optyczne rozchodzące się w kierunku do przodu i do tyłu oddziałują z warstwami o złożonych współczynnikach załamania \(\tilde n(\lambda)=n(\lambda)+ik(\lambda)\). Rozwiązanie pola uwzględnia odbicie i interferencję; lokalna absorpcja jest następnie przekształcana w szybkość generacji przy użyciu energii fotonów. W przypadku absorbentu GaAs, zależny od długości fali współczynnik absorpcji jest kluczowym parametrem wejściowym, ponieważ określa, w których miejscach w warstwie o grubości 2,8 µm absorbowane są fotony o różnych energiach.
Należy zachować szczególną ostrożność podczas interpretacji danych optycznych dotyczących AlGaAs. Model materiału elektrycznego wykorzystuje wyraźnie określoną zawartość stopu, natomiast wybrana baza danych materiałów optycznych jest odwoływana niezależnie. Przed użyciem obliczeń do ilościowego projektowania widmowego należy sprawdzić, czy zestaw danych Inorganic/AlGaAs/Adachi faktycznie dostarcza odpowiednich zależności od składu dla odpowiednich długości fal. Zawartość stopu elektrycznego wynosząca 0,30 nie dowodzi, że każda stała optyczna zmienia się prawidłowo wraz ze zmianą składu.Należy zachować szczególną ostrożność podczas interpretacji danych optycznych dotyczących AlGaAs. Model materiału elektrycznego wykorzystuje wyraźnie określoną zawartość stopu, natomiast wybrana baza danych materiałów optycznych jest odwoływana niezależnie. Przed użyciem obliczeń do ilościowego projektowania widmowego należy sprawdzić, czy zestaw danych Inorganic/AlGaAs/Adachi faktycznie dostarcza odpowiednich zależności od składu dla odpowiednich długości fal. Zawartość stopu elektrycznego wynosząca 0,30 nie dowodzi, że każda stała optyczna zmienia się prawidłowo wraz ze zmianą składu.
Podobnie, symulowana jednowymiarowa struktura upraszcza geometrię optyczną. Nie modeluje ona niezależnie wytworzonej metalowej siatki na przedniej powierzchni, rozprzestrzeniania się prądu w płaszczyźnie, struktury powierzchni, szczegółowych powłok antyrefleksyjnych ani defektów wprowadzonych przez rzeczywisty proces epitaksji, chyba że te elementy zostaną wyraźnie uwzględnione w modelach optycznych i elektrycznych. Analiza optyczna jest zatem szczególnie przydatna do zrozumienia trendów i identyfikacji potencjalnych strat spektralnych, natomiast porównanie z rzeczywistymi ogniwami wymaga starannie dobranej struktury i skalibrowanych danych materiałowych.Podobnie, symulowana jednowymiarowa struktura upraszcza geometrię optyczną. Nie modeluje ona niezależnie wytworzonej metalowej siatki na przedniej powierzchni, rozprzestrzeniania się prądu w płaszczyźnie, struktury powierzchni, szczegółowych powłok antyrefleksyjnych ani defektów wprowadzonych przez rzeczywisty proces epitaksji, chyba że te elementy zostaną wyraźnie uwzględnione w modelach optycznych i elektrycznych. Analiza optyczna jest zatem szczególnie przydatna do zrozumienia trendów i identyfikacji potencjalnych strat spektralnych, natomiast porównanie z rzeczywistymi ogniwami wymaga starannie dobranej struktury i skalibrowanych danych materiałowych.
12. Założenia, weryfikacja i ograniczenia12. Założenia, weryfikacja i ograniczenia
Ten przykład należy traktować jako praktyczną, numeryczną demonstrację heterospójenia AlGaAs/GaAs, a nie jako potwierdzoną prognozę dotyczącą rzeczywistej wydajności konkretnej, wyprodukowanej ogniwa słonecznego. Czterowarstwowa struktura, kontakty i siatka pochodzą z dostarczonej konfiguracji symulacji, a zrzuty ekranu pokazują, że przygotowano kilka obliczeń i wyświetleń wyników. Niezawodna, ilościowa prognoza wymaga dodatkowo niezależnie zweryfikowanych materiałów, parametrów optycznych, współczynników rekombinacji, jakości kontaktów oraz zbieżności numerycznej.Ten przykład należy traktować jako praktyczną, numeryczną demonstrację heterospójenia AlGaAs/GaAs, a nie jako potwierdzoną prognozę dotyczącą rzeczywistej wydajności konkretnej, wyprodukowanej ogniwa słonecznego. Czterowarstwowa struktura, kontakty i siatka pochodzą z dostarczonej konfiguracji symulacji, a zrzuty ekranu pokazują, że przygotowano kilka obliczeń i wyświetleń wyników. Niezawodna, ilościowa prognoza wymaga dodatkowo niezależnie zweryfikowanych materiałów, parametrów optycznych, współczynników rekombinacji, jakości kontaktów oraz zbieżności numerycznej.
Geometria i wymiarowość. Widok 3D to wizualizacja; w tym samouczku zależność elektryczna i optyczna od głębokości jest jednowymiarowa. Efekty, takie jak cieniowanie przedniej siatki, straty na krawędziach, boczne spadki napięcia i nierównomierne styki, wykraczają poza zakres tego podstawowego opisu jednowymiarowego.Geometria i wymiarowość. Widok 3D to wizualizacja; w tym samouczku zależność elektryczna i optyczna od głębokości jest jednowymiarowa. Efekty, takie jak cieniowanie przedniej siatki, straty na krawędziach, boczne spadki napięcia i nierównomierne styki, wykraczają poza zakres tego podstawowego opisu jednowymiarowego.
Spójność materiałowa i modelowa. Energia przerwy i powinowactwo elektronowe w stopie, którego właściwości zależą od składu, powinny być określane za pomocą spójnych modeli. Skrypty Lua dotyczące właściwości elektrycznych oraz zbiory danych optycznych \(n,k\) stanowią odrębne dane wejściowe. Jakakolwiek modyfikacja zastosowana tylko do jednej strony może nieumyślnie spowodować powstanie sztucznych nieciągłości, zwłaszcza na granicy między różnymi materiałami.Spójność materiałowa i modelowa. Energia przerwy i powinowactwo elektronowe w stopie, którego właściwości zależą od składu, powinny być określane za pomocą spójnych modeli. Skrypty Lua dotyczące właściwości elektrycznych oraz zbiory danych optycznych \(n,k\) stanowią odrębne dane wejściowe. Jakakolwiek modyfikacja zastosowana tylko do jednej strony może nieumyślnie spowodować powstanie sztucznych nieciągłości, zwłaszcza na granicy między różnymi materiałami.
Rekombinacja i warstwy graniczne. Obecne są parametry opisujące rekombinację w materiale, ale nie dowodzi to, że rzeczywista warstwa graniczna AlGaAs/GaAs została odwzorowana za pomocą eksperymentalnie potwierdzonej prędkości rekombinacji na tej warstwie. Rzeczywista jakość warstwy granicznej może wpływać na mierzone napięcie i efektywność zbierania nośników. Zmiany w pułapkach oraz procesach radiacyjnych i Augerowskich należy testować pojedynczo i analizować za pomocą profili szybkości.Rekombinacja i warstwy graniczne. Obecne są parametry opisujące rekombinację w materiale, ale nie dowodzi to, że rzeczywista warstwa graniczna AlGaAs/GaAs została odwzorowana za pomocą eksperymentalnie potwierdzonej prędkości rekombinacji na tej warstwie. Rzeczywista jakość warstwy granicznej może wpływać na mierzone napięcie i efektywność zbierania nośników. Zmiany w pułapkach oraz procesach radiacyjnych i Augerowskich należy testować pojedynczo i analizować za pomocą profili szybkości.
Możliwość odtworzenia wyników numerycznych. Sprawdź, czy proces obliczeń zbiega się dla całego zakresu parametrów J–V i natężenia, wykonaj zagęszczenie siatki w pobliżu heterointerfejsu i sprawdź, czy kluczowe parametry fotowoltaiczne i wewnętrzne profile są stabilne. Awaria w konkretnym punkcie o wysokim natężeniu nie jest dowodem na występowanie fizycznego efektu w urządzeniu. Użyj niezmienionej wartości początkowej, aby odróżnić rzeczywiste trendy od wrażliwości algorytmu.Możliwość odtworzenia wyników numerycznych. Sprawdź, czy proces obliczeń zbiega się dla całego zakresu parametrów J–V i natężenia, wykonaj zagęszczenie siatki w pobliżu heterointerfejsu i sprawdź, czy kluczowe parametry fotowoltaiczne i wewnętrzne profile są stabilne. Awaria w konkretnym punkcie o wysokim natężeniu nie jest dowodem na występowanie fizycznego efektu w urządzeniu. Użyj niezmienionej wartości początkowej, aby odróżnić rzeczywiste trendy od wrażliwości algorytmu.
Analiza zrzutów ekranu. Każdy zrzut przedstawia wybrany stan z zapisanych danych symulacji, a numery rysunków w tym samouczku odnoszą się do indeksu zrzutu ekranu, a nie do gwarantowanego przyłożonego napięcia. Przed przypisaniem konkretnej wartości napięcia do danego profilu, należy sprawdzić rzeczywistą wartość napięcia wyświetlaną na wykresie OghmaNano. Należy unikać wyciągania wniosków dotyczących wartości liczbowych \(V_{OC}\), \(J_{SC}\), współczynnika wypełnienia lub wydajności na podstawie małego obrazka, gdy plik z danymi wyjściowymi zawiera dokładne wartości.Analiza zrzutów ekranu. Każdy zrzut przedstawia wybrany stan z zapisanych danych symulacji, a numery rysunków w tym samouczku odnoszą się do indeksu zrzutu ekranu, a nie do gwarantowanego przyłożonego napięcia. Przed przypisaniem konkretnej wartości napięcia do danego profilu, należy sprawdzić rzeczywistą wartość napięcia wyświetlaną na wykresie OghmaNano. Należy unikać wyciągania wniosków dotyczących wartości liczbowych \(V_{OC}\), \(J_{SC}\), współczynnika wypełnienia lub wydajności na podstawie małego obrazka, gdy plik z danymi wyjściowymi zawiera dokładne wartości.