Exact Exchange Hole Research Articles

Solid-state performance of a meta-GGA screened hybrid density functional constructed from Pauli kinetic enhancement factor dependent semilocal exchange hole.

The semilocal form of an exchange hole is highly useful in developing non-local range-separated hybrid density functionals for finite and extended systems. The way to construct the conventional exact exchange hole model is based on either the Taylor series expansion or the reverse engineering technique from the corresponding exchange energy functional. Although the latter technique is quite popular in the context of generalized gradient approximation (GGA) functionals, the same for the meta-GGA functionals is not so much explored. Thus, in this study, we propose a reverse-engineered semilocal exchange hole of a meta-GGA functional, which only depends on the meta-GGA ingredient α (also known as the Pauli kinetic energy enhancement factor). The model is subsequently used to design the short-range-separated meta-GGA hybrid density functional. We show that the present method can be successfully applied for several challenging problems in the context of solids, especially for which the GGA based hybrid fails drastically. This assessment proves that the present functional is quite useful for materials sciences. Finally, we also use this method for several molecular test cases, where the results are also as comparative as its base semilocal functional.

Crystal Structure Evolution of Fluorine under High Pressure

Fluorinated compounds in the last decade were applied as photo-thermo-refractive glasses, high-stress lubricants, and pharmaceutical drugs due to their good mechanical properties and biocompatibility. Although fluorinated materials are largely employed, the possibility of predicting new structures was limited by the impossibility to use density functional theory (DFT) to describe interatomic and intermolecular interactions correctly. This is seen linearly to increase with fluorine concentration. In crystal structure prediction, modern algorithms are usually combined with first-principles methods employed for global solution, which sometimes fail to describe systems as in the case of strongly correlated materials. Fluorine is one of the tricky elements, which is characterized by relativistic effects and no overlap between the DFT exchange hole and the exact exchange hole. Thus, no relativistic exchange–correlation functional was seen to adequately describe fluorine. In this work, we have found an excellent compromise to investigate fluorinated materials using a combination of SCAN (exchange) and rVV10 (correlation) functionals. This was found through the fundamental study of α- and β-fluorine phases, showing α-fluorine as the most stable structure at temperatures lower than 35 K and 0 GPa with respect to β-fluorine. Further, we have computed crystal structure evolution under pressure looking for new stable fluorine allotropes using the USPEX evolutionary algorithm coupled with the SCAN-rVV10 exchange–correlation functional discovering two phase transitions: one from C2/c (i.e., α-fluorine) to Cmca at ∼5.5 GPa and from Cmca to the P4̅21c phase at 220 GPa; all these structures possess metallic behavior. The achievements of this work lie far beyond the thermodynamic of fluorine crystals under pressure. It will give the right instrument to understand the chemical behavior of fluorinated materials under pressure with consequent great speed up to the crystal structure prediction of fluorinated and fluorine-based materials.

Model DFT exchange holes and the exact exchange hole: Similarities and differences.

Exchange and correlation holes are unique quantum concepts for understanding the nature of electron interactions based on quantum conditional probabilities. Among those, the exact exchange hole is of special interest since it is derived rigorously from first principles without approximations and is often modeled by approximate exchange expressions of density functional theory. In this work, the algorithm for the computation of the spherically averaged exact exchange hole for a given reference point is developed and implemented for molecular orbitals in Gaussian basis functions. The formulas include a novel recursive relation for the spherical average of the Bessel function of the first kind and the asymptotic expressions when the exponential factor of the Bessel function becomes large. This new capability is used to explore the extent to which current popular model exchange holes resemble or differ from the exact exchange hole. Point-wise accuracy of the exchange holes for isolated atoms is important in local hybrid schemes, real-space models of static correlation, and others. We find in this vein that among the models tested here, only the BR89 exchange hole seems more or less suitable for that purpose, while better approximations are still very much on demand. Analyzing the deviations of model exchange holes from the exact exchange hole in molecules such as H2 and Cr2 upon bond stretching reveals new aspects of the left-right static correlation.

Communication: Density functional theory model for multi-reference systems based on the exact-exchange hole normalization.

