Fonctionnelle hybride

Fonctionnelle hybride

Une fonctionnelle hybride est un élément d'une classe d'approximations à la fonctionnelle d'échange-corrélation, introduite en 1993 par Axel Becke[1], utilisé au sein de la théorie de la fonctionnelle de la densité (DFT). La caractéristique de ces fonctionnelles est d'avoir une partie d'échange basée sur la méthode de Hartree-Fock (HF) dépendante des orbitales alors que la partie de corrélation est basée sur une autre approche (soit issue de méthodes ab-initio, soit (semi-)empirique). Se basant sur un schéma simple d'amélioration[2], elles sont considérées comme « les fonctionnelles les plus précises disponibles tant que l'énergie est à prendre en compte et la méthode de prédilection de la communauté des chimistes »[3].

Sommaire

Origine

Les fonctionnelles hybrides visent à corriger l'erreur de self-interaction (interaction de l'électron avec lui-même) et la mauvaise représentation des phénomènes de corrélation non-locale des approximations « classiques » de la DFT (LDA et GGA), dans lesquelles elles ne sont pas prises en compte. Ces deux dernières conduisent à une délocalisation électronique trop importante, et donc à la mauvaise représentation des propriétés de certains matériaux comme les oxydes de transition. La construction des fonctionnelles hybrides est basée sur deux considérations. La première est que l'énergie d'échange (Ex) est plus importante que celle de corrélation (Ec), avec : Ex ≈ 10 % Ec (l'énergie d'échange-corrélation, Exc, est égale à Ex+Ec. La seconde est que si la DFT conduit à une surestimation de grandeurs et la méthode HF à une sous-estimation, alors la fonctionnelle exacte d'échange doit contenir une partie de l'échange HF[4]. La première fonctionnelle hybride proposée en se basant sur ces considérations est la fonctionnelle HH (pour half and half)[1], avec :

E_{xc}= \frac{1}{2}(E_x^{HF}+E_c^{DFT})

où \scriptstyle{E_c^{DFT}} indique l'énergie de corrélation obtenue par une fonctionnelle « classique » de la DFT (c'est-à-dire LDA ou GGA).

Fonctionnelles hybrides et connexion adiabatique

L'introduction d'un échange HF peut a priori se faire en substituant la partie d'échange des fonctionnelles par l'expression exacte de cet échange :

E_x^{HF}=\frac{1}{2}\sum_{i,j}\int\int\psi_i^*(\mathbf{r_1})\psi_j^*(\mathbf{r_1})\frac{1}{r_{12}}\psi_i(\mathbf{r_2})\psi_j(\mathbf{r_2})d\mathbf{r_1}d\mathbf{r_2},

