Symulacja diody półprzewodnikowej InP typu PN: dryf i dyfuzja, prostowanie i rekombinacja nośników.Symulacja diody półprzewodnikowej InP typu PN: dryf i dyfuzja, prostowanie i rekombinacja nośników.
1. Wprowadzenie: prostująca struktura półprzewodnikowa na bazie InP.1. Wprowadzenie: prostująca struktura półprzewodnikowa na bazie InP.
Fosforek indu (InP) to półprzewodnik z grupy III–V o bezpośredniej przerwie energetycznej, który stanowi podstawę do tworzenia szybkich urządzeń elektronicznych i zintegrowanych technologii fotonicznych. Diodek PN jest jedną z najprostszych i najbardziej pouczających struktur tego typu: obszar typu p styka się z obszarem typu n, tworząc wewnętrzne pole elektryczne i barierę potencjału, które kontrolują ruch elektronów i dziur. Diodek silnie przewodzi prąd przy polaryzacji przewodzącej i tłumi wstrzykiwanie nośników przy polaryzacji blokującej, co prowadzi do prostowania prądu.Fosforek indu (InP) to półprzewodnik z grupy III–V o bezpośredniej przerwie energetycznej, który stanowi podstawę do tworzenia szybkich urządzeń elektronicznych i zintegrowanych technologii fotonicznych. Diodek PN jest jedną z najprostszych i najbardziej pouczających struktur tego typu: obszar typu p styka się z obszarem typu n, tworząc wewnętrzne pole elektryczne i barierę potencjału, które kontrolują ruch elektronów i dziur. Diodek silnie przewodzi prąd przy polaryzacji przewodzącej i tłumi wstrzykiwanie nośników przy polaryzacji blokującej, co prowadzi do prostowania prądu.
InP ma przerwę energetyczną wynoszącą około 1,34 eV w temperaturze pokojowej oraz stosunkowo wysoką ruchliwość elektronów. Dzięki tym właściwościom, złącza InP znajdują zastosowanie w elektronice opartej na związkach półprzewodnikowych, a także w badaniach nad interfejsami półprzewodnikowymi wykorzystywanymi w urządzeniach opartych na InP. Przykład zawarty w tym samouczku ilustruje działanie elektryczne kompaktowego homozłącza p-InP/n-InP z warstwą kontaktową n⁺.InP ma przerwę energetyczną wynoszącą około 1,34 eV w temperaturze pokojowej oraz stosunkowo wysoką ruchliwość elektronów. Dzięki tym właściwościom, złącza InP znajdują zastosowanie w elektronice opartej na związkach półprzewodnikowych, a także w badaniach nad interfejsami półprzewodnikowymi wykorzystywanymi w urządzeniach opartych na InP. Przykład zawarty w tym samouczku ilustruje działanie elektryczne kompaktowego homozłącza p-InP/n-InP z warstwą kontaktową n⁺.
Zamiast dopasowywać krzywą J–V do samego równania idealnej diody, OghmaNano rozwiązuje sprzężone równania Poissona oraz równania ciągłości dotyczące dryfu i dyfuzji. Otrzymane dane ujawniają zarówno prąd na zaciskach, jak i wewnętrzny rozkład energii pasma przewodnictwa i pasma walencyjnego, gęstości elektronów i dziur, poziomów quasi-Fermiego oraz procesów rekombinacji. Dzięki temu możemy wyjaśnić, dlaczego dioda pełni funkcję prostownika, a nie tylko stwierdzić, że tak się dzieje.Zamiast dopasowywać krzywą J–V do samego równania idealnej diody, OghmaNano rozwiązuje sprzężone równania Poissona oraz równania ciągłości dotyczące dryfu i dyfuzji. Otrzymane dane ujawniają zarówno prąd na zaciskach, jak i wewnętrzny rozkład energii pasma przewodnictwa i pasma walencyjnego, gęstości elektronów i dziur, poziomów quasi-Fermiego oraz procesów rekombinacji. Dzięki temu możemy wyjaśnić, dlaczego dioda pełni funkcję prostownika, a nie tylko stwierdzić, że tak się dzieje.
Główne okno symulacji, przedstawione na rysunku ??, wyświetla strukturę InP ułożoną w warstwy pionowo. Chociaż geometria jest wyświetlana w trzech wymiarach, rozwiązywany tutaj problem elektryczny jest jednowymiarowy i dotyczy kierunku wzrostu.Główne okno symulacji, przedstawione na rysunku ??, wyświetla strukturę InP ułożoną w warstwy pionowo. Chociaż geometria jest wyświetlana w trzech wymiarach, rozwiązywany tutaj problem elektryczny jest jednowymiarowy i dotyczy kierunku wzrostu.
2. Tworzenie nowej symulacji.2. Tworzenie nowej symulacji.
W głównym oknie programu OghmaNano wybierz opcję „Nowa symulacja”. Otworzy się okno wyboru symulacji, w którym wyświetlane są dostępne kategorie przykładów (??). Otwórz sekcję „Przykłady III–V/InP” i wybierz przykład „Diody PN InP” (??). Dokładna lokalizacja przykładu może się różnić w zależności od wersji programu OghmaNano.W głównym oknie programu OghmaNano wybierz opcję „Nowa symulacja”. Otworzy się okno wyboru symulacji, w którym wyświetlane są dostępne kategorie przykładów (??). Otwórz sekcję „Przykłady III–V/InP” i wybierz przykład „Diody PN InP” (??). Dokładna lokalizacja przykładu może się różnić w zależności od wersji programu OghmaNano.
Po otwarciu przykładu wybierz zakładkę „Struktura urządzenia” i otwórz edytor warstw. Przed uruchomieniem symulacji sprawdź nazwy warstw, ich grubości oraz powiązane modele materiałów elektrycznych. W poniższych sekcjach wyjaśniono, w jaki sposób poszczególne obszary wpływają na ogólną charakterystykę diody.Po otwarciu przykładu wybierz zakładkę „Struktura urządzenia” i otwórz edytor warstw. Przed uruchomieniem symulacji sprawdź nazwy warstw, ich grubości oraz powiązane modele materiałów elektrycznych. W poniższych sekcjach wyjaśniono, w jaki sposób poszczególne obszary wpływają na ogólną charakterystykę diody.
3. Analiza struktury urządzenia, połączeń i domieszkowania.3. Analiza struktury urządzenia, połączeń i domieszkowania.
3.1 Struktura warstwy półprzewodnikowej3.1 Struktura warstwy półprzewodnikowej
To urządzenie elektryczne zawiera trzy warstwy InP o łącznej grubości warstwy półprzewodnikowej wynoszącej 700 nm. Warstwy p-InP i n-InP tworzą złącze prostujące, natomiast obszar n⁺-InP stanowi bardziej silnie domieszkowaną warstwę po stronie katody. Rozmieszczenie warstw jest widoczne w edytorze (??).To urządzenie elektryczne zawiera trzy warstwy InP o łącznej grubości warstwy półprzewodnikowej wynoszącej 700 nm. Warstwy p-InP i n-InP tworzą złącze prostujące, natomiast obszar n⁺-InP stanowi bardziej silnie domieszkowaną warstwę po stronie katody. Rozmieszczenie warstw jest widoczne w edytorze (??).
| WarstwaWarstwa | GrubośćGrubość | Doping nominalnyDoping nominalny | RolaRola |
|---|---|---|---|
| p-InPp-InP | 250 nm250 nm | \(N_A=5\times10^{17}\,\mathrm{cm^{-3}}\)\(N_A=5\times10^{17}\,\mathrm{cm^{-3}}\) | Obszar z przewagą otworów i strona anody.Obszar z przewagą otworów i strona anody. |
| n-InPn-InP | 350 nm350 nm | \(N_D=5\times10^{17}\,\mathrm{cm^{-3}}\)\(N_D=5\times10^{17}\,\mathrm{cm^{-3}}\) | Obszar, w którym przeważają elektrony, tworzący złącze.Obszar, w którym przeważają elektrony, tworzący złącze. |
| n⁺-InPn⁺-InP | 100 nanometrów100 nanometrów | \(N_D=1\times10^{18}\,\mathrm{cm^{-3}}\)\(N_D=1\times10^{18}\,\mathrm{cm^{-3}}\) | Warstwa po stronie katody o wyższej koncentracji domieszek.Warstwa po stronie katody o wyższej koncentracji domieszek. |
Warstwa powietrza o grubości 200 nm, widoczna w edytorze struktury, stanowi część wyświetlacza i otaczającej go optycznej geometrii, a nie półprzewodnikowej struktury elektrycznej. W związku z tym jej grubość nie jest uwzględniana w całkowitej grubości urządzenia elektrycznego, która wynosi 700 nm.Warstwa powietrza o grubości 200 nm, widoczna w edytorze struktury, stanowi część wyświetlacza i otaczającej go optycznej geometrii, a nie półprzewodnikowej struktury elektrycznej. W związku z tym jej grubość nie jest uwzględniana w całkowitej grubości urządzenia elektrycznego, która wynosi 700 nm.
3.2 Domieszkowanie i wbudowane ładunki elektrostatyczne3.2 Domieszkowanie i wbudowane ładunki elektrostatyczne
Otwórz „Open Electrical” → „Doping / Jony”, aby zobaczyć stężenia donorów i akceptorów (??). Całkowita zmiana stężenia domieszek zmienia znak na granicy między obszarami p-InP i n-InP, w odległości około 250 nm od górnej warstwy elektrycznie aktywnej struktury półprzewodnikowej. Druga zmiana, między obszarami n-InP i n⁺-InP, modyfikuje gęstość elektronów w kierunku katody, ale nie tworzy kolejnego złącza PN.Otwórz „Open Electrical” → „Doping / Jony”, aby zobaczyć stężenia donorów i akceptorów (??). Całkowita zmiana stężenia domieszek zmienia znak na granicy między obszarami p-InP i n-InP, w odległości około 250 nm od górnej warstwy elektrycznie aktywnej struktury półprzewodnikowej. Druga zmiana, między obszarami n-InP i n⁺-InP, modyfikuje gęstość elektronów w kierunku katody, ale nie tworzy kolejnego złącza PN.
W stanie równowagi elektrony mają tendencję do dyfundowania z obszaru typu n do obszaru typu p, natomiast dziury dyfundują w przeciwnym kierunku. Odsłonięte, zjonizowane domieszki generują pole, które przeciwdziała dalszej dyfuzji. Powstały obszar zubożenia i wewnętrzna bariera potencjału stanowią podstawę działania elementu prostowniczego.W stanie równowagi elektrony mają tendencję do dyfundowania z obszaru typu n do obszaru typu p, natomiast dziury dyfundują w przeciwnym kierunku. Odsłonięte, zjonizowane domieszki generują pole, które przeciwdziała dalszej dyfuzji. Powstały obszar zubożenia i wewnętrzna bariera potencjału stanowią podstawę działania elementu prostowniczego.
\[ \frac{d}{dy}\left(\varepsilon\frac{d\phi}{dy}\right)=-\rho(y), \qquad \rho(y)=q\,[p-n+N_D^+-N_A^-]. \]W tym równaniu \(\phi\) oznacza potencjał elektrostatyczny, \(\varepsilon\) względną przenikalność elektryczną, a \(\rho\) lokalną gęstość ładunku. Zjonizowane domieszki i ruchliwe nośniki ładunku determinują lokalne pole elektryczne, natomiast przyłożone napięcie zmienia efektywną barierę dla wstrzykiwania nośników.W tym równaniu \(\phi\) oznacza potencjał elektrostatyczny, \(\varepsilon\) względną przenikalność elektryczną, a \(\rho\) lokalną gęstość ładunku. Zjonizowane domieszki i ruchliwe nośniki ładunku determinują lokalne pole elektryczne, natomiast przyłożone napięcie zmienia efektywną barierę dla wstrzykiwania nośników.
3.3 Kontakty elektryczne3.3 Kontakty elektryczne
W opisywanym modelu na anodzie występuje warstwa przewodząca głównie przez dziury, a na katodzie – warstwa przewodząca głównie przez elektrony. Takie warunki brzegowe są odpowiednie do badania charakterystyki złącza bez wprowadzania dodatkowej bariery potencjału na stykach. W wytworzonych urządzeniach właściwości metalizacji styków i proces technologiczny decydują o tym, jak blisko styki zbliżają się do idealnego zachowania charakterystycznego dla warstwy przewodzącej.W opisywanym modelu na anodzie występuje warstwa przewodząca głównie przez dziury, a na katodzie – warstwa przewodząca głównie przez elektrony. Takie warunki brzegowe są odpowiednie do badania charakterystyki złącza bez wprowadzania dodatkowej bariery potencjału na stykach. W wytworzonych urządzeniach właściwości metalizacji styków i proces technologiczny decydują o tym, jak blisko styki zbliżają się do idealnego zachowania charakterystycznego dla warstwy przewodzącej.
4. Modele materiałowe i parametry transportu.4. Modele materiałowe i parametry transportu.
4.1 Modele materiałowe w języku Lua: skąd pochodzą parametry InP.4.1 Modele materiałowe w języku Lua: skąd pochodzą parametry InP.
Jedną z ważnych cech tego przykładu diody InP jest to, że parametry elektryczne materiału są definiowane za pomocą modeli materiałowych Lua. Zapewnia to wyraźne powiązanie między właściwościami fizycznymi półprzewodnika a rozwiązaniem opartym na mechanizmie dryfu i dyfuzji. Obszary p-InP, n-InP i n⁺-InP odwołują się do tej samej definicji materiału, Inorganic/InP/InP.lua. Ich odmienne właściwości elektryczne wynikają przede wszystkim z profilu domieszkowania, a nie z przypisywania trzech niezależnych zestawów stałych materiałowych.Jedną z ważnych cech tego przykładu diody InP jest to, że parametry elektryczne materiału są definiowane za pomocą modeli materiałowych Lua. Zapewnia to wyraźne powiązanie między właściwościami fizycznymi półprzewodnika a rozwiązaniem opartym na mechanizmie dryfu i dyfuzji. Obszary p-InP, n-InP i n⁺-InP odwołują się do tej samej definicji materiału, Inorganic/InP/InP.lua. Ich odmienne właściwości elektryczne wynikają przede wszystkim z profilu domieszkowania, a nie z przypisywania trzech niezależnych zestawów stałych materiałowych.
Aby przeanalizować model, otwórz okno Struktura urządzenia → Parametry elektryczne (edytor materiałów elektrycznych) i wybierz jedną z warstw InP. W zakładce Edytor kodu wyświetlany jest skrypt materiału Lua powiązany z danym obszarem; w zakładce Parametry wyświetlane są obliczone właściwości materiału. Należy korzystać z obu tych widoków jednocześnie. Kod definiuje sposób uzyskiwania właściwości, natomiast obliczone wartości pokazują, jakie dane faktycznie otrzymuje solver elektryczny.Aby przeanalizować model, otwórz okno Struktura urządzenia → Parametry elektryczne (edytor materiałów elektrycznych) i wybierz jedną z warstw InP. W zakładce Edytor kodu wyświetlany jest skrypt materiału Lua powiązany z danym obszarem; w zakładce Parametry wyświetlane są obliczone właściwości materiału. Należy korzystać z obu tych widoków jednocześnie. Kod definiuje sposób uzyskiwania właściwości, natomiast obliczone wartości pokazują, jakie dane faktycznie otrzymuje solver elektryczny.
Inorganic/InP/InP.lua.Ten sam schemat pracy edytora kodu można zastosować do innego obszaru półprzewodnikowego. W tym przykładzie z wykorzystaniem InP sprawdź, czy każda aktywna warstwa odnosi się do Inorganic/InP/InP.lua.Obrazy w plikach ?? i ?? pokazują, jak uzyskuje się dostęp do edytora kodu. Konfiguracja warstw symulacji jest podstawą dla przypisanego modelu materiału: każda aktywna область w tym przykładzie jest powiązana z plikiem InP Lua. Nazwy plików obrazów są dziedziczone z sekwencji tworzenia zrzutów ekranu i nie określają, który półprzewodnik jest symulowany.Obrazy w plikach ?? i ?? pokazują, jak uzyskuje się dostęp do edytora kodu. Konfiguracja warstw symulacji jest podstawą dla przypisanego modelu materiału: każda aktywna область w tym przykładzie jest powiązana z plikiem InP Lua. Nazwy plików obrazów są dziedziczone z sekwencji tworzenia zrzutów ekranu i nie określają, który półprzewodnik jest symulowany.
W dostarczonej symulacji przełączniki db_overrides są ustawione na wartość false dla kluczowych parametrów transportu i pasma. W związku z tym OghmaNano jest skonfigurowane tak, aby pobierało odpowiednie wartości z bazy danych materiałów/modelu Lua, zamiast wymuszać użycie wartości domyślnych zapisanych w innym miejscu w pliku JSON. Ma to znaczenie podczas sprawdzania konfiguracji: wyświetlane w zakładce Parametry wartości są bardziej informatywne niż zapisane wartości domyślne w nieaktywnym polu zastępowania.W dostarczonej symulacji przełączniki db_overrides są ustawione na wartość false dla kluczowych parametrów transportu i pasma. W związku z tym OghmaNano jest skonfigurowane tak, aby pobierało odpowiednie wartości z bazy danych materiałów/modelu Lua, zamiast wymuszać użycie wartości domyślnych zapisanych w innym miejscu w pliku JSON. Ma to znaczenie podczas sprawdzania konfiguracji: wyświetlane w zakładce Parametry wartości są bardziej informatywne niż zapisane wartości domyślne w nieaktywnym polu zastępowania.
4.2 Struktura pasmowa, gęstość stanów i ruchliwość.4.2 Struktura pasmowa, gęstość stanów i ruchliwość.
Pierwsze parametry, które należy sprawdzić w zakładce „Parametry obliczone”, to energia przerwy \(E_g\) i powinowactwo elektronowe \(\chi\). Energia przerwy określa różnicę energii między pasmem przewodnictwa i pasmem walencyjnym, natomiast powinowactwo elektronowe wskazuje położenie krawędzi pasma przewodnictwa względem poziomu energii w próżni. W przypadku InP, energia przerwy w temperaturze pokojowej wynosi około \(1.34\,\mathrm{eV}\). Przydatnym wstępnym sprawdzeniem jest porównanie różnicy między \(E_c\) i \(E_v\) na otrzymanym diagramie pasm z obliczoną wartością energii przerwy.Pierwsze parametry, które należy sprawdzić w zakładce „Parametry obliczone”, to energia przerwy \(E_g\) i powinowactwo elektronowe \(\chi\). Energia przerwy określa różnicę energii między pasmem przewodnictwa i pasmem walencyjnym, natomiast powinowactwo elektronowe wskazuje położenie krawędzi pasma przewodnictwa względem poziomu energii w próżni. W przypadku InP, energia przerwy w temperaturze pokojowej wynosi około \(1.34\,\mathrm{eV}\). Przydatnym wstępnym sprawdzeniem jest porównanie różnicy między \(E_c\) i \(E_v\) na otrzymanym diagramie pasm z obliczoną wartością energii przerwy.
Efektywne gęstości stanów w paśmie przewodnictwa i paśmie walencyjnym, oznaczone odpowiednio jako \(N_c\) i \(N_v\), wpływają na zależność między koncentracją nośników ładunku a położeniem quasi-poziomu Fermiego. W związku z tym wpływają one na równowagowe populacje elektronów i dziur oraz na odpowiedź na wstrzykiwanie ładunku przy polaryzacji przewodzącej. Względna przenikalność \(\varepsilon_r\) występuje w równaniu Poissona i wpływa na potencjał elektrostatyczny oraz szerokość warstwy zubożonej. Wielkości te należy analizować jako spójny zestaw, a nie jako niezależne parametry do regulacji.Efektywne gęstości stanów w paśmie przewodnictwa i paśmie walencyjnym, oznaczone odpowiednio jako \(N_c\) i \(N_v\), wpływają na zależność między koncentracją nośników ładunku a położeniem quasi-poziomu Fermiego. W związku z tym wpływają one na równowagowe populacje elektronów i dziur oraz na odpowiedź na wstrzykiwanie ładunku przy polaryzacji przewodzącej. Względna przenikalność \(\varepsilon_r\) występuje w równaniu Poissona i wpływa na potencjał elektrostatyczny oraz szerokość warstwy zubożonej. Wielkości te należy analizować jako spójny zestaw, a nie jako niezależne parametry do regulacji.
Ruchliwość elektronów i dziur, oznaczana odpowiednio jako \(\mu_n\) i \(\mu_p\), określa odpowiedź nośników ładunku na działanie pola elektrycznego; dyfuzja zależy dodatkowo od gradientów gęstości nośników. W modelu jednowymiarowym gęstości prądu zapisuje się w następujący sposób:Ruchliwość elektronów i dziur, oznaczana odpowiednio jako \(\mu_n\) i \(\mu_p\), określa odpowiedź nośników ładunku na działanie pola elektrycznego; dyfuzja zależy dodatkowo od gradientów gęstości nośników. W modelu jednowymiarowym gęstości prądu zapisuje się w następujący sposób:
\[ J_n=q\mu_n nE+qD_n\frac{dn}{dy},\qquad J_p=q\mu_p pE-qD_p\frac{dp}{dy}. \]W tym równaniu \(n\) i \(p\) oznaczają lokalne gęstości elektronów i dziur, \(E\) to natężenie pola elektrycznego, a \(D_n\) i \(D_p\) to współczynniki dyfuzji. Dwie składowe prądu mogą znacznie się różnić w różnych punktach diody, nawet gdy całkowity prąd w stanie ustalonym jest w przybliżeniu niezależny od położenia. Dlatego w dalszych częściach tekstu analizowane są \(J_n\) i \(J_p\) wraz z profilami gęstości nośników.W tym równaniu \(n\) i \(p\) oznaczają lokalne gęstości elektronów i dziur, \(E\) to natężenie pola elektrycznego, a \(D_n\) i \(D_p\) to współczynniki dyfuzji. Dwie składowe prądu mogą znacznie się różnić w różnych punktach diody, nawet gdy całkowity prąd w stanie ustalonym jest w przybliżeniu niezależny od położenia. Dlatego w dalszych częściach tekstu analizowane są \(J_n\) i \(J_p\) wraz z profilami gęstości nośników.
4.3 Jak sprawdzić i porównać ustawienia materiałów.4.3 Jak sprawdzić i porównać ustawienia materiałów.
W edytorze elektrycznym wybierz materiał p-InP i przeanalizuj jego kod Lua oraz wartości parametrów. Następnie przejdź do materiału n-InP, a na końcu do n⁺-InP. Porównaj przypisany model materiału, energię przerwy, powinowactwo elektronowe, ruchliwość, efektywną gęstość stanów, przenikalność elektryczną i współczynniki rekombinacji. W przypadku tego homozłącza we wszystkich trzech obszarach stosuje się ten sam podstawowy model InP; główną celową zmianą między obszarami jest ich koncentracja akceptorów lub donorów.W edytorze elektrycznym wybierz materiał p-InP i przeanalizuj jego kod Lua oraz wartości parametrów. Następnie przejdź do materiału n-InP, a na końcu do n⁺-InP. Porównaj przypisany model materiału, energię przerwy, powinowactwo elektronowe, ruchliwość, efektywną gęstość stanów, przenikalność elektryczną i współczynniki rekombinacji. W przypadku tego homozłącza we wszystkich trzech obszarach stosuje się ten sam podstawowy model InP; główną celową zmianą między obszarami jest ich koncentracja akceptorów lub donorów.
Warto rozróżnić dwie koncepcje: model materiałowy opisuje właściwości InP, natomiast profil domieszkowania określa, które nośniki ładunku dominują lokalnie i gdzie znajduje się złącze. To rozróżnienie staje się szczególnie przydatne przy przejściu od homozłącza PN do heterozłącza opartego na InP, gdzie różnice między modelami materiałowymi mogą same w sobie generować przesunięcia pasm energetycznych. W tej diodzie, zakrzywienie pasma przewodnictwa i pasma walencyjnego, widoczne na zdjęciach, odzwierciedla przede wszystkim elektrostatykę w obszarze przejścia p/n.Warto rozróżnić dwie koncepcje: model materiałowy opisuje właściwości InP, natomiast profil domieszkowania określa, które nośniki ładunku dominują lokalnie i gdzie znajduje się złącze. To rozróżnienie staje się szczególnie przydatne przy przejściu od homozłącza PN do heterozłącza opartego na InP, gdzie różnice między modelami materiałowymi mogą same w sobie generować przesunięcia pasm energetycznych. W tej diodzie, zakrzywienie pasma przewodnictwa i pasma walencyjnego, widoczne na zdjęciach, odzwierciedla przede wszystkim elektrostatykę w obszarze przejścia p/n.
Aby zapewnić powtarzalność badań, należy utrzymać stały model materiału Lua podczas zmiany stężeń domieszek. Z kolei, jeśli analizujemy, jak wynik zależy od współczynnika ruchliwości lub rekombinacji, należy celowo zmieniać ten parametr materiału i porównywać uzyskane krzywe J–V, gęstości nośników i profile prądu wewnętrznego. Edytor parametrów i widok kodu Lua zapewniają połączenie między założeniami dotyczącymi danych wejściowych a symulowanymi wynikami.Aby zapewnić powtarzalność badań, należy utrzymać stały model materiału Lua podczas zmiany stężeń domieszek. Z kolei, jeśli analizujemy, jak wynik zależy od współczynnika ruchliwości lub rekombinacji, należy celowo zmieniać ten parametr materiału i porównywać uzyskane krzywe J–V, gęstości nośników i profile prądu wewnętrznego. Edytor parametrów i widok kodu Lua zapewniają połączenie między założeniami dotyczącymi danych wejściowych a symulowanymi wynikami.
4.4 Parametry rekombinacji4.4 Parametry rekombinacji
Obszary InP obejmują rekombinację Shockleya–Reada–Halla (SRH) oraz kanał rekombinacji nośników swobodnych, uwzględnione w opracowanym modelu materiału. Rekombinacja SRH opisuje proces wychwytywania elektronów i dziur z udziałem defektów; rekombinacja nośników swobodnych jest szczególnie istotna w InP o bezpośredniej przerwie energetycznej, ponieważ przejścia elektron-dziura mogą zachodzić w sposób emisyjny. W niniejszym przykładzie oba procesy są przydatne jako mechanizmy strat nośników, które modyfikują prąd elektryczny.Obszary InP obejmują rekombinację Shockleya–Reada–Halla (SRH) oraz kanał rekombinacji nośników swobodnych, uwzględnione w opracowanym modelu materiału. Rekombinacja SRH opisuje proces wychwytywania elektronów i dziur z udziałem defektów; rekombinacja nośników swobodnych jest szczególnie istotna w InP o bezpośredniej przerwie energetycznej, ponieważ przejścia elektron-dziura mogą zachodzić w sposób emisyjny. W niniejszym przykładzie oba procesy są przydatne jako mechanizmy strat nośników, które modyfikują prąd elektryczny.
\[ R_{\mathrm{rad}}=B\left(np-n_i^2\right), \qquad R_{\mathrm{SRH}}=\frac{np-n_i^2}{\tau_p(n+n_1)+\tau_n(p+p_1)}. \]Współczynnik \(B\), czas życia nośników w pułapkach \(\tau_n\) i \(\tau_p\) oraz dodatkowe gęstości \(n_1\) i \(p_1\) określają intensywność każdego z procesów. Kluczowym elementem analizy w tym samouczku jest to, jak zmienia się przestrzenne rozmieszczenie rekombinacji w zależności od tego, jak polaryzacja powoduje napływ większej liczby elektronów i dziur do tego samego obszaru.Współczynnik \(B\), czas życia nośników w pułapkach \(\tau_n\) i \(\tau_p\) oraz dodatkowe gęstości \(n_1\) i \(p_1\) określają intensywność każdego z procesów. Kluczowym elementem analizy w tym samouczku jest to, jak zmienia się przestrzenne rozmieszczenie rekombinacji w zależności od tego, jak polaryzacja powoduje napływ większej liczby elektronów i dziur do tego samego obszaru.
4.5 Siatka obliczeniowa i rozwiązanie numeryczne4.5 Siatka obliczeniowa i rozwiązanie numeryczne
W obliczeniach elektrycznych wykorzystywana jest jednowymiarowa siatka składająca się z 150 punktów rozmieszczonych na długości 700 nm, co odpowiada średniej odległości około 4,7 nm. Powiązane nieliniowe równania są rozwiązywane za pomocą metody iteracyjnej Newtona oraz rzadkiego liniowego solvera UMFPACK. Ta siatka stanowi użyteczny punkt wyjścia do analizy obszaru złącza i obszarów quasi-neutralnych.W obliczeniach elektrycznych wykorzystywana jest jednowymiarowa siatka składająca się z 150 punktów rozmieszczonych na długości 700 nm, co odpowiada średniej odległości około 4,7 nm. Powiązane nieliniowe równania są rozwiązywane za pomocą metody iteracyjnej Newtona oraz rzadkiego liniowego solvera UMFPACK. Ta siatka stanowi użyteczny punkt wyjścia do analizy obszaru złącza i obszarów quasi-neutralnych.
Wyniki zawierają informacje o zbieżności oraz profile urządzenia zależne od przesunięcia. Dobrym pierwszym krokiem jest sprawdzenie, czy charakterystyka prądowo-napięciowa (J–V) jest płynna, a następnie potwierdzenie, czy prądy elektronów i dziur spełniają ogólne zasady zachowania prądu w stanie ustalonym. Edytor siatki jest wyświetlany w ??.Wyniki zawierają informacje o zbieżności oraz profile urządzenia zależne od przesunięcia. Dobrym pierwszym krokiem jest sprawdzenie, czy charakterystyka prądowo-napięciowa (J–V) jest płynna, a następnie potwierdzenie, czy prądy elektronów i dziur spełniają ogólne zasady zachowania prądu w stanie ustalonym. Edytor siatki jest wyświetlany w ??.
5. Uruchomienie symulacji i analiza wyników w postaci wykresu zależności prądu od napięcia (J–V).5. Uruchomienie symulacji i analiza wyników w postaci wykresu zależności prądu od napięcia (J–V).
Wróć do głównego okna i wybierz opcję „Uruchom symulację”. OghmaNano zmienia wartość przyłożonego napięcia i zapisuje uzyskane rozwiązanie w każdym punkcie. W przykładowym pliku zapisanym w pamięci zastosowano zakres napięcia od 0 do 3 V, z krokiem 25 mV. W oknie terminala rejestrowany jest postęp obliczeń podczas zmiany napięcia, a w zakładce „Wyniki” można uzyskać dostęp do uzyskanych danych.Wróć do głównego okna i wybierz opcję „Uruchom symulację”. OghmaNano zmienia wartość przyłożonego napięcia i zapisuje uzyskane rozwiązanie w każdym punkcie. W przykładowym pliku zapisanym w pamięci zastosowano zakres napięcia od 0 do 3 V, z krokiem 25 mV. W oknie terminala rejestrowany jest postęp obliczeń podczas zmiany napięcia, a w zakładce „Wyniki” można uzyskać dostęp do uzyskanych danych.
Po zakończeniu symulacji wyniki stają się dostępne w zakładce „Wyniki” (??). Otwórz jv.csv, aby wyświetlić wykres gęstości prądu w funkcji przyłożonego napięcia. Na zrzutach ekranu znajdują się zarówno wykres liniowy J–V (??), jak i wykres logarytmiczny prądu (??).Po zakończeniu symulacji wyniki stają się dostępne w zakładce „Wyniki” (??). Otwórz jv.csv, aby wyświetlić wykres gęstości prądu w funkcji przyłożonego napięcia. Na zrzutach ekranu znajdują się zarówno wykres liniowy J–V (??), jak i wykres logarytmiczny prądu (??).
5.1 Zrozumienie charakterystyki.5.1 Zrozumienie charakterystyki.
Przy niskim przyłożonym napięciu wbudowana bariera potencjału znacznie ogranicza wstrzykiwanie nośników. Zwiększenie napięcia przewodzenia zmniejsza tę barierę, dzięki czemu elektrony i dziury łatwiej przechodzą na przeciwną stronę, co powoduje wzrost całkowitego prądu. W pośrednim zakresie napięć prąd często rośnie w przybliżeniu w sposób wykładniczy; przy większych prądach transport i straty rezystywne mogą modyfikować nachylenie krzywej.Przy niskim przyłożonym napięciu wbudowana bariera potencjału znacznie ogranicza wstrzykiwanie nośników. Zwiększenie napięcia przewodzenia zmniejsza tę barierę, dzięki czemu elektrony i dziury łatwiej przechodzą na przeciwną stronę, co powoduje wzrost całkowitego prądu. W pośrednim zakresie napięć prąd często rośnie w przybliżeniu w sposób wykładniczy; przy większych prądach transport i straty rezystywne mogą modyfikować nachylenie krzywej.
\[ J\approx J_0\left[\exp\!\left(\frac{qV}{nk_{\mathrm B}T}\right)-1\right]. \]To przybliżone równanie opisujące idealny diod jest przydatne jako punkt odniesienia przy analizie w przybliżeniu eksponencjalnego przedziału, ale nie wynika ono z modelu dryfu i dyfuzji. Współczynnik idealności \(n\) i współczynnik \(J_0\) powinny być określane tylko w odpowiednim zakresie polaryzacji; mogą się zmieniać w zależności od rekombinacji, rezystancji szeregowej i właściwości kontaktu.To przybliżone równanie opisujące idealny diod jest przydatne jako punkt odniesienia przy analizie w przybliżeniu eksponencjalnego przedziału, ale nie wynika ono z modelu dryfu i dyfuzji. Współczynnik idealności \(n\) i współczynnik \(J_0\) powinny być określane tylko w odpowiednim zakresie polaryzacji; mogą się zmieniać w zależności od rekombinacji, rezystancji szeregowej i właściwości kontaktu.
Wykres w skali logarytmicznej pomaga zidentyfikować moment, w którym pojawia się znaczący prąd, oraz zmiany nachylenia. Wykres liniowy ułatwia interpretację wartości i kształtu charakterystyki przy wysokich prądach. Aby uzyskać jak najpełniejszy obraz działania diody, należy przeanalizować oba wykresy łącznie z wewnętrznymi profilami pasm i gęstości nośników.Wykres w skali logarytmicznej pomaga zidentyfikować moment, w którym pojawia się znaczący prąd, oraz zmiany nachylenia. Wykres liniowy ułatwia interpretację wartości i kształtu charakterystyki przy wysokich prądach. Aby uzyskać jak najpełniejszy obraz działania diody, należy przeanalizować oba wykresy łącznie z wewnętrznymi profilami pasm i gęstości nośników.
6. Budowa wewnętrzna: zasada działania diody.6. Budowa wewnętrzna: zasada działania diody.
Folder „snapshots” zawiera profile przestrzenne zapisane w wybranych punktach podczas pomiaru charakterystyki prądowo-napięciowej (J–V). Dołączone zrzuty ekranu przedstawiają próbki z czterech różnych punktów pomiaru. Pokazują one zmianę od niskiego poziomu wstrzykiwania do wyższego napięcia polaryzacji. Cztery obrazy w każdej grupie odpowiadają tym samym wybranym pozycjom suwaka, co umożliwia bezpośrednie porównanie różnych wielkości fizycznych.Folder „snapshots” zawiera profile przestrzenne zapisane w wybranych punktach podczas pomiaru charakterystyki prądowo-napięciowej (J–V). Dołączone zrzuty ekranu przedstawiają próbki z czterech różnych punktów pomiaru. Pokazują one zmianę od niskiego poziomu wstrzykiwania do wyższego napięcia polaryzacji. Cztery obrazy w każdej grupie odpowiadają tym samym wybranym pozycjom suwaka, co umożliwia bezpośrednie porównanie różnych wielkości fizycznych.
6.1 Pasma przewodnictwa i walencyjne oraz quasi-poziomy Fermiego6.1 Pasma przewodnictwa i walencyjne oraz quasi-poziomy Fermiego
Otwórz program do przeglądania obrazów i wybierz opcje Ec.csv, Ev.csv, Fn.csv i Fp.csv. Oznaczają one odpowiednio: krawędź pasma przewodnictwa, krawędź pasma walencyjnego, quasi-poziom Fermiego dla elektronów i quasi-poziom Fermiego dla dziur. Ich przestrzenne rozdzielenie i nachylenie ujawniają potencjalną barierę, reakcję na przyłożone napięcie oraz odchylenie od stanu równowagi.Otwórz program do przeglądania obrazów i wybierz opcje Ec.csv, Ev.csv, Fn.csv i Fp.csv. Oznaczają one odpowiednio: krawędź pasma przewodnictwa, krawędź pasma walencyjnego, quasi-poziom Fermiego dla elektronów i quasi-poziom Fermiego dla dziur. Ich przestrzenne rozdzielenie i nachylenie ujawniają potencjalną barierę, reakcję na przyłożone napięcie oraz odchylenie od stanu równowagi.
Przy niskim prądzie wtrysku krawędzie pasma energetycznego wyginają się wzdłuż złącza w wyniku wbudowanego elektrostatycznego pola. Wraz ze wzrostem napięcia przewodzenia, efektywna bariera zmniejsza się, a quasi-poziomy Fermiego dla elektronów i dziur oddzielają się, ponieważ nośniki są wprowadzane do urządzenia. Kluczowe jest to, że InP pozostaje tym samym półprzewodnikiem po obu stronach złącza: zmieniające się profile odzwierciedlają zdefiniowaną przez domieszkowanie elektrostatykę i wtrysk nośników.Przy niskim prądzie wtrysku krawędzie pasma energetycznego wyginają się wzdłuż złącza w wyniku wbudowanego elektrostatycznego pola. Wraz ze wzrostem napięcia przewodzenia, efektywna bariera zmniejsza się, a quasi-poziomy Fermiego dla elektronów i dziur oddzielają się, ponieważ nośniki są wprowadzane do urządzenia. Kluczowe jest to, że InP pozostaje tym samym półprzewodnikiem po obu stronach złącza: zmieniające się profile odzwierciedlają zdefiniowaną przez domieszkowanie elektrostatykę i wtrysk nośników.
Porównanie ?? i ?? pokazuje, jak polaryzacja końcowa wpływa na wewnętrzną strukturę pasmową. W szczególności, poziomy quasi-Fermiego pozwalają na bezpośrednią analizę stanu niestabilnego, który powstaje podczas przewodzenia w stanie ustalonym.Porównanie ?? i ?? pokazuje, jak polaryzacja końcowa wpływa na wewnętrzną strukturę pasmową. W szczególności, poziomy quasi-Fermiego pozwalają na bezpośrednią analizę stanu niestabilnego, który powstaje podczas przewodzenia w stanie ustalonym.
6.2 Gęstość elektronów i dziur.6.2 Gęstość elektronów i dziur.
Następnie przedstawiono wykresy Q_nfree.csv i Q_pfree.csv, które ilustrują przestrzenny rozkład swobodnych elektronów i dziur. Przy niskim napięciu polaryzacyjnym elektrony przeważają w obszarze typu n, a dziury w obszarze typu p. Ich gęstość gwałtownie maleje w pobliżu złącza, co odzwierciedla obszar zubożenia.Następnie przedstawiono wykresy Q_nfree.csv i Q_pfree.csv, które ilustrują przestrzenny rozkład swobodnych elektronów i dziur. Przy niskim napięciu polaryzacyjnym elektrony przeważają w obszarze typu n, a dziury w obszarze typu p. Ich gęstość gwałtownie maleje w pobliżu złącza, co odzwierciedla obszar zubożenia.
Wraz ze wzrostem napięcia polaryzacji przewodzącej, nośniki mniejszościowe są wprowadzane do przeciwległych obszarów. Powoduje to zwiększenie nakładania się populacji elektronów i dziur oraz zmianę kształtu obu profili gęstości. Przy większych poziomach wprowadzania nośników, populacje nośników swobodnych po obu stronach złącza znacznie różnią się od wartości obserwowanych przy niskim napięciu polaryzacji.Wraz ze wzrostem napięcia polaryzacji przewodzącej, nośniki mniejszościowe są wprowadzane do przeciwległych obszarów. Powoduje to zwiększenie nakładania się populacji elektronów i dziur oraz zmianę kształtu obu profili gęstości. Przy większych poziomach wprowadzania nośników, populacje nośników swobodnych po obu stronach złącza znacznie różnią się od wartości obserwowanych przy niskim napięciu polaryzacji.
Porównaj punkt zmiany gęstości nośników z położeniem przejścia domieszkowania, a następnie z miejscami, w których występuje silne zakrzywienie pasma energetycznego. Pozwoli to powiązać nośniki dostarczane przez kontakt z barierą elektrostatyczną tworzoną przez złącze PN.Porównaj punkt zmiany gęstości nośników z położeniem przejścia domieszkowania, a następnie z miejscami, w których występuje silne zakrzywienie pasma energetycznego. Pozwoli to powiązać nośniki dostarczane przez kontakt z barierą elektrostatyczną tworzoną przez złącze PN.
6.3 Gęstości prądu elektronów i dziur.6.3 Gęstości prądu elektronów i dziur.
Na wykresach Jn.csv i Jp.csv przedstawiono, jak elektrony i dziury przyczyniają się do przepływu prądu elektrycznego przez urządzenie. Poszczególne składowe mogą się zmieniać w zależności od położenia, ponieważ proces rekombinacji powoduje transfer nośników między populacją elektronów i dziur, ale całkowity prąd w stanie ustalonym powinien pozostać w przybliżeniu stały, gdy nie dochodzi do akumulacji ładunku netto.Na wykresach Jn.csv i Jp.csv przedstawiono, jak elektrony i dziury przyczyniają się do przepływu prądu elektrycznego przez urządzenie. Poszczególne składowe mogą się zmieniać w zależności od położenia, ponieważ proces rekombinacji powoduje transfer nośników między populacją elektronów i dziur, ale całkowity prąd w stanie ustalonym powinien pozostać w przybliżeniu stały, gdy nie dochodzi do akumulacji ładunku netto.
Przy bardzo niskim prądzie, kompensacja pomiędzy składnikami dryfu i dyfuzji może powodować, że poszczególne profile prądowe będą wyglądać na zakłócone lub ostre. Zrzuty ekranu uzyskane przy niskim napięciu pokazują to zjawisko wyraźniej niż zrzuty uzyskane przy wyższym prądzie. Warto przeanalizować profile w odpowiedniej skali i porównać całkowity prąd z wartością J–V na końcu zakresu, zanim zaczniemy interpretować drobne szczegóły numeryczne.Przy bardzo niskim prądzie, kompensacja pomiędzy składnikami dryfu i dyfuzji może powodować, że poszczególne profile prądowe będą wyglądać na zakłócone lub ostre. Zrzuty ekranu uzyskane przy niskim napięciu pokazują to zjawisko wyraźniej niż zrzuty uzyskane przy wyższym prądzie. Warto przeanalizować profile w odpowiedniej skali i porównać całkowity prąd z wartością J–V na końcu zakresu, zanim zaczniemy interpretować drobne szczegóły numeryczne.
6.4 Profil rekombinacji6.4 Profil rekombinacji
Proces rekombinacji decyduje o tym, w których miejscach elektrony i dziury, wprowadzone do materiału, znikają z odpowiednich populacji nośników. Otwórz okno przeglądarki obrazów, używając opcji R_nfree_to_pfree.csv i R_ss_srh.csv, aby porównać odpowiednio bezpośrednie kanały rekombinacji (bez udziału defektów) i kanały rekombinacji wspomagane przez defekty. Oba procesy zależą od dostępności elektronów i dziur, dlatego ich rozkłady przestrzenne zmieniają się wraz ze zmianą populacji nośników.Proces rekombinacji decyduje o tym, w których miejscach elektrony i dziury, wprowadzone do materiału, znikają z odpowiednich populacji nośników. Otwórz okno przeglądarki obrazów, używając opcji R_nfree_to_pfree.csv i R_ss_srh.csv, aby porównać odpowiednio bezpośrednie kanały rekombinacji (bez udziału defektów) i kanały rekombinacji wspomagane przez defekty. Oba procesy zależą od dostępności elektronów i dziur, dlatego ich rozkłady przestrzenne zmieniają się wraz ze zmianą populacji nośników.
Na początku przewodzenia rekombinacja może być skoncentrowana w pobliżu złącza, gdzie nakładają się na siebie populacje elektronów i dziur. Przy silniejszym wstrzykiwaniu nośników w kierunku przewodzenia, obszar nakładania się rozszerza, obejmując większą część półprzewodnika, a profile rekombinacji mogą się poszerzyć. Porównaj położenie i szerokość szczytów rekombinacji z gęstością nośników przy tym samym napięciu polaryzacyjnym.Na początku przewodzenia rekombinacja może być skoncentrowana w pobliżu złącza, gdzie nakładają się na siebie populacje elektronów i dziur. Przy silniejszym wstrzykiwaniu nośników w kierunku przewodzenia, obszar nakładania się rozszerza, obejmując większą część półprzewodnika, a profile rekombinacji mogą się poszerzyć. Porównaj położenie i szerokość szczytów rekombinacji z gęstością nośników przy tym samym napięciu polaryzacyjnym.
Porównując oba kanały, należy wziąć pod uwagę zarówno wartość maksymalną, jak i całkowitą wartość w przestrzeni. Wyraźny lokalny pik nie musi dominować nad całkowitą rekombinacją, jeśli inny kanał jest aktywny na znacznie większej objętości. Jest to szczególnie ważne przy analizie profili rekombinacji w odniesieniu do całkowitego prądu elektrycznego.Porównując oba kanały, należy wziąć pod uwagę zarówno wartość maksymalną, jak i całkowitą wartość w przestrzeni. Wyraźny lokalny pik nie musi dominować nad całkowitą rekombinacją, jeśli inny kanał jest aktywny na znacznie większej objętości. Jest to szczególnie ważne przy analizie profili rekombinacji w odniesieniu do całkowitego prądu elektrycznego.
6.5 Profile rozkładu gęstości elektronów i dziur.6.5 Profile rozkładu gęstości elektronów i dziur.
Zdjęcia oznaczone jako R_n.csv i R_p.csv zawierają dodatkowe informacje na temat współczynników używanych w równaniach ciągłości elektronów i dziur. Pozwalają one zidentyfikować obszary, w których równowaga między lokalnymi składnikami generacji, wychwytu i rekombinacji jest najbardziej wyraźna. Należy interpretować ich znaki zgodnie z konwencją stosowaną przez OghmaNano przy tworzeniu wykresów.Zdjęcia oznaczone jako R_n.csv i R_p.csv zawierają dodatkowe informacje na temat współczynników używanych w równaniach ciągłości elektronów i dziur. Pozwalają one zidentyfikować obszary, w których równowaga między lokalnymi składnikami generacji, wychwytu i rekombinacji jest najbardziej wyraźna. Należy interpretować ich znaki zgodnie z konwencją stosowaną przez OghmaNano przy tworzeniu wykresów.
Uwzględniając wszystkie te elementy – strukturę pasmową, gęstość nośników, prądy i profile zmian – można w pełni wyjaśnić charakterystykę diody: wbudowana bariera kontroluje wstrzykiwanie nośników, przyłożone napięcie modyfikuje wewnętrzną elektrostatykę, a transport i rekombinacja określają wynikowy, stały prąd.Uwzględniając wszystkie te elementy – strukturę pasmową, gęstość nośników, prądy i profile zmian – można w pełni wyjaśnić charakterystykę diody: wbudowana bariera kontroluje wstrzykiwanie nośników, przyłożone napięcie modyfikuje wewnętrzną elektrostatykę, a transport i rekombinacja określają wynikowy, stały prąd.
7. Ćwiczenia i interpretacja.7. Ćwiczenia i interpretacja.
8. Interpretacja modelu i dalsza analiza.8. Interpretacja modelu i dalsza analiza.
W tym przykładzie zastosowano jednowymiarowy, ustalony model dryfu i dyfuzji w arsenku indu (InP), prostą strukturę półprzewodnikową p/n/n⁺ oraz idealizowane styki omowe. Doskonale nadaje się on do zademonstrowania, w jaki sposób pojawia się prostowanie elektryczne na styku i jak to zjawisko wiąże się z przestrzennym rozkładem nośników ładunku. Bardziej rozbudowane badania urządzeń mogą uwzględniać zmierzone charakterystyki styków, udoskonalone parametry materiałowe, zmiany temperatury oraz dodatkowe szczegóły geometryczne.W tym przykładzie zastosowano jednowymiarowy, ustalony model dryfu i dyfuzji w arsenku indu (InP), prostą strukturę półprzewodnikową p/n/n⁺ oraz idealizowane styki omowe. Doskonale nadaje się on do zademonstrowania, w jaki sposób pojawia się prostowanie elektryczne na styku i jak to zjawisko wiąże się z przestrzennym rozkładem nośników ładunku. Bardziej rozbudowane badania urządzeń mogą uwzględniać zmierzone charakterystyki styków, udoskonalone parametry materiałowe, zmiany temperatury oraz dodatkowe szczegóły geometryczne.
Właściwości materiału InP uzyskiwane są z wybranego modelu Lua. W przypadku badań ilościowych należy sprawdzić wyznaczone wartości: szerokości pasma energetycznego, powinowactwa elektronowego, ruchliwości nośników i współczynników rekombinacji w edytorze materiałów lub w danych wyjściowych. Należy również sprawdzić zachowanie prądu i reszty rozwiązań, a następnie powtórzyć obliczenia dla reprezentatywnych punktów polaryzacji, stosując bardziej szczegółową siatkę, aby ustalić stabilność numeryczną.Właściwości materiału InP uzyskiwane są z wybranego modelu Lua. W przypadku badań ilościowych należy sprawdzić wyznaczone wartości: szerokości pasma energetycznego, powinowactwa elektronowego, ruchliwości nośników i współczynników rekombinacji w edytorze materiałów lub w danych wyjściowych. Należy również sprawdzić zachowanie prądu i reszty rozwiązań, a następnie powtórzyć obliczenia dla reprezentatywnych punktów polaryzacji, stosując bardziej szczegółową siatkę, aby ustalić stabilność numeryczną.
W dostarczonym przykładzie standardowy zakres zmian napięcia wynosi od 0 do 3 V. Kolejnym, pouczającym eksperymentem jest uwzględnienie niewielkiego zakresu napięcia wstecznego, co umożliwi bezpośrednie porównanie przewodzenia w kierunku prostym i prądu upływu w kierunku wstecznym. Fizycznie istotny zakres napięcia wstecznego i charakterystyka prądu upływu zależą od danego urządzenia oraz od mechanizmów uwzględnionych w symulacji.W dostarczonym przykładzie standardowy zakres zmian napięcia wynosi od 0 do 3 V. Kolejnym, pouczającym eksperymentem jest uwzględnienie niewielkiego zakresu napięcia wstecznego, co umożliwi bezpośrednie porównanie przewodzenia w kierunku prostym i prądu upływu w kierunku wstecznym. Fizycznie istotny zakres napięcia wstecznego i charakterystyka prądu upływu zależą od danego urządzenia oraz od mechanizmów uwzględnionych w symulacji.
Aby uzyskać więcej informacji, zmieniaj stężenie domieszek pojedynczo. Zwiększenie nominalnego stężenia domieszek zazwyczaj zmniejsza obszar zubożenia i zmienia równowagę elektrostatyczną; zmiana współczynników rekombinacji wpływa na prąd zależny od wstrzykiwania nośników. Są to przydatne, kontrolowane porównania, ponieważ materiały i geometria warstw mogą pozostać niezmienione.Aby uzyskać więcej informacji, zmieniaj stężenie domieszek pojedynczo. Zwiększenie nominalnego stężenia domieszek zazwyczaj zmniejsza obszar zubożenia i zmienia równowagę elektrostatyczną; zmiana współczynników rekombinacji wpływa na prąd zależny od wstrzykiwania nośników. Są to przydatne, kontrolowane porównania, ponieważ materiały i geometria warstw mogą pozostać niezmienione.
9. Kolejne kroki.9. Kolejne kroki.
Dioda InP PN stanowi dobry punkt wyjścia do tworzenia bardziej zaawansowanych struktur półprzewodnikowych na bazie InP. Naturalnym krokiem jest wprowadzenie słabo domieszkowanej warstwy i porównanie właściwości elektrostatycznych oraz parametrów transportowych złącza PN i PIN. Inną możliwością jest badanie heterostruktur zawierających InGaAs/InP, gdzie różnice w składzie powodują powstawanie dodatkowych przesunięć pasm energetycznych i zmieniają sposób wstrzykiwania nośników.Dioda InP PN stanowi dobry punkt wyjścia do tworzenia bardziej zaawansowanych struktur półprzewodnikowych na bazie InP. Naturalnym krokiem jest wprowadzenie słabo domieszkowanej warstwy i porównanie właściwości elektrostatycznych oraz parametrów transportowych złącza PN i PIN. Inną możliwością jest badanie heterostruktur zawierających InGaAs/InP, gdzie różnice w składzie powodują powstawanie dodatkowych przesunięć pasm energetycznych i zmieniają sposób wstrzykiwania nośników.