Scope and limitations of the SCS-MP2 method for stacking and hydrogen bonding interactions.

Fluorobenzenes are pi-acceptor synthons that form pi-stacked structures in molecular crystals as well as in artificial DNAs. We investigate the competition between hydrogen bonding and pi-stacking in dimers consisting of the nucleobase mimic 2-pyridone (2PY) and all fluorobenzenes from 1-fluorobenzene to hexafluorobenzene (n-FB, with n = 1-6). We contrast the results of high level ab initio calculations with those obtained using ultraviolet (UV) and infrared (IR) laser spectroscopy of isolated and supersonically cooled dimers. The 2PY.n-FB complexes with n = 1-5 prefer double hydrogen bonding over pi-stacking, as diagnosed from the UV absorption and IR laser depletion spectra, which both show features characteristic of doubly H-bonded complexes. The 2-pyridone.hexafluorobenzene dimer is the only pi-stacked dimer, exhibiting a homogeneously broadened UV spectrum and no IR bands characteristic for H-bonded species. MP2 (second-order Møller-Plesset perturbation theory) calculations overestimate the pi-stacked dimer binding energies by about 10 kJ/mol and disagree with the experimental observations. In contrast, the MP2 treatment of the H-bonded dimers appears to be quite accurate. Grimme's spin-component-scaled MP2 approach (SCS-MP2) is an improvement over MP2 for the pi-stacked dimers, reducing the binding energy by approximately 10 kJ/mol. When applied to explicitly correlated MP2 theory (SCS-MP2-R12 approach), agreement with the corresponding coupled-cluster binding energies [at the CCSD(T) level] is very good for the pi-stacked dimers, within +/- 1 kJ/mol for the 2PY complexes with 1-fluorobenzene, 1,2-difluorobenzene, 1,2,4,5-tetrafluorobenzene, pentafluorobenzene and hexafluorobenzene. Unfortunately, the SCS-MP2 approach also reduces the binding energy of the H-bonded species, leading to disagreement with both coupled-cluster theory and experiment. The SCS-MP2-R12 binding energies follow the SCS-MP2 binding energies closely, being about 0.5 and 0.7 kJ/mol larger for the H-bonded and pi-stacked forms, respectively, in an augmented correlation-consistent polarized valence quadruple-zeta basis. It seems that the SCS-MP2 and SCS-MP2-R12 methods cannot provide sufficient accuracy to replace the CCSD(T) method for intermolecular interactions where H-bonding and pi-stacking are competitive.

[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]  S. Grimme,et al.  Is spin-component scaled second-order Møller-Plesset perturbation theory an appropriate method for the study of noncovalent interactions in molecules? , 2007, The journal of physical chemistry. A.

[3]  T. Ebata,et al.  Vibrational spectroscopy of 2-pyridone and its clusters in supersonic jets: Structures of the clusters as revealed by characteristic shifts of the NH and C=O bands , 1999 .

[4]  P. Felker,et al.  Stimulated Raman spectroscopy in the ν1 region of isotopically substituted benzene dimers: evidence for symmetrically inequivalent benzene moieties , 1991 .

[5]  Martin Head-Gordon,et al.  Scaled opposite-spin second order Møller-Plesset correlation energy: an economical electronic structure method. , 2004, The Journal of chemical physics.

[6]  S. Leutwyler,et al.  2-pyridone: The role of out-of-plane vibrations on the S1<-->S0 spectra and S1 state reactivity. , 2006, The Journal of chemical physics.

[7]  Florian Weigend,et al.  A fully direct RI-HF algorithm: Implementation, optimised auxiliary basis sets, demonstration of accuracy and efficiency , 2002 .

[8]  Edward F. Valeev.  Improving on the resolution of the identity in linear R12 ab initio theories , 2004 .

[9]  Yasuhito Aoki,et al.  S1–S0 vibronic spectra of benzene clusters revisited. I. The tetramer , 2002 .

[10]  W. Kutzelnigg,et al.  Møller-plesset calculations taking care of the correlation CUSP , 1987 .

[11]  Christof Hättig,et al.  Geometry optimizations with the coupled-cluster model CC2 using the resolution-of-the-identity approximation , 2003 .

[12]  P. Felker,et al.  Mass‐selective ionization‐detected stimulated Raman spectroscopy of benzene trimer and higher clusters , 1993 .

[13]  William Klemperer,et al.  Molecular beam studies of benzene dimer, hexafluorobenzene dimer, and benzene–hexafluorobenzene , 1979 .

