The adiabatic approximation in time-dependent density matrix functional theory: response properties from dynamics of phase-including natural orbitals

The adiabatic approximation is problematic in time-dependent density matrix functional theory. With pure density matrix functionals (invariant under phase change of the natural orbitals) it leads to lack of response in the occupation numbers, hence wrong frequency dependent responses, in particular α(ω→0)≠α(0) (the static polarizability). We propose to relinquish the requirement that the functional must be a pure one-body reduced density matrix (1RDM) functional, and to introduce additional variables which can be interpreted as phases of the one-particle states of the independent particle reference system formed with the natural orbitals, thus obtaining so-called phase-including natural orbital (PINO) functionals. We also stress the importance of the correct choice of the complex conjugation in the two-electron integrals in the commonly used functionals (they should not be of exchange type). We demonstrate with the Lo?wdin-Shull energy expression for two-electron systems, which is an example of a PINO functional, that for two-electron systems exact responses (polarizabilities, excitation energies) are obtained, while writing this energy expression in the usual way as a 1RDM functional yields erroneous responses.     

Adiabatic approximation of time-dependent density matrix functional response theory

Time-dependent density matrix functional theory can be formulated in terms of coupled-perturbed response equations, in which a coupling matrix K(omega) features, analogous to the well-known time-dependent density functional theory (TDDFT) case. An adiabatic approximation is needed to solve these equations, but the adiabatic approximation is much more critical since there is not a good "zero order" as in TDDFT, in which the virtual-occupied Kohn-Sham orbital energy differences serve this purpose. We discuss a simple approximation proposed earlier which uses only results from static calculations, called the static approximation (SA), and show that it is deficient, since it leads to zero response of the natural orbital occupation numbers. This leads to wrong behavior in the omega-->0 limit. An improved adiabatic approximation (AA) is formulated. The two-electron system affords a derivation of exact coupled-perturbed equations for the density matrix response, permitting analytical comparison of the adiabatic approximation with the exact equations. For the two-electron system also, the exact density matrix functional (2-matrix in terms of 1-matrix) is known, enabling testing of the static and adiabatic approximations unobscured by approximations in the functional. The two-electron HeH(+) molecule shows that at the equilibrium distance, SA consistently underestimates the frequency-dependent polarizability alpha(omega), the adiabatic TDDFT overestimates alpha(omega), while AA improves upon SA and, indeed, AA produces the correct alpha(0). For stretched HeH(+), adiabatic density matrix functional theory corrects the too low first excitation energy and overpolarization of adiabatic TDDFT methods and exhibits excellent agreement with high-quality CCSD ("exact") results over a large omega range.     

Excitation energies of pi-conjugated oligomers within time-dependent current-density-functional theory

We study the pi*<--pi singlet excitations of the pi-conjugated oligomers of polyacetylene, polydiacetylene, polybutatriene, polythiophene, poly(para-phenylene vinylene), and the lowest singlet excitations of the hydrogen chain. For this we used time-dependent current-density-functional theory within the Vignale-Kohn and adiabatic local density approximations. By studying the dependence of the excitation spectrum on the chain length we conclude that the reduction of the static polarizability when using the Vignale-Kohn functional has two origins. First, the excitation energies of transitions with a large transition dipole are shifted upward. Second, the character of the transition between the lowest occupied and highest unoccupied molecular orbitals and the oscillator strength of the lowest transition within the adiabatic local density approximation is transferred to higher transitions. The lowest transitions that have a considerable oscillator strength obtained with the Vignale-Kohn functional have excitation energies that are in most cases in better agreement with available reference data than the adiabatic local density approximation.     

Excitation energies from range-separated time-dependent density and density matrix functional theory

