Local hybrid functionals: an assessment for thermochemical kinetics

Local hybrid functionals with position-dependent exact-exchange admixture are a new class of exchange-correlation functionals in density functional theory that promise to advance the available accuracy in many areas of application. Local hybrids with different local mixing functions (LMFs) governing the position dependence are validated for the heats of formation of the extended G3/99 set, and for two sets of barriers of hydrogen-transfer and heavy-atom transfer reactions (HTBH38 and NHTBH38 databases). A simple local hybrid Lh-SVWN with only Slater and exact exchange plus local correlation and a one-parameter LMF, g(r)=b(tau(W)(r)tau(r)), performs best and provides overall mean absolute errors for thermochemistry and kinetics that are a significant improvement over standard state-of-the-art global hybrid functionals. In particular, this local hybrid functional does not suffer from the systematic deterioration that standard functionals exhibit for larger molecules. In contrast, local hybrids based on generalized gradient approximation exchange tend to give rise to nonintuitive LMFs, and no improved functionals have been obtained along this route. The LMF is a real-space function and thus can be analyzed in detail. We use, in particular, graphical analyses to rationalize the performance of different local hybrids for thermochemistry and reaction barriers.     

Importance of the correlation contribution for local hybrid functionals: range separation and self-interaction corrections

Local hybrid functionals with their position-dependent exact-exchange admixture are a conceptually simple and promising extension of the concept of a hybrid functional. Local hybrids based on a simple mixing of the local spin density approximation (LSDA) with exact exchange have been shown to be successful for thermochemistry, reaction barriers, and a range of other properties. So far, the combination of this generation of local hybrids with an LSDA correlation functional has been found to give the most favorable results for atomization energies, for a range of local mixing functions (LMFs) governing the exact-exchange admixture. Here, we show that the choice of correlation functional to be used with local hybrid exchange crucially influences the parameterization also of the exchange part as well as the overall performance. A novel ansatz for the correlation part of local hybrids is suggested based on (i) range-separation of LSDA correlation into short-range (SR) and long-range (LR) parts, and (ii) partial or full elimination of the one-electron self-correlation from the SR part. It is shown that such modified correlation functionals allow overall larger exact exchange admixture in thermochemically competitive local hybrids than before. This results in improvements for reaction barriers and for other properties crucially influenced by self-interaction errors, as demonstrated by a number of examples. Based on the range-separation approach, a fresh view on the breakdown of the correlation energy into dynamical and non-dynamical parts is suggested.     

What can we learn from the adiabatic connection formalism about local hybrid functionals?

Local hybrid functionals with position-dependent exact-exchange admixture are a promising new generation of exchange-correlation functionals for a large variety of applications. So far, the local mixing functions (LMFs) determining the position dependence have been largely constructed in an ad hoc manner, albeit based on physical reasoning. Here the basic formalism of the adiabatic connection is employed to investigate the formal basis of local hybrids and to construct a priori LMFs. Both a local spin density approximation to the LMF (AC-LSDA LMF) and generalized gradient approximation approximations (AC-PW91 LMF and AC-PBE LMF) turn out to provide inferior performance when used in local hybrids to compute atomization energies and reaction barriers compared to previous semiempirical LMFs. This is rationalized by limited flexibility of these first-principles LMFs and some basic limitations of the adiabatic connection formalism in this context. Graphical analyses and formal considerations provide nevertheless important new insight into the physical background of local hybrid functionals.     

Comparative density-functional study of the electron paramagnetic resonance parameters of amavadin

Electronic g tensors and hyperfine coupling tensors have been calculated for amavadin, an unusual eight-coordinate vanadium(IV) complex isolated from Amanita muscaria mushrooms. Different density-functional methods have been compared, ranging from local via gradient-corrected to hybrid functionals with a variable Hartree-Fock exchange admixture. For both electron paramagnetic resonance (EPR) properties, hybrid functionals with an appreciable exact-exchange admixture provide the closest agreement with experimental data. Second-order spin-orbit corrections provide non-negligible contributions to the 51V hyperfine tensor. The orientation of g and A tensors relative to each other also depends on spin-orbit corrections to the A tensor. A rationalization for the close resemblance of the EPR parameters of amavadin to those of the structurally rather different vanadyl complexes is provided, based on the nature of the relevant frontier orbitals.     

A thermochemically competitive local hybrid functional without gradient corrections

Following the suggestion of local hybrid functionals with position-dependent exact-exchange admixture [J. Jaramillo, G. E. Scuseria, and M. Ernzerhof, J. Chem. Phys. 118, 1068 (2003)], a functional that mixes only local and exact exchange plus local correlation has been constructed. With a simple local mixing function for the position dependence, this Lh-SVWN functional provides atomization energies for the G2-1 set that are competitive with currently available state-of-the-art functionals like, e.g., B3LYP. This is achieved without generalized gradient approximations for exchange or correlation.     

Evaluation of a combination of local hybrid functionals with DFT-D3 corrections for the calculation of thermochemical and kinetic data

