The performance and relationship among range-separated schemes for density functional theory

The performance and relationship among different range-separated (RS) hybrid functional schemes are examined using the Coulomb-attenuating method (CAM) with different values for the fractions of exact Hartree-Fock (HF) exchange (α), long-range HF (β), and a range-separation parameter (μ), where the cases of α + β = 1 and α + β = 0 were designated as CA and CA0, respectively. Attenuated PBE exchange-correlation functionals with α = 0.20 and μ = 0.20 (CA-PBE) and α = 0.25 and μ = 0.11 (CA0-PBE) are closely related to the LRC-ωPBEh and HSE functionals, respectively. Time-dependent density functional theory calculations were carried out for a number of classes of molecules with varying degrees of charge-transfer (CT) character to provide an assessment of the accuracy of excitation energies from the CA functionals and a number of other functionals with different exchange hole models. Functionals that provided reasonable estimates for local and short-range CT transitions were found to give large errors for long-range CT excitations. In contrast, functionals that afforded accurate long-range CT excitation energies significantly overestimated energies for short-range CT and local transitions. The effects of exchange hole models and parameters developed for RS functionals for CT excitations were analyzed in detail. The comparative analysis across compound classes provides a useful benchmark for CT excitations.     

Structural and electronic properties of polyacetylene and polyyne from hybrid and Coulomb-attenuated density functionals

The bond length alternation (BLA), the highest-occupied-lowest-unoccupied (HO-LU) orbital energy gap, and the corresponding excitation energy are determined for trans-polyacetylene (PA) and polyyne (PY) using density functional theory. Results from the Coulomb-attenuated CAM-B3LYP functional are compared with those from the conventional BHHLYP and B3LYP hybrid functionals. BLA values and HO-LU gaps are determined using both finite oligomer and infinite chain calculations, subject to periodic boundary conditions. TDDFT excitation energies are determined for the oligomers. The oligomer excitation energies and HO-LU gaps are then used, in conjunction with the infinite chain HO-LU gap, to estimate the infinite chain excitation energy. Overall, BHHLYP and CAM-B3LYP give BLA values and excitation energies that are larger and more accurate than those obtained using B3LYP. The results highlight the degree to which excitation energies can be approximated using the HO-LU gaps-at the infinite limit, this approximation works well for B3LYP, but not for the other functionals, where the HO-LU gap is significantly larger. The study provides further evidence for the high-quality theoretical predictions that can be obtained from the CAM-B3LYP functional.     

Charge-transfer excited states in a pi-stacked adenine dimer, as predicted using long-range-corrected time-dependent density functional theory

The lowest few electronic excitations of a pi-stacked adenine dimer in its B-DNA geometry are investigated, in the gas phase and in a water cluster, using a long-range-corrected version of time-dependent density functional theory (TD-DFT) that asymptotically incorporates Hartree-Fock exchange. Long-range correction is shown to eliminate the catastrophic underestimation of charge-transfer (CT) excitation energies that plagues conventional TD-DFT, at the expense of introducing one adjustable parameter, mu, that determines the length scale on which Hartree-Fock exchange is turned on. This parameter allows us to interpolate smoothly between hybrid density functionals and time-dependent Hartree-Fock theory. Excitation energies for CT states (in which an electron is transferred from one adenine molecule to the other) are found to increase dramatically as a function of mu. Uncorrected hybrid functionals underestimate the CT excitation energies, placing them well below the valence excitations, while time-dependent Hartree-Fock calculations place these states well above the valence states. Values for mu determined from certain benchmark calculations place the CT states well above the valence pipi* and npi* states at the Franck-Condon point.     

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.     

Charge-transfer excitations in uranyl tetrachloride ([UO2Cl4]2-): how reliable are electronic spectra from relativistic time-dependent density functional theory?

Four-component relativistic time-dependent density functional theory (TD-DFT) is used to study charge-transfer (CT) excitation energies of the uranyl molecule as well as the uranyl tetrachloride complex. Adiabatic excitation energies and vibrational frequencies of the excited states are calculated for the lower energy range of the spectrum. The results for TD-DFT with the CAM-B3LYP exchange-correlation functional for the [UO(2)Cl(4)](2-) system are in good agreement with the experimentally observed spectrum of this species and agree also rather well with other theoretical data. Use of the global hybrid B3LYP gives qualitatively correct results, while use of the BLYP functional yields results that are qualitatively wrong due to the too low CT states calculated with this functional. The applicability of the overlap diagnostic of Peach et al. (J. Chem. Phys.2008, 128, 044118) to identify such CT excitations is investigated for a wide range of vertical transitions using results obtained with three different approximate exchange-correlation functionals: BLYP, B3LYP, and CAM-B3LYP.     

Influence of the Delocalization Error and Applicability of Optimal Functional Tuning in Density Functional Calculations of Nonlinear Optical Properties of Organic Donor–Acceptor Chromophores

