Extension of the Coulomb–Hole–Hartree–Fock theory to molecules

The Coulomb–Hole–Hartree–Fock method introduced by E. Clementi in the early 1960s and reparametrized more recently by S. Chakraworty and E. Clementi to compute the correlated electronic energy in atomic systems, is here extended to compute molecules. The new parametrization is obtained empirically by fitting first and second atomic ionization potentials from He to Ca and a few diatomic molecules. The present formulation makes use of either one or more determinants in order to ensure proper dissociation products, following the early proposal of G.C. Lie and E. Clementi in the context of density functional computations for molecular systems. The new formulation is tested against the dissociation energies of a large number of molecules and it is found satisfactory. © 1995 John Wiley & Sons, Inc.     

Coulomb-hole-Hartree-Fock functional for molecular systems

We have extended to molecules a density functional previously parametrized for atomic computations. The Coulombhole-Hartree-Fock functional, introduced by Clementi in 1963, estimated the dynamic correlation energy by the computation of a Hartree-Fock type single-determinant wavefunction, where the Hartree-Fock potential was augmented with an effective potential term, related to a hard Coulomb hole enclosing each electron. The method was later revised by S. Chakravorty and E. Clementi, Phys. Rev. A, 38 (1989) 2290, so that a Yukawa-type soft Coulomb hole replaced the previous hard hole. Atomic correlation energies, computed for atoms with Z = 2 to 54, as well as for a number of excited states, validated the method. In this work we have parametrized for molecules a function which controls the width of the soft Coulomb hole by fitting the first and second atomic ionization potentials of atoms with 1 Z 18 and the binding energies of a few diatomic molecules. The parametrization was successfully validated by computing the dissociation energy for a number of molecules. A few-determinant version of the Coulomb-Hartree-Fock method (CHF-N) necessary to account for the non-dynamic correlation correction and to ensure proper dissociation products, is briefly discussed with reference to a previous proposal by G.C. Lie and E. Clementi, J. Chem. Phys., 60 (1974) 1275 and 60 (1974) 1288.     

Adapted Gaussian basis sets for atoms Cs to Lr based on the generator coordinate Hartree–Fock method

We have applied a discretized version of the generator coordinate Hartree–Fock method to generate adapted Gaussian basis sets for atoms Cs (Z=55) to Lr (Z=103). Our Hartree–Fock total energy results, for all atoms studied, are better than the corresponding Hartree–Fock energy results attained with previous Gaussian basis sets. For the atoms Cs to Lr we have obtained an energy value within the accuracy of 10⁻⁴ to 10⁻³ hartree when compared with the corresponding numerical Hartree–Fock total energy results. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 858–865, 1998     

The perturbation theory of the extended Hartree–Fock approximation for two-electron atoms

The first-order 1/Z perturbation theory of the extended Hartree–Fock approximation for two-electron atoms is described. A number of unexpected features emerge: (a) it is proved that the orbitals must be expanded in powers of Z^?1/2, rather than in Z^?1 as expected; (b) it is shown that the restricted Hartree–Fock and correlation parts of the orbitals can be uncoupled to first order, so that second-order energies are additive; (c) the equation describing the first-order correlation orbital has an infinite number of solutions of all angular symmetries in general, rather than only one of a single symmetry as expected; (d) the first-order correlation equation is a homogeneous linear eigenvalue-type equation with a non-local potential. It involves a parameter μ and an eigenvalue ω(μ) which may be interpreted as the probability amplitude and energy of a virtual correlation state. The second-order correlation energy is 2μ²ω. Numerical solutions for the first-order correlation orbitals, obtained variationally, are presented. The approximate second-order correlation energy is nearly 90% of the exact value. The first-order 1/Z perturbation theory of the natural-orbital expansion is described, and the coupled first-order integro-differential perturbation equations are obtained. The close relationship between the first-order extended Hartree–Fock correlation orbitals and the first-order natural correlation orbitals is discussed. A comparison of the numerical results with those of Kutzelnigg confirms the similarity.     

Soft Coulomb hole method applied to molecules

The soft Coulomb hole method introduces a perturbation operator, defined by ?e/r₁₂ to take into account electron correlation effects, where ω represents the width of the Coulomb hole. A new parametrization for the soft Coulomb hole operator is presented with the purpose of obtaining better molecular geometries than those resulting from Hartree–Fock calculations, as well as correlation energies. The 12 parameters included in ω were determined for a reference set of 12 molecules and applied to a large set of molecules (38 homo‐ and heteronuclear diatomic molecules, and 37 small and medium‐size molecules). For these systems, the optimized geometries were compared with experimental values; correlation energies were compared with results of the MP2, B3LYP, and Gaussian 3 approach. On average, molecular geometries are better than the Hartree–Fock values, and correlation energies yield results halfway between MP2 and B3LYP. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Numerical Hartree–Fock results for atoms Cs through Lr

