Pseudo-potentiel

Les bases d’ondes planes utilisent la même résolution dans chaque région de l’espace de sorte que pour décrire à la fois les cœurs ioniques (c’est-à-dire le noyau entouré du nuage électronique le plus interne) et les états électroniques partiellement localisés autour d’eux, le nombre de vecteurs ${\displaystyle G}$ nécessaires serait relativement prohibitif pour mener à bien la résolution des équations de Kohn-Sham.

Une façon de contourner cette difficulté consiste à utiliser la méthode des pseudopotentiels, appelée encore approximation de cœurs gelés (frozen-core approximation). Cette méthode repose sur l’hypothèse que seuls les électrons de valence (c’est-à-dire les électrons les plus externes) contribuent de façon significative aux propriétés physiques et chimiques d’un système donné alors que les électrons de cœur ne sont pas fortement sensibles à l’environnement chimique. Les cœurs ioniques sont ainsi considérés comme étant “gelés” dans leurs configurations atomiques. La méthode des pseudopotentiels consiste ainsi à ne traiter explicitement que les électrons de valence, qui se déplacent alors dans un potentiel externe effectif, produit par ces cœurs ioniques inertes, appelé pseudopotentiel. Ce pseudopotentiel tente de reproduire l’interaction générée par le vrai potentiel sur les électrons de valence sans inclure explicitement dans le calcul des électrons de cœur.

Le nombre d’électrons apparaissant de manière explicite dans les calculs (et donc le nombre d’équations à résoudre) s’en trouve fortement réduit ; seuls les électrons de valence sont pris en compte, ce qui produit un gain de temps de calcul considérable. Mais, comme en mécanique quantique toutes les fonctions d’onde décrivant les états électroniques (c’est-à-dire les électrons) doivent être orthogonales entre elles. Les électrons de cœur apparaissent toujours de manière implicite. En effet, pour rester orthogonales aux orbitales de cœur, les fonctions d’onde des électrons de valence oscillent rapidement dans la région proche du noyau (Figure 1). Les fonctions d’onde résultant de cette contrait d’orthogonalité restent relativement difficiles à décrire à partir d’une base d’ondes planes (nombre de vecteurs ${\displaystyle G}$ très élevé). On remplace alors la partie de la fonction d’onde proche du cœur ionique par des fonctions d’onde fictives, ou pseudo-fonctions d’onde, qui donnent lieu aux mêmes fonctions d’onde de valence au-delà d’un certain rayon de coupure ${\displaystyle r_{c}}$  (Figure 1). Ces pseudo-fonctions d’onde sont alors plus lisses ou “douces” que les vraies fonctions d’onde et elles peuvent de ce fait être représentées par un nombre de vecteurs  très inférieur à celui nécessité dans le traitement des vraies fonctions d’onde. On dit également que le potentiel très “dur” du cœur ionique est remplacé par un pseudopotentiel plus “doux” (Figure 1).

Le pseudopotentiel se définit donc comme un potentiel ionique effectif relativement doux qui n’agit que sur les électrons de valence.

Si le pseudopotentiel est ajusté de telle manière que la charge intégrée dans la région de cœur ionique correspondant à la pseudo-fonction d’onde soit égale à la charge intégrée associée à la vraie fonction d’onde, le pseudopotentiel est dit à norme conservée. Ainsi, bien que la méthode pseudo-potentielle simplifie fortement la description des électrons de valence, l’utilisation de pseudopotentiels à norme conservée permet de garantir une considération correcte de cette couche électronique externe. De tels pseudopotentiels sont construits de manière à satisfaire un certain nombre de conditions qui les rendent aussi transférables que possible, c’est-à-dire qu’ils puissent être utilisés pour prédire les propriétés chimiques d’un atome dans une large gamme de situations (ex. de l’état massif ou bulk, à la surface). Dans la mesure où les pseudopotentiels à norme conservée ne reposent pas sur la connaissance expérimentale préliminaire d’un élément chimique, on peut les construire pour n’importe quel élément du tableau périodique.

La génération d’un pseudopotentiel se fait à partir d’un atome isolé et on impose que les énergies propres obtenues avec le pseudopotentiel soient égales aux énergies atomiques réelles ou “all electron”. L’hypothèse principale est alors d’admettre que le pseudopotentiel construit pour un atome donné est transférable, c’est-à-dire que les résultats obtenus restent corrects si l’atome est placé dans un certain environnement chimique. La construction d’un tel potentiel est généralement assez délicate, car il faut estimer un rayon de coupure convenable et choisir les électrons pouvant être considérés comme chimiquement inertes (cœur et semi-cœur). On essaie donc de trouver un compromis entre la transférabilité du pseudopotentiel et la diminution de temps de calcul qu’il engendre en autorisant un ${\displaystyle E_{cut}}$ moins important. En général, il est obligatoire de tester un nouveau pseudopotentiel dans des environnements bien connus pour vérifier s’il reproduit bien les résultats auxquels on s’attend.

Plusieurs types de pseudopotentiels ont déjà été développés[1]. Certaines méthodes font appel à des pseudopotentiels ne conservant pas la norme. Cela confère plus de latitude dans leur construction, mais la non conservation de la norme néanmoins des contraintes techniques supplémentaires lors de l’implémentation. Ces pseudopotentiels sont caractérisés par des pseudo-fonctions d’onde arbitrairement lisses dans les régions du cœur (ex. pseudopotentiels de Vanderbilt[2] également appelés pseudopotentiels “ultrasoft”, USPP) permettant de bas .

