InP-PN-Übergangsdiode: Simulation von Drift-Diffusion, Gleichrichtung und Trägerrekombination
1. Einleitung: Eine Gleichrichterdiode aus Indiumphosphid
Indiumphosphid (InP) ist ein Halbleiter der Gruppe III–V mit direkter Bandlücke, der als Grundlage für Hochgeschwindigkeits-Elektronikbauelemente und integrierte photonische Technologien dient. Eine PN-Übergangsdiode ist eine der einfachsten und lehrreichsten Bauelementstrukturen: Eine p-dotierte Region trifft auf eine n-dotierte Region, wodurch ein internes elektrisches Feld und eine Potentialbarriere entstehen, die die Bewegung von Elektronen und Löchern steuern. Die Diode leitet bei Durchlassspannung stark und unterdrückt die Trägerinjektion bei Sperrspannung, was zu einer elektrischen Gleichrichtung führt.
InP besitzt bei Raumtemperatur eine Bandlücke von etwa 1,34 eV und eine vergleichsweise hohe Elektronenbeweglichkeit. Diese Eigenschaften machen den Ladungstransport in InP-Übergängen für die Elektronik von Verbindungshalbleitern relevant und ermöglichen ein besseres Verständnis der Halbleiterschichtübergänge, die in InP-basierten Geräteplattformen verwendet werden. Das Beispiel in diesem Tutorial veranschaulicht den elektrischen Betrieb einer kompakten p-InP/n-InP-Homojunktion mit einer n⁺-Kontaktseitenschicht.
Anstatt die J–V-Kurve allein an eine ideale Diodengleichung anzupassen, löst OghmaNano die gekoppelten Poisson- und Drift-Diffusions-Kontinuitätsgleichungen. Die resultierenden Daten zeigen sowohl den Strom an den Anschlüssen als auch die interne Verteilung der Leitungs- und Valenzbandenergien, der Elektron- und Lochdichten, der quasi-Fermi-Niveaus und der Rekombination. Dies ermöglicht es uns, zu erklären, warum die Diode gleichrichtet, anstatt einfach nur festzustellen, dass sie dies tut.
Das Hauptfenster der Simulation, das in ?? dargestellt ist, zeigt die vertikal geschichtete InP-Komponente. Obwohl die Geometrie in drei Dimensionen dargestellt wird, ist das hier gelöste elektrische Problem eindimensional und bezieht sich auf die Wachstumsrichtung.
2. Erstellen einer neuen Simulation
Im Hauptfenster von OghmaNano wählen Sie „Neue Simulation“ aus. Das Fenster zur Auswahl der Simulationen wird geöffnet und zeigt die verfügbaren Beispielkategorien an (??). Öffnen Sie die Demonstrationen für III–V/InP und wählen Sie das Beispiel „InP-PN-Diode“ aus (??). Die genaue Position des Beispiels kann je nach OghmaNano-Version variieren.
Nachdem das Beispiel geöffnet wurde, wählen Sie die Registerkarte „Gerätestruktur“ aus und öffnen Sie den Layereditor. Überprüfen Sie die Layernamen, deren Dicken und die zugehörigen Modelle für das elektrische Material, bevor Sie die Simulation starten. In den folgenden Abschnitten wird erläutert, welchen Beitrag jede Schicht zur Gesamtcharakteristik der Diode leistet.
3. Untersuchung der Bauelemente, Kontakte und Dotierung.
3.1 Struktur der Halbleiterschicht
Das elektrische Bauelement enthält drei InP-Schichten mit einer Gesamt-Halbleiterdicke von 700 nm. Die p-InP- und n-InP-Schichten bilden den gleichartigen Gleichrichterübergang, während die n⁺-InP-Region eine stärker dotierte Kathodenseite bildet. Die Zuordnung der Schichten ist im Editor sichtbar (??).
| Schicht | Dicke | Nominale Dotierung | Rolle |
|---|---|---|---|
| p-InP | 250 Nanometer | \(N_A=5\times10^{17}\,\mathrm{cm^{-3}}\) | Bereich mit einem hohen Anteil an Löchern und Anodenseite |
| n-dotiertes Indiumphosphid | 350 nm | \(N_D=5\times10^{17}\,\mathrm{cm^{-3}}\) | Bereich mit einem Überschuss an Elektronen, der den Übergangsbereich bildet. |
| n⁺-dotiertes Indiumphosphid | 100 Nanometer | \(N_D=1\times10^{18}\,\mathrm{cm^{-3}}\) | Stärker dotierte Schicht auf der Kathodenseite |
Der in der Strukturansicht dargestellte 200-nm-Luftbereich ist Teil der Anzeige und der umgebenden optischen Geometrie und nicht des elektrischen Halbleiterschichtstapels. Seine Dicke wird daher von der Dicke des elektrischen Bauelements von 700 nm ausgeschlossen.
3.2 Dotierung und integrierte elektrostatische Elemente
Öffnen Sie „Elektrische Eigenschaften“ → „Dotierung / Ionen“, um die Konzentrationen der Donor- und Akzeptor-Atome (??) anzuzeigen. Die Netto-Dotierung ändert ihr Vorzeichen an der Grenzfläche zwischen p-InP und n-InP, etwa 250 nm unterhalb der Oberfläche des elektrisch aktiven Halbleiterschichtstapels. Die zweite Änderung, zwischen n-InP und n⁺-InP, verändert die Elektronendichte in Richtung der Kathode, führt aber nicht zur Bildung einer weiteren PN-Übergangszone.
Im Gleichgewichtszustand neigen Elektronen dazu, von der n-dotierten Region in die p-dotierte Region zu diffundieren, während Löcher sich in die entgegengesetzte Richtung diffundieren. Die freigelegten, ionisierten Dotierstoffe erzeugen ein Feld, das einer weiteren Diffusion entgegenwirkt. Die resultierende Verarmungszone und die interne Potentialbarriere bilden die Grundlage für die Gleichrichterwirkung.
\[ \frac{d}{dy}\left(\varepsilon\frac{d\phi}{dy}\right)=-\rho(y), \qquad \rho(y)=q\,[p-n+N_D^+-N_A^-]. \]Hierbei ist \(\phi\) das elektrostatische Potential, \(\varepsilon\) die dielektrische Permittivität und \(\rho\) die lokale Ladungsdichte. Die ionisierten Dotierstoffe und die beweglichen Ladungsträger bestimmen das lokale elektrische Feld, während die angelegte Spannung die effektive Barriere für die Ladungsträgerinjektion verändert.
3.3 Elektrische Kontakte
Das Modell verwendet eine ohmsche Grenzschicht, in der Löcher überwiegen, an der Anode und eine ohmsche Grenzschicht, in der Elektronen überwiegen, an der Kathode. Diese Randbedingungen sind geeignet, um das Verhalten des Übergangs zu untersuchen, ohne eine zusätzliche Injektionsbarriere an den Kontakten einzuführen. Bei hergestellten Bauelementen bestimmen die Kontaktmetallisierung und die Verarbeitung, wie nahe die Kontakte einem idealen ohmschen Verhalten kommen.
4. Materialmodelle und Transportparameter
4.1 Lua-Materialmodelle: Woher stammen die InP-Parameter?
Eine der wichtigsten Eigenschaften dieses InP-Dioden-Beispiels ist, dass die elektrischen Materialparameter mithilfe von Lua-Materialmodellen definiert werden. Dies schafft eine klare Verbindung zwischen den physikalischen Eigenschaften eines Halbleiters und der Drift-Diffusions-Lösung. Die p-InP-, n-InP- und n⁺-InP-Bereiche beziehen sich alle auf dieselbe Materialdefinition, Inorganic/InP/InP.lua. Ihr unterschiedliches elektrisches Verhalten ergibt sich hauptsächlich aus dem Dotierungsprofil und nicht aus der Verwendung von drei nicht miteinander in Beziehung stehenden Sätzen von Materialkonstanten.
Um das Modell zu untersuchen, öffnen Sie „Gerätestruktur“ → „Elektrische Parameter“ (der Editor für elektrische Materialien) und wählen Sie eine der InP-Schichten aus. Die Registerkarte „Code-Editor“ zeigt das Lua-Skript für das Material an, das dieser Region zugeordnet ist; die Registerkarte „Parameter“ zeigt die berechneten Materialeigenschaften an. Diese beiden Ansichten sollten zusammen verwendet werden. Der Code definiert, wie die Eigenschaften ermittelt werden, während die berechneten Werte zeigen, welche Daten der elektrische Solver tatsächlich erhält.
Inorganic/InP/InP.lua verweist.Die Abbildungen in ?? und ?? zeigen, wie auf den Code-Editor zugegriffen wird. Die Schichtkonfiguration der Simulation ist die maßgebliche Grundlage für das zugewiesene Materialmodell: In diesem Beispiel ist jeder aktive Bereich mit der InP-Lua-Datei verknüpft. Die Dateinamen der Bilder werden aus der Sequenz der Bildschirmaufnahmen übernommen und bestimmen nicht, welcher Halbleiter simuliert wird.
In der bereitgestellten Simulation sind die db_overrides-Schalter auf false eingestellt, um die wichtigsten Transport- und Bandparameter zu steuern. OghmaNano ist daher so konfiguriert, dass es die entsprechenden Werte aus der Materialdatenbank/dem Lua-Modell bezieht, anstatt die in der JSON-Datei an anderer Stelle gespeicherten Standardwerte zu verwenden. Dies ist wichtig, wenn die Konfiguration überprüft wird: Die im Reiter „Parameter“ angezeigten effektiven Werte sind aussagekräftiger als die in einem inaktiven Überschreibungsfeld gespeicherten Standardwerte.
4.2 Bandstruktur, Zustandsdichte und Beweglichkeit
Die ersten Größen, die im Tab „Ausgewertete Parameter“ überprüft werden sollten, sind die Bandlücke \(E_g\) und die Elektronenaffinität \(\chi\). Die Bandlücke bestimmt die Trennung zwischen dem Leitungsband und dem Valenzband, während die Elektronenaffinität die Position des Randes des Leitungsbandes relativ zu einer Vakuumenergie-Referenz angibt. Für InP beträgt die Bandlücke bei Raumtemperatur ungefähr \(1.34\,\mathrm{eV}\). Eine nützliche erste Konsistenzprüfung besteht darin, zu überprüfen, ob die Trennung von \(E_c\) und \(E_v\) im resultierenden Banddiagramm mit der ausgewerteten Bandlücke übereinstimmt.
Die effektiven Zustandsdichten der Leitungs- und Valenzbänder, \(N_c\) bzw. \(N_v\), beeinflussen die Beziehung zwischen Trägerkonzentration und der Position des Quasi-Fermi-Niveaus. Sie tragen daher zu den Gleichgewichtskonzentrationen von Elektronen und Löchern sowie zur Injektionsantwort bei Durchlassvorspannung bei. Die relative Permittivität \(\varepsilon_r\) wird in die Poisson-Gleichung eingesetzt und beeinflusst das elektrostatische Potenzial und die Breite der Verarmungszone. Diese Größen sollten als ein konsistentes Set betrachtet werden, nicht als unabhängige, numerisch einstellbare Parameter.
Die Beweglichkeit von Elektronen und Löchern, \(\mu_n\) bzw. \(\mu_p\), bestimmt die Driftantwort der Ladungsträger auf ein elektrisches Feld; die Diffusion hängt zusätzlich von den Gradienten der Ladungsträgerdichte ab. In dem eindimensionalen Modell werden die Stromdichten wie folgt dargestellt:
\[ J_n=q\mu_n nE+qD_n\frac{dn}{dy},\qquad J_p=q\mu_p pE-qD_p\frac{dp}{dy}. \]Hierbei sind \(n\) und \(p\) die lokalen Elektron- bzw. Lochdichten, \(E\) das elektrische Feld und \(D_n\) und \(D_p\) die Diffusionskoeffizienten. Die beiden Stromkomponenten können innerhalb der Diode erheblich variieren, selbst wenn der gesamte stationäre Strom annähernd positionsunabhängig ist. Aus diesem Grund werden in den folgenden Abschnitten \(J_n\) und \(J_p\) zusammen mit den Trägerdichteprofilen untersucht.
4.3 Wie man die Materialeinstellungen überprüft und vergleicht
Im elektrischen Editor wählen Sie p-InP aus und untersuchen den zugehörigen Lua-Code und die berechneten Parameter. Wechseln Sie dann zu n-InP und schließlich zu n⁺-InP. Vergleichen Sie das zugewiesene Materialmodell, die Bandlücke, die Elektronenaffinität, die Beweglichkeiten, die effektiven Zustandsdichten, die Permittivität und die Rekombinationskoeffizienten. Für diese Homoverbindung wird in allen drei Bereichen dasselbe zugrunde liegende InP-Modell verwendet; die Hauptunterschiede zwischen den Bereichen bestehen in ihrer Akzeptor- oder Donorkonzentration.
Es lohnt sich, zwei Konzepte voneinander zu trennen: Das Materialmodell beschreibt die Eigenschaften von InP, während das Dotierungsprofil festlegt, welche Ladungsträger lokal vorherrschend sind und wo sich der Übergangsbereich befindet. Diese Unterscheidung wird besonders nützlich, wenn man von einer PN-Homojunction zu einer InP-basierten Heterojunction übergeht, wo Unterschiede zwischen den Materialmodellen selbst zu Bandversatz führen können. In dieser Diode spiegelt die in den Momentaufnahmen dargestellte Biegung des Leitungs- und Valenzbandes hauptsächlich die elektrostatischen Eigenschaften über den p/n-Übergang hinweg wider.
Um eine reproduzierbare Studie zu gewährleisten, sollte das Lua-Materialmodell bei Änderungen der Dotierkonzentrationen konstant gehalten werden. Umgekehrt sollte bei der Untersuchung der Empfindlichkeit des Ergebnisses gegenüber einer Beweglichkeits- oder Rekombinationskonstante dieser Materialparameter gezielt geändert und die resultierende J–V-Kurve, die Trägerdichte und die internen Stromprofile verglichen werden. Der Parameter-Editor und die Lua-Code-Ansichten stellen die Verbindung zwischen den physikalischen Eingangsannahmen und der simulierten Ausgabe her.
4.4 Rekombinationsparameter
Die InP-Bereiche umfassen die Shockley-Read-Hall-(SRH-)Rekombination und einen Rekombinationskanal für freie Ladungsträger, wie sie im untersuchten Materialmodell beschrieben werden. Die SRH-Rekombination beschreibt die defektunterstützte Einfangung von Elektronen und Löchern; die Rekombination freier Ladungsträger ist besonders relevant bei InP mit direkter Bandlücke, da Elektronen-Loch-Übergänge auf radiative Weise stattfinden können. In der vorliegenden Übung sind beide Prozesse als Ladungsträgerverlustprozesse nützlich, die den elektrischen Strom beeinflussen.
\[ R_{\mathrm{rad}}=B\left(np-n_i^2\right), \qquad R_{\mathrm{SRH}}=\frac{np-n_i^2}{\tau_p(n+n_1)+\tau_n(p+p_1)}. \]Der Koeffizient \(B\), die Lebensdauer der Ladungsträger in den Fallen \(\tau_n\) und \(\tau_p\) sowie die zusätzlichen Ladungsträgerdichten \(n_1\) und \(p_1\) bestimmen die Größenordnung jedes Prozesses. Der wichtigste Aspekt in diesem Tutorial ist, wie sich die räumliche Verteilung der Rekombination verändert, wenn durch eine Vorspannung mehr Elektronen und Löcher in dieselbe Region getrieben werden.
4.5 Netzgenerierung und numerische Lösung
Für die elektrische Berechnung wird ein eindimensionales Gitter mit 150 Punkten über eine Länge von 700 nm verwendet, was einer durchschnittlichen Punktabstand von etwa 4,7 nm entspricht. Die gekoppelten nichtlinearen Gleichungen werden mit der Newton-Iteration und dem UMFPACK-Sparse-Linearlösungsverfahren gelöst. Dieses Gitter stellt einen nützlichen Ausgangspunkt für die Untersuchung des Übergangsbereichs und der quasi-neutralen Bereiche dar.
Die Ausgabe enthält Informationen zur Konvergenz und gerätespezifische Profile, die von der Vorspannung abhängen. Eine gute erste Überprüfung besteht darin, auf eine glatte J–V-Kurve zu achten und anschließend zu bestätigen, dass die Elektronen- und Lochströme den allgemeinen Anforderungen der Stromerhaltung im stationären Zustand entsprechen. Der Gittereditor wird in ?? dargestellt.
5. Durchführung der Simulation und Auswertung der J–V-Kennlinie
Kehren Sie zum Hauptfenster zurück und wählen Sie „Simulation ausführen“ aus. OghmaNano durchläuft den eingestellten Spannungsbereich und speichert die konvergierte Lösung an jedem Punkt. Die gespeicherte Beispieldatei verwendet einen Spannungsbereich von 0–3 V mit Schritten von 25 mV. Das Terminal protokolliert den Fortschritt des Solvers während des Durchlaufs, während die Registerkarte „Ausgabe“ Zugriff auf die resultierenden Daten bietet.
Sobald die Berechnung abgeschlossen ist, werden die Ergebnisse im Ausgabetab (??) angezeigt. Öffnen Sie jv.csv, um die Stromdichte in Abhängigkeit von der angelegten Spannung darzustellen. Die Screenshots enthalten sowohl eine lineare J–V-Darstellung (??) als auch eine logarithmische Stromdarstellung (??).
5.1 Das Wesen verstehen
Bei niedriger angelegter Spannung schränkt die in der Diode vorhandene Sperrschicht die Injektion von Ladungsträgern stark ein. Eine Erhöhung der Durchlassspannung verringert diese Sperrschicht, wodurch Elektronen und Löcher leichter in die gegenüberliegende Seite gelangen und der Gesamtstrom ansteigt. Im mittleren Bereich ist der Strom oft annähernd exponentiell von der Spannung abhängig; bei höheren Strömen können Transport- und Widerstandsverluste die Steigung verändern.
\[ J\approx J_0\left[\exp\!\left(\frac{qV}{nk_{\mathrm B}T}\right)-1\right]. \]Dieser Ausdruck für die ideale Diode bietet eine nützliche Referenz für die Analyse eines annähernd exponentiellen Intervalls, wird aber nicht durch das Drift-Diffusions-Modell vorgegeben. Der Idealitätsfaktor \(n\) und der Vorfaktor \(J_0\) sollten nur über ein geeignetes Spannungsintervall ermittelt werden; sie können sich in Abhängigkeit von der Rekombination, dem Reihenwiderstand und dem Kontaktverhalten ändern.
Die logarithmische Darstellung hilft dabei, den Beginn eines signifikanten Stromflusses und Änderungen in der Steigung zu erkennen. Die lineare Darstellung erleichtert die Interpretation der Stromstärke und der Form des Stromverlaufs bei hohen Strömen. Für ein möglichst umfassendes Bild des Diodenverhaltens sollten beide Darstellungen zusammen mit den internen Band- und Trägerdichteprofilen betrachtet werden.
6. Interne Funktionsweise: Wie die Diode funktioniert.
Das Verzeichnis „Snapshots“ enthält räumliche Profile, die an ausgewählten Messpunkten während der J–V-Messung gespeichert wurden. Die bereitgestellten Screenshots zeigen vier Positionen innerhalb dieser Messung. Sie veranschaulichen den Übergang von einer geringen Injektion zu einer stärkeren Durchlassvorspannung. Die vier Bilder in jeder Gruppe beziehen sich auf dieselben ausgewählten Positionen des Schiebereglers, wodurch ein direkter Vergleich der verschiedenen physikalischen Größen ermöglicht wird.
6.1 Leitungs- und Valenzbänder sowie Quasi-Fermi-Niveaus
Öffnen Sie die Ansicht für Momentaufnahmen und wählen Sie Ec.csv, Ev.csv, Fn.csv und Fp.csv aus. Dies sind die Energiegrenze des Leitungsbandes, die Energiegrenze des Valenzbandes, das Quasi-Fermi-Niveau der Elektronen und das Quasi-Fermi-Niveau der Löcher. Ihre räumliche Trennung und ihre Steigung geben Aufschluss über die Potentialbarriere, die Reaktion auf die angelegte Vorspannung und die Abweichung vom Gleichgewichtszustand.
Bei geringer Injektion biegen sich die Bandkanten aufgrund des vorhandenen elektrostatischen Feldes über den Übergang. Bei zunehmender Durchlassvorspannung wird die effektive Barrierehöhe verringert, und die quasi-Fermi-Niveaus für Elektronen und Löcher trennen sich, während Ladungsträger durch die Komponente geleitet werden. Der entscheidende Punkt ist, dass InP auf beiden Seiten des Übergangs das gleiche Halbleitermaterial bleibt: Die sich ändernden Profile spiegeln seine durch die Dotierung bestimmte Elektrostatik und die Ladungsträgerinjektion wider.
Der Vergleich zwischen ?? und ?? veranschaulicht, wie sich die terminale Vorspannung auf die interne Bandstruktur auswirkt. Insbesondere liefern die Quasi-Fermi-Niveaus eine direkte Möglichkeit, den Nichtgleichgewichtszustand zu bestimmen, der während des stationären Stromflusses erreicht wird.
6.2 Elektron- und Lochdichten
Als Nächstes werden die Diagramme Q_nfree.csv und Q_pfree.csv dargestellt, die die räumliche Verteilung von freien Elektronen und Löchern zeigen. Bei niedriger Vorspannung überwiegen die Elektronen auf der n-dotierten Seite und die Löcher auf der p-dotierten Seite. Ihre Konzentrationen nehmen in der Nähe des Übergangs stark ab, was die Verarmungszone widerspiegelt.
Wenn die Durchlassspannung zunimmt, werden Minderheitsträger in die gegenüberliegenden Bereiche injiziert. Dadurch erhöht sich die Überlappung der Elektron- und Lochpopulationen, und die Form beider Dichteprofile ändert sich. Bei höheren Injektionspegeln unterscheiden sich die Konzentrationen der freien Ladungsträger über den pn-Übergang erheblich von ihren Werten bei niedriger Durchlassspannung.
Vergleichen Sie den Punkt, an dem sich die Trägerdichte ändert, mit der Position des Dotierungsprofils und anschließend mit den Stellen, an denen eine starke Bandbiegung auftritt. Dadurch wird eine Verbindung zwischen den vom Kontakt gelieferten Ladungsträgern und der elektrostatischen Barriere hergestellt, die durch den PN-Übergang entsteht.
6.3 Stromdichten von Elektronen und Löchern
Stellen Sie die Diagramme Jn.csv und Jp.csv dar. Diese Profile zeigen, wie sich die Elektronen- und Lochbeiträge zum elektrischen Strom durch das Bauelement zusammensetzen. Die einzelnen Komponenten können je nach Position variieren, da Rekombinationsprozesse Ladungsträger zwischen den Elektronen- und Lochpopulationen übertragen. Der gesamte Strom im stationären Zustand sollte jedoch ungefähr konstant bleiben, wenn sich keine Netto-Ladung ansammelt.
Bei sehr geringem Strom können sich die Drift- und Diffusionsanteile gegenseitig aufheben, wodurch einzelne Stromverlaufskurven unregelmäßig oder stark strukturiert erscheinen. Die bei niedriger Vorspannung erstellten Screenshots zeigen dieses Verhalten deutlicher als die bei höherem Strom erstellten Aufnahmen. Es ist sinnvoll, die Verläufe in einer geeigneten Skala zu betrachten und den Gesamtstrom mit dem Endwert der J–V-Kennlinie zu vergleichen, bevor man feine numerische Merkmale interpretiert.
6.4 Rekombinationsprofil
Die Rekombination bestimmt, wo die injizierten Elektronen und Löcher aus ihren jeweiligen Trägerpopulationen entfernt werden. Öffnen Sie die Momentaufnahme-Anzeige mit R_nfree_to_pfree.csv und R_ss_srh.csv, um die direkten, freien Rekombinations- und die defektvermittelten Kanäle zu vergleichen. Beide Prozesse hängen von der Verfügbarkeit von Elektronen und Löchern ab, sodass sich ihre räumliche Verteilung verändert, wenn sich die injizierten Trägerpopulationen ändern.
Unmittelbar nach Beginn der Leitfähigkeit kann die Rekombination in der Nähe des Übergangs konzentriert auftreten, wo sich die Konzentrationen von Elektronen und Löchern überlappen. Bei stärkerer Vorwärtsvorspannung breitet sich diese Überlappung auf einen größeren Bereich des Halbleiters aus, und die Rekombinationsprofile können breiter werden. Vergleichen Sie die Position und Breite der Rekombinationspeaks mit den Trägerdichten bei derselben Vorspannung.
Beim Vergleich der beiden Kanäle sollten sowohl der Spitzenwert als auch das räumliche Integral berücksichtigt werden. Ein ausgeprägter lokaler Spitzenwert muss nicht unbedingt die gesamte Rekombinationsrate bestimmen, wenn ein anderer Kanal über ein viel größeres Volumen aktiv ist. Dies ist besonders wichtig, wenn die Rekombinationsprofile in Beziehung zum gesamten elektrischen Strom gesetzt werden.
6.5 Profile der resultierenden Elektron- und Lochdichte
Die mit R_n.csv und R_p.csv gekennzeichneten Abbildungen liefern weitere Informationen über die in die Kontinuitätsgleichungen für Elektronen und Löcher einfließenden Raten. Sie zeigen, wo das Gleichgewicht der lokalen Erzeugungs-, Einfangs- und Rekombinationsgrößen am stärksten ausgeprägt ist. Interpretieren Sie die Vorzeichen entsprechend der in OghmaNano verwendeten Darstellungskonvention.
Zusammenfassend lässt sich sagen, dass die Bandstruktur, die Trägerdichte, die Ströme und die zeitlichen Verläufe eine vollständige Erklärung für das Verhalten der Diode liefern: Die interne Barriere steuert die Injektion, die angelegte Vorspannung verändert die interne Elektrostatik, und Transport- sowie Rekombinationsprozesse bestimmen den resultierenden stationären Strom.
7. Übungen und Interpretation
8. Interpretation des Modells und weitere Analyse
Dieses Beispiel verwendet ein eindimensionales, stationäres Drift-Diffusions-Modell für InP, einen kompakten p/n/n⁺-Halbleiterstapel und idealisierte ohmsche Anschlüsse. Es eignet sich gut, um zu veranschaulichen, wie elektrische Gleichrichtung aus der Übergangssperrschicht entsteht und wie dieses Verhalten mit der räumlichen Ladungsträgerverteilung zusammenhängt. Umfangreichere Geräteuntersuchungen können gemessene Kontaktcharakteristiken, verfeinerte Materialparameter, Temperaturvariation und zusätzliche geometrische Details einbeziehen.
Die Materialeigenschaften von InP werden aus dem ausgewählten Lua-Modell entnommen. Für quantitative Studien sollten die ermittelten Werte für Bandlücke, Elektronenaffinität, Beweglichkeit und Rekombinationskoeffizienten im Materialeditor oder in den Ausgabedaten überprüft werden. Überprüfen Sie die Erhaltung des Stroms und die Solver-Residuen und wiederholen Sie repräsentative Arbeitspunkte mit einem feineren Netz, um die numerische Stabilität sicherzustellen.
Der in dem bereitgestellten Beispiel verwendete Standard-Messbereich erstreckt sich von 0 bis 3 V. Ein aufschlussreicher nächster Schritt wäre, einen moderaten Bereich für die Sperrspannung einzubeziehen, um einen direkten Vergleich der Durchlassleitung und der Sperrleckströme zu ermöglichen. Der physikalisch relevante Bereich für die Sperrspannung und die Leckcharakteristiken hängen vom Bauelement und von den in der Simulation berücksichtigten Mechanismen ab.
Um weitere Erkenntnisse zu gewinnen, sollten Sie die Dotierkonzentrationen nacheinander variieren. Eine Erhöhung der nominalen Dotierkonzentration verringert im Allgemeinen die Verarmungszone und verändert die elektrostatischen Gleichgewichtseigenschaften; eine Änderung der Rekombinationskoeffizienten beeinflusst den injektionsabhängigen Strom. Dies sind nützliche, kontrollierte Vergleichsversuche, da die Schichtmaterialien und die Geometrie unverändert bleiben können.
9. Nächste Schritte
Die InP-PN-Diode stellt einen nützlichen Ausgangspunkt für komplexere, auf InP basierende Halbleiterstrukturen dar. Eine naheliegende Weiterentwicklung besteht darin, eine schwach dotierte Region einzuführen und das elektrostatische und Transportverhalten einer PN- und einer PIN-Diode zu vergleichen. Eine weitere Möglichkeit besteht darin, Heterostrukturen mit InGaAs/InP zu untersuchen, wobei Unterschiede in der Zusammensetzung zusätzliche Bandversätze und ein neues Verhalten bei der Ladungsträgerinjektion verursachen.