Time-dependent Density Functional Theory Excitation Research Articles

Electronic structure and optical spectral analysis of the MnO42- anion with consideration of site and Jahn-Teller distortion.

The ground state of 3d1 MnO42- was studied by density functional theory (DFT) and complete active space self-consistent field (CASSCF) methods in terms of a variety of molecular point group structures to ascertain the site and Jahn-Teller (JT) distortion effect. Modeling results from UB3LYP/6-31+G(d) calculations with natural bond orbital analysis show the four Mn-O bonds are coordinate covalent. The one-electron matrix elements from CASSCF(AILFT) (abinitio ligand field theory) with a second order perturbation treatment were used to calculate the parameters of the angular overlap model. These allowed the calculation of JT stabilization energies and first and second order JT coupling constants for MnO42- with C2v and D2d symmetry. Absorption spectra and excited state transition energies were calculated assuming the lattice distorted geometry was C2v with time-dependent DFT (TDDFT) using the unrestricted coulomb-attenuating UCAM-B3LYP density functional and with the spin-adapted spin-flip DFT using the UBH&HLYP density functional, both with the AUG-CC-PVTZ basis set. The mean absolute deviation of the calculations from experimental excited states for the ligand field and ligand-to-metal charge transfer (LMCT) bands was better than 0.1eV. Several new assignments for LMCT excited states were made on the basis of the TDDFT excitation results.

Revisiting the Performance of Time-Dependent Density Functional Theory for Electronic Excitations: Assessment of 43 Popular and Recently Developed Functionals from Rungs One to Four.

In this paper, the performance of more than 40 popular or recently developed density functionals is assessed for the calculation of 463 vertical excitation energies against the large and accurate QuestDB benchmark set. For this purpose, the Tamm-Dancoff approximation offers a good balance between computational efficiency and accuracy. The functionals ωB97X-D and BMK are found to offer the best performance overall with a root-mean square error (RMSE) of around 0.27 eV, better than the computationally more demanding CIS(D) wave function method with a RMSE of 0.36 eV. The results also suggest that Jacob's ladder still holds for time-dependent density functional theory excitation energies, though hybrid meta generalized-gradient approximations (meta-GGAs) are not generally better than hybrid GGAs. Effects of basis set convergence, gauge invariance correction to meta-GGAs, and nonlocal correlation (VV10) are also studied, and practical basis set recommendations are provided.

The TDDFT Excitation Energies of the BODIPYs; The DFT and TDDFT Challenge Continues.

The derivatives of 4,4-difluoro-4-bora-3a,4a-diaza-s-indacene (BODIPY) are pivotal ingredients for a large number of functional, stimuli-responsive materials and therapeutic molecules based on their photophysical properties, and there is a urgent need to understand and predict their optical traits prior to investing a large amount of resources in preparing them. Density functional theory (DFT) and time-dependent DFT (TDDFT) computations were performed to calculate the excitation energies of the lowest-energy singlet excited state of a large series of common BODIPY derivatives employing various functional aiming at the best possible combination providing the least deviations from the experimental values. Using the common “fudge” correction, a series of combinations was investigated, and a methodology is proposed offering equal or better performances than what is reported in the literature.

Self-interaction correction, electrostatic, and structural influences on time-dependent density functional theory excitations of bacteriochlorophylls from the light-harvesting complex 2.

First-principles calculations offer the chance to obtain a microscopic understanding of light-harvesting processes. Time-dependent density functional theory can have the computational efficiency to allow for such calculations. However, the (semi-)local exchange-correlation approximations that are computationally most efficient fail to describe charge-transfer excitations reliably. We here investigate whether the inexpensive average density self-interaction correction (ADSIC) remedies the problem. For the systems that we study, ADSIC is even more prone to the charge-transfer problem than the local density approximation. We further explore the recently reported finding that the electrostatic potential associated with the chromophores' protein environment in the light-harvesting complex 2 beneficially shifts spurious excitations. We find a great sensitivity on the chromophores' atomistic structure in this problem. Geometries obtained from classical molecular dynamics are more strongly affected by the spurious charge-transfer problem than the ones obtained from crystallography or density functional theory. For crystal structure geometries and density-functional theory optimized ones, our calculations confirm that the electrostatic potential shifts the spurious excitations out of the energetic range that is most relevant for electronic coupling.

