GaAs/AlGaAs-Quantenpotentialtopf: Einführung in das 8-Band-Zinkblende-k·p-Modell

1. Einleitung

In diesem Tutorial werden wir die elektronischen und optischen Eigenschaften einer GaAs/AlGaAs-Quantenstruktur mithilfe des Acht-Band-Zinkblende-k·p-Modells von OghmaNano berechnen. Wir beginnen mit einem vorgefertigten Beispiel, untersuchen die Schichtstruktur, führen die Simulation durch und analysieren die resultierenden Leitungs- und Valenzteilbandstrukturen. Abschließend betrachten wir die aus den berechneten Zuständen ermittelten optischen Spektren.

Das Tutorial zum parabolischen Quantenpotential behandelt Elektronen und Löcher mithilfe vergleichsweise einfacher effektiver Massebandstrukturen. Hierbei enthält der Hamilton-Operator eine Leitbandsbasis sowie die Basen für das schwere Loch (HH), das leichte Loch (LH) und das abgespaltene Valenzband (SO), jeweils mit zwei Spin-Komponenten, was insgesamt acht Basiszustände ergibt. Die Kopplung zwischen diesen Zuständen ermöglicht es dem Modell, Effekte zu beschreiben, die durch unabhängige Parabeln nicht erfasst werden können, insbesondere die Mischung von Valenzbändern und Änderungen der Krümmung in Abhängigkeit vom Wellenvektor in der Ebene.

2. Öffnen des Beispiels mit acht Frequenzbändern

In OghmaNano klicken Sie im Menüband „Datei“ auf „Neue Simulation“. Doppelklicken Sie auf den Ordner „GaAs-Beispiele“ (??). Wählen Sie aus den Beispielen für Quantenpotentiale „GaAs/AlGaAs QW 8 Band k.p“ aus, indem Sie darauf doppelklicken (??). Dadurch wird eine voreingestellte Simulation geöffnet; es ist nicht erforderlich, eine 8x8-Hamilton-Matrix manuell einzugeben.

3. Quantenwellenstruktur und Materialparameter

Das Hauptfenster zeigt eine dreischichtige AlGaAs/GaAs/AlGaAs-Struktur (??). Öffnen Sie den Schichteneditor mit der Schaltfläche auf der linken Seite. Das mitgelieferte Beispiel besteht aus zwei 40 nm dicken AlGaAs-Barrieren, die eine 10 nm dicke GaAs-Schicht umschließen:

Die Bandkanten des GaAs liegen innerhalb der umgebenden Bandlücke des AlGaAs, was eine räumliche Begrenzung in Wachstrichtung der Schicht ermöglicht. Im Schichteditor ist die zentrale GaAs-Schicht speziell als Quantenpotential (??) gekennzeichnet. Die optischen Materialeinträge und Schichttypen sind in diesem Beispiel bereits konfiguriert.

Der Solver für acht Bänder benötigt mehr Materialinformationen als ein skalares parabolisches Bandmodell: Zusätzlich zu den Bandrändern und den effektiven Massen werden Parameter verwendet, die die interbandische Kopplung, die Aufspaltungsenergie und (falls zutreffend) die Dehnung beschreiben, und diese werden in den Hamilton-Operator eingegeben. Sie können die GaAs-Funktionen über die elektrischen Parameter und den Code-Editor des Materials (??) einsehen. Der Screenshot zeigt die Verformungspotenziale des Valenzbandes b und d; diese müssen für diese erste Testausführung nicht geändert werden.

Welche Änderungen gibt es bei der Berechnung mit den acht Frequenzbändern?

Ein parabolisches Modell ordnet jeder Bandstruktur eine eigene Krümmung zu. Der k·p-Ansatz konstruiert stattdessen eine Matrix, deren nicht-diagonale Elemente die zugrunde liegenden Basiszustände miteinander koppeln. Schematisch dargestellt, ist das Problem der eingeschränkten Zustände:

\[ H_{8\times8}(k_x,k_y,-i\partial_z)\,\Psi_n(z) =E_n(k_x,k_y)\,\Psi_n(z). \]

Hierbei bezeichnet \(z\) die Richtung der Einschränkung; \(k_x\) und \(k_y\) beschreiben die Bewegung parallel zu den Schichten. Die Hüllwellenfunktion \(\Psi_n\) hat acht Komponenten, die die Beiträge der verschiedenen Basiszustände beschreiben. Für dieses Tutorial konzentrieren wir uns auf die berechneten Energien und nicht auf die vollständige Matrixalgebra. Eine detaillierte Erklärung finden Sie im Kapitel Theorie der achtfachen Bandstruktur von Zinkblende.

4. Durchführung der Simulation

Kehren Sie zur Registerkarte „Datei“ zurück und klicken Sie auf die blaue Schaltfläche „Simulation starten“. Das Acht-Band-Hamilton-System wird zusammengestellt und für eine Reihe von Wellenvektoren in der Ebene gelöst. Sie können den Berechnungsprozess auf der Registerkarte „Terminal“ verfolgen (??).

