Gradient expansion of the kinetic energy density functional: Local behavior of the kinetic energy density

Calculations with Hartree—Fock electron densities for the rare gas atoms He through Xe show that the gradient expansion for the kinetic energy functional, T[] = T₀[] + T₂[] + T₄[] + … = ∫t() dτ, approximates the kinetic energy by averaging over the shell structure present in the true local kinetic energy density, t(), and that the accuracy of the gradient expansion improves with increasing atomic number. Components of t(), t₀(), t₂() and t₄(), are exhibited and discussed. The defined function t() is everywhere positive.     

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.     

Chemical bonding in view of electron charge density and kinetic energy density descriptors

Stalke's dilemma, stating that different chemical interpretations are obtained when one and the same density is interpreted either by means of natural bond orbital (NBO) and subsequent natural resonance theory (NRT) application or by the quantum theory of atoms in molecules (QTAIM), is reinvestigated. It is shown that within the framework of QTAIM, the question as to whether for a given molecule two atoms are bonded or not is only meaningful in the context of a well‐defined reference geometry. The localized‐orbital‐locator (LOL) is applied to map out patterns in covalent bonding interaction, and produces results that are consistent for a variety of reference geometries. Furthermore, LOL interpretations are in accord with NBO/NRT, and assist in an interpretation in terms of covalent bonding. © 2008 Wiley Periodicals, Inc.J Comput Chem, 2009.     

Kinetic energy density study of some representative semilocal kinetic energy functionals

There is a number of explicit kinetic energy density functionals for noninteracting electron systems that are obtained in terms of the electron density and its derivatives. These semilocal functionals have been widely used in the literature. In this work, we present a comparative study of the kinetic energy density of these semilocal functionals, stressing the importance of the local behavior to assess the quality of the functionals. We propose a quality factor that measures the local differences between the usual orbital-based kinetic energy density distributions and the approximated ones, allowing us to ensure if the good results obtained for the total kinetic energies with these semilocal functionals are due to their correct local performance or to error cancellations. We have also included contributions coming from the Laplacian of the electron density to work with an infinite set of kinetic energy densities. For all but one of the functionals, we have found that their success in the evaluation of the total kinetic energy is due to global error cancellations, whereas the local behavior of their kinetic energy density becomes worse than that corresponding to the Thomas-Fermi functional.     

The carbonyl stretch of furfural: Fermi resonance and rotational isomerism

The i.r. region of the carbonyl stretch of furfural is interpreted in terms of Fermi resonance and rotational isomerism. The Fermi doublets for each rotamer are assigned and the unperturbed frequencies and Fermi coupling coefficients W, are evaluated by the method of solvent variation.     

Structure of isolated and solvated peroxyl radicals

We have investigated the structure of HO₂ and a series of alkyl peroxyl radicals ROO using a variety of quantum mechanical methods. We first compute the geometries, vibrational frequencies, electronic charge distributions, and spin densities for the series of radicals considered in the gas phase. Significant differences with respect to previous calculations have been pointed out in a few cases. In particular, we show the fundamental importance of electronic correlation when computing net atomic charges and spin densities, which have generally been estimated in the litterature by means of Hartree–Fock SCF electronic densities. Solvation effects on the geometry and electronic structure have been estimated by carrying out self-consistent reaction field computations in a polarizable continuum environment with relative dielectric permittivity equal to that of liquid water. Large electronic polarization is predicted in such conditions. This may be important in order to understand reactive properties of the radicals in different media. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 1039–1048, 1999     

Constraint-based analysis of a physics-guided kinetic energy density expansion

An approach guided by physical consistency in determining the general forms of D-dimensional kinetic energy density functionals (KEDF) has been demonstrated previously, producing an expansion which contains the majority of the known one-point KEDF forms. It has also been shown that any noninteracting KEDF must necessarily have a homogeneity degree of 2 in coordinate scaling, and that the ratio of the collective KED to electron density must approach the ionization energy as $r\to \infty$. This article demonstrates that the scaling condition is already satisfied in the general expansion despite not being conceived with the scaling as a constraint, and that the second condition places a restriction on the expansion terms of the KED. The discussion is extended as well for some known KEDs for comparison.     

Theoretical study of atoms by the electronic kinetic energy density and stress tensor density

The electronic structure of atoms in the first, second, and third periods were analyzed using the electronic kinetic energy density and stress tensor density, which are local quantities motivated by quantum field theoretic consideration, specifically the rigged quantum electrodynamics. The zero surfaces of the electronic kinetic energy density, which are called as the electronic interfaces, of the atoms were computed. It was found that their sizes exhibited clear periodicity and were comparable to the conventional atomic and ionic radii. The electronic stress tensor density and its divergence, tension density, of the atoms, were also computed and how their electronic structures were characterized by them was discussed. © 2016 Wiley Periodicals, Inc.     

Approximate density matrices and wigner distribution functions from density,kinetic energy density,and idempotency constraints

