Exchange and correlation energy in density functional theory

Several different versions of density functional theory (DFT) that satisfy Hohenberg–Kohn theorems are characterized by different definitions of a reference or model state determined by an N‐electron ground state. A common formalism is developed in which exact Kohn–Sham equations are derived for standard Kohn–Sham theory, for reference‐state density functional theory, and for unrestricted Hartree–Fock (UHF) theory considered as an exactly soluble model Hohenberg–Kohn theory. A natural definition of exchange and correlation energy functionals is shown to be valid for all such theories. An easily computed necessary condition for the locality of exchange and correlation potentials is derived. While it is shown that in the UHF model of DFT the optimized effective potential (OEP) exchange satisfies this condition by construction, the derivation shows that this condition is not, in general, sufficient to define an exact local exchange potential. It serves as a test to eliminate proposed local potentials that are not exact for ground states. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 77: 521–525, 2000     

Investigation of dephasing in an open quantum system under chaotic influence via a fractional Kohn–Sham scheme

In this work, the dynamics of dephasing (without relaxation) in the presence of a chaotic oscillator is theoretically investigated. The time‐dependent density functional theory framework was used in tandem with the Lindblad master equation approach for modeling the dissipative dynamics. Using the Kohn–Sham (K–S) scheme under certain approximations, the exact model for the potentials was acquired. In addition, a space‐fractional K–S scheme was developed (using the modified Riemann–Liouville operator) for modeling the dephasing phenomenon. Extensive analyses and comparative studies were then done on the results obtained using the space‐fractional K–S system and the conventional K–S system. © 2014 Wiley Periodicals, Inc.     

Analysis of self‐interaction correction for describing core excited states

Core‐excitation energies are calculated by the self‐interaction‐corrected time‐dependent density functional theory (SIC‐TDDFT) and SIC‐delta‐self‐consistent field (SIC‐ΔSCF) methods. For carbon monoxide, SIC‐TDDFT severely overestimates core‐excitation energies, while the SIC‐ΔSCF method using Kohn–Sham density functional theory (KS‐DFT) slightly overestimates. These behaviors are attributed to the fact that the self‐interaction errors in the total and orbital energies considerably differ. We evaluate the difference of the self‐interaction errors for the Slater exchange functional. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Examination of several density functionals in numerical Kohn–Sham calculations for atoms

We studied several exchange‐only and exchange–correlation energy density functionals in numerical, i.e., basis‐set‐free, nonrelativistic Kohn–Sham calculations for closed‐shell ¹S states of atoms and atomic ions with N electrons, where 2≤N≤120. Accurate total energies are presented to serve as reference data for algebraic approaches, as do the numerical Hartree–Fock results, which are also provided. Gradient‐corrected exchange‐only functionals considerably improve the total energies obtained from the usual local density approximation, when compared to the Hartree–Fock results. Such an improvement due to gradient corrections is not seen in general for highest orbital energies, neither for exchange‐only results (to be compared with Hartree–Fock results), nor for exchange–correlation results (to be compared with experimental ionization energies). © 2001 John Wiley & Sons, Inc. Int J Quant Chem 82: 227–241, 2001     

Orbital functional theory of linear response and excitation

Orbital functional theory (OFT) is based on a rule that determines a single‐determinant reference state Φ for any exact N‐electron eigenstate Ψ. An OFT model postulates an explicit correlation energy functional E_c of occupied orbital functions {?_i} and occupation numbers {n_i}. The orbital Euler–Lagrange equations are analogous to Kohn–Sham equations, but do not in general contain local potential functions. Time‐dependent Hartree–Fock theory is generalized in OFT to a formally exact linear response theory that includes electronic correlation. In the exchange‐only limit, the theory reduces to the random‐phase approximation of many‐body theory. The formalism determines excitation energies. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

Density functional theory for molecular and periodic systems using density fitting and continuous fast multipole method: Analytical gradients

A full implementation of analytical energy gradients for molecular and periodic systems is reported in the TURBOMOLE program package within the framework of Kohn–Sham density functional theory using Gaussian‐type orbitals as basis functions. Its key component is a combination of density fitting (DF) approximation and continuous fast multipole method (CFMM) that allows for an efficient calculation of the Coulomb energy gradient. For exchange‐correlation part the hierarchical numerical integration scheme (Burow and Sierka, Journal of Chemical Theory and Computation 2011, 7, 3097) is extended to energy gradients. Computational efficiency and asymptotic O(N) scaling behavior of the implementation is demonstrated for various molecular and periodic model systems, with the largest unit cell of hematite containing 640 atoms and 19,072 basis functions. The overall computational effort of energy gradient is comparable to that of the Kohn–Sham matrix formation. © 2016 Wiley Periodicals, Inc.     

Exact relations between the electron density and external potential for systems of interacting and noninteracting electrons

