Non-local relation between kinetic and exchange energy densities in Hartree–Fock theory

A non-local generalization K( r, r' ) of the kinetic energy t( r ) such that t( r ) = ∫K( r, r' ) dr' is defined using the idempotency property of the Hartree–Fock first-order density matrix. This is, in turn, related by means of an explicit differential equation to the non-local exchange energy density X( r, r' ). The relationship is illustrated for a couple of examples: with the Fermi-hole in a uniform electron gas, of importance in the local density version of density functional theory, and with inhomogeneous electron systems.     

Some questions on the exchange contribution to the effective potential of the Kohn–Sham theory

The current success of Density Functional Theory applications hinges upon the availability of explicitly density-dependent functionals to self-consistently solve a set of one-electron equations, the Kohn–Sham (KS) equations, which determine the occupied orbitals and its associated electronic density. In KS theory, a local exchange potential is proposed as part of an effective potential. This potential is compared to the exchange operator of the Hartree–Fock theory, which is of a non-local nature. The present paper discusses the variational framework of the KS equations, and the equivalence between both exchange potentials within a correlation-free theory. The common difficulties of current local exchange functionals to correctly simulate the non-locality of the exchange energy density in chemical systems are also analyzed and explained through an exactly solvable model. We give then numerical arguments and conclude by analyzing the performance of various commonly used approximations to exchange functionals.     

Semiempirical double-hybrid density functional with improved description of long-range correlation

The recently proposed new family of "double-hybrid" density functionals [Grimme, S. J. Chem. Phys. 2006, 124, 34108] replaces a fraction of the semi-local correlation energy by a non-local correlation energy expression that employs the Kohn-Sham orbitals in second-order many-body perturbation theory. These functionals have provided results of high accuracy over a wide range of properties but fail to accurately describe long-range van der Waals interactions. In this work, a distance-dependent scaling factor for the non-local correlation energy is introduced to address this problem, and two new double-hybrid density functionals are proposed. The new functionals are optimized with the finite cc-pVTZ basis on training sets of atomization energies and intermolecular interaction energies. They are compared against (scaled) second-order M?ller-Plesset perturbation theories and popular density functionals including the hybrid-GGA functional B3-LYP and the first double-hybrid functional (B2-PLYP). Tests are performed on an extensive set including reaction energies, barrier heights, weakly interacting complexes, transition-metal systems, molecular geometries, and harmonic vibrational frequencies. Within the cc-pVTZ atomic orbital basis, we have demonstrated the ability to find a parametrization scheme which is simultaneously able to describe thermochemistry and weakly bound systems with a satisfactory degree of accuracy.     

Successful a priori modeling of CO adsorption on Pt(111) using periodic hybrid density functional theory

The adsorption of CO on the surface of metals such as Pt(111) is of great interest owing to the industrial importance of the catalytic oxidation of pollutant CO. To date, reliable high-level calculations of this process have not been possible, a situation often referred to as the "CO/Pt(111) puzzle". Standard generalized-gradient-approximation density functional theory approaches fail to capture key details of the binding, such as the location of the adsorption site, while cluster approaches using alternative methods show some but insufficient improvement. Using a new computational methodology combining hybrid density functionals containing non-local Hartree-Fock exchange with periodic imaging plane-wave-based techniques, we demonstrate that key aspects of the adsorption of CO on Pt(111), including the identification of the absorption site and CO frequency change, can now be adequately modeled. The binding is dominated by both CO dative covalent bonding and metal-to-molecule pi back-bonding, effects requiring realistic alignment of both the molecular HOMO and LUMO orbitals with respect to the metal Fermi energy.     

An improved model potential for liquid lithium

An improvement in the computed structure of liquid lithium was obtained by adding the Xα Slater exchange energy to the electron pseudopotential of Ashcroft, Singwi, Sjölander et al. The remaining discrepancies can only be corrected by using a fully non-local pseudopotential for lithium, and perhaps by improving experimental accuracies of structure factors.     

