GaInNAs – półprzewodnik z niską zawartością azotka: kwantowa studnia o strukturze cynkowej – przewodnik po metodzie k·p z wykorzystaniem 10 pasm energetycznych.
1. Wprowadzenie
W tym samouczku zasymulujemy strukturę GaInNAs, czyli warstwę półprzewodnikową o niskiej zawartości azotu, wykorzystując dziesięciopasmowy model k·p dla struktury cynkowej, dostępny w programie OghmaNano. Otworzymy istniejący przykład, przeanalizujemy strukturę warstw, uruchomimy obliczenia dotyczące stanów elektronowych i zbadamy przewidywane widma optyczne. W przykładzie znajduje się warstwa GaInNAs o grubości 7 nm, umieszczona w większej strukturze GaAs/AlGaAs, a nie jako izolowana, trz warstwowa struktura.
Dlaczego przydatny jest model dziesięciopasmowy? Zwykła ośmiopasmowa teoria k·p opisuje pasmo przewodnictwa oraz pasma walencyjne ciężkich dziur, lekkich dziur i rozszczepionych dziur, uwzględniając spin. Wprowadzenie niewielkiej koncentracji azotu do stopu III–V może stworzyć rezonansowy poziom związany z azotem, który szczególnie silnie oddziałuje z pasmem przewodnictwa materiału bazowego. Opis dziesięciopasmowy dodaje dwie składowe spinowe tego stanu związanego z azotem. To oddziaływanie może znacząco obniżyć podstawową energię przejścia i zmienić krzywiznę pasma przewodnictwa. Dzięki temu rozcieńczone azotki stają się interesujące w kontekście urządzeń optycznych działających w zakresie dłuższych fal.
2. Otwieranie przykładowego pliku z dziesięcioma pasmami.
Wybierz Plik → Nowa symulacja. W oknie nowej symulacji kliknij dwukrotnie GaAs demos (??). W następnym oknie znajdują się różne przykłady struktur kwantowych. Wybierz AlGaAs, rozcieńczony azotek, 10 pasm k.p (??). Spowoduje to załadowanie wstępnie skonfigurowanej symulacji z dziesięcioma pasmami, w tym odpowiednich definicji warstw i plików z parametrami materiałowymi.
3. Budowa urządzenia i parametry materiałowe.
Główne okno wyświetla wielowarstwową strukturę (??). W przeciwieństwie do prostych struktur GaAs/AlGaAs, stosowanych w poprzednich samouczkach, w tym przykładzie zastosowano dodatkowe warstwy ochronne, warstwy o stopniowo zmieniającej się grubości oraz warstwy barierowe, otaczające warstwę azotkowego półprzewodnika. Aby sprawdzić geometrię, kliknij opcję „Edytor warstw” w lewym pasku narzędzi.
Edytor warstw (??) wyświetla osiem wpisów, zaczynając od p-GaAs, przez AlGaAs i warstwę aktywną, aż do n-GaAs. W szczególności, centralna warstwa GaInNAs (QW) ma grubość 7,0 nm, a jej typ warstwy jest ustawiony na „Quantum well” (studnia kwantowa). Widoczny przykład zawiera również warstwy AlGaAs o grubości 0,16 µm po obu stronach warstwy aktywnej oraz warstwę barierową/SCL o grubości 0,02 µm. Korzystając z dostarczonego szablonu, nie trzeba ręcznie odtwarzać struktury warstw.
3.1 Sprawdzanie definicji materiałów.
Kliknij „Parametry elektryczne”, wybierz zakładkę „GaInNAs (QW)” i otwórz edytor kodu (??). Zrzut ekranu ilustruje, w jaki sposób plik materiałowy zawiera wielkości fizyczne, takie jak efektywna gęstość stanów pasma walencyjnego i mobilność elektryczna. Na przykład, obserwowana mobilność elektronów przy niskim polu wynosi 0.05 m² V⁻¹ s⁻¹, co odpowiada 500 cm² V⁻¹ s⁻¹. Są to dane wejściowe dla modelu materiałowego, a nie wartości wyliczone na podstawie wykresu dyspersji studni kwantowej.
W obliczeniach dotyczących roztworów azotków szczególnie ważne jest, aby skład, energie na krawędzi pasma i parametry oddziaływania azotu były wzajemnie spójne. W tym pierwszym etapie pozostaw wartości domyślne bez zmian; omówimy przydatne eksperymenty dotyczące składu na końcu strony.
4. Jakie zmiany zachodzą, gdy wprowadzamy azot?
Zanim uruchomimy model, warto zrozumieć, dlaczego hamiltonian ma dziesięć pasm, a nie osiem. Wyobraźmy sobie, że krawędź pasma przewodnictwa półprzewodnika bazowego przesuwa się w górę w skali energii wraz ze wzrostem wektora fal elektronu w płaszczyźnie. Stan rezonansowy związany z azotem ma swoją charakterystyczną energię. Gdyby oba stany nie oddziaływały ze sobą, ich energie mogłyby się zbliżać i przecinać. Jednak w rozcieńczonym azotku azot oddziałuje ze stanami przewodnictwa półprzewodnika bazowego, co zmienia wynikowe poziomy energii.
Dobrym przykładem ilustrującym tę zasadę fizyki jest macierz dwustanowego przecięcia się pasm (BAC):
\[ H_{\mathrm{BAC}}(k)= \begin{pmatrix} E_C(k) & V_N\\ V_N & E_N \end{pmatrix}, \qquad V_N=\beta_N\sqrt{x_N}. \]
W tym przypadku \(E_C(k)\) oznacza rozkład pasma przewodnictwa, \(E_N\) – energię poziomu związanego z azotem, \(x_N\) – udział azotu, a \(\beta_N\) kontroluje siłę sprzężenia. Diagonalizacja tej małej macierzy daje dwie gałęzie:
\[ E_{\pm}(k)=\frac{E_C(k)+E_N}{2} \pm\frac{1}{2}\sqrt{\big[E_C(k)-E_N\big]^2+4V_N^2}. \]
Wyraz zawierający pierwiastek kwadratowy zapobiega przecięciu się rozgałęzień, gdy \(V_N\neq0\). W szczególności, dolna gałąź o charakterze przewodzącym może przemieszczać się w dół w stosunku do pierwotnej gałęzi przewodzącej, co zmniejsza podstawową przerwę energetyczną. Ponieważ oddziaływanie zależy od odległości między \(E_C(k)\) a \(E_N\), dyspersja może również przybrać silnie nieparaboliczny charakter.
Ten prosty przykład (macierz 2x2) ilustruje podstawową ideę; nie przedstawia on jednak pełnego rozwiązania hamiltonianu dla studni kwantowej, które jest omawiane w tym samouczku. Dziesięciopasmowy model k·p zachowuje szerszą bazę przewodnictwa/walencji charakterystyczną dla struktury cynkowej, uwzględniając spin, i dodaje komponenty związane z azotem. Następnie uwzględnia się kwantowe ograniczenie wzdłuż kierunku wzrostu. Pełne opracowanie teoretyczne można znaleźć w rozdziale poświęconym dziesięciopasmowemu modelowi struktury cynkowej.
5. Uruchamianie symulacji.
Wróć do zakładki „Plik” i naciśnij „Uruchom symulację”. Wybierz zakładkę „Terminal”, aby śledzić obliczenia (??). Program oblicza rzadką hamiltonianę, określa ograniczone stany własne przewodnictwa i walencyjności, a następnie powtarza obliczenia dla sekwencji wektorów fal w płaszczyźnie \(k_{\parallel}\).
Na zrzucie ekranu widać, że w dzienniku zidentyfikowano solver KP10, wyświetlono obliczone poziomy energii i zgłoszono, że proces zbiegł się po określonej liczbie iteracji. Pary blisko położonych poziomów energii mogą obejmować degeneracje związane ze spinem; nie należy ich automatycznie interpretować jako odrębnych gałęzi antyprzecięcia azotu. Poczekaj na zakończenie obliczeń, zanim przejdziesz do analizy plików wynikowych.
6. Zrozumienie obliczonej struktury pasmowej.
Otrzymane wykresy przedstawiają zależność energii od wektora falowego w płaszczyźnie, a nie energii od głębokości wewnątrz urządzenia. Potencjał studni kwantowej ogranicza nośniki ładunku wzdłuż kierunku wzrostu warstwy, tworząc dyskretne stany; ruch równoległy do warstw powoduje, że każdy z ograniczonych stanów ma określoną dyspersję \(E_n(k_{\parallel})\). W prostym parabolicznym modelu ta dyspersja jest kwadratowa. Zastosowanie sprzężonego hamiltonianu dziesięciu pasm pozwala na odchylenie od paraboli.
6.1 Stany przewodnictwa i zmniejszona przerwa energetyczna
Wykres przewodnictwa (??) przedstawia dwie obliczone gałęzie, rozpoczynające się przy różnych poziomach energii. Najniższy stan przewodnictwa zaczyna się przy około 0,89 eV w wybranym punkcie odniesienia dla energii w symulacji. Nie należy mylić tej wartości z przerwą energetyczną materiału GaInNAs: jest to energia własna dla danej, ograniczonej struktury. Zauważ, że dwie przedstawione krzywe nie mają identycznego przebiegu. Jest to jeden z elementów, które należy zbadać podczas porównywania obliczeń z wykorzystaniem ośmiu pasm z modelem związków azotkowych.
Uproszczone powyżej równania BAC pomagają wyjaśnić, dlaczego rozproszenie przewodnictwa może się zmieniać w materiale zawierającym azot. Same w sobie nie dowodzą jednak, jaki jest udział azotu w kształtowaniu krzywych o określonym kolorze. Aby to ustalić, należy przeanalizować odpowiednie wektory własne lub zmieniać wartość \(x_N\) i obserwować, jak zmieniają się stany.
6.2 Stany walencyjne
Diagram zależności energii walencyjnej (??) zawiera kilka gałęzi o wyraźnie nieparabolicznym kształcie. Charakter ciężkich i lekkich dziur może się mieszać przy skończonej wartości wektora fal w płaszczyźnie; na kształt krzywych wpływają energie związane z ograniczeniem, oddziaływanie spinowo-orbitalne oraz sprzężenie między stanami bazowymi. Wznoszący się odcinek lub zmiana krzywizny niekoniecznie wskazują na błąd numeryczny. Z kolei gładko narysowana krzywa sama w sobie nie jest dowodem na to, że parametry materiału są poprawne.
Najlepiej traktować poszczególne kolory jako oddzielne gałęzie stanów własnych, zamiast przypisywać im na stałe oznaczenia HH lub LH w całym zakresie wektora fal, bez sprawdzania ich charakteru.
6.3 Analiza obu zespołów łącznie
Włączenie skumulowanego rozpraszania (??) powoduje umieszczenie stanów elektronowych i walencyjnych na tej samej skali energii. Podstawowe przejście między pasmami jest związane z ich rozdzieleniem, a nie z absolutną energią któregokolwiek z przedstawionych stanów:
\[ \hbar\omega_{nm}(k_{\parallel}) =E_{c,n}(k_{\parallel})-E_{v,m}(k_{\parallel}). \]
To wyrażenie opisuje bezpośrednią, w przybliżeniu pionową, optyczną zmianę stanu w płaszczyźnie wektora falowego: pęd fotonu jest mały w porównaniu z pędem elektronów w sieci krystalicznej. Obserwowany widm jest jednak zazwyczaj złożony i nie ogranicza się do jednej ostrej linii. Wiele dozwolonych przejść przyczynia się do powstania widma, a ich intensywność i prawdopodobieństwo wystąpienia różnią się w zależności od stanu.
7. Od stanów elektronicznych do widm optycznych.
Dzięki możliwości definiowania określonych stanów, OghmaNano może generować widma optyczne. Kluczowe pytanie z zakresu fizyki nie dotyczy jedynie tego, czy elektron i dziura mają odpowiednią różnicę energii, ale także tego, jak silnie oba stany oddziałują z polem elektromagnetycznym oraz czy ich stany początkowe i końcowe są odpowiednio wypełnione.
7.1 Dlaczego niektóre przejścia mają większy wpływ niż inne?
Zasada Fermiego stanowi użyteczny punkt wyjścia do zrozumienia szybkości przejść. W przypadku przejścia z początkowego stanu elektronowego \(i\) do stanu końcowego \(f\), wywołanego słabym zakłóceniem optycznym, jej schematyczny zapis ma postać:
\[ W_{i\rightarrow f}=\frac{2\pi}{\hbar} \left|\langle f|H_{\mathrm{int}}|i\rangle\right|^2 \delta(E_f-E_i-\hbar\omega). \]
Element macierzy \(\langle f|H_{\mathrm{int}}|i\rangle\) opisuje oddziaływanie stanów z światłem. Zależy on między innymi od ich charakteru orbitalnego, funkcji obwiedni i polaryzacji światła. Funkcja delta zapewnia zachowanie energii. W numerycznym widmie idealizowane, ostre przejście jest reprezentowane za pomocą odpowiedniej funkcji poszerzającej linie. Po zsumowaniu wartości dla dozwolonych podpasm i wektorów fal w płaszczyźnie, wynikiem jest widmo, a nie pojedyncza energia przejścia.
W przypadku nośników ładunku rozróżnia się absorpcję i emisję. Jeśli \(f_c\) oznacza zajętość stanu przewodnictwa, a \(f_v\) zajętość stanu elektronu walencyjnego, charakterystyczne czynniki to:
\[ \begin{aligned} \text{absorption:}\quad &(1-f_c)f_v,\\ \text{stimulated emission:}\quad &f_c(1-f_v),\\ \text{net stimulated contribution:}\quad &f_c-f_v. \end{aligned} \]
Ostatnie wyrażenie to różnica między procesami optycznymi zachodzącymi w kierunku wzrostu i spadku. Dodatni wkład w wzmocnienie wymaga wystarczającej inwersji populacji. Czynniki te wyjaśniają interpretację przedstawionych poniżej wykresów; więcej szczegółów można znaleźć w sekcji teoria przejść optycznych.
7.2 Emisja i absorpcja wymuszona
Na dwóch wykresach (?? i ??) przedstawiono konkurujące ze sobą procesy optyczne zachodzące w tej samej strukturze pasmowej studni kwantowej. Liczne krzywe ilustrują sekwencję warunków zastosowanych w przykładzie; ich względne wartości zmieniają się, ponieważ to poziom wypełnienia nośników decyduje o tym, które przejścia są możliwe.
Widma obejmują zakres bliskiej podczerwieni, mniej więcej od 1,1 do 1,5 µm. Dla orientacji, energia fotonu i długość fali w próżni są powiązane następującym wzorem:
\[ E_{\mathrm{photon}}=\frac{hc}{\lambda} \simeq \frac{1240\ \mathrm{eV\,nm}}{\lambda\,[\mathrm{nm}]}. \]
Zatem 1300 nm odpowiada w przybliżeniu 0,954 eV. Pozwala to na sprawdzenie spójności między różnicą energii w połączonym wykresie pasm a skalą widm optycznych, przy uwzględnieniu efektów kwantowych, zakresu uwzględnionych przejść oraz wybranej wartości odniesienia energii.
7.3 Zysk netto: występuje, gdy emisja przewyższa absorpcję.
W przypadku ośrodka wzmacniającego istotna jest równowaga między emisją wymuszoną a absorpcją. Przy zastosowaniu przyjętej konwencji znaków, można to przedstawić schematycznie w następujący sposób:
\[ g_{\mathrm{net}}(\lambda) =g_{\mathrm{stim}}(\lambda)-\alpha(\lambda). \]
Na wykresie przedstawiającym obliczoną wartość zysku netto (??) znajdują się zarówno krzywe o wartościach ujemnych, jak i dodatnich. Gdy absorpcja jest dominująca, \(g_{\mathrm{net}}\) ma wartość ujemną. Wraz ze zmianą koncentracji elektronów i dziur, stymulowana emisja może stać się na tyle silna, że w pewnych zakresach widma pojawi się dodatkowy zysk materiału. Przedstawiony na wykresie obszar przejścia charakteryzuje się wyraźną strukturą, a nie jest to idealnie gładki pojedynczy pik, ponieważ w procesie uczestniczy wiele przejść między stanami o ograniczonych energiach.
Ważne: sam fakt uzyskania dodatkowych korzyści materialnych nie jest wystarczającym wskaźnikiem, że w pełni sprawny laser będzie działał powyżej progu. Należy również uwzględnić straty w rezonatorze, ograniczenie optyczne, współczynnik odbicia lusterek oraz wstrzykiwanie nośników.
8. Lokalizowanie wyników symulacji.
bands i gain.Po zakończeniu obliczeń otwórz zakładkę „Wyniki” (??). Na zrzucie ekranu widać główne pliki wynikowe, w tym: Ec.csv (gałęzie przewodnictwa), Ev.csv (gałęzie walencyjne), Ec_Ev.csv (połączone pasma), gain.csv, absorption.csv, net_gain.csv i spontaneous.csv. Dodatkowe wyniki znajdują się w katalogach bands i gain. Kliknij dwukrotnie plik CSV, aby wyświetlić jego zawartość w postaci wykresu, lub wykorzystaj zapisane dane liczbowe do dalszej analizy.
Chociaż w tym samouczku przedstawiono procesy emisji wymuszonej, absorpcji i wzmocnienia netto, przykład ilustruje również emisję spontaniczną. Emisja spontaniczna to proces radiacyjnej rekombinacji, który nie wymaga obecności przychodzącego fotonu, a jej charakterystyczny współczynnik zajętości wynosi \(f_c(1-f_v)\). Jej bezwzględne współczynniki spektralne różnią się od współczynników emisji wymuszonej.
9. Wypróbuj to samodzielnie: przeprowadź badania dotyczące oddziaływania azotu.
Proponowane ćwiczenie: przed rozpoczęciem edycji zapisz oddzielną kopię tego przykładowego pliku. Następnie przeanalizuj stężenie azotu w definicji materiału GaInNAs, zmieniając jeden parametr na raz, przy jednoczesnym zachowaniu stałej geometrii warstwy o grubości 7 nm.
- Zapisz najniższą wartość energii własnej elektronu w punkcie \(k_{\parallel}=0\).
- Zmień wartość udziału azotu, tak aby mieściła się w odpowiednim, fizycznie uzasadnionym zakresie dla wybranego modelu.
- Ponownie oblicz rozproszenie przewodnictwa i porównaj kształt najniższej gałęzi.
- Porównaj położenie krawędzi absorpcji z głównymi cechami wzmocnienia.
Ilustracja BAC przewiduje, że zwiększenie oddziaływania azotu może prowadzić do większego rozdzielenia połączonych gałęzi i zmiany dolnej energii przewodnictwa. Jednak w przypadku pełnych obliczeń dotyczących stopu, skład może wpływać również na inne parametry materiału. Dlatego należy interpretować wszelkie wynikowe zmiany, korzystając z rzeczywistych funkcji materiałowych dostępnych w symulacji, zamiast zakładać, że proste dwustanowe wzory stanowią pełny model stopu.
✅ Na co zwrócić uwagę
Zbadaj, czy zmienia się położenie krawędzi pasma przewodnictwa i czy zmienia się jej krzywizna w zależności od stężenia azotu. Porównaj widma optyczne przy identycznych ustawieniach wypełnienia, aby zmiany w energiach przejść nie były mylone z niezwiązanymi zmianami w populacji nośników. Jeśli masz możliwość analizy składu wektorów własnych, szczególnie pouczające jest śledzenie zmian wkładu pasma przewodnictwa gospodarza oraz wkładu związanego z azotem w stanach.
W drugim eksperymencie należy utrzymać stałe parametry materiału i zmienić szerokość warstwy kwantowej z 7 nm na 6 nm. Zwiększenie stopnia kwantowego uwięzienia zazwyczaj powoduje podniesienie najniższego poziomu energii elektronu i zmienia dozwolone energie przejść między pasmami. Wielkość przesunięcia jest określana przez złożone obliczenia uwzględniające dziesięć pasm, a nie przez proste przypisanie jednej stałej efektywnej masy do każdego stanu.