Effect of modified virtual orbitals on local Coulomb correlation holes in FCN

Local Coulomb correlation hole distribution functions may be used to assess the extent to which electron correlation effects are present in large scale SCF + CI wave functions. From a set of modified virtual orbitals, ordered according to their interaction with the SCF configuration, we have constructed a limited SCF + CI wave function with improved convergence characteristics with respect to that formed from the canonical virtual orbital set. These wave functions, of the same size yet with different energies, have been used to examine the range and depth of local Coulomb correlation holes in FCN. In all cases, the depth of the local Coulomb hole is no more than 10% or so of that of the corresponding Fermi hole, and the range Fermi correlation is generally less than that of Fermi correlation. This is particularly marked in the high density regions around the nuclei. The significance of our results is discussed in relation to a recent proposal for the incorporation of Coulomb correlation into the local exchange method.     

Statistical electron correlation — coefficients and — holes in molecules

Summary Electron correlation in the H₂, LiH and BH molecules has been analyzed in terms of the statistical correlation coefficients introduced by Kutzelnigg, Del Re, and Berthier. Angular, radial (in-out), longitudinal (left-right) and transverse correlation coefficients have been evaluated from both self-consistent-field (SCF) and configuration interaction (CI) wave functions. It has been found that these coefficients reflect fairly well the correlation behavior in the molecular system. The lack of spherical symmetry in molecular densities adds new features to these correlation coefficients and this information can be useful for the study of electronic structure in molecules. The correlation hole function, Fermi and Coulomb holes in these systems have also been calculated and discussed.Dedicated to Professor Werner Kutzelnigg on the occasion of his sixtieth birthday     

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     

Role of electron correlation in the quantum-mechanical calculations of the coulomb interaction energy in the DNA base pairs

The intermolecular electronic correlation contributions to the Coulomb component of the nucleic acid base interaction energy are estimated. The Coulomb energy is evaluated with the use of atomic monopoles, which are determined from the π-electronic densities calculated by the SCF method and by employing partially or completely optimized APSG wave functions. When the correlation is thus taken into account, a systematic decrease in atomic charges occurs; this effect is considerable only if an optimized orbital set is used. As a result, the Coulomb interaction energy due to the π-electronic atoms decreases from ?1.13 to ?0.85 kcal/mol for the AT pair and from ?7.15 to ?4.61 kcal/mol for the GC pair.     

Electron correlation: the many-body problem at the heart of chemistry

The physical interactions among electrons and nuclei, responsible for the chemistry of atoms and molecules, is well described by quantum mechanics and chemistry is therefore fully described by the solutions of the Schr?dinger equation. In all but the simplest systems we must be content with approximate solutions, the principal difficulty being the treatment of the correlation between the motions of the many electrons, arising from their mutual repulsion. This article aims to provide a clear understanding of the physical concept of electron correlation and the modern methods used for its approximation. Using helium as a simple case study and beginning with an uncorrelated orbital picture of electronic motion, we first introduce Fermi correlation, arising from the symmetry requirements of the exact wave function, and then consider the Coulomb correlation arising from the mutual Coulomb repulsion between the electrons. Finally, we briefly discuss the general treatment of electron correlation in modern electronic-structure theory, focussing on the Hartree-Fock and coupled-cluster methods and addressing static and dynamical Coulomb correlation.     

Angular aspects of exchange correlation and the fermi hole

The complete (nonreduced) αα probability density functions evaluated from the Hartree–Fock and simple Hartree product wavefunctions have been used to elucidate the angular features of spin correlation and the Fermi hole in the 2³S state of helium and the ground state of beryllium. This approach shows that the local Fermi holes in these two cases are very similar and that the Fermi hole is essentially spherically symmetric when the reference electron is close to the nucleus. As the reference electron is removed to larger radial distances, appreciable polarization of the Fermi hole is observed. The polarization is greater in the direction of the nucleus than away from the nucleus, contrary to the situation in the Coulomb hole of the helium ground state where the polarization is greater away from the nucleus than toward the nucleus. Several other differences between the He 2³S Fermi hole and the He 1¹S Coulomb hole are noted.     

Small-angle electron scattering by formaldehyde and ketene: Effects of electron correlation and chemical binding