The differential virial theorem (DVT) is an explicit relation between the electron density ρ( r ), the external potential, kinetic energy density tensor, and (for interacting electrons) the pair function. The time‐dependent generalization of this relation also involves the paramagnetic current density. We present a detailed unified derivation of all known variants of the DVT starting from a modified equation of motion for the current density. To emphasize the practical significance of the theorem for noninteracting electrons, we cast it in a form best suited for recovering the Kohn–Sham effective potential v_s( r ) from a given electron density. The resulting expression contains only ρ( r ), v_s( r ), kinetic energy density, and a new orbital‐dependent ingredient containing only occupied Kohn–Sham orbitals. Other possible applications of the theorem are also briefly discussed. © 2012 Wiley Periodicals, Inc.     

Accurate Kohn–Sham ionization potentials from scaled‐opposite‐spin second‐order optimized effective potential methods

One important property of Kohn–Sham (KS) density functional theory is the exact equality of the energy of the highest occupied KS orbital (HOMO) with the negative ionization potential of the system. This exact feature is out of reach for standard density‐dependent semilocal functionals. Conversely, accurate results can be obtained using orbital‐dependent functionals in the optimized effective potential (OEP) approach. In this article, we investigate the performance, in this context, of some advanced OEP methods, with special emphasis on the recently proposed scaled‐opposite‐spin OEP functional. Moreover, we analyze the impact of the so‐called HOMO condition on the final quality of the HOMO energy. Results are compared to reference data obtained at the CCSD(T) level of theory. © 2016 Wiley Periodicals, Inc.     

Complexes of 4‐substituted phenolates with HF and HCN: Energy decomposition and electronic structure analyses of hydrogen bonding

We have computationally studied para‐X‐substituted phenols and phenolates (X = NO, NO₂, CHO, COMe, COOH, CONH₂, Cl, F, H, Me, OMe, and OH) and their hydrogen‐bonded complexes with B^? and HB (B = F and CN), respectively, at B3LYP/6‐311++G** and BLYP‐D/QZ4P levels of theory. Our purpose is to explore the structures and stabilities of these complexes. Moreover, to understand the emerging trends, we have analyzed the bonding mechanisms using the natural bond orbital scheme as well as Kohn–Sham molecular orbital (MO) theory in combination with quantitative energy decomposition analyses [energy decomposition analysis (EDA), extended transition state‐natural orbitals for chemical valence (ETS‐NOCV)]. These quantitative analyses allow for the construction of a simple physical model that explains all computational observations. © 2012 Wiley Periodicals, Inc.     

Density functional theory for molecular and periodic systems using density fitting and continuous fast multipole method: Stress tensor

A full implementation of the analytical stress tensor for periodic systems is reported in the TURBOMOLE program package within the framework of Kohn–Sham density functional theory using Gaussian-type orbitals as basis functions. It is the extension of the implementation of analytical energy gradients (Lazarski et al., Journal of Computational Chemistry 2016, 37, 2518–2526) to the stress tensor for the purpose of optimization of lattice vectors. Its key component is the efficient calculation of the Coulomb contribution by combining density fitting approximation and continuous fast multipole method. For the exchange-correlation (XC) part the hierarchical numerical integration scheme (Burow and Sierka, Journal of Chemical Theory and Computation 2011, 7, 3097–3104) is extended to XC weight derivatives and stress tensor. The computational efficiency and favorable scaling behavior of the stress tensor implementation are demonstrated for various model systems. The overall computational effort for energy gradient and stress tensor for the largest systems investigated is shown to be at most two and a half times the computational effort for the Kohn–Sham matrix formation. © 2019 Wiley Periodicals, Inc.     

Representation of Kohn–Sham free atom eigenfunctions by Slater‐type orbitals

Slater‐type orbitals are applied to represent the numerically obtained Kohn–Sham eigenfunction of free atom. The algorithm evaluating the nonlinear expansion coefficients of this approximation is described. Standard iterative solution of Kohn–Sham equation to obtain the nonlinear expansion coefficients is avoided and replaced by the projection method. First, the eigenfunction is obtained in the B‐spline space based on the Galerkin formulation of the finite element method. Then, based on the density functional theory, the conditions are formulated, which leads to the set of nonlinear equations. The proposed algorithm is general and can be applied for any atomic Kohn–Sham eigenfunction. As an examplary application of the proposed algorithm, the set of nonlinear equations is derived for occupied states of N, Al, Ga, and In atoms. The expansion coefficients, obtained for these atoms, are evaluated numerically by Newton procedure and listed in the tables. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

Energy‐surfaces from the upper bound of the Pauli kinetic energy

