Step structure in the atomic Kohn-Sham potential

In this work we analyze the exchange-correlation potentialv _xc within the Kohn-Sham approach to density functional theory for the case of atomic systems. The exchange-correlation potential is written as the sum of two potentials. One of these potentialsv _xc,scr is the long-range. Coulombic potential of the coupling constant integrated exchange-correlation hole which represents the screening of the two-particle interactions due to exchange-correlation effects. The other potentialv _xc,scr ^resp contains the functional derivative with respect to the electron density of the coupling constant integrated pair-correlation function representing the sensitivity of this exchange-correlation screening to density variations. As explicit expression of the exchange-part of this functional derivative is derived using an approximation for the Greens function of the Kohn-Sham system and is shown to display a distinct atomic shell structure. The corresponding potentialv _xc,scr ^resp has a clear step structure and is constant within the atomic shells and changes rapidly at the atomic shell boundaries. Numerical examples are presented for the Be and Kr atoms using the Optimized Potential Model (OPM).     

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.     

Atomic shell structure based on inhomogeneity measures of the electron density

A measure of electron density inhomogeneity is proposed based on the distance of function values to the average value within a chosen region. The choice of the examined regions follows the approach of restricted space partitioning with fixed inhomogeneity as the restriction. The integration of electron density over the regions of the partitioning yields discrete distribution of charges, which is analyzed. The resulting functionals depend on the particular definition of distance. Two possibilities are selected for the distance definition, and the functionals are applied to examine the shell structure of the atoms Li to Xe.     

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.     

Gradient expansion of the atomic kinetic energy functional

It is demonstrated that the ground-state atomic kinetic energy functional T[?], where ? is the electron density, can be computed to surprising accuracy from the truncated gradient expansion: T[?] = + T₂[?] + T₄[?], with T_o[?] = $\text{3}\text{10}$(3π²)$\text{2}\text{3}$ ∫ ?$\text{5}\text{3}$ d_τ, T₂ [?] = $\text{1}\text{72}$ ∫ (??)²?^?1 d_τ, and T₄ [?] given by the formula of Hodges. Calculations of T₀, T₂ and T₄ are reported for He with ? both the Hartree—Fock and a very accurate density, and for Ne, Ar and Kr with ? the Hartree—Fock density. For Kr, T₀ + T₂ + T₄ is within 0.3% of the exact Hartree—Fock T, with T₂/T₀ = 0.05, T₄/T₂ = 0.17.     

Topological analysis of electron densities from Kohn-Sham and subsystem density functional theory

In this study, we compare the electron densities for a set of hydrogen-bonded complexes obtained with either conventional Kohn-Sham density functional theory (DFT) calculations or with the frozen-density embedding (FDE) method, which is a subsystem approach to DFT. For a detailed analysis of the differences between these two methods, we compare the topology of the electron densities obtained from Kohn-Sham DFT and FDE in terms of deformation densities, bond critical points, and the negative Laplacian of the electron density. Different kinetic-energy functionals as needed for the frozen-density embedding method are tested and compared to a purely electrostatic embedding. It is shown that FDE is able to reproduce the characteristics of the density in the bonding region even in systems such as the F-H-F(-) molecule, which contains one of the strongest hydrogen bonds. Basis functions on the frozen system are usually required to accurately reproduce the electron densities of supermolecular calculations. However, it is shown here that it is in general sufficient to provide just a few basis functions in the boundary region between the two subsystems so that the use of the full supermolecular basis set can be avoided. It also turns out that electron-density deformations upon bonding predicted by FDE lack directionality with currently available functionals for the nonadditive kinetic-energy contribution.     

Properties of the atomic form factor: Applications and shell structure

A relationship between the atomic form factor, F(k), and its first derivative, both at the origin, is presented. A new function f(k), related in a simple way to F(k), has been studied and some applications have been performed. They led us to find lower bounds to F(k) for all k and to other quantities such as the charge density at the origin, ρ(0), and radial expectation values. Finally, interesting effects on the Laplacian of f(k) due to the atomic shell structure were found. © 1995 John Wiley & Sons, Inc.     

Dynamic evolution of Kohn-Sham electron density in the real-time domain with finite basis expansion.

We present a novel method for time-dependent density functional theory calculations on dynamic linear response and electron density evolution in the real-time domain with the finite basis expansion approach of conventional quantum chemistry. To demonstrate the validity and efficiency of this method, dynamic polarizabilities of a water chain and diphenylene molecules are computed by employing the Chebyshev interpolation algorithm, which was developed by Baer and co-workers. The calculated dynamic polarizabilities show good agreement with those obtained from conventional linear response calculations. The density evolution in the real-time domain with application of a long-duration electric field gives electronic conduction in molecules, where a dynamic process of charge transfer is observed with the snapshots of response density in real time. Charge transfer oscillating with the frequency of the highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) gap is shown in a diphenylene molecule while there is little change in time for a water chain.     

The structure of the ion: An ion kinetic energy spectrometry study

The kinetic energy release (T) values for the loss of CO from ions from five compounds have been obtained and are consistent with at least two different structures for the ion. This result is supported by the T value obtained for the decomposition of the molecular ions.     

Structure of the electron momentum density of atomic systems

