Direct minimization on the complex Stiefel manifold in Kohn-Sham density functional theory for finite and extended systems

Abstract

Direct minimization method on the complex Stiefel manifold in Kohn-Sham density functional theory is formulated to treat both finite and extended systems in a unified manner. This formulation is well-suited for scenarios where straightforward iterative diagonalization becomes challenging, especially when the Aufbau principle is not applicable. We present the theoretical foundation and numerical implementation of the Riemannian conjugate gradient (RCG) within a localized non-orthogonal basis set. Riemannian Broyden-Fletcher-Goldfarb-Shanno (RBFGS) method is tentatively implemented. Extensive testing compares the performance of the proposed methods and highlights that the quasi-Newton method is more efficient. However, for extended systems, the computational time required grows rapidly with respect to the number of $\mathbf{k}$-points.

20 pages, 5 figures

References in corpus (8)

• Optimization Algorithm for the Generation of ONCV Pseudopotentials
• Multiresolution analysis of electronic structure: semicardinal and wavelet bases
• DFT-FE -- A massively parallel adaptive finite-element code for large-scale density functional theory calculations
• Reduced Density Matrix Functional for Many-Electron Systems
• Systematically improvable optimized atomic basis sets for {\it ab inito} calculations
• New Algebraic Formulation of Density Functional Calculation
• Steering Orbital Optimization out of Local Minima and Saddle Points Toward Lower Energy
• Direct Energy Minimization Based on Exponential Transformation in Density Functional Calculations of Finite and Extended Systems