Local effective potential theory: nonuniqueness of potential and wave function

In local effective potential energy theories such as the Hohenberg-Kohn-Sham density functional theory (HKS-DFT) and quantal density functional theory (Q-DFT), electronic systems in their ground or excited states are mapped to model systems of noninteracting fermions with equivalent density. From these models, the equivalent total energy and ionization potential are also obtained. This paper concerns (i) the nonuniqueness of the local effective potential energy function of the model system in the mapping from a nondegenerate ground state, (ii) the nonuniqueness of the local effective potential energy function in the mapping from a nondegenerate excited state, and (iii) in the mapping to a model system in an excited state, the nonuniqueness of the model system wave function. According to nondegenerate ground state HKS-DFT, there exists only one local effective potential energy function, obtained as the functional derivative of the unique ground state energy functional, that can generate the ground state density. Since the theorems of ground state HKS-DFT cannot be generalized to nondegenerate excited states, there could exist different local potential energy functions that generate the excited state density. The constrained-search version of HKS-DFT selects one of these functions as the functional derivative of a bidensity energy functional. In this paper, the authors show via Q-DFT that there exist an infinite number of local potential energy functions that can generate both the nondegenerate ground and excited state densities of an interacting system. This is accomplished by constructing model systems in configurations different from those of the interacting system. Further, they prove that the difference between the various potential energy functions lies solely in their correlation-kinetic contributions. The component of these functions due to the Pauli exclusion principle and Coulomb repulsion remains the same. The existence of the different potential energy functions as viewed from the perspective of Q-DFT reaffirms that there can be no equivalent to the ground state HKS-DFT theorems for excited states. Additionally, the lack of such theorems for excited states is attributable to correlation-kinetic effects. Finally, they show that in the mapping to a model system in an excited state, there is a nonuniqueness of the model system wave function. Different wave functions lead to the same density, each thereby satisfying the sole requirement of reproducing the interacting system density. Examples of the nonuniqueness of the potential energy functions for the mapping from both ground and excited states and the nonuniqueness of the wave function are provided for the exactly solvable Hooke's atom. The work of others is also discussed.     

Pair Function and Structure Factor of an Interacting Electron Liquid

Abstract

Recent work on small angle scattering from liquid metals has caused renewed interest in the electron pair function in the uniform interacting electron liquid, jellium. Therefore we have re-examined this problem, starting from an analysis of the exchange hole, in which the only correlations are due to the Pauli Principle and solely therefore between parallel spin electrons. The pair function g(r) of noninteracting Fermions is expressed in terms of the density of the p-component in the free electron density matrix. This motivates the treatment of the Coulomb repulsion via a potential energy V(r) To close the theory, one must either invoke self-consistency to determine V(r), or relate it to the (direct) correlation function c(r) as in classical liquids. Both methods are briefly considered; the second has the advantage that here the collective plasma oscillations can be introduced through their zero-point energy.     

Theoretical studies on the ionization potential of interacting atoms at large separations

The theory of the one-particle Green's function is applied to calculations of the ionization potential of interacting atoms which are at large separations. Equations for the ionization potential involve terms which relate to Van der Waals interactions between separated atoms and long-range interactions between an atom and an ion. Numerical calculations of the ionization potential of two hydrogen atoms and two helium atoms at large separations are performed. Applications to the ionization potentials of weakly-interacting Van der Waals molecules (NeAr, NeKr, NeXe) are also reported.     

Koopmans' theorem for large molecular systems within density functional theory

It is shown that in density functional theory (DFT), Koopmans' theorem for a large molecular system can be stated as follows: The ionization energy of the system equals the negative of the highest occupied molecular orbital (HOMO) energy plus the Coulomb electrostatic energy of removing an electron from the system, or equivalently, the ionization energy of an N-electron system is the negative of the arithmetic average of the HOMO energy of this system and the lowest unoccupied molecular orbital (LUMO) energy of the (N - 1)-electron system. Relations between this DFT Koopmans' theorem and its existing counterparts in the literature are discussed. Some of the previous results are generalized and some are simplified. DFT calculation results of a fullerene molecule, a finite single-walled carbon nanotube and a finite boron nitride nanotube are presented, indicating that this Koopmans' theorem approximately holds, even if the orbital relaxation is taken into consideration.     

