Self-Consistent Constricted Variational Theory RSCF-CV(∞)-DFT and Its Restrictions To Obtain a Numerically Stable ΔSCF-DFT-like Method: Theory and Calculations for Triplet States

Ensemble Density Functional Theory of Neutral and Charged Excitations : Exact Formulations, Standard Approximations, and Open Questions

Recent progress in the field of (time-independent) ensemble density-functional theory (DFT) for excited states are reviewed. Both Gross-Oliveira-Kohn (GOK) and N-centered ensemble formalisms, which are mathematically very similar and allow for an in-principle-exact description of neutral and charged electronic excitations, respectively, are discussed. Key exact results, for example, the equivalence between the infamous derivative discontinuity problem and the description of weight dependencies in the ensemble exchange-correlation density functional, are highlighted. The variational evaluation of orbital-dependent ensemble Hartree-exchange (Hx) energies is discussed in detail. We show in passing that state-averaging individual exact Hx energies can lead to severe (although solvable) v-representability issues. Finally, we explore the possibility of using the concept of density-driven correlation, which has been introduced recently and does not exist in regular ground-state DFT, for improving state-of-the-art correlation density-functional approximations for ensembles. The present review reflects the efforts of a growing community to turn ensemble DFT into a rigorous and reliable low-cost computational method for excited states. We hope that, in the near future, this contribution will stimulate new formal and practical developments in the field.

Method for Calculating Excited Electronic States Using Density Functionals and Direct Orbital Optimization with Real Space Grid or Plane-Wave Basis Set

A direct orbital optimization method is presented for density functional calculations of excited electronic states using either a real space grid or a plane-wave basis set. The method is variational, provides atomic forces in the excited states, and can be applied to Kohn-Sham (KS) functionals as well as orbital-density-dependent (ODD) functionals including explicit self-interaction correction. The implementation for KS functionals involves two nested loops: (1) An inner loop for finding a stationary point in a subspace spanned by the occupied and a few virtual orbitals corresponding to the excited state; (2) an outer loop for minimizing the energy in a tangential direction in the space of the orbitals. For ODD functionals, a third loop is used to find the unitary transformation that minimizes the energy functional among occupied orbitals only. Combined with the maximum overlap method, the algorithm converges in challenging cases where conventional self-consistent field algorithms tend to fail. The benchmark tests presented include two charge-transfer excitations in nitrobenzene and an excitation of CO to degenerate π* orbitals where the importance of complex orbitals is illustrated. The application of this method to several metal-to-ligand charge-transfer and metal-centered excited states of an Fe^II photosensitizer complex is described, and the results are compared to reported experimental estimates. This method is also used to study the effect of the Perdew-Zunger self-interaction correction on valence and Rydberg excited states of several molecules, both singlet and triplet states, and the performance compared to semilocal and hybrid functionals.

Variational calculations of excited states via direct optimization of the orbitals in DFT

A direct optimization method for obtaining excited electronic states using density functionals is presented.

Half-Projected σ Self-Consistent Field For Electronic Excited States

Fully self-consistent mean-field solutions of electronic excited states have been much less accessible compared to ground state solutions (e.g., Hartree-Fock). The main reason for this is that most excited states are energy saddle points, and hence energy-based optimization methods such as Δ-SCF often collapse to the ground state. Recently, our research group has developed a new method, σ-SCF [ J. Chem. Phys. 2017 , 147 , 214104 ], that successfully solves the "variational collapse" problem of energy-based methods. Despite the success, σ-SCF solutions are often spin-contaminated for open-shell states due to the single-determinant nature; unphysical behaviors such as disappearing solutions and discontinuous first-order energy derivatives are also observed along with the spontaneous breaking of spin or spatial symmetries. In this work, we tackle these problems by partially restoring the broken spin-symmetry of a σ-SCF solution through an approximate spin-projection scheme called half-projection. Orbitals of the projected wave function are optimized in a variation-after-projection (VAP) manner. The resulting theory, which we term half-projected (HP) σ-SCF, brings substantial improvement to the description of singlet and triplet excitations of the original σ-SCF method. Numerical simulations on small molecules suggest that HP σ-SCF delivers high-quality excited-state solutions that exist in a wide range of geometries with smooth potential energy surfaces. We also show that local excitations in HP σ-SCF can be size-intensive.

Constricted variational density functional theory for spatially clearly separated charge-transfer excitations

Constricted Variational Density Functional Theory (CV-DFT) is known to be one of the successful methods in predicting charge-transfer excitation energies. In this paper, we apply the CV-DFT method to the well-known model systems ethylene-tetrafluoroethylene (C₂H₄ × C₂F₄) and the zincbacteriochlorin-bacteriochlorin complex (ZnBC-BC). The analysis of the CV-DFT energies enables us to understand the -1/R charge-transfer behaviour in CV-DFT for large separation distances R. With this we discuss the importance of orbital relaxations using the relaxed version of CV(∞)-DFT, the R-CV(∞)-DFT method. Possible effects of the optimization of the transition matrix for the relaxed self-consistent field version of CV(∞)-DFT, RSCF-CV(∞)-DFT in the case of large fragment separations are shown and we introduce two possible gradient restrictions to avoid the unwanted admixing of other transitions.