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     

On the exchange-hole model of London dispersion forces

First-principles derivation is given for the heuristic exchange-hole model of London dispersion forces by Becke and Johnson [J. Chem. Phys. 122, 154104 (2005)]. A one-term approximation is used for the dynamic charge density response function, and it is shown that a central nonempirical ingredient of the approximate nonexpanded dispersion energy is the charge density autocorrelation function, a two-particle property, related to the exchange-correlation hole. In the framework of a dipolar approximation of the Coulomb interaction around the molecular origin, one obtains the so-called Salem-Tang-Karplus approximation to the C(6) dispersion coefficient. Alternatively, by expanding the Coulomb interaction around the center of charge (centroid) of the exchange-correlation hole associated with each point in the molecular volume, a multicenter expansion is obtained around the centroids of electron localization domains, always in terms of the exchange-correlation hole. In order to get a formula analogous to that of Becke and Johnson, which involves the exchange-hole only, further assumptions are needed, related to the difficulties of obtaining the expectation value of a two-electron operator from a single determinant. Thus a connection could be established between the conventional fluctuating charge density model of London dispersion forces and the notion of the "exchange-hole dipole moment" shedding some light on the true nature of the approximations implicit in the Becke-Johnson model.     

Efficient evaluation of Coulomb integrals in a mixed Gaussian and plane-wave basis using the density fitting and Cholesky decomposition

The objective of this paper is to show that the density fitting (resolution of the identity approximation) can also be applied to Coulomb integrals of the type (k(1)(1)k(2)(1)|g(1)(2)g(2)(2)), where k and g symbols refer to plane-wave functions and gaussians, respectively. We have shown how to achieve the accuracy of these integrals that is needed in wave-function MO and density functional theory-type calculations using mixed Gaussian and plane-wave basis sets. The crucial issues for achieving such a high accuracy are application of constraints for conservation of the number electrons and components of the dipole moment, optimization of the auxiliary basis set, and elimination of round-off errors in the matrix inversion.     

Efficient evaluation of short-range Hartree-Fock exchange in large molecules and periodic systems

We present an efficient algorithm for the evaluation of short-range Hartree-Fock exchange energies and geometry gradients in Gaussian basis sets. Our method uses a hierarchy of screening levels to eliminate negligible two-electron integrals whose evaluation is the fundamental computational bottleneck of the procedure. By applying our screening technique to the Heyd-Scuseria-Ernzerhof [J. Chem. Phys. 118, 8207 (2003)] short-range Coulomb hybrid density functional, we achieve a computational efficiency comparable with that of standard nonhybrid density functional calculations.     

Significant reduction of on-site Coulomb energy U due to short-range correlation in an organic Mott insulator.

By carrying out a first-principles T-matrix calculation on multiple scatterings between electrons, we show that the intramolecular electron-electron interaction energy U, of a Mott insulator of the organic radical 1,3,5-trithia-2,4,6-triazapentalenyl (TTTA) is significantly reduced from the naive expectation value of the Coulomb interaction (7.3 eV and 5.3 eV, respectively, for the bare and screened Coulomb interactions) to 2.9 eV due to the short-range correlation. This result together with the intermolecular interaction energy D=1.3 eV explains the experimental optical gap (1.5 eV). The associated two-particle wavefunction clearly shows the Coulomb hole indicating that two electrons with antiparallel spins cannot approach because of the Coulomb repulsion. We also discuss the energetics and magnetics of this system.     

Fundamental importance of the Coulomb hole sum rule to the understanding of the Colle-Salvetti wave function functional

In this paper we consider the general form of the correlated-determinantal wave function functional of Colle and Salvetti (CS) for the He atom. The specific form employed by CS is the basis for the widely used CS correlation energy formula and the Lee-Yang-Parr correlation energy density functional of Kohn-Sham density functional theory. We show the following: (i) The key assumption of CS for the determination of this wave function functional, viz., that the resulting single-particle density matrix and the Hartree-Fock theory Dirac density matrix are the same, is equivalent to the satisfaction of the Coulomb hole sum rule for each electron position. The specific wave function functional derived by CS does not satisfy this sum rule for any electron position. (ii) Application of the theorem on the one-to-one correspondence between the Coulomb hole sum rule for each electron position and the constraint of normalization for approximate wave functions then proves that the wave function derived by CS violates charge conservation. (iii) Finally, employing the general form of the CS wave function functional, the exact satisfaction of the Coulomb hole sum rule at each electron position then leads to a wave function that is normalized. The structure of the resulting approximate Coulomb holes is reasonably accurate, reproducing both the short- and the long-range behavior of the hole for this atom. Thus, the satisfaction of the Coulomb hole sum rule by an approximate wave function is a necessary condition for constructing wave functions in which electron-electron repulsion is represented reasonably accurately.     

The Coulomb Hole of the Ne Atom