Time-dependent multiconfiguration theory for electronic dynamics of molecules in intense laser fields: a description in terms of numerical orbital functions

The equations of motion (EOMs) for spin orbitals in the coordinate representation are derived within the framework of the time-dependent multiconfiguration theory developed for electronic dynamics of molecules in intense laser fields. We then tailor the EOMs for diatomic (or linear) molecules to apply the theory to the electronic dynamics of a hydrogen molecule in an intense, near-infrared laser field. Numerical results are presented to demonstrate that the time-dependent numerical multiconfiguration wave function is able to describe the correlated electron motions as well as the ionization processes of a molecule in intense laser fields.     

Effective potential energy curves of H2 molecule evolving in a strong time‐dependent magnetic field

Evolution of hydrogen molecule, starting initially from its field‐free ground state, in a time‐dependent (TD) magnetic field of order 10¹¹ G is presented in a parallel internuclear axis and magnetic field‐axis configuration. Effective potential energy curves (EPECs), in terms of exchange and correlation energy, of the hydrogen molecule as a function of TD magnetic‐field strength, are analyzed through TD density functional computations based on a quantum fluid dynamics approach. The numerical computations are performed for internuclear separation R ranging from 0.1 to 14.0 a.u. The EPECs exhibit field‐dependent significant potential‐well minima both at large internuclear separations and at short internuclear separations with a considerable increase in the exchange and correlation energy of the hydrogen molecule. The results, when compared with the time‐independent (TI) studies involving static TI magnetic fields, reveal TD behavior of field‐dependent crossovers between different spin‐states of hydrogen molecule as indicated by the TI investigations in static magnetic fields. Besides this, present work reveals interesting dynamics in the TD total‐electronic charge‐density distribution of the hydrogen molecule. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Normalization and Fermi–Coulomb and Coulomb hole sum rules for approximate wave functions

For approximate wave functions, we prove the theorem that there is a one‐to‐one correspondence between the constraints of normalization and of the Fermi–Coulomb and Coulomb hole charge sum rules at each electron position. This correspondence is surprising in light of the fact that normalization depends on the probability of finding an electron at some position. In contrast, the Fermi–Coulomb hole sum rule depends on the probability of two electrons staying apart because of correlations due to the Pauli exclusion principle and Coulomb repulsion, while the Coulomb hole sum rule depends on Coulomb repulsion. We demonstrate the theorem for the ground state of the He atom by the use of two different approximate wave functions that are functionals rather than functions. The first of these wave function functionals is constructed to satisfy the constraint of normalization, and the second that of the Coulomb hole sum rule for each electron position. Each is then shown to satisfy the other corresponding sum rule. The significance of the theorem for the construction of approximate “exchange‐correlation” and “correlation” energy functionals of density functional theory is also discussed. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Electronic energy functionals: Levy–Lieb principle within the ground state path integral quantum Monte Carlo

We propose a theoretical/computational protocol based on the use of the Ground State Path Integral Quantum Monte Carlo for the calculation of the kinetic and Coulomb energy density for a system of N interacting electrons in an external potential. The idea is based on the derivation of the energy densities via the (N ? 1)‐conditional probability density within the framework of the Levy–Lieb constrained search principle. The consequences for the development of energy functionals within the context of density functional theory are discussed. We propose also the possibility of going beyond the energy densities and extend this idea to a computational procedure where the (N ? 1)‐conditional probability is an implicit functional of the electron density, independently from the external potential. In principle, such a procedure paves the way for an on‐the‐fly determination of the energy functional for any system. © 2012 Wiley Periodicals, Inc.     

Multipole polarizabilities of spherically confined hydrogen atom and positronium in screened Coulomb potential

The static multipole polarizabilities of the hydrogen atom interacting with a screened Coulomb potential confined in an impenetrable spherical box are calculated in the sum-over-states formalism. The system eigenenergies and wave functions are produced by employing the generalized pseudospectral method. High-precision results are obtained to resolve the discrepancies between previous predictions, and used to analyze the critical behavior of the system with varying the confinement radius and screening parameter. A scaling law of the multipole polarizabilities with respect to the nuclear charge and electron reduced mass is proposed. Based on the scaling law we extend the investigation to the spatially confined positronium atom in a screened Coulomb potential by utilizing a reduced one-electron model. The results are in good agreement with previous estimates. However, it is conjectured that when the confinement radius is small the confined positronium atom should be treated more reasonably by an effective two-electron model.     

