Pourquoi la résolution itérative des équations de Hartree-Fock entraîne-t-elle une convergence?

Dans la méthode de champ auto-cohérent de Hartree-Fock pour résoudre l'équation électronique de Schroedinger indépendante du temps, nous cherchons à minimiser l'énergie d'état fondamental, , d'un système d'électrons dans un champ externe par rapport au choix des orbitales de spin, { χ i } .$E_{0}$$\{\chi_{i}\}$

Nous faisons cela en résolvant de manière itérative les équations 1 à électrons où x i est le spin / coordonnée spatiale des électrons i , ε est la valeur propre orbital et f i est l'opérateur Fock (un opérateur d'électrons 1), avec la forme f i = -

$$\hat{f}_{i}\chi(\mathbf{x}_{i})=\varepsilon\chi(\mathbf{x}_{i})$$ $\mathbf{x}_{i}$$i$$\varepsilon$$\hat{f}_{i}$ (la sommation s'étend ici sur les noyaux,ZAétant la charge nucléaire sur le noyau A etriAétant la distance entre l'électroniet le noyauA). V H F i est le potentiel moyen ressenti par l'électronidû à tous les autres électrons du système. PuisqueV H F i dépend des orbitales de spin,χ $$\hat{f}_{i} = -\frac{1}{2}\nabla^{2}_{i}-\sum_{A=1}^{M}\frac{Z_{A}}{r_{iA}}+V^{\mathrm{HF}}_{i}$$ $Z_{A}$$r_{iA}$$i$$A$$V^{\mathrm{HF}}_{i}$$i$$V_{i}^{\mathrm{HF}}$$\chi_{j}$, des autres électrons, on peut dire que l'opérateur de Fock dépend de ses fonctions propres. Dans "Modern Quantum Chemistry" de A. Szabo et N. Ostlund, pp. 54 (la première édition), ils écrivent que "l'équation de Hartree-Fock (2.52) est non linéaire et doit être résolue de manière itérative" . J'ai étudié les détails de cette solution itérative dans le cadre de mes recherches, mais pour cette question, je pense qu'ils sont sans importance, sauf pour énoncer la structure de base de la méthode, qui est:

1. Faites une première estimation des spin-orbitales, et calculez V H F i .$\{\chi_{i}\}$$V_{i}^{\mathrm{HF}}$
2. Résolvez l'équation des valeurs propres ci-dessus pour ces orbitales de spin et obtenez de nouvelles orbitales de spin.
3. Répétez le processus avec vos nouvelles orbitales de spin jusqu'à ce que l'auto-cohérence soit atteinte.

$V_{i}^{\mathrm{HF}}$

Ma question est la suivante: comment savoir que cette convergence se produira? Pourquoi les fonctions propres des solutions itératives successives "s'améliorent" en quelque sorte vers le cas convergent? N'est-il pas possible que la solution diverge? Je ne vois pas comment cela est empêché.

Comme autre question, je serais intéressé de savoir pourquoi les fonctions propres convergentes (orbitales de spin) donnent la meilleure (c'est-à-dire la plus basse) énergie d'état fondamental. Il me semble que la solution itérative de l'équation a en quelque sorte une convergence et une minimisation d'énergie "intégrées". Peut-être y a-t-il une contrainte intégrée dans les équations qui assure cette convergence?

Post-cross du Physics Stack Exchange: https://physics.stackexchange.com/q/20703/why-does-iterativement-solving-the-hartree-fock-equations-result-in-convergence

Les équations de Hartree-Fock sont le résultat d'une minimisation contrainte de Newton-Raphson de l'énergie par rapport à l'espace des paramètres des déterminants de Slater (je n'ai pas ma copie de Szabo-Ostlund à portée de main, mais je crois que cela est souligné dans la dérivation). Par conséquent, HF-SCF convergera si votre estimation de départ se situe dans une région convexe autour d'un minimum. Ailleurs, il peut ou non converger. La convergence SCF échoue tout le temps.

La théorie fonctionnelle de la densité (DFT) utilise également une approche à une particule similaire à Hartree-Fock, bien que le potentiel effectif soit un peu plus impliqué. Pour atteindre un minimum global, le problème est abordé comme un problème de point fixe non linéaire qui, comme l' a dit Deathbreath , peut être résolu via une minimisation de Newton-Raphson contrainte . Une approche courante dans la communauté DFT consiste à utiliser la méthode de Broyden qui, si elle est correctement organisée ( J Phys A 17 (1984) L317 ), ne nécessite que deux vecteurs: l'entrée et la sortie actuelles. (Voir Singh et Nordstrom , p. 91-92, pour un aperçu rapide de cette méthode, ou Martin, Annexe L, pour un aperçu plus complet des techniques connexes.) Une technique plus récente utilisée dans Wien2k tente de surmonter les difficultés de convergence avec la méthode de Broyden en utilisant une méthode multi-sécante ( PRB 78 (2008) 075114 , arXiv: 0801.3098 ).

On peut utiliser l'algorithme d'amortissement optimal ODA dans le cycle SCF pour obtenir un algorithme de minimisation réel. Ensuite, il converge toujours. (Les articles connexes d'Eric Cancès méritent également d'être lus.)