Kohn-Sham Kinetic Energy Research Articles

Reworking the Tao-Mo exchange-correlation functional. I. Reconsideration and simplification.

The revised, regularized Tao-Mo (rregTM) exchange-correlation density functional approximation (DFA) [A. Patra, S. Jana, and P. Samal, J. Chem. Phys. 153, 184112 (2020) and Jana et al., J. Chem. Phys. 155, 024103 (2021)] resolves the order-of-limits problem in the original TM formulation while preserving its valuable essential behaviors. Those include performance on standard thermochemistry and solid data sets that is competitive with that of the most widely explored meta-generalized-gradient-approximation DFAs (SCAN and r2SCAN) while also providing superior performance on elemental solid magnetization. Puzzlingly however, rregTM proved to be intractable for de-orbitalization via the approach of Mejía-Rodríguez and Trickey [Phys. Rev. A 96, 052512 (2017)]. We report investigation that leads to diagnosis of how the regularization in rregTM of the z indicator functions (z = the ratio of the von-Weizsäcker and Kohn-Sham kinetic energy densities) leads to non-physical behavior. We propose a simpler regularization that eliminates those oddities and that can be calibrated to reproduce the good error patterns of rregTM. We denote this version as simplified, regularized Tao-Mo, sregTM. We also show that it is unnecessary to use rregTM correlation with sregTM exchange: Perdew-Burke-Ernzerhof correlation is sufficient. The subsequent paper shows how sregTM enables some progress on de-orbitalization.

Self-consistent implementation of locally scaled self-interaction-correction method.

Recently proposed local self-interaction correction (LSIC) method [Zope et al., J. Chem. Phys. 151, 214108 (2019)] is a one-electron self-interaction-correction (SIC) method that uses an iso-orbital indicator to apply the SIC at each point in space by scaling the exchange-correlation and Coulomb energy densities. The LSIC method is exact for the one-electron densities, also recovers the uniform electron gas limit of the uncorrected density functional approximation, and reduces to the well-known Perdew-Zunger SIC (PZSIC) method as a special case. This article presents the self-consistent implementation of the LSIC method using the ratio of Weizsäcker and Kohn-Sham kinetic energy densities as an iso-orbital indicator. The atomic forces as well as the forces on the Fermi-Löwdin orbitals are also implemented for the LSIC energy functional. Results show that LSIC with the simplest local spin density functional predicts atomization energies of the AE6 dataset better than some of the most widely used generalized-gradient-approximation (GGA) functional [e.g., Perdew-Burke-Ernzerhof (PBE)] and barrier heights of the BH6 database better than some of the most widely used hybrid functionals (e.g., PBE0 and B3LYP). The LSIC method [a mean absolute error (MAE) of 0.008 Å] predicts bond lengths of a small set of molecules better than the PZSIC-LSDA (MAE 0.042 Å) and LSDA (0.011 Å). This work shows that accurate results can be obtained from the simplest density functional by removing the self-interaction-errors using an appropriately designed SIC method.

Restoring the iso-orbital limit of the kinetic energy density in relativistic density functional theory.

In contrast to nonrelativistic density functional theory, the ratio between the von Weizsäcker and the Kohn-Sham kinetic energy density, commonly used as iso-orbital indicator t within exchange-correlation functionals beyond the generalized-gradient level, violates the exact iso-orbital limit and the appropriate parameter range, 0 ≤ t ≤ 1, in relativistic density functional theory. Based on the exact decoupling procedure within the infinite-order two-component method and the Cauchy-Schwarz inequality, we present corrections to the relativistic and the picture-change-transformed nonrelativistic kinetic energy density that restores these exact constraints. We discuss the origin of the new correction terms and illustrate the effectiveness of the current approach for several representative cases. The proposed generalized iso-orbital indicator tλ is expected to be a useful ingredient for the development of relativistic exchange-correlation functionals.

On the many-electron self-interaction error of the semilocal exchange hole based meta-GGA level range-separated hybrid with the B88 hybrids

Delocalization error or many electron self-interaction error related to the fractional occupation number is studied within the framework of the density matrix expansion based meta-GGA level range-separated hybrid functional. The present study mainly focuses on the performance of that functional with the BECKE88 exchange based local, global and range-separated hybrids. It is observed that using the Kohn–Sham kinetic energy dependent exchange hole, several properties relating to the fractional occupation number are significantly improved without hindering the atomization energy accuracy. The present systematic assessment and functional performance for G2/148 test set and (non-) hydrogen barrier heights indicate that beyond the BECKE88 based range-separated functional, present meta-GGA level range-separated functional can be used for further molecular simulations.