The Wigner distribution function and the corresponding density matrix are calculated using a form for the distribution function suggested by maximization of the entropy. Wigner functions and density matrices are determined by imposing conditions of idempotency on the density matrix. Exchange energies and Compton profiles calculated from density matrices obtained by imposing the idempotency constraints are compared with the results of calculations using the Hartree–Fock density matrix and a Gaussian approximation for the density matrix for H and the noble gases He through Xe. Compton profiles from Wigner functions with idempotency constraints show improvements over the Gaussian approximation for the lighter atoms, but do not show significant changes for the heavier atoms. Exchange energies from density matrices with idempotency constraints show improvements over the Gaussian approximation except for the heavier atoms Kr and Xe.     

On the evaluation of the non-interacting kinetic energy in density functional theory

The utility of both an orbital-free and a single-orbital expression for computing the non-interacting kinetic energy in density functional theory is investigated for simple atomic systems. The accuracy of both expressions is governed by the extent to which the Kohn-Sham equation is solved for the given exchange-correlation functional and so special attention is paid to the influence of finite Gaussian basis sets. The orbital-free expression is a statement of the virial theorem and its accuracy is quantified. The accuracy of the single-orbital expression is sensitive to the choice of Kohn-Sham orbital. The use of particularly compact orbitals is problematic because the failure to solve the Kohn-Sham equation exactly in regions where the orbital has decayed to near-zero leads to unphysical behaviour in regions that contribute to the kinetic energy, rendering it inaccurate. This problem is particularly severe for core orbitals, which would otherwise appear attractive due to their formally nodeless nature. The most accurate results from the single-orbital expression are obtained using the relatively diffuse, highest occupied orbitals, although special care is required at orbital nodes.     

Theoretical study of lithium ionic conductors by electronic stress tensor density and electronic kinetic energy density

We analyze the electronic structure of lithium ionic conductors, and , using the electronic stress tensor density and kinetic energy density with special focus on the ionic bonds among them. We find that, as long as we examine the pattern of the eigenvalues of the electronic stress tensor density, we cannot distinguish between the ionic bonds and bonds among metalloid atoms. We then show that they can be distinguished by looking at the morphology of the electronic interface, the zero surface of the electronic kinetic energy density. © 2016 Wiley Periodicals, Inc.     

Interatomic interactions in momentum space. Momentum density and kinetic energy in chemical bonding

On the basis of the virial theorem for a uniform scaling process of a polyatomic system, the total energy and its gradient are quantitatively related with the behavior of the electron density in momentum space through the kinetic energy of the system. For attractive and repulsive interactions, the behavior of the momentum density distribution and its effect on the stabilization energy and the interatomic force are examined. Some guiding principles are deduced for their interrelation. The results are used to clarify the role of kinetic energy in chemical bonding. Possible energy partitioning in this approach is also mentioned.     

Effective homogeneity of the exchange-correlation and non-interacting kinetic energy functionals under density scaling

Correlated electron densities, experimental ionisation potentials, and experimental electron affinities are used to investigate the homogeneity of the exchange-correlation and non-interacting kinetic energy functionals of Kohn-Sham density functional theory under density scaling. Results are presented for atoms and small molecules, paying attention to the influence of the integer discontinuity and the choice of the electron affinity. For the exchange-correlation functional, effective homogeneities are highly system-dependent on either side of the integer discontinuity. By contrast, the average homogeneity-associated with the potential that averages over the discontinuity-is generally close to 4/3 when the discontinuity is computed using positive affinities for systems that do bind an excess electron and negative affinities for those that do not. The proximity to 4/3 becomes increasingly pronounced with increasing atomic number. Evaluating the discontinuity using a zero affinity in systems that do not bind an excess electron instead leads to effective homogeneities on the electron abundant side that are close to 4/3. For the non-interacting kinetic energy functional, the effective homogeneities are less system-dependent and the effect of the integer discontinuity is less pronounced. Average values are uniformly below 5/3. The study provides information that may aid the development of improved exchange-correlation and non-interacting kinetic energy functionals.     

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.     

Non-Gaussian statistics of amide I mode frequency fluctuation of N-methylacetamide in methanol solution: linear and nonlinear vibrational spectra

By carrying out molecular dynamics simulations of an N-methylacetamide (NMA) in methanol solution, the amide I mode frequency fluctuation and hydrogen bonding dynamics were theoretically investigated. Combining an extrapolation formula developed from systematic ab initio calculation studies of NMA-(CH3OH)n clusters with a classical molecular dynamics simulation method, we were able to quantitatively describe the solvatochromic vibrational frequency shift induced by the hydrogen-bonding interaction between NMA and solvent methanol. It was found that the fluctuating amide I mode frequency distribution is notably non-Gaussian and it can be decomposed into two Gaussian peaks that are associated with two distinctively different solvation structures. The ensemble-average-calculated linear response function associated with the IR absorption is found to be oscillating, which is in turn related to the doublet amide I band shape. Numerically calculated infrared absorption spectra are directly compared with experiment and the agreement was found to be excellent. By using the Onsager's regression hypothesis, the rate constants of the interconversion process between the two solvation structures were obtained. Then, the nonlinear response functions associated with two-dimensional infrared pump-probe spectroscopy were simulated. The physics behind the two-dimensional line shape and origin of the cross peaks in the time-resolved pump-probe spectra is explained and the result is compared with 2D spectra experimentally measured recently by Woutersen et al.     

IR spectral fingerprints of carbonyl groups in graphite oxide: a theoretical study

1. Download : Download high-res image (134KB)
2. Download : Download full-size image