Efficient time-dependent density functional theory approximations for hybrid density functionals: analytical gradients and parallelization

In this paper, we present the implementation of efficient approximations to time-dependent density functional theory (TDDFT) within the Tamm-Dancoff approximation (TDA) for hybrid density functionals. For the calculation of the TDDFT/TDA excitation energies and analytical gradients, we combine the resolution of identity (RI-J) algorithm for the computation of the Coulomb terms and the recently introduced "chain of spheres exchange" (COSX) algorithm for the calculation of the exchange terms. It is shown that for extended basis sets, the RIJCOSX approximation leads to speedups of up to 2 orders of magnitude compared to traditional methods, as demonstrated for hydrocarbon chains. The accuracy of the adiabatic transition energies, excited state structures, and vibrational frequencies is assessed on a set of 27 excited states for 25 molecules with the configuration interaction singles and hybrid TDDFT/TDA methods using various basis sets. Compared to the canonical values, the typical error in transition energies is of the order of 0.01 eV. Similar to the ground-state results, excited state equilibrium geometries differ by less than 0.3 pm in the bond distances and 0.5° in the bond angles from the canonical values. The typical error in the calculated excited state normal coordinate displacements is of the order of 0.01, and relative error in the calculated excited state vibrational frequencies is less than 1%. The errors introduced by the RIJCOSX approximation are, thus, insignificant compared to the errors related to the approximate nature of the TDDFT methods and basis set truncation. For TDDFT/TDA energy and gradient calculations on Ag-TB2-helicate (156 atoms, 2732 basis functions), it is demonstrated that the COSX algorithm parallelizes almost perfectly (speedup ~26-29 for 30 processors). The exchange-correlation terms also parallelize well (speedup ~27-29 for 30 processors). The solution of the Z-vector equations shows a speedup of ~24 on 30 processors. The parallelization efficiency for the Coulomb terms can be somewhat smaller (speedup ~15-25 for 30 processors), but their contribution to the total calculation time is small. Thus, the parallel program completes a Becke3-Lee-Yang-Parr energy and gradient calculation on the Ag-TB2-helicate in less than 4 h on 30 processors. We also present the necessary extension of the Lagrangian formalism, which enables the calculation of the TDDFT excited state properties in the frozen-core approximation. The algorithms described in this work are implemented into the ORCA electronic structure system.     

Comparing modern density functionals for conjugated polymer band structures: screened hybrid, Minnesota, and Rung 3.5 approximations

Semiconducting polymers with π-conjugated backbones show promise in fields such as photovoltaics. Practical applications of conjugated polymers require precise control over the polymer's electronic band structure. Several new classes of density functional approximation, including screened hybrids, semilocal Minnesota functionals, and Rung 3.5 functionals, show potential for improved predictions of conjugated polymer band structures. This work compares these methods to standard global hybrid density functionals for bandgaps and band structures of representative conjugated polymers. The new methods exhibit particular promise for modeling three-dimensionally periodic bulk polymers, which can be problematic for global hybrids.     

Density functionals for exchange and correlation energies: Exact conditions and comparison of approximations

A set of exact conditions is compiled for the purpose of developing and testing approximations for the exchange-correlation energy as a functional of the electron density. Special emphasis is placed upon recently developed density-scaling relationships. Commonly used generalized gradient approximations are compared against several of these conditions. A direct tabular comparison of these functionals (not of calculated properties) with one another is also made. © 1994 John Wiley & Sons, Inc.     

Nature of one-dimensional short hydrogen bonding: bond distances, bond energies, and solvent effects