Semi-local machine-learned kinetic energy density functional with third-order gradients of electron density.

A semi-local kinetic energy density functional (KEDF) was constructed based on machine learning (ML). The present scheme adopts electron densities and their gradients up to third-order as the explanatory variables for ML and the Kohn-Sham (KS) kinetic energy density as the response variable in atoms and molecules. Numerical assessments of the present scheme were performed in atomic and molecular systems, including first- and second-period elements. The results of 37 conventional KEDFs with explicit formulae were also compared with those of the ML KEDF with an implicit formula. The inclusion of the higher order gradients reduces the deviation of the total kinetic energies from the KS calculations in a stepwise manner. Furthermore, our scheme with the third-order gradient resulted in the closest kinetic energies to the KS calculations out of the presented functionals.

Exploration of near the origin and the asymptotic behaviors of the Kohn-Sham kinetic energy density for two-dimensional quantum dot systems with parabolic confinement.

The behaviors of the positive definite Kohn-Sham kinetic energy density near the origin and at the asymptotic region play a major role in designing meta-generalized gradient approximations (meta-GGAs) for exchange in low-dimensional quantum systems. It is shown that near the origin of the parabolic quantum dot, the Kohn-Sham kinetic energy differs from its von Weizsäcker counterpart due to the p orbital contributions, whereas in the asymptotic region, the difference between the above two kinetic energy densities goes as ∼ρ(r)r2. All these behaviors have been explored using the two-dimensional isotropic quantum harmonic oscillator as a test case. Several meta-GGA ingredients are then studied by making use of the above findings. Also, the asymptotic conditions for the exchange energy density and the potential at the meta-GGA level are proposed using the corresponding behaviors of the two kinetic energy densities.

Kinetic‐energy‐density dependent semilocal exchange‐correlation functionals

AbstractWe present the theory of semilocal exchange‐correlation (XC) energy functionals which depend on the Kohn–Sham kinetic energy density (KED), including the relevant class of meta‐generalized gradient approximation (meta‐GGA) functionals. Thanks to the KED ingredient, meta‐GGA functionals can satisfy different exact constraints for XC energy and can be made one‐electron self‐correlation free. This leads to a better accuracy over a wider range of properties with respect to GGAs, often reaching the accuracy of hybrid functionals, but at much reduced computational cost. An extensive survey of the relevant literature on existing KED dependent XC functionals is provided, considering nonempirical, semi‐empirical, and fully empirical ones. A deeper analysis and a wide benchmark are presented for functionals derived considering only exact constraints and parameters obtained from model and/or atomic systems.

Visualization and analysis of the Kohn-Sham kinetic energy density and its orbital-free description in molecules.

We visualize the Kohn-Sham kinetic energy density (KED) and the ingredients--the electron density, its gradient, and Laplacian--used to construct orbital-free models of it, for the AE6 test set of molecules. These are compared to related quantities used in metaGGA's, to characterize two important limits--the gradient expansion and the localized-electron limit typified by the covalent bond. We find the second-order gradient expansion of the KED to be a surprisingly successful predictor of the exact KED, particularly at low densities where this approximation fails for exchange. This contradicts the conjointness conjecture that the optimal enhancement factors for orbital-free kinetic and exchange energy functionals are closely similar in form. In addition we find significant problems with a recent metaGGA-level orbital-free KED, especially for regions of strong electron localization. We define an orbital-free description of electron localization and a revised metaGGA that improves upon atomization energies significantly.

Kinetic Energy of Hydrocarbons as a Function of Electron Density and Convolutional Neural Networks.

We demonstrate a convolutional neural network trained to reproduce the Kohn-Sham kinetic energy of hydrocarbons from an input electron density. The output of the network is used as a nonlocal correction to conventional local and semilocal kinetic functionals. We show that this approximation qualitatively reproduces Kohn-Sham potential energy surfaces when used with conventional exchange correlation functionals. The density which minimizes the total energy given by the functional is examined in detail. We identify several avenues to improve on this exploratory work, by reducing numerical noise and changing the structure of our functional. Finally we examine the features in the density learned by the neural network to anticipate the prospects of generalizing these models.

Subsystem density functional theory with meta-generalized gradient approximation exchange-correlation functionals.