The correct description of multi-reference electronic ground states within Kohn-Sham density functional theory (DFT) requires an ensemble-state representation, employing fractionally occupied orbitals. However, the use of fractional orbital occupation leads to non-normalized exact-exchange holes, resulting in large fractional-spin errors for conventional approximative density functionals. In this communication, we present a simple approach to directly include the exact-exchange-hole normalization into DFT. Compared to conventional functionals, our model strongly improves the description for multi-reference systems, while preserving the accuracy in the single-reference case. We analyze the performance of our proposed method at the example of spin-averaged atoms and spin-restricted bond dissociation energy surfaces.

Testing Exact Upper Bounds to Exact Exchange.

The exact exchange energy and its energy density are useful but computationally expensive ingredients in density functional approximations for Kohn-Sham density functional theory. We present detailed tests of some exact nonempirical upper bounds to exact exchange. These "Rung 3.5" upper bounds contract the Kohn-Sham one-particle density matrix with model density matrices used to construct semilocal model exchange holes and invoke the Cauchy-Schwarz inequality. The contraction automatically eliminates the computationally expensive long-range component of the exact exchange hole. Numerical tests show that the exchange upper bounds underestimate total exchange energies while predicting other properties with accuracy approaching standard hybrid approximations.

Semilocal exchange hole with an application to range-separated density functionals

Exchange-correlation hole is a central concept in density functional theory. It not only provides justification for an exchange-correlation energy functional, but also serves as a local ingredient in nonlocal range-separation density functional. However, due to the nonlocal nature, modelig the conventional exact exchange hole presents a great challenge to density functional theory. In this work, we propose a semilocal exchange hole underlying the Tao-Perdew-Staroverov-Scuseria (TPSS) meta-GGA functional. The present model is distinct from previous models at small separation between an electron and the hole around the electron. It is also different in the way it interpolates between the rapidly varying iso-orbital density and the slowly varying density, which is determined by the wave vector analysis based on the exactly solvable infinite barrier model for jellium surface. Our numerical tests show that the exchange hole generated from this model mimics the conventional exact exchange hole quite well for atoms. Finally, as a simple application, we apply the hole model to construct a TPSS-based range-separation functional. Our tests show that this TPSS-based range-separation functional can substantially improve TPSS band gaps and barrier heights, without losing much accuracy of molecular atomization energies.

Communication: Near-locality of exchange and correlation density functionals for 1- and 2-electron systems.

The uniform electron gas and the hydrogen atom play fundamental roles in condensed matter physics and quantum chemistry. The former has an infinite number of electrons uniformly distributed over the neutralizing positively charged background, and the latter only one electron bound to the proton. The uniform electron gas was used to derive the local spin density approximation to the exchange-correlation functional that undergirds the development of the Kohn-Sham density functional theory. We show here that the ground-state exchange-correlation energies of the hydrogen atom and many other 1- and 2-electron systems are modeled surprisingly well by a different local spin density approximation (LSDA0). LSDA0 is constructed to satisfy exact constraints but agrees surprisingly well with the exact results for a uniform two-electron density in a finite, curved three-dimensional space. We also apply LSDA0 to excited or noded 1-electron densities, where it works less well. Furthermore, we show that the localization of the exact exchange hole for a 1- or 2-electron ground state can be measured by the ratio of the exact exchange energy to its optimal lower bound.

How close are the Slater and Becke-Roussel potentials in solids?

The Becke–Roussel (BR) potential [Phys. Rev. A1989, 39, 3761] was proposed as an approximation to the Slater potential, which is the Coulomb potential generated by the exact exchange hole. In the present work, a detailed comparison between the Slater and BR potentials in solids is presented. It is shown that the two potentials usually lead to very similar results for the electronic structure; however, in a few cases, e.g., Si, Ge, or strongly correlated systems like NiO, the fundamental band gap or magnetic properties can differ markedly. Such differences should not be neglected when the computationally expensive Slater potential is replaced by the cheap semilocal BR potential in approximations to the exact-exchange Kohn–Sham potential, such as the one proposed by Becke and Johnson [J. Chem. Phys.2006, 124, 221101].

Halogen Bonding from Dispersion-Corrected Density-Functional Theory: The Role of Delocalization Error.