Density approximation to the average Hartree-Fock exchange potential for atoms

A non-local density-based approximation to the average Hartree-Fock (HF) exchange potential is developed. The new potential is formulated within the spin-dependent version of the weighted density approximation, and is based on a novel form of the (exchangeonly) pair-correlation function for electrons in finite systems. The results for total energies and one-electron orbital energies of atoms are reasonably accurate in comparison with those obtained using the exact average HF potential, or the exact orbital-dependent HF potential.     

The importance of electronic exchange-correlation in non-self-consistent frozen density approximation

The effects of different treatments for the exchange-correlation energy on the accuracy of non-self-consistent frozen density approximation (FDA) are discussed. Local spin density approximation (LSDA) and non-local spin density approximation (NLSDA) are employed, respectively. Corresponding results obtained by using full-self-consistent density functional theory (DFT) are also given for the purpose of comparison. Explicit calculations for hydrogen bonds, covalent bonds and ionic bonds indicate that, comparing with LSDA, NLSDA can improve the accuracy of FDA as well as that of DFT. This improvement attributed to the refinements in the treatment for the electronic exchange-correlation energy may help to extend the application of FDA.     

Towards chemical accuracy for the thermodynamics of large molecules: new hybrid density functionals including non-local correlation effects

Two hybrid density functionals that include a second-order perturbation correction for non-local correlation effects are tested for the full G3/05 test set. Very large AO basis sets including core-polarization/correlation functions have been employed that yield for the first time results quite close to the basis set limit for this set. The B2-PLYP functional and the new mPW2-PLYP approach with a modified exchange part give by far the lowest MAD over the whole G3/05 set ever reported for a DFT method (2.5 and 2.1 kcal mol(-1), respectively). The big improvement compared to common density functionals is further demonstrated by the reduction of the maximum and minimum errors (outliers) and by much smaller errors for complicated molecular systems.     

Design of an exchange functional with correct asymptotics

Knowledge of asymptotic conditions on exchange allows for a better design of exchange energy expressions in density functional theory. By working inside an exchange-only framework, the fulfillment of such conditions by some of the most widely used exchange functionals is discussed. In turn, we propose a model expression which partially meets the energetics and asymptotics of both the exchange energy density and potential. Improvement upon the local spin density approximation without the use of generalized gradient corrections is also presented. Hartree-Fock orbitals are employed to build electron densities. © 1997 John Wiley & Sons, Inc.     

Local and non-local properties controlling the strength of Brønsted acids and bases

A general formula giving both the energy of heterolytic dissociation of Brønsted acids and the proton affinity of Brønsted bases is proposed, thus providing a rationale for the proton exchange between acids and bases. The formula improves and extends previous ones proposed in the past (DeKock and Barbachyn, Cerofolini and Re, this collaboration) and establishes a correlation between, on one side, the proton affinity of the (neutral or negatively charged) base and, on the other side, the local properties of the proton-hosting site and non-local properties of the base, like its size and ionization energy. Our model involves intuitively appealing chemical concepts and can be very useful to close the gap between the experimental view and the ab initio calculations approach.     

A simplified implementation of van der Waals density functionals for first-principles molecular dynamics applications

We present a simplified implementation of the non-local van der Waals correlation functional introduced by Dion et al. [Phys. Rev. Lett. 92, 246401 (2004)] and reformulated by Román-Pérez et al. [Phys. Rev. Lett. 103, 096102 (2009)]. The proposed numerical approach removes the logarithmic singularity of the kernel function. Complete expressions of the self-consistent correlation potential and of the stress tensor are given. Combined with various choices of exchange functionals, five versions of van der Waals density functionals are implemented. Applications to the computation of the interaction energy of the benzene-water complex and to the computation of the equilibrium cell parameters of the benzene crystal are presented. As an example of crystal structure calculation involving a mixture of hydrogen bonding and dispersion interactions, we compute the equilibrium structure of two polymorphs of aspirin (2-acetoxybenzoic acid, C(9)H(8)O(4)) in the P2(1)/c monoclinic structure.     