The present paper addresses the controversial problem on the nonmonotonic behavior of the spherically-averaged momentum density γ(p) observed previously for some ground-state atoms based on the Roothaan-Hartree-Fock (RHF) wave functions of Clementi and Roetti. Highly accurate RHF wave functions of Koga et al. are used to study the existence of extrema in the momentum density γ(p) of all the neutral atoms from hydrogen to xenon. Three groups of atoms are clearly identified according to the nonmonotonicity parameter μ, whose value is either equal to, larger, or smaller than unity. Additionally, it is found that the function p^?α γ(p) is (i) monotonically decreasing from the origin for α ≥ 0.75, (ii) convex for α ≥ 1.35, and (iii) logarithmically convex for α ≥ 3.64 for all the neutral atoms with nuclear charges Z = 1–54. Finally, these monotonicity properties are applied to derive simple yet general inequalities which involve three momentum moments 〈p ^t≥. These inequalities not only generalize similar inequalities reported so far but also allow us to correlate some fundamental atomic quantities, such as the electron-electron repulsion energy and the peak height of Compton profile, in a simple manner.     

The effect of electron correlation on the shell structure of atoms

The effect of electron correlation on the spherically averaged electron density (r) and the radial density distribution function D(r), exhibiting the shell structure of atoms, is studied, using basis sets consisting of contracted Gaussians, for the atoms Li through Ar. These effects were shown to be of smaller importance as compared to the effect associated to basis set extension.     

Introduction of the shell structure into the Thomas—Fermi energy density functional for neutral atoms

The Thomas—Fermi energy density functional is constructed for the Na atom using H-like wavefunctions. Minimization of the functional with respect to the Z in the wavefunctions leads to the energy value ot ?166.35 e_O²/a_O, which compares favorably with the Hartree—Fock value of ?161.8 e_O²/a_O.     

Towards an exact orbital-free single-particle kinetic energy density for the inhomogeneous electron liquid in the Be atom

Holas and March [Phys. Rev. A51, 2040 (1995)] wrote the gradient of the one-body potential V(r) in terms of low-order derivatives of the idempotent Dirac density matrix built from a single Slater determinant of Kohn–Sham orbitals. Here, this is first combined with the study of Dawson and March [J. Chem. Phys. 81, 5850 (1984)] to express the single-particle kinetic energy density of the Be atom ground-state in terms of both the electron density n(r) and potential V(r). While this is the more compact formulation, we then, by removing V(r), demonstrate that the ratio t(r)/n(r) depends, though non-locally, only on the single variable n′(r)/n(r), no high-order gradients entering for the spherical Be atom.     

Zero electron kinetic energy spectroscopy of the para-fluorotoluene cation

Zero electron kinetic energy (ZEKE) spectroscopy is employed to gain information on the vibrational energy levels of the para-fluorotoluene (pFT) cation. Vibrationally resolved spectra are obtained following excitation through a range of intermediate vibrational energy levels in the S1 state. These spectra allow the observation of different cationic vibrational modes, whose assignment is achieved both from a knowledge of the S1 vibrational states and also by comparison with density functional calculations. In one notable case, clean ZEKE spectra were obtained from two overlapped S1 features. From the authors' data, the adiabatic ionization energy of pFT was derived as 70,946+/-4 cm(-1). The information on the cationic energy levels obtained will be useful in untangling the intramolecular vibrational redistribution dynamics of pFT in the S1 state.     

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.     

Charge localization: The ion kinetic energy spectra of pyrazine,pyrimidine and pyridazine

The ion kinetic energy (IKE) spectra of pyrazine, pyrimidine, pyridazine, pyrazine-d₄, pyrimidine-d₄ and pyridazine-d₄ have shown that these compounds can more readily be distinguished on the basis of IKE spectra than mass spectra. Studies of the fragmentation of various doublycharged ions have given information about the localization of charge on the heteroatoms and on the structures of the ions themselves. Transition states from which fragmentation occurs for doublycharged ions with loss of either [H]⁺ or [D]⁺ have shown large differences of the order 1 Å between the intercharge distances in the two cases.     

Vibrationally promoted electron emission at a metal surface: electron kinetic energy distributions

We report the first direct measurement of the kinetic energy of exoelectrons produced by collisions of vibrationally excited molecules with a low work function metal surface exhibiting electron excitations of 64% (most probable) and 95% (maximum) of the initial vibrational energy. This remarkable efficiency for vibrational-to-electronic energy transfer is in good agreement with previous results suggesting the coupling of multiple vibrational quanta to a single electron.     

The effect of electron correlation on the topological and atomic properties of the electron density distributions of molecules

The topological properties of the electron density and the properties of an atom in a molecule are calculated by means of second-order Møller-Plesset perturbation theory (MP2) and compared with the results of configuration interaction calculations (C12) which include all single and double substitutions from the Hartree-Fock reference configuration. A software package for analyzing the effects of electron correlation on the topological properties of the electron density of molecules is described. H₂CO is used to provide a numerical example and to indicate that the number of bond critical points is unaffected by the inclusion of electron correlation. Correlation leads to only a small shift in the positions of bond critical points and a small change in the electron density at bond critical points. It is further shown that the energy of an atom in a molecule can be calculated to an accuracy of 1 kcal/mol and the electron population of an atom to about 0.001e. A statistical method is used to show that the deviation of the MP2 correlation correction relative to the CI2 correlation correction for a variety of atomic properties is about 25%.