Relationships between Orbital Energies, Optical and Fundamental Gaps, and Exciton Shifts in Approximate Density Functional Theory and Quasiparticle Theory.

The relationships between Kohn-Sham (KS) and generalized KS (GKS) density functional orbital energies and fundamental gaps or optical gaps raise many interesting questions including the physical meanings of KS and GKS orbital energies when computed with currently available approximate density functionals (ADFs). In this work, by examining three diverse databases with various ADFs, we examine such relations from the point of view of the exciton shift of quasiparticle theory. We start by calculating a large number of excitation energies by time-dependent density functional theory (TDDFT) with a large number of ADFs. To relate the exciton shift implicit in TDDFT to the exciton shift that is explicit in Green's function theory, we define the exciton shift in TDDFT as the difference of the response shift and the quasiparticle shift. We found a strong correlation between the response shift and the amount of Hartree-Fock exchange included in the density functional, with the response shift varying between -1 and 5 eV. This range is an order of magnitude larger than the mean errors of the TDDFT excitation energies. This result suggests that, with currently available functionals, the KS or GKS orbital energies should be treated as intermediate mathematical variables in the calculation of excitation energies rather than as the energies of independent-particle reference states for quasiparticle theory.

Benchmarking Quantum Chemistry Methods in Calculations of Electronic Excitations

Quantum chemistry methods are applied to obtain numerical solutions of the Schr¨odinger equation for molecular systems. Calculations of transitions between electronic states of large molecules present one of the greatest challenges in this field which require the use of supercomputer resources. In this work we describe the results of benchmark calculations of electronic excitation in the protein domains which were designed to engineer novel fluorescent markers operating in the near-infrared region. We demonstrate that such complex systems can be efficiently modeled with the hybrid qunatum mechanics/molecular mechanics approach (QM/MM) using the modern supercomputers. More specifically, the time-dependent density functional theory (TD-DFT) method was primarily tested with respect to its performance and accuracy. GAMESS (US) and NWChem software were benchmarked in direct and storage-based TDDFT calculations with the hybrid B3LYP density functional, both showing good scaling up to 32 nodes. We note that conventional SCF calculations greatly outperform direct SCF calculations for our test system. Accuracy of TD-DFT excitation energies was estimated by a comparison to the more accurate ab initio XMCQDPT2 method.

Valence and charge-transfer optical properties for some SinCm (m, n ≤ 12) clusters: Comparing TD-DFT, complete-basis-limit EOMCC, and benchmarks from spectroscopy.

Accurate optical characterization of the closo-Si12C12 molecule is important to guide experimental efforts toward the synthesis of nano-wires, cyclic nano-arrays, and related array structures, which are anticipated to be robust and efficient exciton materials for opto-electronic devices. Working toward calibrated methods for the description of closo-Si12C12 oligomers, various electronic structure approaches are evaluated for their ability to reproduce measured optical transitions of the SiC2, Si2Cn (n = 1-3), and Si3Cn (n = 1, 2) clusters reported earlier by Steglich and Maier [Astrophys. J. 801, 119 (2015)]. Complete-basis-limit equation-of-motion coupled-cluster (EOMCC) results are presented and a comparison is made between perturbative and renormalized non-iterative triples corrections. The effect of adding a renormalized correction for quadruples is also tested. Benchmark test sets derived from both measurement and high-level EOMCC calculations are then used to evaluate the performance of a variety of density functionals within the time-dependent density functional theory (TD-DFT) framework. The best-performing functionals are subsequently applied to predict valence TD-DFT excitation energies for the lowest-energy isomers of SinC and Sin-1C7-n (n = 4-6). TD-DFT approaches are then applied to the SinCn (n = 4-12) clusters and unique spectroscopic signatures of closo-Si12C12 are discussed. Finally, various long-range corrected density functionals, including those from the CAM-QTP family, are applied to a charge-transfer excitation in a cyclic (Si4C4)4 oligomer. Approaches for gauging the extent of charge-transfer character are also tested and EOMCC results are used to benchmark functionals and make recommendations.

