Calcul Hartree Fock à la main
Estimez pas à pas l’énergie d’un système héliumoïde 1s² avec l’approximation Hartree-Fock minimale. Cet outil pédagogique utilise un paramètre orbital effectif zeta, sépare les contributions cinétique, attraction noyau-électron et répulsion électronique, puis trace la courbe d’énergie en fonction de zeta.
Comprendre le calcul Hartree Fock à la main
Le calcul Hartree Fock à la main est souvent présenté comme une étape charnière entre la mécanique quantique élémentaire et la chimie quantique computationnelle moderne. En pratique, il ne s’agit pas de refaire à la main l’ensemble d’un programme de structure électronique complet, mais plutôt de comprendre la logique de l’approximation du champ moyen, la construction des orbitales, l’énergie électronique et l’auto-cohérence. Pour un étudiant, un enseignant ou un praticien qui souhaite réviser les fondements, ce type de calcul permet d’identifier très clairement ce qui relève de l’approximation physique et ce qui relève ensuite de l’algorithme numérique.
Dans sa forme générale, la méthode Hartree Fock suppose que la fonction d’onde électronique peut être approchée par un seul déterminant de Slater. Chaque électron évolue dans un champ moyen créé par le noyau et les autres électrons. Cela remplace un problème de corrélation électronique à N corps, insoluble analytiquement dans le cas général, par un ensemble d’équations orbitalaires couplées. À la main, on simplifie énormément le cadre pour rendre les intégrales manipulables. C’est pour cela que les exemples les plus fréquents portent sur l’atome d’hélium, les ions héliumoïdes ou les systèmes minimaux avec très peu de fonctions de base.
Pourquoi apprendre Hartree Fock sans ordinateur
Apprendre Hartree Fock à la main présente plusieurs avantages pédagogiques. D’abord, cela clarifie la structure de l’Hamiltonien électronique non relativiste. Ensuite, cela permet de voir d’où viennent les termes d’énergie. Enfin, cela évite l’effet de boîte noire fréquemment rencontré lorsqu’on se contente de lancer un calcul dans un logiciel sans comprendre les matrices ni les intégrales sous-jacentes.
- Vous comprenez le rôle de l’énergie cinétique électronique.
- Vous visualisez l’attraction noyau-électron et sa dépendance en Z.
- Vous distinguez la répulsion coulombienne moyenne et les effets d’échange.
- Vous saisissez la logique du principe variationnel.
- Vous voyez pourquoi l’énergie Hartree Fock reste au-dessus de l’énergie exacte.
Le cas le plus simple: l’ion héliumoïde 1s²
Le calculateur ci-dessus se concentre sur un système héliumoïde à deux électrons dans une orbitale spatiale de type 1s. C’est un cas classique, car il permet de condenser l’idée de Hartree Fock en une seule variable, le paramètre effectif zeta, qui joue le rôle d’une charge nucléaire ressentie par les électrons après écrantage partiel. Avec une orbitale de type Slater 1s normalisée, on obtient une expression analytique de l’énergie totale. Elle est suffisamment simple pour être dérivée, minimisée et interprétée à la main.
Pour cette approximation, l’énergie totale en unités atomiques s’écrit:
E(zeta) = zeta² – 2Zzeta + 5zeta/8
Les trois blocs physiques sont immédiatement visibles:
- zeta² représente la contribution cinétique totale des deux électrons.
- -2Zzeta correspond à l’attraction noyau-électron.
- 5zeta/8 représente la répulsion électronique moyenne dans cette base très simple.
Le minimum énergétique est obtenu en annulant la dérivée:
dE/dzeta = 2zeta – 2Z + 5/8 = 0, donc zeta optimum = Z – 5/16.
Exemple rapide pour l’hélium
Pour l’hélium, Z = 2. On trouve alors zeta optimum = 2 – 0,3125 = 1,6875. Si on remplace cette valeur dans l’expression de l’énergie, on obtient une estimation proche de l’énergie Hartree Fock minimale dans ce modèle à un seul paramètre. Cette énergie n’est pas l’énergie exacte de l’hélium, mais elle a une valeur pédagogique considérable, car elle montre comment l’écrantage électronique réduit la charge nucléaire effectivement ressentie par chaque électron.
Étapes d’un calcul Hartree Fock à la main
Même si les détails changent selon le système étudié, la logique générale reste la même. Voici la méthode canonique que l’on suit habituellement.
- Écrire l’Hamiltonien électronique dans l’approximation de Born Oppenheimer, en séparant énergie cinétique, attraction noyau-électron et répulsion interélectronique.
- Choisir une forme d’orbitales d’essai, souvent des orbitales de Slater ou une base minimale gaussienne dans les traitements plus modernes.
- Construire le déterminant de Slater afin d’assurer l’antisymétrie de la fonction d’onde.
- Calculer l’énergie moyenne à l’aide des intégrales une et deux électrons.
- Appliquer le principe variationnel pour minimiser l’énergie par rapport aux paramètres libres ou aux coefficients de la base.
- Interpréter les orbitales et l’énergie en séparant contributions directes, échange et éventuelle erreur due à l’absence de corrélation électronique.
Que signifie vraiment “à la main” en chimie quantique
Dans l’usage universitaire, “faire Hartree Fock à la main” signifie souvent l’une des trois choses suivantes. Premièrement, dériver l’énergie variationnelle pour un atome simple. Deuxièmement, établir les équations de Roothaan pour une très petite base. Troisièmement, résoudre une ou deux itérations auto-cohérentes sur papier avec des nombres simplifiés. Le but n’est donc pas l’exhaustivité numérique, mais la transparence conceptuelle.
Version analytique contre version matricielle
Dans les systèmes atomiques les plus simples, on travaille directement avec un paramètre comme zeta. Dans les molécules ou les bases gaussiennes, on passe à la forme matricielle: matrice de recouvrement S, matrice du cœur H, matrices coulombiennes et d’échange, puis diagonalisation de la matrice de Fock. La version analytique est excellente pour construire l’intuition. La version matricielle est indispensable dès que le problème se complexifie.
| Méthode | Principe | Coût conceptuel | Précision typique | Usage pédagogique |
|---|---|---|---|---|
| Hartree simple | Champ moyen sans antisymétrie complète | Faible | Limitée | Introduction historique |
| Hartree Fock | Déterminant de Slater et échange exact au niveau mono-déterminant | Moyen | Bonne base qualitative | Référence conceptuelle essentielle |
| Post-HF | Ajout explicite de corrélation électronique | Élevé | Très élevée selon la méthode | Niveau avancé |
| DFT | Fonctionnelle de la densité | Moyen | Souvent excellente en pratique | Comparaison moderne |
Statistiques et ordres de grandeur utiles
Pour donner un cadre quantitatif, il est utile de rappeler quelques nombres standards. Un hartree vaut environ 27,2114 eV, ce qui permet de convertir directement les résultats de nombreux calculs atomiques vers une unité plus familière en spectroscopie ou en chimie. Pour l’hélium, l’énergie non relativiste exacte du fondamental est proche de -2,9037 hartree, alors que la limite Hartree Fock est aux environs de -2,8617 hartree. La différence, d’environ 0,042 hartree, est l’énergie de corrélation. Elle paraît petite en valeur absolue, mais elle est conceptuellement fondamentale, car elle quantifie ce que Hartree Fock ne capture pas dans une description mono-déterminant.
| Système | Valeur ou statistique | Ordre de grandeur | Commentaire |
|---|---|---|---|
| 1 hartree | 27,2114 eV | Échelle énergétique atomique | Conversion standard utilisée dans les logiciels et les cours |
| Hélium exact non relativiste | Environ -2,9037 Eh | Référence benchmark | Valeur modèle pour comparer les méthodes approchées |
| Hélium Hartree Fock limite | Environ -2,8617 Eh | Très proche mais non exact | Montre l’importance de la corrélation électronique |
| Énergie de corrélation de He | Environ -0,042 Eh | Près de -1,14 eV | Petite en apparence, décisive en précision chimique |
Interprétation du paramètre zeta
Le paramètre zeta peut être compris comme une mesure de compaction orbitale. Quand zeta augmente, l’orbitale 1s se contracte autour du noyau. Cela augmente en général l’énergie cinétique, mais renforce aussi l’attraction noyau-électron. La répulsion entre électrons augmente également, car les densités se superposent davantage. Le minimum de la courbe E(zeta) traduit l’équilibre entre ces effets concurrents. C’est exactement le type de compromis que l’approximation variationnelle sait identifier avec élégance.
Pourquoi la courbe d’énergie a un minimum
Si zeta devient trop petit, les électrons sont trop diffus et n’exploitent pas assez l’attraction du noyau. Si zeta devient trop grand, la contraction excessive pénalise le système via l’énergie cinétique et la répulsion mutuelle. Le résultat est une courbe convexe, avec un point bas bien défini. Le graphe du calculateur rend ce comportement visible en temps réel.
Différence entre énergie totale et énergie orbitale
Une confusion fréquente chez les débutants consiste à assimiler la somme des énergies orbitalaires à l’énergie totale du système. En Hartree Fock, cette identité n’est pas vraie directement, car les interactions électron-électron seraient alors comptées deux fois. Il faut toujours distinguer la valeur propre de l’opérateur de Fock associée à une orbitale et l’énergie électronique totale calculée à partir de l’expression complète. Dans le cas très simple traité ici, nous affichons surtout l’énergie totale et ses composantes physiques, ce qui est généralement le plus clair pour un apprentissage progressif.
Erreurs fréquentes dans un calcul Hartree Fock à la main
- Oublier de normaliser correctement l’orbitale d’essai.
- Confondre l’énergie d’une orbitale avec l’énergie totale du système.
- Négliger l’antisymétrie et donc les termes d’échange dans les cas où ils sont pertinents.
- Utiliser une unité incohérente entre hartree, rydberg et électron-volt.
- Perdre le facteur 2 lié au nombre d’électrons dans une couche fermée.
- Conclure trop vite que Hartree Fock est “faux” alors qu’il s’agit d’une approximation systématique, souvent excellente comme point de départ.
Comment relier ce calcul à la pratique moderne
Les programmes de chimie quantique modernes utilisent des bases gaussiennes, des algorithmes SCF robustes, des stratégies de convergence avancées et parfois des traitements relativistes. Pourtant, la logique reste la même que dans le petit calcul à la main: on écrit une énergie, on construit une matrice de Fock, on résout pour obtenir des orbitales, puis on boucle jusqu’à auto-cohérence. Maîtriser la version simplifiée rend la version logicielle beaucoup plus intelligible.
Quand Hartree Fock suffit-il
Hartree Fock est souvent qualitativement très utile pour:
- obtenir une première structure orbitale,
- préparer un calcul post-HF,
- interpréter des tendances périodiques,
- enseigner l’échange et l’antisymétrie.
En revanche, pour des énergies de réaction fines, des interactions faibles, des états excités délicats ou des systèmes fortement corrélés, il faut généralement aller au-delà.
Ressources académiques et institutionnelles fiables
Pour approfondir la théorie et vérifier les constantes physiques, vous pouvez consulter des sources institutionnelles reconnues:
- NIST Physics Laboratory pour les constantes, unités atomiques et références physiques.
- LibreTexts Chemistry pour des explications universitaires détaillées sur Hartree, Hartree Fock et la chimie quantique.
- NIST Chemistry WebBook pour des données thermochimiques et spectroscopiques utiles en mise en contexte.
Comment utiliser intelligemment le calculateur ci-dessus
Commencez avec Z = 2 pour l’hélium. Laissez l’option optimisée activée, puis observez le zeta obtenu et l’énergie totale. Ensuite, passez en mode manuel et entrez des valeurs de zeta plus petites et plus grandes pour voir comment la courbe réagit. Faites enfin varier Z vers 3 ou 4 pour représenter Li+ ou Be2+. Vous verrez que l’orbital se contracte à mesure que la charge nucléaire augmente, et que l’énergie totale devient plus négative, ce qui traduit une liaison électronique plus forte au noyau.
Cette exploration développe une intuition très précieuse: en chimie quantique, les formules ne servent pas seulement à donner un nombre final. Elles décrivent des compromis physiques entre localisation, attraction coulombienne, répulsion interélectronique et structure antisymétrique de la fonction d’onde. C’est exactement cette intuition qui distingue un simple utilisateur de logiciel d’un vrai praticien de la structure électronique.