Analytical time-dependent density functional derivative methods within the RI-J approximation, an approach to excited states of large molecules

Time-dependent density functional theory (TDDFT) is now well established as an efficient method for molecular excited state treatments. In this work, we introduce the resolution of the identity approximation for the Coulomb energy (RI-J) to excited state gradient calculations. In combination with nonhybrid functionals, the RI-J approximation leads to speed ups in total timings of an order of magnitude compared to the conventional method; this is demonstrated for oligothiophenes with up to 40 monomeric units and adamantane clusters. We assess the accuracy of the computed adiabatic excitation energies, excited state structures, and vibrational frequencies on a set of 36 excited states. The error introduced by the RI-J approximation is found to be negligible compared to deficiencies of standard basis sets and functionals. Auxiliary basis sets optimized for ground states are suitable for excited state calculations with small modifications. In conclusion, the RI-J approximation significantly extends the scope of applications of analytical TDDFT derivative methods in photophysics and photochemistry.     

Pseudospectral time-dependent density functional theory

Time-dependent density functional theory (TDDFT) is implemented within the Tamm-Dancoff approximation (TDA) using a pseudospectral approach to evaluate two-electron repulsion integrals. The pseudospectral approximation uses a split representation with both spectral basis functions and a physical space grid to achieve a reduction in the scaling behavior of electronic structure methods. We demonstrate here that exceptionally sparse grids may be used in the excitation energy calculation, following earlier work employing the pseudospectral approximation for determining correlation energies in wavefunction-based methods with similar conclusions. The pseudospectral TDA-TDDFT method is shown to be up to ten times faster than a conventional algorithm for hybrid functionals without sacrificing chemical accuracy.     

Time-dependent density functional theory based upon the fragment molecular orbital method

Time-dependent density functional theory (TDDFT) was combined with the two-body fragment molecular orbital method (FMO2). In this FMO2-TDDFT scheme, the system is divided into fragments, and the electron density for fragments is determined self-consistently. Consequently, only one main fragment of interest and several fragment pairs including it are calculated by TDDFT. To demonstrate the accuracy of FMO2-TDDFT, we computed several low-lying singlet and triplet excited states of solvated phenol and polyalanine using our method and the standard TDDFT for the full system. The BLYP functional with the long-range correction (LC-BLYP) was employed with the 6-31G(*) basis set (some tests were also performed with 6-311G(*), as well as with B3LYP and time-dependent Hartree-Fock). Typically, FMO2-TDDFT reproduced the full TDDFT excitation energies within 0.1 eV, and for one excited state the error was about 0.2 eV. Beside the accurate reproduction of the TDDFT excitation energies, we also automatically get an excitation energy decomposition analysis, which provides the contributions of individual fragments. Finally, the efficiency of our approach was exemplified on the LC-BLYP6-31G(*) calculation of the lowest singlet excitation of the photoactive yellow protein which consists of 1931 atoms, and the obtained value of 3.1 eV is in agreement with the experimental value of 2.8 eV.     

Analytical Hessian of electronic excited states in time-dependent density functional theory with Tamm-Dancoff approximation

We present the analytical expression and computer implementation for the second-order energy derivatives of the electronic excited state with respect to the nuclear coordinates in the time-dependent density functional theory (TDDFT) with Gaussian atomic orbital basis sets. Here, the Tamm-Dancoff approximation to the full TDDFT is adopted, and therefore the formulation process of TDDFT excited-state Hessian is similar to that of configuration interaction singles (CIS) Hessian. However, due to the replacement of the Hartree-Fock exchange integrals in CIS with the exchange-correlation kernels in TDDFT, many quantitative changes in the derived equations are arisen. The replacement also causes additional technical difficulties associated with the calculation of a large number of multiple-order functional derivatives with respect to the density variables and the nuclear coordinates. Numerical tests on a set of test molecules are performed. The simulated excited-state vibrational frequencies by the analytical Hessian approach are compared with those computed by CIS and the finite-difference method. It is found that the analytical Hessian method is superior to the finite-difference method in terms of the computational accuracy and efficiency. The numerical differentiation can be difficult due to root flipping for excited states that are close in energy. TDDFT yields more exact excited-state vibrational frequencies than CIS, which usually overestimates the values.     

