Half-projected and projected Hartree-Fock calculations for singlet ground states. II. Lithium hydride

The half-projected Hartree–Fock function (HPHF ) for singlet states is defined as a linear combination of two Slater determinants which contains only spin eigenstates with even spin quantum numbers. The possible uses of such an approach for determining molecular properties are investigated computing the potential energy curve, binding energy, force constant, and dipole moment variation corresponding to the lithium hydride ground state. Full projected and restricted Hartree–Fock calculations (PHF and RHF ) are performed simultaneously for comparison purposes. It is found that the HPHF model yields very satisfactory results, very close to those of the PHF scheme. Both models predict properly the molecular behavior as a function of nuclear separation, whereas the RHF one fails. A discussion is given in terms of configuration equivalents. It is concluded that the HPHF scheme seems to be useful for determining molecular properties specially in the case of large systems in which the more sophisticated methods are unmanageable.     

Spin-unrestricted calculations in quantum chemistry

The unrestricted complete active space self-consistent field (UCASSCF ) function is defined, and a proof that a UCASSCF eigenfunction of the spin operator S ² is a CASSCF function is given. The spin-contamination for an unrestricted Hartree–Fock (UHF ) function is evaluated by using Araki angle operators, and the UHF function is then projected on the restricted open-shell Hartree–Fock (ROHF ) space. The present analysis has deep consequences since it implies that the only non-spin-contaminated UHF functions are the ROHF functions. This is illustrated in a calculation of the spin density of He. © 1993 John Wiley & Sons, Inc.     

Spin-Extended Hartree–Fock calculation for hydrocarbon π-radicals

An iterative method of solving the Mestechkin spin-extended Hartree–Fock equations is proposed. On this basis projected and spin-extended values of energies and spin densities for some hydrocarbon π-radicals are calculated. Relations among UHF , ASA , AA , projected, and spin-extended energies are discussed.     

Communication: ROHF theory made simple

Restricted open-shell Hartree-Fock (ROHF) theory is formulated as a projected self-consistent unrestricted HF (UHF) model by mathematically constraining spin density eigenvalues. This constrained UHF (CUHF) wave function is identical to that obtained from Roothaan's effective Fock operator. The α and β CUHF Fock operators are parameter-free and have eigenvalues (orbital energies) that are physically meaningful as in UHF, except for eliminating spin contamination. This new way of solving ROHF leads to orbitals that turn out to be identical to semicanonical orbitals. The present approach removes ambiguities in ROHF orbital energies.     

Electronic structure of long cumulene chains

The analysis of the equations of the unrestricted Hartree–Fock (UHF) method for polyenes C_NH_N+2 with even and odd N » 1 is carried out. The equations of the UHF method are shown to be the same in both cases. The comparison of the UHF method with the extended Hartree–Fock (EHF) method applied to large systems is performed. The ground state and π-electron spectra of long cumulene chains C_NH₄ are treated by the EHF Method. The end effects are taken into consideration. It is shown that the EHF method gives a finite value of the first optical transition frequency and, at the same time, zero value of torsion barrier of end CH₂–groups in long cumulene chains (N → ) in contrast to previous calculations of cumulenes by the Huckel method and the restricted Hartree–Fock method.     

On the nature of electron correlation in C60

The ground state restricted Hartree Fock (RHF) wave function of C(60) is found to be unstable with respect to spin symmetry breaking, and further minimization leads to a significantly spin contaminated unrestricted Hartree Fock (UHF) solution ( = 7.5, 9.6 for singlet and triplet, respectively). The nature of the symmetry breaking in C(60) relative to the radicaloid fullerene, C(36), is assessed by energy lowering of the UHF solution, , and the unpaired electron number. We conclude that the high value of each of these measures in C(60) is not attributable to strong correlation behavior as is the case for C(36). Instead, their origin is from the collective effect of relatively weak, global correlations present in the π space of both fullerenes. Second order perturbation (MP2) calculations of the singlet triplet gap are significantly more accurate with RHF orbitals than UHF orbitals, while orbital optimized opposite spin second order correlation (O2) performs even better.     

Spin-projected EHF method. IV. Comparison of potential curves given by different one-electron methods

Some aspects of the computer realization of the spin-projected extended Hartree–Fock (EHF ) method at the ab initio level are briefly discussed for the algorithm of solution developed in the previous papers of this series. Calculations have been performed for the BH molecule by using a small basis of contracted Gaussian lobes with the purpose of comparing the potential curves given by the different one-electron methods RHF , UHF , UHF with subsequent spin projection, and EHF . It is concluded that the UHF and, in particular, the EHF methods give a qualitatively correct shape of the potential curve; the RHF method shows the known incorrect dissociation behavior while the potential curve obtained by subsequent spin projection of the UHF wave function exhibits spurious extrema at intermediate internuclear separations.     

Corrections to coupled Hartree–Fock theory: the rearrangement effect

In this communication we discuss a method of incorporating corrections to the coupled Hartree–Fock (CHF ) formalism by introducing the so-called “rearrangement effect.” In this we take account of the relaxation of the core orbitals when excitations from a starting Hartree–Fock wave function occur. The magnitude of this correction numerically is found to be quite significant for the polarizabilities of two-electron atomic systems, results for which are reported.     

Ab initio study of organic mixed valency