Do fractionally incremented nuclear charges improve time-dependent density functional theory excitation energies as reliably as fractional orbital populations?

Gaiduk et al. (Phys Rev Lett 108:253005, 2012) showed that one can improve local, semilocal, and hybrid approximations to the Kohn–Sham effective potentials of atoms and molecules by removing a system-independent fraction of electron charge from the highest occupied molecular orbital (HOMO); if the corrected Kohn–Sham potential is used for adiabatic linear-response time-dependent density functional theory (TDDFT) calculations, accurate Rydberg excitation energies are obtained. One may ask whether the same effect could also be achieved by fractionally increasing the positive charges of the nuclei. We investigate this question and find that a small increase in nuclear charges can indeed substantially reduce errors in TDDFT Rydberg excitation energies. However, the optimal magnitude of the charge increase is system-dependent. In addition, the procedure is ambiguous for molecules, where one has to decide how to distribute the additional charge among individual nuclei. These two drawbacks of the fractional nuclear charge method make it disadvantageous compared to the HOMO depopulation technique.

On low-lying excited states of extended nanographenes.

Low-lying excited states of planarly extended nanographenes are investigated using the long-range corrected (LC) density functional theory (DFT) and the spin-flip (SF) time-dependent density functional theory (TDDFT) by exploring the long-range exchange and double-excitation correlation effects on the excitation energies, band gaps, and exciton binding energies. Optimizing the geometries of the nanographenes indicates that the long-range exchange interaction significantly improves the CC bond lengths and amplify their bond length alternations with overall shortening the bond lengths. The calculated TDDFT excitation energies show that long-range exchange interaction is crucial to provide accurate excitation energies of small nanographenes and dominate the exciton binding energies in the excited states of nanographenes. It is, however, also found that the present long-range correction may cause the overestimation of the excitation energy for the infinitely wide graphene due to the discrepancy between the calculated band gaps and vertical ionization potential (IP) minus electron affinity (EA) values. Contrasting to the long-range exchange effects, the SF-TDDFT calculations show that the double-excitation correlation effects are negligible in the low-lying excitations of nanographenes, although this effect is large in the lowest excitation of benzene molecule. It is, therefore, concluded that long-range exchange interactions should be incorporated in TDDFT calculations to quantitatively investigate the excited states of graphenes, although TDDFT using a present LC functional may provide a considerable excitation energy for the infinitely wide graphene mainly due to the discrepancy between the calculated band gaps and IP-EA values. © 2017 Wiley Periodicals, Inc.

On the calculation of [formula omitted] for electronic excitations in time-dependent density-functional theory

Excited states are often treated within the context of time-dependent (TD) density-functional theory (DFT), making it important to be able to assign the excited spin-state symmetry. While there is universal agreement on how Δ〈Sˆ2〉, the difference between 〈Sˆ2〉 for ground and excited states, should be calculated in a wave-function-like formalism such as the Tamm–Dancoff approximation (TDA), confusion persists as to how to determine the spin-state symmetry of excited states in TD-DFT. We try to clarify the origins of this confusion by examining various possibilities for the parameters (σ1,σ2) in the formula Δ〈Sˆ2〉=[Δ〈SˆTDA2〉(X→)+Δ〈SˆTDA2〉(Y→∗)+σ1Δ〈Sˆmixed2〉(X→,Y→∗)]/(X→†X→+σ2Y→†Y→), where X→ is the particle–hole part and Y→ is the hole–particle part of the response theory vector. A first principles derivation leads directly to (σ1,σ2)=(+1,−1) which we argue is the best simple formula linking spin with energy, albeit approximately. On the other hand, if the desire is to recover wave-function-like values of Δ〈Sˆ2〉, then we argue that the choice (σ1,σ2)=(+1,+1) should be made. Additional examples are offered to justify that the choice of σ1=0 should also be made when seeking wave-function-like values of Δ〈Sˆ2〉.

Generalized Energy-Based Fragmentation Approach for Localized Excited States of Large Systems.