où \scriptstyle{\psi} est une fonction d'onde, \scriptstyle{\mathbf{r}} un vecteur position et \scriptstyle{r_{12}} la distance entre le point 1 et le point, après substitution des fonctions d'onde classiques par les fonctions d'onde Kohn-Sham (KS) afin d'obtenir l'échange exact Kohn-Sham. Cette substitution a été proposée par Kohn et Sham dès 1965[réf. nécessaire] qui avaient relevé l'intérêt de l'introduction de l'échange exact, la partie de corrélation étant celle de la LDA. Bien que la forme utilisée convienne (en particulier, le potentiel effectif adopte un comportement asymptotique en 1/r à l'infini) et soit efficace pour des atomes seuls, elle s'avéra moins pertinente que la GGA pour des systèmes plus complexes dans la majorité des cas testés. Cette description erronée est imputable à la séparation artificielle des contributions des termes d'échange et de corrélation dans l'énergie d'échange-corrélation. En effet, cette configuration combine un trou d'échange non-local (contribution HF) avec un trou de corrélation local (contribution DFT-LDA) qui fait que la description locale du trou n'est plus possible.

Afin de résoudre ce type de problème, A. Becke a mis en avant[4] l'utilité de la formule de connexion adiabatique[5],[6],[7],[8] qui s'écrit :

E_{xc}=\int_0^1 U_{xc}^{\lambda}d\lambda

où \scriptstyle{\lambda} est un paramètre de couplage inter-électronique permettant de mettre en jeu de manière plus ou moins prononcée la répulsion coulombienne électronique en \scriptstyle{\frac{1}{r_{12}}} et \scriptstyle{U_{xc}^{\lambda}} l'énergie potentielle d'échange et corrélation à \scriptstyle{\lambda} donné. Cette formule permet de connecter le système KS sans interactions (\scriptstyle{\lambda}=0) incluant l'échange exact (les électrons sont considérés comme des fermions non chargés) au système réel (\scriptstyle{\lambda}=1) avec interactions de manière continue. Tous ces intermédiaires sont à même densité électronique, celle du système réel. La fonction \scriptstyle{U_{xc}^{\lambda}} pouvant être considérée comme une fonction linéaire de la variable \scriptstyle{\lambda} en première approximation, une solution simple[1] est donnée par :

E_{xc}^{HH}=\frac{1}{2}E_{x}^{HF}+\frac{1}{2}E_{c}^{LDA},

c'est-à-dire l'expression de l'énergie d'échange et corrélation de la fonctionnelle hybride « half-and-half ». Cette dernière a des performances comparables à la GGA aux tests G2[9],[10], et a de fait ouvert la voie aux fonctionnelles « adiabatiquement connectées » (adiabatically connected functionals en anglais).

Au-delà de la fonctionnelle hybride HH

Fonctionnelles hybrides paramétrées

Peu de temps après avoir proposé la fonctionnelle hybride HH, A. Becke, dans son article Density‐functional thermochemistry. III. The role of exact exchange[4] - l'article plus cité de la littérature scientifique en chimie en 2004[11] et figurant dans les 10 premiers de 1999 à 2005[12] - propose une fonctionnelle hybride paramétrée contenant des éléments de fonctionnelles « classiques », dont l'énergie d'échange-corrélation s'exprime de la manière suivante :

E_{xc}^{B3PW91}=E_{xc}^{LSDA}+a_0(E^{exact}_x-E^{LSDA}_x)+a_x\Delta E^{B88}_x+a_c\Delta E^{PW91}_c.

Dans cette expression, \scriptstyle{E_{xc}^{LSDA}} est l'approximation de densité de spin locale (LSDA) dont la partie de corrélation est celle proposée par J.P. Perdew et Y. Wang[13], \scriptstyle{E^{exact}_x} l'énergie d'échange exact, \scriptstyle{\Delta E^{B88}_x} la correction de gradient pour l'échange à la LSDA proposée par A. Becke en 1988[14] et \scriptstyle{\Delta E^{PW91}_c} la correction de gradient pour la corrélation de J.P. Perdew et Y. Wang de 1991. Les paramètres \scriptstyle{a_0}, \scriptstyle{a_x} et \scriptstyle{a_c} sont ajustés sur la base G1[15] par une méthode des moindres carrés. Leurs valeurs sont, respectivement de 0,20, 0,72 et 0,81. Cette fonctionnelle hybride, nommée B3PW91, permet d'améliorer les résultats obtenus dans le test G2 par rapport à la fonctionnelle hybride HH.

L'une des fonctionnelles hybrides les plus utilisées est directement dérivée de cette fonctionne B3PW91. La fonctionnelle B3LYP[16] est obtenue par substitution dans la formule de B3PW91 de la fonctionnelle de corrélation de Perdew-Wang par la fonctionnelle de corrélation de Lee-Yang-Parr[17] et de la fonctionnelle d'échange et corrélation LSDA par la fonctionnelle LDA de Vosko-Wilk-Nusair[18]. La B3LYP permet de réduire encore l'écart à la déviation chimique dans les tests G2. Cette fonctionnelle hybride a également été reparamétrée (en faisant varier la proportion d'échange exact) pour différentes applications. On citera par exemple la fonctionnelle B3LYP* de Reiher et collaborateurs[19] contenant 15 % d'échange HF et non plus 20 %.

Génération théorique de fonctionnelles hybrides