On the basis of recently synthesized calix[4]hydroquinone (CHQ) nanotubes which were self-assembled with infinitely long one-dimensional (1-D) short hydrogen bonds (SHB), we have investigated the nature of 1-D SHB using first-principles calculations for all the systems including the solvent water. The H-bonds relay (i.e., contiguous H-bonds) effect in CHQs shortens the H...O bond distances significantly (by more than 0.2 A) and increases the bond dissociation energy to a large extent (by more than approximately 4 kcal/mol) due to the highly enhanced polarization effect along the H-bond relay chain. The H-bonds relay effect shows a large increase in the chemical shift associated with the SHB. The average binding energies for the infinite 1-D H-bond arrays of dioles and dions increase by approximately 4 and approximately 9 kcal/mol per H-bond, respectively. The solvent effect (due to nonbridging water molecules) has been studied by explicitly adding water molecules in the CHQ tube crystals. This effect is found to be small with slight weakening of the SHB strength; the H...O bond distance increases only by 0.02 A, and the average binding energy decreases by approximately 1 kcal/mol per H-bond. All these results based on the first-principles calculations are the first detailed analysis of energy gain by SHB and energy loss by solvent effect, based on a partitioning scheme of the interaction energy components. These reliable results elucidate not only the self-assembly phenomena based on the H-bond relay but also the solvent effect on the SHB strength.     

On the accuracy of frozen density embedding calculations with hybrid and orbital-dependent functionals for non-bonded interaction energies

We analyze the accuracy of the frozen density embedding (FDE) method, with hybrid and orbital-dependent exchange-correlation functionals, for the calculation of the total interaction energies of weakly interacting systems. Our investigation is motivated by the fact that these approaches require, in addition to the non-additive kinetic energy approximation, also approximate non-additive exact-exchange energies. Despite this further approximation, we find that the hybrid/orbital-dependent FDE approaches can reproduce the total energies with the same accuracy (about 1 mHa) as the one of conventional semi-local functionals. In many cases, thanks to error cancellation effects, hybrid/orbital-dependent approaches yield even the smallest error. A detailed energy-decomposition investigation is presented. Finally, the Becke-exchange functional is found to reproduce accurately the non-additive exact-exchange energies also for non-equilibrium geometries. These performances are rationalized in terms of a reduced-gradient decomposition of the non-additive exchange energy.     

Assessment of new meta and hybrid meta density functionals for predicting the geometry and binding energy of a challenging system: the dimer of H2S and benzene

Noncovalent interactions of a hydrogen bond donor with an aromatic pi system present a challenge for density functional theory, and most density functionals do not perform well for this kind of interaction. Here we test seven recent density functionals from our research group, along with the popular B3LYP functional, for the dimer of H 2S with benzene. The functionals considered include the four new meta and hybrid meta density functionals of the M06 suite, three slightly older hybrid meta functionals, and the B3LYP hybrid functional, and they were tested for their abilities to predict the dissociation energies of three conformations of the H 2S-benzene dimer and to reproduce the key geometric parameters of the equilibrium conformation of this dimer. All of the functionals tested except B3LYP correctly predict which of the three conformations of the dimer is the most stable. The functionals that are best able to reproduce the geometry of the equilibrium conformation of the dimer with a polarized triple-zeta basis set are M06-L, PWB6K, and MPWB1K, each having a mean unsigned relative error across the two experimentally verifiable geometric parameters of only 8%. The success of M06-L is very encouraging because it is a local functional, which reduces the cost for large simulations. The M05-2X functional yields the most accurate binding energy of a conformation of the dimer for which a binding energy calculated at the CCSD(T) level of theory is available; M05-2X gives a binding energy for the system with a difference of merely 0.02 kcal/mol from that obtained by the CCSD(T) calculation. The M06 functional performs well in both categories by yielding a good representation of the geometry of the equilibrium structure and by giving a binding energy that is only 0.19 kcal/mol different from that calculated by CCSD(T). We conclude that the new generation of density functionals should be useful for a variety of problems in biochemistry and materials where aromatic functional groups can serve as hydrogen bond acceptors.     

Density functional resonance theory: complex density functions, convergence, orbital energies, and functionals

