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.

An Efficient Step-Forward Way To Solve The Schrodinger Eigenvalue Equation In Self-Consistent Calculations

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

Self-consistent Density Functional Calculations without Diagonalization

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

Adaptive Finite Differencing in High Accuracy Electronic Structure Calculations

Efficient iterative method for solving the Dirac–Kohn–Sham density functional theory

Accelerating Eigenvalue Computation for Nuclear Structure Calculations via Perturbative Corrections

Efficient Spectral and Spectral Element Methods for Eigenvalue Problems of Schrödinger Equations with an Inverse Square Potential

Linear Scaling Calculations of Excitation Energies with Active-Space Particle-Particle Random Phase Approximation

Variational projector-augmented wave method: A full-potential approach for electronic structure calculations in solid-state physics

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

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

Eigenstate calculation in the state-averaged (multi-layer) multi-configurational time-dependent Hartree approach

Robust Determination Of The Chemical Potential In The Pole Expansion And Selected Inversion Method For Solving Kohn-Sham Density Functional Theory

A Simple and Efficient Lattice Summation Method for Metallic Electrodes in Constant Potential Molecular Dynamics Simulation

Parallel implementation of electronic structure eigensolver using a partitioned folded spectrum method

New Algorithms For Iterative Matrix-Free Eigensolvers In Quantum Chemistry

Diagonalization-free self-consistent field approach with localized molecular orbitals

A novel splitting strategy to accelerate solving generalized eigenvalue problem from Kohn--Sham density functional theory

Parallel Self-Consistent-Field Calculations via Chebyshev-Filtered Subspace Acceleration

Gaussian Statistics of the Cosmic Microwave Background: Correlation of Temperature Extrema In the COBE DMR Two-Year Sky Maps