¿Por qué la resolución iterativa de las ecuaciones de Hartree-Fock resulta en convergencia?

En el método de campo autoconsistente de Hartree-Fock para resolver la ecuación electrónica de Schroedinger independiente del tiempo, buscamos minimizar la energía del estado fundamental, , de un sistema de electrones en un campo externo con respecto a la elección de los orbitales giratorios, { χ i } .$E_{0}$$\{\chi_{i}\}$

Hacemos esto mediante la resolución iterativa los 1-electrón ecuaciones donde x i es el giro / espacial de coordenadas del electrón i , ε es el valor propio orbital y f i es el operador de Fock (un operador de 1-electrón), con la forma f i = -

$$\hat{f}_{i}\chi(\mathbf{x}_{i})=\varepsilon\chi(\mathbf{x}_{i})$$ $\mathbf{x}_{i}$$i$$\varepsilon$$\hat{f}_{i}$ (la suma se ejecuta sobre los núcleos, aquí, conZAsiendo la carga nuclear en el núcleo A yriAsiendo la distancia entre el electróniy el núcleoA). V H F i es el potencial promedio que siente el electrónidebido a todos los demás electrones del sistema. ComoV H F i depende de los orbitales giratorios,χ $$\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}$, de los otros electrones, podemos decir que el operador Fock depende de sus funciones propias. En "Modern Quantum Chemistry" de A. Szabo y N. Ostlund, págs. 54 (la primera edición) escriben que "la ecuación de Hartree-Fock (2.52) no es lineal y debe resolverse de forma iterativa" . He estudiado los detalles de esta solución iterativa como parte de mi investigación, pero para esta pregunta creo que no son importantes, excepto para establecer la estructura básica del método, que es:

1. Haga una conjetura inicial de los orbitales giratorios, y calcule V H F i .$\{\chi_{i}\}$$V_{i}^{\mathrm{HF}}$
2. Resuelva la ecuación de valor propio anterior para estos orbitales giratorios y obtenga nuevos orbitales giratorios.
3. Repita el proceso con sus nuevos orbitales giratorios hasta alcanzar la autoconsistencia.

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

Mi pregunta es esta: ¿cómo podemos saber que ocurrirá esta convergencia? ¿Por qué las funciones propias de las soluciones iterativas sucesivas en algún sentido "mejoran" hacia el caso convergente? ¿No es posible que la solución pueda divergir? No veo cómo se previene esto.

Como otra pregunta, me interesaría saber por qué las funciones propias convergentes (orbitales giratorios) dan la mejor (es decir, la más baja) energía del estado fundamental. Me parece que la solución iterativa de la ecuación de alguna manera tiene convergencia y minimización de energía "incorporada". ¿Quizás hay alguna restricción incorporada en las ecuaciones que asegura esta convergencia?

Publicación cruzada de Physics Stack Exchange: https://physics.stackexchange.com/q/20703/why-does-iteratively-solving-the-hartree-fock-equations-result-in-convergence

Las ecuaciones de Hartree-Fock son el resultado de realizar una minimización restringida de Newton-Raphson de la energía con respecto al espacio de parámetros de los determinantes Slater (no tengo a mano mi copia de Szabo-Ostlund, pero creo que esto se señala en la derivación) Por lo tanto, HF-SCF convergerá si su conjetura inicial está en una región convexa alrededor de un mínimo. En otros lugares, puede o no converger. La convergencia SCF falla todo el tiempo.

La teoría funcional de la densidad (DFT) también utiliza un enfoque de una partícula similar a Hartree-Fock, aunque el potencial efectivo es un poco más complicado. Para lograr un mínimo global, el problema se aborda como un problema de punto fijo no lineal que, como dijo Deathbreath , se puede resolver a través de una minimización restringida de Newton-Raphson . Un enfoque común en la comunidad DFT es utilizar el Método de Broyden, que si se organiza correctamente ( J Phys. A 17 (1984) L317 ) requiere solo dos vectores: la entrada y la salida actuales. (Ver Singh y Nordstrom , p. 91-92, para una descripción rápida de este método, o Martin, Apéndice L, para una descripción más completa de las técnicas relacionadas.) Una técnica más reciente utilizada en Wien2k intenta superar las dificultades de convergencia con el método Broyden mediante el uso de un método multi-secante. ( PRB 78 (2008) 075114 , arXiv: 0801.3098 )

Se puede usar el algoritmo de amortiguación óptimo ODA en el ciclo SCF para obtener un algoritmo de minimización real. Entonces siempre converge. (También vale la pena leer los documentos relacionados de Eric Cancès).