Due to their position-dependent exact exchange admixture, local hybrid functionals offer a higher flexibility and thus the potential for more universal and accurate exchange correlation functionals compared to global hybrids with a constant admixture, as has been demonstrated in previous work. Yet, the local hybrid constructions used so far do not account for the inclusion of dispersion-type interactions. As a first exploratory step toward a more general approach that includes van der Waals-type interactions with local hybrids, the present work has added DFT-D3-type corrections to a number of simple local hybrid functionals. Optimization of only the s(8) and s(r,6) parameters for the S22 set provides good results for weak interaction energies but deteriorates the excellent performance of the local hybrids for G3 atomization energies and for classical reaction barriers. A combined optimization of the two DFT-D3 parameters with one of the two parameters of the spin-polarized local mixing function (LMF) of a local hybrid for a more general optimization set provides simultaneously accurate dispersion energies, improved atomization energies, and accurate reaction barriers, as well as excellent alkane protobranching ratios. For other LMFs, the improvements of such a combined optimization for the S22 energies have been less satisfactory. The most notable advantage of the dispersion-corrected local hybrids over, for example, a B3LYP-D3 approach, is in the much more accurate reaction barriers.     

Hybrid exchange correlation functionals and potentials: Concept elaboration

This paper deals with hybrid functionals that contain exact exchange energy and are the most popular and effective functionals in modern density functional theory. Emphasis is laid on generalization of the notion of a “hybrid functional,” which arises from the introduction of the spatial dependence of the exact exchange admixture (local hybrid functionals). Problems inherent in hybrid functionals are considered along with problems inherent in a wider class of so-called orbital-dependent functionals. In particular, the technique for constructing the local and multiplicative potentials, including the optimized effective potential method, is considered in detail. The theoretical approaches under study are illustrated by calculations of atomization molecular energies and magnetic resonance parameters.     

Effective local potentials for excited states

The constrained variational Hartree-Fock method for excited states of the same symmetry as the ground state [Chem. Phys. Lett. 287, 189 (1998)] is combined with the effective local potential (ELP) method [J. Chem. Phys. 125, 081104 (2006)] to generate Kohn-Sham-type exact-exchange potentials for singly excited states of many-electron systems. Illustrative examples include the three lowest (2)S states of the Li and Na atoms and the three lowest (3)S states of He and Be. For the systems studied, excited-state ELPs differ from the corresponding ground-state potentials in two respects: They are less negative and have small additional "bumps" in the outer electron region. The technique is general and can be used to approximate excited-state exchange-correlation potentials for other orbital-dependent functionals.     

Structure and electron paramagnetic resonance parameters of the manganese site of concanavalin A studied by density functional methods

The EPR parameters of the manganese site in the saccharide-binding protein concanavalin A have been studied by density functional methods, with an emphasis on metal (⁵⁵Mn) and ligand (¹H and ¹⁷O) hyperfine couplings, in comparison with high-field EPR and ENDOR data. Results for gradient-corrected and hybrid functionals with different exact-exchange admixture have been compared with experiment for the ⁵⁵Mn and the ¹H ligand hyperfine coupling and have been predicted for ¹⁷O hyperfine coupling based on comparison with experiment for the related [Mn(H₂O)₆]²⁺. Appreciable exact-exchange admixture in the hybrid functional is needed to obtain an adequate spin-density distribution and thus near-quantitative agreement with experimental EPR parameters. The common use of experimental proton hyperfine coupling tensors together with the point-dipole approximation for determination of bond lengths is evaluated by explicit calculations.     

Self-consistent effective local potentials

An effective local potential (ELP) is a multiplicative operator whose deviation from a given nonlocal potential has the smallest variance evaluated with a prescribed single-determinant wave function. ELPs are useful in density functional theory as alternatives to optimized effective potentials (OEPs) because they do not require special treatment in finite basis set calculations as OEPs do. We generalize the idea of variance-minimizing potentials by introducing the concept of a self-consistent ELP (SCELP), a local potential whose deviation from its nonlocal counterpart has the smallest variance in terms of its own Kohn-Sham orbitals. A semi-analytical method for computing SCELPs is presented. The OEP, ELP, and SCELP techniques are applied to the exact-exchange-only Kohn-Sham problem and are found to produce similar results for many-electron atoms.     

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.     

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.     

Structure and properties of metal-exchanged zeolites studied using gradient-corrected and hybrid functionals. II. Electronic structure and photoluminescence spectra

The influence of the choice of the exchange-correlation functional (semilocal gradient corrected or hybrid functionals) on the electronic properties of metal-exchanged zeolites has been investigated for Cu- and Co-exchanged chabazite. The admixture of exact exchange in hybrid functionals increases the fundamental gap of purely siliceous chabazite, leading to better agreement with experiment and many-body perturbation theory for close-packed SiO(2) polymorphs where detailed experimental information is available. For the metal-exchanged chabazite the increased exchange splitting strongly influences the position of the cation states relative to the framework bands-in general, gradient-corrected functionals locate the occupied cation states close to the valence-band maximum of the framework, while hybrid functionals shift the occupied cation states to larger binding energies and the empty states to higher energies within the fundamental gap. The photoluminescence spectra have been analyzed using fixed-moment total-energy calculations for excited spin states in structurally relaxed and frozen geometries. The geometrical relaxation of the excited states leads to large differences in excitation and emission energies which are more pronounced in calculations using hybrid functionals. Due to the stronger relaxation effects calculated with hybrid functionals, the large differences in the electronic spectra calculated with both types of functionals are not fully reflected in the photoluminescence spectra.     

