MOS-Kondensator (NMOS): Elektrostatik, Verarmungszone und Bandbiegung
1. Einleitung
Ein MOSFET (Metall-Oxid-Halbleiter-Feldeffekttransistor) funktioniert, indem eine Gate-Spannung verwendet wird, um die Trägerdichte an einer Halbleiteroberfläche elektrostatisch zu steuern. Das einfachste „Atom“ der MOSFET-Physik ist der MOS-Kondensator: ein Metall-Gate, eine isolierende Oxidschicht und ein dotiertes Halbleitermaterial. Wenn Sie den MOS-Kondensator verstehen, verstehen Sie auch die Elektrostatik, die einen MOSFET-Schalter ermöglicht (siehe ??).
Ein NMOS-Bauelement wird als „NMOS“ bezeichnet, da es in einer Transistorgeometrie, wenn es eingeschaltet ist, einen n-leitenden Inversionskanal bildet – die leitfähigen Kanalträger sind Elektronen. In der Praxis wird dies üblicherweise mit einem p-leitenden Siliziumkörper (akzeptordotiertes Si) und einer ausreichend positiven Gate-Spannung realisiert. Schon bevor Quelle und Drain hinzugefügt werden, ist die gleiche Gate-Elektrostatik im MOS-Kondensator sichtbar: Eine positive Gate-Spannung zieht Elektronen in Richtung der Si/SiO2-Grenzfläche (Oberflächenakkumulation), während eine negative Gate-Spannung Löcher abstößt und eine Region mit festen, ionisierten Akzeptoren (Verarmungszone) hinterlässt, wie in ?? dargestellt.
Ein PMOS-Bauelement stellt den komplementären Fall dar: Es bildet einen p-Kanal (Löcher als leitfähige Ladungsträger), wobei üblicherweise ein n-dotiertes Substrat bei negativer Gate-Spannung verwendet wird. NMOS- und PMOS-Bauelemente bilden zusammen das grundlegende komplementäre Paar, das in der CMOS-Logikschaltung eingesetzt wird.
In diesem Tutorial erstellen und untersuchen wir eine zweidimensionale NMOS-Kondensatorstruktur in OghmaNano. Wir konzentrieren uns gezielt auf die Anreicherungs- und Verarmungsbereiche, in denen die Physik hauptsächlich von elektrostatischen Effekten bestimmt wird und die Ergebnisse leicht zu interpretieren sind. Die starke Inversion (ein echter leitfähiger Kanal) wird auf Teil B und spätere Tutorials verschoben.
2. OghmaNano starten
Starten Sie OghmaNano über das Windows-Startmenü. Das Hauptfenster von OghmaNano wird wie in ?? dargestellt angezeigt. Klicken Sie auf „Neue Simulation“. Dadurch wird die Bibliothek der verfügbaren Gerätetypen geöffnet, wie in ?? dargestellt. Doppelklicken Sie auf „Si Demos“, um die Silizium-Beispielsammlung zu öffnen, und wählen Sie dann das nMOS-Kondensator-Beispiel aus, wie in ?? dargestellt. Wenn Sie dazu aufgefordert werden, speichern Sie die Simulation in einem lokalen Ordner mit Schreibzugriff. Nach dem Speichern wird das Hauptfenster der Simulation geöffnet (siehe ??).
💡 Tipp: Speichern Sie Simulationen auf einem lokalen Laufwerk (z. B. C:\), um die beste Leistung zu erzielen. Ordner, die auf einem Netzwerk oder in der Cloud gespeichert sind, können wiederholte Berechnungen verlangsamen.
3. Überprüfen Sie die Struktur des Geräts (Schichten und Materialien).
Nachdem das Projekt erstellt wurde, öffnet sich das Hauptfenster der Simulation, das eine bewusst vereinfachte NMOS-MOS-Kondensatorstruktur anzeigt: ein Metallgitter, ein SiO2-Gate-Oxid und ein p-dotiertes Siliziumsubstrat mit einem Body-Kontakt. Es wird als NMOS-Kondensator bezeichnet, da der Halbleiterkörper p-dotiert ist und eine positive Gate-Vorspannung die Oberflächenelektrostatik in Richtung Elektronenanreicherung lenkt. Diese Struktur stellt die typische NMOS-elektrostatische Konfiguration dar und wird häufig verwendet, um die Bandbiegung, Anreicherung und Verarmung an der Si/SiO2-Grenzfläche zu untersuchen.
In praktischen Bauelementen ist der MOS-Kondensator oft in komplexere Strukturen eingebettet oder mit zusätzlichen Kontakten kombiniert. Hier werden diese Elemente jedoch bewusst weggelassen, damit der Fokus vollständig auf der vertikalen Elektrostatik des NMOS-Kondensators selbst liegt. Durch die Eliminierung der lateralen Ladungsträgerbewegung und der geometrischen Komplexität kann der physikalische Ursprung der Ladungsverteilung und der Variation des elektrostatischen Potenzials im Halbleiter direkt und eindeutig untersucht werden.
Die Schichtstruktur kann im Schichteditor überprüft werden. Um die Darstellung zu vereinfachen, wird in diesem Tutorial eine etwas größere Oxidschichtdicke gewählt, als es in modernen Siliziumtechnologien üblich ist, sodass der resultierende Potentialabfall, die Akkumulationsschicht und die Verarmungszone in den Ausgabediagrammen deutlich sichtbar sind. In modernen NMOS-Kondensatoren liegen die Oxidschichtdicken typischerweise im Bereich von 10 nm oder darunter. Nachdem Sie das Tutorial abgeschlossen haben, werden Sie ermutigt, die Schichtdicken anzupassen und zu untersuchen, wie sich die elektrostatische Reaktion ändert, wenn sich die Struktur den realistischen Abmessungen eines Bauelements nähert.
4. Überprüfen Sie das elektrische Band und das zweidimensionale Netz.
Der größte Teil des Arbeitsablaufs für NMOS-Kondensatoren findet im Bereich „Elektrisch“ statt. Hier können Sie auswählen, welches elektrische Modell aktiv ist, Parameter bearbeiten, die Dotierung überprüfen und das Netz konfigurieren. Durch Klicken auf „Elektrisches Netz“ wird der Netzeditor geöffnet (siehe ??).
In diesem Tutorial wird ein zweidimensionales Netz verwendet: Das Bauelement wird in zwei Dimensionen aufgelöst, um die Bandbiegung und die Verarmungszone in der Nähe der Grenzfläche zu erfassen, während gleichzeitig die Berechnungszeit kurz gehalten wird. In dem gezeigten Beispiel gibt es nur wenige Punkte entlang einer Richtung (lateral) und eine dichtere Abtastung durch die Oxidschicht und die Siliziumschicht (vertikal), wo die elektrostatischen Gradienten auftreten. Sie können die Anzahl der Netzpunkte erhöhen, aber die Rechenkosten steigen in etwa proportional zur Anzahl der Unbekannten (und bei nichtlinearen Problemen zur Anzahl der Iterationen pro Arbeitspunkt).
5. Editor für elektrische Parameter (Oxid im Vergleich zu Silizium)
Wenn Sie auf „Elektrische Parameter“ klicken, wird der Editor für elektrische Parameter geöffnet. Dieser Editor ist nach Materialien geordnet: Auf der Registerkarte „SiO2“ finden Sie die Parameter für das Oxid, während auf der Registerkarte „Si“ die Parameter für Silizium angezeigt werden. Für die Siliziumschicht verwenden wir eine physikalisch sinnvolle Parametersammlung: eine relative Permittivität von etwa 11,7, Standardwerte für Silizium in Bezug auf Elektronenaffinität und Bandlücke, endliche Elektronen- und Lochmobilitäten sowie realistische effektive Zustandsdichten. Wichtig ist, dass wir die Fermi-Dirac-Statistik anstelle der Maxwell-Boltzmann-Statistik verwenden, da bei angelegter Gate-Spannung die Trägerdichte in der Nähe der Si/SiO2-Grenzfläche so hoch werden kann, dass Degenerierungseffekte eine Rolle spielen.
Im Gegensatz dazu ist bei der Oxidschicht die relative Permittivität der einzige Parameter, der in diesem Tutorial eine zentrale Rolle spielt. Das SiO2 wird als idealer Isolator behandelt: Es ermöglicht ein elektrisches Feld und einen Spannungsabfall, unterstützt aber keine bewegliche Ladung. Infolgedessen wird das elektrostatische Verhalten des NMOS-Kondensators dadurch bestimmt, wie sich die angelegte Gate-Spannung zwischen der Oxidschicht und dem Silizium aufteilt, was ausschließlich durch die Dicke und die Permittivität der Oxidschicht bestimmt wird.
Diese Wahl des Modells spiegelt sich in den Schaltflächen oben im Editor wider. In der SiO2-Schicht ist die Drift-Diffusions-Funktion deaktiviert, sodass dort nur die Poisson-Gleichung gelöst wird. Im Gegensatz dazu ist in der Siliziumschicht die Drift-Diffusions-Funktion aktiviert, damit die Umlagerung von Elektronen und Löchern als Reaktion auf die Gate-Vorspannung korrekt erfasst werden kann. Diese Trennung spiegelt die physikalischen Funktionen der beiden Materialien in einem NMOS-Kondensator wider: Das Oxid leitet das elektrische Feld weiter, während das Silizium die beweglichen Ladungsträger aufnimmt.
Möglicherweise fällt Ihnen auf, dass Parameter wie die Elektronenaffinität und die Bandlücke, die für SiO2 angegeben sind, denen von Silizium ähneln. Im Rahmen dieses Tutorials ist dies beabsichtigt und unbedenklich. Da in der Oxidschicht keine Gleichungen zur Beschreibung des Ladungstransports gelöst werden, haben diese Bandparameter keinen Einfluss auf die Lösung: Es gibt keinen Mechanismus, durch den Elektronen oder Löcher in das Modell eindringen und in die Oxidschicht gelangen könnten. Daher wirkt sich SiO2 nur über seine dielektrische Reaktion aus, nicht über seine elektronische Bandstruktur.
6. Dotierungsprofil-Editor (Dotierungsprofil des Substrats)
Öffnen Sie abschließend den Doping-Editor im Bereich „Elektrisch“ und überprüfen Sie die Dotierstoffverteilung (siehe ??). Bei einem NMOS-Kondensator ist das Siliziumsubstrat typischerweise p-dotiert (akzeptordotiert). Bei einer Null- oder negativen Gate-Vorspannung sind Löcher die vorherrschenden Ladungsträger in der Nähe der Si/SiO2-Grenzfläche. Wenn die Gate-Vorspannung positiv wird, stößt das elektrische Feld am Gate die Löcher von der Oberfläche weg, wodurch ionisierte Akzeptoren (fixierte negative Ladung) zurückbleiben und eine Verarmungszone entsteht.
Die Dotierstoffkonzentration beeinflusst maßgeblich das elektrostatische Verhalten des Kondensators. Eine höhere Dotierung vom p-Typ führt bei einem gegebenen Oberflächenpotenzial zu einer dünneren Verarmungszone, während eine geringere Dotierung eine breitere Verarmungszone und eine größere räumliche Ausdehnung der Bandbiegung zur Folge hat. In Teil B werden wir diesen Editor verwenden, um die Dotierung des Substrats gezielt zu verändern und die daraus resultierenden Änderungen in der Breite der Verarmungszone, dem elektrostatischen Potenzial und der Bandstruktur zu beobachten.
7. Führen Sie die Simulation aus (Ausgabe im Terminal).
Klicken Sie im Hauptfenster auf „Simulation starten“ (der blaue Pfeil). OghmaNano wechselt zum Tab „Terminal“ und startet den Solver. Die erste Phase einer typischen Ausführung wird in ?? dargestellt. Dort können Sie die Befehlszeile sehen, die zum Starten des Hauptsolvers verwendet wird, gefolgt vom Laden des Projekts und der Initialisierung des Modells. Die wiederholten Zeilen „Generiere…“ zeigen, dass OghmaNano Nachschlagetabellen für die Fermi-Dirac-Statistik erstellt (für Elektronen und Löcher in den von Ihnen ausgewählten Schichten). Die genaue Berechnung von Fermi-Dirac-Integralen ist aufwendig, daher berechnet der Code diese Werte einmal und verwendet sie dann während der nichtlinearen Iterationen wieder. Wenn Sie „Datei wird geladen“ sehen, bedeutet dies, dass die zwischengespeicherten Tabellen (oder andere zwischengespeicherte Zwischenergebnisse) neu geladen werden, um zu vermeiden, dass sie bei jeder Ausführung neu generiert werden müssen.
Eine zweite wichtige Phase ist die schrittweise Anpassung der Vorspannung. Bevor die Hauptberechnung beginnt, erhöht der Algorithmus die Gate-Spannung schrittweise von einem bekannten, einfachen Ausgangspunkt (in der Regel in der Nähe des Gleichgewichtszustands) auf den ersten angeforderten Vorspannungspunkt (hier -2 V). Dies ist eine Strategie, die auf der Kontinuität basiert: Newton-ähnliche Algorithmen konvergieren nur dann zuverlässig, wenn die anfängliche Schätzung ausreichend nahe an der Lösung liegt. Ein direkter Sprung zu einer hohen Vorspannung kann dazu führen, dass die anfängliche Schätzung „zu weit entfernt“ ist, was zu einer Divergenz oder einer langsamen, instabilen Konvergenz führt. Durch die schrittweise Erhöhung der Gate-Spannung liefert jede Lösung eine gute anfängliche Schätzung für die nächste.
Sobald das Gate die Startspannung erreicht hat, setzt der Solver den Gate-Sweep fort (die „J–V“-/Bias-Kurve). Der Körper des Bauelements wird auf 0 V gehalten (wie in der Terminalausgabe angegeben). Möglicherweise stellen Sie fest, dass die gedruckte Stromdichte in Vorzeichen und Größe schwankt, obwohl eine Kapazität im Wesentlichen einen konstanten Gleichstrom von Null haben sollte. Dieses Verhalten ist hier normal: Der Solver meldet sehr kleine Restströme in der Nähe des numerischen Rauschpegels der nichtlinearen Berechnung. Mit anderen Worten, die „Spitzen“ stellen keinen physikalischen Gleichstrom dar; sie sind ein Zeichen dafür, dass versucht wird, eine Größe zu berechnen, die tatsächlich nahezu 0 ist, und zwar mit einer hohen relativen Genauigkeit.
Die als f(...)= angegebene Größe ist ein nichtlineares Rest-/Fehlmaß (eine Konvergenzmetrik). Wenn dieser Wert sehr klein ist, wird davon ausgegangen, dass die elektrostatische Lösung (Potential, Ladung, Bandbiegung) konvergiert ist. Die rechts angegebene Zeit (in Millisekunden) ist die Laufzeit pro Bias-Schritt/Iterationsblock und dient zur Leistungsanalyse und zur Beurteilung der Kosten für eine Erhöhung der Netzauflösung.
f(...) stellt die Solver-Residuen dar (Konvergenzmaß).8. Ausgabedateien und aktuelle Kontaktinformationen
Sobald die Simulation abgeschlossen ist, wechseln Sie zum Reiter „Ausgabe“, um die von der Simulation erzeugten Dateien anzuzeigen (siehe ??). Sie können auf Einträge wie iv_contact_0.csv und iv_contact_1.csv doppelklicken, um die Diagramme der Kontaktstromwerte zu öffnen.
Bei einem idealen MOS-Kondensator ist der stationäre Gleichstrom im Wesentlichen null, da die Struktur rein elektrostatisch ist. Die Diagramme der Kontaktströme sind daher hauptsächlich zur Information und Vollständigkeit enthalten sowie als einfache Plausibilitätsprüfung, um sicherzustellen, dass keine ungewollten Strompfade vorhanden sind. Die Kurven sollten über den gesamten Spannungsbereich nahe bei Null bleiben; kleine Schwankungen sind normal und spiegeln lediglich die Tatsache wider, dass der tatsächliche Strom extrem gering ist.
Im nächsten Abschnitt werden wir uns von diesen zusätzlichen Ausgaben abwenden und uns auf die physikalisch relevanten Ergebnisse der Simulation konzentrieren: das elektrostatische Potenzial, die Bandkanten, die Fermi-Niveaus und die Ladungsträgerverteilungen innerhalb des NMOS-Kondensators.
💡 Was Sie bisher gelernt haben: In diesem Tutorial richten Sie eine NMOS-Kondensatorsimulation ein und untersuchen diese. Dabei werden folgende Aspekte behandelt: Bauelementstruktur, Netz, Materialparameter und Dotierung. Teil B Darauf aufbauend wenden wir eine Gate-Vorspannung an und interpretieren die Simulationsergebnisse – Bandbiegung, elektrostatische Potentiale, Trägerdichten und Oberflächenverarmung.