Failure of time-dependent density functional theory for long-range charge-transfer excited states: the zincbacteriochlorin-bacteriochlorin and bacteriochlorophyll-spheroidene complexes

It is well-known that time-dependent density functional theory (TDDFT) yields substantial errors for the excitation energies of charge-transfer (CT) excited states, when approximate standard exchange-correlation (xc) functionals are used, for example, SVWN, BLYP, or B3LYP. Also, the correct 1/R asymptotic behavior of CT states with respect to a distance coordinate R between the separated charges of the CT state is not reproduced by TDDFT employing these xc-functionals. Here, we demonstrate by analysis of the TDDFT equations that the first failure is due to the self-interaction error in the orbital energies from the ground-state DFT calculation, while the latter is a similar self-interaction error in TDDFT arising through the electron transfer in the CT state. Possible correction schemes, such as inclusion of exact Hartree-Fock or exact Kohn-Sham exchange, as well as aspects of the exact xc-functional are discussed in this context. Furthermore, a practical approach is proposed which combines the benefits of TDDFT and configuration interaction singles (CIS) and which does not suffer from electron-transfer self-interaction. The latter approach is applied to a (1,4)-phenylene-linked zincbacteriochlorin-bacteriochlorin complex and to a bacteriochlorophyll-spheroidene complex, in which CT states may play important roles in energy and electron-transfer processes. The errors of TDDFT alone for the CT states are demonstrated, and reasonable estimates for the true excitation energies of these states are given.     

Highly efficient implementation of pseudospectral time‐dependent density‐functional theory for the calculation of excitation energies of large molecules

We have developed and implemented pseudospectral time‐dependent density‐functional theory (TDDFT) in the quantum mechanics package Jaguar to calculate restricted singlet and restricted triplet, as well as unrestricted excitation energies with either full linear response (FLR) or the Tamm–Dancoff approximation (TDA) with the pseudospectral length scales, pseudospectral atomic corrections, and pseudospectral multigrid strategy included in the implementations to improve the chemical accuracy and to speed the pseudospectral calculations. The calculations based on pseudospectral time‐dependent density‐functional theory with full linear response (PS‐FLR‐TDDFT) and within the Tamm–Dancoff approximation (PS‐TDA‐TDDFT) for G2 set molecules using B3LYP/6‐31G*^* show mean and maximum absolute deviations of 0.0015 eV and 0.0081 eV, 0.0007 eV and 0.0064 eV, 0.0004 eV and 0.0022 eV for restricted singlet excitation energies, restricted triplet excitation energies, and unrestricted excitation energies, respectively; compared with the results calculated from the conventional spectral method. The application of PS‐FLR‐TDDFT to OLED molecules and organic dyes, as well as the comparisons for results calculated from PS‐FLR‐TDDFT and best estimations demonstrate that the accuracy of both PS‐FLR‐TDDFT and PS‐TDA‐TDDFT. Calculations for a set of medium‐sized molecules, including C_n fullerenes and nanotubes, using the B3LYP functional and 6‐31G^** basis set show PS‐TDA‐TDDFT provides 19‐ to 34‐fold speedups for C_n fullerenes with 450–1470 basis functions, 11‐ to 32‐fold speedups for nanotubes with 660–3180 basis functions, and 9‐ to 16‐fold speedups for organic molecules with 540–1340 basis functions compared to fully analytic calculations without sacrificing chemical accuracy. The calculations on a set of larger molecules, including the antibiotic drug Ramoplanin, the 46‐residue crambin protein, fullerenes up to C₅₄₀ and nanotubes up to 14×(6,6), using the B3LYP functional and 6‐31G^** basis set with up to 8100 basis functions show that PS‐FLR‐TDDFT CPU time scales as N^2.05 with the number of basis functions. © 2016 Wiley Periodicals, Inc.     

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)相似文献   

