The Liu‐Parr power series expansion of the Pauli kinetic energy functional with the incorporation of shell‐inducing traits: Atoms

An approximate expression for the Pauli kinetic energy functional T_p is advanced in terms of the Liu‐Parr expansion [S. Liu, R.G. Parr, Phys. Rev. A 1997 , 55, 1792] which involves a power series of the one‐electron density. We use this explicit functional for T_p to compute the value of the noninteracting kinetic energy functional T_s of 34 atoms, from Li to Kr (and their positive and negative monoions). In particular, we examine the effect that a shell‐by‐shell mean‐square optimization of the expansion coefficients has on the kinetic energy values and explore the effect that the size of the expansion, given by the parameter n, has on the accuracy of the approximation. The results yield a mean absolute percent error for 34 neutral atoms of 0.15, 0.08, 0.04, 0.03, and 0.01 for expansions with n = 3, 4, 5, 6, and 7, respectively (where ). We show that these results, which are the most accurate ones obtained to date for the representation of the noninteracting kinetic energy functional, stem from the imposition of shell‐inducing traits. We also compare these Liu‐Parr functionals with the exact but nonexplicit functional generated in the local‐scaling transformation version of DFT.     

Static and Frequency‐Dependent Dipole–Dipole Polarizabilities of All Closed‐Shell Atoms up to Radium: A Four‐Component Relativistic DFT Study

We test the performance of four‐component relativistic density functional theory by calculating the static and frequency‐dependent electric dipole–dipole polarizabilities of all (ground‐state) closed‐shell atoms up to Ra. We consider 12 nonrelativistic functionals, including three asymptotically shape‐corrected functionals, by using two smooth interpolation schemes introduced by the Baerends group: the gradient‐regulated asymptotic connection (GRAC) procedure and the statistical averaging of (model) orbital potentials (SAOP). Basis sets of doubly augmented triple‐zeta quality are used. The results are compared to experimental data or to accurate ab initio results. The reference static electric dipole polarizability of palladium has been obtained by finite‐field calculations using the coupled‐cluster singles, doubles, and perturbative triples method within this work. The best overall performance is obtained using hybrid functionals and their GRAC shape‐corrected versions. The performance of SAOP is among the best for nonhybrid functionals for Group 18 atoms but its precision degrades when considering the full set of atoms. In general, we find that conclusions based on results obtained for the rare‐gas atoms are not necessarily representative of the complete set of atoms. GRAC cannot be used with effective core potentials since the asymptotic correction is switched on in the core region.     

Interesting properties of Thomas–Fermi kinetic and Parr electron–electron‐repulsion DFT energy functional generated compact one‐electron density approximation for ground‐state electronic energy of molecular systems

The reduction of the electronic Schrodinger equation or its calculating algorithm from 4N‐dimensions to a (nonlinear, approximate) density functional of three spatial dimension one‐electron density for an N‐electron system, which is tractable in the practice, is a long desired goal in electronic structure calculation. If the Thomas‐Fermi kinetic energy (～∫ρ^5/3d r ₁) and Parr electron–electron repulsion energy (～∫ρ^4/3d r ₁) main‐term functionals are accepted, and they should, the later described, compact one‐electron density approximation for calculating ground state electronic energy from the 2nd Hohenberg–Kohn theorem is also noticeable, because it is a certain consequence of the aforementioned two basic functionals. Its two parameters have been fitted to neutral and ionic atoms, which are transferable to molecules when one uses it for estimating ground‐state electronic energy. The convergence is proportional to the number of nuclei (M) needing low disc space usage and numerical integration. Its properties are discussed and compared with known ab initio methods, and for energy differences (here atomic ionization potentials) it is comparable or sometimes gives better result than those. It does not reach the chemical accuracy for total electronic energy, but beside its amusing simplicity, it is interesting in theoretical point of view, and can serve as generator function for more accurate one‐electron density models. © 2008 Wiley Periodicals, Inc. J Comput Chem 2009     

Analytical nuclear excited‐state gradients for the Tamm–Dancoff approximation using uncoupled frozen‐density embedding

