Electrostatic interactions between molecules from relaxed one-electron density matrices of the coupled cluster singles and doubles model

The influence of electron correlation on the electrostatic interaction between closed shell molecules is studied using the relaxed electron densities of the coupled cluster singles and doubles (CCSD) model. The corresponding CCSD one-electron density matrices are efficiently computed without full four-index transformation by employing the generalized exchange and Coulomb operator technique. Using several representative van der Waals and hydrogen bonded complexes it was found that in most cases the convergence of the Møller-Plesset expansion of the electrostatic energy, restricted to single, double and quadruple excitations, is satisfactory and the fourth-order triple excitation term is more important than the sum of the fifth- and higher-order contributions from CCSD theory. The importance of the CCSD correlation correction to the electrostatic energy was gauged by comparison of the total interaction energy computed by symmetry-adapted perturbation theory (SAPT) and by the super-molecular CCSD(T) approach (coupled cluster singles and doubles model with a non-iterative inclusion of triple excitations). Except for the CO and N2 dimers, very good agreement between the two sets of results is observed. For the difficult case of the CO dimer the difference between the SAPT and CCSD(T) results can be explained by the truncation of the SAPT expansion for the dispersion energy at second order in the intramonomer correlation operator.

[1]  Hans-Joachim Werner,et al.  A comparison of the efficiency and accuracy of the quadratic configuration interaction (QCISD), coupled cluster (CCSD), and Brueckner coupled cluster (BCCD) methods , 1992 .

[2]  C. E. Dykstra Intermolecular electrical interaction: a key ingredient in hydrogen bonding , 1988 .

[3]  Tatiana Korona,et al.  Anisotropic intermolecular interactions in van der Waals and hydrogen-bonded complexes: What can we get from density functional calculations? , 1999 .

[4]  K. Szalewicz,et al.  Many‐body theory of exchange effects in intermolecular interactions. Second‐quantization approach and comparison with full configuration interaction results , 1994 .

[5]  Guntram Rauhut,et al.  Analytical energy gradients for local second-order Mo/ller–Plesset perturbation theory , 1998 .

[6]  Peter J. Knowles,et al.  Perturbative corrections to account for triple excitations in closed and open shell coupled cluster theories , 1994 .

[7]  Anthony J. Stone,et al.  Distributed multipole analysis, or how to describe a molecular charge distribution , 1981 .

[8]  S. F. Boys,et al.  The calculation of small molecular interactions by the differences of separate total energies. Some procedures with reduced errors , 1970 .

[9]  K. Szalewicz,et al.  Analytic first-order properties from explicitly correlated many-body perturbation theory and Gaussian geminal basis , 1998 .

[10]  K. Szalewicz,et al.  Effects of monomer geometry and basis set saturation on computed depth of water dimer potential , 1996 .

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

[12]  P. Wormer,et al.  Reply to the Comment on "The importance of high-order correlation effects for the CO-CO interaction potential" - [Chem. Phys. Lett. 314 (1999) 326] , 2001 .

[13]  H. Koch,et al.  Comment on ``The importance of high-order correlation effects for the CO CO interaction potential'' [Chem. Phys. Lett. 314 (1999) 326] , 2001 .

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

[15]  Jiří Čížek,et al.  Direct calculation of the Hartree–Fock interaction energy via exchange–perturbation expansion. The He … He interaction , 1987 .

[16]  P. Jørgensen,et al.  First-order one-electron properties in the integral-direct coupled cluster singles and doubles model , 1997 .

[17]  A. Stone,et al.  Electrostatic predictions of shapes and properties of Van der Waals molecules , 1986 .

[18]  K. Szalewicz,et al.  Symmetry-adapted perturbation theory calculation of the He-HF intermolecular potential energy surface , 1993 .

[19]  Hideto Kanamori,et al.  Ab initio MO studies of van der Waals molecule (N2)2: Potential energy surface and internal motion , 1998 .

[20]  Peter Pulay,et al.  An efficient reformulation of the closed‐shell self‐consistent electron pair theory , 1984 .