Time-dependent density functional theory (TD-DFT) in the adiabatic formulation exhibits known failures when applied to predicting excitation energies. One of them is the lack of the doubly excited configurations. On the other hand, the time-dependent theory based on a one-electron reduced density matrix functional (time-dependent density matrix functional theory, TD-DMFT) has proven accurate in determining single and double excitations of H(2) molecule if the exact functional is employed in the adiabatic approximation. We propose a new approach for computing excited state energies that relies on functionals of electron density and one-electron reduced density matrix, where the latter is applied in the long-range region of electron-electron interactions. A similar approach has been recently successfully employed in predicting ground state potential energy curves of diatomic molecules even in the dissociation limit, where static correlation effects are dominating. In the paper, a time-dependent functional theory based on the range-separation of electronic interaction operator is rigorously formulated. To turn the approach into a practical scheme the adiabatic approximation is proposed for the short- and long-range components of the coupling matrix present in the linear response equations. In the end, the problem of finding excitation energies is turned into an eigenproblem for a symmetric matrix. Assignment of obtained excitations is discussed and it is shown how to identify double excitations from the analysis of approximate transition density matrix elements. The proposed method used with the short-range local density approximation (srLDA) and the long-range Buijse-Baerends density matrix functional (lrBB) is applied to H(2) molecule (at equilibrium geometry and in the dissociation limit) and to Be atom. The method accounts for double excitations in the investigated systems but, unfortunately, the accuracy of some of them is poor. The quality of the other excitations is in general much better than that offered by TD-DFT-LDA or TD-DMFT-BB approximations if the range-separation parameter is properly chosen. The latter remains an open problem.     

Excited state geometries within time-dependent and restricted open-shell density functional theories

Singlet excited state geometries of a set of medium sized molecules with different characteristic lowest excitations are studied. Geometry optimizations of excited states are performed with two closely related restricted open-shell Kohn–Sham methods and within linear response to time-dependent density functional theory. The results are compared to wave-function based methods. Excitation energies (vertical and adiabatic) calculated from the open-shell methods show systematic errors depending on the type of excitation. However, for all states accessible by the restricted methods a good agreement for the geometries with time-dependent density functional theory and wave-function based methods is found. An analysis of the energy with respect to the mixing angle for the singly occupied orbitals reveals that some states (mostly [n→π^*]) are stable when symmetry constraints are relaxed and others (mostly [π→π^*]) are instable. This has major implications on the applicability of the restricted open-shell methods in molecular dynamics simulations.     

Perspectives on double-excitations in TDDFT

The adiabatic approximation in time-dependent density functional theory (TDDFT) yields reliable excitation spectra with great efficiency in many cases, but fundamentally fails for states of double-excitation character. We discuss how double-excitations are at the root of some of the most challenging problems for TDDFT today. We then present new results for (i) the calculation of autoionizing resonances in the helium atom, (ii) understanding the nature of the double excitations appearing in the quadratic response function, and (iii) retrieving double-excitations through a real-time semiclassical approach to correlation in a model quantum dot.     

Relationship between orbital energy gaps and excitation energies for long‐chain systems

The difference between the excitation energies and corresponding orbital energy gaps, the exciton binding energy, is investigated based on time‐dependent (TD) density functional theory (DFT) for long‐chain systems: all‐trans polyacetylenes and linear oligoacenes. The optimized geometries of these systems indicate that bond length alternations significantly depend on long‐range exchange interactions. In TDDFT formalism, the exciton binding energy comes from the two‐electron interactions between occupied and unoccupied orbitals through the Coulomb‐exchange‐correlation integral kernels. TDDFT calculations show that the exciton binding energy is significant when long‐range exchange interactions are involved. Spin‐flip (SF) TDDFT calculations are then carried out to clarify double‐excitation effects in these excitation energies. The calculated SF‐TDDFT results indicate that double‐excitation effects significantly contribute to the excitations of long‐chain systems. The discrepancies between the vertical ionization potential minus electron affinity (IP–EA) values and the HOMO–LUMO excitation energies are also evaluated for the infinitely long polyacetylene and oligoacene using the least‐square fits to estimate the exciton binding energy of infinitely long systems. It is found that long‐range exchange interactions are required to give the exciton binding energy of the infinitely long systems. Consequently, it is concluded that long‐range exchange interactions neglected in many DFT calculations play a crucial role in the exciton binding energies of long‐chain systems, while double‐excitation correlation effects are also significant to hold the energy balance of the excitations. © 2016 Wiley Periodicals, Inc.     

