Nonequilibrium Green's Function Formalism Applicable to Discrete Impurities in Semiconductor Nanostructures

Cet article présente un nouveau cadre de fonctions de Green hors équilibre qui prend explicitement en compte la nature discrète des impuretés dans les nanostructures semi-conductrices en séparant la diffusion à courte portée et les potentiels de Hartree à longue portée, révélant ainsi la non-localité intrinsèque de la diffusion par les impuretés et démontrant son impact significatif sur les propriétés de transport dans les systèmes quasi unidimensionnels.

L'essentiel

Imaginez une nanostructure semi-conductrice (comme un minuscule fil utilisé dans les puces informatiques futures) comme un long couloir étroit. À l'intérieur de ce couloir, les électrons tentent de courir d'une extrémité à l'autre pour transporter un courant électrique. Cependant, le couloir n'est pas vide ; il est rempli d'« impuretés » — de la saleté ou des débris laissés accidentellement lors de la fabrication. Ces impuretés sont en réalité des atomes individuels (des dopants) qui agissent comme des obstacles.

Pendant des décennies, les scientifiques ont modélisé ces obstacles en faisant semblant qu'ils formaient un brouillard lisse et invisible réparti uniformément dans tout le couloir. Ils supposaient que, parce qu'il y avait tant d'obstacles, les électrons ne verraient qu'un « nuage » moyen de résistance. Cela fonctionnait bien pour les couloirs larges et grands.

Mais dans les fils minuscules et ultra-fins de la technologie moderne, cette idée de « brouillard » s'effondre. Le couloir est si étroit que l'emplacement spécifique de chaque grain de poussière compte. Si un grain est juste au milieu du chemin, il bloque l'électron. S'il est sur le côté, l'électron peut peut-être le contourner. L'ancien modèle de « brouillard » manque ce détail crucial.

Le Nouveau Cadre : Deux Types de Tracas

Ce papier, par Nobuyuki Sano, propose une nouvelle façon de calculer comment les électrons se déplacent dans ces fils minuscules en traitant les impuretés comme des points distincts et individuels plutôt que comme un brouillard. L'auteur divise le problème d'une impureté en deux parties, en utilisant une analogie ingénieuse :

1. La Partie « Longue Portée » (L'Effet de Quartier) : Imaginez qu'une impureté est une personne debout dans le couloir. Même si vous ne les touchez pas, leur présence modifie légèrement l'atmosphère. Ils peuvent repousser les gens ou les attirer de plus loin. En physique, c'est le champ électrique « longue portée ». Le papier traite cela comme un potentiel d'arrière-plan lisse et auto-cohérent (comme une pente douce dans le couloir) qui affecte tout le monde.
2. La Partie « Courte Portée » (Le Risque de Trébuchement) : C'est le heurt immédiat et net contre lequel vous trébuchez si vous marchez exactement sur l'impureté. C'est la diffusion « courte portée ». Le papier traite cela comme une collision spécifique et localisée qui ne se produit que lorsqu'un électron s'approche très près d'un atome d'impureté spécifique.

Le Système de Coordonnées « Fantôme »

La découverte la plus surprenante du papier concerne où ces collisions se produisent.

En physique traditionnelle, nous pensons qu'une collision se produit à un endroit spécifique sur une carte (l'Espace Réel). Si une impureté est à la position X, la collision se produit à X.

Cependant, ce papier montre que dans le monde quantique de ces fils minuscules, l'« emplacement » d'une collision est en fait un mélange de l'endroit où l'électron était et de l'endroit où il va. L'auteur utilise un outil mathématique appelé coordonnées de Wigner (spécifiquement le « centre de masse » du trajet de l'électron) pour décrire cela.

L'Analogie :
Pensez à un flou de mouvement. Si vous prenez une photo d'une voiture en mouvement rapide, vous ne la voyez pas à un point exact ; vous voyez une traînée. Le papier soutient que le « taux de diffusion » (la probabilité qu'un électron rebondisse sur une impureté) n'est pas lié à un seul point sur la carte. Au lieu de cela, il est lié à la position moyenne du voyage de l'électron.