Aspects of density functional resonance theory (DFRT) [D. L. Whitenack and A. Wasserman, Phys. Rev. Lett. 107, 163002 (2011)], a recently developed complex-scaled version of ground-state density functional theory (DFT), are studied in detail. The asymptotic behavior of the complex density function is related to the complex resonance energy and system's threshold energy, and the function's local oscillatory behavior is connected with preferential directions of electron decay. Practical considerations for implementation of the theory are addressed including sensitivity to the complex-scaling parameter, θ. In Kohn-Sham DFRT, it is shown that almost all θ-dependence in the calculated energies and lifetimes can be extinguished via use of a proper basis set or fine grid. The highest occupied Kohn-Sham orbital energy and lifetime are related to physical affinity and width, and the threshold energy of the Kohn-Sham system is shown to be equal to the threshold energy of the interacting system shifted by a well-defined functional. Finally, various complex-scaling conditions are derived which relate the functionals of ground-state DFT to those of DFRT via proper scaling factors and a non-Hermitian coupling-constant system.     

Asymptotic properties of the exchange energy density and the exchange potential of finite systems: relevance for generalized gradient approximations

It is shown that generalized gradient approximations (GGAs) for exchange only, due to their very limited form, quite generally can not simultaneously reproduce both the asymptotic forms of the exchange energy density and the exchange potential of finite systems. Furthermore, mechanisms making GGAs formally approach at least one of these asymptotic forms do not improve the corresponding quantity in the relevant part of the asymptotic regime of atoms. By constructing a GGA which leads to superior atomic exchange energies compared to all GGAs heretofore but does not reproduce the asymptotic form of the exact exchange energy density it is demonstrated that this property is not important for obtaining extremely accurate atomic exchange energies. We conclude that GGAs by their very concept are not suited to reproduce these asymptotic properties of finite systems. As a byproduct of our discussion we present a particularly simple and direct proof of the well known asymptotic structure of the exchange potential of finite spherical systems.     

Calculation of electronic g-tensors for transition metal complexes using hybrid density functionals and atomic meanfield spin-orbit operators

We report the first implementation of the calculation of electronic g-tensors by density functional methods with hybrid functionals. Spin-orbit coupling is treated by the atomic meanfield approximation. g-Tensors for a set of small main group radicals and for a series of ten 3d and two 4d transition metal complexes have been compared using the local density approximation (VWN functional), the generalized gradient approximation (BP86 functional), as well as B3-type (B3PW91) and BH-type (BHPW91) hybrid functionals. For main group radicals, the effect of exact-exchange mixing is small. In contrast, significant differences between the various functionals arise for transition metal complexes. As has been shown previously, local and in particular gradient-corrected functionals tend to underestimate the "paramagnetic" contributions to the g-tensors in these cases and thereby recover only about 40-50% of the range of experimental g-tensor components. This is improved to ca. 60% by the B3PW91 functional, which also gives slightly reduced standard deviations. The range increases to almost 100% using the half-and-half functional BHPW91. However, the quality of the correlation with experimental data worsens due to a significant overestimate of some intermediate g-tensor values. The worse performance of the BHPW91 functional in these cases is accompanied by spin contamination. Although none of the functionals tested thus appears to be ideal for the treatment of electronic g-tensors in transition metal complexes, the B3PW91 hybrid functional exhibited the overall most satisfactory performance. Apart from the validation of hybrid functionals, some aspects in the treatment of spin-orbit contributions to the g-tensor are discussed.     

A variational principle for transition and density matrices and approximations to the equations-of-motion

A general variational principle for transition and density matrices is proposed. The principle is closely related to Rowe's variational treatment of the equations-of-motion method. It permits the simultaneous construction of coupled approximations for two eigenstates, and it is a straightforward extension of the usual variational method.     