Assessment of dressed time-dependent density-functional theory for the low-lying valence states of 28 organic chromophores

Almost all time-dependent density-functional theory (TDDFT) calculations of excited states make use of the adiabatic approximation, which implies a frequency-independent exchange-correlation kernel that limits applications to one-hole/one-particle states. To remedy this problem, Maitra et al. [N.T. Maitra, F. Zhang, R.J. Cave, K. Burke, Double excitations within time-dependent density functional theory linear response theory, J. Chem. Phys. 120 (2004) 5932 ] proposed dressed TDDFT (D-TDDFT), which includes explicit two-hole/two-particle states by adding a frequency-dependent term to adiabatic TDDFT. This paper offers the first extensive test of D-TDDFT, and its ability to represent excitation energies in a general fashion. We present D-TDDFT excited states for 28 chromophores and compare them with the benchmark results of Schreiber et al. [M. Schreiber, M.R. Silva-Junior, S.P.A. Sauer, W. Thiel, Benchmarks for electronically excited states: CASPT2, CC2, CCSD, and CC3, J. Chem. Phys. 128 (2008) 134110]. We find the choice of functional used for the A-TDDFT step to be critical for positioning the 1h1p states with respect to the 2h2p states. We observe that D-TDDFT without HF exchange increases the error in excitations already underestimated by A-TDDFT. This problem is largely remedied by implementation of D-TDDFT including Hartree-Fock exchange.     

A DFT/TDDFT interpretation of the ground and excited states of porphyrin and porphyrazine complexes

A comprehensive treatment is given of the electronic excitation spectra of Mg, Zn and Ni complexes of porphyrin and porphyrazine using time-dependent density functional theory (TDDFT). It is emphasized that the Kohn–Sham (KS) molecular orbital (MO) method, which is the basis for the TDDFT calculations, affords a MO interpretation of the ground state electronic structure and of the nature of the excitations. This implies that a direct connection can be made to many previous MO treatments of the title compounds. We discuss in particular, how the original explanations of the intensity distribution over the lowest excitations (the Q and B bands) in terms of a cyclic polyene model, or even a free-electron model, can be reconciled with the actual molecular and electronic structure of these compounds being much more complicated than these simple models. A fragment approach is used, building the porphyrin ring from pyrrole rings and CH or N bridges. This leads directly to a simple interpretation of the orbitals of Gouterman's four-orbital model, which are responsible for the Q and B bands. It also leads to additional occupied π-orbitals which are absent in the cyclic polyene model and which need to be invoked to understand other features of the electronic spectra such as the origin of the N, L and M bands. Considerable attention is given to the intensities of the various transitions, explaining why the transitions within the so-called four-orbital model of Gouterman have large transition dipoles, why transitions from additional occupied π-orbitals have relatively small transition dipoles.     

Performance of the M11 and M11-L density functionals for calculations of electronic excitation energies by adiabatic time-dependent density functional theory

Adiabatic time-dependent density functional theory is a powerful method for calculating electronic excitation energies of complex systems, but the quality of the results depends on the choice of approximate density functional. In this article we test two promising new density functionals, M11 and M11-L, against databases of 214 diverse electronic excitation energies, and we compare the results to those for 16 other density functionals of various kinds and to time-dependent Hartree-Fock. Charge transfer excitations are well known to be the hardest challenge for TDDFT. M11 is a long-range-corrected hybrid meta-GGA, and it shows better performance for charge transfer excitations than any of the other functionals except M06-HF, which is a specialized functional that does not do well for valence excitations. Several other long-range-corrected hybrid functionals also do well, and we especially recommend M11, ωB97X, and M06-2X for general spectroscopic applications because they do exceptionally well on ground-state properties as well as excitation energies. Local functionals are preferred for many applications to extended systems because of their significant cost advantage for large systems. M11-L is a dual-range local functional and-unlike all previous local functionals-it has good performance for Rydberg states as well as for valence states. Thus it is highly recommended for excitation energy calculations on extended systems.     

Long-range excitations in time-dependent density functional theory

