Difference projection-after-variation double-hybrid density functional theory applied to the calculation of vertical excitation energies

Role of Exact Exchange in Difference Projected Double-Hybrid Density Functional Theory for Treatment of Local, Charge Transfer, and Rydberg Excitations

Difference approaches to the study of excited states have undergone a renaissance in recent years, with the development of a plethora of algorithms for locating self-consistent field approximations to excited states. Density functional theory is likely to offer the best balance of cost and accuracy for difference approaches, and yet there has been little investigation of how the parametrization of density functional approximations affects performance. In this work, we aim to explore the role of the global Hartree-Fock exchange parameter in tuning accuracy of different excitation types within the framework of the recently introduced difference projected double-hybrid density functional theory approach and contrast the performance with conventional time-dependent double-hybrid density functional theory. Difference projected double-hybrid density functional theory was demonstrated to give vertical excitation energies with average error and standard deviation with respect to multireference perturbation theory comparable to more expensive linear-response coupled cluster approaches ( J. Chem. Phys.2020, 153, 074103). However, despite benchmarking of local excitations, there has been no investigation of the methods performance for charge transfer or Rydberg excitations. In this work we report a new benchmark of charge transfer, Rydberg, and local excited state vertical excitation energies and examine how the exact Hartree-Fock exchange affects the benchmark performance to provide a deeper understanding of how projection and nonlocal correlation balance differing sources of error in the ground and excited states.

Orbital optimization in nonorthogonal multiconfigurational self-consistent field applied to the study of conical intersections and avoided crossings

Nonorthogonal approaches to electronic structure methods have recently received renewed attention, with the hope that new forms of nonorthogonal wavefunction Ansätze may circumvent the computational bottleneck of orthogonal-based methods. The basis in which nonorthogonal configuration interaction is performed defines the compactness of the wavefunction description and hence the efficiency of the method. Within a molecular orbital approach, nonorthogonal configuration interaction is defined by a "different orbitals for different configurations" picture, with different methods being defined by their choice of determinant basis functions. However, identification of a suitable determinant basis is complicated, in practice, by (i) exponential scaling of the determinant space from which a suitable basis must be extracted, (ii) possible linear dependencies in the determinant basis, and (iii) inconsistent behavior in the determinant basis, such as disappearing or coalescing solutions, as a result of external perturbations, such as geometry change. An approach that avoids the aforementioned issues is to allow for basis determinant optimization starting from an arbitrarily constructed initial determinant set. In this work, we derive the equations required for performing such an optimization, extending previous work by accounting for changes in the orthogonality level (defined as the dimension of the orbital overlap kernel between two determinants) as a result of orbital perturbations. The performance of the resulting wavefunction for studying avoided crossings and conical intersections where strong correlation plays an important role is examined.