Accurate nonrelativistic numerical Hartree–Fock results are reported for the heavy atoms Cs (Z = 55) through Lr (Z = 103) in their ground states. © 1995 John Wiley & Sons, Inc.     

Examination of several density functionals in numerical Kohn–Sham calculations for atoms

We studied several exchange‐only and exchange–correlation energy density functionals in numerical, i.e., basis‐set‐free, nonrelativistic Kohn–Sham calculations for closed‐shell ¹S states of atoms and atomic ions with N electrons, where 2≤N≤120. Accurate total energies are presented to serve as reference data for algebraic approaches, as do the numerical Hartree–Fock results, which are also provided. Gradient‐corrected exchange‐only functionals considerably improve the total energies obtained from the usual local density approximation, when compared to the Hartree–Fock results. Such an improvement due to gradient corrections is not seen in general for highest orbital energies, neither for exchange‐only results (to be compared with Hartree–Fock results), nor for exchange–correlation results (to be compared with experimental ionization energies). © 2001 John Wiley & Sons, Inc. Int J Quant Chem 82: 227–241, 2001     

Relativistic theory of an inhomogeneous electron liquid in relation to atomic binding energies

Experimentally determined ionization potentials in the literature are used to plot the binding energies for neutral atoms as a function of atomic number Z for Z?=?2–30, 32, 36, 42. From this pretty smooth plot we have subtracted non-relativistic Hartree–Fock binding energies, using both available numerical values and the almost analytical result, based on the non-relativistic Thomas–Fermi statistical theory valid for large Z. The difference is still relatively smooth. For Mo, with Z?=?42, the difference is about 70 atomic units. This difference is then analyzed using first relativistic theory of an inhomogeneous electron liquid and then the Local Density Approximation (LDA), and for Mo their results yield approximately 88 and 67 atomic units respectively. We infer that a highly accurate relativistic many-electron theory will therefore be needed before reliable electron correlation energies can be extracted from the experimental binding energies for atoms heavier than Argon. This fact has prompted us to use available LDA calculations to confront three theoretical predictions of the Z dependence of non-relativistic electron correlation energies at large Z.     

Near-Dirac–Hartree–Fock results for first-row atoms calculated with GTO basis sets

Relativistic basis sets for first-row atoms have been constructed by using the near-Hartree–Fock (nonrelativistic) eigenvectors calculated by Partridge. These bases generate results of near-Dirac–Hartree–Fock quality. Relativistic total and orbital energies, relativistic corrections to the total energy, and magnetic interaction energies for the first-row atoms have been presented. The smallest Gaussian expansions (13s8 p expansions) yield Dirac–Hartree–Fock total energies accurate through six significant digits, while the largest expansions (18s13p expansions) give these energies accurate through seven significant digits. These results are more accurate than some of the results reported earlier, particularly for the open-shell atoms, indicating that the basis employed is reasonably economical for relativistic calculations. © 1995 John Wiley & Sons, Inc.     

Analytical density-dependent representation of Hartree—Fock atomic potentials

A procedure to represent atomic electron charge densities [L. Fernandez Pacios, J. Phys. Chem., 95 , 10653 (1991); J. Phys. Chem., 96 , 7294 (1992)] is here generalized to obtain simple analytical functions for potential energy contributions. Based upon suitable functions to describe atomic electron densities in a physically meaningful form, the procedure is developed to define density-dependent analytical expressions for the electrostatic (classical) and exchange (quantum) potentials by means of proper approximate functionals. Calculations of correlation energies by using various density-functional approaches are also performed. The whole scheme is used to represent Hartree–Fock limit atomic wave functions by Clementi–Roetti. This way, a set of analytically simple, nonbasis set-dependent functions are defined with the aim to be further implemented in energy decomposition schemes for molecular interactions studies using atomic instead of electronic building blocks. © 1993 John Wiley & Sons, Inc.     

Atomic Xα calculations based on EHFS(α) = Eexp

The total energies and one-electron energies for first- and second-row atoms were calculated by using the Hartree–Fock and the Hartree–Fock-Slater Hamiltonian with X_α orbitals, u_i(α_exp); α was parametrized from E^HFS (α_exp) = E^exp. The E^HF (α_exp) total energies are always higher than the Hartree–Fock energies for the atoms. The relation of the calculated ionization potential to the experimental ionization potential depends on the α used to define u_i(α), α_exp, or α_HF.     

Use of noninteger n‐Slater type orbitals in combined Hartree–Fock–Roothaan theory for calculation of isoelectronic series of atoms Be to Ne