Long-range and short-range Coulomb correlation effects as simulated by Hartree–Fock,local density approximation,and generalized gradient approximation exchange functionals

 Exchange functionals used in density functional theory (DFT) are generally considered to simulate long-range electron correlation effects. It is shown that these effects can be traced back to the self-interaction error (SIE) of approximate exchange functionals. An analysis of the SIE with the help of the exchange hole reveals that both short-range (dynamic) and long-range (nondynamic) electron correlation effects are simulated by DFT exchange where the local density approximation (LDA) accounts for stronger effects than the generalized gradient expansion (GGA). This is a result of the fact that the GGA exchange hole describes the exact exchange hole close to the reference electron more accurately than the LDA hole does. The LDA hole is more diffuse, thus leading to an underestimation of exchange and stronger SIE effects, where the magnitude of the SIE energy is primarily due to the contribution of the core orbitals. The GGA exchange hole is more compact, which leads to an exaggeration of exchange in the bond and the nonbonding region and negative SIE contributions. Partitioning of the SIE into intra-/interelectronic and individual orbital contributions makes it possible to test the performance of a given exchange functional in different regions of the molecule. It is shown that Hartree–Fock exchange always covers some long-range effects via interelectronic exchange while self-interaction-corrected DFT is lacking these effects. Received: 25 May 2002 / Accepted: 7 October 2002 / Published online: 21 January 2003 Correspondence to: E. Kraka e-mail: kraka@theoc.gn.se Acknowledgements. This work was supported financially by the Swedish Natural Science Research Council (NFR). Calculations were done on the supercomputers of Nationellt Superdatorcentrum (NSC), Link?ping, Sweden. The authors thank the NSC for a generous allotment of computer time.     

Analytic derivatives for perturbatively corrected "double hybrid" density functionals: theory, implementation, and applications

A recently proposed new family of density functionals [S. Grimme, J. Chem. Phys. 124, 34108 (2006)] adds a fraction of nonlocal correlation as a new ingredient to density functional theory (DFT). This fractional correlation energy is calculated at the level of second-order many-body perturbation theory (PT2) and replaces some of the semilocal DFT correlation of standard hybrid DFT methods. The new "double hybrid" functionals (termed, e.g., B2-PLYP) contain only two empirical parameters that have been adjusted in thermochemical calculations on parts of the G2/3 benchmark set. The methods have provided the lowest errors ever obtained by any DFT method for the full G3 set of molecules. In this work, the applicability of the new functionals is extended to the exploration of potential energy surfaces with analytic gradients. The theory of the analytic gradient largely follows the standard theory of PT2 gradients with some additional subtleties due to the presence of the exchange-correlation terms in the self-consistent field operator. An implementation is reported for closed-shell as well as spin-unrestricted reference determinants. Furthermore, the implementation includes external point charge fields and also accommodates continuum solvation models at the level of the conductor like screening model. The density fitting resolution of the identity (RI) approximation can be applied to the evaluation of the PT2 part with large gains in computational efficiency. For systems with approximately 500-600 basis functions the evaluation of the double hybrid gradient is approximately four times more expensive than the calculation of the standard hybrid DFT gradient. Extensive test calculations are provided for main group elements and transition metal containing species. The results reveal that the B2-PLYP functional provides excellent molecular geometries that are superior compared to those from standard DFT and MP2.     

Higher approximations for transition matrices and their application to the calculation of atomic spectra

Theoretical and Experimental Chemistry -     

The M06 suite of density functionals for main group thermochemistry, thermochemical kinetics, noncovalent interactions, excited states, and transition elements: two new functionals and systematic testing of four M06-class functionals and 12 other functionals

We present two new hybrid meta exchange- correlation functionals, called M06 and M06-2X. The M06 functional is parametrized including both transition metals and nonmetals, whereas the M06-2X functional is a high-nonlocality functional with double the amount of nonlocal exchange (2X), and it is parametrized only for nonmetals.The functionals, along with the previously published M06-L local functional and the M06-HF full-Hartree–Fock functionals, constitute the M06 suite of complementary functionals. We assess these four functionals by comparing their performance to that of 12 other functionals and Hartree–Fock theory for 403 energetic data in 29 diverse databases, including ten databases for thermochemistry, four databases for kinetics, eight databases for noncovalent interactions, three databases for transition metal bonding, one database for metal atom excitation energies, and three databases for molecular excitation energies. We also illustrate the performance of these 17 methods for three databases containing 40 bond lengths and for databases containing 38 vibrational frequencies and 15 vibrational zero point energies. We recommend the M06-2X functional for applications involving main-group thermochemistry, kinetics, noncovalent interactions, and electronic excitation energies to valence and Rydberg states. We recommend the M06 functional for application in organometallic and inorganometallic chemistry and for noncovalent interactions. Electronic Supplementary Material The online version of this article (doi:) contains supplementary material, which is available to authorized users. Contribution to the Mark S. Gordon 65th Birthday Festschrift Issue.     