多组态含时自旋非限制HartreeFock理论

基于自旋非限制HartreeFock理论,发展了自旋非限制多组态含时HartreeFock理论方法来研究激光场中的多电子相关动力学．自旋向上和自旋向下的自旋轨道分别在他们各自的子空间内传播;并通过约化密度矩阵和平均场算符相互作用．分别利用了自旋限制和非限制的多组态含时HartreeFock方法虚时和实时传播计算氦原子基态能量和电离几率．自旋非限制的计算结果与其他报道相吻合．     

Atomic and electronic structures of neutral and charged Pbn clusters (n=2-15): theoretical investigation based on density functional theory

The geometric and electronic structures of the Pbn+ clusters (n=2-15) have been investigated and compared with neutral clusters. The search for several low-lying isomers was carried out under the framework of the density functional theory formalism using the generalized gradient approximation for the exchange correlation energy. The wave functions were expanded using a plane wave basis set and the electron-ion interactions have been described by the projector augmented wave method. The ground state geometries of the singly positively charged Pbn+ clusters showed compact growth pattern as those observed for neutrals with small local distortions. Based on the total energy of the lowest energy isomers, a systematic analysis was carried out to obtain the physicochemical properties, viz., binding energy, second order difference in energy, and fragmentation behavior. It is found that n=4, 7, 10, and 13 clusters are more stable than their neighbors, reflecting good agreement with experimental observation. The chemical stability of these clusters was analyzed by evaluating their energy gap between the highest occupied and lowest unoccupied molecular orbitals and adiabatic ionization potentials. The results revealed that, although Pb13 showed higher stability from the total energy analysis, its energy gap and ionization potential do not follow the trend. Albeit of higher stability in terms of binding energy, the lower ionization potential of Pb13 is interesting which has been explained based on its electronic structure through the density of states and electron shell filling model of spherical clusters.     

Nonlinear response properties of atomic hydrogen under quantum plasma environment: A time-dependent variation perturbation study on hyperpolarizability and two-photon excitations

Pilot calculations on the frequency-dependent nonlinear response property, viz. the electric dipole hyperpolarizability of atomic hydrogen under quantum plasma environment, have been performed using an external oscillatory electric field. Fourth-order perturbation theory within a variational scheme is adopted to obtain the hyperpolarizability within and beyond normal dispersion region. Two-photon absorption from the ground state is explicitly obtained from the pole positions of nonlinear response of the system and studied up to principal quantum number n = 4 . Ground and perturbed wave functions of appropriate symmetries are represented by linear combination of Slater-type orbitals. Exponential cosine-screened Coulomb potential is used to simulate the quantum plasma environment. With respect to plasma strength, the nonlinear response properties are considerably enhanced. Results are compared with those under classical plasma environment represented by screened Coulomb potential. Departure from Coulomb potential results in lifting of the accidental degeneracy in the respective two-photon excited states beyond n = 2 . For free hydrogen atom, the transition energies and the radial density profiles of the respective two-photon excited states match exactly with those obtained from analytical wave functions.     

Towards a force field based on density fitting

Total intermolecular interaction energies are determined with a first version of the Gaussian electrostatic model (GEM-0), a force field based on a density fitting approach using s-type Gaussian functions. The total interaction energy is computed in the spirit of the sum of interacting fragment ab initio (SIBFA) force field by separately evaluating each one of its components: electrostatic (Coulomb), exchange repulsion, polarization, and charge transfer intermolecular interaction energies, in order to reproduce reference constrained space orbital variation (CSOV) energy decomposition calculations at the B3LYP/aug-cc-pVTZ level. The use of an auxiliary basis set restricted to spherical Gaussian functions facilitates the rotation of the fitted densities of rigid fragments and enables a fast and accurate density fitting evaluation of Coulomb and exchange-repulsion energy, the latter using the overlap model introduced by Wheatley and Price [Mol. Phys. 69, 50718 (1990)]. The SIBFA energy scheme for polarization and charge transfer has been implemented using the electric fields and electrostatic potentials generated by the fitted densities. GEM-0 has been tested on ten stationary points of the water dimer potential energy surface and on three water clusters (n = 16,20,64). The results show very good agreement with density functional theory calculations, reproducing the individual CSOV energy contributions for a given interaction as well as the B3LYP total interaction energies with errors below kBT at room temperature. Preliminary results for Coulomb and exchange-repulsion energies of metal cation complexes and coupled cluster singles doubles electron densities are discussed.     