L'utilisation de la connexion adiabatique, et plus particulièrement du paramètre \scriptstyle{\lambda} de couplage inter-électronique, a été effectuée par J.P. Perdew et collaborateurs dans l'optique de fournir une fonctionnelle hybride dont la forme est étayée théoriquement, et donc se passant d'un paramétrage sur une base de données[3]. La fonctionnelle hybride PBE0[2], de forme :

E_{xc}=E_{x}^{DFT}+\frac{1}{4}(E_x^{HF}-E_{xc}^{DFT}),

dans laquelle la fonctionnelle d'échange et corrélation \scriptstyle{E_{x}^{DFT}} est la fonctionnelle GGA PBE[20] fut proposée en 1996, à partir de considérations sur la variation de \scriptstyle{U_{xc}^{\lambda}} avec \scriptstyle{\lambda}. J. Perdew et collaborateurs donnèrent une assise théorique pour les valeurs de paramètres et variations qui avaient été obtenus par ajustement sur des bases de données.

La fonctionnelle hybride PBE0 fournit de bons résultats sur un nombre important et varié de données[21],[22].

Liste de fonctionnelles hybrides

  • B1B95
  • B1LYP
  • B3PW91
  • B3LYP
  • B3LYP5
  • B3P86
  • BHandH
  • BHandHLYP
  • BMK
  • BP
  • EDF1
  • EDF2
  • M05
  • M05-2X
  • M06
  • M06-2X
  • M06-HF
  • M06-L
  • MPW1PW91
  • O3LYP
  • PBE1PBE
  • ωB97
  • ωB97X
  • ωB97X-DBP

Références et notes