The total (elastic plus inelastic) intensities of 51 keV electrons scattered by H₂CO and H₂CCO have been measured over a range of K = (4π/λ) sin(θ/2) = 1–9.5 Å^?1 and compared with the theoretical intensities calculated with SCF and CI wave functions. Significant discrepancies are found between the experimental intensities and the theoretical ones based on the SCF wave functions. Most of the chemical binding and electron correlation effects observed in the total scattered intensities are reproduced by the theoretical intensities based on the CI wave functions calculated with the basis set including polarization functions on all atoms. © 1992 John Wiley & Sons, Inc.     

Fermi and Coulomb holes of molecules in excited states

Correlation holes of electrons with the same (Fermi hole) and different (Coulomb hole) spins in the ground (X¹Σ⁺), first (A¹Σ⁺) and second (B¹II) excited states of LiH were constructed from full configuration interaction (CI ) wave functions. It was found that the shapes of both the Fermi and Coulomb holes in these states are dependent on the location of the reference electron. When the reference electron is chosen to be close to the Li nucleus, the Fermi correlation results in a large negative hole for all three states. However, the A¹Σ⁺ excited state is further characterized by displaying a second hole around the H nucleus, and in the B¹II state, the hole is elongated along the molecular axis. Coulomb correlation shows up strongly in the A¹Σ⁺ state and, in addition, there is clearly correlation of electrons at the two nuclei. These features of the correlation holes were compared with those from a two-Slater-determinant model wave function. The Hartree, Fermi, and Coulomb screening potentials in these states were also studied in the light of possible modeling of the correlation functionals for the excited states. © 1995 John Wiley & Sons, Inc.     

A new look at correlations in atomic and molecular systems. I. Application of fermion monte carlo variational method

We are engaged in research directed toward the development of compact and accurate correlation functions for many-electron systems. Our computational tool is the variational method in which the many-electron integrals are calculated by Monte Carlo using the fermion Metropolis sampling algorithm. That is, a many-fermion system is simulated by sampling the square of a correlated antisymmetric wave function. The principal advantage of the method is that interelectronic distance r_ij may be included directly in the wave function without adding significant computational complexity. In addition, other quantities of physical and theoretical interest such as electron correlation functions and representations of Coulomb and Fermi “holes” are very easily obtained. Preliminary results are reported for He, H₂, and Li₂.     

The calculation of frequency-dependent hyperpolarizabilities including electron correlation effects

The theory for the calculation of the frequency-dependent hyperpolarizabilities β(?2ω; 0, ω), β(?ω; 0, ω), and β(0; ω, ?ω) is discussed. New relations between these tensors are derived for those wave functions that obey the time-dependent Hellmann–Feynman theorem (e.g., the self-consistent field [SCF] or the exact wave function). Using second-order Møller–Plesset perturbation theory (MP2), expressions are obtained for the hyperpolarizabilities in terms of derivatives of appropriately defined linear polarizability tensors with respect to a static electric field. Results are presented for ammonia and formaldehyde for the optical Kerr effect and for secondharmonic generation. These results indicate that it is desirable to determine the frequency-dependent contribution to the hyperopolarizability at the MP2 rather than the SCF level of theory, in cases where the static hyperpolarizability has a large contribution from electron correlation and/or where the frequency-dependent contribution may be more significant, such as for secondharmonic generation.     

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.     

On the SCF calculation of excited states: Singlet states in the two-electron problem

The problem of determining SCF wave functions for excited electronic states is examined for singlet states of two-electron systems using a Lowdin natural orbital transformation of the full CI wave function. This analysis facilitates the comparison of various SCF methods with one another. The distribution of the full CI states among the natural orbital MCSCF states is obtained for the S states of helium using a modest Gaussian basis set. For SCF methods that are not equivalent to the full CI wave functions, it is shown that the Hartree-Fock plus all single excitation wave functions are equivalent to that of Hartree-Fock plus one single excitation. It is further shown that these wave functions are equivalent to the perfect pair or TCSCF wave functions in which the CI expansion coefficients are restricted to have opposite signs. The case of the natural orbital MCSCF wave function for two orbitals is examined in greater detail. It is shown that the first excited state must always be found on the lower natural orbital MCSCF CI root, thus precluding the use of the Hylleras-Undeim-MacDonald (HUM) theorem in locating this state. It is finally demonstrated that the solution obtained by applying the HUM theorem (minimizing the upper MCSCF CI root with respect to orbital mixing parameters) is an artifact of the MCSCF method and does not correspond to any of the full CI states.     