We have extended the generalized energy-based fragmentation (GEBF) approach to localized excited states of large systems. In this approach, the excited-state energy of a large system could be expressed as the combination of the excited-state energies of "active subsystems", which contains the chromophore center, and the ground-state energies of "inactive subsystems". The GEBF approach has been implemented at the levels of time-dependent density functional theory (TDDFT) and approximate coupled cluster singles and doubles (CC2) method. Our results show that GEBF-TDDFT can reproduce the TDDFT excitation energies and solvatochromic shifts for large systems and that GEBF-CC2 could be used to validate GEBF-TDDFT result (with different functionals). The GEBF-TDDFT method is found to be able to provide satisfactory or reasonable descriptions on the experimental solvatochromic shifts for the n → π* transitions of acetone in various solutions, and the lowest π → π* transitions of pyridine and uracil in aqueous solutions.

Excited States of Xanthene Analogues: Photofragmentation and Calculations by CC2 and Time-Dependent Density Functional Theory.

Action spectroscopy has emerged as an analytical tool to probe excited states in the gas phase. Although comparison of gas-phase absorption properties with quantum-chemical calculations is, in principle, straightforward, popular methods often fail to describe many molecules of interest-such as xanthene analogues. We, therefore, face their nano- and picosecond laser-induced photofragmentation with excited-state computations by using the CC2 method and time-dependent density functional theory (TDDFT). Whereas the extracted absorption maxima agree with CC2 predictions, the TDDFT excitation energies are blueshifted. Lowering the amount of Hartree-Fock exchange in the DFT functional can reduce this shift but at the cost of changing the nature of the excited state. Additional bandwidth observed in the photofragmentation spectra is rationalized in terms of multiphoton processes. Observed fragmentation from higher-lying excited states conforms to intense excited-to-excited state transitions calculated with CC2. The CC2 method is thus suitable for the comparison with photofragmentation in xanthene analogues.

Solvent Effects on Electronic Excitations of an Organic Chromophore.

In this work we study the solvatochromic shift of a selected low-energy excited state of alizarin in water by using a linear-scaling implementation of large-scale time-dependent density functional theory (TDDFT). While alizarin, a small organic dye, is chosen as a simple example of solute-solvent interactions, the findings presented here have wider ramifications for the realistic modeling of dyes, paints, and pigment-protein complexes. We find that about 380 molecules of explicit water need to be considered in order to yield an accurate representation of the solute-solvent interaction and a reliable solvatochromic shift. By using a novel method of constraining the TDDFT excitation vector, we confirm that the origin of the slow convergence of the solvatochromic shift with system size is due to two different effects. The first factor is a strong redshift of the excitation due to an explicit delocalization of a small fraction of the electron and the hole from the alizarin onto the water, which is mainly confined to within a distance of 7 Å from the alizarin molecule. The second factor can be identified as long-range electrostatic influences of water molecules beyond the 7 Å region on the ground-state properties of alizarin. We also show that these electrostatic influences are not well reproduced by a QM/MM model, suggesting that full QM studies of relatively large systems may be necessary in order to obtain reliable results.

Why do TD-DFT excitation energies of BODIPY/Aza-BODIPY families largely deviate from experiment? Answers from electron correlated and multireference methods.

The vertical excitation energies of 17 boron-dipyrromethene (BODIPY) core structures with a variety of substituents and ring sizes are benchmarked using time-dependent density functional theory (TD-DFT) with nine different functionals combined with the cc-pVTZ basis set. When compared to experimental measurements, all functionals provide mean absolute errors (mean AEs) greater than 0.3 eV, larger than the 0.1-0.3 eV differences typically expected from TD-DFT. Due to the high linear correlation of TD-DFT results with experiment, most functionals can be used to predict excitation energies if corrected empirically. Using the CAM-B3LYP functional, 0-0 transition energies are determined, and while the absolute difference is improved (mean AE = 0.478 eV compared to 0.579 eV), the correlation diminishes substantially (R(2) = 0.961 to 0.862). Two very recently introduced charge transfer (CT) indices, q(CT) and d(CT), and electron density difference (EDD) plots demonstrate that CT does not play a significant role for most of the BODIPYs examined and, thus, cannot be the source of error in TD-DFT. To assess TD-DFT methods, vertical excitation energies are determined utilizing TD-HF, configuration interaction CIS and CIS(D), equation of motion EOM-CCSD, SAC-CI, and Laplace-transform based local coupled-cluster singles and approximate doubles LCC2* methods. Moreover, multireference CASSCF and CASPT2 vertical excitation energies were also obtained for all species (except CASPT2 was not feasible for the four largest systems). The SAC-CI/cc-pVDZ, LCC2*/cc-pVDZ, and CASPT2/cc-pVDZ approaches are shown to have the smallest mean AEs of 0.154, 0.109, and 0.100 eV, respectively; the utility of the LCC2* approach is demonstrated for eight extended BODIPYs and aza-BODIPYs. We found that the problems with TD-DFT arise from difficulties in dealing with the differential electron correlation (as assessed by comparing CCS, CC2, LR-CCSD, CCSDR(T), and CCSDR(3) vertical excitation energies for five compounds) and from contributions of multireference character and double excitations (from analysis of the CASSCF wave functions).