We analyze the Coulomb hole of Ne from highly-accurate CISD wave functions obtained from optimized even-tempered basis sets. Using a two-fold extrapolation procedure we obtain highly accurate results that recover 97 % of the correlation energy. We confirm the existence of a shoulder in the short-range region of the Coulomb hole of the Ne atom, which is due to an internal reorganization of the K-shell caused by electron correlation of the core electrons. The feature is very sensitive to the quality of the basis set in the core region and it is not exclusive to Ne, being also present in most of second-row atoms, thus confirming that it is due to K-shell correlation effects.     

The properties of off-shell coulomb radial wavefunctions

Approximations to exact wave functions for the scattering of few-particle systems often involve components corresponding to the interaction of two of the particles “off the energy-shell”. Several examples arising in the collision of ions and photons with atoms are given. An expansion in partial waves leads to an off-shell radial wave function. The defining differential equation is solved here numerically with particular emphasis on the behaviour arising from two-body potentials of long-range Coulomb form. The transition to shell of the radial wave functions, Jost functions and solutions andT-matrix elements is discussed for both short-range and Coulomb potentials. It is shown that the approximation of a Coulomb potential by a shorter-range form involves little error when sufficiently far off the energy shell.     

Unpaired electrons,spin polarization,and bond orders in radicals from the 2‐RDM in orbital spaces: Basic notions and testing calculations

Adopting the second‐order reduced density matrix level, the conventional α‐ and β‐spin populations in radicals are split into paired and unpaired or electropon (referring to the simultaneous occurrence of an electron and a hole of opposite spins in an orbital) populations. This analysis gives the possibility to distinguish the (un)favorable for chemical bonding electronic interactions by means of positive or negative Coulomb and/or Fermi correlations of two electropons. To overcome the conceptual difficulties originated from the subtle superposition of unpaired electrons due to spin density and those responsible for chemical bonding, we use the notion of properly unpaired electrons. The quantity describing this notion provides a global picture for the ability of electrons of a given orbital to form covalent bonds with the electrons of all remaining orbitals. More detailed information, concerning the behavior of electrons in two distinct target orbitals, is obtained by means of the two‐electropon correlations. As shown, the boundary values of the used quantities are physically meaningful, and the whole theory is tested from various points of view concerning: localized and delocalized radical centers, orthogonal and nonorthogonal orbitals, uncorrelated and correlated levels, Coulomb and Fermi correlations. We also check the electropon based analysis by investigating the spin polarization effects and bond orders in radicals. The tests are achieved for well‐known radicals, and to preserve the stability of the numerical results and the invariance of the obtained conceptual pictures, we used natural basis sets introduced within the natural bond orbital methodology. © 2014 Wiley Periodicals, Inc.     

Kinetic Energy Density of the Two-Electron Hookean Atom in Terms of the Ground-State Electron Density

An explicit expression is derived for the kinetic energy density, including the correlation contribution, in terms of the ground-state electron density for the two-electron Hookean atom. This model atom has the merit that while the electrons are tied to an origin by springs, the Coulomb interaction energy between the two electrons is fully incorporated.     

Electron correlation effects in SCF+CI wave functions—Fermi and coulomb correlation holes in HCN

The Coulomb correlation hole distribution function has been computed with respect to various reference centers in the HCN molecule, using standard SCF +CI type wave functions. The extent to which statistical correlation between unlike-spin electrons is introduced into an SCF wave function through the inclusion of configuration interaction has been assessed by an examination of the range and depth of such holes, and compared with the behavior of analogous Fermi distribution functions. Our results show that the range of Fermi correlation is consistently longer than that of the corresponding Coulomb correlation.     

On the electron-electron counterbalance hole

When a many-electron system has spatial inversion symmetry, the electron-electron counterbalance hole implies that two electrons with parallel spins cannot be at opposite positions with respect to the inversion center, and its presence was pointed out in the literature [T. Koga, J. Chem. Phys. 108, 2515 (1998)] for any pairs of Hartree-Fock orbitals with the same inversion parity. We report here a generalized result that in all two-electron systems with spatial inversion symmetry, the electron-electron counterbalance hole always exists for any approximate and exact wave functions with even inversion parity. The same is also true in momentum space. An extension of the hole to systems with three or more electrons is discussed.     

Long-range energy transfer and ionization in extended quantum systems driven by ultrashort spatially shaped laser pulses

The processes of ionization and energy transfer in a quantum system composed of two distant H atoms with an initial internuclear separation of 100 atomic units (5.29 nm) have been studied by the numerical solution of the time-dependent Schr?dinger equation beyond the Born-Oppenheimer approximation. Thereby it has been assumed that only one of the two H atoms was excited by temporally and spatially shaped laser pulses at various laser carrier frequencies. The quantum dynamics of the extended H-H system, which was taken to be initially either in an unentangled or an entangled ground state, has been explored within a linear three-dimensional model, including the two z coordinates of the electrons and the internuclear distance R. An efficient energy transfer from the laser-excited H atom (atom A) to the other H atom (atom B) and the ionization of the latter have been found. It has been shown that the physical mechanisms of the energy transfer as well as of the ionization of atom B are the Coulomb attraction of the laser driven electron of atom A by the proton of atom B and a short-range Coulomb repulsion of the two electrons when their wave functions strongly overlap in the domain of atom B.     