Références et notes citées
  1. ↑ a, b et c (en) A.D. Becke, « A new mixing of Hartree-Fock and local density-functional theories », dans J. Chem. Phys., vol. 98, 1993, p. 1372–1377 [lien DOI] 
  2. ↑ a et b (en) John P. Perdew, Matthias Ernzerhof et Kieron Burke, « Rationale for mixing exact exchange with density functional approximations », dans J. Chem. Phys., vol. 105, 1996, p. 9982–9985 [texte intégral [PDF], lien DOI (pages consultées le 9 août 2010)] 
  3. ↑ a et b (en) Richard M. Martin, Electronic Structure : Basic Theory and Practical Methods, Cambridge University Press, 2004 (réimpr. 2005), 624 p. (ISBN 0 521 78285 6), p. 165-169 
  4. ↑ a, b et c (en) A.D. Becke, « Density‐functional thermochemistry. III. The role of exact exchange », dans J. Chem. Phys., vol. 98, 1993, p. 5648 (ISSN 0021-9606) [texte intégral, lien DOI] 
  5. ↑ (en) J. Harris et R. O. Jones, « The surface energy of a bounded electron gas », dans Journal of Physics F: Metal Physics, vol. 4, 1974, p. 1170 [texte intégral, lien DOI] 
  6. ↑ (en) O. Gunnarsson et B. I. Lundqvist, « Exchange and correlation in atoms, molecules, and solids by the spin-density-functional formalism », dans Phys. Rev. B, vol. 13, 1976, p. 4274 [texte intégral, lien DOI] 
  7. ↑ (en) David C. Langreth et John P. Perdew, « Exchange-correlation energy of a metallic surface: Wave-vector analysis », dans Phys. Rev. B, vol. 15, 1977, p. 2884 [texte intégral, lien DOI] 
  8. ↑ (en) J. Harris, « Adiabatic-connection approach to Kohn-Sham theory », dans Phys. Rev. A, vol. 29, 1984, p. 1648 [texte intégral, lien DOI] 
  9. ↑ (en) Oscar N. Ventura et Martina Kieninger, « Computational chemistry as an analytical tool: thermochemical examples in atmospheric chemistry », dans Pure & Appl. Chem, IUPAC, vol. 70, no 12, 1998, p. 2301 [texte intégral (page consultée le 9 août 2010)] 
  10. ↑ Les tests G2 sont une comparaison des résultats obtenus à la base de données utilisées par J. Pople et collaborateurs pour construire leur méthode Gaussian-2.
  11. ↑ (en) Most Cited Journal Articles 2004-Chemistry and Related, American Chemical Society. Consulté le 10 août 2010
  12. ↑ (en) CAS Science Spotlight Most Cited Chemistry and Related Science, American Chemical Society. Consulté le 10 août 2010
  13. ↑ (en) J. P. Perdew et Y. Wang, « Accurate and simple analytic representation of the electron-gas correlation energy », dans Phys. Rev. B, vol. 45, 1992, p. 13244 [texte intégral, lien DOI] .
  14. ↑ (en) A.D. Becke, « Density-functional exchange-energy approximation with correct asymptotic behavior », dans Phys. Rev. A, American Physical Society, vol. 38, no 6, 1988, p. 3098-3100 [texte intégral, lien DOI] 
  15. ↑ Base de données comprenant 56 énergies d'atomisation, 42 potentiels d'ionisation, 8 affinités protoniques et les énergies atomiques des 10 atomes des deux premières rangées du tableau périodique, ayant servi à J. Pople et collaborateurs pour construire leur méthode Gaussian-1.
  16. ↑ (en) P.J. Stephens, F. J. Devlin, C. F. Chabalowski and M. J. Frisch, « Ab Initio Calculation of Vibrational Absorption and Circular Dichroism Spectra Using Density Functional Force Fields », dans J. Phys. Chem., vol. 98, 1994, p. 11623–11627 [lien DOI] 
  17. ↑ (en) Chengteh Lee, Weitao Yang et Robert G. Parr, « Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density », dans Phys. Rev. B, American Physical Society, vol. 37, no 2, 1988, p. 785-789 [texte intégral, lien DOI] 
  18. ↑ (en) S. H. Vosko, L. Wilk et M. Nusair, « Accurate spin-dependent electron liquid correlation energies for local spin density calculations: a critical analysis », dans Can. J. Phys., vol. 58, 1980, p. 1200 [lien DOI] 
  19. ↑ (en) Markus Reiher, Oliver Salomon, Bernd Artur Hess, « Reparameterization of hybrid functionals based on energy differences of states of different multiplicity », dans Theoretical Chemistry Accounts: Theory, Computation, and Modeling (Theoretica Chimica Acta), Springer Berlin / Heidelberg, vol. 107, no 1, 2001, p. 48-55 (ISSN 1432-881X) [texte intégral, lien DOI] 
  20. ↑ (en) John P. Perdew, Kieron Burke, et Matthias Ernzerhof, « Generalized Gradient Approximation Made Simple », dans Phys. Rev. Lett., American Physical Society, vol. 77, no 18, 1996, p. 3865-3868 [texte intégral, lien DOI] 
  21. ↑ (en) Carlo Adamo, Gustavo E. Scuseria et Vincenzo Barone, « Accurate excitation energies from time-dependent density functional theory: Assessing the PBE0 model », dans The Journal of Chemical Physics, AIP, vol. 111, no 7, 1999, p. 2889-2899 [texte intégral, lien DOI] 
  22. ↑ (en) "Carlo Adamo, Maurizio Cossi, Giovanni Scalmani et Vincenzo Barone, « Accurate static polarizabilities by density functional theory: assessment of the PBE0 model », dans Chemical Physics Letters, vol. 307, no 3-4, 1999, p. 265 - 271 (ISSN 0009-2614) [texte intégral, lien DOI] 
