Abstract:Variational calculations of excited electronic states are carried out by finding saddle points on the surface that describes how the energy of the system varies as a function of the electronic degrees of freedom. This approach has several advantages over commonly used methods especially in the context of density functional calculations, as collapse to the ground state is avoided and yet, the orbitals are variationally optimized for the excited state. This optimization makes it possible to describe excitations with large charge transfer where calculations based on ground state orbitals are problematic, as in linear response time-dependent density functional theory. A generalized mode following method is presented where an $n^{\text{th}}$-order saddle point is found by inverting the components of the gradient in the direction of the eigenvectors of the $n$ lowest eigenvalues of the electronic Hessian matrix. This approach has the distinct advantage of following a chosen excited state through atomic configurations where the symmetry of the single determinant wave function is broken, as demonstrated in calculations of potential energy curves for nuclear motion in the ethylene and dihydrogen molecules. The method is implemented using a generalized Davidson algorithm and an exponential transformation for updating the orbitals within a generalized gradient approximation of the energy functional. Convergence is found to be more robust than for a direct optimization approach previously shown to outperform standard self-consistent field approaches, as illustrated here for charge transfer excitations in nitrobenzene and N-phenylpyrrole, involving calculations of $4^{\text{th}}$- and $6^{\text{th}}$-order saddle points, respectively. Finally, calculations of a diplatinum and silver complex are presented, illustrating the applicability of the method to excited state energy curves of large molecules.

Gradient Type Optimization Methods for Electronic Structure Calculations

Self-Consistent Solution Of Kohn-Sham Equation By Real-Space Finite-Difference Method

A Proximal Gradient Method for Ensemble Density Functional Theory

Efficient simulations of Hartree--Fock equations by an accelerated gradient descent method

Convergence of density functional iterative procedures with a Newton-Raphson algorithm

Adaptive Regularized Self-Consistent Field Iteration with Exact Hessian for Electronic Structure Calculation.

A Convergence Analysis of a Structure-Preserving Gradient Flow Method for the All-Electron Kohn-Sham Model

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

Excited state orbital optimization via minimizing the square of the gradient: General approach and application to singly and doubly excited states via density functional theory

Convergence Analysis of Gradient Algorithms on Riemannian Manifolds Without Curvature Constraints and Application to Riemannian Mass

A Time-Dependent Random State Approach for Large-Scale Density Functional Calculations

Efficient Treatment of Relativistic Effects with Periodic Density Functional Methods: Energies, Gradients, and Stress Tensors

Efficient numerical methods for computing the stationary states of phase field crystal models

An Efficient Implementation of Analytical Nuclear Gradients for Linear-Response Time-Dependent Density Functional Theory in the Plane Wave Basis

Structure Optimization with Stochastic Density Functional Theory

Projection Gradient Method for Energy Functional Minimization with A Constraint and Its Application to Computing the Ground State of Spin-Orbit-Coupled Bose-Einstein Condensates

Calculations of Excited Electronic States by Converging on Saddle Points Using Generalized Mode Following

Direct Unconstrained Optimization of Molecular Orbital Coefficients in Density Functional Theory

On the Convergence of the Self-Consistent Field Iteration in Kohn--Sham Density Functional Theory

Diagonalization without Diagonalization: A Direct Optimization Approach for Solid-State Density Functional Theory

Gradient Optimization of Fermionic Projected Entangled Pair States on Directed Lattices