Das Beispielterminal zeigt, wie die dünnbesetzte Matrix aufgebaut wird, gefolgt von den Energiewerten der Eigenzustände und Meldungen wie block converged after 12 iterations. Leitungs- und Valenzzustände werden bei aufeinanderfolgenden Wellenvektoren separat berechnet. Warten Sie, bis die Berechnung abgeschlossen ist und die Schaltfläche wieder den Zustand „Simulation ausführen“ anzeigt, bevor Sie die Ergebnisse überprüfen.

5. Das berechnete Bandmodell verstehen

Der Solver berechnet für jeden in der Ebene liegenden Wellenvektor \(k_{\parallel}\) eine diskrete Menge von räumlich begrenzten Eigenzuständen. Die Grafik zeigt daher die Subband-Dispersion und nicht das räumliche Bandprofil über die Schichten hinweg. Öffnen Sie die Registerkarte „Ausgabe“ und suchen Sie im Ordner bands oder öffnen Sie die zusammenfassenden Banddateien, die in ?? angezeigt werden.

5.1 Leitungsunterbänder

Die Darstellung des Leitungsbandes (??) zeigt zwei nach oben gebogene Zweige. An der Stelle \(k_{\parallel}=0\) entsteht die vertikale Trennung aufgrund der unterschiedlichen räumlich begrenzten Zustände. Die Zunahme der Energie bei \(k_{\parallel}\) spiegelt eine Bewegung parallel zur Potentialmulde wider, wobei die Ladungsträger nicht räumlich begrenzt sind.

In der Nähe des Minimums eines Leitungsband-Nebenbandes kann die Dispersionsbeziehung durch folgende Formel angenähert werden:

\[ E_n(k_{\parallel})\simeq E_n(0) +\frac{\hbar^2 k_{\parallel}^{2}}{2m_{n,\parallel}^{*}}. \]

Im Mehrbandmodell handelt es sich hierbei um eine lokale Approximation: Die Krümmung und damit die scheinbare effektive Masse müssen nicht über den gesamten dargestellten Bereich der Wellenvektoren konstant bleiben. Die Kopplung zwischen Leitungs- und Valenzband trägt zu diesem Verhalten bei.

5.2 Valenzteilbandstrukturen: Mischung und Krümmung

Das Diagramm des Valenzbandes ist interessanter (??). Seine beiden dargestellten Zweige beginnen bei unterschiedlichen Energien, biegen sich nach unten und scheinen sich bei einem endlichen Wert von \(k_{\parallel}\) zu schneiden. Im Zinkblende-Valenzband verhalten sich die Zustände der schweren und leichten Löcher nicht wie zwei vollständig unabhängige Parabeln: bei einem endlichen Wellenvektor können die Terme in der Hamilton-Funktion ihre Eigenschaften vermischen.

Eine sinnvolle Methode, dies auszudrücken, besteht darin, zwei benachbarte Zustände zu betrachten, die durch ein Matrixelement \(V\) miteinander verbunden sind.

\[ H_{\mathrm{illustrative}}= \begin{pmatrix}E_{\mathrm{HH}}(k)&V(k)\\ V^{*}(k)&E_{\mathrm{LH}}(k)\end{pmatrix}. \]

Wenn sie sich paaren dürfen, stoßen sich ihre Energien in der Nähe einer Begegnung ab, und ihr HH/LH-Charakter kann zwischen den Zweigen ausgetauscht werden. Eine echte Kreuzung ist auch möglich, bei der die betreffende Paarung aufgrund von Symmetrie verschwindet. Allein anhand des vorliegenden Screenshots lassen sich diese Fälle nicht unterscheiden: Um den Zustandscharakter zu bestimmen, müssen die entsprechenden Eigenvektoren untersucht und gegebenenfalls weitere Punkte in der Nähe des scheinbaren Schnittpunkts berechnet werden. Die oben dargestellte 2x2-Matrix ist lediglich eine Illustration und nicht der vollständige Acht-Band-Hamilton-Operator.

5.3 Gemeinsame Betrachtung von Leitungs- und Valenzzuständen

Die kombinierte Darstellung (??) ordnet die Elektron- und Valenzelektronenzustände einer gemeinsamen Energieskala zu. In diesem Beispiel liegen die Valenzzustände nahe dem Energienullpunkt, während die räumlich begrenzten Leitungszustände bei etwa 1,5 eV und darüber liegen. Die Energie­differenz zwischen einem bestimmten Paar von Zuständen entspricht der Energie eines möglichen vertikalen optischen Übergangs, wobei das jeweilige optische Matrixelement und die Besetzung der Trägerzustände berücksichtigt werden müssen.

Da der Impuls eines Photons relativ gering ist, ist die erste sinnvolle Näherung, dass ein optischer Übergang den Wellenvektor in der Ebene beibehält:

\[ \hbar\omega=E_{c,n}(k_{\parallel})-E_{v,m}(k_{\parallel}), \qquad \lambda=\frac{hc}{\hbar\omega}. \]