We report the derivation and implementation of analytical nuclear gradients for excited states using time‐dependent density functional theory using the Tamm–Dancoff approximation combined with uncoupled frozen‐density embedding using density fitting. Explicit equations are presented and discussed. The implementation is able to treat singlet as well as triplet states and functionals using the local density approximation, the generalized gradient approximation, combinations with Hartree–Fock exchange (hybrids), and range‐separated functionals such as CAM‐B3LYP. The new method is benchmarked against supermolecule calculations in two case studies: The solvatochromic shift of the (vertical) fluorescence energy of 4‐aminophthalimide on solvation, and the first local excitation of the benzonitrile dimer. Whereas for the 4‐aminophthalimide–water complex deviations of about 0.2 eV are obtained to supermolecular calculations, for the benzonitrile dimer the maximum error for adiabatic excitation energies is below 0.01 eV due to a weak coupling of the subsystems. © 2017 Wiley Periodicals, Inc.     

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.     

Alternative Representations of the Correlation Energy in Density‐Functional Theory: A Kinetic‐Energy Based Adiabatic Connection

The adiabatic‐connection framework has been widely used to explore the properties of the correlation energy in density‐functional theory. The integrand in this formula may be expressed in terms of the electron–electron interactions directly, involving intrinsically two‐particle expectation values. Alternatively, it may be expressed in terms of the kinetic energy, involving only one‐particle quantities. In this work, we explore this alternative representation for the correlation energy and highlight some of its potential for the construction of new density functional approximations. The kinetic‐energy based integrand is effective in concentrating static correlation effects to the low interaction strength regime and approaches zero asymptotically, offering interesting new possibilities for modeling the correlation energy in density‐functional theory     

Applications of density functional theory methods in millimeter‐wave spectroscopy

Density functional theory (DFT), using the most common functionals, and ab initio quantum chemistry methods are used to calculate the rotational constants and dipole moments of the astrophysically important molecules HCN, CH₃CN, CH₃CNH⁺, HCCCN, and HCCNC. As far as millimeter‐wave spectroscopy is of interest the DFT methods performed well with most functionals, giving results within ±1% of experiments for rotational constants and ±3% for dipole moments. Analyzing the results obtained with all theoretical models, it may be concluded that the Becke's three‐parameter exchange functional and the gradient‐corrected functional of Lee, Yang, and Paar (B3LYP) and Becke's three‐parameter functional with Perdew–Wang correlational functional [B3PW91/6‐31G(d, p)] give the best performances. A detailed analysis of the electron correlation effects shows that HCCCN is more stable than is HCCNC, by 1.16 eV, with important contribution arising from triple excitations. This result is also compared with those obtained with DFT methods. Despite occasional difficulties, DFT with the currently available functionals are of great utility in quickly assessing spectroscopic parameters of astrophysical interest. © 2002 Wiley Periodicals, Inc. Int J Quantum Chem, 2003     

libKEDF: An accelerated library of kinetic energy density functionals

Kinetic energy density functionals (KEDFs) approximate the kinetic energy of a system of electrons directly from its electron density. They are used in electronic structure methods that lack direct access to orbitals, for example, orbital‐free density functional theory (OFDFT) and certain embedding schemes. In this contribution, we introduce libKEDF, an accelerated library of modern KEDF implementations that emphasizes nonlocal KEDFs. We discuss implementation details and assess the performance of the KEDF implementations for large numbers of atoms. We show that using libKEDF, a single computing node or (GPU) accelerator can provide easy computational access to mesoscale chemical and materials science phenomena using OFDFT algorithms. © 2017 Wiley Periodicals, Inc.     

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.     

First-order gradient correction for the exchange-energy density functional for atoms

Summary Spurred by earlier discoveries by Deb and others that a first-order correction in the kinetic energy functional leads to an improved kinetic energy, a first-order gradient term is studied as a component of the gradient-corrected functional for atomic exchange energy. This term is shown to improve the local density approximation to the exchange energy more than does the usual second-order gradient correction. The imperative of systematically deriving this gradient correction is discussed but not resolved.     

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.     

Performance of DFT hybrid functionals in the theoretical treatment of H‐bonds: Analysis term‐by‐term

