Analytical gradient of the linear combination of Gaussian‐type orbitals—local spin density energy

An expression has been derived for the analytical evaluation of the energy gradient within the linear combination of Gaussian‐type orbitals—local spin density method. This expression is valid for any exchange‐correlation energy functional which can be represented in a density gradient expansion. In practice, because the exchange‐correlation terms are fitted with auxiliary functions, one has to introduce an approximation. Results are reported of tests on diatomics that show that it is possible to attain a typical accuracy of ±0.01 a.u. on equilibrium distances, relative to the energy minimum. The formulas for molecular integral derivatives that we implemented are based on the highly efficient recurrence formulas of Obara and Saika. We report here an additional formula for angular momentum transfer which is very useful for efficient programming of the gradient. In all cases studied, the time required to compute the gradient is a fraction of the time spent to solve the self‐consistent‐field Kohn–Sham equations.

[1]  D. Salahub,et al.  Surface extended energy loss fine structure and local spin density investigation of carbidic carbon on the Ni(100) surface , 1988 .

[2]  M. Dupuis.  Energy derivatives for configuration interaction wave functions , 1981 .

[3]  R. Bartlett,et al.  Analytical gradients for the coupled-cluster method† , 1984 .

[4]  Bernard R. Brooks,et al.  Analytic gradients from correlated wave functions via the two‐particle density matrix and the unitary group approach , 1980 .

[5]  D. Salahub,et al.  Model potential for As and LCGTO MP LSD calculation of geometry, electronic structure and photoelectron spectra for As2 and As4 , 1987 .

[6]  Brett I. Dunlap.  Geometry optimization using local density functional methods , 1986 .

[7]  Michael C. Zerner,et al.  The Challenge of d and f electrons : theory and computation , 1989 .

[8]  N. H. March,et al.  Theory of the inhomogeneous electron gas , 1983 .

[9]  John R. Sabin,et al.  On some approximations in applications of Xα theory , 1979 .

[10]  R. Bartlett,et al.  Analytical gradient evaluation in coupled-cluster theory , 1985 .

[11]  T. Carrington,et al.  Vinylidene: Potential energy surface and unimolecular reaction dynamics , 1984 .

[12]  G. Scuseria,et al.  A systematic theoretical study of harmonic vibrational frequencies: The single and double excitation coupled cluster (CCSD) method , 1988 .

[13]  Timothy J. Lee.  Theory for externally contracted configuration interaction energy gradients , 1987 .

[14]  Dennis R. Salahub,et al.  Compact basis sets for LCAO‐LSD calculations. Part I: Method and bases for Sc to Zn , 1985 .

[15]  Tom Ziegler,et al.  The determination of molecular structures by density functional theory. The evaluation of analytical energy gradients by numerical integration , 1988 .

[16]  R. Bartlett.  Electronic Structure Methods. (Book Reviews: Comparison of Ab Initio Quantum Chemistry with Experiment for Small Molecules) , 1985 .

[17]  A. D. McLean,et al.  f-Type functions in the orbital basis for calculating molecular interactions involving d electrons Cr2 and Mo2 , 1983 .

[18]  H. Schaefer,et al.  A New dimension to quantum chemistry: Theoretical methods for the analytic evaluation of first, second, and third derivatives of the molecular electronic energy with respect to nuclear coordinates , 1986 .

[19]  M. Karplus,et al.  Analysis of Charge Distributions: Hydrogen Fluoride , 1964 .

[20]  John Arents,et al.  Atomic Structure Calculations , 1964 .

[21]  P. Pulay.  Comment on ‘‘Force in SCF theories’’ , 1983 .

[22]  Shigeru Obara,et al.  General recurrence formulas for molecular integrals over Cartesian Gaussian functions , 1988 .

[23]  J. Pople,et al.  Derivative studies in configuration–interaction theory , 1980 .

[24]  Poul Jørgensen,et al.  Geometrical derivatives of energy surfaces and molecular properties , 1986 .

[25]  Analytic energy derivative methods for excited singlet states of the same symmetry as the electronic ground state , 1985 .

[26]  D. Salahub,et al.  Molecular and electronic structure of ozone and thiozone from LCAO local density calculations , 1985 .

[27]  Peter Pulay,et al.  Ab initio calculation of force constants and equilibrium geometries in polyatomic molecules , 1969 .

[28]  Keiji Morokuma,et al.  Energy gradient in a multi-configurational SCF formalism and its application to geometry optimization of trimethylene diradicals , 1979 .

[29]  Ronald H. Felton,et al.  A new computational approach to Slater’s SCF–Xα equation , 1975 .

[30]  D. Salahub,et al.  Model potential calculations for second‐row transition metal molecules within the local‐spin‐density method , 1985 .

[31]  L. Hedin,et al.  A local exchange-correlation potential for the spin polarized case. i , 1972 .

[32]  W. Fink,et al.  Numerical Test of the Integral Hellmann—Feynman Theorem , 1967 .

[33]  Henry F. Schaefer,et al.  Gradient techniques for open‐shell restricted Hartree–Fock and multiconfiguration self‐consistent‐field methods , 1979 .

[34]  Shigeru Obara,et al.  Efficient recursive computation of molecular integrals over Cartesian Gaussian functions , 1986 .

[35]  P. Hohenberg,et al.  Inhomogeneous Electron Gas , 1964 .

[36]  L. Curtiss,et al.  Theoretical study of structures and energies of acetylene, ethylene, and vinyl radical and cation , 1988 .