Cela signifie que la diffusion est non locale. L'électron « ressent » l'impureté non seulement lorsqu'il la touche, mais en fonction d'une relation plus large et floue entre ses positions passées et futures. C'est comme si l'électron avait un sens « fantomatique » de l'obstacle qui s'étend au-delà du point de contact physique.

Ce Qui Se Passe Lorsque Vous Utilisez le Nouveau Modèle ?

L'auteur a appliqué ces nouvelles mathématiques pour simuler un fil cylindrique (un nano-fil) et l'a comparé aux anciens modèles de « brouillard » :

• L'Ancien Modèle (Local/Diagonal) : Il suppose que la diffusion se produit à un point unique et agit comme un mur simple. Ce modèle tend à surestimer la vitesse à laquelle les électrons peuvent se déplacer (la mobilité). Il pense que les électrons sont moins « confus » par les obstacles qu'ils ne le sont réellement.
• Le Nouveau Modèle (Non Local/Hors Diagonale) : Parce qu'il prend en compte la nature « floue » de la collision, il montre que les électrons perdent leur « cohérence de phase » (leur rythme synchronisé) beaucoup plus rapidement. Ils se confondent et sont diffusés plus facilement.
• Le Résultat : Le nouveau modèle prédit que le courant électrique et la mobilité sont en réalité plus faibles que ce que les anciens modèles suggéraient, en particulier lorsque le nombre d'impuretés est modéré (ni trop peu, ni trop).

La Surprise de l'« Auto-Moyennage »

Le papier a également trouvé quelque chose d'intéressant concernant les moyennes. Si vous prenez de nombreuses dispositions aléatoires différentes d'impuretés et que vous les moyennez (comme regarder une foule de loin), le nouveau modèle « non local » correspond toujours étonnamment bien aux résultats de l'ancien modèle de « brouillard ».

Cependant, si vous regardez un seul fil spécifique avec une disposition spécifique d'impuretés, l'ancien modèle échoue complètement. Il manque les variations sauvages de performance qui se produisent d'un fil minuscule à l'autre, simplement parce que les grains de poussière sont tombés à des endroits légèrement différents.

En Résumé

Ce papier fournit une « carte » plus précise pour naviguer dans le monde quantique des fils minuscules. Il nous dit que nous ne pouvons pas simplement traiter les impuretés comme un brouillard lisse ou comme de simples bosses ponctuelles. Nous devons reconnaître que dans le domaine quantique, l'« emplacement » d'une collision est un peu flou et dépend de tout le trajet de l'électron. En faisant cela, nous obtenons une image plus vraie de la vitesse à laquelle l'électricité peut réellement circuler dans la prochaine génération de puces informatiques, révélant qu'elles pourraient être légèrement plus lentes (et plus variables) que prévu précédemment.

Résumé technique

1. Énoncé du problème

Dans les nanostructures semi-conductrices dopées, les propriétés de transport électronique présentent une variabilité significative d'un dispositif à l'autre, même à température ambiante. Cette variabilité découle de la nature discrète des impuretés (dopants), que les approches de modélisation traditionnelles ne parviennent pas à capturer avec précision.

• Limites du formalisme NEGF traditionnel : Les schémas standards de fonctions de Green hors équilibre (NEGF) supposent généralement un « auto-moyennage », où les impuretés sont traitées comme un fond continu de type « jellium ». Cela suppose que le dispositif est beaucoup plus grand que la longueur de cohérence de phase, moyennant ainsi les détails microscopiques.
• Limites des approches alternatives : Traiter les impuretés comme des charges ponctuelles explicites dans le hamiltonien (non perturbé) conduit à des potentiels singuliers, provoquant des instabilités numériques dans les méthodes aux différences finies.
• Le vide théorique : Il existe un besoin d'un cadre théorique rigoureux qui intègre la distribution spatiale discrète des impuretés dans l'auto-énergie NEGF sans singularités numériques, tout en distinguant correctement les fluctuations de potentiel à longue portée et la diffusion à courte portée.