Adiabatic time-dependent density functional theory fails for excitations of a heteroatomic molecule composed of two open-shell fragments at large separation. Strong frequency dependence of the exchange-correlation kernel is necessary for both local and charge-transfer excitations. The root of this is the static correlation created by the step in the exact Kohn-Sham ground-state potential between the two fragments. An approximate nonempirical kernel is derived for excited molecular dissociation curves at large separation. Our result is also relevant when the usual local and semilocal approximations are used for the ground-state potential, as static correlation there arises from the coalescence of the highest occupied and lowest unoccupied orbital energies as the molecule dissociates.     

Electronic spectrum of UO2(2+) and [UO2Cl4]2- calculated with time-dependent density functional theory

The electronic spectra of UO(2) (2+) and [UO(2)Cl(4)](2-) are calculated with a recently proposed relativistic time-dependent density functional theory method based on the two-component zeroth-order regular approximation for the inclusion of spin-orbit coupling and a noncollinear exchange-correlation functional. All excitations out of the bonding sigma(u) (+) orbital into the nonbonding delta(u) or phi(u) orbitals for UO(2) (2+) and the corresponding excitations for [UO(2)Cl(4)](2-) are considered. Scalar relativistic vertical excitation energies are compared to values from previous calculations with the CASPT2 method. Two-component adiabatic excitation energies, U-O equilibrium distances, and symmetric stretching frequencies are compared to CASPT2 and combined configuration-interaction and spin-orbit coupling results, as well as to experimental data. The composition of the excited states in terms of the spin-orbit free states is analyzed. The results point to a significant effect of the chlorine ligands on the electronic spectrum, thereby confirming the CASPT2 results: The excitation energies are shifted and a different luminescent state is found.     

Computational study of iron(II) systems containing ligands with nitrogen heterocyclic groups

Density functional theory (DFT) calculations show the higher energy HOMO (highest occupied molecular orbital) orbitals of four iron(II) diimine complexes are metal centered and the lower energy LUMO (lowest unoccupied molecular orbitals) are ligand centered. The energy of the orbitals correlates with electrochemical redox potentials of the complexes. Time-dependent density functional theory (TDDFT) calculations reveal ligand centered (LC) and metal-to-ligand charge transfer (MLCT) at higher energy than experimentally observed. TDDFT calculations also reveal the presence of d-d transitions which are buried under the MLCT and LC transitions. The difference in chemical and photophysical behavior of the iron complexes compared to that of their ruthenium analogues is also addressed.     

Theoretical study of the electronic spectroscopy of CO adsorbed on Pt(111)

The excited states of CO adsorbed on the Pt(111) surface are studied using a time-dependent density functional theory formalism. To reduce the computational cost, electronic excitations are computed within a reduced single excitation space. Using cluster models of the surface, excitation energies are computed for CO in the on-top, threefold, and bridge binding sites. On adsorption, there is a lowering of the 5sigma orbital energy. This leads to a large blueshift in the 5sigma- -> pi(CO*) excitation energy for all adsorption sites. The 1pi and 4sigma orbital energies are lowered to a lesser extent, and smaller shifts in the corresponding excitation energies are predicted. For the larger clusters, pi* excitations at lower energies are observed. These transitions correspond to excitations to virtual orbitals of pi* character which lie below the pi* orbitals of gas phase CO. These orbitals are associated predominantly with the metal atoms of the cluster. The excitation energies are also found to be sensitive to changes in the adsorption geometry. The electronic spectrum of CO on Pt(111) is simulated and the assignment of the bands observed in experimental electron energy loss spectroscopy discussed.     

A long-range-corrected time-dependent density functional theory

We apply the long-range correction (LC) scheme for exchange functionals of density functional theory to time-dependent density functional theory (TDDFT) and examine its efficiency in dealing with the serious problems of TDDFT, i.e., the underestimations of Rydberg excitation energies, oscillator strengths, and charge-transfer excitation energies. By calculating vertical excitation energies of typical molecules, it was found that LC-TDDFT gives accurate excitation energies, within an error of 0.5 eV, and reasonable oscillator strengths, while TDDFT employing a pure functional provides 1.5 eV lower excitation energies and two orders of magnitude lower oscillator strengths for the Rydberg excitations. It was also found that LC-TDDFT clearly reproduces the correct asymptotic behavior of the charge-transfer excitation energy of ethylene-tetrafluoroethylene dimer for the long intramolecular distance, unlike a conventional far-nucleus asymptotic correction scheme. It is, therefore, presumed that poor TDDFT results for pure functionals may be due to their lack of a long-range orbital-orbital interaction.     