Phenomenological chemical reactivity theory for mobile electrons

We present herein a model to deal with the chemical reactivity, selectivity and site activation concepts of π electron systems derived by merging the classical Coulson–Longuet-Higgins response function theory based on the Hückel molecular orbital (HMO) theory and the conceptual density functional theory. HMO-like expressions for the electronic chemical potential, chemical hardness and softness, including their local counterparts, atomic and bond Fukui functions and non-local response functions are derived. It is shown that sophisticated non-local concepts as site activation may be cast into deeper physical grounds by introducing a simplified version of static response functions. In this way, useful quantities such as self and mutual polarizabilities originally defined through the HMO parameters can be redefined as self and mutual softnesses. The model is illustrated by discussing the classical Hammett free energy relationship describing inductive substituent effects on the reactivity of benzoic acids.     

Does the gradient-regulated connection improve the description of correlated metal bond properties?

Gradient-regulated connection (GRAC) is a generalized gradient approximation exchange density functional designed by combining the revPBE and PW91 exchange functionals to impose their behaviors in the slowly- and fast-varying density regions, respectively. Such a construction allows one single density functional to accurately estimate both covalent and weak interactions occurring in main-group-based molecular systems. For the first time, the assessment of the performance of the GRAC exchange functional is extended to the modeling of various metal bond energy and structure properties. This assessment shows that when GRAC is coupled with the Perdew, Burke, Ernzerhof (PBE) correlation, the resulting exchange-correlation density functional is an excellent alternative to global hybrids to model bond dissociation energy, atomic electronic excitation energy, and bond length structure properties of single-reference metal bonds. It also shows that coupling with the Tognetti, Cortona, Adamo (TCA) correlation constitutes a robust approach to tackle energy bond properties of organometallic complexes with multi-reference character.     

Benchmark results for empirical post-GGA functionals: difficult exchange problems and independent tests

Many of the most promising new density functionals have improved the treatment of non-local exchange effects with the help of semi-empirical information and more sophisticated recipes for combining Hartree-Fock and local exchange approximations. In order to quantify recent advancements and identify directions for improvement, we have examined a broad spectrum of test problems. We evaluate the performance of several new hybrid density functionals (ωB97, ωB97X, ωB97X-D, LRC-ωPBEh, M06, M06-2X, and M06-HF) on a variety of chemical problems, some sensitive to the treatment of exact exchange (which we have hoped to systematically improve) and some which require a balanced treatment of correlation. Since all of the functionals under consideration are parameterized with ground-state thermochemical data, the benchmark aims to determine the applicability of the new density functionals to cases that have not been considered in the optimization of the semi-empirical parameters. The first class of benchmarks includes the excitation energies of 21 molecules (83 states) primarily from a recent benchmark conducted by Tozer and co-workers, with some additional references from data made available from the groups of Thiel and Truhlar. We briefly examine the conformational preferences of a small peptide and complete our study with two recently published sets of data that have shown large, systematic errors in simple alkane thermochemistry. While our results indicate that the more general hybrids currently under development perform well for problems outside of their parameterization and improve over the standard hybrid density functionals in an essentially systematic way, there is still a significant self-interaction error in the more difficult cases. Functionals based on a range-separation of exchange and functionals depending on the kinetic-energy density both perform comparably, and there is evidence for complementary strengths.     

Hybrid density functional theory for pi-stacking interactions: application to benzenes, pyridines, and DNA bases