2. Méthodologie

L'auteur propose un nouveau cadre théorique qui couple le schéma NEGF à l'équation de Poisson de manière auto-cohérente, spécifiquement conçu pour gérer les impuretés discrètes.

A. Séparation du potentiel d'impureté

Le potentiel électrostatique d'une impureté discrète est divisé en deux composantes :

1. Composante à longue portée ($V_{long}$) : Traitée comme un potentiel de Hartree auto-cohérent. Elle est résolue via l'équation de Poisson couplée à la densité de porteurs. Cela rend compte des fluctuations de potentiel causées par l'écran collectif des impuretés.
2. Composante à courte portée ($V_{short}$) : Traitée comme un potentiel de diffusion inclus dans l'auto-énergie ($\Sigma_{imp}$). Cela représente les événements de diffusion localisés.

B. Nouvelle formulation de l'auto-énergie

L'innovation centrale réside dans la dérivation de l'auto-énergie de diffusion par impureté ($\Sigma_{imp}$) :

• Coordonnées de Wigner : L'auto-énergie est dérivée en utilisant les coordonnées de Wigner ($X = (r+r')/2$, $\xi = r-r'$) plutôt que les coordonnées d'espace réel standard.
• Non-localité : La matrice d'auto-énergie résultante est intrinsèquement non diagonale dans l'espace réel. La dépendance spatiale du taux de diffusion est régie par la coordonnée du « centre de masse » ($X$) dans la représentation de Wigner, et non par la coordonnée locale d'espace réel ($r$).
• Forme mathématique : L'auto-énergie est exprimée comme une somme sur les positions discrètes des impuretés ($R_l$) impliquant une fonction de corrélation $C_X(\xi)$, qui est la transformée de Fourier du potentiel écranté au carré.
  $$\Sigma^\alpha_{imp}(r, r'; E) \approx \sum_l \delta(X - R_l) C_X(\xi) G^\alpha(X, \xi; E)$$

C. Réduction au taux de diffusion

L'article démontre que l'application d'un « auto-moyennage » (moyenne d'ensemble) sur les distributions d'impuretés à ce nouveau formalisme retrouve le taux de diffusion standard dérivé de la règle d'or de Fermi. Cependant, il clarifie que la dépendance spatiale dans le taux de diffusion est fondamentalement liée à la coordonnée de Wigner ($X$), et non à la coordonnée locale ($r$).

D. Modèle d'application

Le cadre est appliqué à des fils minces cylindriques (rayon $\rho_s = 4$ nm, longueur $L_c = 30$ nm) dans le cadre d'une approximation quasi-1D :

• L'équation de Poisson 3D fournit un potentiel à longue portée « figé » le long de l'axe du fil.
• Le calcul NEGF est effectué en 1D en utilisant l'approximation de la sous-bande la plus basse.
• La matrice de taux de diffusion est construite explicitement sur la base des positions discrètes des impuretés.

3. Contributions clés

1. Dérivation rigoureuse de la diffusion dépendante de la position : L'article dérive systématiquement la dépendance spatiale des taux de diffusion par impureté sous des profils d'impuretés inhomogènes, prouvant que le taux de diffusion dépend de la coordonnée du centre de masse (Wigner) plutôt que de la coordonnée locale d'espace réel.
2. Non-localité intrinsèque : Il établit que la diffusion par impuretés est intrinsèquement non locale dans l'espace. La matrice de taux de diffusion est non diagonale, reflétant l'interférence de phase quantique, contrairement aux matrices diagonales utilisées dans les approximations « locales » traditionnelles.
3. Cadre unifié : Il comble avec succès le fossé entre la modélisation des impuretés discrètes et les modèles traditionnels moyennés sur l'ensemble. Le nouveau formalisme se réduit au modèle « jellium » standard lors de l'auto-moyennage, mais conserve la capacité de simuler des configurations discrètes spécifiques.
4. Résolution des singularités : En traitant la partie à courte portée comme un potentiel de diffusion dans l'auto-énergie (plutôt que comme un potentiel singulier dans le hamiltonien), la méthode évite les instabilités numériques associées aux simulations de charges ponctuelles.

