Convergence of LCAO density matrices in Hartree—Fock calculations on extended model chains

Some mathematical implications of the extended nature of model chains are reviewed to call attention to their importance in electronic structure calculations at the restricted Hartree-Fock Roothaan level.     

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.     

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.     

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.

     

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.     

Post Hartree–Fock and density functional theory studies on Di-Protonated Allopurinol

Post Hartree–Fock and density functional theory (DFT) methods have been employed to study the molecular properties of Di-Protonated Allopurinol²⁺ tautomers in gaseous and aqueous phase environments. The tautomers in gaseous phase have been optimized at MP2/6-311G(2d,2p) and B3LYP/6-311G(2d,2p) levels of theory. The self-consistent reaction field theory (SCRF) has been employed to optimize the tautomers in aqueous phase (ε = 78.5) at B3LYP/6-311G(2d,2p) level of theory and the solvent effect has been studied. The structure, energetics and relative stabilities of the tautomers have been analyzed both in gaseous and aqueous phases. The principle of maximum hardness (MHP) has been tested at B3LYP/6-311G(2d,2p) level of theory. The condensed Fukui functions have been calculated using the atomic charges obtained through Natural population analysis to identify the relative change in the most reactive site of the optimized structures. NMR studies have been carried out, on the basis of Cheeseman coworker’s method, to analyze the molecular environment as well as the delocalization activities of electron clouds.     

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     

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.     

On the connections between Brueckner–coupled-cluster,density-dependent Hartree–Fock,and density functional theory

The existence of an effective one-particle Hamiltonian in the Brueckner coupled cluster model naturally leads to the definition of an effective interaction G, which is a function of the T₂ amplitudes. Two types of approximations to G are proposed: One is purely phenomenological, while the other is based on approximations to the Brueckner T₂ equation. In both cases, the resulting effective interaction may be viewed as electron-density-dependent. Generalizing Hartree–Fock theory to accommodate density-dependent interactions (DDHF ), a method is obtained that is capable of accounting for correlation effects in an independent particle framework. The heuristic Skyrme force, successfully used in nuclear physics to model nucleon–nucleon interactions, is presented here as an example of an effective electron–electron correlation interaction. Due to the δ-function character of the Skyrme force, it is possible to express the energy in this model by an integral over an energy density, thus formally providing a connection between DDHF and density functional theory for this particular case. An approximation to the Brueckner T₂ equation is also proposed in the coordinate representation. In this model, the density-matrix dependence of T₂ is reduced to a nonlocal electron density dependence by means of an expansion which introduces terms that depend on the gradient of the density. The first term in this expansion amounts to a “local density approximation” to Brueckner coupled cluster theory. © 1995 John Wiley & Sons, Inc.     

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.     

Description of core excitations by time-dependent density functional theory with local density approximation, generalized gradient approximation, meta-generalized gradient approximation, and hybrid functionals

Time-dependent density functional theory (TDDFT) is employed to investigate exchange-correlation-functional dependence of the vertical core-excitation energies of several molecules including H, C, N, O, and F atoms. For the local density approximation (LDA), generalized gradient approximation (GGA), and meta-GGA, the calculated X1s-->pi* excitation energies (X = C, N, O, and F) are severely underestimated by more than 13 eV. On the other hand, time-dependent Hartree-Fock (TDHF) overestimates the excitation energies by more than 6 eV. The hybrid functionals perform better than pure TDDFT because HF exchange remedies the underestimation of pure TDDFT. Among these hybrid functionals, the Becke-Half-and-Half-Lee-Yang-Parr (BHHLYP) functional including 50% HF exchange provides the smallest error for core excitations. We have also discovered the systematic trend that the deviations of TDHF and TDDFT with the LDA, GGA, and meta-GGA functionals show a strong atom-dependence. Namely, their deviations become larger for heavier atoms, while the hybrid functionals are significantly less atom-dependent.     

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     

He2+2: A comparison between Hartree–Fock and density functional techniques

A comparison of Roothaan-Hartree–Fock methods (both restricted and unrestricted) with density functional ones (LCAO -Xα and cellular MS -Xαβ) is made using as test case the He²₂ + molecular ion. It is shown the analogy that exists between RHF and symmetry-adapted LCAO -Xα potential energy curves, as well as between UHF and symmetry-unconstrained LCAO -Xα ones. The influence of symmetry adaptation on the overall behavior of the potential energy curve is also discussed. Finally, the difference in the behavior of the LCAO -Xα and cellular MS -Xαβ calculations is explained as an artifact of the space partitioning in the latter technique. It is concluded that LCAO -Xα method is superior to cellular MS -Xαβ because it requires less effort to reach the same results and that the general behavior is similar to UHF , although the former affords a better equilibrium bond distance and a worse energy barrier than the latter.     

