Por que resolver iterativamente as equações de Hartree-Fock resulta em convergência?

No método de campo auto-consistente de Hartree-Fock para resolver a equação eletrônica de Schroedinger independente do tempo, procuramos minimizar a energia do estado fundamental, , de um sistema de elétrons em um campo externo com relação à escolha dos orbitais de rotação, { χ i } .$E_{0}$$\{\chi_{i}\}$

Fazemos isso iterativamente resolvendo as equações de electrões onde x i é o spin / coordenada espacial de electrões i , ε é o valor próprio e f i é o operador Fock (um operador 1-electrão), com a forma f i = -

$$\hat{f}_{i}\chi(\mathbf{x}_{i})=\varepsilon\chi(\mathbf{x}_{i})$$ $\mathbf{x}_{i}$$i$$\varepsilon$$\hat{f}_{i}$ (as pistas soma sobre os núcleos, aqui, comZUmasendo a carga nuclear sobre um núcleo eriUmasendo a distância entre electrõesie núcleoA). V H F i é o potencial médio sentido pelo elétronidevido a todos os outros elétrons no sistema. ComoV H F i é dependente dos orbitais de rotação,χ $$\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}$, dos outros elétrons, podemos dizer que o operador Fock depende de suas próprias funções. Em "Modern Quantum Chemistry", de A. Szabo e N. Ostlund, pp. 54 (primeira edição), eles escrevem que "a equação de Hartree-Fock (2,52) é não-linear e deve ser resolvida iterativamente" . Estudei os detalhes dessa solução iterativa como parte de minha pesquisa, mas, para esta questão, acho que eles não são importantes, exceto para indicar a estrutura básica do método, que é:

1. Faça um palpite inicial dos orbitais giratórios e calcule V H F i .$\{\chi_{i}\}$$V_{i}^{\mathrm{HF}}$
2. Resolva a equação de autovalor acima para esses orbitais de rotação e obtenha novos orbitais de rotação.
3. Repita o processo com seus novos orbitais de rotação até que a consistência seja atingida.

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

Minha pergunta é a seguinte: como podemos saber que essa convergência ocorrerá? Por que as funções próprias das soluções iterativas sucessivas, em certo sentido, "melhoram" em direção a casos convergentes? Não é possível que a solução possa divergir? Não vejo como isso é evitado.

Como uma pergunta adicional, eu estaria interessado em saber por que as funções próprias convergentes (orbitais de rotação) fornecem a melhor (ou seja, mais baixa) energia do estado fundamental. Parece-me que a solução iterativa da equação de alguma forma tem convergência e minimização de energia "embutidas". Talvez haja alguma restrição embutida nas equações que assegure essa convergência?

Postagens cruzadas do Physics Stack Exchange: https://physics.stackexchange.com/q/20703/why-does-iteratively-solving-the-hartree-fock-equations-result-in-convergence

As equações de Hartree-Fock são o resultado da realização de minimização restrita de Newton-Raphson da energia em relação ao espaço de parâmetros dos determinantes de Slater (não tenho minha cópia de Szabo-Ostlund em mãos, mas acredito que isso seja apontado em a derivação). Portanto, o HF-SCF convergirá se o seu palpite inicial estiver em uma região convexa em torno de um mínimo. Em outros lugares, pode ou não convergir. A convergência SCF falha o tempo todo.

A teoria funcional da densidade (DFT) também usa uma abordagem de uma partícula semelhante à Hartree-Fock, embora o potencial efetivo esteja um pouco mais envolvido. Para atingir um mínimo global, o problema é abordado como um problema de ponto fixo não linear que, como disse Deathbreath , pode ser resolvido através de uma minimização restrita de Newton-Raphson . Uma abordagem comum na comunidade DFT é usar o método de Broyden que, se organizado corretamente ( J Phys A 17 (1984) L317 ), requer apenas dois vetores: a entrada e a saída atuais. (Veja Singh e Nordstrom , p. 91-92, para uma rápida visão geral desse método, ou Martin, Apêndice L, para uma visão geral mais completa das técnicas relacionadas.) Uma técnica mais recente usada no Wien2k tenta superar as dificuldades de convergência com o método Broyden, empregando um método multissecante ( PRB 78 (2008) 075114 , arXiv: 0801.3098 ).

Pode-se usar o algoritmo de amortecimento ideal ODA no ciclo SCF para obter um algoritmo de minimização real. Então sempre converge. (Artigos relacionados de Eric Cancès também merecem ser lidos.)