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.

Efficient Softest Mode Finding in Transition States Calculations

Constrained Broyden Minimization Combined with the Dimer Method for Locating Transition State of Complex Reactions

Force reversed method for locating transition states

Iterative Minimization Algorithm for Efficient Calculations of Transition States

Locating Transition States by Variational Reaction Path Optimization with an Energy-Derivative-Free Objective Function

Efficient Chain-of-States Approach for Locating Transition State via Spherical Optimization.

Optimization-based Shrinking Dimer Method for Finding Transition States

Adaptive nudged elastic band approach for transition state calculation

A Simple Density Functional Fractional Occupation Number Procedure to Determine the Low Energy Transition Region of Spin-Flip Reactions

TRANSITION PATH THEORY FOR LANGEVIN DYNAMICS ON MANIFOLDS: OPTIMAL CONTROL AND DATA-DRIVEN SOLVER

Gradient Type Optimization Methods for Electronic Structure Calculations

Identification of mechanisms of magnetic transitions using an efficient method for converging on first order saddle points

Accurate transition state generation with an object-aware equivariant elementary reaction diffusion model

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

Force-free identification of minimum-energy pathways and transition states for stochastic electronic structure theories

Crossing The Gap Using Variational Quantum Eigensolver: A Comparative Study

Identifying a Feasible Transition Pathway between Two Conformational States for a Protein

Analytical ab initio hessian from a deep learning potential for transition state optimization

Biomolecular transitions: efficient computation of pathways, free energies, and rates

Data-driven Computation of Molecular Reaction Coordinates

An Efficient Algorithm for Constrained CASSCF(1,2) and CASSCF(3,2) Simulations as Relevant to Electron and Hole Transfer Problems