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.     

Extended Hartree–Fock calculations for the helium ground state

The extended Hartree–Fock (EHF) wave function of an n-electron system is defined (Löwdin, Phys. Rev. 97 , 1509 (1955)) as the best Slater determinant built on one-electron spin orbitals having a complete flexibility and projected onto an appropriate symmetry subspace. The configuration interaction equivalent to such a wavefunction for the ¹S state of a two-electron atom is discussed. It is shown that there is in this case an infinite number of solutions to the variational problem with energies lower than that of the usual Hartree–Fock function, and with spin orbitals satisfying all the extremum conditions. Two procedures for obtaining EHF spin orbitals are presented. An application to the ground state of Helium within a basic set made up of 4(s), 3(p₀), 2(d₀) and 1 (f₀) Slater orbitals has produced 90% of the correlation energy.     

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.     

Exact and approximate Hartree–Fock calculations for extended metallic systems

Starting out with the electron gas, we make a survey of the reasons for the singularity in the derivative of the orbital energy with respect to the wave number at the Fermi level for a realistic extended metallic system. Some properties of the occupation function are reviewed and it is pointed out that the direct reason for the singularity resides in a divergent lattice sum originating in the exchange part of the orbital energy. Numerical aspects are discussed, in particular with reference to the difficulty in detecting this singularity in actual computations.     

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.     

Electronic structure calculations of 1,3-dipolar cycloadditions using density functional and Hartree–Fock methods

Local spin density (LSD ) methods were used to study the concerted 1,3-dipolar cycloadditions for fulminic acid plus acetylene, fulminic acid plus ethylene, and nitrone plus ethylene. Cartesian Gaussian double-zeta split-valence basis sets augmented with one set of polarization functions (DZVP ) were used for the LSD calculations. The LSD calculations were performed with the LSD exchange functional (Dirac) and with the Vosko, Wilk, and Nusair correlation energy functional (VWN ). Nonlocal spin-density corrections (NLSD ) were estimated with the exchange functional of Becke and the correlation energy functional of Perdew (VWN + BP ) and Becke, and the correlation energy functional of Lee, Yang, and Parr (B-LYP ). Vibrational frequencies were computed at the VWN and B-LYP levels by numerical differentiation of the analytical first derivatives of the energy. Each of these reactions was examined using Hartree–Fock and Møller–Plesset perturbation theory for comparison. Geometry optimizations were carried out at the Hartree–Fock level with the 6-311G (d,p) basis set, and correlation energies were computed up to the MP 4SDTQ /6-311G (d,p) level of theory. For the reactions of fulminic acid plus acetylene, fulminic acid plus ethylene, and nitrone plus ethylene, our best estimated density functional barrier heights are 7.8 ± 1.5, 8.9 ± 0.3, and 11.05 ± 1.9 kcal/mol, respectively. These results are in reasonable agreement with the correlated wave-function calculations and experimental estimates. © 1994 John Wiley & Sons, Inc.     

Parallel pseudospectral electronic structure: I. Hartree–Fock calculations

We present an outline of the parallel implementation of our pseudospectral electronic structure program, Jaguar, including the algorithm and timings for the Hartree–Fock and analytic gradient portions of the program. We also present the parallel algorithm and timings for our Lanczos eigenvector refinement code and demonstrate that its performance is superior to the ScaLAPACK diagonalization routines. The overall efficiency of our code increases as the size of the calculation is increased, demonstrating actual as well as theoretical scalability. For our largest test system, alanine pentapeptide [818 basis functions in the cc-pVTZ(-f) basis set], our Fock matrix assembly procedure has an efficiency of nearly 90% on a 16-processor SP2 partition. The SCF portion for this case (including eigenvector refinement) has an overall efficiency of 87% on a partition of 8 processors and 74% on a partition of 16 processors. Finally, our parallel gradient calculations have a parallel efficiency of 84% on 8 processors for porphine (430 basis functions). © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 1017–1029, 1998     

Computation of Fourier transform quantities in Hartree–Fock calculations for simple crystals

An efficient method is presented to compute the Fourier transforms of lattice sums over Slater-type orbital products that arise in crystal Hartree–Fock calculations. Introduction of one-dimensional integral representations for the Fourier transforms of the STO 'S enables a separation of the three infinite sums under the (two-dimensional) integral sign. This, in turn, makes it possible to calculate and store quantities associated with the three Fourier transform vector components separately. The orbital transform integrations are then performed numerically. The method is very advantageous when lattice sum transforms are needed for a large number of transform vectors. Formulas and results are presented for simple-cubic crystals when 1s, 2s, and 2p Slater orbitals are involved.     

Improved uncoupled Hartree–Fock (IUCHF) perturbation methods and bounds for the second-order energy in coupled Hartree–Fock perturbation theory