Molecular structure, IR spectra of 2-mercaptobenzothiazole and 2-mercaptobenzoxazole by density functional theory and ab initio Hartree–Fock calculations

Quantum chemistry calculations have been performed using Gaussian03 program to compute optimized geometry, harmonic vibrational frequency along with intensities in IR and Raman spectra and atomic charges at RHF/6-31+G*, B3LYP/6-31+G* and B3LYP/6-31++G* levels for 2-mercaptobenzothiazole (MBT, C₇H₅NS₂) and 2-mercaptobenzoxazole (MBO, C₇H₅NOS) in the ground state. The scaled harmonic vibrational frequencies have been compared with experimental FT-IR and FT-Raman spectra. The results show that the scaled theoretical vibrational frequencies is very good agreement with the experimental values. A detailed interpretation of the infrared and Raman spectra of 2-mercaptobenzothiazole and 2-mercaptobenzoxazole was reported. Comparison of calculated spectra with the experimental spectra provides important information about the ability of the computational method to describe the vibrational modes.     

New method for approximate Hartree–Fock calculations using density approximations and coulomb field corrections. I

A method is proposed that reduces the computational effort of HF calculations considerably by reducing the number of two-electron integrals that have to be calculated. The following concepts are used: (i) approximation of the electron density by only few functions for the Coulomb part of the HF matrix; (ii) modification of this approximate density, to improve its Coulomb field; (iii) in the exchange part, a basis function χ is replaced by a function \documentclass{article}\pagestyle{empty}\begin{document}$ \tilde \chi $\end{document} consisting of fewer Gaussian lobes; (iv) the error caused by this replacement is reduced by a modification of the densities \documentclass{article}\pagestyle{empty}\begin{document}$ \tilde \chi _i \tilde \chi _j $\end{document} in the exchange integrals. The computation time of the integral part is reduced by a factor 6 for molecules containing five first-row atoms as, e.g., CF₄, if one uses a 7S/3P basis set contracted to (5, 1, 1/3). The integral time increases roughly with n³, if n is the number of Gaussian lobes.     

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.     

Long-range and short-range Coulomb correlation effects as simulated by Hartree–Fock,local density approximation,and generalized gradient approximation exchange functionals

 Exchange functionals used in density functional theory (DFT) are generally considered to simulate long-range electron correlation effects. It is shown that these effects can be traced back to the self-interaction error (SIE) of approximate exchange functionals. An analysis of the SIE with the help of the exchange hole reveals that both short-range (dynamic) and long-range (nondynamic) electron correlation effects are simulated by DFT exchange where the local density approximation (LDA) accounts for stronger effects than the generalized gradient expansion (GGA). This is a result of the fact that the GGA exchange hole describes the exact exchange hole close to the reference electron more accurately than the LDA hole does. The LDA hole is more diffuse, thus leading to an underestimation of exchange and stronger SIE effects, where the magnitude of the SIE energy is primarily due to the contribution of the core orbitals. The GGA exchange hole is more compact, which leads to an exaggeration of exchange in the bond and the nonbonding region and negative SIE contributions. Partitioning of the SIE into intra-/interelectronic and individual orbital contributions makes it possible to test the performance of a given exchange functional in different regions of the molecule. It is shown that Hartree–Fock exchange always covers some long-range effects via interelectronic exchange while self-interaction-corrected DFT is lacking these effects. Received: 25 May 2002 / Accepted: 7 October 2002 / Published online: 21 January 2003 Correspondence to: E. Kraka e-mail: kraka@theoc.gn.se Acknowledgements. This work was supported financially by the Swedish Natural Science Research Council (NFR). Calculations were done on the supercomputers of Nationellt Superdatorcentrum (NSC), Link?ping, Sweden. The authors thank the NSC for a generous allotment of computer time.     

Molecular structure and vibrational spectra of indole and 5-aminoindole by density functional theory and ab initio Hartree–Fock calculations

The molecular geometry and vibrational frequencies of indole and 5-aminoindole in the ground state have been calculated by using the Hartree–Fock and density functional method with 6-311++G(d,p) basis set. The optimized geometrical parameters obtained by DFT calculations are in good agreement with the experimental values. Comparison of the observed fundamental vibrational frequencies of indole and 5-aminoindole with the calculated results by density functional and Hartree–Fock methods indicates that B3LYP is superior to the scaled Hartree–Fock approach for molecular vibrational problems. The theoretical spectrograms for FT-IR spectrum of 5-aminoindole have been constructed.