Origin of substituent effects in the absorption spectra of peroxy radicals: time dependent density functional theory calculations.

We investigate the assignment of electronic transitions in alkyl peroxy radicals. Past experimental work has shown that the phenyl peroxy radical exhibits a transition in the visible region; however, previous high level calculations have not reproduced this observed absorption. We use time dependent density functional theory (TDDFT) to characterize the electronic excitations of the phenyl peroxy radical as well as other hydrocarbon substituted peroxy radicals. TDDFT calculations of the phenyl peroxy radical support an excitation in the visible spectrum. Further, we investigate the nature of this visible absorption using electron attachment/detachment density diagrams of the peroxy radicals and present a qualitative picture of the origin of the visible absorption based on molecular orbital perturbations. The peroxy radical substituent is also compared against isoelectronic radical groups. The visible absorption is determined to be dependent on mixing of the alkyl and radical substituent orbitals.     

Time-dependent density functional theory gradients in the Amsterdam density functional package: geometry optimizations of spin-flip excitations

An implementation of time-dependent density functional theory (TDDFT) energy gradients into the Amsterdam density functional theory program package (ADF) is described. The special challenges presented by Slater-type orbitals in quantum chemical calculation are outlined with particular emphasis on details that are important for TDDFT gradients. Equations for the gradients of spin-flip TDDFT excitation energies are derived. Example calculations utilizing the new implementation are presented. The results of standard calculations agree well with previous results. It is shown that starting from a triplet reference, spin-flip TDDFT can successfully optimize the geometry of the four lowest singlet states of CH₂ and three other isovalent species. Spin-flip TDDFT is used to calculate the potential energy curve of the breaking of the C?CC bond of ethane. The curve obtained is superior to that from a restricted density functional theory calculation, while at the same time the problems with spin contamination exhibited by unrestricted density functional theory calculations are avoided.     

A TDDFT study of the optical response of DNA bases, base pairs, and their tautomers in the gas phase

We present calculations of the optical response of the DNA bases and base pairs both in their normal and tautomeric forms in the gas phase, using time-dependent density functional theory (TDDFT). These calculations are performed in real time within the adiabatic approximation with a basis of local orbitals. Our results for the individual bases are in good agreement with experiment and computationally more demanding calculations of chemical accuracy. The optical response of base pairs indicates that the differences between normal and tautomeric forms in certain cases are significant enough to provide a means of identification.     

Asymptotic correction of the exchange-correlation kernel of time-dependent density functional theory for long-range charge-transfer excitations

Time-dependent density functional theory (TDDFT) calculations of charge-transfer excitation energies omegaCT are significantly in error when the adiabatic local density approximation (ALDA) is employed for the exchange-correlation kernel fxc. We relate the error to the physical meaning of the orbital energy of the Kohn-Sham lowest unoccupied molecular orbital (LUMO). The LUMO orbital energy in Kohn-Sham DFT--in contrast to the Hartree-Fock model--approximates an excited electron, which is correct for excitations in compact molecules. In CT transitions the energy of the LUMO of the acceptor molecule should instead describe an added electron, i.e., approximate the electron affinity. To obtain a contribution that compensates for the difference, a specific divergence of fxc is required in rigorous TDDFT, and a suitable asymptotically correct form of the kernel fxc(asymp) is proposed. The importance of the asymptotic correction of fxc is demonstrated with the calculation of omegaCT(R) for the prototype diatomic system HeBe at various separations R(He-Be). The TDDFT-ALDA curve omegaCT(R) roughly resembles the benchmark ab initio curve omegaCT CISD(R) of a configuration interaction calculation with single and double excitations in the region R=1-1.5 A, where a sizable He-Be interaction exists, but exhibits the wrong behavior omegaCT(R)相似文献