The suitability of a hybrid density functional to qualitatively reproduce geometric and energetic details of parallel pi-stacked aromatic complexes is presented. The hybrid functional includes an ad hoc mixture of half the exact (HF) exchange with half of the uniform electron gas exchange, plus Lee, Yang, and Parr's expression for correlation energy. This functional, in combination with polarized, diffuse basis sets, gives a binding energy for the parallel-displaced benzene dimer in good agreement with the best available high-level calculations reported in the literature, and qualitatively reproduces the local MP2 potential energy surface of the parallel-displaced benzene dimer. This method was further critically compared to high-level calculations recently reported in the literature for a range of pi-stacked complexes, including monosubstituted benzene-benzene dimers, along with DNA and RNA bases, and generally agrees with MP2 and/or CCSD(T) results to within +/-2 kJ mol(-1). We also show that the resulting BH&H binding energy is closely related to the electron density in the intermolecular region. The net result is that the BH&H functional, presumably due to fortuitous cancellation of errors, provides a pragmatic, computationally efficient quantum mechanical tool for the study of large pi-stacked systems such as DNA.     

元素电负性和硬度的密度泛函理论研究

应用密度泛函理论的DFT LDA、DFT LDA/NL和改进的Slater过渡态方法,把元素的电离能和电子亲合能的计算扩展到周期表的103种元素.并用有限差分方法计算了这103种元素的电负性和硬度.计算中考虑了相对论效应.计算结果比以前Robles等用密度泛函理论的XGL和Xα近似的交换相关泛函的计算结果有所改进,更接近实验值.     

Regional self-interaction correction of density functional theory

We propose a new simple scheme for self-interaction correction (SIC) of exchange functionals in the density functional theory. In the new scheme, exchange energies are corrected by substituting exchange self-interactions for exchange functionals in regions of self-interaction. To classify the regions of self-interaction, we take advantage of the property of the total kinetic energy density approaching the Weizs?cker density in the case of electrons in isolated orbitals. The scheme differs from conventional SIC methods in that it produces optimized molecular structures. Applying the scheme to the calculation of reaction energy barriers showed that it provides a clear improvement in cases where the barriers are underestimated by conventional "pure" functionals. In particular, we found that this scheme even reproduces a transition state that is not given by pure functionals.     

Hartree–Fock exchange energy of an inhomogeneous electron gas

We examine the short-range behavior of the spherically averaged Hartree–Fock exchange charge density by performing a simple Taylor expansion. On the basis of this expansion, a theoretical model is constructed that generates gradient correction terms to the local density approximation for the exchange energy of an inhomogeneous electron gas. In particular, we derive the Xαβ exchange energy functional and a theoretical value for the parameter β. Our value for β agrees well with previous empirical estimates, and with empirical calculations in the present work.     

A generalized gradient approximation for exchange derived from the model potential of van Leeuwen and Baerends

The common way to obtain energies from Kohn-Sham exchange potentials is by using the Levy-Perdew virial relation. For potentials that are not functional derivatives (i.e., nearly all model exchange potentials in existence), this approach leads to energy expressions that lack translational and rotational invariance. We propose a method for constructing potential-based energy functionals that are free from these artifacts. It relies on the same line-integration technique that gives rise to the Levy-Perdew relation, but uses density scaling instead of coordinate scaling. The method is applicable to any exchange or correlation potential that depends on the density explicitly, and correctly recovers the parent energy functional from a functional derivative. To illustrate our approach we develop a properly invariant generalized gradient approximation for exchange starting from the model potential of van Leeuwen and Baerends.     

Development of density functionals for thermochemical kinetics

A density functional theory exchange-correlation functional for the exploration of reaction mechanisms is proposed. This functional, denoted BMK (Boese-Martin for Kinetics), has an accuracy in the 2 kcal/mol range for transition state barriers but, unlike previous attempts at such a functional, this improved accuracy does not come at the expense of equilibrium properties. This makes it a general-purpose functional whose domain of applicability has been extended to transition states, rather than a specialized functional for kinetics. The improvement in BMK rests on the inclusion of the kinetic energy density together with a large value of the exact exchange mixing coefficient. For this functional, the kinetic energy density appears to correct "back" the excess exact exchange mixing for ground-state properties, possibly simulating variable exchange.