[21]  Rodney J. Bartlett,et al.  COUPLED-CLUSTER THEORY: AN OVERVIEW OF RECENT DEVELOPMENTS , 1995 .

[22]  K. Szalewicz,et al.  On the effectiveness of monomer‐, dimer‐, and bond‐centered basis functions in calculations of intermolecular interaction energies , 1995 .

[23]  Robert Moszynski,et al.  Perturbation Theory Approach to Intermolecular Potential Energy Surfaces of van der Waals Complexes , 1994 .

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

[25]  Sl,et al.  Many‐body theory of intermolecular induction interactions , 1994 .

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

[27]  K. Szalewicz,et al.  Symmetry-adapted double-perturbation analysis of intramolecular correlation effects in weak intermolecular interactions , 1979 .

[28]  Per Jensen,et al.  Computational molecular spectroscopy , 2000, Nature Reviews Methods Primers.

[29]  Raymond A. Poirier,et al.  Cumulative atomic multipole representation of the molecular charge distribution and its basis set dependence , 1983 .

[30]  A. van der Avoird,et al.  Symmetry‐adapted perturbation theory of nonadditive three‐body interactions in van der Waals molecules. I. General theory , 1995 .

[31]  P. Wormer,et al.  Ab Initio Potential Energy Surface and Infrared Spectrum of the Ne−CO Complex , 1997 .

[32]  Donald C. Comeau,et al.  The equation-of-motion coupled-cluster method. Applications to open- and closed-shell reference states , 1993 .

[33]  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 .

[34]  K. Szalewicz,et al.  Many‐body theory of exchange effects in intermolecular interactions. Density matrix approach and applications to He–F−, He–HF, H2–HF, and Ar–H2 dimers , 1994 .

[35]  Stanisl,et al.  Many‐body perturbation theory of electrostatic interactions between molecules: Comparison with full configuration interaction for four‐electron dimers , 1993 .

[36]  U. Kaldor,et al.  Many-Body Methods in Quantum Chemistry , 1989 .

[37]  Patrick W. Fowler,et al.  Theoretical studies of van der Waals molecules and intermolecular forces , 1988 .

[38]  A. J. Sadlej Perturbation theory of the electron correlation effects for atomic and molecular properties , 1981 .

[39]  Stanisl,et al.  Many‐body symmetry‐adapted perturbation theory of intermolecular interactions. H2O and HF dimers , 1991 .

[40]  T. Heijmen,et al.  Symmetry-adapted perturbation theory for the calculation of Hartree-Fock interaction energies , 1996 .

[41]  K. Szalewicz,et al.  Many-Body Theory of Van der Waals Interactions , 1989 .

[42]  D. Yarkony,et al.  Modern Electronic Structure Theory: Part I , 1995 .

[43]  K. Szalewicz,et al.  Perturbation theory calculations of intermolecular interaction energies , 1991 .

[44]  T. Heijmen,et al.  Ab initio potential-energy surface and rotationally inelastic integral cross sections of the Ar–CH4 complex , 1997 .

[45]  K. Szalewicz,et al.  Møller–Plesset expansion of the dispersion energy in the ring approximation , 1993 .

[46]  John F. Stanton,et al.  The ACES II program system , 1992 .

[47]  Jeppe Olsen,et al.  On the divergent behavior of Møller–Plesset perturbation theory for the molecular electric dipole moment , 2000 .

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

[49]  J. V. Lenthe,et al.  State of the Art in Counterpoise Theory , 1994 .

[50]  Rodney J. Bartlett,et al.  Theory and application of MBPT(3) gradients: The density approach , 1987 .

[51]  P. Wormer,et al.  The importance of high-order correlation effects for the CO–CO interaction potential , 1999 .

[52]  V. Kellö,et al.  Medium-size polarized basis sets for high-level-correlated calculations of molecular electric properties , 1991 .

[53]  J. Olsen,et al.  Divergence in Møller–Plesset theory: A simple explanation based on a two-state model , 2000 .

[54]  K. Szalewicz,et al.  Helium dimer potential from symmetry-adapted perturbation theory calculations using large Gaussian geminal and orbital basis sets , 1997 .