First‐order correlation‐kinetic contribution to Kohn–Sham exchange charge density function in atoms,using quantal density functional theory approach

Using the static exchange‐correlation charge density concept, the total integrated exchange‐charge density function is calculated within the nonrelativistic spin‐restricted exchange‐only (i) optimized effective potential model, and (ii) nonvariational local potential derived from the exchange‐only work potential within the quantal density functional theory, for the ground‐state isoelectronic series: Ga⁺, Zn, Cu^?; In⁺, Cd, Ag^?; and Tl⁺, Hg, Au^?. The difference between the exchange charge density function derived from these potentials is employed to evaluate the first‐order correlation‐kinetic contribution to the integrated exchange charge density. This contribution is found to be important for both the intra‐ and inter‐shell regions. Screening effects on the contribution due to the nd¹⁰ (n = 3–5) subshells are discussed through comparisons with similar calculations on Ca, Sr, and Ba, wherein nd¹⁰ electrons are absent. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

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     

Electron correlations in Kohn–Sham exchange‐only theory

We provide an interpretation for the “exchange” energy and potential of Kohn–Sham exchange‐only theory, or equivalently that of the optimized potential method (OPM), which shows that in addition to contribution due to the Pauli exclusion principle, there is a kinetic component to these properties. The interpretation is in terms of a conservative field R ^OPM( r ), which is a sum of two fields, one representative of Pauli electron correlations and the other of kinetic effects. The OPM exchange potential is derived via the differential virial theorem to be the work done to move an electron in the field R ^OPM( r ). The OPM exchange energy is then expressed via the integral virial theorem in terms of this field. A similar interpretation for the energy and potential may also be derived directly from the OPM integral equation. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 71:473–480, 1999     

Formulation of Strutinsky's method for atomic systems in the extended Kohn–Sham scheme

Strutinsky's standard averaging method is formulated in the framework of the extended Kohn–Sham scheme and a two‐step procedure permitting the application of the method is proposed. A Taylor‐series expansion of the ground‐state energy‐function of the occupation numbers is derived, which involves the averaged energy as the leading term and shell corrections as smaller terms. Numerical applications for atoms and ions from Be through Ar are presented and discussed. © 2002 Wiley Periodicals, Inc. Int J Quantum Chem, 2002     

Fitting elephants in the density functionals zoo: Statistical criteria for the evaluation of density functional theory methods as a suitable replacement for counting parameters

Counting parameters has become customary in the density functional theory community as a way to infer the transferability of popular approximations to the exchange‐correlation functionals. Recent work in data science, however, has demonstrated that the number of parameters of a fitted model is not related to the complexity of the model itself, nor to its eventual overfitting. Using similar arguments, here, we show that it is possible to represent every modern exchange‐correlation functional approximations using just one single parameter. This procedure proves the futility of the number of parameters as a measure of transferability. To counteract this shortcoming, we introduce and analyze the performance of three statistical criteria for the evaluation of the transferability of exchange‐correlation functionals. The three criteria are called Akaike information criterion, Vapnik‐Chervonenkis criterion, and cross‐validation criterion and are used in a preliminary assessment to rank 60 exchange‐correlation functional approximations using the ASCDB database of chemical data.     

Kinetic energy density for orbital‐free density functional calculations by axiomatic approach

An axiomatic approach is herein used to determine the physically acceptable forms for general D‐dimensional kinetic energy density functionals (KEDF). The resulted expansion captures most of the known forms of one‐point KEDFs. By statistically training the KEDF forms on a model problem of noninteracting kinetic energy in 1D (six terms only), the mean relative accuracy for 1000 randomly generated potentials is found to be better than the standard KEDF by several orders of magnitudes. The accuracy improves with the number of occupied states and was found to be better than for a system with four occupied states. Furthermore, we show that free fitting of the coefficients associated with known KEDFs approaches the exactly analytic values. The presented approach can open a new route to search for physically acceptable kinetic energy density functionals and provide an essential step toward more accurate large‐scale orbital free density functional theory calculations.     

Numerical integration of the Kohn–Sham equations: Integral method for bound states

A numerical method is presented that solves the multicenter Kohn–Sham equations. The method couples the resolution of the integral form of the equation at a given energy with an iterative search for the eigenvalues. The validity of the method is checked by comparing some test calculations for diatomics with results in the literature from other numerical methods. For these calculations the wave functions are expanded in partial waves either on one center or on two centers with the help of the partitioning of space in fuzzy cells. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 74: 49–54, 1999     

Self-consistent field calculations using two-body density functionals for correlation energy component: I. Atomic systems

Self-consistent field calculations are done using two-body density functionals for the correlation energy. The corresponding functional derivatives are obtained and used in pseudo-eigenvalue equations analogous to the Kohn–Sham ones. The examples studied include atomic systems from He to Ar. The values obtained for ionization potentials, electron affinities, dipole polarizabilities, and virial ratios from these calculations are given, and the effect of exchange is addressed. The results obtained are in good agreement with experimental values, and are of the same quality as those given by accurate exchange-correlation functionals. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 1887–1898, 1998     

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.     