4. Résultats et discussion

Le cadre a été appliqué à des nanofils de silicium avec des densités de donneurs variables ($\bar{N}_d^+$ de $10^{18}$ à $10^{20}$ cm$^{-3}$).

• Fluctuations de potentiel : Les impuretés discrètes créent des fluctuations de potentiel de l'ordre de dizaines de mV (jusqu'à ~100 mV à fort dopage). Ces fluctuations sont suffisamment importantes pour provoquer une localisation électronique locale même à température ambiante.
• Matrices de taux de diffusion :
  • Cas discret : Les matrices sont non diagonales, symétriques et localisées autour des positions des impuretés. Les éléments hors diagonale peuvent être négatifs, indiquant une forte interférence de phase (similaire aux fonctions de Wigner).
  • Cas auto-moyenné (non local) : Même après moyennage, la matrice reste non diagonale en raison de la dépendance à la coordonnée de Wigner.
  • Approximation locale (traditionnelle) : Suppose une matrice diagonale, ce qui échoue à capturer les effets de déphasage.
• Propriétés de transport (Courant et mobilité) :
  • Surestimation dans les modèles traditionnels : Les modèles traditionnels « locaux » auto-moyennés surestiment considérablement la mobilité et le courant électroniques car ils sous-estiment le déphasage causé par la diffusion élastique.
  • Discret vs Ensemble : Dans le régime de dopage modéré ($10^{18} \leq \bar{N}_d^+ \leq 10^{19}$ cm$^{-3}$), les propriétés de transport d'un dispositif spécifique varient considérablement selon la configuration exacte des impuretés. Le modèle « auto-moyenné non local » fournit une bonne approximation de la moyenne d'ensemble, tandis que le modèle « local » ne le fait pas.
  • Données quantitatives : Pour $\bar{N}_d^+ = 10^{19}$ cm$^{-3}$, le modèle local a prédit une mobilité de 215 cm$^2$/V/s, tandis que le modèle auto-moyenné non local plus précis a prédit 100 cm$^2$/V/s, et la simulation discrète a donné 95 cm$^2$/V/s.

5. Importance

• Variabilité des dispositifs : Ce travail fournit une base théorique pour comprendre et prédire la variabilité statistique dans les dispositifs semi-conducteurs de nouvelle génération (par exemple, FinFET, FET à nanofils) où le nombre de dopants est faible et leurs positions aléatoires.
• Précision de la modélisation : Il met en évidence que l'approximation « locale » pour les taux de diffusion est insuffisante pour les nanostructures. Ignorer la nature non locale de la diffusion (dépendance de la coordonnée de Wigner) conduit à des prédictions erronées de la mobilité et du courant.
• Perspectives futures : Le cadre est généralisable à d'autres mécanismes de diffusion (par exemple, diffusion par phonons, rugosité de surface), suggérant que la non-localité est une caractéristique fondamentale du transport quantique dans les systèmes désordonnés qui doit être prise en compte dans la simulation précise des dispositifs.

En conclusion, Sano présente un formalisme NEGF robuste et auto-cohérent qui traite correctement la nature discrète des impuretés, révélant que la diffusion par impuretés est intrinsèquement non locale et que les approximations locales traditionnelles surestiment considérablement les performances de transport dans les dispositifs à l'échelle nanométrique.