Different kinds of improved uncoupled Hartree–Fock methods are proposed for the calculation of second-order perturbation energies. Using these methods inequalities are derived for the error of the uncoupled procedure with geometric approximation.     

B‐spline one‐center method for molecular Hartree–Fock calculations

We introduce one‐center method in spherical coordinates to carry out Hartree–Fock calculations. Both the radial wave function and the angular wave function are expanded by B‐splines to deal with electron‐nucleus cusps, which results in the improved convergence for several typical closed‐shell diatomic molecules with moderate basis numbers. B‐splines could represent both the bound state and continuum state wave functions properly, and the present approach has been applied to investigate the ionization dynamics for H₂ in the intense laser field adopting single‐active‐electron model. © 2013 Wiley Periodicals, Inc.     

Uncoupled Hartree–Fock calculations of the polarizability and hyperpolarizabilities of nitrophenols

The polarizability and hyperpolarizabilities of nitrophenols as model compounds for studying nonlinear optics have been investigated at the Hartree-Fock level of approximation by means of the Dalgarno Uncoupled Hartree-Fock (DUHF) or Sum Over Orbitals (SOO) method. The additive character and the charge transfer effects in α, β, and γ have been analyzed in terms of the σ and π molecular orbital contributions, the contribution of the individual π molecular orbitals, and the contribution of the highest occupied and the lowest unoccupied molecular orbitals. Within the SOO approach, the reliability of the Two-Level Model has been tested and the influence of the rotation of the nitro group and of the presence of the intramolecular hydrogen bonding in ortho-nitrophenol have been studied. The results show that the present method is a reliable and efficient tool for the prediction of trends in the molecular polarizability and hyperpolarizabilities of large molecules. © 1995 by John Wiley & Sons, Inc.     

Three-dimensional Hartree–Fock crystal-orbital calculations on conducting polymers: trans-polyacetylene and polythiophene

Ab initio Hartree–Fock crystal orbital quantum chemical calculations on various structural models of three-dimensional crystals of trans-polyacetylene and polythiophene were carried out. The results provide insight into the actual structure and symmetry of the crystalline polymers under study, which are not easily amenable to experimental determination. Both conducting polymers under study were calculated to form crystals with monoclinic unit cells. A possible molecular basis of the specific electronic properties of these polymers is discussed. © 1998 John Wiley & Sons, Inc. Int J Quant Chem 68: 421–429, 1998     

Exact-exchange Hartree–Fock calculations for periodic systems. I. Illustration of the method

A scheme is presented for performing linear-combination-of-atomic-orbitals (LCAO ) self-consistent-field (SCF ) ab initio Hartree–Fock calculations of the electronic structure of periodic systems. The main aspects which characterize the present method are (i) a thorough discussion of both translational and local symmetry properties and the derivation of general formulas for the transformation of all the relevant monoelectronic and bielectronic terms under symmetry operators. (ii) The use of general yet powerful criteria for the truncation of infinite sums; in particular, the Coulomb electron–electron interactions are subdivided into terms corresponding to intersecting or nonintersecting charge distributions; the latter are grouped into shell contributions and the interaction is evaluated by multipolar expansions; the exchange interaction may be evaluated with great precision by retaining a relatively small number of two-electron integrals according to a truncation criterion which fully preserves its nonlocal character. (iii) The use of a procedure for performing integrals over k , as needed in the evaluation of the Fermi energy and in the reconstruction of the Fock matrix, which is particularly simple because it employs partially intersecting small spheres as integration subdomains where linear extrapolation is admitted. A comparison is finally made of our fundamental equations in the critical SCF stage with those obtainable by a recent proposal which uses Fourier transforms to express Coulomb and exchange integrals.     

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     

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.     

Triazene proton affinities: A comparison between density functional,Hartree–Fock,and post-Hartree–Fock methods

The consistency of three density functional computational implementations (DMol, DGauss, and deMon) are compared with high-level Hartree–Fock and Møller–Plesset (MP) calculations for triazene (HN?NNH₂) and formyl triazene (HN?NNHCOH). Proton affinities on all electronegative sites are investigated as well as the geometries of the neutral and protonated species. Density functional calculations employing the nonlocal gradient corrections show agreement with MP calculations for both proton affinities and geometries of neutral and protonated triazenes. Local spin density approximation DMol calculations using numerical basis sets must employ an extended basis to agree with other density functional codes using analytic Gaussian basis sets. The lowest energy conformation of triazene was found to be nonplanar; however, the degree of nonplanarity, as well as some bond lengths, is dependent on the basis set, electron correlation treatment, and methods used for the calculation. © 1994 by John Wiley & Sons, Inc.

1 This article is a U.S. Government work and, as such, is in the public domain in the United States of America.