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

Résumé Technique : Origine Électrostatique Nonlocale de la Variabilité de la Barrière de Schottky dans les Contacts 2D

Énoncé du Problème
Les contacts électriques limitent fréquemment les performances des dispositifs à semi-conducteurs de faible épaisseur atomique. Un défi persistant dans ce domaine est la grande variabilité des valeurs de hauteur de barrière de Schottky (SBH) rapportées pour des paires métal/semi-conducteur 2D identiques. Par exemple, les SBH rapportées pour les contacts Au–MoS2 varient de 0,05 à 0,95 eV, et pour Ti–MoS2 de 0,05 à 0,46 eV. Les cadres théoriques conventionnels traitent la SBH comme une propriété d'interface locale, attribuant les écarts par rapport à la limite de Schottky–Mott aux états de gap d'interface et au piégeage du niveau de Fermi. Dans cette perspective, le désordre microscopique est censé être moyenné par les réservoirs de charges d'interface ou contraint par l'effet de « pincement » (pinch-off) des régions environnantes à haute barrière. Cependant, l'étalement géant de valeurs de SBH observé expérimentalement ne peut être réconcilié avec un modèle où la barrière est déterminée uniquement par des propriétés d'interface locales.

Méthodologie
Les auteurs proposent un modèle électrostatique nonlocal où la SBH effective dans les dispositifs 2D à contact supérieur est traitée comme une observable de bord pilotée par des perturbations électrostatiques distantes, plutôt que comme une distribution spatiale de barrières locales.

1. Cadre Théorique :

   • L'étude considère une géométrie de contact supérieur avec une région couverte ($x < 0$), un canal non couvert ($x > 0$), et un bord de contact à $x = 0$.
   • La SBH latérale ($\Phi^{\text{lat}}_B$) est définie comme le saut de la bande de conduction à travers le bord. Le modèle postule que le décalage de la barrière est contrôlé par le changement de l'alignement au bord ($\Delta J$), qui est linéairement couplé au potentiel électrostatique au bord, $u(0)$.
   • En utilisant une fonctionnelle d'énergie libre quadratique pour le potentiel électrostatique, les auteurs dérivent un noyau de réponse qui cartographie une source de défaut distant $s_d(x)$ vers le potentiel de bord.
   • Le décalage de barrière résultant dû à un défaut ponctuel à une distance $x_d$ est donné par :
     $$\Delta \Phi^{\text{lat}}_B (x_d) = \chi_e S_d \frac{\lambda_l}{\kappa + K_e \lambda_l} e^{-x_d/\lambda_l}$$
     où $\lambda_l$ est la longueur de propagation déterminée par la rigidité électrostatique en plan du canal ($\kappa$) et la capacité de rappel ($C_{ch}$), et l'amplitude est déterminée par des paramètres spécifiques au contact (susceptibilité de bord $\chi_e$ et rigidité de bord $K_e$).

2. Validation Computationnelle :

   • Le modèle est validé par des calculs de premier principe utilisant la Théorie de la Fonctionnelle de la Densité (DFT) combinée à la méthode des Fonctions de Green hors équilibre (NEGF).
   • Les systèmes étudiés incluent des contacts supérieurs Ti–MoS2 et Au–MoS2.
   • Des lacunes de soufre sont introduites dans le canal MoS2 à des distances variables ($0d$ à $10d$) du bord du contact.
   • Les SBH effectives sont extraites de la densité d'états locale (LDOS) près du bord du contact, plus précisément en mesurant le saut de la bande de conduction.

Contributions Clés et Résultats

• Mécanisme Nonlocal : L'article établit que des défauts situés loin du bord du contact peuvent moduler de manière significative la SBH mesurée via un couplage électrostatique nonlocal. L'influence d'un défaut décroît exponentiellement avec la distance, régie par une échelle de longueur caractéristique $\lambda_l$.
• Découplage des Rôles : Le modèle révèle une division claire des rôles :
  • Le canal semi-conducteur détermine la portée de l'effet nonlocal ($\lambda_l$), qui dépend de la propagation électrostatique en plan et de l'écrantage.
  • Le contact métallique détermine l'amplitude du décalage de la SBH, régi par la susceptibilité du bord aux changements de potentiel et l'efficacité de la conversion des décalages de potentiel de bord en modulation d'alignement.
• Accord Quantitatif :
  • Les calculs DFT-NEGF pour Ti–MoS2 et Au–MoS2 montrent que la SBH latérale suit la décroissance exponentielle prédite en fonction de la distance du défaut.
  • Une longueur d'interaction commune de $\lambda_l \approx 1,1$ nm s'ajuste aux données pour les deux métaux, confirmant que cette échelle de longueur est principalement une propriété de l'électrostatique du canal MoS2 plutôt que du métal spécifique.
  • Les variations de SBH calculées couvrent une plage qui rend compte d'une part substantielle de la variabilité rapportée expérimentalement pour Ti–MoS2 et d'une part significative de l'étalement pour Au–MoS2.
• Explication Unifiée : Les résultats fournissent une explication unifiée de la variabilité de SBH longtemps observée, démontrant que des contacts nominalement similaires peuvent produire des hauteurs de barrière très différentes selon la distribution spatiale des défauts par rapport au bord du contact.

Signification
L'article prétend résoudre une incohérence majeure dans la littérature des contacts 2D en démontrant que la barrière de Schottky mesurée n'est pas une propriété d'interface purement locale, mais une quantité de bord contrôlée de manière nonlocale. En identifiant l'environnement de bord et la distribution spatiale des défauts proches comme des variables de conception clés, ce travail suggère que l'injection de porteurs dans les dispositifs 2D est gouvernée par des environnements de désordre mésoscopiques plutôt que par de simples structures d'interface à l'échelle atomique. Ce cadre dépasse les modèles classiques de Schottky inhomogène et de pincement (pinch-off), offrant un mécanisme pour l'« étalement géant » des valeurs de SBH et établissant l'électrostatique nonlocale médiée par les défauts proches du bord comme un facteur fondamental de la performance des contacts 2D.