Search

We present an algorithm for efficient calculation of analytic nonadiabatic derivative couplings between spin-adiabatic, time-dependent density functional theory states within the Tamm-Dancoff approximation. Our derivation is based on the direct differentiation of the Kohn-Sham pseudowavefunction using the framework of Ou et al. Our implementation is limited to the case of a system with an even number of electrons in a closed shell ground state, and we validate our algorithm against finite difference at an S1/T2 crossing of benzaldehyde. Through the introduction of a magnetic field spin-coupling operator, we break time-reversal symmetry to generate complex valued nonadiabatic derivative couplings. Although the nonadiabatic derivative couplings are complex valued, we find that a phase rotation can generate an almost entirely real-valued derivative coupling vector for the case of benzaldehyde. [ABSTRACT FROM AUTHOR]

Working within the Tamm-Dancoff approximation, we calculate the derivative couplings between time-dependent density-functional theory excited states by assuming that the Kohn-Sham superposition of singly excited determinants represents a true electronic wavefunction. All Pulay terms are included in our derivative coupling expression. The reasonability of our approach can be established by noting that, for closely separated electronic states in the infinite basis limit, our final expression agrees exactly with the Chernyak-Mukamel expression (with transition densities from response theory). Finally, we also validate our approach empirically by analyzing the behavior of the derivative couplings around the T1/T2 conical intersection of benzaldehyde. [ABSTRACT FROM AUTHOR]

We demonstrate that Boys-localized diabatic states do indeed exhibit small derivative couplings, as is required of quasidiabatic states. In doing so, we present a general formalism for calculating derivative couplings and analytic gradients for diabatic states. We then develop additional equations specific to the case of Boys-localized configuration-interaction singles (CIS)-in particular, the analytic gradient of the CIS dipole matrix-and we validate our implementation against finite-difference results. In a forthcoming paper, we will publish additional algorithmic and computational details and apply our method to the Closs energy-transfer systems as a further test of the validity of Boys-localized diabatic states. [ABSTRACT FROM AUTHOR]

We present a method for analytically calculating the derivative couplings between a pair of configuration-interaction-singles (CIS) excited states obtained in an atom-centered basis. Our theory is exact and has been derived using two completely independent approaches: one inspired by the Hellmann-Feynman theorem and the other following from direct differentiation. (The former is new, while the latter is in the spirit of existing approaches in the literature.) Our expression for the derivative couplings incorporates all Pulay effects associated with the use of an atom-centered basis, and the computational cost is minimal, roughly comparable to that of a single CIS energy gradient. We have validated our method against CIS finite-difference results and have applied it to the lowest lying excited states of naphthalene; we find that naphthalene derivative couplings include Pulay contributions sufficient to have a qualitative effect. Going beyond standard problems in analytic gradient theory, we have also constructed a correction, based on perturbative electron-translation factors, for including electronic momentum and eliminating spurious components of the derivative couplings that break translational symmetry. This correction is general and can be applied to any level of electronic structure theory. [ABSTRACT FROM AUTHOR]

We study the relative ability of several models of x-ray absorption spectra to capture the Franck–Condon structure apparent from an experiment on gaseous nitrogen. In doing so, we adopt the Born–Oppenheimer approximation and a constrained density functional theory method for computing the energies of the x-ray-excited molecule. Starting from an otherwise classical model for the spectrum, we systematically introduce more realistic physics, first by substituting the quantum mechanical nuclear radial density in the bond separation R for the classical radial density, then by adding the effect of zero-point energy and other level shifts, and finally by including explicit rovibrational quantization of both the ground and excited states. The quantization is determined exactly, using a discrete variable representation (DVR). We show that the near-edge x-ray absorption fine structure (NEXAFS) spectrum can be predicted semiquantitatively within this framework. We also address the possibility of non-trivial temperature dependence in the spectrum. By using constrained density functional theory in combination with more accurate potentials, we demonstrate that it is possible to improve the predicted spectrum. Ultimately, we establish the predictive limits of our method with respect to vibrational fine structure in NEXAFS spectra. [ABSTRACT FROM AUTHOR]

In a recent article, we showed that configuration interaction singles (CIS) has a systematic bias against charge-transfer (CT) states: CT vertical excitation energies are consistently too high (by 1-2 eV) as compared with non-CT energies [J. E. Subotnik, J. Chem. Phys. 137, 071104 (2011)]. We now show that this CIS error can be corrected approximately by performing a single Newton-Raphson step to reoptimize orbitals, thus establishing a new set of orbitals which better balances ground and excited state energies. The computational cost of this correction is exactly that of one coupled-perturbed Hartree-Fock calculation, which is effectively the cost of the CIS calculation itself. In other words, for twice the computational cost of a standard CIS calculation, or roughly the same cost as a linear-response time-dependent Hartree-Fock calculation, one can achieve a balanced, size-consistent description of CT versus non-CT energies, ideally with the accuracy of a much more expensive doubles CIS(D) calculation. [ABSTRACT FROM AUTHOR]

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear-electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an "open teamware" model and an increasingly modular design.

Working within the Tamm-Dancoff approximation, we calculate the derivative couplings between time-dependent density-functional theory excited states by assuming that the Kohn-Sham superposition of singly excited determinants represents a true electronic wavefunction. All Pulay terms are included in our derivative coupling expression. The reasonability of our approach can be established by noting that, for closely separated electronic states in the infinite basis limit, our final expression agrees exactly with the Chernyak-Mukamel expression (with transition densities from response theory). Finally, we also validate our approach empirically by analyzing the behavior of the derivative couplings around the T1/T2 conical intersection of benzaldehyde.