Direct Energy Minimization Based on Exponential Transformation in Density Functional Calculations of Finite and Extended Systems

Abstract

The energy minimization involved in density functional calculations of electronic systems can be carried out using an exponential transformation that preserves the orthonormality of the orbitals. The energy of the system is then represented as a function of the elements of a skew-Hermitian matrix that can be optimized directly using unconstrained minimization methods. An implementation based on the limited memory Broyden-Fletcher-Goldfarb-Shanno approach with inexact line search and a preconditioner is presented and the performance compared with that of the commonly used self-consistent field approach. Results are presented for the G2 set of 148 molecules, liquid water configurations with up to 576 molecules and some insulating crystals. A general preconditioner is presented that is applicable to systems with fractional orbital occupation as is, for example, needed in the k-point sampling for periodic systems. This exponential transformation direct minimization approach is found to outperform the standard implementation of the self-consistent field approach in that all the calculations converge with the same set of parameter values and it requires less computational effort on average. The formulation of the exponential transformation and the gradients of the energy presented here are quite general and can be applied to energy functionals that are not unitary invariant such as self-interaction corrected functionals.

References in corpus (2)

• Localized atomic basis set in the projector augmented wave method
• Steering Orbital Optimization out of Local Minima and Saddle Points Toward Lower Energy

Cited by in corpus (9)

• GPAW: An open Python package for electronic-structure calculations
• Orbital-optimized versus time-dependent density functional calculations of intramolecular charge transfer excited states
• Mn Dimer can be Described Accurately with Density Functional Calculations when Self-interaction Correction is Applied
• Economical Quasi-Newton Self Consistent Field Solver
• Freeze-and-release direct optimization method for variational calculations of excited electronic states
• Direct minimization on the complex Stiefel manifold in Kohn-Sham density functional theory for finite and extended systems
• Variational Density Functional Calculations of Excited States: Conical Intersection and Avoided Crossing in Ethylene Bond Twisting
• Method for Calculating Excited Electronic States Using Density Functionals and Direct Orbital Optimization with Real Space Grid or Plane Wave Basis Set
• Orbital Optimization and Neural-Network-Assisted Configuration Interaction Calculations of Rydberg States