Thomas–Fermi theory from a perturbative treatment of the hydrogenic limit energy density functional: A possible link to local density functional theory

A perturbation expansion which connects the hydrogenic limit energy density functional to the Thomas–Fermi functional is discussed. This perturbation series, where the Coulomb energy density functional is treated as the perturbation to the hydrogenic limit functional, is, in fact, the q = (N/Z) expansion of Thomas–Fermi theory. A truncated form of the first-order correction to the functional provides further insight into the model which treats the ground state energy as a local functional of the electron density.     

A density matrix variational calculation for atomic Be

The ground-state energy of the beryllium atom is calculated using a variational procedure in which the elements of the two-body reduced density matrix (particle–particle matrix) are the variational parameters. It is shown that, for this problem and with the limited number of spin-orbitals used, the trace condition and the simultaneous nonnegativity conditions on the particle–particle, the particle–hole, and the hole–hole matrices form a complete solution to the N-representability problem. The energy obtained is – 14.61425 a.u., practically identical to the value given by a configuration interaction calculation which uses the same states. The effects of weakening the nonnegativity conditions on each of the matrices in turn were also explored.     

One‐particle density matrix functional for correlation in molecular systems

Based on the analysis of the general properties for the one‐ and two‐particle reduced density matrices, a new natural orbital functional is obtained. It is shown that by partitioning the two‐particle reduced density matrix in an antisymmeterized product of one‐particle reduced density matrices and a correction Γ^c we can derive a corrected Hartree–Fock theory. The spin structure of the correction term from the improved Bardeen–Cooper–Schrieffer theory is considered to take into account the correlation between pairs of electrons with antiparallel spins. The analysis affords a nonidempotent condition for the one‐particle reduced density matrix. Test calculations of the correlation energy and the dipole moment of several molecules in the ground state demonstrate the reliability of the formalism. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem 94: 317–323, 2003     

Direct minimization of the energy in density functional theory