Based on the Kohn–Sham Pauli potential and the Kohn–Sham electron density, the upper bound of the Pauli kinetic energy is tested as a suitable replacement for the exact Pauli kinetic energy for application in orbital‐free density functional calculations. It is found that bond lengths for strong and moderately bound systems can be qualitatively predicted, but with a systematic shift toward larger bond distances with a relative error of 6% up to 30%. Angular dependence of the energy‐surface cannot be modeled with the proposed functional. Therefore, the upper bound model is the first parameter‐free functional expression for the kinetic energy that is able to qualitatively reproduce binding curves with respect to bond distortions. © 2016 Wiley Periodicals, Inc.     

About the compatibility between ansatzes and constraints for a local formulation of orbital‐free density functional theory

Functional properties that are exact for the Hohenberg–Kohn functional may turn into mutually exclusive constraints at a given level of ansatz. This is exemplarily shown for the local density approximation. Nevertheless, it is possible to reach exactly the Kohn–Sham data from an orbital‐free density functional framework based on simple one‐point functionals by starting from the Levy–Perdew–Sahni formulation. The energy value is obtained from the density‐potential pair, and therefore does not refer to the functional dependence of the potential expression. Consequently, the potential expression can be obtained from any suitable model and is not required to follow proper scaling behavior.     

Analytic derivatives for the XYG3 type of doubly hybrid density functionals: Theory,implementation, and assessment

We present a theoretical development of the equations required to perform an analytic geometry optimization of a molecular system using the XYG3 type of doubly hybrid (xDH) functionals. In contrast to the well‐established B2PLYP type of DH functionals, the energy expressions in the xDH functionals are constructed by using density and orbital information from another standard Kohn–Sham (KS) functional (e.g., B3LYP) for doing the self‐consistent field calculations. Thus, the xDH functionals are nonvariational in both the hybrid density functional part and the second‐order perturbation part, each of which requires formally to solve a coupled‐perturbed KS equation. An implementation is reported here which combines the two parts by defining a total Lagrangian such that only a single set of the Z‐vector equations need to be solved. The computational cost with our implementation is of the same order as those for the conventional Møller–Plesset theory to the second order (MP2) and B2PLYP. Systematic test calculations are provided for covalently bonded molecules as well as compounds involving the intramolecular nonbonded interactions for the main group elements. Satisfactory performance of the xDH functionals demonstrates that the extra computer time on top of the conventional KS procedure is well‐invested, in particular, when the standard KS functionals and MP2 as well, are problematic. © 2013 Wiley Periodicals, Inc.     

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.     

Charge‐Shift Bonding: A New and Unique Form of Bonding

Charge‐shift bonds (CSBs) constitute a new class of bonds different than covalent/polar‐covalent and ionic bonds. Bonding in CSBs does not arise from either the covalent or the ionic structures of the bond, but rather from the resonance interaction between the structures. This Essay describes the reasons why the CSB family was overlooked by valence‐bond pioneers and then demonstrates that the unique status of CSBs is not theory‐dependent. Thus, valence bond (VB), molecular orbital (MO), and energy decomposition analysis (EDA), as well as a variety of electron density theories all show the distinction of CSBs vis‐à‐vis covalent and ionic bonds. Furthermore, the covalent–ionic resonance energy can be quantified from experiment, and hence has the same essential status as resonance energies of organic molecules, e.g., benzene. The Essay ends by arguing that CSBs are a distinct family of bonding, with a potential to bring about a Renaissance in the mental map of the chemical bond, and to contribute to productive chemical diversity.     

The Stabilizing Effects in Gold Carbene Complexes

Bonding and stabilizing effects in gold carbene complexes are investigated by using Kohn–Sham density functional theory (DFT) and the intrinsic bond orbital (IBO) approach. The π‐stabilizing effects of organic substituents at the carbene carbon atom coordinated to the gold atom are evaluated for a series of recently isolated and characterized complexes, as well as intermediates of prototypical 1,6‐enyne cyclization reactions. The results indicate that these effects are of particular importance for gold complexes especially because of the low π‐backbonding contribution from the gold atom.     

New Insights into the Band‐Gap Narrowing of (N,P)‐Codoped TiO2 from Hybrid Density Functional Theory Calculations

The electronic properties of anatase‐TiO₂ codoped by N and P at different concentrations have been investigated via generalized Kohn–Sham theory with the Heyd–Scuseria–Ernzerhof (HSE06) hybrid functional for exchange‐correlation in the context of density functional theory. At high doping concentrations, we find that the high photocatalytic activity of (N, P)‐codoped anatase TiO₂ vis‐à‐vis the N‐monodoped case can be rationalized by a double‐hole‐mediated coupling mechanism [Yin et al., Phys. Rev. Lett. 2011, 106, 066801] via the formation of an effective N? P bond. On the other hand, Ti³⁺ and Ti⁴⁺ ions’ spin double‐exchange results in more substantial gap narrowing for larger separations between N and P atoms. At low doping concentrations, double‐hole‐coupling is dominant, regardless of the N? P distance.