Halogen bonds are formed when a Lewis base interacts with a halogen atom in a different molecule, which acts as an electron acceptor. Due to its charge transfer component, halogen bonding is difficult to model using many common density-functional approximations because they spuriously overstabilize halogen-bonded dimers. It has been suggested that dispersion-corrected density functionals are inadequate to describe halogen bonding. In this work, we show that the exchange-hole dipole moment (XDM) dispersion correction coupled with functionals that minimize delocalization error (for instance, BH&HLYP, but also other half-and-half functionals) accurately model halogen-bonded interactions, with average errors similar to other noncovalent dimers with less charge-transfer effects. The performance of XDM is evaluated for three previously proposed benchmarks (XB18 and XB51 by Kozuch and Martin, and the set proposed by Bauzá et al.) spanning a range of binding energies up to ∼50 kcal/mol. The good performance of BH&HLYP-XDM is comparable to M06-2X, and extends to the "extreme" cases in the Bauzá set. This set contains anionic electron donors where charge transfer occurs even at infinite separation, as well as other charge transfer dimers belonging to the pnictogen and chalcogen bonding classes. We also show that functional delocalization error results in an overly delocalized electron density and exact-exchange hole. We propose intermolecular Bader delocalization indices as an indicator of both the donor-acceptor character of an intermolecular interaction and the delocalization error coming from the underlying functional.

Approximating the exchange energy through the nonempirical exchange-factor approach

Various widely used approaches for the exchange-correlation energy of Kohn-Sham theory, such as the generalized gradient approximation and hybrid methods, suffer from the one-electron self-interaction error. This error manifests itself in one-electron or two-electron spin-unpolarized systems. We propose new approximations to the exchange energy that are exact for the mentioned problems. Starting from the weighted density approximation (WDA), we introduce several weighting factors (that we refer to as exchange factors) which are designed to reproduce various properties of the exact exchange hole. Furthermore, a new exact condition for the exchange factor is presented that concerns the asymptotic behavior of the exchange hole as one electron coordinate moves out to infinity. Observance of this and other conditions results in exchange functionals that improve upon the WDA.

Counterintuitive Coulomb hole around the bond midplane

The Coulomb hole does not have its largest depth around an electron in or near the bond midplane. It splits into two parts, localized on both nuclear sites forming the bond. Even counterintuitive positive values of the "hole" around such a position may be observed. This happens when the Fermi hole is deeper than the total exchange-correlation hole at the reference electron position. This Coulomb "heap" is shown to arise from correlation effects on the one-electron density rather than correlation effects in the pair density. Left-right correlation tends to enhance the effect of the nuclear attraction, contracting the electron density around the nuclear positions and depleting the bond center region. Possible alternative definitions of the Coulomb hole are discussed, including one based on the exact Kohn-Sham exchange hole. Approximate density functional theory methods (generalized gradient approximation) are not accurate enough to realize the advantages of this definition.

Position-dependent exact-exchange energy for slabs and semi-infinite jellium

The position-dependent exact-exchange energy per particle $\varepsilon_x(z)$ (defined as the interaction between a given electron at $z$ and its exact-exchange hole) at metal surfaces is investigated, by using either jellium slabs or the semi-infinite (SI) jellium model. For jellium slabs, we prove analytically and numerically that in the vacuum region far away from the surface $\varepsilon_{x}^{\text{Slab}}(z \to \infty) \to - e^{2}/2z$, {\it independent} of the bulk electron density, which is exactly half the corresponding exact-exchange potential $V_{x}(z \to \infty) \to - e^2/z$ [Phys. Rev. Lett. {\bf 97}, 026802 (2006)] of density-functional theory, as occurs in the case of finite systems. The fitting of $\varepsilon_{x}^{\text{Slab}}(z)$ to a physically motivated image-like expression is feasible, but the resulting location of the image plane shows strong finite-size oscillations every time a slab discrete energy level becomes occupied. For a semi-infinite jellium, the asymptotic behavior of $\varepsilon_{x}^{\text{SI}}(z)$ is somehow different. As in the case of jellium slabs $\varepsilon_{x}^{\text{SI}}(z \to \infty)$ has an image-like behavior of the form $\propto - e^2/z$, but now with a density-dependent coefficient that in general differs from the slab universal coefficient 1/2. Our numerical estimates for this coefficient agree with two previous analytical estimates for the same. For an arbitrary finite thickness of a jellium slab, we find that the asymptotic limits of $\varepsilon_{x}^{\text{Slab}}(z)$ and $\varepsilon_{x}^{\text{SI}}(z)$ only coincide in the low-density limit ($r_s \to \infty$), where the density-dependent coefficient of the semi-infinite jellium approaches the slab {\it universal} coefficient 1/2.