Verschiedene Paare von Elektron- und Lochzuständen tragen daher bei unterschiedlichen Wellenlängen bei. Die mehrbandigen Wellenfunktionen bestimmen auch die Stärke und Polarisation jeder Übergangsreaktion, und zwar in einem Umfang, der über das hinausgeht, was allein aus den Energiewerten abgelesen werden kann.

6. Optische Eigenschaften des Quantenbrunnens

Die Berechnung mit acht Bändern liefert elektronische Energien und Wellenfunktionen für gemischte Zustände für das optische Modell. Das Beispiel erzeugt separate Spektren für stimulierte Emission, Absorption, Nettoverstärkung und spontane Emission. Die einzelnen Kurven in jeder Grafik entsprechen der Reihenfolge der Trägerbedingungen, die in der bereitgestellten Simulation verwendet werden; es handelt sich nicht um acht einzelne Hamilton-Bänder.

6.1 Anregte Emission und Absorption

Die stimulierte Emission tritt auf, wenn ein einfallendes Photon ein Elektron dazu anregt, einen Übergang von einem besetzten Leitungszustand zu einem freien Valenzzustand zu vollziehen, wobei Energie in Form von Strahlung abgegeben wird. Die Absorption ist der umgekehrte Prozess: Ein Valenzelektron nimmt ein Photon auf und bewegt sich in einen freien Leitungszustand. Ihre Besetzungszahlen sind schematisch dargestellt als:

\[ \begin{aligned} \text{stimulated emission:}&\quad f_c(1-f_v),\\ \text{absorption:}&\quad f_v(1-f_c). \end{aligned} \]

Die Kurven für die stimulierte Emission (??) und die Absorptionskurven (??) haben die gleichen zugrunde liegenden Übergangenergien, aber ihre Intensitäten unterscheiden sich, je nachdem, wie sich die Besetzung ändert. In diesen Ergebnissen liegt ein Großteil der spektralen Struktur um 790–850 nm; die genauen Merkmale hängen von den konfigurierten Subbändern, den Matrixelementen, den Trägerkonzentrationen und der Linienverbreiterung ab.

6.2 Netto-Emissionsgewinn: tritt auf, wenn die Emission die Absorption übersteigt.

Die Grafik der Netto-Verstärkung (??) fasst die Beiträge der konkurrierenden stimulierten Emission und Absorption zusammen. Wenn man von diesen Beiträgen ihre jeweiligen Besetzungswahrscheinlichkeiten abzieht, erhält man ein besonders einfaches Ergebnis:

\[ f_c(1-f_v)-f_v(1-f_c)=f_c-f_v. \]

Demnach führt der Übergang zu einer Verstärkung, wenn der beteiligte Leitfähigkeitszustand stärker besetzt ist als der entsprechende Valenzelektronenzustand. Negative Nettoverstärkungswerte deuten auf Bedingungen hin, die von Absorption dominiert werden; positive Werte entsprechen einer stimulierten Emission, die die Absorption übersteigt. Der hier dargestellte Nettoverstärkungskoeffizient ist das Ergebnis einer Quantenwellenoptik; er sollte nicht automatisch als die Modenverstärkung einer vollständigen Laserkavität interpretiert werden, die auch von der Einschränkung und anderen optischen Verlusten abhängt.

6.3 Spontane Emission

Spontane Emission benötigt kein einfallendes Photon, um den Übergang auszulösen. Sie tritt auf, wenn sich ein besetzter Elektronenenergiezustand durch radiative Rekombination in einen verfügbaren Valenzzustand umwandelt. Der entsprechende Belegungsfaktor ist wiederum \(f_c(1-f_v)\), aber spontane und stimulierte Emission haben unterschiedliche spektrale Vorfaktoren und physikalische Interpretationen. Die Spektren in ?? verwenden eine logarithmische vertikale Skala; daher bedeutet eine gleiche vertikale Abfolge nicht, dass die absolute Zunahme der emittierten Photonenrate gleich ist.

Eine ausführlichere Erklärung der Übergangsstärken, der Trägerberufe und der Linienverbreiterung finden Sie auf den Seiten optische Übergänge und Verstärkungstheorie.

7. Die Ausgabedateien finden

Nachdem die Simulation abgeschlossen ist, wählen Sie im Hauptfenster die Option „Ausgabe“ (??). Zu den angezeigten Ergebnissen gehören Ec.csv, Ev.csv und Ec_Ev.csv für die Banddiagramme sowie gain.csv, absorption.csv, net_gain.csv und spontaneous.csv. Die Ordner bands und gain enthalten weitere Ausgabedaten des Solvers. Doppelklicken Sie auf eine CSV-Datei, um das entsprechende Diagramm in OghmaNano zu öffnen, oder exportieren Sie die numerischen Daten für die weitere Analyse.

8. Probieren Sie es selbst aus: Vergleichen Sie die Ergebnisse mit dem parabolischen Modell.

Als zweites Experiment verändern Sie die Dicke der GaAs-Schicht von 10 nm auf 8 nm, wobei die Dicke der beiden 40 nm dicken AlGaAs-Barrieren unverändert bleibt. Führen Sie das Modell erneut aus und vergleichen Sie die berechneten Eigenenergien und den Wellenlängenbereich der optischen Übergänge.