Samouczek dotyczący parabolicznej studni kwantowej GaAs/AlGaAs
1. Wprowadzenie
W tym samouczku obliczymy właściwości elektroniczne i optyczne prostej struktury GaAs/AlGaAs, wykorzystując model parabolicznych pasm w programie OghmaNano. Cienka warstwa GaAs znajduje się pomiędzy dwiema szerszymi barierami AlGaAs. Zmiana energii na krawędzi pasma ogranicza ruch elektronów i dziur w kierunku wzrostu, podczas gdy pozostają one swobodne w ruchu równoległym do warstw. Powoduje to powstawanie dyskretnych energii ograniczenia, które przekształcają się w dwuwymiarowe podpasma elektroniczne.
Zaczniemy od gotowego modelu symulacyjnego, przeanalizujemy strukturę warstw, uruchomimy model i przyjrzymy się obliczonym pasmom przewodnictwa i walencyjnym. Następnie zbadamy, w jaki sposób te stany wpływają na absorpcję, emisję wymuszoną, całkowite wzmocnienie optyczne i emisję spontaniczną. Jeśli dopiero zaczynasz przygodę z fizyką, strona parabolic-band theory page zawiera bardziej szczegółowe wprowadzenie do tematu.
2. Otwieranie przykładowego pliku dotyczącego studni kwantowej.
Otwórz OghmaNano i wybierz opcję „Nowa symulacja” z menu „Plik”. W oknie wyboru urządzenia kliknij dwukrotnie przykładowe pliki GaAs (??). Spowoduje to otwarcie listy dostępnych przykładów dotyczących półprzewodników i struktur kwantowych. Kliknij dwukrotnie przykładowy plik GaAs/AlGaAs QW parabolic (??). W głównym oknie symulacji otworzy się wcześniej skonfigurowany przykład.
3. Poznanie budowy urządzenia.
Główne okno wyświetla strukturę składającą się z trzech warstw (??). Warstwa GaAs stanowi studnię, a dwie otaczające ją warstwy AlGaAs pełnią funkcję barier. Aby sprawdzić wymiary, kliknij przycisk „Edytor warstw” znajdujący się po lewej stronie głównego okna. W dostarczonym przykładzie zastosowano następującą strukturę:
| Warstwa | Grubość | Cel |
|---|---|---|
| AlGaAs (stop galowo-arsenkowy) | 40 nm | Pierwsza przeszkoda |
| GaAs (struktura wielowarstwowa) | 10 nanometrów | Studnia kwantowa |
| AlGaAs (stop galowo-arsenkowy) | 40 nm | Druga bariera |
Zwróć uwagę, że środkowej warstwie przypisano w sposób wyraźny typ warstwy „Quantum well”, a nie „Active”. Dostępne opcje i wybrane materiały optyczne można zobaczyć w ??. Warstwy o większej przerwie energetycznej tworzą potencjały ograniczające; dokładna wartość przesunięcia pasm zależy od parametrów materiałowych użytych w przykładzie.
Dlaczego zmiana szerokości odwiertu ma znaczenie?
Ograniczenie wzdłuż kierunku wzrostu pozwala na występowanie tylko określonych funkcji fal elektronów i dziur. W najprostszym, idealnym modelu jednowymiarowego nieskończonego potencjału o szerokości \(L\), energie ograniczenia elektronów, mierzone od dna potencjału, wynoszą:
\[ E_n=\frac{n^2\pi^2\hbar^2}{2m^{*}L^2},\qquad n=1,2,3,\ldots \]
W tym wzorze \(m^{*}\) oznacza efektywną masę nośników, a \(\hbar\) – zredukowaną stałą Plancka. W rzeczywistych przykładach bariery AlGaAs mają skończoną grubość, dlatego OghmaNano rozwiązuje równanie dla stanów skwantowanych, wykorzystując rzeczywisty profil materiału, zamiast zakładać model nieskończonej studni potencjału. Niemniej jednak, równanie to pozwala na uzyskanie użytecznych przewidywań: zmniejszenie szerokości warstwy GaAs zazwyczaj zwiększa energię skwantowania elektronów. Dziury mają własną efektywną masę i odpowiadające jej energie skwantowania.
Można przeglądać lub modyfikować podstawowe funkcje materiału za pomocą przycisku „Parametry elektryczne”. Na zrzucie ekranu widać, że definicja materiału AlGaAs zawiera obliczenia dotyczące szerokości pasma, zależne od składu i temperatury (??). Aby ukończyć ten wstępny samouczek, nie trzeba modyfikować kodu materiału.
4. Uruchamianie symulacji.
Wróć do zakładki „Plik” i naciśnij przycisk „Uruchom symulację” (niebieski trójkątny przycisk). Na karcie „Terminal” wyświetlane są komunikaty generowane przez program podczas działania. W podanym przykładzie w terminalu wyświetlane są wartości energii stanów własnych oraz powtarzające się komunikaty block converged after ... iterations (??). Poczekaj, aż symulacja się zakończy, zanim otworzysz wykresy.
Rozwiązywacz oblicza skwantowane stany w poszczególnych warstwach, a następnie analizuje, jak zmieniają się one w zależności od wektora fal w płaszczyźnie, \(k_{\parallel}\). Jeśli skonfigurowano obliczenia optyczne, generuje on również widma dla różnych warunków nośnika w danym przykładzie.
5. Elektroniczne podpasma: interpretacja wykresu pasmowego.
Otwórz wykres rozkładu energii w pasmach (??). Widoczne są dwie wyraźnie odrębne grupy krzywych: pasma przewodnictwa odpowiadające dodatnim wartościom energii oraz pasma walencyjne znajdujące się w pobliżu i poniżej poziomu odniesienia pasma walencyjnego. Pionowa odległość między nimi odzwierciedla szerokość przerwy energetycznej półprzewodnika oraz energię elektronów i dziur w obszarze potencjału. Różne wartości energii początkowej w obrębie każdej grupy wynikają z skwantowanych stanów w potencjalnej studni.
Choć ruch wewnątrz studni jest ograniczony, nośniki mogą nadal przemieszczać się wzdłuż warstw. Dlatego pojedyncza, zlokalizowana energia, po przedstawieniu jej w formie wykresu w odniesieniu do równoległego wektora falowego, przyjmuje kształt paraboli:
\[ E_{e,n}(k_{\parallel})=E_{e,n}(0) +\frac{\hbar^2 k_{\parallel}^2}{2m_e^{*}}. \]
W tym miejscu \(E_{e,n}(0)\) przedstawiono energię na dolnej granicy pasma elektronowego \(n\). W przypadku odpowiadających mu pasm elektronów walencyjnych, prosta, paraboliczna krzywa wygina się w dół:
\[ E_{v,m}(k_{\parallel})=E_{v,m}(0) -\frac{\hbar^2 k_{\parallel}^2}{2m_h^{*}}. \]
To jest kluczowe przybliżenie w tym samouczku: w kierunku ograniczenia stany mają skwantowane energie, ale w płaszczyźnie zachowują dyspersję związaną z efektywną masą. Bardziej zaawansowane wielopasmowe modele k·p pozwalają na mieszanie się pasm i sprawiają, że stają się one nieparaboliczne; strona 8-band zincblende theory page wyjaśnia, dlaczego ma to znaczenie.
6. Od pasm częstotliwości do widm optycznych.
Elektrony mogą przechodzić między pasmami walencyjnymi i przewodnictwa w wyniku oddziaływania ze światłem. Przy normalnym padaniu promieni światła, pęd fotonu jest mały w porównaniu z wektorem falowym elektronu, dlatego prosty model uwzględnia tylko przejścia, w których zachowywany jest wektor falowy w płaszczyźnie. W centrum rozproszenia, przybliżona energia fotonu dla przejścia między zlokalizowanym poziomem elektronowym \(n\) a poziomem dziury \(m\) wynosi:
\[ \hbar\omega\simeq E_g+E_{e,n}^{\mathrm{conf}}+E_{h,m}^{\mathrm{conf}}. \]
Obie energie związane z ograniczeniem są mierzone w górę od odpowiednich krawędzi pasm energetycznych. Dostępne przejścia, zajętość stanów początkowych i końcowych oraz wszelkie poszerzenia widmowe decydują o kształcie obliczonej odpowiedzi optycznej. Można dokładniej zapoznać się z tymi koncepcjami na stronie poświęconej teorii przejść optycznych strona o teorii przejść optycznych.
6.1 Emisja wymuszona i całkowite wzmocnienie
Wykres przedstawiający emisję wymuszoną (??) pokazuje obliczony udział emisji w zależności od długości fali dla serii warunków nośnych zastosowanych w przykładzie. Wyraźne skoki i ramiona na spektrum odpowiadają energiom, przy których zachodzą różne przejścia między podpasmami. Oddzielny, skumulowany współczynnik wzmocnienia (??) uwzględnia konkurencyjny udział absorpcji wymuszonej i dlatego może być dodatni lub ujemny.
Schematycznie, dla stanu przewodnictwa, w którym prawdopodobieństwo wystąpienia wynosi \(f_c\), oraz dla stanu elektronu walencyjnego, w którym prawdopodobieństwo wystąpienia wynosi \(f_v\), całkowity wkład wynikający ze stymulacji zawiera współczynnik zajętości.
\[ \underbrace{f_c(1-f_v)}_{\text{stimulated emission}} -\underbrace{f_v(1-f_c)}_{\text{absorption}} =f_c-f_v. \]
Dodatni wkład netto w proces stymulowanej emisji wymaga odpowiedniego odwrócenia stanów, które biorą udział w przejściu. Całkowity rejestrowany widm wymaga dodatkowo uwzględnienia sił przejściowych oraz liczby dostępnych stanów.
6.2 Absorpcja i spontaniczna emisja.
Wykres absorpcji (??) rejestruje udział procesów absorpcji w dozwolonych przejściach. Absorpcja wymaga zajętego początkowego stanu walencyjnego i dostępnego końcowego stanu przewodnictwa. Wykres emisji spontanicznej (??) przedstawia radiacyjną rekombinację zajętych stanów przewodnictwa w puste stany walencyjne, bez konieczności użycia padającego fotonu do zainicjowania przejścia. W związku z tym, współczynnik zajętości emisji spontanicznej zawiera \(f_c(1-f_v)\).
Zauważ, że wyraźne cechy widma występują w podobnych zakresach długości fal na wszystkich wykresach. Jest to zgodne z oczekiwaniami: ta sama podstawowa struktura podpasm elektronów i dziur określa możliwe energie przejść optycznych, natomiast rozkład nośników ładunku odróżnia absorpcję od różnych rodzajów emisji. Na wykresie przedstawiającym emisję spontaniczną zastosowano logarytmiczną skalę pionową, dlatego należy zachować ostrożność podczas porównywania wysokości krzywych.
7. Lokalizowanie wygenerowanych plików.
Po zakończeniu symulacji otwórz zakładkę „Wyniki” w głównym oknie (??). W przykładzie generowane są oddzielne pliki CSV zawierające wyniki różnych obliczeń spektralnych. Widoczne nazwy plików zawierają wpisy zgodne ze wzorcem gain_00.csv, z dodatkiem przyrostków wskazujących na absorpcję, zysk netto, emisję spontaniczną i wymuszoną, a następnie pliki o podobnych nazwach dla kolejnych obliczonych warunków. Kliknij dwukrotnie wynik, aby wyświetlić go w oknie graficznym. Ponieważ wyniki są zapisywane jako pliki CSV, można je również analizować lub ponownie wyświetlać w programach zewnętrznych, poza OghmaNano.
Ponumerowane pliki pozwalają śledzić zmiany obliczonego widma w różnych warunkach pracy urządzenia, a nie tylko analizować pojedyncze obliczenia optyczne.
8. Wypróbuj to sam: zmień szerokość odwiertu.
Proponowane doświadczenie: Wróć do edytora warstw i zmień grubość warstwy GaAs (QW) z 10 nm na 8 nm, pozostawiając warstwy barierowe AlGaAs bez zmian. Uruchom ponownie symulację i porównaj uzyskane energie podpasm oraz widma optyczne z pierwotnym przykładem.
- Czy najniższy poziom ograniczenia elektronów przesuwa się w górę w stosunku do krawędzi pasma przewodnictwa GaAs?
- Jak zmieniają się energie przejścia elektronów i dziur?
- Czy początek widma optycznego przesuwa się w kierunku krótszych czy dłuższych długości fal?
✅ Oczekiwany trend w zakresie parametrów fizycznych
Zazwyczaj, im węższa studnia potencjału, tym większa energia ograniczenia. W związku z tym najniższa energia przejścia między pasmami ma tendencję do przesuwania się w kierunku wyższej energii fotonu, co odpowiada krótszej długości fali, chociaż dokładne przesunięcie zależy od wysokości barier i efektywnych mas elektronów i dziur. Jak wykazano w \(\lambda=hc/E_{\mathrm{photon}}\), wyższa energia przejścia oznacza mniejszą długość fali emisji lub absorpcji.
Gdy już dobrze zapoznasz się z tym przykładem, naturalnym następnym krokiem będzie porównanie go z ośmiopasmową wersją modelu k·p dla struktury GaAs/AlGaAs, opisującej kwantową studnię. Zachowuje ona efekt kwantowego uwięzienia, ale wykracza poza niezależne paraboliczne pasma, aby opisać oddziaływania i mieszanie się stanów przewodnictwa i walencyjności.