Analytical approach for the excited-state Hessian in time-dependent density functional theory: formalism, implementation, and performance

The paper presents the formalism, implementation, and performance of the analytical approach for the excited-state Hessian in the time-dependent density functional theory (TDDFT) that extends our previous work [J. Liu and W. Z. Liang, J. Chem. Phys. 135, 014113 (2011)] on the analytical Hessian in TDDFT within Tamm-Dancoff approximation (TDA) to full TDDFT. In contrast to TDA-TDDFT, an appreciable advantage of full TDDFT is that it maintains the oscillator strength sum rule, and therefore yields more precise results for the oscillator strength and other related physical quantities. For the excited-state harmonic vibrational frequency calculation, however, full TDDFT does not seem to be advantageous since the numerical tests demonstrate that the accuracy of TDDFT with and without TDA are comparable to each other. As a common practice, the computed harmonic vibrational frequencies are scaled by a suitable scale factor to yield good agreement with the experimental fundamental frequencies. Here we apply both the optimized ground-state and excited-state scale factors to scale the calculated excited-state harmonic frequencies and find that the scaling decreases the root-mean-square errors. The optimized scale factors derived from the excited-state calculations are slightly smaller than those from the ground-state calculations.     

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.     

Time-dependent density functional theory: past, present, and future

Time-dependent density functional theory (TDDFT) is presently enjoying enormous popularity in quantum chemistry, as a useful tool for extracting electronic excited state energies. This article discusses how TDDFT is much broader in scope, and yields predictions for many more properties. We discuss some of the challenges involved in making accurate predictions for these properties.     

A simplified relativistic time-dependent density-functional theory formalism for the calculations of excitation energies including spin-orbit coupling effect

In the present work we have proposed an approximate time-dependent density-functional theory (TDDFT) formalism to deal with the influence of spin-orbit coupling effect on the excitation energies for closed-shell systems. In this formalism scalar relativistic TDDFT calculations are first performed to determine the lowest single-group excited states and the spin-orbit coupling operator is applied to these single-group excited states to obtain the excitation energies with spin-orbit coupling effects included. The computational effort of the present method is much smaller than that of the two-component TDDFT formalism and this method can be applied to medium-size systems containing heavy elements. The compositions of the double-group excited states in terms of single-group singlet and triplet excited states are obtained automatically from the calculations. The calculated excitation energies based on the present formalism show that this formalism affords reasonable excitation energies for transitions not involving 5p and 6p orbitals. For transitions involving 5p orbitals, one can still obtain acceptable results for excitations with a small truncation error, while the formalism will fail for transitions involving 6p orbitals, especially 6p1/2 spinors.     

Performance of parallel TURBOMOLE for density functional calculations

The parallelization of density functional treatments of molecular electronic energy and first-order gradients is described, and the performance is documented. The quadrature required for exchange correlation terms and the treatment of exact Coulomb interaction scales virtually linearly up to 100 nodes. The RI-J technique to approximate Coulomb interactions (by means of an auxiliary basis set approximation for the electron density) even shows superlinear speedup on distributed memory architectures. The bottleneck is then linear algebra. Demonstrative application examples include molecules with up to 300 atoms and 3000 basis functions that can now be treated in a few hours per geometry optimization cycle in C₁ symmetry. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 1746–1757, 1998     

Theoretical study of the electronic spectra of square-planar platinum (II) complexes based on the two-component relativistic time-dependent density-functional theory

In the present work the electronic spectra of [PtCl(4)](2-), [PtBr(4)](2-), and [Pt(CN)(4)](2-) are studied with a recently proposed relativistic time-dependent density-functional theory (TDDFT) based on the two-component zeroth-order regular approximation and a noncollinear exchange-correlation (XC) functional. The contribution to the double group excited states in terms of singlet and triplet single group excited states is estimated through the inner product of the transition density matrix obtained from two-component and scalar relativistic TDDFT calculations to better understand the double group excited states. Spin-orbital coupling effects are found to be very important in order to simulate the electronic spectra of these complexes. The results show that the two-component TDDFT formalism can afford excitation energies with high accuracy for the transition-metal systems studied here when use is made of a proper XC potential.     

