Analytic evaluation of energy gradients for the singles and doubles coupled cluster method including perturbative triple excitations: Theory and applications to FOOF and Cr2

The analytic energy gradient for the singles and doubles coupled cluster method including a perturbative correction due to triple excitations [CCSD(T)] is formulated and computationally implemented. Encouraged by the recent success in reproducing the experimental equilibrium structure and vibrational frequencies of ozone, the new CCSD(T) gradient method is tested with two other ‘‘difficult’’ quantum chemistry problems: FOOF and Cr2. With the largest basis set employed in this work [triple zeta plus two sets of polarization functions (TZ2Pf)] at the CCSD(T) level of theory, the predictions for the O–O and O–F bond lengths in FOOF are 1.218 and 1.589 A, respectively. These figures are in good agreement with the experimental values 1.216 and 1.575 A. Based on CCSD calculations with even larger basis sets, it is concluded that the error of 0.014 A in the O–F bond length at the TZ2Pf/CCSD(T) level of theory is due to the remaining basis set deficiency. On the other hand, the CCSD(T) prediction for the equilibr...

[1]  G. Scuseria,et al.  Comparison of coupled-cluster methods which include the effects of connected triple excitations , 1990 .

[2]  M. Plesset,et al.  Note on an Approximation Treatment for Many-Electron Systems , 1934 .

[3]  Krishnan Raghavachari,et al.  Highly correlated systems: Structure, binding energy and harmonic vibrational frequencies of ozone , 1989 .

[4]  B. Delley,et al.  Metal-metal bonding in Cr-Cr and Mo-Mo dimers: Another success of local spin-density theory , 1983 .

[5]  P. Taylor,et al.  Vibrations in small Mg clusters , 1990 .

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

[7]  H. Schaefer,et al.  Fluorine peroxide (FOOF): A continuing problem for normally reliable theoretical methods , 1987 .

[8]  M. Head‐Gordon,et al.  A fifth-order perturbation comparison of electron correlation theories , 1989 .

[9]  Henry F. Schaefer,et al.  Ordering of the O-O stretching vibrational frequencies in ozone , 1989 .

[10]  P. Taylor,et al.  The structure of F2O2: Theoretical predictions and comparisons with F2 and F2O , 1982 .

[11]  Gustavo E. Scuseria,et al.  The vibrational frequencies of ozone , 1990 .

[12]  J. Bernholc,et al.  Local Spin-Density Description of Multiple Metal-Metal Bonding: Mo 2 and Cr 2 , 1983 .

[13]  R. Bartlett,et al.  The full CCSDT model for molecular electronic structure , 1987 .

[14]  Trygve Helgaker,et al.  Molecular Hessians for large‐scale MCSCF wave functions , 1986 .

[15]  Henry F. Schaefer,et al.  A new implementation of the full CCSDT model for molecular electronic structure , 1988 .

[16]  R. Bartlett,et al.  A coupled cluster approach with triple excitations , 1984 .

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

[18]  M. Page,et al.  Multireference CI Gradients and MCSCF Second Derivatives. , 1984 .

[19]  K. Richman,et al.  Numerical Hartree–Fock calculations on diatomic chromium , 1987 .

[20]  R. Ahlrichs,et al.  Many body perturbation calculations and coupled electron pair models , 1979 .

[21]  Rodney J. Bartlett,et al.  Analytic energy derivatives in many‐body methods. I. First derivatives , 1989 .

[22]  R. O. Jones,et al.  Density functional theory and molecular bonding. III. Iron‐series dimers , 1979 .

[23]  Ron Shepard,et al.  Geometrical energy derivative evaluation with MRCI wave functions , 1987 .

[24]  W. D. Allen,et al.  The analytic evaluation of energy first derivatives for two‐configuration self‐consistent‐field configuration interaction (TCSCF‐CI) wave functions. Application to ozone and ethylene , 1987 .

[25]  G. Scuseria,et al.  Analytic evaluation of energy gradients for the single, double and linearized triple excitation coupled-cluster CCSDT-1 wavefunction: Theory and applications , 1988 .

[26]  G. Scuseria,et al.  Is coupled cluster singles and doubles (CCSD) more computationally intensive than quadratic configuration interaction (QCISD) , 1989 .

[27]  J. Gauss,et al.  Analytical differentiation of the energy contribution due to triple excitations in fourth-order Møller-Plesset perturbation theory , 1988 .

[28]  Peter J. Knowles,et al.  Studies using the CASSCF wavefunction , 1982 .

[29]  R. Bartlett,et al.  The equilibrium structure and harmonic vibrational frequencies of ozone: coupled cluster results including triple excitations , 1989 .

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

