Simulation einer AlGaAs/GaAs-Heterostrukturdiode: Bandanordnung, J–V-Kennlinien und Ladungsträgertransport
1. Einleitung: Warum sollte man eine AlGaAs/GaAs-Heteroverbindung simulieren?
Eine AlGaAs/GaAs-Heteroverbindungdiode ist eine Halbleiterdiode, die an der Grenzfläche zweier unterschiedlicher III–V-Halbleitermaterialien entsteht: Aluminiumgalliumarsenid (AlxGa1−xAs) und Galliumarsenid (GaAs). Im Gegensatz zu einer herkömmlichen GaAs-PN-Diode kombiniert dieses Bauelement eine Änderung der Dotierung mit einer Änderung der Bandstruktur des Halbleiters. Die Grenzfläche weist daher sowohl ein eingebautes elektrostatischs Potential als auch Versetzungen zwischen dem Leitungs- und dem Valenzband auf.
In diesem Tutorial untersuchen wir einen p-Al0,30Ga0,70As-/n-GaAs-Übergang mithilfe des eindimensionalen Poisson- und Drift-Diffusions-Solvers von OghmaNano. Der AlGaAs-Bereich mit seiner größeren Bandlücke trifft auf den GaAs-Bereich mit seiner kleineren Bandlücke; ihre Grenzfläche verändert die Injektion von Elektronen und Löchern und bietet eine einfache Einführung in die Bandanordnung von Heterostrukturübergängen. Die gesamte Halbleiterschicht ist nur 700 nm dick, was die Visualisierung interner Transportprofile erleichtert.
AlGaAs/GaAs-Heterostrukturen sind wichtige Bausteine in der III–V-Elektronik und -Optoelektronik, einschließlich Heteroübergangskomponenten, photoelektrischer Schichtmaterialien und lichtemittierenden Strukturen. Dieses Beispiel dient speziell als Demonstration einer Diode/Gleichrichterfunktion: Es wird nicht als optimierte Photodiode, Laserdiode oder als kommerzielle Hochgeschwindigkeitskomponente präsentiert. Wir werden das Modell verwenden, um eine messbare Kennlinie für den Dunkelstrom in Abhängigkeit von der Spannung (J–V) mit dem mikroskopischen Banddiagramm und den Rekombinationsprozessen in Beziehung zu setzen.
Die Hauptansicht wird in ?? angezeigt. Der Rest der Anleitung folgt dem gleichen Ablauf wie die anderen Beispiele für OghmaNano-Dioden: Laden Sie ein Bauelement, überprüfen Sie dessen Eingänge, führen Sie eine Simulation durch und interpretieren Sie die Ergebnisse.
2. Erstellen einer neuen Simulation
Wählen Sie im Hauptfenster von OghmaNano die Option „Neue Simulation“ aus. Der Demonstrationsbrowser gruppiert verwandte Beispiele für Halbleiter. Öffnen Sie die GaAs/III–V-Demonstrationen und wählen Sie das Beispiel für die AlGaAs/GaAs-Heterostrukturdiode aus, oder laden Sie das mitgelieferte Beispielprojekt direkt, falls sich die Menübezeichnung zwischen den Versionen unterscheidet. Die beiden Phasen des Browsers sind in ?? und ?? dargestellt.
Sobald das Beispiel geöffnet ist, untersuchen Sie die vertikale Anordnung in der Geräteansicht. Obwohl die Ansicht dreidimensional ist, behandelt dieses Tutorial die eindimensionale elektrische Leitfähigkeit durch die Schichtdicke. Das Modell berücksichtigt weder die laterale Stromkonzentration noch die tatsächlichen Geometrien der Metallkontakte.
3. Untersuchung der Schichtstruktur, der Legierungszusammensetzung und der Dotierung
3.1 Schichtaufbau
Öffnen Sie den Layer-Editor über die Symbolleiste des Geräts. Die aktive Diode besteht aus drei Halbleiterschichten, wobei sich oberhalb dieser eine Luftschicht befindet, die für die Darstellung des Geräts verwendet wird. Das elektrische Problem betrifft den 700 nm dicken Halbleiterschichtstapel und nicht die 200 nm dicke Luftschicht.
| Region (von oben nach unten) | Dicke | Doping beim Model | Zweck |
|---|---|---|---|
| Luft (nicht elektrisch leitfähig) | 200 Nanometer | — | Umgebung über dem Gerät |
| p-Al0,30Ga0,70As | 200 Nanometer | NA = 5 × 1017 cm−3 | Weitbandabstand-p-Seite einer Heteroverbindung |
| n-dotiertes Galliumarsenid | 400 nm | ND = 5 × 1017 cm−3 | n-Dotierte Schicht mit geringerer Bandlücke |
| n+-Galliumarsenid | 100 Nanometer | ND = 1 × 1018 cm−3 | Stärker dotierte Region auf der Kathodenseite |
Die p-AlGaAs/n-GaAs-Grenzschicht stellt sowohl den metallurgischen PN-Übergang als auch die Materialgrenzfläche dar. Die n-GaAs/n+-GaAs-Grenzschicht ist ein Dotierungsprofil innerhalb von GaAs und keine zweite Heteroverbindung. Der Schichteneditor (??) und die Kontaktkonfiguration (??) werden im Folgenden zusammen dargestellt.
3.2 Dotierung und Aluminiumzusammensetzung
Öffnen Sie den Doping-/Ionen-Editor über die Registerkarte „Elektrisch“. Die mit Akzeptoren dotierte AlGaAs-Seite und die mit Donatoren dotierte GaAs-Seite bilden den Gleichrichterübergang (??). Die Dotierung ändert ihr Vorzeichen an der Materialgrenzfläche, wodurch ein eingebautes elektrisches Feld entsteht. Die n+-Region erhöht die Elektronendichte in der Nähe der Kathode.
Öffnen Sie als Nächstes den Editor für die Zusammensetzung/den Gradienten der Legierung. Die AlGaAs-Legierung hat einen Aluminium-Molenbruch von x = 0,30, während GaAs einem Wert von x = 0 entspricht. Der abrupte Zusammensetzungsübergang (??) ist einer der Hauptgründe, warum das Energiediagramm neben der üblichen elektrostatischen Bandbiegung auch einen Sprung aufweist.
3.3 Warum die Benutzeroberfläche wichtig ist
Bei Raumtemperatur hat GaAs eine Bandlücke von etwa 1,42 eV, während Al0,30Ga0,70As eine größere Bandlücke aufweist, die je nach Legierungsparametrisierung oft bei etwa 1,8 eV liegt. In einer Heterostruktur wird die gesamte Bandlückenunterschied auf die Leitungs- und Valenzbandverschiebung verteilt:
\[\Delta E_g = \Delta E_c + \Delta E_v\]Hier bezeichnen die Verschiebungssymbole die Größen der entsprechenden Sprungstellen. Ihre numerischen Werte werden durch das Materialmodell und die Annahmen zur Bandanordnung bestimmt und nicht allein durch das Dotierungsprofil. Bei der späteren Analyse der simulierten Bandkanten ist zu unterscheiden zwischen einer ausgeprägten, materialbedingten Verschiebung und der allmählichen Bandbiegung, die durch die Raumladung verursacht wird.
4. Elektrische Materialparameter, Schnittstellen und Rekombination
4.1 Die Materialdatenbank und die Lua-Modelle
Öffnen Sie den Editor für die elektrischen Parameter für jede Schicht. In diesem Beispiel wird Inorganic/AlGaAs/AlGaAs.lua für die p-dotierte Legierung und Inorganic/GaAs/GaAs.lua für beide n-dotierten GaAs-Bereiche verwendet. Das AlGaAs-Lua-Modell (??) und seine Parameter (??) können mit dem GaAs-Modell (??) und seinen Parametern (??) verglichen werden.
Die Simulation ist so konfiguriert, dass sie wichtige Materialwerte – einschließlich Bandlücke, Elektronenaffinität, Trägermobilität, effektive Zustandsdichte und Rekombinationsparameter – aus der Lua-Materialdatenbank bezieht, anstatt lokale Überschreibungen zu erzwingen. Dies ist besonders wichtig für AlGaAs, da seine Parameter von der Legierungszusammensetzung abhängen. Lesen Sie die auf dem Bildschirm angezeigten Werte, die vom Modell zurückgegeben werden, anstatt beliebige inaktive Standardwerte in der Projektdatei als tatsächliche Materialeigenschaften zu interpretieren.
4.2 Bandanpassung und Schnittstelleneinstellungen
Öffnen Sie den Schnittstellen-Editor, der in ?? dargestellt ist. Dieses Beispiel zeigt eine abrupte Schnittstelle. Bei einer grundlegenden Drift-Diffusions-Betrachtung hängt die Trägerinjektion von den Leitungs- und Valenzband-Diskontinuitäten, der Trägerstatistik und der Behandlung der Schnittstellengrenze ab. Reale Proben können zusätzlich Schnittstellenzustände, Legierungsstörungen, thermionische Emission, Tunnelung oder durch die Verarbeitung verursachte Defekte aufweisen; diese Effekte sollten nicht einfach deshalb als berücksichtigt angenommen werden, weil eine Heteroverbindung vorhanden ist.
4.3 Rekombination durch Strahlung und durch Defekte begünstigte Prozesse
GaAs und Al0,30Ga0,70As sind bei dieser Zusammensetzung Materialien mit direkter Bandlücke, sodass unter Trägerinjektion eine strahlungsinduzierte Rekombination ohne Energieverlust stattfinden kann. Ein gängiges Modell lautet:
\[R_{\mathrm{rad}} = B\left(np-n_i^2\right)\]wobei B ein vom Material abhängiger Koeffizient ist und n und p die lokalen Konzentrationen freier Ladungsträger darstellen. Die gewählte Konfiguration ermöglicht außerdem eine einstufige Shockley-Read-Hall-(SRH-)Rekombination in den drei aktiven Halbleiterbereichen. Eine idealisierte Formel lautet:
\[R_{\mathrm{SRH}}=\frac{np-n_i^2}{\tau_p(n+n_1)+\tau_n(p+p_1)}\]SRH beschreibt die durch Fallen unterstützte Rekombination und ist empfindlich gegenüber den durch die gewählten Modelle bereitgestellten Fallenparametern. In diesem Beispiel ist die Auger-Rekombination deaktiviert. Man sollte nicht den Schluss ziehen, dass ein einzelner Rekombinationsprozess vorherrscht, ohne dessen Stärke mit der der anderen aktivierten Prozesse unter den gleichen Vorspannungsbedingungen zu vergleichen.
Die Anode und die Kathode verwenden idealisierte ohmsche Randbedingungen. Diese helfen dabei, die physikalischen Eigenschaften des Heteroübergangs zu isolieren, ermöglichen aber keine Vorhersage des Verhaltens einer bestimmten Metallisierungsmethode oder ihres Kontaktwiderstands.
5. Durchführung der Simulation mit dem dunklen J–V-Modell
Kehren Sie zum Hauptfenster der Simulation zurück und klicken Sie auf „Simulation starten“. Die Eingabedatei legt eine Spannungsänderung von -0,1 V bis +2,5 V in Schritten von 0,025 V fest, ohne dass eine gezielte optische Beleuchtung erfolgt. OghmaNano löst die gekoppelten Poisson- und Kontinuitätsgleichungen für Elektronen und Löcher an aufeinanderfolgenden Arbeitspunkten.
Während der Berechnung können Sie zum Tab „Terminal“ wechseln, um den Fortschritt des Solvers zu beobachten. Nach Abschluss der Berechnung können Sie dann die Ausgabe öffnen. Der Ausgabebrowser (??) und das Terminal (??) veranschaulichen diese Phasen. Die Screenshots des Terminals werden zu unterschiedlichen Zeitpunkten aufgenommen, nicht bei festen elektrischen Vorspannungen.
Doppelklicken Sie auf jv.csv, um die simulierte Stromdichte-Spannungskurve (??) zu öffnen. Eine Gleichrichtercharakteristik sollte bei moderater Sperrspannung einen relativ schwachen Strom und bei zunehmender Injektion einen immer größeren Vorwärtsstrom aufweisen. Bei dieser speziellen Heteroverbindung hängen der Beginn und die detaillierte Form von den Bandversatzwerten, den Materialparametern und den Rekombinationsraten ab – und nicht nur von der Bandlücke des GaAs.
Die Formel für die ideale Diode kann als nützliches Hilfsmittel für die lokale Interpretation dienen:
\[J \approx J_0\left[\exp\!\left(\frac{qV}{nk_{\mathrm B}T}\right)-1\right]\]Allerdings sollte ein angepasster Idealitätsfaktor n oder ein Sättigungsstrom J0 nur aus einem geeigneten Spannungsbereich ermittelt werden. Reihenwiderstand, hohe Injektion und Heteroflächen-Transport können zu Abweichungen von einem einfachen exponentiellen Gesetz führen. Der beigefügte Screenshot dient lediglich der Veranschaulichung und stellt keine unabhängige experimentelle Bestätigung dar.
6. Analyse von Simulationsmomentaufnahmen: Was geschieht im Inneren der Diode?
OghmaNano speichert während der J–V-Messung räumlich aufgelöste elektrische Ergebnisse. Öffnen Sie die Ansicht für Momentaufnahmen und verwenden Sie den Schieberegler, um ausgewählte Bias-Positionen zu vergleichen. Die vier Bilder in jeder Zeile unten entsprechen dem Anfang, ungefähr einem Drittel, ungefähr zwei Dritteln und dem Ende der gespeicherten Messung. Die Beschriftungen 0–3 geben die Positionen der Screenshots an, nicht die Spannungen in Volt; verwenden Sie die Bias-Anzeige des Betrachters, um die genauen Werte für eine bestimmte Messreihe zu erhalten.
6.1 Leitungsband, Valenzband und Quasi-Fermi-Niveaus
Stellen Sie die Diagramme Ec.csv, Ev.csv, Fn.csv und Fp.csv zusammen dar. Der erste Screenshot des Banddiagramms (??) zeigt den unteren Rand des Leitungsbandes (Ec) und den oberen Rand des Valenzbandes (Ev) und veranschaulicht sowohl das gleichmäßige elektrostatische Potenzial als auch die abrupte Banddiskontinuität an der Grenzfläche zwischen AlGaAs und GaAs. Die quasi-Fermi-Niveaus für Elektronen und Löcher zeigen Abweichungen vom Gleichgewichtszustand und die Verteilung der injizierten Ladungsträger.
Vergleichen Sie das erste und das letzte Bild. Bei Durchlassvorspannung verändert sich die Sperrschicht und die quasi-Fermi-Niveaus werden getrennt. Verwechseln Sie einen Bandkanten-Sprung an der Grenzfläche nicht mit numerischem Rauschen: eine echte Heteroverbindung kann auch dann diskontinuierliche Bandkanten aufweisen, wenn das elektrostatische Potenzial kontinuierlich bleibt.
6.2 Elektron- und Lochdichten
Stellen Sie in diesem Projekt die Ausgabewerte für die Konzentration freier Elektronen und freier Löcher, Q_nfree.csv bzw. Q_pfree.csv, grafisch dar. Die p-AlGaAs-Seite besteht hauptsächlich aus Lochträgermaterial, während die GaAs-Seite hauptsächlich aus Elektronenträgermaterial besteht. Eine Erhöhung der Durchlassspannung führt dazu, dass Minderheitsträger über den Übergang injiziert werden, wodurch sich die räumliche Verteilung der Träger ändert. Wie die Ergebnisse zur Trägerdichte (??) zeigen, verändert eine abrupte Materialgrenze die verfügbaren elektronischen Zustände und die Bandanordnung, sodass das Trägerdichteprofil nicht unbedingt dem einer homogenen GaAs-PN-Verbindung ähneln muss.
6.3 Elektronen- und Lochstrom
Stellen Sie die Diagramme Jn.csv und Jp.csv so auf, dass die Beiträge von Elektronen und Löchern zum Gesamtstrom getrennt dargestellt werden. Im stationären Zustand können sich die beiden Komponenten in Abhängigkeit von der Position neu verteilen, da Elektronen und Löcher rekombinieren, während die Summe der vorzeichenbehafteten elektrischen Stromdichten in Abwesenheit zusätzlicher Stromquellen ungefähr konstant bleiben sollte. Lokale Spitzen oder rasche Änderungen in der Nähe der Heteroverbindung in den Diagrammen des Trägerstroms (??) sollten anhand des Gitters und des Banddiagramms überprüft werden.
Bei der Analyse von Stromverläufen mit geringen Stromstärken ist Vorsicht geboten. Beurteilen Sie den allgemeinen Trend und die Kontinuität des Stromverlaufs, anstatt jede kleine Oszillation auf einen physikalischen Transportmechanismus zurückzuführen. Wenn Sie große Abweichungen feststellen, verfeinern Sie das Netz und prüfen Sie die Konvergenz des Solvers, bevor Sie Schlussfolgerungen ziehen.
6.4 Strahlungsbedingte und SRH-Rekombination
Untersuchen Sie R_nfree_to_pfree.csv (freie-zu-freie radiative Rekombination) und R_ss_srh.csv (Rekombination auf einer einzelnen Energieniveau, SRH-Rekombination). Die Screenshots der Rekombinationsprozesse (??) zeigen diese Kanäle zusammen. Unter Injektionsbedingungen kann die absolute Rekombinationsrate stark ansteigen; ihre Position wird durch die Überlappung freier Ladungsträger, die Materialeigenschaften und die gewählten Fallenparameter bestimmt.
Achten Sie besonders auf die Eigenschaften an der p-AlGaAs/n-GaAs-Grenzfläche. Eine lokale Eigenschaft kann physikalisch bedeutsam sein, aber ihre Größe allein beweist nicht, dass es sich um eine Ladungsträgerspeicherung an der Grenzfläche handelt: In der vorliegenden Übung werden Modelle für das Volumenmaterial und die konfigurierten SRH-Parameter verwendet. Vergleichen Sie den integrierten Beitrag und die räumliche Verteilung beider Prozesse, bevor Sie einen dominanten Verlustmechanismus festlegen.
6.5 Profile der Elektronen- und Lochbeweglichkeitswerte
Die getrennten Ergebnisse für R_n.csv und R_p.csv zeigen die für die Simulation verwendeten Gleichgewichtsbeziehungen zwischen Elektronen- und Lochkonzentrationen. Diese Diagramme (??) ermöglichen eine zusätzliche Überprüfung, wo Ladungsträger während des Messvorgangs entfernt oder umverteilt werden. Interpretieren Sie die angezeigten Werte und die grafisch dargestellten Komponenten anhand der jeweiligen Datenüberschriften und gehen Sie nicht davon aus, dass alle Kurven eine einzige positive Rekombinationsrate darstellen.
6.6 Trägermobilitäten
Die Mobilitätsausgabe ermöglicht die Analyse von Transportparametern in Abhängigkeit von der Position. Änderungen an der Grenzfläche sind zu erwarten, da die AlGaAs- und GaAs-Modelle nicht unbedingt die gleiche Mobilität liefern müssen. Das Beispiel für einen Mobilitäts-Screenshot (??) zeigt eine gespeicherte Anzeigeposition; es handelt sich nicht um einen Beweis für einen bestimmten feldabhängigen Mobilitätseffekt.
7. Übungen: Zusammenhang zwischen den Bandversatzwerten und dem Verhalten von Dioden herstellen
8. Annahmen und Einschränkungen des Beispiels
Dies ist ein idealisiertes Lehrbeispiel für eine eindimensionale Struktur. Die Halbleiterdicken und Dotierkonzentrationen sind für eine epitaktisch aufgebrachte III–V-Schichtstruktur plausibel, aber das Beispiel allein führt nicht zu einem herstellbaren oder hinsichtlich der Leistung optimierten Produkt. Die Berechnung geht von idealisierten ohmschen Kontakten, einer lateral gleichförmigen Struktur, einer abrupten Materialgrenzfläche und den durch die ausgewählten Modelle bereitgestellten Transport- und Rekombinationsmechanismen aus.
Die AlGaAs/GaAs-Grenzfläche kann in realen Bauelementen zusätzliche Transportphänomene aufweisen; praktische Kontakte, Defekte, Oberflächenrekombination, Legierungsstörungen, Selbstaufheizung und Substrateffekte können ebenfalls die gemessenen Eigenschaften verändern. In dieser Eingabe ist die Auger-Rekombination deaktiviert, und die Simulation ist eine elektrische Messung im Dunkeln, nicht eine kalibrierte LED-, Photodioden- oder Solarzellen-Effizienzvorhersage. Eine quantitative Validierung erfordert experimentelle Daten und unabhängig überprüfte Materialparameter.
9. Nächste Schritte
Nachdem die J–V-Kurve reproduziert und die Heteroverbindung im Banddiagramm lokalisiert wurde, ist die nächste aussagekräftige Berechnung ein kontrollierter Vergleich der Legierungszusammensetzung: Variieren Sie den Aluminiumanteil, halten Sie die Dotierung und Dicke konstant und verfolgen Sie den Sprung in der Leitbandkante, den Sprung in der Valenzbandkante und den Strom bei Durchlassspannung. Eine Berechnung für eine GaAs-Homoverbindung liefert eine nützliche Referenz.
Weitere Informationen darüber, wie OghmaNano Gerätekonfigurationen und Materialmodellparameter speichert, finden Sie unter Dokumentation zum erweiterten Oghma-Dateiformat.