The RPA Atomization Energy Puzzle.

There is current interest in the random phase approximation (RPA), a "fifth-rung" density functional for the exchange-correlation energy. RPA has full exact exchange and constructs the correlation with the help of the unoccupied Kohn-Sham orbitals. In many cases (uniform electron gas, jellium surface, and free atom), the correction to RPA is a short-ranged effect that is captured by a local spin density approximation (LSDA) or a generalized gradient approximation (GGA). Nonempirical density functionals for the correction to RPA were constructed earlier at the LSDA and GGA levels (RPA+), but they are constructed here at the fully nonlocal level (RPA++), using the van der Waals density functional (vdW-DF) of Langreth, Lundqvist, and collaborators. While they make important and helpful corrections to RPA total and ionization energies of free atoms, they correct the RPA atomization energies of molecules by only about 1 kcal/mol. Thus, it is puzzling that RPA atomization energies are, on average, about 10 kcal/mol lower than those of accurate values from experiment. We find here that a hybrid of 50% Perdew-Burke-Ernzerhof GGA with 50% RPA+ yields atomization energies much more accurate than either one does alone. This suggests a solution to the puzzle: While the proper correction to RPA is short-ranged in some systems, its contribution to the correlation hole can spread out in a molecule with multiple atomic centers, canceling part of the spread of the exact exchange hole (more so than in RPA or RPA+), making the true exchange-correlation hole more localized than in RPA or RPA+. This effect is not captured even by the vdW-DF nonlocality, but it requires the different kind of full nonlocality present in a hybrid functional.

Local hybrid functionals based on density matrix products.

We present a novel similarity metric comparing exact and semilocal density functional theory (DFT) exchange holes in real space. This metric is obtained from the product of the one-particle density matrix and the uniform electron gas model density matrix. The metric is bound between 0 and 1, 1 in the uniform electron gas, 0 in regions asymptotically far from finite systems, and can detect delocalization of the exact exchange hole and effective fractional occupations. We also present a parameter-free local hybrid functional that uses this similarity metric to locally mix exact and semilocal DFT exchange energy densities. The resulting functional gives better thermochemistry and reaction barrier heights than our original local hybrids [Jaramillo et al., J. Chem. Phys. 118, 1068 (2003)], while retaining moderate accuracy for symmetric radical cation dimers.

Improved electronic structure and optical properties of sp-hybridized semiconductors using LDA+U SIC

We propose the local density approximation (LDA) plus an on-site Coulomb self-interaction-like correction (SIC) potential for describing sp-hybridized bonds in semiconductors and insulators. We motivate the present LDA+U SIC scheme by comparing the exact exchange (EXX) hole with the LDA exchange hole. The LDA+U SIC method yields good band-gap energies Eg and dielectric constants epsilon(omega ~ 0) of Si, Ge, GaAs, and ZnSe. We also show that LDA consistently underestimates the gamma-point effective electron m c and light-hole m lh masses, and the underlying reason for this is a too strong light-hole-electron coupling within LDA. The advantages of the LDA+U SIC approach are a computational time of the same order as the ordinary LDA, the orbital dependent LDA+U SIC exchange-correlation interaction is asymmetric analogously to the EXX potential, and the method can be used for materials and compounds involving localized d- and f-orbitals.

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.

Why the generalized gradient approximation works and how to go beyond it

The local spin density LSD approximation, while of only moderate accuracy, has proven extremely reliable over three decades of use. We argue that any gradient-corrected functional should preserve the correct features of LSD, even if the system under study . contains no regions of small density gradient. The Perdew)Wang 1991 PW91 functional . respects this condition, while, e.g., the Lee)Yang)Parr LYP correlation functional violates it. We extend this idea to the next generation of density functionals, those which . incorporate exact exchange via the optimized effective potential OEP , with a model in which the correlation hole is constructed from the exact exchange hole. The resulting exchange-correlation hole is deeper and less diffuse than the exact exchange hole. We denote such a functional as ''locally correlated Hartree)Fock'' and list a variety of conditions such a functional should satisfy. We demonstrate the promise of this approach with a crude simple model. Q 1997 John Wiley & Sons, Inc.