The ground state calculations in the combined Hartree–Fock–Roothaan approach are performed for the neutral and the first 20 cationic members of the isoelectronic series of atoms from Be to Ne using noninteger n‐Slater type orbitals. For the total energies obtained, only a small deviation has been found. At the same time, the size of the present noninteger n‐Slater type orbitals is smaller than that of the usual extended integer n‐Slater functions in literature. All of the nonlinear parameters are fully optimized. The relationship between optimized parameters and atomic number Z is also investigated. For each atom, the total energies are given in tables. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

Perturbation treatment of the Hartree–Fock equations of 1s2p3s 4Pθ state of the lithium isoelectronic sequence

The first order Hartree–Fock equations of the 1s2p3s ⁴P⁰ state of the three-electron atomic systems have been solved exactly. These solutions are used to evaluate Hartree–Fock energy up to third order with high accuracy. The third order Hartree–Fock energies for Li to Ne⁷⁺ are compared with those derived from experiment and other theoretical calculations.     

The modified Hartree–Fock self-consistent field equation

A one-electron correlation operator is introduced into the Hartree–Fock self-consistent field equation. The correlation operator is derived from the second-order perturbation theory. Energies of atomic and molecular systems calculated from this modified Hartree–Fock equation are equal to that from second-order perturbation of Hartree–Fock equation. The modified equation can also be solved self-consistently by the LCAO approximation. We also presented the modified expressions for other operators.     

The Hartree–Fock ground state of atomic two-electron systems and the Wilson–Silverstone 1s orbital

For the Hartree–Fock ground state of atomic two-electron systems, the variational function of Wilson and Silverstone, ?(r) = (a + kr)^?1 exp(-kr) / (4π)^1/2, can be optimized in two complementary ways. For small values of the atomic number Z, all intergrals have been calculated numerically and optimization can be performed accurately. However, as Z increases, loss of significant figures is increasingly detrimental to the optimization process. For sufficiently large values of Z, the integrals may be replaced by asymptotic expansions in terms of (2a)^?. As a result of optimization, the parameters and expectation values can be given as expansions in terms of (32Z)^?1/2. Both methods yield good results for Z ≈ 25, so that the whole range of Z can be treated accurately. The results have been compared with those derived from other analytical two-parameter functions. It is found that ?(r) is indeed the outstanding two-parameter function, at least for small and intermediate values of Z.     

Hartree–Fock instabilities and electronic properties

Hartree–Fock instabilities are investigated for about 80 compounds, from acetylene to mivazerol (27 atoms) and a cluster of 18 water molecules, within a double ζ basis set. For most conjugated systems, the restricted Hartree–Fock wave function of the singlet fundamental state presents an external or so‐called triplet instability. This behavior is studied in relation with the electronic correlation, the vicinity of the triplet and singlet excited states, the electronic delocalization linked with resonance, the nature of eventual heteroatoms, and the size of the systems. The case of antiaromatic systems is different, because they may present a very large internal Hartree–Fock instability. Furthermore, the violation of Hund's rule, observed for these compounds, is put in relation with the fact that the high symmetry structure in its singlet state has no feature of a diradical‐like species. It appears that the triplet Hartree–Fock instability is directly related with the spin properties of nonnull orbital angular momentum electronic systems. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 483–504, 2000     

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.     

Spin–other–orbit and spin–orbit interactions in ƒ2 and ƒ3 electron configurations

Expressions of the matrix elements of the spin–other–orbit and spin–orbit interactions for the various multiplets of all the states of ?²- and ?³-electron configurations are reported and used to evaluate the Hartree–Fock values of these interactions in the neutral atoms Ce(4?²), Pr(4?³), Ho(4?¹¹) and Er(4?¹²). The required values of the spin–spin parameters M, and the spin-orbit parameter ζ for these atoms were obtained using numerical Hartree–Fock wave functions.     

Hartree-Fock soft Coulomb hole to estimate correlation energies in atoms, molecules and polymers

We have shown that the empirical correction introduced into the Hartree-Fock method to calculate correlation energies for atoms and therefore to remove the error caused by the so-called Coulomb hole can be extended from atoms to molecules and polymers. A reformulation was required of the necessary parameter representation. The reparametrization has been performed staying as close as possible to the original expressions for atoms reported by Chakravorty and Clementi (S.J. Chakravorty and E. Clementi, Phys. Rev. A, 39 (1989) 2290). In addition to their work, where the correlation energy has been calculated with the self-consistent Hartree-Fock wavefunction and the correction integrals, we have performed investigations, including the perturbation operator in the Fock operator, so that the total energy also contains the correlation energy. The applications of this approach to atoms and molecules show that the total electron correlation energies and ionization potentials calculated as differences of total energies can be obtained very satisfactorily. On the basis of the reported calculations it turns out that one obtains better agreement with reference values of more sophisticated calculations when the correction integrals are used to build up the Fock matrix. Furthermore we have found that the magnitude of the correlation energy depends only weakly on the size of the basis sets, which makes this empirical method very attractive for its application to large molecular and polymeric systems.