Le concept de pseudopotentiel a été introduit dans les années 1930 par Fermi[3]. Par la suite, Hellmann utilise cette notion pour le calcul des niveaux énergétiques de métaux alcalins[4] - [5]. Ces premiers pseudopotentiels sont qualifiés d'empiriques; ce qui signifie qu'ils ne sont pas obtenus par calcul mais paramétrés pour reproduire au mieux des résultats expérimentaux de référence. L'utilisation d'un tel type de pseudopotentiel se fonde sur deux constatations. Tout d'abord, s'il était certain à cette époque que l'on pouvait obtenir de manière exacte et par calcul des pseudopotentiels cela passait par la résolution de calculs complexes (impliquant les fonctions d'onde atomiques) impossible à résoudre sans ressources informatiques. L'utilisation d'une méthode empirique beaucoup plus simple allait donc de soi. Ensuite, on peut constater qu'un certain nombre d'éléments peuvent être décrits par des pseudopotentiels paramétrés tout en fournissant une approximation acceptable de l'interaction électron de valence - ion de cœur. Cela était vrai notamment pour les métaux alcalins, les métaux "simples" tels que l'aluminium ainsi que les semi-conducteurs. L'utilisation de ces pseudopotentiels va permettre en une dizaine d'années d'augmenter le champ des connaissances dans le domaine de l'état solide avant d'être remplacés par des pseudopotentiels ab initio plus efficaces.

La méthode des ondes planes orthogonales (OPW pour orthogonalized planes waves) a été introduite par Convers Herring dans les années 1940[6]. La méthode permit de mieux comprendre la nature de structure de bande de matériaux semiconducteurs tels que le silicium et le germanium et fut la première à expliquer de manière théorique que le silicium est un matériau à gap indirect[7]. Le développement de cette méthode mérite d’être introduit car celle-ci est l’ancêtre direct de la notion de pseudopotentiel.

Formalisme mathématique

De manière concrète, la méthode OPW est une approche générale qui vise à construire des fonctions de bases pour la description des états de valence. Ces fonctions sont définies de la manière suivante :

${\displaystyle \chi _{q}^{OPW}(r)={\frac {1}{\Omega }}[e^{iqr}-\sum _{j}<u_{j}|q>u_{j}(r)]}$

Les fonctions uj sont arbitraires mais nécessitent d’être localisées autour des noyaux. De la définition précédente, il s’ensuite que ${\displaystyle \chi _{q}^{OPW}}$ est bien orthogonale à toutes les fonctions uj c’est-à-dire que pour tout uj :

${\displaystyle <u_{j}|\chi _{q}^{OPW}>=0}$

Si les fonctions uj sont correctement choisies, l’expression (X) peut alors être vue comme étant la somme de deux contributions ; une partie adoucie (soft), c’est-à-dire ne comportant pas de nœuds, et une partie localisée. La partie adoucie peut être représenté aisément par une combinaison d’ondes planes ce qui était l’objectif d’Herring comme il le précise lui-même[6]

« This suggest that it would be practical to try to approximate [the eigenfunction] by a linear combinaison of a few planes waves, plus a linear combinaison of a few function ${\displaystyle u_{j}}$ localized about each nucleus and obeying wave equations of the form :

${\displaystyle \nabla ^{2}u_{j}+(E_{j}-V_{j})u_{j}=0}$ »

Le potentiel Vj et les fonctions uj qui apparaissent dans l’équation doivent être choisie de manière optimales. Il est important de préciser que, les états de valence étant étiquetés, par leurs nombres quantiques l et m, il est normal que les fonctions que l’on ajoute possèdent également les mêmes valeurs de l et m. Des précédentes définitions, les fonctions d’onde de valence peuvent être exprimées comme :

C'est par une reformulation de la méthode OPW d'Herring qu'en 1959 Phillips et Kleinman développent la première approche formelle (non empirique) de la notion de pseudopotentiel[8]. Antoncik, de manière indépendante, publie la même année une approche similaire[9]. La méthode du pseudopotentiel de Phillips-Kleinman-Antoncik (PKA) est la première à montrer que la condition d'orthogonalité dans la région de cœur entre les états de cœur et de valence agit comme un potentiel répulsif qui tend à s'opposer au potentiel nucléaire attractif ressenti par les électrons de valence. Le plus souvent, ces deux effets se combinent pour former un potentiel faiblement répulsif, le pseudopotentiel[10].

L'introduction du principe de pseudopotentiel à norme conservée représente la percée la plus significative dans le traitement des électrons de cœur. L'approche fut développée par Hamann, Schlüter et Chiang[11] et fut suivie peu de temps après par une méthode très similaire mais établie de manière indépendante par Kerker[12].

Méthodologie

La première étape dans la génération d'un pseudopotentiel à norme conservée consiste à réaliser un calcul atomique ab initio tous-électrons. La configuration électronique de l'atome peut être choisie arbitrairement. Généralement il s'agit de l'atome à l'état neutre. Le calcul se fait au départ de l'équation de Kohn-Sham écrite sous sa forme radiale :

${\displaystyle \left[-{\frac {1}{2}}{\frac {d^{2}}{dr^{2}}}+{\frac {l(l+1)}{2r^{2}}}+{\frac {-Z}{r}}+V_{H}+V_{XC}\right]rR_{nl}(r)=\epsilon _{nl}rR_{nl}(r)}$

La fonction d'onde réelle est ensuite remplacée par une pseudofonction d'onde à laquelle est associée un potentiel modèle (le pseudopotentiel) choisi pour reproduire correctement les propriétés des états de valence. Le schéma mathématique qui permet de générer la pseudofonction d'onde n'est pas unique. Il existe du point de vue mathématique une certaine liberté dans le choix de la méthode comme en témoignent les nombreuses publications qui présentent différentes manière de faire[11] - [12] - [13] - [14] - [15] - [16] - [17] - [18].