We analyze the methodology and the performance of subsystem density functional theory (DFT) with meta-generalized gradient approximation (meta-GGA) exchange-correlation functionals for non-bonded molecular systems. Meta-GGA functionals depend on the Kohn-Sham kinetic energy density (KED), which is not known as an explicit functional of the density. Therefore, they cannot be directly applied in subsystem DFT calculations. We propose a Laplacian-level approximation to the KED which overcomes this limitation and provides a simple and accurate way to apply meta-GGA exchange-correlation functionals in subsystem DFT calculations. The so obtained density and energy errors, with respect to the corresponding supermolecular calculations, are comparable with conventional approaches, depending almost exclusively on the approximations in the non-additive kinetic embedding term. An embedding energy error decomposition explains the accuracy of our method.

Kohn-Sham kinetic energy density in the nuclear and asymptotic regions: Deviations from the von Weizsäcker behavior and applications to density functionals

We show that the Kohn-Sham positive-definite kinetic energy (KE) density significantly differs from the von Weizs\"acker (VW) one at the nuclear cusp as well as in the asymptotic region. At the nuclear cusp, the VW functional is shown to be linear and the contribution of p-type orbitals to the KE density is theoretically derived and numerically demonstrated in the limit of infinite nuclear charge, as well in the semiclassical limit of neutral large atoms. In the latter case, it reaches 12 of the KE density. In the asymptotic region we find new exact constraints for meta Generalized Gradient Approximation (meta-GGA) exchange functionals: with an exchange enhancement factor proportional to $\sqrt{\alpha}$, where $\alpha$ is the common meta-GGA ingredient, both the exchange energy density and the potential are proportional to the exact ones. In addition, this describes exactly the large-gradient limit of quasi-two dimensional systems.

Single-point kinetic energy density functionals: A pointwise kinetic energy density analysis and numerical convergence investigation

We present a comprehensive study of single-point kinetic energy density functionals (KEDFs) to be used in orbital-free density functional theory (DFT) calculations. We first propose a form of KEDFs based on a pointwise Kohn-Sham (KS) kinetic energy density (KED) and electron localization function (ELF) analysis. We find that the ELF and modified enhancement factor have a very strong and transferable correlation with the reduced density in various bulk metals. The non-self-consistent kinetic energy errors predicted by our KEDF models are decreased greatly compared to previously reported generalized gradient approximation (GGA) KEDFs. Second, we perform self-consistent calculations with various single-point KEDFs and investigate their numerical convergence behavior. We find striking numerical instabilities for previous GGA KEDFs; most of the GGA KEDFs fail to converge and show unphysical densities during the optimization. In contrast, our KEDFs demonstrate stable convergence, and their self-consistent results of various bulk properties agree reasonably well with KSDFT. A further detailed KED analysis reveals an interesting bifurcation phenomenon in defective metals and alloys, which may shed light on directions for future KEDF development.

One-electron self-interaction and the asymptotics of the Kohn–Sham potential: an impaired relation

One-electron self-interaction and an incorrect asymptotic behavior of the Kohn–Sham exchange–correlation potential are among the most prominent limitations of many present-day density functionals. However, a one-electron self-interaction-free energy does not necessarily lead to the correct long-range potential. This is shown here explicitly for local hybrid functionals. Furthermore, carefully studying the ratio of the von Weizsäcker kinetic energy density to the (positive) Kohn–Sham kinetic energy density, τW/τ, reveals that this ratio, which frequently serves as an iso-orbital indicator and is used to eliminate one-electron self-interaction effects in meta-generalized-gradient approximations and local hybrid functionals, can fail to approach its expected value in the vicinity of orbital nodal planes. This perspective article suggests that the nature and consequences of one-electron self-interaction and some of the strategies for its correction need to be reconsidered.

ACKS2: Atom-condensed Kohn-Sham DFT approximated to second order

A new polarizable force field (PFF), namely atom-condensed Kohn-Sham density functional theory approximated to second order (ACKS2), is proposed for the efficient computation of atomic charges and linear response properties of extended molecular systems. It is derived from Kohn-Sham density functional theory (KS-DFT), making use of two novel ingredients in the context of PFFs: (i) constrained atomic populations and (ii) the Legendre transform of the Kohn-Sham kinetic energy. ACKS2 is essentially an extension of the Electronegativity Equalization Method (EEM) [W. J. Mortier, S. K. Ghosh, and S. Shankar, J. Am. Chem. Soc. 108, 4315 (1986)] in which two major EEM shortcomings are fixed: ACKS2 predicts a linear size-dependence of the dipole polarizability in the macroscopic limit and correctly describes the charge distribution when a molecule dissociates. All ACKS2 parameters are defined as atoms-in-molecules expectation values. The implementation of ACKS2 is very similar to that of EEM, with only a small increase in computational cost.