Time-dependent density functional theory based on a noncollinear formulation of the exchange-correlation potential

In this study we have introduced a formulation of time-dependent density functional theory (TDDFT) based on a noncollinear exchange-correlation potential. This formulation is a generalization of conventional TDDFT. The form of this formulation is exactly the same as that of the conventional TDDFT for the excitation energies of transitions that do not involve spin flips. In addition, this noncollinear TDDFT formulation allows for spin-flip transitions. This feature makes it possible to resolve more fully excited state spin multiplets, while for closed-shell systems, the spin-flip transitions will result in singlet-triplet excitations and this excitation energy calculated from this formulation of TDDFT is exactly the same as that from ordinary TDDFT. This formulation is applied to the dissociation of H(2) in its (1)Sigma(g) (+) ground state and (1)Sigma(u) (+) and (3)Sigma(u) (-) excited states with (3)Sigma(u) (-) (M(s)=+1) as the reference state and the multiplets splitting of some atoms.     

Coherent control and time-dependent density functional theory: towards creation of wave packets by ultrashort laser pulses

Explicitly time-dependent density functional theory (TDDFT) is a formally exact theory, which can treat very large systems. However, in practice it is used almost exclusively in the adiabatic approximation and with standard ground state functionals. Therefore, if combined with coherent control theory, it is not clear which control tasks can be achieved reliably, and how this depends on the functionals. In this paper, we continue earlier work in order to establish rules that answer these questions. Specifically, we look at the creation of wave packets by ultrashort laser pulses that contain several excited states. We find that (i) adiabatic TDDFT only works if the system is not driven too far from the ground state, (ii) the permanent dipole moments involved should not differ too much, and (iii) these results are independent of the functional used. Additionally, we find an artifact that produces fluence-dependent excitation energies.     

An overlap fitted chain of spheres exchange method

The "chain of spheres" (COS) algorithm, as part of the RIJCOSX SCF procedure, approximates the exchange term by performing analytic integration with respect to the coordinates of only one of the two electrons, whereas for the remaining coordinates, integration is carried out numerically. In the present work, we attempt to enhance the efficiency of the method by minimizing numerical errors in the COS procedure. The main idea is based on the work of Friesner and consists of finding a fitting matrix, Q, which leads the numerical and analytically evaluated overlap matrices to coincide. Using Q, the evaluation of exchange integrals can indeed be improved. Improved results and timings are obtained with the present default grid setup for both single point calculations and geometry optimizations. The fitting procedure results in a reduction of grid sizes necessary for achieving chemical accuracy. We demonstrate this by testing a number of grids and comparing results to the fully analytic and the earlier COS approximations. This turns out to be favourable for total and reaction energies, for which chemical accuracy can now be reached with a corresponding ~30% speedup over the original RIJCOSX procedure for single point energies. Results are slightly less favourable for the accuracy of geometry optimizations, but the procedure is still shown to yield geometries with errors well below the method inherent errors of the employed theoretical framework.     

Benchmarking the performance of time-dependent density functional methods

The performance of 24 density functionals, including 14 meta-generalized gradient approximation (mGGA) functionals, is assessed for the calculation of vertical excitation energies against an experimental benchmark set comprising 14 small- to medium-sized compounds with 101 total excited states. The experimental benchmark set consists of singlet, triplet, valence, and Rydberg excited states. The global-hybrid (GH) version of the Perdew-Burke-Ernzerhoff GGA density functional (PBE0) is found to offer the best overall performance with a mean absolute error (MAE) of 0.28 eV. The GH-mGGA Minnesota 2006 density functional with 54% Hartree-Fock exchange (M06-2X) gives a lower MAE of 0.26 eV, but this functional encounters some convergence problems in the ground state. The local density approximation functional consisting of the Slater exchange and Volk-Wilk-Nusair correlation functional (SVWN) outperformed all non-GH GGAs tested. The best pure density functional performance is obtained with the local version of the Minnesota 2006 mGGA density functional (M06-L) with an MAE of 0.41 eV.