Autres références utilisées
  • Cécil Combelles, Modélisation ab-initio appliquée à la conception de nouvelles batteries Li-Ion, Université Montpellier II - Sciences et Techniques du Languedoc, 2009, 208 p. [lire en ligne (page consultée le 8 août 2010)], p. 35-37.
    thèse de doctorat sous la direction de Marie-Liesse Doublet
     
  • Leonard Jeloaica, Etude Ab initio des mécanismes réactionnels dans la phase initiale du dépôt par couches atomiques des oxydes à moyenne et forte permittivité sur silicium, Université Paul Sabatier - Toulouse III, 2006, 141 p. [lire en ligne (page consultée le 8 août 2010)], p. 41.
    thèse de doctorat sous la direction de Mehdi Djafari Rouhani
     
  • Antony Fouqueau, Étude théorique de complexes de Fer par la théorie de la fonctionnelle de la densité: Application au problème de piégeage de spin, Université Joseph-Fourier - Grenoble I, 2005, 264 p. [lire en ligne (page consultée le 8 août 2010)], p. 65-66.
    thèse de doctorat sous la direction de Mark E. Casida
     
  • Fayçal Allouti, Chemins réactionnels conduisant à la formation des oxydes des métaux de transition 3dn (n>5) : structure électronique des composés antiferromagnétiques M2O2 (M=Cr, Mn, Fe, Co, Ni, Cu), Université Pierre et Marie Curie - Paris VI, 2006, 214 p. [lire en ligne (page consultée le 8 août 2010)], p. 48-50.
    thèse de doctorat sous la direction de Mohammad Esmaïl Alikhani
     

Voir aussi



Wikimedia Foundation. 2010.

Contenu soumis à la licence CC-BY-SA. Source : Article Fonctionnelle hybride de Wikipédia en français (auteurs)

Игры ⚽ Поможем сделать НИР

Regardez d'autres dictionnaires:

  • Théorie de la fonctionnelle de la densité — Pour les articles homonymes, voir DFT. La théorie de la fonctionnelle de la densité (pour Density Functional Theory, sous entendu électronique : DFT) constitue au début du XXIe siècle l une des méthodes les plus utilisées dans les… …   Wikipédia en Français

  • Theorie de la fonctionnelle de la densite — Théorie de la fonctionnelle de la densité Méthodes numériques pour le calcul de la structure électronique Hartree Fock Théorie de la perturbation de Møller Plesset Interaction de configuration Méthode du cluster couplé Champ multi… …   Wikipédia en Français

  • Théorie de la fonctionnelle densité — Théorie de la fonctionnelle de la densité Méthodes numériques pour le calcul de la structure électronique Hartree Fock Théorie de la perturbation de Møller Plesset Interaction de configuration Méthode du cluster couplé Champ multi… …   Wikipédia en Français

  • Méthode composite de chimie quantique — Les méthodes composites de chimie quantique sont des méthodes de chimie numérique dont le but est d obtenir une grande précision par combinaison des résultats de plusieurs calculs. Cette combinaison est constituée à partir de méthodes d un haut… …   Wikipédia en Français

  • FUSIONS CELLULAIRES — Depuis le début des années soixante, il est possible de fusionner des cellules somatiques entre elles. Ce mode de croisement est bien distinct de celui qui s’opère lors de la reproduction sexuelle. Il constitue néanmoins une approche alternative… …   Encyclopédie Universelle

  • Proteomique — Protéomique La protéomique désigne la science qui étudie les protéomes, c est à dire l ensemble des protéines d une cellule, organelle, tissu, organe ou organisme à un moment donné et sous des conditions données. Dans la pratique, la protéomique… …   Wikipédia en Français

  • Protéomique — La protéomique désigne la science qui étudie les protéomes, c est à dire l ensemble des protéines d une cellule, d un organite, d un tissu, d un organe ou d un organisme à un moment donné et sous des conditions données. Dans la pratique, la… …   Wikipédia en Français

  • Saab 9-X BioHybrid — Saab 9X BioHybrid Constructeur Saab Classe Compacte Usines d’assemblage Trollhättan …   Wikipédia en Français

  • Champ de force (chimie) — Pour les articles homonymes, voir champ de force. Un champ de force peut par exemple être utilisé afin de minimiser l énergie d étirement de cette molécule d éthane. Dans le cadre de la m …   Wikipédia en Français

  • CONCEPTION ASSISTÉE PAR ORDINATEUR — La conception assistée par ordinateur, ou C.A.O., peut être considérée comme l’ensemble des aides informatiques aux bureaux d’études et aux bureaux des méthodes. Elle concerne donc le processus conception fabrication depuis l’élaboration du… …   Encyclopédie Universelle

Share the article and excerpts

Direct link
… Do a right-click on the link above
and select “Copy Link”