Ground-state and some excited states of Li2 by the half-projected Hartree–Fock method

The half-projected Hartree–Fock wave function (HPHF ) is one of simplest models for introducing some electronic correlation effects. In this model, the wave function is built up with only two Slater determinants. This simple form suggests its application for the direct determination of singlet excited states. On the other hand, because the HPHF model does not mix singlet and triplet states with M_s = 0, it can be used for determining independently singlet and triplet states without any mutual contamination. In the present work, we applied this model to determine nine electronic states of the lithium molecule; one of them exhibits even the same symmetry of the fundamental one. For this purpose, the 6-311G (d) basis was used. Potential energy curves were determined and some spectroscopic constants derived. The numerical results were compared with the available experimental data, as well as with other theoretical values. © 1995 John Wiley & Sons, Inc.     

Coulomb-Hole–Hartree–Fock functional

We have extended to molecules a density functional previously parametrized for atomic computations. The Coulomb-hole–Hartree–Fock functional, introduced by Clementi in 1963, estimates the dynamical correlation energy by the computations of a Hartree–Fock-type single-determinant wave function, 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 revisited by S. Chakravorty and E. Clementi [Phys. Rev. A 39 , 2290 (1989)], where a Yukawa-type soft Coulomb hole replaced the previous hard hole; atomic correlation energies, computed for atoms with Z = 2 to Z = 54 as well as for a number of excited states, validated the method. In this article, we parametrized a function, which controls the width of the soft Coulomb hole, by fitting the first and second atomic ionization potentials of the atoms with 1 ? Z ? 18. The parametrization has been preliminarily validated by computing the dissociation energy for a number of molecules. A few-determinant version of the Coulomb-hole–Hartree–Fock method, necessary to account for the nondynamic correlation corrections, is briefly discussed. © 1994 John Wiley & Sons, Inc.     

Comparative study of the Hartree–Fock and the Hartree–Fock–Slater method

The Hartree–Fock method (standard Roothaan closed-shell HF –LCAO theory) and the Hartree–Fock–Slater method (restricted HFS –LCAO –DV method developed by Baerends and Ros) have been compared with emphasis on the respective one-electron equations and on the matrix elements of the respective Fock operators. Using the same STO basis in the two cases, the matrix elements of the Fock operators and of their separate one-electron, Coulomb, and exchange contributions have been calculated for the same orbitals and density of the ground state of the diatomic molecule ZnO. The effects of methodical (exchange potential) and numerical (DV method, density fit) differences between the HF and HFS methods on the various matrix elements have been analyzed. As expected the methodical effect prevails and is responsible for the higher (less negative) values of the matrix elements of the HFS Fock operator compared to those of the HF Fock operator. Numerical effects are observable also and are caused by the difference in integration procedures (DV method), not by the density fit.     

Simple derivation of conditions for instability in the Hartree–Fock and projected Hartree–Fock schemes

The conditions for instability of solutions of Hartree–Fock and projected Hartree–Fock equations are derived in a form involving finite real symmetric matrices. These conditions are also expressed in terms of the Fock–Dirac density matrix, both at the spin–orbital and at the orbital level. The particular variations which give rise to the so-called singlet and triplet instabilities are described.     

Gaussian overlap approximation in the projected Hartree–Fock method

A method for the approximate calculation of matrix elements with respect to projected Hartree–Fock wave functions is proposed. The method is tested on some calculations in the many-parameter AMO method. It is found that the approximation reduces the amount of work, involved in the evaluation of the energy, by a factor of five and that it reproduces the exact values to within a few per cent.     

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     

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.     

Multiple time scale Hartree–Fock molecular dynamics

This is the first application of a rigorous, established multiple time-step method to ab initio molecular dynamics. The resulting algorithm is conceptually simple and easy to implement, but very effective. It translates the large mass differences present in ab initio molecular dynamics into substantial savings in computer time while retaining high accuracy. This transforms ab initio molecular dynamics from a desirable but prohibitively expensive possibility into a viable method, at least for short-time phenomena in small systems or for otherwise inaccessibly complicated potential energy surfaces. © 1993 John Wiley & Sons, Inc.     

Spline methods for multiconfiguration Hartree–Fock calculations

The earlier numerical multiconfiguration Hartree–Fock atomic structure package was not designed with high-performance computers in mind. In this paper, some new algorithms based on spline–Galerkin methods are described that are appropriate for concurrent/vector architectures. The goal is to improve the level of numerical accuracy by several orders of magnitude using fewer basis functions than points in a numerical grid. Of critical importance is the robustness of the code: The most serious problems in the numerical implementation were associated with orthogonality constraints. In a spline basis approach, the orthogonality requirements can be integrated into quadratically convergent update procedures. These procedures are evaluated for a number of cases.     

Equations de Hartree–Fock dependant du temps