Nonempirically tuned hybrid density functionals with range‐separated exchange are applied to calculations of the first hyperpolarizability (β_∥) and charge‐transfer (CT) excitations of linear “push–pull” donor–acceptor‐substituted organic molecules with extended π‐conjugated bridges. An unphysical delocalization with increasing chain length in density functional calculations can be reduced significantly by enforcing an asymptotically correct exchange‐correlation potential adjusted to give frontier orbital energies representing ionization potentials. The delocalization error for a number of donor–acceptor systems is quantified by calculations with fractional electron numbers and from orbital localizations. Optimally tuned hybrid variants of the PBE functional incorporating range‐separated exchange can produce similar magnitudes for β_∥ as Møller–Plesset second‐order perturbation (MP2) correlated calculations. Improvements are also found for CT excitation energies, with results similar to an approximate coupled‐cluster model (CC2).     

Assessment of the accuracy of density functionals for prediction of relative energies and geometries of low-lying isomers of water hexamers

Water hexamers provide a critical testing ground for validating potential energy surface predictions because they contain structural motifs not present in smaller clusters. We tested the ability of 11 density functionals (four of which are local and seven of which are nonlocal) to accurately predict the relative energies of a series of low-lying water hexamers, relative to the CCSD(T)/aug'-cc-pVTZ level of theory, where CCSD(T) denotes coupled cluster theory with an interative treatment of single and double excitations and a quasi-perturbative treatment of connected triple excitations. Five of the density functionals were tested with two different basis sets, making a total of 16 levels of density functional theory (DFT) tested. When single-point energy calculations are carried out on geometries obtained with second-order M?ller-Plesset perturbation theory (MP2), only three density functionals, M06-L, M05-2X, and M06-2X, are able to correctly predict the relative energy ordering of the hexamers. These three functionals predict that the range of energies spanned by the six isomers is 3.2-5.6 kcal/mol, whereas the other eight functionals predict ranges of 1.0-2.4 kcal/mol; the benchmark value for this range is 3.1 kcal/mol. When the hexamers are optimized at each level of theory, all methods are able to reproduce the MP2 geometries well for all isomers except the boat and bag isomers, and DFT optimization changes the energy ordering for seven of the 16 methods tested. The addition of zero-point energy changes the energy ordering for all of the density functionals studied except for M05-2X and M06-2X. The variation in relative energies predicted by the different methods highlights the necessity for exercising caution in the choice of density functionals used in future studies. Of the 11 density functionals tested, the most accurate results for energies were obtained with the PWB6K, MPWB1K, and M05-2X functionals.     

Away from generalized gradient approximation: orbital-dependent exchange-correlation functionals

The local-density approximation of density functional theory (DFT) is remarkably accurate, for instance, for geometries and frequencies, and the generalized gradient approximations have also made bond energies quite reliable. Sometimes, however, one meets with failure in individual cases. One of the possible routes towards better functionals would be the incorporation of orbital dependence (which is an implicit density dependency) in the functionals. We discuss this approach both for energies and for response properties. One possibility is the use of the Hartree-Fock-type exchange energy expression as orbital-dependent functional. We will argue that in spite of the increasing popularity of this approach, it does not offer any advantage over Hartree-Fock for energies. We will advocate not to apply the separation of exchange and correlation, which is so ingrained in quantum chemistry, but to model both simultaneously. For response properties the energies and shapes of the virtual orbitals are crucial. We will discuss the benefits that Kohn-Sham potentials can offer which are derived from either an orbital-dependent energy functional, including the exact-exchange functional, or which can be obtained directly as orbital-dependent functional. We highlight the similarity of the Hartree-Fock and Kohn-Sham occupied orbitals and orbital energies, and the essentially different meanings the virtual orbitals and orbital energies have in these two models. We will show that these differences are beneficial for DFT in the case of localized excitations (in a small molecule or in a fragment), but are detrimental for charge-transfer excitations. Again, orbital dependency, in this case in the exchange-correlation kernel, offers a solution.     

Performance assessment of density-functional methods for study of charge-transfer complexes

Various density functionals are applied to a number of weakly bound intermolecular pi-pi charge-transfer (CT) complexes. Most functionals, including the recently developed mPWPW91 and mPW1PW91, grossly underestimate experimental excitation energies; good agreement is obtained only with the half-and-half hybrid BH&HLYP functional. PW91PW91 provides the best agreement with intermolecular distances measured in crystal, while the BH&HLYP values are about 0.1 A too long. Various hybrid functionals with nonlocal exchange correction provide binding energies that compare favorably with the experimental heats of formation measured in solution.     

Description of core excitations by time-dependent density functional theory with local density approximation, generalized gradient approximation, meta-generalized gradient approximation, and hybrid functionals