Resonant processes and Coulomb interactions in (C59N)2

We have determined the on-site molecular Coulomb interaction energy U of a (C59N)2 bulk film and find values ranging from 1.10+/-0.10 eV for the highest occupied molecular orbital to 1.35+/-0.10 eV for the deeper lying orbitals, comparable to values found in C60. The on-site Coulomb interaction between a carbon core hole and valence electrons, Uc, is, however, substantially lower than in C60 at 1.35+/-0.07 eV. Resonant photoemission (RESPES) results show a weakened participator decay channel, especially around the N 1s threshold, where resonance of the highest occupied molecular orbital shoulder is absent. Near-edge x-ray absorption fine structure and constant initial state measurements, taken in parallel with the RESPES data, indicate, however, that matrix element effects cannot be ruled out.     

A short-range correlation energy density functional with multi-determinantal reference

We introduce a short-range correlation density functional defined with respect to a multi-determinantal reference which is meant to be used in a multi-determinantal extension of the Kohn–Sham scheme of density functional theory based on a long-range/short-range decomposition of the Coulomb electron–electron interaction. We construct the local density approximation for this functional and discuss its performance on the He atom.     

Review of theories on ionization in fast ion-atom collisions with prospects for applications to hadron therapy

This work reviews quantum-mechanical four-body distorted wave theories for double electron capture in collisions between fast heavy multiply charged ions and heliumlike atomic systems. The widely used distorted wave methods of the first- and second-order in the pertinent perturbation series expansions are compared with each other. This tests the presumed importance of double continuum intermediate states of two electrons. Further, the relative performance is evaluated of the second-order theories with and without the eikonalization of the two-electron Coulomb wave functions for double continuum intermediate states. This checks the correctness and usefulness of the eikonalized Coulomb waves when two electrons participate actively to the transition from the initial to the final state of the entire system. We also analyze the significance of the contributions from excited heliumlike states especially in comparison between theory and measurement. The overall goal of the present study is to determine how much of the unprecedented experience gained over several decades in studying high-energy theories of pure three-body charge exchange could be exported directly to four-body double-electron capture without much of additional and essential eleaborations, besides the naturally increased computational demand. In particular, we address the unexpected breakdown of the continuum distorted wave eikonal initial state approximation and the anticipated success of continuum distorted wave theory for double charge exchange in ion-atom collisions at high impact energies.     

Adiabatic approximation of time-dependent density matrix functional response theory

Time-dependent density matrix functional theory can be formulated in terms of coupled-perturbed response equations, in which a coupling matrix K(omega) features, analogous to the well-known time-dependent density functional theory (TDDFT) case. An adiabatic approximation is needed to solve these equations, but the adiabatic approximation is much more critical since there is not a good "zero order" as in TDDFT, in which the virtual-occupied Kohn-Sham orbital energy differences serve this purpose. We discuss a simple approximation proposed earlier which uses only results from static calculations, called the static approximation (SA), and show that it is deficient, since it leads to zero response of the natural orbital occupation numbers. This leads to wrong behavior in the omega-->0 limit. An improved adiabatic approximation (AA) is formulated. The two-electron system affords a derivation of exact coupled-perturbed equations for the density matrix response, permitting analytical comparison of the adiabatic approximation with the exact equations. For the two-electron system also, the exact density matrix functional (2-matrix in terms of 1-matrix) is known, enabling testing of the static and adiabatic approximations unobscured by approximations in the functional. The two-electron HeH(+) molecule shows that at the equilibrium distance, SA consistently underestimates the frequency-dependent polarizability alpha(omega), the adiabatic TDDFT overestimates alpha(omega), while AA improves upon SA and, indeed, AA produces the correct alpha(0). For stretched HeH(+), adiabatic density matrix functional theory corrects the too low first excitation energy and overpolarization of adiabatic TDDFT methods and exhibits excellent agreement with high-quality CCSD ("exact") results over a large omega range.     

Linear scaling computation of the Fock matrix. VII. Periodic density functional theory at the Gamma point

Linear scaling quantum chemical methods for density functional theory are extended to the condensed phase at the Gamma point. For the two-electron Coulomb matrix, this is achieved with a tree-code algorithm for fast Coulomb summation [M. Challacombe and E. Schwegler, J. Chem. Phys. 106, 5526 (1997)], together with multipole representation of the crystal field [M. Challacombe, C. White, and M. Head-Gordon, J. Chem. Phys. 107, 10131 (1997)]. A periodic version of the hierarchical cubature algorithm [M. Challacombe, J. Chem. Phys. 113, 10037 (2000)], which builds a telescoping adaptive grid for numerical integration of the exchange-correlation matrix, is shown to be efficient when the problem is posed as integration over the unit cell. Commonalities between the Coulomb and exchange-correlation algorithms are discussed, with an emphasis on achieving linear scaling through the use of modern data structures. With these developments, convergence of the Gamma-point supercell approximation to the k-space integration limit is demonstrated for MgO and NaCl. Linear scaling construction of the Fockian and control of error is demonstrated for RBLYP6-21G* diamond up to 512 atoms.