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.     

New link between conceptual density functional theory and electron delocalization

In this paper we give a new definition of the softness kernel based on the exchange-correlation density. This new kernel is shown to correspond to the change of electron fluctuation upon external perturbation, thus helping to bridge the gap between conceptual density functional theory and some tools describing electron localization in molecules. With the aid of a few computational calculations on diatomics we illustrate the performance of this new computational tool.     

Generalized-gradient exchange-correlation hole obtained from a correlation factor ansatz

The Perdew-Burke-Ernzerhof (PBE) approximation to the exchange-correlation energy is employed as reference point for the construction of an angle-averaged exchange-correlation hole. First, we develop a new model for the PBE exchange hole. In contrast to the previous model [Ernzerhof and Perdew, J. Chem. Phys. 109, 3313 (1998)], it contains an atomic exchange hole, similar to the Becke-Roussel model [Becke and Roussel, Phys. Rev. A 39, 3761 (1989)]. A correlation factor, i.e., a function multiplying the exchange hole, is proposed that turns the exchange into an exchange-correlation hole. The correlation factor has a simple form and is determined through a number of known conditions that should be satisfied by a generalized-gradient exchange-correlation hole.     

Interatomic exchange-correlation interaction energy from a measure of quantum theory of atoms in molecules topological bonding: A diatomic case

The potential relations between the measure of topological interatomic bonding—integrals of electron density with respect to internuclear axis over the corresponding quantum theory of atoms in molecules (QTAIM)-defined interatomic surface (IAS)—and interatomic exchange-correlation contributions from the interacting quantum atoms approach are discussed. The quantum chemical computations of 38 equilibrium diatomic systems at different levels of theory (HF, MP2, MP4SDQ, and CCSD) are invoked to support abstract considerations. Parameters of excellent correlations between IAS integrals and interatomic exchange-correlation energy are found by the optimization. The performance of these trends depends on the accuracy of the electronic correlation treatment. The resulting trends are a unique feature of equilibrium states, whereas more complicated dependencies are explored for several systems at non-equilibrium conditions. The relations of established trends with other IAS-based estimations of strength of bonding interactions between topological atoms and issues explored for multiatomic systems are briefly discussed.     

Electron pairing and chemical bonds. Electron fluctuation and pair localization in ELF domains

This article reports the numerical comparison of the quantities characterizing the extent of electron fluctuation and pair localization in the domains determined by the direct minimization of electron fluctuation with the domains resulting from the partitioning of the molecules based on the topological analysis of the so-called electron localization function (ELF). Such a comparison demonstrates that the ELF partitioning can be regarded as a feasible alternative to computationally much more demanding direct optimization of minimum fluctuation domains. This opened the possibility of the systematic scrutiny of the electron pair model of the chemical bond, and as it was demonstrated, the previous pessimistic claims about the applicability of this model are not completely justified.     

Real-space analysis of the exchange-correlation energy

The exchange-correlation energy of a many-electron system may be written as the electrostatic interaction between the electron density at position r and the density of the exchange-correlation hole at position r + u. If we average the hole over the entire system, we find that the energy is uniquely decomposed into contributions from various electronic separations u. We may also decompose the hole into contributions from parallel and antiparallel spins. We give several exact conditions which this system-averaged, spindecomposed exchange-correlation hole satisfies. Local spin density (LSD ) and generalized gradient approximations (GGAS ), are more appropriate for u → 0 than for large u and more trustworthy for antiparallel spins than for parallel spins. We illustrate how good LSD is as u = 0 with explicit examples, but also note that, contrary to expectation, LSD is not exact for u=0, except in certain limiting cases. We show that the dramatic failure of the second-order gradient expansion for large u can be cured by a real-space cutoff procedure which generates a nonempirical GGA, the Pw91 functional. We conclude with some thoughts about the search for greater accuracy in the next 30 years of density functional theory. © 1995 John Wiley & Sons, Inc.     

Electron delocalization and aromaticity variations in the stacked nucleic acid base pairs

A newly developed exchange-correlation functional (MPWB1K) in density functional theory has been applied to evaluate the electron delocalization of individual fragments in the stacking interaction between nucleic acid bases (NABs). Electronically and structural-based indices have been employed to investigate the aromaticity variation during stacking interaction. A quantitative study of NABs in their isolated and stacked forms reveals that stacking interaction causes a decrease in bond delocalization. It is shown that the decrease in the aromaticity is accompanied by local decrease in two-center delocalization indices within the pyrimidine rings. We found that the aromaticity exhibits a similar trend for NABs in both their isolated and stacked forms. Moreover, it is indicated that aromatic fluctuation index is more sensible index to delineate the aromaticity variation during stacking interaction.

     

Orbital energies and negative electron affinities from density functional theory: Insight from the integer discontinuity