Insights into Magnetically Induced Current Pathways and Optical Properties of Isophlorins

The magnetically induced current density of tetraoxa-isophlorin and dioxa-dithia-isophlorin have been studied at the density functional theory (DFT) level using the gauge including magnetically induced current method (GIMIC). The current density calculations show that the studied isophlorins with formally 28 π electrons are strongly antiaromatic, sustaining paratropic ring currents of -48.5 and -58.3 nA/T, respectively. All chemical bonds of the porphyrinoid macroring participate in the transport of the paratropic ring current. Calculations of excitation energies at the time-dependent density functional theory (TDDFT) level and at correlated ab initio levels show that tetraoxa-isophlorin and dithia-isophlorin have small optical gaps of 0.82-1.34 and 0.90-1.25 eV, respectively. The transition to the lowest excited state, which belongs to the Bg irreducible representation, is dipole forbidden and has therefore not been observed in the spectroscopical studies. For dioxa-dithia-isophlorin, the excitation energies of the Q and B bands calculated at the TDDFT level agree rather well with experimental values with deviations in the range of [-0.15, +0.27] eV, whereas corresponding excitation energies obtained at the ab initio levels are 0.14-0.51 eV too large as compared to experimental values. For tetraoxa-isophlorin, the deviations of the TDDFT excitation energies from experimental values are in the range [-0.10, +0.50] eV and the ab initio values are 0.53-0.97 eV larger than the experimental values due to the significant double excitation character of the excited states.

Permanent dipole moments and energies of excited states from density functional theory compared with coupled cluster predictions: Case of para-nitroaniline

Different ways to extract properties of excited states from time-dependent density functional theory (TD-DFT) calculations are compared to ab initio results obtained with the Equation of Motion Coupled Cluster approach. The recently proposed a posteriori Tamm–Dancoff approximation (ATDA) predicts the permanent dipole moments to be underestimated by 25% on average, close to the results of the relaxed density TD-DFT formalism, quadratic response formalism, and numerical energy derivatives, while the unrelaxed density approximation results are less accurate (40% overestimate). We also propose a correction for TD-DFT excitation energies, which are known to be problematic for charge transfer states. The static DFT energies evaluated on the relaxed densities of the excited states are found to be more accurate than TD-DFT excitation energies (RMSD is 0.7eV vs. 1.1eV, while maximum deviation is −1.0eV vs. −2.0eV). This validates ATDA for description of nonlinear optical properties of donor–acceptor molecules, exemplified by para-nitroaniline, and extends this method to improve the excitation energy predictions.

Electronic Spectroscopy and Computational Studies of Glutathionylco(III)balamin

