Nonlocal Electrostatic Origin of Schottky-Barrier Variability in 2D Contacts

Technisches Resümee: Nichtlokaler elektrostatischer Ursprung der Variabilität der Schottky-Barriere in 2D-Kontakten

Problemstellung
Elektrische Kontakte begrenzen häufig die Leistungsfähigkeit von atomar dünnen Halbleiterbauelementen. Eine anhaltende Herausforderung auf diesem Gebiet ist die weite Variabilität der berichteten Schottky-Barrierenhöhen (SBH) für identische Metall/2D-Halbleiter-Paare. Beispielsweise liegen die berichteten SBH-Werte für Au–MoS₂-Kontakte zwischen 0,05 und 0,95 eV und für Ti–MoS₂ zwischen 0,05 und 0,46 eV. Konventionelle theoretische Rahmenwerke behandeln die SBH als eine lokale Grenzflächeneigenschaft und führen Abweichungen vom Schottky–Mott-Limit auf Grenzflächenzustände und Fermi-Level-Pinning zurück. In dieser Sichtweise wird erwartet, dass mikroskopische Unordnung durch Grenzflächen-Ladungsreservoirs gemittelt oder durch den „Pinch-off“-Effekt umliegender Hochbarriere-Regionen begrenzt wird. Die experimentell beobachtete „riesige Streuung“ der SBH-Werte lässt sich jedoch nicht mit einem Modell vereinbaren, bei dem die Barriere ausschließlich durch lokale Grenzflächeneigenschaften bestimmt wird.

Methodik
Die Autoren schlagen ein nichtlokales elektrostatisches Modell vor, bei dem die effektive SBH in top-kontaktierten 2D-Bauelementen als eine Kanten-Observierbare behandelt wird, die durch ferne elektrostatische Störungen getrieben wird, anstatt als räumliche Verteilung lokaler Barrieren.

1. Theoretischer Rahmen:

   • Die Studie betrachtet eine Top-Kontakt-Geometrie mit einem bedeckten Bereich ($x < 0$), einem unbedeckten Kanal ($x > 0$) und einer Kontaktkante bei $x = 0$.
   • Die laterale SBH ($\Phi^{\text{lat}}_B$) wird als Leitungsband-Stufe über die Kante definiert. Das Modell postuliert, dass die Barriereverschiebung durch die Änderung der Kanten-Linien-Anordnung ($\Delta J$) gesteuert wird, welche linear an das elektrostatische Potenzial an der Kante, $u(0)$, gekoppelt ist.
   • Unter Verwendung eines quadratischen Freie-Energie-Funktionals für das elektrostatische Potenzial leiten die Autoren einen Antwort-Kernel ab, der eine ferne Defektquelle $s_d(x)$ auf das Kantenpotenzial abbildet.
   • Die resultierende Barriereverschiebung aufgrund eines Punktdefekts in einem Abstand $x_d$ ist gegeben durch:
     $$\Delta \Phi^{\text{lat}}_B (x_d) = \chi_e S_d \frac{\lambda_l}{\kappa + K_e \lambda_l} e^{-x_d/\lambda_l}$$
     wobei $\lambda_l$ die Propagationslänge ist, die durch die in-plane elektrostatische Steifigkeit ($\kappa$) des Kanals und die Rückstellkapazität ($C_{ch}$) bestimmt wird, und die Amplitude durch kontaktspezifische Parameter (Kanten-Suszeptibilität $\chi_e$ und Kanten-Steifigkeit $K_e$) bestimmt wird.

2. Komputationelle Validierung:

   • Das Modell wird mittels First-Principles-Berechnungen mittels Dichtefunktionaltheorie (DFT) in Kombination mit der Nichtgleichgewichts-Green-Funktions-Methode (NEGF) validiert.
   • Untersuchte Systeme sind Ti–MoS₂- und Au–MoS₂-Top-Kontakte.
   • Schwefelvakanzen werden im MoS₂-Kanal in variierenden Abständen ($0d$ bis $10d$) zur Kontaktkante eingeführt.
   • Effektive SBH-Werte werden aus der lokalen Zustandsdichte (LDOS) nahe der Kontaktkante extrahiert, spezifisch durch Messung der Leitungsband-Stufe.

Zentrale Beiträge und Ergebnisse

• Nichtlokaler Mechanismus: Die Arbeit stellt fest, dass Defekte, die fernab der Kontaktkante liegen, die gemessene SBH signifikant durch nichtlokale elektrostatische Kopplung modulieren können. Der Einfluss eines Defekts nimmt exponentiell mit dem Abstand ab, gesteuert durch eine charakteristische Längenskala $\lambda_l$.
• Trennung der Rollen: Das Modell zeigt eine klare Aufteilung der Rollen auf:
  • Der Halbleiterkanal bestimmt den Reichweite des nichtlokalen Effekts ($\lambda_l$), welche von der in-plane elektrostatischen Ausbreitung und Abschirmung abhängt.
  • Der Metallkontakt bestimmt die Amplitude der SBH-Verschiebung, gesteuert durch die Suszeptibilität der Kante gegenüber Potenzialänderungen und die Effizienz der Umwandlung von Kantenpotenzialverschiebungen in eine Linien-Anordnungsmodulation.
• Quantitative Übereinstimmung:
  • DFT-NEGF-Berechnungen für sowohl Ti–MoS₂ als auch Au–MoS₂ zeigen, dass die laterale SBH dem vorhergesagten exponentiellen Abfall mit dem Defektabstand folgt.
  • Eine gemeinsame Interaktionslänge von $\lambda_l \approx 1,1$ nm passt zu den Daten für beide Metalle, was bestätigt, dass diese Längenskala primär eine Eigenschaft der MoS₂-Kanal-Elektrostatik und nicht des spezifischen Metalls ist.
  • Die berechneten SBH-Variationen decken einen Bereich ab, der einen wesentlichen Teil der experimentell berichteten Variabilität für Ti–MoS₂ erklärt und einen signifikanten Teil der Streuung für Au–MoS₂ abdeckt.
• Vereinheitlichte Erklärung: Die Ergebnisse liefern eine einheitliche Erklärung für die langjährige Variabilität der extrahierten SBH-Werte, indem sie zeigen, dass nominell ähnliche Kontakte je nach räumlicher Verteilung von Defekten relativ zur Kontaktkante sehr unterschiedliche Barrierenhöhen liefern können.

Bedeutung
Die Arbeit beansprucht, eine große Inkonsistenz in der 2D-Kontaktliteratur gelöst zu haben, indem sie zeigt, dass die gemessene Schottky-Barriere keine rein lokale Grenzflächeneigenschaft ist, sondern eine nichtlokale, kantengesteuerte Größe. Durch die Identifizierung der Kantenumgebung und der räumlichen Verteilung nahegelegener Defekte als Schlüsseldesignvariablen legt die Arbeit nahe, dass der Ladungsträgerinjektion in 2D-Bauelementen durch mesoskopische Unordnungsumgebungen und nicht nur durch atomare Grenzflächenstrukturen unterliegt. Dieser Rahmen geht über klassische inhomogene-Schottky- und Pinch-off-Modelle hinaus und bietet einen Mechanismus für die „riesige Streuung“ der SBH-Werte sowie etabliert nichtlokale Elektrostatik, vermittelt durch kantennahe Defekte, als einen fundamentalen Faktor für die 2D-Kontaktleistung.