Orbital energies in Kohn-Sham density functional theory (DFT) are investigated, paying attention to the role of the integer discontinuity in the exact exchange-correlation potential. A series of closed-shell molecules are considered, comprising some that vertically bind an excess electron and others that do not. High-level ab initio electron densities are used to calculate accurate orbital energy differences, Deltavarepsilon, between the lowest unoccupied molecular orbital (LUMO) and the highest occupied molecular orbital (HOMO), using the same potential for both. They are combined with accurate vertical ionization potentials, I(0), and electron affinities, A(0), to determine accurate "average" orbital energies. These are the orbital energies associated with an exchange-correlation potential that averages over a constant jump in the accurate potential, of magnitude Delta(XC)=(I(0)-A(0))-Deltavarepsilon, as given by the discontinuity analysis. Local functional HOMO energies are shown to be almost an order of magnitude closer to these average values than to -I(0), with typical discrepancies of just 0.02 a.u. For systems that do not bind an excess electron, this level of agreement is only achieved when A(0) is set equal to the negative experimental affinity from electron transmission spectroscopy (ETS); it degrades notably when the zero ground state affinity is instead used. Analogous observations are made for the local functional LUMO energies, although the need to use the ETS affinities is less pronounced for systems where the ETS values are very negative. The application of an asymptotic correction recovers the preference, leading to positive LUMO energies (but bound orbitals) for these systems, consistent with the behavior of the average energies. The asymptotically corrected LUMO energies typically agree with the average values to within 0.02 a.u., comparable to that observed with the HOMOs. The study provides numerical support for the view that local functionals exhibit a near-average behavior based on a constant jump of magnitude Delta(XC). It illustrates why a recently proposed DFT expression involving local functional frontier orbital energies and ionization potential yields reasonable estimates of negative ETS affinities and is consistent with earlier work on the failure of DFT for charge-transfer excited states. The near-average behavior of the exchange-correlation potential is explicitly illustrated for selected systems. The nature of hybrid functional orbital energies is also mentioned, and the results of the study are discussed in terms of the variation in electronic energy as a function of electron number. The nature of DFT orbital energies is of great importance in chemistry; this study contributes to the understanding of these quantities.     

Asymptotic correction of the exchange-correlation kernel of time-dependent density functional theory for long-range charge-transfer excitations

Time-dependent density functional theory (TDDFT) calculations of charge-transfer excitation energies omegaCT are significantly in error when the adiabatic local density approximation (ALDA) is employed for the exchange-correlation kernel fxc. We relate the error to the physical meaning of the orbital energy of the Kohn-Sham lowest unoccupied molecular orbital (LUMO). The LUMO orbital energy in Kohn-Sham DFT--in contrast to the Hartree-Fock model--approximates an excited electron, which is correct for excitations in compact molecules. In CT transitions the energy of the LUMO of the acceptor molecule should instead describe an added electron, i.e., approximate the electron affinity. To obtain a contribution that compensates for the difference, a specific divergence of fxc is required in rigorous TDDFT, and a suitable asymptotically correct form of the kernel fxc(asymp) is proposed. The importance of the asymptotic correction of fxc is demonstrated with the calculation of omegaCT(R) for the prototype diatomic system HeBe at various separations R(He-Be). The TDDFT-ALDA curve omegaCT(R) roughly resembles the benchmark ab initio curve omegaCT CISD(R) of a configuration interaction calculation with single and double excitations in the region R=1-1.5 A, where a sizable He-Be interaction exists, but exhibits the wrong behavior omegaCT(R)相似文献   

Derivative discontinuity, bandgap and lowest unoccupied molecular orbital in density functional theory