[14]  Frederick R. Manby,et al.  R12 methods in explicitly correlated molecular electronic structure theory , 2006 .

[15]  Christof Hättig,et al.  Optimization of auxiliary basis sets for RI-MP2 and RI-CC2 calculations: Core–valence and quintuple-ζ basis sets for H to Ar and QZVPP basis sets for Li to Kr , 2005 .

[16]  Krzysztof Szalewicz,et al.  Potential energy surface for the benzene dimer and perturbational analysis of π-π interactions , 2006 .

[17]  Christof Hättig,et al.  Quintuple-ζ quality coupled-cluster correlation energies with triple-ζ basis sets , 2007 .

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

[19]  Seiichiro Ten-no,et al.  Initiation of explicitly correlated Slater-type geminal theory , 2004 .

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

[21]  T. Takatani,et al.  Performance of spin-component-scaled Møller-Plesset theory (SCS-MP2) for potential energy curves of noncovalent interactions. , 2007, Physical chemistry chemical physics : PCCP.

[22]  T. Dunning,et al.  Electron affinities of the first‐row atoms revisited. Systematic basis sets and wave functions , 1992 .

[23]  Roland Lindh,et al.  The reduced multiplication scheme of the Rys quadrature and new recurrence relations for auxiliary function based two‐electron integral evaluation , 1991 .

[24]  Martin Head-Gordon,et al.  Optimized spin-component scaled second-order Møller-Plesset perturbation theory for intermolecular interaction energies , 2007 .

[25]  Hans-Joachim Werner,et al.  Calculation of intermolecular interactions in the benzene dimer using coupled-cluster and local electron correlation methods. , 2006, Physical chemistry chemical physics : PCCP.

[26]  Roman Leist,et al.  Hydrogen bonding of the nucleobase mimic 2-pyridone to fluorobenzenes: an ab initio investigation. , 2006, The journal of physical chemistry. A.

[27]  Frederick R. Manby,et al.  Density fitting in second-order linear-r12 Møller–Plesset perturbation theory , 2003 .

[28]  Masuhiro Mikami,et al.  Energy profile of the interconversion path between T-shape and slipped-parallel benzene dimers , 2002 .

[29]  W. Klopper,et al.  Nucleobase-fluorobenzene interactions: hydrogen bonding wins over pi stacking. , 2007, Angewandte Chemie.

[30]  C. David Sherrill,et al.  High-Accuracy Quantum Mechanical Studies of π−π Interactions in Benzene Dimers , 2006 .

[31]  Hans W. Horn,et al.  ELECTRONIC STRUCTURE CALCULATIONS ON WORKSTATION COMPUTERS: THE PROGRAM SYSTEM TURBOMOLE , 1989 .

[32]  James A Platts,et al.  Spin-Component Scaling Methods for Weak and Stacking Interactions. , 2007, Journal of chemical theory and computation.

[33]  T. Ebata,et al.  Population labeling spectroscopy for the electronic and the vibrational transitions of 2-pyridone and its hydrogen-bonded clusters , 2000 .

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

[35]  Elangannan Arunan,et al.  The rotational spectrum, structure and dynamics of a benzene dimer , 1993 .

[36]  P. Felker,et al.  The Raman and vibronic activity of intermolecular vibrations in aromatic‐containing complexes and clusters , 1994 .

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

[38]  P. Felker,et al.  RAMAN-VIBRONIC DOUBLE RESONANCE SPECTROSCOPY OF BENZENE DIMER , 1992 .

[39]  Christof Hättig,et al.  CC2 excitation energy calculations on large molecules using the resolution of the identity approximation , 2000 .

[40]  F. Weigend,et al.  Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: Design and assessment of accuracy. , 2005, Physical chemistry chemical physics : PCCP.

[41]  S. Grimme.  Improved second-order Møller–Plesset perturbation theory by separate scaling of parallel- and antiparallel-spin pair correlation energies , 2003 .

[42]  J. Noga,et al.  Alternative formulation of the matrix elements in MP2‐R12 theory , 2005 .

[43]  D. Tew,et al.  New correlation factors for explicitly correlated electronic wave functions. , 2005, The Journal of chemical physics.

[44]  Seiichiro Ten-no,et al.  Explicitly correlated second order perturbation theory: introduction of a rational generator and numerical quadratures. , 2004, The Journal of chemical physics.

[45]  Wim Klopper,et al.  Explicitly correlated second-order Møller–Plesset methods with auxiliary basis sets , 2002 .