Direct orbital optimization methods for variational density functional calculations of excited electronic states
ORAL
Abstract
Excited states can be calculated as stationary states of the energy expressed as a density functional. This variational approach provides atomic forces, can describe doubly excited states and, thanks to state-specific orbital relaxation, gives better approximations of charge-transfer, Rydberg and core-level excitations compared to the widely used linear-response time-dependent density functional theory (TDDFT) formalism. However, widespread application in excited-state geometry optimization and molecular dynamics is hindered by limitations of standard self-consistent field (SCF) algorithms. Such algorithms perform well in calculations of non-degenerate ground states (minima of the energy functionals) but commonly fail to converge excited-state solutions because the latter are typically saddle points and are often nearly degenerate to each other or the ground state. To overcome this impediment we develop robust direct orbital optimization methods for variational excited-state DFT calculations [1-3]. It will be shown how a novel approach using quasi-Newton algorithms that can develop negative Hessian eigenvalues in combination with the maximum overlap method can converge excited states in molecules where conventional SCF approaches usually fail, such as charge-transfer excitations in nitrobenzene [1-3]. The new method can be used to calculate potential energy surfaces (PESs) of challenging electronically degenerate systems, as it will be shown for the double bond breaking and conical intesection in ethylene. There, unrestricted broken-symmetry solutions can provide PESs in agreement with multireference calculations. Challenges for application in excited-state structural optimization or dynamics arise due to the presence of lower symmetry solutions with unphysical PESs. It will be illustrated how these issues can be overcome using strategies based on following the lowest eigenmodes of the electronic Hessian to converge on the target excited-state solutions [4].
–
Publication: [1] Ivanov, A. V., Levi, G., Jónsson, E. Ö. & Jónsson, H. Method for Calculating Excited Electronic States Using Density Functionals and Direct Orbital Optimization with Real Space Grid or Plane-Wave Basis Set. J. Chem. Theory Comput. 17, 5034–5049 (2021).<br>[2] Levi, G., Ivanov, A. V. & Jónsson, H. Variational Density Functional Calculations of Excited States via Direct Optimization. J. Chem. Theory Comput. 16, 6968–6982 (2020).<br>[3] Levi, G., Ivanov, A. V. & Jonsson, H. Variational Calculations of Excited States Via Direct Optimization of Orbitals in DFT. Faraday Discuss. 224, 448 (2020).<br>[4] Schmerwitz, L. A. Y., Ivanov, A. V., Jónsson, E. Ö., Jónsson, H. & Levi, G. Orbital optimized density functional theory description of potential energy surfaces in presence of electronic degeneracies: Double bond breaking and conical intersection in ethylene. In preparation (2021).