The conventional analysis of Perdew and Levy, and Sham and Schlu?ter shows that the functional derivative discontinuity of the exchange-correlation density functional plays a critical role in the correct prediction of bandgaps, or the chemical hardness. In a recent work by the present authors, explicit expressions for bandgap prediction with all common types of exchange-correlation functionals have been derived without invoking the concept of exchange-correlation energy functional derivative discontinuity at all. We here analyze the two approaches and establish their connection and difference. The present analysis further leads to several important results: (1) The lowest unoccupied molecular orbital (LUMO) in DFT has as much meaning in describing electron addition as the highest occupied molecular orbital (HOMO) in describing electron removal. (2) Every term in the total energy functional contributes to the energy gap because of the discontinuity of the derivative of the density (or density matrix) with respect to the number of electrons, ((?ρ(s)(r('),r))/?N)(v(s) ), at integers. (3) Consistent with the Perdew-Levy-Sham-Schlu?ter conclusion that the exact Kohn-Sham energy gap differs from the fundamental bandgap by a finite correction due to the functional derivative discontinuity of the exchange-correlation energy, we show that the exchange-correlation functional cannot be an explicit and differentiable functional of the electron density, either local or nonlocal. The last result is further strengthened when we consider Mott insulators. There, the exact exchange-correlation functional needs to have an explicitly discontinuous (nondifferentiable) dependence on the density or the density matrix. (4) We obtain exact conditions on the derivatives of total energy with respect to the spin-up and spin-down number of electrons.     

Electron density, exchange-correlation density, and bond characterization from the perspective of the valence-bond theory. I. Two simple analytical cases

In this work, using a valence-bond wave function we obtain analytical expressions for the first- and second-order reduced density matrices of two simple, but quite representative, cases of diatomic molecular systems, namely, H2 and LiH. A detailed study of their exchange-correlation density is performed for both equilibrium and nonequilibrium internuclear distances, discriminating the parallel- and antiparallel-spin contributions. The results show that the behavior of the exchange-correlation density clearly changes with the character of the bond, making it possible to obtain a good deal of information regarding the type of the bond interaction.     

Diminished gradient dependence of density functionals: constraint satisfaction and self-interaction correction

The Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation for the exchange-correlation energy functional has two nonempirical constructions, based on satisfaction of universal exact constraints on the hole density or on the energy. We show here that, by identifying one possible free parameter in exchange and a second in correlation, we can continue to satisfy these constraints while diminishing the gradient dependence almost to zero (i.e., almost recovering the local spin density approximation or LSDA). This points out the important role played by the Perdew-Wang 1991 nonempirical hole construction in shaping PBE and later constructions. Only the undiminished PBE is good for atoms and molecules, for reasons we present, but a somewhat diminished PBE could be useful for solids; in particular, the surface energies of solids could be improved. Even for atoms and molecules, a strongly diminished PBE works well when combined with a scaled-down self-interaction correction (although perhaps not significantly better than LSDA). This shows that the undiminished gradient dependence of PBE and related functionals works somewhat like a scaled-down self-interaction correction to LSDA.     

Heteronuclear diatomic force constants clarified through perturbation theory II

A simple relationship between the heteronuclear diatomic force constant (K(AB)) and the homonuclear diatomic force constants (K(AA), K(BB)), which was proposed in a previous report, has been improved through the second-order perturbation theory as K(AB) = zeta3(K(AA) x K(BB))(1/2); zeta = (R(AA) x R(BB))(1/2)/R(AB) where zeta denotes the correction factor in which R(AB), R(AA), and R(BB) are the equilibrium internuclear distances of diatomic molecules AB, AA, and BB, respectively. To test the above expression, a large number of heteronuclear diatomic force constants have been calculated and compared with those obtained from normal coordinate analyses as well as ab initio quantum mechanical methods (Gaussian 98W). We have found that the above modified expression better reproduces the force constants of most heteronuclear diatomic molecules than the previous expression. It is therefore expected that the expression may also be applied to the prediction of stretching force constants between heteronuclear diatomics in various polyatomic molecules.     

Ionization potentials and electron affinities in the Perdew-Zunger self-interaction corrected density-functional theory

Using a self-consistent implementation of the Perdew-Zunger self-interaction corrected (PZ-SIC) density-functional theory, we have calculated ionization potentials (IP) and electron affinities (EA) of first- and second-row atoms and a set of small molecules. Several exchange-correlation functionals were tested. IPs and EAs were obtained by two methods: as the difference in self-consistent field (SCF) energies of neutrals and ions (deltaSCF) and as negatives of highest-occupied orbital energies. We found that, except for local spin-density approximation, PZ-SIC worsens DeltaSCF IPs and EAs. On the other hand, PZ-SIC brings orbital eigenvalues into much better agreement with electron removal energies. The Perdew-Zunger SIC seems to over-correct many-electron systems; for molecules it performs worse than for atoms. We also discuss several common approximations to PZ-SIC such as spherical averaging of orbital densities in atoms.     

Electronic spectroscopy of UO2(2+), NUO(+) and NUN: an evaluation of time-dependent density functional theory for actinides

The performance of the time-dependent density functional theory (TDDFT) approach has been evaluated for the electronic spectrum of the UO(2)(2+), NUO(+) and NUN molecules. Different exchange-correlation functionals (LDA, PBE, BLYP, B3LYP, PBE0, M06, M06-L, M06-2X, CAM-B3LYP) and the SAOP model potential have been investigated, as has the relative importance of the adiabatic local density approximation (ALDA) to the exchange-correlation kernel. The vertical excitation energies have been compared with reference data obtained using accurate wave-function theory (WFT) methods.     

Some theoretical aspects of vibrational fine structure accompanying core ionizations in N2 and CO

Ab initio calculations have been carried out on CO and N₂ and the relevant core hole states with different basis sets to investigate differences in geometries and force constants. From these calculations vibrational band profiles of the core level ESCA spectra for these molecules have been interpreted, obviating the need to rely on data pertaining to the equivalent core species. The agreement with experimental profiles is excellent. The O_1s level of CO which has not been subjected to detailed theoretical analysis previously, is predicted to show substantial vibrational structure in excellent agreement with recently acquired experimental data. The effect of temperature on the band profiles has also been considered. Theoretically derived core binding and relaxation energies of these systems have been investigated both as a function of basis set, and of internuclear distance. Density difference contours have been computed and give a straightforward pictorial representation of the substantial electron reorganizations accompanying core ionizations. Small basis sets with valence exponents appropriate to the equivalent core species when used in hole state calculations describe bond lengths, force constants, core binding energies and relaxation energies with an accuracy comparable to that appropriate to the corresponding extended basis set calculations.