A Proximal Gradient Method for Ensemble Density Functional Theory

Michael Ulbrich,Zaiwen Wen,Chao Yang,Dennis Klöckner,Zhaosong Lu
DOI: https://doi.org/10.1137/14098973X
IF: 2.968
2015-01-01
SIAM Journal on Scientific Computing
Abstract:The ensemble density functional theory (E-DFT) is valuable for simulations of metallic systems due to the absence of a gap in the spectrum of the Hamiltonian matrices. Although the widely used self-consistent field (SCF) iteration method can be extended to solve the minimization of the total energy functional with respect to orthogonality constraints, there is no theoretical guarantee on the convergence of these algorithms. In this paper, we consider an equivalent model with a single variable and a single spherical constraint by eliminating the dependence on the fractional occupancies. A proximal gradient method is developed by keeping the entropy term but linearizing all other terms in the total energy functional. Convergence to the stationary point is established. Numerical results using the MATLAB toolbox KSSOLV can outperform SCF consistently on many metallic systems.
What problem does this paper attempt to address?