Where is the spin? Understanding electronic structure and g-tensors for ruthenium complexes with redox-active quinonoid ligands

Understanding the bonding in transition metal complexes with redox-active ligands is a major challenge, for example in redox catalysis or in bioinorganic chemistry. In this work, electronic g-tensors, spin-density distributions, and electronic structure have been studied by different density functional methods for an extended series of complexes [Ru(acac)2(L)]n (n = -1, 0, +1; L = redox-active o-quinonoid ligand). Comparison is made with experimental g-tensors and g-tensor-based oxidation-state assignments for a number of experimentally studied examples, using both gradient-corrected (BP86) and hybrid functionals (B3LYP, BHLYP) representing a range of exact-exchange admixtures. Reasonable, albeit not perfect, agreement with experimental g-tensors is obtained in one-component DFT calculations with hybrid functionals. Analyses of spin densities confirm the assignment of the cationic complexes as predominantly d5-Ru(III) with a neutral quinonoid ligand. However, this conclusion is obtained only after inclusion of the appreciable spin polarization of the unrestricted determinant, while the singly occupied molecular orbital (SOMO) is localized more on the acac ligands. The anionic complexes turn out to be approximately halfway between a d6-Ru(II)/semiquinone and a d5-Ru(III)/catecholate formulation, but again only after taking into account the extensive spin polarization. Even the previous assignment of the neutral parent systems as d5-Ru(III)/semiquinone is not accurate, as a d6-Ru(II)/quinone resonance structure contributes to some extent. Very unusual trends in the spin contamination of the Kohn-Sham determinant with increasing exact-exchange admixture in some of the cationic complexes have been traced to an interplay between spin delocalization and spin polarization.     

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.     

Molecular dynamics simulation of liquid water: hybrid density functionals

The structure, dynamical, and electronic properties of liquid water utilizing different hybrid density functionals were tested within the plane wave framework of first-principles molecular dynamics simulations. The computational approach, which employs modified functionals with short-ranged Hartree-Fock exchange, was first tested in calculations of the structural and bonding properties of the water dimer and cyclic water trimer. Liquid water simulations were performed at the state point of 350 K at the experimental density. Simulations included three different hybrid functionals, a meta-functional, four gradient-corrected functionals, and the local density and Hartree-Fock approximations. It is found that hybrid functionals are superior in reproducing the experimental structure and dynamical properties as measured by the radial distribution function and self-diffusion constant when compared to the pure density functionals. The local density and Hartree-Fock approximations show strongly over- and understructured liquids, respectively. Hydrogen bond analysis shows that the hybrid functionals give slightly smaller average numbers of hydrogen bonds than pure density functionals but similar hydrogen bond populations. The average molecular dipole moments in the liquid from the three hybrid functionals are lower than those of the corresponding pure density functionals.     

The effective local potential method: implementation for molecules and relation to approximate optimized effective potential techniques

We have recently formulated a new approach, named the effective local potential (ELP) method, for calculating local exchange-correlation potentials for orbital-dependent functionals based on minimizing the variance of the difference between a given nonlocal potential and its desired local counterpart [V. N. Staroverov et al., J. Chem. Phys. 125, 081104 (2006)]. Here we show that under a mildly simplifying assumption of frozen molecular orbitals, the equation defining the ELP has a unique analytic solution which is identical with the expression arising in the localized Hartree-Fock (LHF) and common energy denominator approximations (CEDA) to the optimized effective potential. The ELP procedure differs from the CEDA and LHF in that it yields the target potential as an expansion in auxiliary basis functions. We report extensive calculations of atomic and molecular properties using the frozen-orbital ELP method and its iterative generalization to prove that ELP results agree with the corresponding LHF and CEDA values, as they should. Finally, we make the case for extending the iterative frozen-orbital ELP method to full orbital relaxation.     

Seminumerical calculation of the Hartree-Fock exchange matrix: application to two-component procedures and efficient evaluation of local hybrid density functionals

A two-component extension of the seminumerical procedure for the calculation of the Hartree-Fock (HF) exchange matrix recently presented by Neese et al. (Chem Phys 2009, 356, 98) was implemented into the program system TURBOMOLE. It is demonstrated that this allows for efficient self-consistent treatment of spin-orbit coupling at HF and hybrid density functional theory level. One-component HF calculations were performed to study the accuracy of integration grids and the exploitation of the molecular point group symmetry. The efficiency was tested, and for one-component calculations compared to the implementation realized by Neese. It was further demonstrated that local hybrid density functionals can be evaluated with this technique. The "prototype" of this class of functionals, Lh-BLYP, was applied to an organic molecule with more than 150 atoms.