The performance of DFT to reproduce noncovalent interactions like H‐bonds was assessed. Three DFT hybrid functionals (B3LYP, BHandHLYP from Jaguar5.5, and BHandHLYP from G03) were used to calculate the interaction energies and H‐bond distances of several host–guest complexes theoretically designed. Two reference systems (whose experimental data were available) were also calculated at the same levels of theory. In all the cases B3LYP and BHandHLYP from G03 gave rather the same results of interaction energy and distances, whereas the functional BHandHLYP from Jaguar5.5 overestimated the interaction energies and underestimated the H‐bond distances. Considering the construction of each functional, the terms responsible for such differences are the gradient correction to the exchange functional (Becke88) and the correlation functional (VWN, LYP) and not the Hartree–Fock contribution. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

Determining the role of the underlying orbital‐dependence of PBE0‐DH and PBE‐QIDH double‐hybrid density functionals

We study the orbital‐dependence of three (parameter‐free) double‐hybrid density functionals, namely the PBE0‐DH, the PBE‐QIDH models, and the SOS1‐PBE‐QIDH spin‐opposite‐scaled variant of the latter. To do it, we feed all their energy terms with different sets of orbitals obtained previously from self‐consistent density functional theory calculations using several exchange‐correlation functionals (e.g., PBE, PBE0, PBEH&H), or directly with HF‐PBE orbitals, to see their effect on selected datasets for atomization and reaction energies, the latter proned to marked self‐interaction errors. We find that the PBE‐QIDH double‐hybrid model shows a great consistency, as the best results are always obtained for the set of orbitals corresponding to its hybrid scheme, which prompts us to recommend this model without any other fitting or reparameterization. © 2017 Wiley Periodicals, Inc.     

Parallel implementation of efficient charge–charge interaction evaluation scheme in periodic divide‐and‐conquer density‐functional tight‐binding calculations

A low‐computational‐cost algorithm and its parallel implementation for periodic divide‐and‐conquer density‐functional tight‐binding (DC‐DFTB) calculations are presented. The developed algorithm enables rapid computation of the interaction between atomic partial charges, which is the bottleneck for applications to large systems, by means of multipole‐ and interpolation‐based approaches for long‐ and short‐range contributions. The numerical errors of energy and forces with respect to the conventional Ewald‐based technique can be under the control of the multipole expansion order, level of unit cell replication, and interpolation grid size. The parallel performance of four different evaluation schemes combining previous approaches and the proposed one are assessed using test calculations of a cubic water box on the K computer. The largest benchmark system consisted of 3,295,500 atoms. DC‐DFTB energy and forces for this system were obtained in only a few minutes when the proposed algorithm was activated and parallelized over 16,000 nodes in the K computer. The high performance using a single node workstation was also confirmed. In addition to liquid water systems, the feasibility of the present method was examined by testing solid systems such as diamond form of carbon, face‐centered cubic form of copper, and rock salt form of sodium chloride. © 2017 Wiley Periodicals, Inc.     

Grid‐based density functional calculations of many‐electron systems

Exploratory variational pseudopotential density functional calculations are performed for the electronic properties of many‐electron systems in the 3D cartesian coordinate grid (CCG). The atom‐centered localized gaussian basis set, electronic density, and the two‐body potentials are set up in the 3D cubic box. The classical Hartree potential is calculated accurately and efficiently through a Fourier convolution technique. As a first step, simple local density functionals of homogeneous electron gas are used for the exchange‐correlation potential, while Hay‐Wadt‐type effective core potentials are employed to eliminate the core electrons. No auxiliary basis set is invoked. Preliminary illustrative calculations on total energies, individual energy components, eigenvalues, potential energy curves, ionization energies, and atomization energies of a set of 12 molecules show excellent agreement with the corresponding reference values of atom‐centered grid as well as the grid‐free calculation. Results for three atoms are also given. Combination of CCG and the convolution procedure used for classical Coulomb potential can provide reasonably accurate and reliable results for many‐electron systems. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

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.     

Three‐body expansion of the fragment molecular orbital method combined with density‐functional tight‐binding

The three‐body fragment molecular orbital (FMO3) method is formulated for density‐functional tight‐binding (DFTB). The energy, analytic gradient, and Hessian are derived in the gas phase, and the energy and analytic gradient are also derived for polarizable continuum model. The accuracy of FMO3‐DFTB is evaluated for five proteins, sodium cation in explicit solvent, and three isomers of polyalanine. It is shown that FMO3‐DFTB is considerably more accurate than FMO2‐DFTB. Molecular dynamics simulations for sodium cation in water are performed for 100 ps, yielding radial distribution functions and coordination numbers. © 2017 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.