Méthodes Ab-Initio
Méthodes Ab-Initio
Méthodes Ab-Initio
Les principales méthodes ab-initio utilisés pour les calculs des structures électroniques
peuvent être classés en trois catégories : la méthode Hartree-Fock, les méthodes post Hatree-
Fock et les méthodes basées sur la théorie de la fonctionnelle de la densité. Dans les méthodes
Hatree-Fock et post Hartree-Fock les principales quantités considérées pour les calculs sont
les fonctions d’ondes. Dans la théorie de la fonctionnelle de la densité, on considère plutôt la
densité électronique.
Les propriétés physiques d’un système solide, illustré par l’image de négatons légers
en mouvement autour de noyaux lourds, dépendent de sa structure électronique et
géométrique. Ces propriétés physiques se déduisent de l’équation de Schrödinger
indépendante du temps.
HΨ=EΨ
On a explicitement
Avec :
𝑻𝒆 l’Opérateur énergie cinétique des négatons
𝑻𝒏𝒐𝒚 l’Opérateur énergie cinétique des noyaux
𝑽𝒏𝒐𝒚-𝒏𝒆𝒈 Interaction noyau-négaton
𝑽𝒏𝒆𝒈-𝒏𝒆𝒈 Interaction négaton-négaton
𝑽𝒏𝒐𝒚-𝒏𝒐𝒚 Interaction noyau-noyau
Les deux premiers termes 𝑻𝒆 et 𝑻𝒏𝒐𝒚 désignent les énergies cinétiques du système et
les termes 𝑽𝒏𝒐𝒚-𝒏𝒆𝒈, 𝑽𝒏𝒆𝒈-𝒏𝒆𝒈 et 𝑽𝒏𝒐𝒚-𝒏𝒐𝒚 désignent les énergies potentielles.
Approximation Born-Oppenheimer
[𝑇𝑒 + 𝑉𝑛𝑒𝑔 − 𝑛𝑒𝑔({𝑟 𝑖}) + 𝑉𝑖𝑜𝑛 − 𝑛𝑒𝑔(𝑅⃗ I , r i )]𝛷({r i}) = 𝐸(𝑅⃗ I)𝛷({r i}, {𝑅⃗ I })
[𝑇𝑖𝑜𝑛𝑠 + 𝑉𝑖𝑜𝑛 − 𝑖𝑜𝑛({𝑅⃗ I }) + 𝐸({𝑅⃗ I})]𝛤({𝑅𝑅⃗ I }) = 𝐸𝑡𝑜𝑡 ({𝑅⃗ I })𝛤({𝑅⃗ I})
𝐸 = 𝐸(𝑅⃗ 𝐼 ) est l’énergie des négatons pour une configuration donnée des noyaux.
C’est une valeur propre de l’opérateur hamiltonien négatonique 𝑯 qui est donné par
H = 𝑻𝒆 +𝑽𝒏𝒐𝒚-𝒏𝒆𝒈 +𝑽𝒏𝒆𝒈-𝒏𝒆𝒈
⃗⃗⃗ 𝐼 }) = Φ(r1 , r2 , … , r𝑛 , {R
Φ({r⃗⃗⃗ 𝑖 }, {R ⃗⃗⃗ 𝐼 })
⃗⃗⃗ 𝐼 })Φ2 ( r2 , {R
= Φ1 ( r1 , {R ⃗⃗⃗ 𝐼 }) … Φ𝑛 ( r𝑛 , {R
⃗⃗⃗ 𝐼 })
Et chaque fonction d’ondes 𝛷𝑖 (𝑟𝑖 ) d’un négaton obéit à une équation de Schrödinger
de la forme
ℏ2
[− (𝛻⃗⃗⃗𝑟𝑖 )2 + 𝑈(𝑟𝑖 , {R
⃗⃗⃗ 𝐼 })]Φ𝑖 = 𝜖𝑖 Φ𝑖
2𝑚𝑒
⃗⃗⃗ 𝐼 }) est écrit sous forme
𝑈(𝑟𝑖 , {R
𝑁
𝑒2 𝑍𝐼
⃗⃗⃗ 𝐼 }) =
𝑈(𝑟𝑖 , {R ∑ ⃗⃗⃗ 𝐼 })
+ 𝑉(𝑟𝑖 , {R
4𝜋𝜀0 ⃗
||𝑅𝐼 − 𝑟𝑖 || 𝐼=1
avec 𝑉(𝑟 , {𝑅 𝐼 }) l’énergie potentielle associée à l’action des autres négatons sur un
électron 𝑖 considéré.
𝐹𝑖 𝛷 𝑖 (𝑟 𝑖, 𝜗𝑖 , {𝑅⃗𝐼 }) = 𝜖𝑖 𝛷 𝑖 (𝑟 𝑖, 𝜗𝑖 , {𝑅⃗𝐼 })
Les équations de Hartree ou de Hartree-Fock (ou Hartree) peuvent être déduite d’un
principe variationnel, en considérant la valeur moyenne de l’énergie comme une fonctionnelle
de l’état (ou de manière équivalente une fonctionnelle de la fonction d’onde) : 〈𝐸〉 = 〈𝐸〉(|Φ⟩
) = ⟨Φ|𝑯|Φ⟩ et en la minimisant.
Mais ses méthodes sont toujours aussi inefficaces quand on a rencontré plus tard des
systèmes plus complexe.
Cette approche présente plusieurs avantages qui facilitent son utilisation dans
différents domaines, surtout en chimie et science des matériaux.
• D’abord, elle inclut dans son formalisme une grande part de la corrélation
électronique.
• Elle peut être appliqué a des systèmes aussi différents que les matériaux covalents,
ioniques.
• L’aspect mono électronique permet de donner une interprétation « chimique » de la
fonction d’onde.
Dans la théorie de la DFT, l’énergie c’est une fonction de la densité électronique ρ du
système.
Densité électronique
Dans les sections précédentes, nous avions défini les négatons comme étant des
particules indissociables. Or, en tant que particule isolée, le négaton ne peut pas être localisé :
nous parlerons plutôt de sa probabilité de présence dans un élément de volume dr, définissant
ainsi la densité électronique. En fait, les négatons sont considérés dans leur ensemble (nuage
électronique) et la densité électronique localisera les régions de l’espace ou les négatons sont
les plus probables. Cette probabilité de présence sera définie comme étant l’intégrale multiple
du carré de la fonction d’onde. Il s’agit d’une fonction positive uniquement des trois
coordonnées spatiales (x, y, z). Cette quantité s’annule à l’infini, intégrée sur tout l’espace,
elle est égale au nombre total de négatons N :
𝜌(𝑟 → ∞) = 0,
∫ 𝜌(𝑟)𝑑𝑟 = 𝑁
Contrairement à la fonction d’onde, la densité électronique est une observable qui peut
être mesurée diffraction X et neutronique. De ce fait, la densité électronique est riche
d’informations : « la connaissance de la densité électronique est tout ce dont nous avons
besoin pour une détermination complète des propriétés moléculaires » (E. Bright Wilson,
1965).
Méthode de Thomas-Fermi
Le premier terme est l’énergie cinétique des négatons. Le second terme est l’énergie
potentielle qui résulte de l’interaction du système avec un potentiel extérieur 𝑉𝑒𝑥𝑡(𝑟 ) qui peut
être le potentiel créé par des noyaux atomiques ou des ions. Le dernier terme correspond aux
interactions coulombiennes entre les négatons eux-mêmes. L’expression de l’énergie
cinétique est déduite de la statistique de Fermi-Dirac d’un système de négatons libre : c’est
donc une approximation.
La méthode de Thomas-Fermi donne des résultats utiles dans certains cas mais elle est
basée sur plusieurs approximations qui limitent ses domaines de validités pratiques.
E=E ( 𝝆𝟎 )
E0( 𝝆𝟎 )=min E( 𝝆𝟎 )
Les conséquences des deux théorèmes fondamentaux de la DFT sont que la variation
du potentiel externe implique une variation de la densité. L’équation d’énergie des fonctions
mono électroniques est alors décrite par l’expression :
𝐄[𝛒] = 𝐅𝐇𝐊 (𝛒) + ∫(𝐫)𝐕𝐞𝐱𝐭 𝐝𝟑 𝐫
La fonctionnelle 𝐅𝐇𝐊 (𝛒) est universelle pour n’importe quel système à plusieurs
négatons. Si la fonctionnelle 𝐅𝐇𝐊 (𝛒) est connue, alors il sera relativement facile d’uriliser le
principe variationnel pour déterminer l’énergie totale et la densité électronique de l’état
fondamental pour un potentiel extérieur donné. Malheureusement, le théorème de Hohenberg
et Kohn ne donne aucune indice sur la forme de 𝐅𝐇𝐊 (𝛒) .
Théorème de Kohn-Sham
Cette approche a pour but de déterminer les propriétés exactes d’un système à
plusieurs particules en utilisant des méthodes à particules indépendantes. En pratique, cette
révolution en la matière a permis d’effectuer certaines approximations qui se sont révélées
très satisfaisantes. L’approche de Kohn et Sham remplace le système à particules interagissant
entre elles et qui obéit à l’Hamiltonien, par un système moins complexe facilement résolu.
Cette approche assure que la densité à l’état fondamental du système en question est égale à
celle de certains systèmes à particules non interagissant entre elles.
L’énergie cinétique d’un gaz de négatons en interaction est inconnue. L’idée de Kohn
et de Sham est de remplacer le terme 𝑇’[𝜌(𝑟)] par le terme 𝑇 [ 𝜌(𝑟)] représentant l’énergie
cinétique du gaz de négaton sans interaction. Cette correspondance entre un système de
négatons en interaction et un système sans interaction a en fait plusieurs conséquences, parmi
lesquelles le passage d’une description, basée sur la fonction d’onde à N négatons, à une autre
avec N fonctions d’ondes mono-électroniques.
Étant donné que la forme de 𝑇[𝜌] et de 𝐸𝑛𝑐𝑙 [𝜌] est inconnue, il est impossible de
résoudre analytiquement l’équation de Schrödinger multi électronique liée au système
d’électrons en interaction. En 1965, Kohn et Sham proposent de réécrire la fonctionnelle E [ρ]
de la manière suivante :
𝑇0[ρ] est l’énergie cinétique d’un système d’électrons indépendants, et J[ρ] est la
contribution classique à la fonctionnelle de l’interaction électron-électron. Évidemment 𝑇0
n’est pas équivalente à l’énergie cinétique du système en interaction ; la différence des deux et
la contribution non classique à l’interaction électron-électron composent un terme appelé
énergie d’échange- corrélation 𝐸XC :
𝑀
𝑍 1 𝑁
𝐸 [𝜌] = 𝛴 1𝑁 ∫ 𝛴𝐴=1 2
|Ѱ𝑖(𝑟1 )| 𝑑𝑟1 − 𝛴 𝑖=1 < Ѱ𝑖|𝛥2|Ѱ𝑖
𝑟1 2
𝑁 𝑁 | 2
1
> + 𝛴𝑖=1 𝛴𝑖<𝑗 Ѱ𝑖 (𝑟1 )| + 𝑑𝑟 𝑑𝑟 + 𝐸𝑋𝐶 [𝜌]
𝑟12 1 2
En appliquant le principe variationnel par rapport aux fonctions d’onde ψi, on arrive à
un ensemble d’équations de Schrödinger mono électroniques appelées équations de Kohn et
Sham :
𝛻
(− + 𝑣𝑠 (𝑟))𝜓𝑖 = 𝜀𝑖 𝜓𝑖
2
𝑣𝑠 Représente le potentiel local fictif de Kohn et Sham, dans lequel se déplacent les
fermions indépendants :
𝑧𝐴
𝑣𝑠 (𝑟1) = ∫ 𝜌(𝑟2 )𝑟12 𝑑𝑟12 𝑑𝑟2 + 𝑣𝑋𝐶 (𝑟1) −
𝑟1𝐴
Les pseudo-potentiels
Les premiers pseudo-potentiels ont été mis en place en 1958 par Phillips et Kleiman.
Les Pseudo-potentiels utilisés actuellement sont déterminés à partir de calculs « tous négatons
», ce qui rend la méthode plus précise. Les pseudo-potentiels associés à des énergies de
coupure (Ecut-off) élevées sont appelés « hard », par comparaison aux pseudos potentiels
appelés « soft ». Vanderbilt a développé une base d’ondes plus réduite que ces pseudo-
potentiels traditionnels, avec des énergies de coupure encore plus basses. Ces pseudo-
potentiels sont appelés « ultra soft ». Ceci a permis de considérer des systèmes plus
complexes et/ou de grande taille.
Les fonctions d’ondes des négatons dans un cristal peuvent être décomposées en une
série d’ondes planes. Comme celle –ci possède a priori un nombre infini d’éléments, on est
obligé de tronquer ce développement en série, et, pour garder des temps de calcul
raisonnables, on doit souvent travailler avec une base minimale. Malheureusement, on besoin
d’un très grand nombre d’ondes planes pour décrire les négatons de cœur fortement liés aux
noyaux ainsi que les fonctions d’ondes des négatons de valence qui oscillent rapidement dans
la région centrale d’un atome.
Le fait d’utiliser des pseudo-potentiels permet d’utiliser une base plus restreinte ce qui
réduit les temps de calculs. Cette approximation est basée sur les idées suivantes :
-les propriétés chimiques des atomes sont principalement déterminées par les négatons
de valence tandis que les négatons de cœur restent chimiquement inertes. On n’a donc
pas besoin de recalculer leurs fonctions d’ondes lorsque l’atome est placé dans un
cristal (approximation des cœurs gelés).
-les fonctions d’ondes des négatons de valence oscillent rapidement dans la région
centrale de l’atome. Ces oscillations assurent l’orthogonalité aux orbitales de cœur
mais n’influencent que très peu les propriétés chimiques de l’atome. Le potentiel dû au
noyau et aux négatons profonds peut donc être remplacé par un pseudo-potentiel
variant plus lentement que le potentiel réel ce qui permet d’éviter ces nœuds dans la
partie radicale de la fonction d’onde. Un pseudo-potentiel tient compte de toutes les
interactions entre les négatons de valence et le cœur ionique (l’interaction
coulombienne avec le noyau écranté par les négatons de cœur, répulsion due au
principe de Pauli, interaction d’échange-corrélation entre négatons de cœur et négatons
de valence). Il facilite la résolution numérique de l’équation de Schrödinger et abaisse
le nombre d’ondes planes nécessaires pour représenter correctement les fonctions du
Bloch. En plus, l’approximation des cœurs gelés diminue le nombre de négatons pour
lesquels l’équation de Kohn-Sham doit être résolue et permet de se débarrasser
souvent des effets relativistes.
Pour les négatons de valence d’un solide, on peut négliger les effets relativistes. Ceci
n’est plus valable pour les négatons des couches profondes pour lesquels ces effets jouent un
rôle très important. On doit par conséquent en tenir compte lors de la construction du pseudo
potentiel. Celui-ci reprend donc tous les éléments du tableau périodique, des pseudo potentiels
ont été générés par différents groupes de recherche comme par exemple le groupe de
C.Hartwigsen, S.Goedecker et J.Hutter (HGH) . Mais en utilisant ces pseudos potentiels, on
est souvent très limité en ce qui concerne le nombre de négatons de valence ou
l’approximation de l’énergie d’échange-corrélation (la plupart de ces pseudo potentiels ne
permettent que l’utilisation de la LDA). Pour cela, nous avons préféré générer les pseudo-
potentiels utilisés dans nos calculs à l’aide du programme ATOM qui accompagne le code
SIESTA qui sera exposé au chapitre IV. En plus, nous avons pu utiliser soit la LDA soit la
GGA comme approximation de l’énergie d’échange-corrélation.
Pour des fermions chargés tels que les négatons, la probabilité de trouver le négaton 1
en r1 est largement dépendante de la position du négaton 2 : nous dirons que leurs
mouvements sont corrélés. Mais la nature des fermions fait que les négatons se repoussent à
cause de la répulsion colombienne. Le principe de Pauli ainsi que la corrélation inter-
électronique interdisent la présence des négatons de spin (𝜎 = ±1 /2) dans la même région
de l’espace au moyen du trou d’échange-corrélation (XC). Le terme « trou » représente la
diminution de la densité électronique en r2 liée aux interactions négaton-négaton non
classique. La probabilité d’occupation spatiale s’écrit alors :
𝜌2 (𝑟2 𝜎2 , 𝑟1 𝜎1 ) = 𝜌(𝑟2/1 )
Et se développe en :
𝑡𝑟𝑜𝑢 𝑡𝑟𝑜𝑢
𝜌(𝑟2/1 ) = 𝜌(𝑟1 )[𝜌(𝑟2 ) + 𝜌𝑥𝑐 ] = 𝜌(𝑟1 )𝜌(𝑟2 ) + 𝜌(𝑟𝑖 )𝜌𝑥𝑐
1 𝜌(𝑟1 , 𝑟2 )
𝑉𝑒𝑒 = ∫ ∫ 𝑑 3 𝑟1 𝑑 3 𝑟2
2 |𝑟1 − 𝑟2 |
1 𝜌(𝑟1 ,𝑟2 )
= 2 ∫ ∫ 𝑑 3 𝑟1 𝑑 3 𝑟2 |𝑟 −𝑟 |
1 2
( 𝑡𝑟𝑜𝑢
1 3 3
𝜌 𝑟1 )𝜌 𝑥𝑐 𝜌(𝑟1 , 𝑟2 )
+ ∫ ∫ 𝑑 𝑟1 𝑑 𝑟2
2 |𝑟1 − 𝑟2 |
𝑉𝑒𝑒 = 𝑈 + 𝐸𝑥𝑐
Dans leurs articles originaux, Kohn et Sham ont souligné le fait que l’on peut
considérer les solides très proches d’un gaz de négatons homogène. Dans cette limite, il est
soutenu que les effets d’échange-corrélation ont un caractère local. Les deux auteurs ont
proposé l’utilisation de l’approximation de la densité locale (LDA), dans laquelle l’énergie
𝐿𝐷𝐴
d’échange-Corrélation 𝐸𝑋𝐶 [𝜌] n’est autre qu’une intégrale sur tout l’espace.
En d’autres termes, elle repose sur l’hypothèse que les termes d’échange-corrélation
ne dépendent que de la valeur locale de 𝜌(r) ; c'est-à-dire qu’elle traite un système non
homogène comme étant localement homogène. L’énergie d’échange-corrélation s’exprime
alors de la manière suivante :
𝐿𝐷𝐴 [ ] 𝐿𝐷𝐴
𝐸𝑋𝐶 𝜌 = ∫ 𝜌(𝑟)𝜀𝑋𝐶 [𝜌(𝑟)]𝑑 3 r
𝐿𝐷𝐴
avec 𝜀𝑋𝐶 [𝜌(𝑟)] ont été paramétrées à partir des calculs Monte Carlo ,choisi comme
localement égal à celui d’un gaz de négatons homogène
𝐿𝐷𝐴 [ ( )] 𝑔𝑎𝑧−ℎ𝑜𝑚
𝜀𝑋𝐶 𝜌 𝑟 = 𝜀𝑋𝐶 [𝜌(𝑟)]
𝐿𝐷𝐴
𝜀𝑋𝐶 [𝜌(𝑟)] est souvent écrit comme somme de deux termes : un terme d’échange 𝜀𝑋𝐿𝐷𝐴 [𝜌(𝑟)]
et un terme de corrélation 𝜀 𝐿𝐷𝐴
𝐶 [𝜌(𝑟)]
𝐿𝐷𝐴 [ ( )]
𝜀𝑋𝐶 𝜌 𝑟 = 𝜀𝑋𝐿𝐷𝐴 [𝜌(𝑟)] + 𝜀𝐶𝐿𝐷𝐴 [𝜌(𝑟)]
𝐿𝐷𝐴
𝑉𝑋𝐶 = 𝛿𝐸𝑋𝐶 [𝜌(𝑟)]/𝛿𝜌(𝑟)
Kohn et Sham ont aussi utilisés cette approximation, pour les systèmes magnétique
polarisation de spin, ils ont renommé l’approximation de la densité locale de spin (LSPA) ou
l’énergie d’échange et de corrélation devient une fonctionnelle des deux densités de spin haut
et spin bas :
𝐿𝐷𝐴 [ 𝐿𝐷𝐴
𝐸𝑋𝐶 𝜌↓ , 𝜌↑ ] = ∫ 𝜌↓ (𝑟)𝜌↑ (𝑟)𝜀𝑋𝐶 [𝜌↓ (𝑟)𝜌↑ (𝑟)]𝑑 3 r
Par sa construction, on s’attend à ce que la LDA soit une bonne approximation pour
des densités qui varient lentement dans l’espace. Bien que cette condition ne soit pas toujours
vérifiée pour des systèmes réels, la LDA s’est avérée remarquablement appropriée pour une
grande variété de systèmes. Néanmoins, la LDA donne une très mauvaise estimation du gap
des isolants et des semi-conducteurs, ce qui n’est très surprenant, cette grandeur partiellement
de l’état fondamental. Elle permet, par exemple, de décrire la liaison covalente mais ne peut
pas en principe être utilisée pour étudier les systèmes dans lesquels les corrélations non
locales sont non négligeables. C’est le cas par exemple des systèmes dans lesquels les
interactions de Van der Walls sont importantes. Il a été également noté que les énergies de
cohésion des solides sont systématiques surestimées, et que l’erreur augmente au fur et à
mesure que la taille et/ou la dimensionnalité du système diminuent. En plus les distances de
liaison à l’équilibre sont souvent faiblement sous-estimées, mais l’erreur peut atteindre 10°/°
dans les petits systèmes.
Approximation de la GGA (Generalized Gradient Approximation)
𝐺𝐺𝐴 [ ] 𝐺𝐺𝐴
𝐸𝑋𝐶 𝜌 = ∫ 𝜌(𝑟)𝜀𝑋𝐶 [𝜌(𝑟) ∇𝜌 (𝑟)]𝑑 3 r
De façon générale, la GGA est amélioré par rapport à la LDA dans un certain nombre
de propriétés, comme l’énergie de liaison et l’énergie de cohésion. Par exemple, la LSDA a
montré que l’état fondamental du fer est non magnétique de structure à faces centrées. Ce
résultat est en contradiction avec les observations expérimentales qui montre que l’état
fondamental du fer est ferromagnétique de structure centrée (bcc). Par contre la GGA donne
l’état fondamental correct donné expérimentalement.