[31]  David H. Magers,et al.  Highly correlated single‐reference studies of the O3 potential surface. I. Effects of high order excitations on the equilibrium structure and harmonic force field of ozone , 1989 .

[32]  Henry F. Schaefer,et al.  On the evaluation of analytic energy derivatives for correlated wave functions , 1984 .

[33]  Josef Paldus,et al.  Coupled cluster approach or quadratic configuration interaction , 1989 .

[34]  N. Handy,et al.  The analytic gradient for the coupled pair functional method: Formula and application for HCl, H2CO, and the dimer H2CO⋅⋅⋅HCl , 1988 .

[35]  R. Kok,et al.  A theoretical investigation of the bond length of dichromium , 1983 .

[36]  T. H. Dunning.  Gaussian Basis Functions for Use in Molecular Calculations. III. Contraction of (10s6p) Atomic Basis Sets for the First‐Row Atoms , 1970 .

[37]  M. Guest,et al.  Theoretical study of the electronic structure of the transition metal dimers, Sc2, Cr2, Mo2 and Ni2 , 1980 .

[38]  D. Cremer,et al.  Analytical evaluation of energy gradients in quadratic configuration interaction theory , 1988 .

[39]  F. Cotton,et al.  Dimolybdenum: nature of the sextuple bond , 1980 .

[40]  W. Goddard,et al.  Nature of Mo-Mo and Cr-Cr Multiple Bonds: A Challenge for the Local-Density Approximation , 1982 .

[41]  D. Salahub,et al.  LCAO local-spin-density and Xα calculations for Cr2 and Mo2 , 1984 .

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

[43]  Julia E. Rice,et al.  Analytic evaluation of energy gradients for the single and double excitation coupled cluster (CCSD) wave function: Theory and application , 1987 .

[44]  S. J. Cole,et al.  Towards a full CCSDT model for electron correlation , 1985 .

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

[46]  Analytical MBPT(4) gradients , 1988 .

[47]  B. Roos,et al.  Theoretical evidence for multiple 3d bondig in the V2 and Cr2 molecules , 1983 .

[48]  P. Jeffrey Hay,et al.  The structure and bonding of F2O2 , 1987 .

[49]  K. Raghavachari,et al.  Highly correlated systems. Structure and harmonic force field of F2O2 , 1989 .

[50]  H. Schaefer,et al.  The analytic configuration interaction gradient method: Application to the cyclic and open isomers of the S3 molecule , 1986 .

[51]  A. Wachters,et al.  Gaussian Basis Set for Molecular Wavefunctions Containing Third‐Row Atoms , 1970 .

[52]  B. Dunlap.  Xα,Cr2, and the symmetry dilemma , 1983 .

[53]  Peter R. Taylor,et al.  Analytical MCSCF energy gradients: Treatment of symmetry and CASSCF applications to propadienone , 1984 .

[54]  Michael A. Robb,et al.  MC SCF gradient optimization of the H2CO→H2 + CO transition structure , 1982 .

[55]  Curtis L. Janssen,et al.  An efficient reformulation of the closed‐shell coupled cluster single and double excitation (CCSD) equations , 1988 .

[56]  R. Bartlett,et al.  Analytic energy gradients for general coupled‐cluster methods and fourth‐order many‐body perturbation theory , 1986 .

[57]  Alistair P. Rendell,et al.  Comparison of the quadratic configuration interaction and coupled cluster approaches to electron correlation including the effect of triple excitations , 1990 .

[58]  Peter R. Taylor,et al.  General contraction of Gaussian basis sets. I. Atomic natural orbitals for first‐ and second‐row atoms , 1987 .

[59]  H. Schaefer,et al.  Analytic evaluation of infrared intensities and polarizabilities by two-configuration self-consistent field wave functions , 1986 .

[60]  A. Wolf,et al.  Nonempirical calculations on diatomic transition metals. II. RHF investigation of lowest closed‐shell states of homonuclear 3d transition‐metal dimers , 1980 .

[61]  H. Schaefer,et al.  Theoretical investigations of molecules composed only of fluorine, oxygen and nitrogen: determination of the equilibrium structures of FOOF, (NO)2 and FNNF and the transition state structure for FNNF cis-trans isomerization , 1989 .

[62]  I. H. Hillier,et al.  Correlation effects and the bonding in Mo2 and Cr2 , 1982 .

[63]  Charles W. Bauschlicher,et al.  Theoretical calculation of ozone vibrational infrared intensities , 1985 .

[64]  Martin Head-Gordon,et al.  Quadratic configuration interaction. A general technique for determining electron correlation energies , 1987 .

[65]  H. Oberhammer,et al.  An improved ab initio structure for fluorine peroxide (FOOF) , 1988 .

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