Multiconfiguration electron density function for the ATSP2K-package

A new ATSP2K module is presented for evaluating the electron density function of any multiconfiguration Hartree-Fock or configuration interaction wave function in the non relativistic or relativistic Breit-Pauli approximation. It is first stressed that the density function is not a priori spherically symmetric in the general open shell case. Ways of building it as a spherical symmetric function are discussed, from which the radial electron density function emerges. This function is written in second quantized coupled tensorial form for exploring the atomic spherical symmetry. The calculation of its expectation value is performed using the angular momentum theory in orbital, spin, and quasispin spaces, adopting a generalized graphical technique. The natural orbitals are evaluated from the diagonalization of the density matrix

Complexity analysis of Klein-Gordon single-particle systems

The Fisher-Shannon complexity is used to quantitatively estimate the contribution of relativistic effects to on the internal disorder of Klein-Gordon single-particle Coulomb systems which is manifest in the rich variety of three-dimensional geometries of its corresponding quantum-mechanical probability density. It is observed that, contrary to the non-relativistic case, the Fisher-Shannon complexity of these relativistic systems does depend on the potential strength (nuclear charge). This is numerically illustrated for pionic atoms. Moreover, its variation with the quantum numbers (n, l, m) is analysed in various ground and excited states. It is found that the relativistic effects enhance when n and/or l are decreasing.Comment: 4 pages, 3 figures, Accepted in EPL (Europhysics Letters

Density Scaling of Noninteracting Kinetic Energy Functionals

The influence of imposing an approximate density scaling condition on a noninteracting kinetic energy functional is investigated. A simple generalized gradient approximation (GGA) is presented, which satisfies both the density scaling condition and the usual coordinate scaling condition; the remaining multiplicative constant is determined from an energy criterion. In post-Kohn–Sham calculations, noninteracting kinetic energies of the closed-shell molecules of the G1 set determined using the GGA are a modest improvement over those determined using the corresponding local functional, which does not satisfy the density scaling condition. Potential energy curves of CO, F2, and P2 exhibit binding with the GGA, compared to purely repulsive curves with the local functional. Adjusting the exponent in the GGA form in order to optimize energy accuracy violates the density scaling condition, and two of the diatomics no longer exhibit binding. Results are compared with those from other local/GGA functionals in the literature

Revisiting the density scaling of the non-interacting kinetic energy

Scaling relations play an important role in the understanding and development of approximate functionals in density functional theory. Recently, a number of these relationships have been redefined in terms of the Kohn–Sham orbitals [Calderín, Phys. Rev. A: At., Mol., Opt. Phys., 2013, 86, 032510]. For density scaling the author proposed a procedure involving a multiplicative scaling of the Kohn–Sham orbitals whilst keeping their occupation numbers fixed. In the present work, the differences between this scaling with fixed occupation numbers and that of previous studies, where the particle number change implied by the scaling was accommodated through the use of the grand canonical ensemble, are examined. We introduce the terms orbital and ensemble density scaling for these approaches, respectively. The natural ambiguity of the density scaling of the non-interacting kinetic energy functional is examined and the ancillary definitions implicit in each approach are highlighted and compared. As a consequence of these differences, Calderín recovered a homogeneity of degree 1 for the non-interacting kinetic energy functional under orbital scaling, contrasting recent work by the present authors [J. Chem. Phys., 2012, 136, 034101] where the functional was found to be inhomogeneous under ensemble density scaling. Furthermore, we show that the orbital scaling result follows directly from the linearity and the single-particle nature of the kinetic energy operator. The inhomogeneity of the non-interacting kinetic energy functional under ensemble density scaling can be quantified by defining an effective homogeneity. This quantity is shown to recover the homogeneity values for important approximate forms that are exact for limiting cases such as the uniform electron gas and one-electron systems. We argue that the ensemble density scaling provides more insight into the development of new functional forms

Molecular Binding in Post-Kohn-Sham Orbital-Free DFT

Molecular binding in post-Kohn–Sham orbital-free DFT is investigated, using noninteracting kinetic energy functionals that satisfy the uniform electron gas condition and which are inhomogeneous under density scaling. A parameter is introduced that quantifies binding, and a series of functionals are determined from fits to near-exact effective homogeneities and/or Kohn–Sham noninteracting kinetic energies. These are then used to investigate the relationship between binding and the accuracy of the effective homogeneity and noninteracting kinetic energy at the equilibrium geometry. For a series of 11 molecules, the binding broadly improves as the effective homogeneity improves, although the extent to which it improves is dependent on the accuracy of the noninteracting kinetic energy; optimal binding appears to require both to be accurate simultaneously. The use of a Thomas–Fermi–von Weizsäcker form, augmented with a second gradient correction, goes some way toward achieving this, exhibiting molecular binding on average. The findings are discussed in terms of the noninteracting kinetic potential and the Hellmann–Feynman theorem. The extent to which the functionals can reproduce the system-dependence of the near-exact effective homogeneity is quantified, and potential energy curves are presented for selected molecules. The study provides impetus for including density scaling homogeneity considerations in the design of noninteracting kinetic energy functionals

Electron density and Fisher information of Dirac-Fock atoms

Numerical calculations on the gradient and Laplacian forms of the position space Fisher information measure are reported using the 1-normalised Dirac-Fock densities (shape function), σ(r), for atoms H-Lr. It is shown that the difference in effective electrostatic potentials, corresponding to the gradient and the Laplacian form of Fishers' information, is completely defined by the shape function (the density per particle) at the nucleus, σ(r=0). The influence of relativistic effects on the Fisher information is recovered for the first time

Characterization of the Chandrasekhar correlated two-electron wavefunction using Fisher, Shannon and statistical complexity information measures

The three-parameter two-electron correlated analytic wavefunction of Chandrasekhar, parametrized using three physically meaningful conditions on the electron density, is assessed using several information theory measures, the Fisher information, Shannon entropy, and statistical complexity, and compared to the results for the same measures using hydrogenic and Hartree-Fock (HF) orbitals

Comparative characterization of two-electron wavefunctions using information-theory measures

Information-theory measures, in particular the Shannon entropy, Fisher information and statistical complexity, are used to discuss the variations among several commonly encountered model two-electron correlated wavefunctions. The Hookean, Moshinsky, and three-parameter Chandrasekhar wavefunctions are considered in real and momentum space, with further comparisons to the Hookean-Hartree-Fock (HF) wavefunction of Ragot, the numerical HF limit, and the hydrogenic (pure Coulomb) limit. The purpose of the study is to quantitatively analyze the effect of different models for inclusion of electron-electron correlation on information-theoretical measures, including statistical complexity, which characterize the electron distribution in position and momentum space

N-derivative of Shannon entropy of shape function for atoms

The Shannon entropy of the ratio of electron density and the number of electrons, shape function entropy, is reported for the atoms He-Ac within the non-relativistic exchange-only optimized effective potential model. The derivative of the shape function entropy with electron number at constant external potential is related to an integral containing the difference between the average Fukui function and the shape function weighted by the logarithm of electron density. The trends in the shape function entropy, its spin analogue and the corresponding derivatives with electron number reveal interesting periodic behaviour