Second order perturbation corrections to singles and doubles coupled-cluster methods: General theory and application to the valence optimized doubles model

We present a general perturbative method for correcting a singles and doubles coupled-cluster energy. The coupled-cluster wave function is used to define a similarity-transformed Hamiltonian, which is partitioned into a zeroth-order part that the reference problem solves exactly plus a first-order perturbation. Standard perturbation theory through second-order provides the leading correction. Applied to the valence optimized doubles ~VOD! approximation to the full-valence complete active space self-consistent field method, the second-order correction, which we call ~2!, captures dynamical correlation effects through external single, double, and semi-internal triple and quadruple substitutions. A factorization approximation reduces the cost of the quadruple substitutions to only sixth order in the size of the molecule. A series of numerical tests are presented showing that VOD~2! is stable and well-behaved provided that the VOD reference is also stable. The second-order correction is also general to standard unwindowed coupled-cluster energies such as the coupled-cluster singles and doubles ~CCSD! method itself, and the equations presented here fully define the corresponding CCSD~2! energy. © 2000 American Institute of Physics. @S0021-9606~00!30130-1#

[1]  R. Bartlett,et al.  An efficient way to include connected quadruple contributions into the coupled cluster method , 1998 .

[2]  John F. Stanton,et al.  INVESTIGATION OF AN ASYMMETRIC TRIPLE-EXCITATION CORRECTION FOR COUPLED-CLUSTER ENERGIES , 1998 .

[3]  Nicholas C. Handy,et al.  Size-consistent Brueckner theory limited to double substitutions , 1989 .

[4]  John F. Stanton,et al.  The equation of motion coupled‐cluster method. A systematic biorthogonal approach to molecular excitation energies, transition probabilities, and excited state properties , 1993 .

[5]  M. Head‐Gordon,et al.  An accurate local model for triple substitutions in fourth order Møller–Plesset theory and in perturbative corrections to singles and doubles coupled cluster methods , 2000 .

[6]  Anna I. Krylov,et al.  Energies and analytic gradients for a coupled-cluster doubles model using variational Brueckner orbitals: Application to symmetry breaking in O4+ , 1998 .

[7]  John F. Stanton,et al.  A simple correction to final state energies of doublet radicals described by equation-of-motion coupled cluster theory in the singles and doubles approximation , 1997 .

[8]  R. Bartlett,et al.  GEOMETRY AND HARMONIC FREQUENCY OF N2 WITH COUPLED CLUSTER METHODS THAT INCLUDE CONNECTED QUADRUPLE EXCITATIONS , 1999 .

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

[10]  J. Stanton Why CCSD(T) works: a different perspective , 1997 .

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

[12]  John F. Stanton,et al.  Many-body perturbation theory with a restricted open-shell Hartree—Fock reference , 1991 .

[13]  Martin Head-Gordon,et al.  A tensor formulation of many-electron theory in a nonorthogonal single-particle basis , 1998 .

[14]  G. Scuseria,et al.  The optimization of molecular orbitals for coupled cluster wavefunctions , 1987 .

[15]  J. Stanton,et al.  Application of an equation-of-motion coupled cluster method including higher-order corrections to potential energy surfaces of radicals , 1999 .

[16]  B. Roos,et al.  Lecture notes in quantum chemistry , 1992 .

[17]  R. Bartlett,et al.  Connected quadruples for the frequencies of O3 , 1999 .

[18]  Björn O. Roos,et al.  Second-order perturbation theory with a complete active space self-consistent field reference function , 1992 .

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

[20]  Rodney J. Bartlett,et al.  Full configuration-interaction and state of the art correlation calculations on water in a valence double-zeta basis with polarization functions , 1996 .

[21]  R. Bartlett,et al.  A full coupled‐cluster singles and doubles model: The inclusion of disconnected triples , 1982 .

[22]  Kerstin Andersson,et al.  Second-order perturbation theory with a CASSCF reference function , 1990 .

[23]  H. Koch,et al.  Gauge invariant coupled cluster response theory , 1999 .

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

[25]  P. Löwdin Studies in Perturbation Theory. IV. Solution of Eigenvalue Problem by Projection Operator Formalism , 1962 .

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

[27]  Martin Head-Gordon,et al.  A second-order correction to singles and doubles coupled-cluster methods based on a perturbative expansion of a similarity-transformed Hamiltonian , 2000 .

[28]  T. H. Dunning Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen , 1989 .

[29]  Rodney J. Bartlett,et al.  Noniterative energy corrections through fifth-order to the coupled cluster singles and doubles method , 1998 .

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

[31]  N. Handy,et al.  Dynamical and Nondynamical Correlation , 1996 .

[32]  Anna I. Krylov,et al.  Size-consistent wave functions for nondynamical correlation energy: The valence active space optimized orbital coupled-cluster doubles model , 1998 .