Time-dependent density functional theory (TDDFT) is employed to investigate exchange-correlation-functional dependence of the vertical core-excitation energies of several molecules including H, C, N, O, and F atoms. For the local density approximation (LDA), generalized gradient approximation (GGA), and meta-GGA, the calculated X1s-->pi* excitation energies (X = C, N, O, and F) are severely underestimated by more than 13 eV. On the other hand, time-dependent Hartree-Fock (TDHF) overestimates the excitation energies by more than 6 eV. The hybrid functionals perform better than pure TDDFT because HF exchange remedies the underestimation of pure TDDFT. Among these hybrid functionals, the Becke-Half-and-Half-Lee-Yang-Parr (BHHLYP) functional including 50% HF exchange provides the smallest error for core excitations. We have also discovered the systematic trend that the deviations of TDHF and TDDFT with the LDA, GGA, and meta-GGA functionals show a strong atom-dependence. Namely, their deviations become larger for heavier atoms, while the hybrid functionals are significantly less atom-dependent.     

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.     

Benchmarking two-photon absorption with CC3 quadratic response theory, and comparison with density-functional response theory

We present a detailed study of the effects of electron correlation on two-photon absorption calculated by coupled cluster quadratic response theory. The hierarchy of coupled cluster models CCS, CC2, CCSD, and CC3 has been used to investigate the effects of electron correlation on the two-photon absorption cross sections of formaldehyde (CH2O), diacetylene (C4H2), and water (H2O). In particular, the effects of triple excitations on two-photon transition cross sections are determined for the first time. In addition, we present a detailed comparison of the coupled cluster results with those obtained from Hartree-Fock and density-functional response theories. We have investigated the local-density approximation, the pure Becke-Lee-Yang-Parr (BLYP) functional, the hybrid Becke-3-parameter-Lee-Yang-Parr (B3LYP), and the Coulomb-attenuated B3LYP (CAM-B3LYP) functionals. Our results show that the CAM-B3LYP functional, when used in conjuction with a one-particle basis-set containing diffuse functions, has much promise; however, care must still be exercised for diffuse Rydberg-type states.     

Optical excitations of defects in realistic nanoscale silica clusters: comparing the performance of density functional theory using hybrid functionals with correlated wavefunction methods

Optical excitations of low energy silica (SiO(2))(4) clusters obtained by global optimization, as opposed to constructed by hand, are studied using a range of theoretical methods. By focusing on the lowest energy silica clusters we hope to capture at least some of the characteristic ways by which the dry surfaces of silica nanosystems preferentially terminate. Employing the six lowest energy (SiO(2))(4) cluster isomers, we show that they exhibit a surprisingly wide range of geometries, defects, and associated optical excitations. Some of the clusters show excitations localized on isolated defects, which are known from previous studies using hydrogen-terminated versions of the defect in question. Other clusters, however, exhibit novel charge-transfer excitations in which an electron transfers between two spatially separated defects. In these cases, because of the inherent proximity of the constituent defects due to the small cluster dimensions, the excitation spectrum is found to be very different from that of the same defects in isolation. Excitation spectra of all clusters were calculated using time-dependent density functional theory (TD-DFT) and delta-SCF DFT (DeltaDFT) methods employing two different hybrid density functionals (B3LYP and BB1K) differing essentially in the amount of incorporated Hartree-Fock-like exchange (HFLE). In all cases the results were compared with CASPT2 calculated values which are taken as a benchmark standard. In line with previous work, the spatially localized excitations are found to be well described by TD-DFT/B3LYP but which gives excitation energies that are significantly underestimated in the case of the charge-transfer excitations. The TD-DFT/BB1K combination in contrast is found to give generally good excitation energies for the lowest excited states of both localized and charge-transfer excitations. Finally, our calculations suggest that the increased quality of the predicted excitation spectra by adding larger amounts of HFLE is mainly due to an increased localization of the excited state associated with the elimination of spurious self-interaction inherent to (semi-)local DFT functionals.     

Accurate evaluation of valence and low-lying Rydberg states with standard time-dependent density functional theory

Using a standard exchange-correlation functional, namely, PBE0, the basis set dependence of time-dependent density functional theory (TD-DFT) calculations has been explored using 33 different bases and five organic molecules as test cases. The results obtained show that this functional can provide accurate (i.e., at convergence) results for both valence and low-lying Rydberg excitations if at least one diffuse function for the heavy atoms is included in the basis set. Furthermore, these results are in fairly good agreement with the experimental data and with those delivered by other functionals specifically designed to yield correct asymptotic/long-range behavior. More generally, the PBE0 calculations show that a greater accuracy can be obtained for both Rydberg and valence excitations if they occur at energies below the epsilonHOMO + 1 eV threshold. This latter value is proposed as a thumb rule to verify the accuracy of TD-DFT/PBE0 applications.     

Functional dependence of core-excitation energies

We examine in depth the functional dependence of computed core-electron binding and excitation energies based on a total-energy difference approach within Kohn-Sham density functional theory. Twenty-seven functional combinations were studied using a database of reliable experimental data on 18 molecules. The computed core-electron binding energies are largely dependent on the choice of exchange functional. The term value of the first resonant excited state and energy differences between the lowest core-excited states are, however, quite insensitive to the choice of functionals since the errors due to the core-region cancel out. Using these results we define a different exchange functional, which mixes two functionals designed by Perdew and Wang (PD86 and PD91), with the best results for both excitation and binding energies obtained for a mixing ratio 60:40 between these. We also reexamine the relativistic corrections for inner-shell excitations.