Einführung in parabolische Quantenpotentialtöpfe aus GaAs/AlGaAs
1. Einleitung
In diesem Tutorial werden wir die elektronischen und optischen Eigenschaften eines einfachen GaAs/AlGaAs-Quantenraums mithilfe des parabolischen Bandmodells in OghmaNano berechnen. Eine dünne GaAs-Schicht wird zwischen zwei breiteren AlGaAs-Barrieren angeordnet. Die Änderung der Bandkantenenergie begrenzt die Bewegung von Elektronen und Löchern in Wachstumsrichtung, während sie sich weiterhin frei parallel zu den Schichten bewegen können. Dadurch entstehen diskrete Einschlussenergien, die sich zu zweidimensionalen elektronischen Unterschichten entwickeln.
Wir werden mit einer vorgefertigten Simulation beginnen, die Schichtstruktur untersuchen, das Modell ausführen und die berechneten Leitungs- und Valenzunterbänder betrachten. Anschließend werden wir untersuchen, wie diese Zustände zur Absorption, zur stimulierten Emission, zur Netto-optischen Verstärkung und zur spontanen Emission beitragen. Wenn Sie mit der Physik noch nicht vertraut sind, bietet die Seite parabolische-Band-Theorie eine detailliertere Einführung.
2. Öffnen des Beispiels für eine Quantenpotentialstruktur
Öffnen Sie OghmaNano und wählen Sie im Menü „Datei“ die Option „Neue Simulation“ aus. Im Geräteauswahlfenster doppelklicken Sie auf „GaAs-Beispiele“ (??). Dadurch wird die Auswahl an Halbleiter- und Quantenwellenbeispielen geöffnet. Doppelklicken Sie auf „GaAs/AlGaAs-Quantenwelle, parabolisch“ (??). Das voreingestellte Beispiel wird im Hauptfenster der Simulation geöffnet.
3. Die Struktur des Geräts verstehen
Das Hauptfenster zeigt die dreischichtige Struktur (??). Die GaAs-Schicht bildet den aktiven Bereich; die beiden umgebenden AlGaAs-Schichten bilden die Barrieren. Klicken Sie auf „Schichteditor“ auf der linken Seite des Hauptfensters, um die Abmessungen zu überprüfen. Das mitgelieferte Beispiel verwendet die folgende Struktur:
| Schicht | Dicke | Zweck |
|---|---|---|
| Aluminiumgalliumarsenid | 40 Nanometer | Erste Hürde |
| GaAs (Quantenbrunnen) | 10 Seemeilen | Quantenpotentialtopf |
| Aluminiumgalliumarsenid | 40 Nanometer | Zweite Hürde |
Beachten Sie, dass der mittleren Schicht explizit der Quantenwellenschicht-Typ zugewiesen wird und nicht der aktive Schicht-Typ. Die Einträge und die Auswahl der optischen Materialien sind in ?? sichtbar. Die breiterbandigen Barrieren erzeugen die Einschluss-Potenziale; die genauen Bandversätze hängen von den in dem Beispiel verwendeten Materialparametern ab.
Warum ist es wichtig, die Breite des Bohrlochs zu verändern?
Die räumliche Begrenzung entlang der Wachstumsrichtung erlaubt nur bestimmte Elektronen- und Loch-Wellenfunktionen. In der einfachsten idealisierten Darstellung eines eindimensionalen, unendlich tiefen Potentials mit der Breite \(L\) sind die Elektronen-Einkapselungsenergien, gemessen vom tiefsten Punkt des Potentials,
\[ E_n=\frac{n^2\pi^2\hbar^2}{2m^{*}L^2},\qquad n=1,2,3,\ldots \]
Hierbei ist \(m^{*}\) die effektive Masse der Ladungsträger und \(\hbar\) die reduzierte Planck-Konstante. In der realen Anwendung gibt es endliche AlGaAs-Barrieren, daher berechnet OghmaNano die eingeschränkten Zustände unter Verwendung des tatsächlichen Materialprofils und nicht unter der Annahme eines unendlichen Potentialtopfs. Dennoch liefert die Gleichung eine nützliche Vorhersage: Eine Verringerung der GaAs-Schichtdicke führt im Allgemeinen zu einer Erhöhung der Energie der eingeschränkten Elektronen. Löcher haben ihre eigene effektive Masse und die entsprechenden Einschlussenergien.
Sie können die zugrunde liegenden Materialfunktionen über die Schaltfläche „Elektrische Parameter“ untersuchen oder bearbeiten. In dem Screenshot ist zu sehen, dass die Definition des AlGaAs-Materials eine berechnete Bandlücke enthält, die von der Zusammensetzung und der Temperatur abhängt (??). Um dieses Einführungstutorial abzuschließen, müssen Sie den Materialcode nicht bearbeiten.
4. Durchführung der Simulation
Kehren Sie zur Registerkarte „Datei“ zurück und klicken Sie auf „Simulation starten“ (der blaue, dreieckige Button). In der Registerkarte „Terminal“ werden die Meldungen des Solvers während der Ausführung angezeigt. Im bereitgestellten Beispiel werden im Terminal die Eigenzustandsenergien und wiederholte block converged after ... iterations-Meldungen (??) angezeigt. Warten Sie, bis die Ausführung abgeschlossen ist, bevor Sie die Diagramme öffnen.
Der Solver berechnet die quantisierten Zustände über die verschiedenen Schichten und untersucht anschließend deren Veränderung in Abhängigkeit vom Wellenvektor in der Ebene, \(k_{\parallel}\). Wenn die optische Berechnung aktiviert ist, werden außerdem Spektren für die verschiedenen Trägerzustände im Beispiel generiert.
5. Elektronische Nebenbänder: Interpretation des Banddiagramms
Öffnen Sie die berechnete Darstellung der Bandstruktur (??). Es lassen sich zwei deutlich unterscheidbare Kurvengruppen erkennen: die Leitungsbänder bei positiven Energien und die Valenzbänder in der Nähe und unterhalb der Valenzbandreferenz. Die vertikale Trennung spiegelt die Bandlücke des Halbleiters sowie die Energien der eingeschlossenen Elektronen und Löcher wider. Die unterschiedlichen Anfangsenergien innerhalb jeder Gruppe ergeben sich aus den quantisierten Zuständen des Potentialsprungs.
Obwohl die Bewegung innerhalb des Potentials begrenzt ist, können sich die Ladungsträger dennoch entlang der Schichten bewegen. Aus diesem Grund wird eine einzelne, lokalisierte Energieform, wenn sie gegen den parallelen Wellenvektor aufgetragen wird, zu einer Parabel:
\[ E_{e,n}(k_{\parallel})=E_{e,n}(0) +\frac{\hbar^2 k_{\parallel}^2}{2m_e^{*}}. \]
Hierbei ist \(E_{e,n}(0)\) die Energie am unteren Ende des Elektronen-Unterbandes \(n\). Für die entsprechenden Valenzelektronen-Unterbänder nimmt die einfache parabolische Form eine nach unten gerichtete Krümmung an:
\[ E_{v,m}(k_{\parallel})=E_{v,m}(0) -\frac{\hbar^2 k_{\parallel}^2}{2m_h^{*}}. \]
Dies ist die entscheidende Annäherung in diesem Tutorial: Die Energieniveaus der Zustände sind in Richtung der Einschränkung quantisiert, behalten aber in der Ebene eine effektive Massenabhängigkeit bei. Fortschrittlichere Mehrband-k·p-Modelle ermöglichen es, dass sich die Bänder vermischen und nicht mehr parabolisch sind; die Seite 8-Band-Zinkblende-Theorie erklärt, warum dies von Bedeutung sein kann.
6. Von den Subbändern zu den optischen Spektren
Elektronen können optische Übergänge zwischen Valenz- und Leitungsunterbändern durchlaufen. Bei normalem Einfall von Licht ist der Impuls des Photons im Vergleich zum elektronischen Wellenvektor gering, daher berücksichtigt ein einfaches Modell Übergänge, bei denen der Wellenvektor in der Ebene erhalten bleibt. Im Zentrum der Dispersionskurve beträgt die angenäherte Photonenenergie für einen Übergang zwischen einem lokalisierten Elektronenenergieniveau \(n\) und einem Loch-Energieniveau \(m\)
\[ \hbar\omega\simeq E_g+E_{e,n}^{\mathrm{conf}}+E_{h,m}^{\mathrm{conf}}. \]
Die beiden Einschlussenergien werden von ihren jeweiligen Bandkanten aus nach oben gemessen. Die verfügbaren Übergänge, die Besetzung ihrer Anfangs- und Endzustände sowie eine eventuelle spektrale Verbreiterung bestimmen die Form der berechneten optischen Antwort. Diese Konzepte können Sie auf der Seite Theorie der optischen Übergänge weiter untersuchen.
6.1 Stimulierte Emission und Nettoverstärkung
Die Grafik zur stimulierten Emission (??) zeigt den berechneten Emissionsbeitrag in Abhängigkeit von der Wellenlänge für die in dem Beispiel verwendeten Trägerzustände. Die deutlichen spektralen Stufen und Maxima hängen mit den Energien zusammen, bei denen verschiedene Übergänge zwischen den Nebenbändern stattfinden. Die separate Grafik zur Nettoverstärkung (??) berücksichtigt den konkurrierenden Beitrag der stimulierten Absorption und kann daher positiv oder negativ sein.
Vereinfacht dargestellt enthält der resultierende, durch Anregung hervorgerufene Beitrag für einen Leitfähigkeitszustand, der mit der Wahrscheinlichkeit \(f_c\) besetzt ist, und einen Valenzelektronenzustand, der mit der Wahrscheinlichkeit \(f_v\) besetzt ist, den Belegungsfaktor.
\[ \underbrace{f_c(1-f_v)}_{\text{stimulated emission}} -\underbrace{f_v(1-f_c)}_{\text{absorption}} =f_c-f_v. \]
Ein positiver Netto-Stimulationsbeitrag erfordert eine ausreichende Inversion der an der Übergangsreaktion beteiligten Zustände. Das insgesamt dargestellte Spektrum hängt zusätzlich von den optischen Übergangsstärken und der Anzahl der verfügbaren Zustände ab.
6.2 Absorption und spontane Emission
Das Absorptionsdiagramm (??) erfasst den Beitrag der erlaubten Übergänge zur Absorption. Die Absorption erfordert einen besetzten Anfangszustand im Valenzband und einen verfügbaren Endzustand im Leitungsband. Das Diagramm für die spontane Emission (??) stellt die radiative Rekombination besetzter Zustände im Leitungsband in leere Zustände im Valenzband dar, ohne dass ein einfallendes Photon erforderlich ist, um den Übergang anzuregen. Entsprechend enthält der Faktor für die Besetzung bei spontaner Emission \(f_c(1-f_v)\).
Beachten Sie, dass die ausgeprägten spektralen Merkmale in ähnlichen Wellenlängenbereichen in allen Diagrammen auftreten. Das ist zu erwarten: Die gleiche zugrunde liegende Struktur von Elektronen- und Loch-Subbändern liefert die möglichen Energien für optische Übergänge, während die Besetzung der Ladungsträger die Absorption von den verschiedenen Arten der Emission unterscheidet. Das Diagramm für die spontane Emission verwendet eine logarithmische vertikale Skala. Achten Sie daher beim Vergleichen der Kurvenhöhen darauf.
7. Die berechneten Dateien finden
Öffnen Sie nach Abschluss der Simulation das Register „Ausgabe“ im Hauptfenster (??). Das Beispiel erstellt separate CSV-Dateien für die verschiedenen spektralen Berechnungen. Zu den sichtbaren Dateinamen gehören Einträge, die dem Muster gain_00.csv mit den Suffixen „Absorption“, „Nettoverstärkung“, „spontan“ und „angeregt“ folgen, gefolgt von ähnlich benannten Dateien für nachfolgend berechnete Bedingungen. Doppelklicken Sie auf ein Ergebnis, um es im Grafikfenster anzuzeigen. Da die Ergebnisse als CSV-Dateien gespeichert werden, können sie auch außerhalb von OghmaNano analysiert oder erneut grafisch dargestellt werden.
Die nummerierten Dateien ermöglichen es Ihnen, die Veränderungen des berechneten Spektrums über die gesamte Reihe der Trägerbedingungen hinweg zu verfolgen, anstatt sich nur ein einzelnes optisches Berechnungsergebnis anzusehen.
8. Probieren Sie es selbst aus: Ändern Sie die Breite des Brunnens.
Vorschlag für ein Experiment: Kehren Sie zum Layer-Editor zurück und ändern Sie die Dicke der GaAs-(Quantenwellen-)Schicht von 10 nm auf 8 nm, wobei die AlGaAs-Barriere-Schichten unverändert bleiben. Führen Sie die Simulation erneut aus und vergleichen Sie die resultierenden Teilbandenergien und optischen Spektren mit dem ursprünglichen Beispiel.
- Verschiebt sich das niedrigste Energieniveau, in dem sich Elektronen aufhalten können, im Verhältnis zum Leitungsbandrand von GaAs nach oben?
- Wie verändern sich die Übergangsenergien von Elektronen und Löchern?
- Verschiebt sich der Beginn des optischen Spektrums zu kürzeren oder längeren Wellenlängen?
✅ Erwartete Entwicklung der physischen Eigenschaften
Im Allgemeinen führt eine schmalere Quantenstruktur zu einer Zunahme der Einschlussenergien. Daher verschiebt sich die niedrigste interbandige Übergangsanregung tendenziell zu einer höheren Photonenenergie, was einer kürzeren Wellenlänge entspricht, obwohl die genaue Verschiebung von den endlichen Barrierenhöhen und den effektiven Massen von Elektronen und Löchern abhängt. Seit \(\lambda=hc/E_{\mathrm{photon}}\) bedeutet eine höhere Übergangsanregungsenergie eine kleinere Emissions- oder Absorptionswellenlänge.
Sobald Sie sich mit diesem Beispiel vertraut gemacht haben, liegt es auf der Hand, dass der nächste Vergleich die 8-Band-Zinkblende-k·p-Version eines GaAs/AlGaAs-Quantenraums ist. Diese behält die Quanteneinschränkung bei, geht aber über unabhängige parabolische Bänder hinaus, um die Wechselwirkungen und das Mischen zwischen den Leitungs- und Valenzzuständen zu beschreiben.