The minimization of the energy functional of the first-order density matrix γ( r , r ') is achieved using unitary transformations applied to γ. Equivalently, such transformations can be carried out also on one-electron orbitals (natural orbitals) and their occupation (integer or non-integer) numbers. The conventional local density approximation based on the electron density p( r ) is then considered as a special case. The direct minimization of the energy functional of p with respect to the parameters of the unitary transformation leads to stationary conditions that are all equivalent to the Kohn–Sham equations. Preliminary numerical tests show that the proposed algorithms for the direct minimization of the energy work in a satisfactory manner. © John Wiley & Sons, Inc.     

Kinetic energy density in terms of electron density for closed-shell atoms in a bare Coulomb field

A long-term aim in density functional theory is to obtain the kinetic energy density t(r) in terms of the ground-state electron density ρ(r). Here, t(r) is written explicitly in terms of ρ(r) for an arbitrary number 𝒩 of closed shells in a bare Coulomb field. In the limit as 𝒩→∞, closed results for t(r) follow from the earlier analysis of ρ(r) by Heilmann and Lieb. [Phys. Rev. A 52 , 3628 (1995)]. © 1998 John Wiley & Sons, Inc. Int J Quant Chem 66 : 281–283, 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.     

Structure of the ammonia dimer studied by density functional theory

Self-consistent Kohn–Sham density functional calculations have been carried out to study the structure of the ammonia dimer. The local-density approximation yields unusually large binding energy and short internitrogen distance compared with the experimental and more accurate theoretical data. The results from the Becke–Perdew gradient-corrected functionals are generally in good agreement with those at the SCF MP 2 level when the geometry is fully optimized with various large basis sets. With our best estimation, the staggered quasi-linear structure (C_s) is 0.6 kcal/mol lower in energy than the symmetric cyclic one (C_2h). The hydrogen-bonded N—H bond in the staggered quasi-linear structure is found to be 0.008 Å longer than the N—H bond in ammonia. In our calculations, we could not find the minima on the energy surface corresponding to the two asymmetric cyclic structures suggested by microwave spectra and coupled pair functional calculations. © 1994 John Wiley & Sons, Inc.     

Simple tests for density functional methods

The performance of the currently used generalized gradient approximation density functionals is analyzed using several simple, yet critical requirements. We analyze the effects of the self-interaction error, the inclusion of the exact exchange, and the parameter settings used in the popular three-parameter hybrid density functionals. The results show that the elimination of the self-interaction error from the current density functionals lead to very poor results for H₂. The inclusion of the exact exchange does not significantly influence the self-interaction corrected results. The variation of the A, B, and C parameters of a hybrid DFT method influences the H(SINGLE BOND)H equilibrium bond length through a very simple linear equation, and it is possible to reproduce the experimental H(SINGLE BOND)H distance with appropriate selection of these parameters, although an infinite number of solutions exists. Similar results were obtained for the total energy and the electron density along the internuclear axis. The analysis of the exact KS potential at the bond critical point of the dissociating H₂ molecule shows that, for this property, the second order Moller–Plesset perturbation theory yields a better potential than the density functionals studied in this article. © 1997 John Wiley & Sons, Inc. J Comput Chem 18 : 1534–1545, 1997     

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     

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.     

Correlation of electron density and bond length to band gap for binary oxides and halides

The best quantity correlated to the electronic energy band gap is found for alkali and alkaline-earth metal oxides and halides with face centered cubic (fcc) structure based on density functional theory and Bader's atom-in-molecule theory. Previous studies show the correlation of the band gap to the ground state electron density at the bond critical point (BCP). Whereas, in quantum mechanics, the gap between the energy levels of one dimensional square well potential is inversely proportional to the square of the width of the well which is the metal–nonmetal chemical bond length in our case. These motivate the proposition of a new quantity Q, the ratio of the density at the BCP to the square of the bond length. Our study reveals that, for the aforementioned materials, the band gap has a strong correlation to Q when they are multiplied by the density at the BCP.     

Energy band structure of (SN)x chain: Unrestricted Hartree–Fock and charge density wave solutions

The ab initio Hartree–Fock crystal orbital method is used for the calculation of the energy band structure of a one-dimensional model of (SN)_x. Two energy band structures are described corresponding to the self-consistent spin density wave (SDW ) and the self-consistent charge density wave (CDW ) solution, respectively.     

Performance of parallel TURBOMOLE for density functional calculations

The parallelization of density functional treatments of molecular electronic energy and first-order gradients is described, and the performance is documented. The quadrature required for exchange correlation terms and the treatment of exact Coulomb interaction scales virtually linearly up to 100 nodes. The RI-J technique to approximate Coulomb interactions (by means of an auxiliary basis set approximation for the electron density) even shows superlinear speedup on distributed memory architectures. The bottleneck is then linear algebra. Demonstrative application examples include molecules with up to 300 atoms and 3000 basis functions that can now be treated in a few hours per geometry optimization cycle in C₁ symmetry. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 1746–1757, 1998     

On the connections between Brueckner–coupled-cluster,density-dependent Hartree–Fock,and density functional theory

The existence of an effective one-particle Hamiltonian in the Brueckner coupled cluster model naturally leads to the definition of an effective interaction G, which is a function of the T₂ amplitudes. Two types of approximations to G are proposed: One is purely phenomenological, while the other is based on approximations to the Brueckner T₂ equation. In both cases, the resulting effective interaction may be viewed as electron-density-dependent. Generalizing Hartree–Fock theory to accommodate density-dependent interactions (DDHF ), a method is obtained that is capable of accounting for correlation effects in an independent particle framework. The heuristic Skyrme force, successfully used in nuclear physics to model nucleon–nucleon interactions, is presented here as an example of an effective electron–electron correlation interaction. Due to the δ-function character of the Skyrme force, it is possible to express the energy in this model by an integral over an energy density, thus formally providing a connection between DDHF and density functional theory for this particular case. An approximation to the Brueckner T₂ equation is also proposed in the coordinate representation. In this model, the density-matrix dependence of T₂ is reduced to a nonlocal electron density dependence by means of an expansion which introduces terms that depend on the gradient of the density. The first term in this expansion amounts to a “local density approximation” to Brueckner coupled cluster theory. © 1995 John Wiley & Sons, Inc.     

Improved local density functional approach for atomic systems

Through a new local density approximation to the kinetic energy density functional introduced by us recently, a simple Thomas–Fermi-like scheme for the direct calculation of electron density in atoms is proposed. The calculated density is nonsingular at the nucleus and the energy values are in very good agreement with the corresponding Hartree–Fock results for atoms. © 1994 John Wiley & Sons, Inc.     

Small-angle x-ray scattering study of density fluctuation in pressure-densified polystyrene glasses

Densified polystyrene glasses, prepared by cooling from the liquid state under elevated pressure, were studied by small-angle x-ray scattering at ambient pressure. The density fluctuation, determined from the x-ray data, showed a decrease with increasing pressure up to about 1.5 kbar, and then leveled off to a fairly constant value. The reduction in the density fluctuation produced by the pressure is much greater than the associated decrease in the specific volume. The observed change in density fluctuation is consistent with the view that the density fluctuation in glassy polymers consists of dynamic and quasistatic components and that the first of these can be correlated with the compressibility of the glass. The present data on the density fluctuation, in conjunction with the available data on volume and enthalpy, can be interpreted to mean that in pressure-densified glasses unfavorable chain configurations are trapped in local energy minima, and the strain energy thus stored can promote segmental motion leading to volume expansion at temperatures far below T_g. Some preliminary evidence indicating the formation of microcavities in these pressure-densified glasses is also presented.     

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     

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     

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.