A series of six radical cations of the type (D L D)⁺ was investigated at the ab initio unrestricted Hartree–Fock level. One localized and one delocalized conformation were systematically searched by full geometry optimization. At both nuclear arrangements, mostly found as being minima in the symmetry‐restrained Hartree–Fock framework, excitation energies were calculated through the expansion of the wave function on single electronic excitations of the Hartree–Fock fundamental determinant and at the unrestricted Hartree–Fock or at the multiconfigurational self consistent field levels. Few calculations were also performed by taking into account some part of the electronic correlation. Except for N,N,N′,N′‐tetramethyl p‐phenylenediamine, all the studied compounds are localized stable cations, at the symmetry‐restrained Hartree–Fock level. However, the reoptimization of their wave function changes this observation since only three of them seem to conserve a localized stable conformation. Most of the studied systems are characterized by one or two excited electronic states very close to the fundamental one and should thus present an unresolved broadened first absorption band in the near‐infrared region. These features are in agreement with the available experimental data. Strong Hartree–Fock instabilities are found for the delocalized structure and put in relation with the existence of the large nonadiabatic coupling in this conformational region. The solvent influence is discussed in the Onsager dipolar reaction field framework. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 76: 552–573, 2000     

Bond length alternation in cyclic polyenes. II. Unrestricted hartree–fock method

The problem of bond length alternation in cyclic polyene models as described by the Pariser–Parr–Pople π-electron Hamiltonian, together with an empirical quasi harmonic σ-core potential is investigated using the unrestricted Hartree–Fock wave function employing different spatial orbitals for different spins. It is shown that in contrast to the restricted Hartree–Fock method, which favors bond alternation in large cyclic polyenes, the unrestricted Hartree–Fock method stabilizes the symmetric structures with equidistant internuclear separation. An assessment of the amount of correlation error recovered by the unrestricted Hartree–Fock procedure is examined and the qualitatively different behavior of the cyclic polyene models when described by restricted and unrestricted Hartree–Fock wave functions is discussed from this viewpoint.     

Formulation of unrestricted and restricted Hartree–Fock–Bogoliubov equations

The Hartree–Fock–Bogoliubov (HFB) method, dealing with Bogoliubov orbitals, which consist of particle and hole part, can provide states with pair correlations associated with Cooper pairs. The dimension of HFB Fock matrices can be reduced by restrictions of spin states of Bogoliubov orbitals similarly to ordinary Hartree–Fock (HF) equations such as restricted HF (RHF), unrestricted HF (UHF), and generalized HF (GHF). However, there are few studies of moderate restricted HFB equations such as UHF‐based HFB equations. In this article, formulation and calculations of restricted HFB equations are described. The solutions of general and restricted HFB equations are compared. Pair correlations taking account of restricted and general HFB equations are discussed. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem, 2004     

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     

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.     

Unrestricted-Hartree–Fock wave functions approximating low-lying covalent states of ring π systems

We show that spin projected unrestricted-Hartree–Fock (PUHF ) wave functions are able to well approximate some low-lying covalent states of ring π systems. The UHF wave functions belong to either the axial or torsional spin density wave class. Their spin structures are found to be approximations for the spin correlation structures of the corresponding exact wave functions. The PUHF wave functions become close to exact eigenstates with homopolar valence bond characters in the strong correlation limit.     

Time-dependent Hartree–Fock calculations in the Pariser–Parr–Pople model. Applications to aniline,azulene and pyridine

Propagator or Green's function methods are used to analyze the time-dependent Hartree–Fock model. The non-hermitian matrix problem for the time-dependent Hartree–Fock solution is reduced to a problem related only to hermitian matrices. Particular attention is given to the calculation of oscillator strength in different approximations. The connection between the stability of the Hartree–Fock solution and the solution of the time-dependent Hartree–Fock problem is demonstrated. The results of numerical calculations are given for aniline, azulene and pyridine.     

Broken symmetry solutions of the extended Hubbard model

This paper describes a group theoretical classification of the unrestricted Hartree–Fock (UHF ) solution of the extended Hubbard model including nearest-neighbor repulsion and exchange interactions. We derived systematically all types of the broken symmetry solution with a single-order parameter for the ordering vector Q = (π,π). It is shown that there are four nonmagnetic states: a charge density wave (CDW), a charge current wave (CCW), and two bond-order waves (BOW), and five magnetic states: an ordinary antiferromagnetic state (AF), a spin current wave (SCW), two axial spin bond-order waves (ASBOW), and an helical spin bond-order wave (HSBOW). We showed that these phases can be derived naturally from the irreducible decomposition of the Fourier transformed interaction and the self-consistent field (SCF ) condition. Some numerical results are given for the half-filled case.     

Quantum chemistry of many-electron systems: high-priority aspects

The review generalizes the studies devoted to the development of a new quantum chemistry method representing an alternative to the Hartree–Fock approximation. Based on the hypothesis of prohibition of equipotential surfaces, which clarifies the physical sense of the Pauli exclusion principle, and taking account of the condition for antisymmetrical wave function of the triplet state (3S) of He atom, the Hartree–Fock approximation is inappropriate for a priori determination of the nodal surfaces of many-electron wave functions (MWFs) for the test systems traditionally used in quantum chemistry, namely, excited triplet state of H₂ molecule and the ground electronic states of Li atom and LiH molecule. The nodal surfaces of the wave functions corresponding to the minimum basis set of Slater orbitals in the Hartree–Fock approximation are constructed and analyzed. An alternative to the Hartree–Fock approximation is provided by the MWF quantum chemical method being developed by the authors. In the MWF method, the nodal surfaces for H₂(³Σ _u ^v ) and Li(²S) are specified a priori. Some aspects of geometric interpretation of the Pauli exclusion principle are discussed. Unlike the MWF method, the Hartree–Fock approximation is unsuitable for taking account of the dependence of the MWF nodal surfaces on the nuclear charges and on correlation effects related to the motion of electrons with antiparallel spins because such nodal surfaces are predefined by the mathematical properties of Slater determinants rather than by physically clear and more practically valuable algebraic products of electrostatic potential differences.