Evaluation of range-separated hybrid density functionals for the prediction of vibrational frequencies, infrared intensities, and Raman activities

We present an assessment of different density functionals, with emphasis on range-separated hybrids, for the prediction of fundamental and harmonic vibrational frequencies, infrared intensities, and Raman activities. Additionally, we discuss the basis set convergence of vibrational properties of H2O with long-range corrected hybrids. Our results show that B3LYP is the best functional for predicting vibrational frequencies (both fundamental and harmonic); the screened-PBE hybrid (HSE) density functional works best for infrared intensities, and the long-range corrected PBE (LC-omegaPBE), M06-HF, and M06-L density functionals are almost as good as MP2 for predicting Raman activities. We show the predicted Raman spectrum of adenine as an example of a medium-size molecule where a DFT/Sadlej pVTZ calculation is affordable and compare our results against the experimental spectrum.     

TiCl, TiH, and TiH+ bond energies: a test of a correlation-consistent Ti basis set

Correlation-consistent basis sets are developed for the Ti atom. The polarization functions are optimized for the average of the ³F and ⁵F states. One series of correlation-consistent basis sets is for 3d and 4s correlation, while the second series includes 3s and 3p correlation as well as 3d and 4s correlation. These basis sets are tested using the Ti ³F–⁵F separation and the dissociation energies of TiCl X⁴Φ, TiH X⁴Φ, and TiH⁺ X³Φ. The CCSD(T) complete basis set limit values are determined by extrapolation. The Douglas–Kroll approach is used to compute the scalar relativistic effect. Spin-orbit effects are taken from experiment and/or are computed at the CASSCF level. The Ti³F–⁵F separation is in excellent agreement with experiment, while the TiCl, TiH, and TiH⁺ bond energies are in good agreement with experiment. Extrapolation with the valence basis set is consistent with other atoms, while including 3s and 3p correlation appears to make extrapolation more difficult. Received: 20 January 1999 / Accepted: 26 February 1999 / Published online: 7 June 1999     

Local and density fitting approximations within the short-range/long-range hybrid scheme: application to large non-bonded complexes

We introduce local and density fitting approximations into a hybrid approach which couples short-range density functionals with long-range wave-functions. The use of density fitting methods within the local-correlation approximation makes it possible to tackle bigger systems than without these methods, so that we are now able to treat systems of biological interest within the DFT/ab initio scheme. In first benchmark calculations, we apply it to the S22 database of Hobza and co-workers for binding energies of weakly bound molecular clusters.     

Assessment of density functional theory optimized basis sets for gradient corrected functionals to transition metal systems: the case of small Nin (n<or=5) clusters

Density functional calculations have been performed for small nickel clusters, Ni(n), Ni(n) (+), and Ni(n)(-) (n相似文献   

A new local density functional for main-group thermochemistry, transition metal bonding, thermochemical kinetics, and noncovalent interactions

We present a new local density functional, called M06-L, for main-group and transition element thermochemistry, thermochemical kinetics, and noncovalent interactions. The functional is designed to capture the main dependence of the exchange-correlation energy on local spin density, spin density gradient, and spin kinetic energy density, and it is parametrized to satisfy the uniform-electron-gas limit and to have good performance for both main-group chemistry and transition metal chemistry. The M06-L functional and 14 other functionals have been comparatively assessed against 22 energetic databases. Among the tested functionals, which include the popular B3LYP, BLYP, and BP86 functionals as well as our previous M05 functional, the M06-L functional gives the best overall performance for a combination of main-group thermochemistry, thermochemical kinetics, and organometallic, inorganometallic, biological, and noncovalent interactions. It also does very well for predicting geometries and vibrational frequencies. Because of the computational advantages of local functionals, the present functional should be very useful for many applications in chemistry, especially for simulations on moderate-sized and large systems and when long time scales must be addressed.