Visualization of electron correlation in a series of helium S states

Electron correlation has been studied for a series of helium S states represented by a variety of wavefunctions, the best of which are accurate Hylleraas—Kinoshita functions. The states studied are the ground state, the lowest excited ¹S and ³S states, and the (2s)² and (2p)² doubly-excited ¹S states. Primary data is obtained from graphs of the conditional probability density as a function of the radial distance r₂ and the interelectronic angle θ₁₂, given that r₁ is fixed at various distances. Such graphs make clear the extent to which characteristics such as angular and radial correlations, and Fermi and Coulomb holes, are consequences of the relative motion of electrons in two-electron atoms.     

Accurate predictions of the energetics of silicon compounds using the multireference correlation consistent composite approach

Theoretical studies, using the multireference correlation consistent composite approach (MR-ccCA), have been carried out on the ground and lowest lying spin-forbidden excited states of a series of silicon-containing systems. The MR-ccCA method is the multireference equivalent of the successful single reference ccCA method that has been shown to produce chemically accurate (within ±1.0 kcal mol(-1) of reliable, well-established experiment) results. The percentage contributions of the SCF configurations to complete active space self-consistent field wave functions together with the Frobenius norm of the t(1) vectors and related D(1) diagnostics of the coupled-cluster single double wave function with the cc-pVTZ basis set have been utilized to illustrate the multi-configurational characteristics of the compounds considered. MR-ccCA incorporates additive terms to account for relativistic effects, atomic spin-orbit coupling, scalar relativistic effects, and core-valence correlation. MR-ccCA has been utilized to predict the atomization energies, enthalpies of formation, and the lowest energy spin-forbidden transitions for Si(n)X(m) (2 ≤ n + m ≥ 3 where n ≠ 0 and X = B, C, N, Al, P), silicon hydrides, and analogous compounds of carbon. The energetics of small silicon aluminides and phosphorides are predicted for the first time.     

Characterization of bond and current correlation structures in a molecular many-electron wave function

A method has been developed to analyzed the bond and current correlation structures of a molecular many-electron wave function. It is shown that the second order density matrix contains information about the bond and current correlations in its off-diagonal components with respect to the indices of orbital basis functions. We break down the off-diagonal correlation functions into five kinds: charge, spin scalar, spin quadrupole, charge spin, and spin polar correlation functions. For a real wave function, the four correlation functions, except for the spin polar one, have only symmetric–symmetric and antisymmetric–antisymmetric components. The former components give site–bond and bond–bond correlations of charges and spins, while the latter components give current–current correlations of charges and spins. The spin polar correlation function has only symmetric–antisymmetric components that give site–current and bond–current correlations of spins. The five off-diagonal correlation functions are expressed in terms of the off-diagonal components of the second order density matrix. The linked off-diagonal correlation functions are defined in that they give dynamical bond and current correlations. The method is applied to the analyses of the bond and current correlations in the low lying exact eigenstates of the PPP Hamiltonian of benzene.     

Electron correlation in pericyclic reactivity: A similarity approach

The recently proposed topological approach to chemical reactivity in terms of the secondorder similarity index was applied to the detailed analysis of correlation effects in pericyclic reactivity. This analysis implies that (i) the electron correlation is generally more important in forbidden reactions than in the allowed ones; (ii) in contributing to the discrimination between the allowed and forbidden reactions, the Fermi correlation acts in parallel with the Coulomb one; and (iii) the so-called multibond reactions are more sensitive to electron correlation than are the electrocyclic ones.     

Electronic correlation studies. III. Self-correlated field method. Application to 2S ground state and 2P excited state of three-electron atomic systems

The self-correlated field method is based on the insertion in the group product wave function of pair functions built upon a set of correlated “local” functions and of “nonlocal” functions. This work is an application to three-electron systems. The effects of the outer electron on the inner pair are studied. The total electronic energy and some intermediary results such as pair energies, Coulomb and exchange “correlated” integrals, are given. The results are always better than those given by conventional SCF computations and reach the same level of accuracy as those given by more laborious methods used in correlation studies.