Role of electron–electron coalescence density in density functional theory

An expression for the evaluation of electron–electron coalescence density as a functional of the density for any electron system is proposed. The formula, clarifies previously advanced upper bounds for this quantity and provides a method to independently estimate the system‐averaged on‐top exchange–correlation hole. The relationship with the on‐top pair density shows that producing the true electron–electron coalescense should be considered as a leading physical requirement for trial wave functions in any energy minimization scheme. © 2002 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

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.     

Self-consistent field calculations using two-body density functionals for correlation energy component: II. Small molecules

In part I of this series, self-consistent calculations using two-body density functionals for correlation energy were done and applied to atomic systems, giving very good results. We now apply the same scheme to small molecules. The examples studied include diatomic (H₂, Li₂, B₂, C₂, N₂, O₂, F₂, HLi, HBe, HB, HF, and HCl) as well as polyatomic (H₂O, NH₃, H₂O₂, and O₃) molecules at their ground states. The values reported for equilibrium geometries, atomization energies, vibrational frequencies, and dipole moments are compared with experimental and other theoretical calculations, with good agreement in most cases. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 1899–1908, 1998     

Standard grids for high‐precision integration of modern density functionals: SG‐2 and SG‐3

Density‐functional approximations developed in the past decade necessitate the use of quadrature grids that are far more dense than those required to integrate older generations of functionals. This category of difficult‐to‐integrate functionals includes meta‐generalized gradient approximations, which depend on orbital gradients and/or the Laplacian of the density, as well as functionals based on B97 and the popular “Minnesota” class of functionals, each of which contain complicated and/or oscillatory expressions for the exchange inhomogeneity factor. Following a strategy introduced previously by Gill and co‐workers to develop the relatively sparse “SG‐0” and “SG‐1” standard quadrature grids, we introduce two higher‐quality grids that we designate SG‐2 and SG‐3, obtained by systematically “pruning” medium‐ and high‐quality atom‐centered grids. The pruning procedure affords computational speedups approaching a factor of two for hybrid functionals applied to systems of atoms, without significant loss of accuracy. The grid dependence of several popular density functionals is characterized for various properties. © 2017 Wiley Periodicals, Inc.     

Novel method of self‐interaction corrections in density functional calculations

It is demonstrated that the commonly applied self‐interaction correction (SIC) used in density functional theory does not remove all self‐interaction. We present as an alternative a novel method that, by construction, is totally free from self‐interaction. The method has the correct asymptotic 1/r dependence. We apply the new theory to localized f electrons in praseodymium and compare with the old version of SIC, the local density approximation (LDA) and with an atomic Hartree–Fock calculation. The results show a lowering of the f level, a contraction of the f electron cloud and a lowering of the total energy by 13 eV per 4 f electron compared to LDA. The equilibrium volume of the new SIC method is close to the ones given by LDA and the older SIC method and is in good agreement with experiment. The experimental cohesive energy is in better agreement using the new SIC method, both compared to LDA and another SIC method. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 81: 247–252, 2001     

Coupled‐cluster reaction barriers of : An application of the coupled‐cluster//Kohn–Sham density functional theory model chemistry

In this work, we report a theoretical investigation concerning the use of the popular coupled‐cluster//Kohn‐Sham density functional theory (CC//KS‐DFT) model chemistry, here applied to study the entrance channel of the reaction, namely by comparing CC//KS‐DFT calculations with KS‐DFT, MRPT2//CASSCF, and CC//CASSCF results from our previous investigations. This was done by performing single point energy calculations employing several coupled cluster methods and using KS‐DFT geometries optimized with six different functionals, while conducting a detailed analysis of the barrier heights and topological features of the curves and surfaces here obtained. The quality of this model chemistry is critically discussed in the context of the title reaction and also in a wider context. © 2013 Wiley Periodicals, Inc.     

Correlation energy,correlated electron density,and exchange‐correlation potential in some spherically confined atoms

We report correlation energies, electron densities, and exchange‐correlation potentials obtained from configuration interaction and density functional calculations on spherically confined He, Be, Be²⁺, and Ne atoms. The variation of the correlation energy with the confinement radius R_c is relatively small for the He, Be²⁺, and Ne systems. Curiously, the Lee–Yang–Parr (LYP) functional works well for weak confinements but fails completely for small R_c. However, in the neutral beryllium atom the CI correlation energy increases markedly with decreasing R_c. This effect is less pronounced at the density‐functional theory level. The LYP functional performs very well for the unconfined Be atom, but fails badly for small R_c. The standard exchange‐correlation potentials exhibit significant deviation from the “exact” potential obtained by inversion of Kohn–Sham equation. The LYP correlation potential behaves erratically at strong confinements. © 2016 Wiley Periodicals, Inc.     

Hierarchy of equations in the generalized density functional theory

Functional relations and equations of hierarchy in the generalized density functional theory (DFT) are derived from coordinate scaling and adiabatic connection. Local and nonlocal solutions for the noninteracting kinetic energy, exchange energy, correlation energy, and the kinetic energy correction functionals are presented. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006