Convergence analysis of direct minimization and self-consistent iterations

Abstract

This article is concerned with the numerical solution of subspace optimization problems, consisting of minimizing a smooth functional over the set of orthogonal projectors of fixed rank. Such problems are encountered in particular in electronic structure calculation (Hartree-Fock and Kohn-Sham Density Functional Theory -DFT- models). We compare from a numerical analysis perspective two simple representatives, the damped self-consistent field (SCF) iterations and the gradient descent algorithm, of the two classes of methods competing in the field: SCF and direct minimization methods. We derive asymptotic rates of convergence for these algorithms and analyze their dependence on the spectral gap and other properties of the problem. Our theoretical results are complemented by numerical simulations on a variety of examples, from toy models with tunable parameters to realistic Kohn-Sham computations. We also provide an example of chaotic behavior of the simple SCF iterations for a nonquadratic functional.

References in corpus (12)

• Separable Dual Space Gaussian Pseudo-potentials
• Mathematical theory and numerical methods for Bose-Einstein condensation
• Ensemble density-functional theory for ab-initio molecular dynamics of metals and finite-temperature insulators
• Efficient spectral computation of the stationary states of rotating Bose-Einstein condensates by the preconditioned nonlinear conjugate gradient method
• There are No Unfilled Shells in Hartree-Fock Theory
• Computation of Ground States of the Gross-Pitaevskii Functional via Riemannian Optimization
• Robust Mixing for Ab-Initio Quantum Mechanical Calculations
• Computing the Self-Consistent Field in Kohn-Sham Density Functional Theory
• A Projected Preconditioned Conjugate Gradient Algorithm for Computing Many Extreme Eigenpairs of a Hermitian Matrix
• Difficulties in applying pure Kohn-Sham density functional theory electronic structure methods to protein molecules
• Convergence analysis of adaptive DIIS algorithms with application to electronic ground state calculations
• On the Nature of Self-Consistency in Density Functional Theory