Conditions on the Kohn–Sham kinetic energy and associated density

AbstractEffective, explicitly density‐dependent (i.e., orbital‐free) approximations to the Kohn–Sham kinetic energy functional Ts remain elusive. In the course of developing such functionals in collaboration with Frank Harris, various issues have arisen that have gone unaddressed. We consider four of them here: positivity of the KS kinetic energy density, supposed requirements on its functional dependence, the use of cutoffs to insure positivity, and the role of the Coulomb virial theorem in the asymptotic analysis of Ts. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2009

Kinetic energy density functionals from the Airy gas with an application to the atomization kinetic energies of molecules

We construct and study several semilocal density-functional approximations for the positive Kohn-Sham kinetic energy density. These functionals fit the kinetic energy density of the Airy gas and they can be accurate for integrated kinetic energies of atoms, molecules, jellium clusters, and jellium surfaces. We find that these functionals are the most accurate ones for atomization kinetic energies of molecules and for fragmentation of jellium clusters. We also report that local and semilocal kinetic energy functionals can show ``binding'' when the density of a spin-unrestricted Kohn-Sham calculation is used.

Condition on the Kohn–Sham kinetic energy and modern parametrization of the Thomas–Fermi density

We study the asymptotic expansion of the neutral-atom energy as the atomic number Z-->infinity, presenting a new method to extract the coefficients from oscillating numerical data. Recovery of the correct expansion yields a condition on the Kohn-Sham kinetic energy that is important for the accuracy of approximate kinetic energy functionals for atoms, molecules, and solids. For example, this determines the small gradient limit of any generalized gradient approximation and conflicts somewhat with the standard gradient expansion. Tests are performed on atoms, molecules, and jellium clusters using densities constructed from Kohn-Sham orbitals. We also give a modern, highly accurate parametrization of the Thomas-Fermi density of neutral atoms.

On the local representation of the electronic momentum operator in atomic systems

The local quantum theory is applied to the study of the momentum operator in atomic systems. Consequently, a quantum-based local momentum expression in terms of the single-electron density is determined. The limiting values of this function correctly obey two fundamental theorems: Kato's cusp condition and the Hoffmann-Ostenhof and Hoffmann-Ostenhof exponential decay. The local momentum also depicts the electron shell structure in atoms as given by its local maxima and inflection points. The integration of the electron density in a shell gives electron populations that are in agreement with the ones expected from the Periodic Table of the elements. The shell structure obtained is in agreement with the higher level of theory computations, which include the Kohn-Sham kinetic energy density. The average of the local kinetic energy associated with the local momentum is the Weizsacker kinetic energy. In conclusion, the local representation of the momentum operator provides relevant information about the electronic properties of the atom at any distance from the nucleus.

The Kohn-Sham kinetic energy density as indicator of the electron localization: Atomic shell structure

In this report, it is shown that the Kohn-Sham (KS) kinetic energy density (KED) contains the average local electrostatic potential (ALEP) and the average local ionization energy (ALIE); the shell structure in atomic systems is presented as one application of the KS-KED. By writing the KS-KED from the KS equations, this quantity was divided in three contributions: orbital, Coulomb, and exchange correlation. By studying several closed and open shell atoms, the shell structure was established by the maxima presented by the Coulomb contribution and the minima in the orbital contribution of the KS-KED. The exchange-correlation contribution to the KS-KED does not show maxima or minima, but this quantity shows bumps where the division between shells is expected. The results obtained in this work were compared with other shell structure indicators such as the electron localization function, the ALEP, the ALIE, and the radial distribution function. The most important result in this work is related to the fact that even when the ALEP and the ALIE functions were built with different arguments to each other, they are contained in the KS-KED. In this way, the KS-KED shows its importance to reveal the electron localization in atomic systems.

Using the Kohn–Sham formalism in pair density-functional theories

There has been some recent interest in ‘generalized’ density functional theories, in which the fundamental variable is the density of electron pairs (pair density), instead of the electron density. Much of the work along these lines has focused on Weisacker-type functionals of the density and has not explicitly considered the N-representability problem. As is discussed here, such theories cannot be accurate. Some new functionals for the kinetic energy, which may address the problem, are proposed, and the correlation component is discussed as a relatively small correction to the Kohn–Sham kinetic energy.