We have studied glutathionylcobalamin (GS-Cbl) by optical spectroscopy and with density functional theory (DFT) and time-dependent DFT (TD-DFT) electronic structure methods of truncated geometric models. We examined the geometric structure of the models by comparison of DFT calculations with recent high-resolution experimental X-ray structure data ( Hannibal, L. et al. Inorg. Chem. 2010, 49, 9921) for GS-Cbl, and we examined the TD-DFT excitation simulations by comparison of the models with measured optical spectra. The calculations employed the B3LYP hybrid functional and the nonhybrid BP86 functional in both vacuum and water (conductor polarized continuum model (cpcm)) with the 6-311G(d,p) basis set. The optimized geometric structure calculations for six truncated models were made by varying the chemical structure, solvent model, and the two DFT functionals. All showed similar geometry. Charge decomposition analysis (CDA) and extended charge decomposition analysis (ECDA), especially with BP86 shows the similar charge transfer nature of the Co-S bond in GS-Cbl and the Co-C bond in CH(3)Cbl. Mayer and Wiberg bond orders illustrate the similar covalent nature of the two bonds. Finally, absolute optical spectral simulation calculations were compared with the experimental UV-visible extinction spectrum and the electronic circular dichroism (ECD) differential extinction spectrum. The BP86 method shows more spectral features, and the best fit was found for a GS-Cbl model with 5,6-dimethylbenzimidazole at the BP86/6-311G(d,p) level with a water cpcm solvent model. The excited state transitions were investigated with Martin's natural transition orbitals (NTOs). The BP86 calculations also showed π bonding interactions between Co and the axial S of the GS- ligand in the molecular orbitals (MOs) and NTOs.

Response calculations based on an independent particle system with the exact one-particle density matrix: Excitation energies

Adiabatic response time-dependent density functional theory (TDDFT) suffers from the restriction to basically an occupied → virtual single excitation formulation. Adiabatic time-dependent density matrix functional theory allows to break away from this restriction. Problematic excitations for TDDFT, viz. bonding-antibonding, double, charge transfer, and higher excitations, are calculated along the bond-dissociation coordinate of the prototype molecules H(2) and HeH(+) using the recently developed adiabatic linear response phase-including (PI) natural orbital theory (PINO). The possibility to systematically increase the scope of the calculation from excitations out of (strongly) occupied into weakly occupied ("virtual") natural orbitals to larger ranges of excitations is explored. The quality of the PINO response calculations is already much improved over TDDFT even when the severest restriction is made, to virtually the size of the TDDFT diagonalization problem (only single excitation out of occupied orbitals plus all diagonal doubles). Further marked improvement is obtained with moderate extension to allow for excitation out of the lumo and lumo+1, which become fractionally occupied in particular at longer distances due to left-right correlation effects. In the second place the interpretation of density matrix response calculations is elucidated. The one-particle reduced density matrix response for an excitation is related to the transition density matrix to the corresponding excited state. The interpretation of the transition density matrix in terms of the familiar excitation character (single excitations, double excitations of various types, etc.) is detailed. The adiabatic PINO theory is shown to successfully resolve the problematic cases of adiabatic TDDFT when it uses a proper PI orbital functional such as the PILS functional.

Discontinuities of the exchange-correlation kernel and charge-transfer excitations in time-dependent density-functional theory

We identify the key property that the exchange-correlation (XC) kernel of time-dependent density-functional theory must have in order to describe long-range charge-transfer excitations. We show that the discontinuity of the XC potential as a function of particle number induces a space- and frequency-dependent discontinuity of the XC kernel that diverges as $r\ensuremath{\rightarrow}\ensuremath{\infty}$. In a combined donor-acceptor system, the same discontinuity compensates for the vanishing overlap between the acceptor and donor orbitals, thereby yielding a finite correction to the Kohn-Sham eigenvalue differences. This mechanism is illustrated to first order in the Coulomb interaction.

A block variational procedure for the iterative diagonalization of non-Hermitian random-phase approximation matrices

We present a technique for the iterative diagonalization of random-phase approximation (RPA) matrices, which are encountered in the framework of time-dependent density-functional theory (TDDFT) and the Bethe-Salpeter equation. The non-Hermitian character of these matrices does not permit a straightforward application of standard iterative techniques used, i.e., for the diagonalization of ground state Hamiltonians. We first introduce a new block variational principle for RPA matrices. We then develop an algorithm for the simultaneous calculation of multiple eigenvalues and eigenvectors, with convergence and stability properties similar to techniques used to iteratively diagonalize Hermitian matrices. The algorithm is validated for simple systems (Na(2) and Na(4)) and then used to compute multiple low-lying TDDFT excitation energies of the benzene molecule.