Variation–iteration method for one-dimensional two-electron systems

A quantum chemical study of the two low-lying quartet states of seven model compound I iron–porphyrin complexes with varying axial ligands has been carried out using the INDO method. The varying axial ligands included in this study are five that are models for those in the intact enzymes: imidazole and imidazolate (model peroxidase HRP and CCP), CH₃CONH₂ (Gln175 mutant of CCP), PhO^?1 (catalase), CH₃S^?1 (P450), and two that have been used in biomimetics of these enzymes: Cl^?1 (hemin) and PhS^?1 (model P450s). The purpose of these studies was to determine the role of the axial ligands in determining (i) the relative energies of the two nearly degenerate quartet electronic states of compound I, involved either as an a_1u or a_2u porphyrin π cation radical and (ii) the electron and spin distributions in the a_1u and a_2u radical cations of compound I. For most of the model complexes, including both HRP-I and CAT-I, a moderate effect of the axial ligand on the relative energy of these two states was observed and the a_1u radical cation was found to be the ground state. The energy order of these two radical cations, however, was reversed in the P450-I model complexes, indicating an association of the unique property of the Fe?O bond breaking with an a_2u radical cation. The symmetry-allowed overlap between the Fe?O and 3a_2u orbitals may lower the activation energy for the Fe?O bond cleavage in P450-I. However, the calculated electronic and spin properties, including the unpaired spin and net charge on the oxygen and the Fe?O bond overlap density, important determinants of the reactivity of this complex in the ligand–Fe?O region, are very similar for all complexes and in both cation states. © 1992 John Wiley & Sons, Inc.     

Half-projected Hartree–Fock natural orbitals for defining CAS–SCF active spaces

We have extended the range of systems to which the half-projected Hartree–Fock (HPHF ) method has been applied, including the triplet state of the wave function. In our implementation, DIIS overcomes the convergence difficulties reported in earlier studies. HPHF allows generation of a symmetry-broken wave function in regions of the potential energy surface where the RHF wave function is triplet-stable. The fractionally occupied natural orbitals (FONOS ) of the HPHF wave function are good starting vectors for CAS –SCF calculations. A CAS –SCF in the space defined by the HPHF FONOS should be used instead of the unrestricted natural orbital CAS –SCF method in regions of triplet stability and for small active space problems. We draw extensive comparisons between the results of both the UNO –CAS and HPNO –CAS methods and those of full CAS –SCF calculations. © 1993 John Wiley & Sons, Inc.     

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.     

Perturbative calculation of the Hartree–Fock interaction energy using orthogonalized orbitals

A many‐body perturbation theory based on the partitioning of the dimer Hamiltonian, formulated in an orthogonalized basis set, is used for the calculation of interaction energies at the Hartree–Fock (HF) level. Numerical results for the (HF)₂ and (H₂O)₂ systems in selected geometries are presented. The interaction‐energy components are compared with the results obtained from the standard supermolecular approach and the intermolecular perturbation theory based on the biorthogonal basis set. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 75: 81–88, 1999     

Numerical solution of Hartree–Fock equations for a polyatomic molecule: Linear H3 in momentum space

We give an example of multidimensional numerical solutions in momentum space, for the Hartree–Fock equations of the simplest molecule with more than two electrons and two nuclei.     

Ground-state multiplicity by spin-extended Hartree–Fock method

New possibilities of the spin-extended Hartree–Fock method in determining the sequence of energy levels are analyzed and demonstrated by the example of some oxygen-contained compounds.     

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.     

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.     

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.     

Decoupled Hartree–Fock methods. I. Calculations for atoms with orthogonalized orbitals

By a proper approximation of the interaction term in a many-electron Hamiltonian the Hartree-Fock equations are decoupled. Making use of this simplification one obtains a good initial guess for the wave function with minimal computational work. Refining the procedure, the exact HF limit can be achieved. Total energies, ionization potentials and excitation energies for light atoms are calculated.     

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.     

Local Hartree–Fock orbitals using a three‐level optimization strategy for the energy

Using the three‐level energy optimization procedure combined with a refined version of the least‐change strategy for the orbitals—where an explicit localization is performed at the valence basis level—it is shown how to more efficiently determine a set of local Hartree–Fock orbitals. Further, a core–valence separation of the least‐change occupied orbital space is introduced. Numerical results comparing valence basis localized orbitals and canonical molecular orbitals as starting guesses for the full basis localization are presented. The results show that the localization of the occupied orbitals may be performed at a small computational cost if valence basis localized orbitals are used as a starting guess. For the unoccupied space, about half the number of iterations are required if valence localized orbitals are used as a starting guess compared to a canonical set of unoccupied Hartree–Fock orbitals. Different local minima may be obtained when different starting guesses are used. However, the different minima all correspond to orbitals with approximately the same locality. © 2013 Wiley Periodicals, 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.     

Form factors directly determined from momentum space Hartree—Fock solutions: H2 and Li2

The characteristics of the theoretical form factors computed from momentum space wave functions of quantum chemistry (i.e., vectorial quantities including appropriate phase factors) are examined in connection with the so-called atomicity assumption of crystallography. Using molecular orbitals previously computed by direct numerical integration of momentum space Hartree-Fock equations for the H₂ and Li₂ molecules, the patterns of parallel and perpendicular form factors with respect to bond axes are shown to be more significant than are the corresponding spherical averages. The case of the triatomic system H₃ is briefly considered. © 1997 John Wiley & Sons, Inc. Int J Quant Chem 63: 451–457, 1997     

Spin-Adapted restricted Hartree–Fock reference coupled-cluster theory for open-shell systems: Noniterative triples for noncanonical orbitals

The coupled clusters singles and doubles (CCSD ) method for calculations of open-shell systems with the single restricted Hartree–Fock (ROHF ) reference determinant is extended by the noniterative triples to give CCSD(T) . Our approach profits from the fact that (a) single- and double-excitation amplitudes are spin-adapted, which directly leads to a computationally less demanding algorithm than are nonadapted procedures and produces the spin-adapted CCSD wave function and (b) triple excitations calculated from converged spin-adapted (SA ) CCSD amplitudes are also obtained more effectively. Altogether, computer demands of our SA CCSD(T) approach, applicable to high-spin open-shell cases which are well represented by a single-determinant reference is comparable to that for closed-shell systems. Our approach is not based on semicanonical orbitals, applied by Bartlett's group. However, we compare some other possible choices of ROHF orbitals to this “standard.” Numerical results for a series of atoms and molecules demonstrate little sensitivity to this selection. © 1995 John Wiley & Sons, Inc.