Bosonised DFT potential estimated from QMC calculations of the ground-state density for the inhomogeneous electron liquid in Be

To avoid the solution of numerous Kohn–Sham one-body potential equations for wave functions in density functional theory, various groups independently proposed the use of Pauli potential to bosonise the customary one-body potential theory. Here, we utilise our recent quantum Monte Carlo calculations of the ground-state electron density of the Be atom to estimate the bosonised one-body potential V_B(r) and hence extract the Pauli potential for this atom.     

Approximating the Pauli Potential in Bound Coulomb Systems

It is shown that the Pauli potential in bound Coulomb systems can in good approximation be composed from the corresponding atomic fragments. This provides a simple and fast procedure how to generate the Pauli potential in bound systems, which is needed to perform an orbital‐free density functional calculation. The method is applicable to molecules and solids. © 2016 Wiley Periodicals, Inc.     

The Relationship of the Energy of Interaction of Alkali Metal Cations with an Aprotic Solvent Molecule with Quantum Topological Electron Density Characteristics

The optimal geometry and wave functions of the complexes [M(Solv)]⁺ (M = Li, Na, K; Solv is an aprotic solvent molecule) were calculated and the topological characteristics of the electron density distribution at the (3,–1) critical points (CP) of ion–molecule bonds were analyzed by the density functional theory in the B3LYP/6-31+G(d, p) approximation. The parametric dependences for the energy of ion–molecule bonds in terms of the local kinetic and potential electron energy densities at the bond CTs were proposed.     

Exact electronic properties in the classically forbidden region of a metal surface

We derive exact properties of the inhomogeneous electron gas in the asymptotic classically forbidden region at a metal–vacuum interface within the framework of local effective potential energy theory. We derive a new expression for the asymptotic structure of the Kohn–Sham density functional theory (KS‐DFT) exchange‐correlation potential energy v_xc(r) in terms of the irreducible electron self‐energy. We also derive the exact asymptotic structure of the orbitals, density, the Dirac density matrix, the kinetic energy density, and KS exchange energy density. We further obtain the exact expression for the Fermi hole and demonstrate its structure in this asymptotic limit. The exchange‐correlation potential energy is derived to be v_xc(z → ∞) = ?α_KS,xc/z, and its exchange and correlation components to be v_x(z → ∞) = ?α_KS,x/z and v_c(z → ∞) = ?α_KS,c/z, respectively. The analytical expressions for the coefficients α_KS,xc and α_KS,x show them to be dependent on the bulk‐metal Wigner–Seitz radius and the barrier height at the surface. The coefficient α_KS,c = 1/4 is determined in the plasmon‐pole approximation and is independent of these metal parameters. Thus, the asymptotic structure of v_xc(z) in the vacuum region is image‐potential‐like but not the commonly accepted one of ?1/4z. Furthermore, this structure depends on the properties of the metal. Additionally, an analysis of these results via quantal density functional theory (Q‐DFT) shows that both the Pauli W_x(z → ∞) and lowest‐order correlation‐kinetic W(z → ∞) components of the exchange potential energy v_x(z → ∞), and the Coulomb W_c(z → ∞) and higher‐order correlation‐kinetic components of the correlation potential energy v_c(z → ∞), all contribute terms of O(1/z) to the structure. Hence correlations attributable to the Pauli exclusion principle, Coulomb repulsion, and correlation‐kinetic effects all contribute to the asymptotic structure of the effective potential energy at a metal surface. The relevance of the results derived to the theory of image states and to KS‐DFT is also discussed. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005