The time-dependent Hartree–Fock (TDHF ) equations are derived up to the second order when the system is perturbed by a monochromatic plane wave. The solutions of the equations are subjected to the orthonormalization conditions satisfied by the orbitals. In the equations, these conditions are expressed by the appearing of coefficients λ playing the part of Lagrangian multipliers. Relations between the coefficients λ are established. These relations are equivalent to the above-mentioned orthonormalization conditions. This equivalence enables us to substitute for the solution of an integrodifferential equation system subject to constraint conditions, that of a free system. The TDHF equations obtained determine the first- and second-order orbital perturbations, which no doubt verify the orthonormalization conditions. These orbitals can be used in the calculation, up to second order, of different nonlinear optical effects.     

Variation–iteration method in momentum space: Determination of Hartree–Fock atomic orbitals

A method using the Svartholm iterative procedure to solve atomic Hartree–Fock equations in momentum space is defined and applied to the ground states of Be and B⁺. The calculated atomic orbital properties follow a monotonic and stable convergence, but with rates of convergence depending on each property. The evolution of the orbitals during the iterations is explained by the combined actions of the variational principle, the Svartholm iterative procedure, and the momentum space representation. © 1993 John Wiley & Sons, Inc.     

Molecular symmetry. IV. The coupled perturbed Hartree–Fock method

Symmetry methods employed in the ab initio polyatomic program HONDO are extended to the coupled perturbed Hartree–Fock (CPHF) formalism, a key step in the analytical computation of energy first derivatives for configuration interaction (CI) wavefunctions, and energy second derivatives for Hartree–Fock (HF) wavefunctions. One possible computational strategy is to construct Fock-like matrices for each nuclear coordinate in which the one- and two-electron integrals of the usual Fock matrix are replaced by the integral first derivatives. “Skeleton” matrices are constructed from the unique blocks of electron-repulsion integral derivatives. The correct matrices are generated by applying a symmetrization operator. The analysis is valid for many wavefunctions, including closed- or open-shell spin-restricted and spin-unrestricted HF wavefunctions. To illustrate the method, we compare the computer time required for setting up the coupled perturbed HF equations for eclipsed ethane using D_3h symmetry point group and various subgroups of D_3h. Computational times are roughly inversely proportional to the order of the point group.     

Calculation of molecular electronic spectra by projected-unrestricted Hartree–Fock theory

This paper is concerned with a new application of projected-unrestricted Hartree–Fock theory, namely, the calculation of electronic spectra for symmetric molecules. The excited electronic state is represented by a single determinant whose unrestricted nature allows for orbital rearrangement relative to the self-consistent ground state. The self-consistent calculation must be followed by spin projection to obtain appropriate spin eigenstates. It was necessary to develop modified procedures for portions of the spin projection calculation because our method of constructing the wave functions produces degeneracies among the natural orbitals. Illustrative calculations using the all-valence-electron INDO approximations produced results which compared favorably with configuration-interaction treatments. The method described here should be most useful, however, in conjunction with ab initio calculations using flexible basis sets.     

Study of potential curves by UHF-type methods. VI. Spin-extended Hartree–Fock method

An improved computational procedure of the spin-extended Hartree–Fock method is presented. The main features of the EHF potential curves for the dissociation of the LiH, BeH, and Li₂ molecules are discussed. The results are compared with curves calculated by the perturbation or CI methods using the same basis set.     

Convergence properties of Hartree–Fock SCF molecular calculations

Hartree-Fock equations are viewed as nonlinear algebraic equations that can be solved iteratively. Provided we assume the existence of a solution, valuable properties of convergence may be assessed. The close connection between convergence of the SCF procedure and stability properties of the solution is shown from a nonapproximate standpoint. The convergence features of level-shifting convergence-forcing techniques are analyzed. The connection between this nonlinear algebraic approach and the related gap equation is displayed and the example of the restricted Hartree-Fock hydrogen molecule is discussed.     

Some remarks on multiconfiguration time-dependent Hartree–Fock theory

Recent progress in multiconfiguration time-dependent Hartree–Fock (MC –TDHF ) theory is reviewed. The generality of a formalism first put forward by Moccia [Int. J. Quantum Chemistry 8 , 293 (1974)] is emphasized, and illustrative applications are discussed. A calculation of the frequency-dependent electric-dipole polarizability of the helium atom shows that even with a modest amount of CI, it is possible to achieve an accuracy comparable to that obtained by use of Hylleraas-type wave functions.     

Hartree–Fock MO–LCAO equations with charge-dependent atomic integrals

The Hartree–Fock equations are derived in the MO -LCAO approximation for the case when the integrals (except overlap integrals) over the atomic orbitals are charge-dependent. It is shown that inclusion of the overlap matrix in the iterative procedure gives equations which are too complicated for the simple model under consideration. The approach is applied to the VESCF method in the PPP scheme.