Hydrogen bonding of the nucleobase mimic 2-pyridone to fluorobenzenes: an ab initio investigation

The hydrogen-bonded complexes of the nucleobase mimic 2-pyridone (2PY) with seven different fluorinated benzenes (1-, 1,2-, 1,4-, 1,2,3-, 1,3,5-, 1,2,3,4-, and 1,2,4,5-fluorobenzene) are important model systems for investigating the relative importance of hydrogen bonding versus pi-stacking interactions in DNA. We have shown by supersonic-jet spectroscopy that these dimers are hydrogen bonded and not pi-stacked at low temperature (Leist, R.; Frey, J. A.; Leutwyler, S. J. Phys. Chem. A 2006, 110, 4180). Their geometries and binding energies D(e) were calculated using the resolution of identity (RI) M?ller-Plesset second-order perturbation theory method (RIMP2). The most stable dimers are bound by antiparallel N-H...F-C and C-H...O=C hydrogen bonds. The binding energies are extrapolated to the complete basis set (CBS) limit, , using the aug-cc-pVXZ basis set series. The CBS binding energies range from -D(e,CBS) = 6.4-6.9 kcal/mol and the respective dissociation energies from -D(0,CBS) = 5.9-6.3 kcal/mol. In combination with experiment, the latter represent upper limits to the dissociation energies of the pi-stacked isomers (which are not observed experimentally). The individual C-H...O=C and N-H...F-C contributions to D(e) can be approximately separated. They are nearly equal for 2PY.fluorobenzene; each additional F atom strengthens the C-H...O=C hydrogen bond by approximately 0.5 kcal/mol and weakens the C-F...H-N hydrogen bond by approximately 0.3 kcal/mol. The single H-bond strengths and lengths correlate with the gas-phase acid-base properties of the C-H and C-F groups of the fluorobenzenes.     

Is spin-component scaled second-order Møller-Plesset perturbation theory an appropriate method for the study of noncovalent interactions in molecules?

Testing of the spin-component scaled second-order M?ller-Plesset (SCS-MP2) method for the computation of noncovalent interaction energies is done with a database of 165 biologically relevant complexes. The effects of the spin-scaling procedure (i.e., MP2 vs SCS-MP2), the basis set size, and the corrections for basis set superposition error (BSSE) are systematically examined. When using two-point basis set extrapolations for the correlation energy, augmentation of the atomic orbital basis with computationally costly diffuse functions is found to be obsolete. In general, SCS-MP2 also improves results for noncovalent interactions statistically on MP2, and significant outliers are removed. Moreover, it is shown that effects of BSSE and one-particle basis set incompleteness almost cancel each other in the case of triple-zeta sets (SCS-MP2/TZVPP or SCS-MP2/cc-pVTZ without counterpoise correction), which opens a practical route to efficient computations for large systems. We recommend SCS-MP2 as the preferred quantum chemical wave function based method for the noncovalent interactions in large biologically relevant systems when reasonable coupled-cluster with single and double and perturbative triple excitations (CCSD(T)) calculations cannot be performed anymore. A comparison to MP2 and CCSD(T) interaction energies for n-alkane dimers, however, indicates (and this also holds to a lesser extent for hydrogen-bonded systems) limitations of SCS-MP2 when treating chemically "saturated" interactions. The different behavior of second-order perturbation theory for saturated and for stacked pi-systems is discussed.     

Intermolecular interactions in nitrogen-containing aromatic systems

π–π and CH···N interactions are vital in biological systems. In this study, stacking and hydrogen-bonded interactions in pyrazine and triazine dimers were investigated by density functional theory combined with symmetry-adapted perturbation theory (DFT-SAPT) and counterpoise (CP)-corrected supermolecular MP2, SCS-MP2, B3LYP-D and CCSD(T) calculations. All interaction energies were computed using the optimized structures at the CP-corrected SCS/aug-cc-pVDZ level, which gave 1–2 kJ/mol lower interaction energies than the ones computed at the MP2 level. For both dimers, doubly hydrogen-bonded and cross-(displaced) stacked orientations were found to be the lowest energy ones. The reference CCSD(T) calculations favored the former structure in both dimer systems, whereas MP2 and SCS-MP2 located the latter as the lowest energy isomer. In particular, the former was found to be lower in energy than the latter by 2.28 and 1.01 kJ/mol at the CCSD(T)/aug-cc-pVDZ level for pyrazine and triazine, respectively. B3LYP-D produced interaction energies in agreement with the CCSD(T) at the equilibrium geometries, but it overestimates them at the short range and underestimates at the long intermonomer separations. Furthermore, it tends to give smaller equilibrium distances compared to the CCSD(T). DFT-SAPT method was in a good agreement with the reference CCSD(T) calculations. This suggests that DFT-SAPT can be employed to compute the full potential energy surface of these dimers. Moreover, DFT-SAPT calculations showed that the electrostatic and dispersion contributions are the most important energy components stabilizing these dimers. The present study aims to show which theoretical method is the most promising one for the investigation of intermolecular interactions dominated by π–π and CH···N. Therefore, the findings obtained in this study can be used to unravel the structures of nucleic acid bases and other systems stabilized by π–π and CH···N interactions.     

Physical nature of silane⋯carbene dimers revealed by state-of-the-art ab initio calculations

Using the SAPT2 + 3(CCD)δMP2 method in complete basis set (CBS) limit, it is shown that the interactions in the recently studied silane⋯carbene dimers are mainly dispersive in nature. Consequently, slow convergence of dispersion energy also forces slow convergence of the interaction energy. Therefore, obtaining very accurate values requires extrapolation of the correlation part to the CBS limit. The most accurate values obtained at the CCSD(T)/CBS level of theory show that the studied silane⋯carbene dimers are rather weakly bound, with interaction energies ranging from about −1.9 to −1.3 kcal/mol. Comparing to CCSD(T)/CBS, it will be shown that SCS-MP2 and MP2C methods clearly underestimate and methods based on SAPT2+ and having some third-order corrections, as well as the MP2 method, overestimate values of interaction energies. Popular SAPT(DFT) method performs better than SCS-MP2 and MP2C; nevertheless, underestimation is still considerable. The underestimation is slightly quenched if third-order dispersion energy and its exchange counterpart is added to the SAPT(DFT). The closest value of CCSD(T)/CBS has been given by the SAPT2 + (3)(CCD)δMP2 method in quadruple-ζ basis set. © 2019 Wiley Periodicals, Inc.     

Probing phenylalanine/adenine pi-stacking interactions in protein complexes with explicitly correlated and CCSD(T) computations

To examine the effects of pi-stacking interactions between aromatic amino acid side chains and adenine bearing ligands in crystalline protein structures, 26 toluene/(N9-methyl)adenine model configurations have been constructed from protein/ligand crystal structures. Full geometry optimizations with the MP2 method cause the 26 crystal structures to collapse to six unique structures. The complete basis set (CBS) limit of the CCSD(T) interaction energies has been determined for all 32 structures by combining explicitly correlated MP2-R12 computations with a correction for higher-order correlation effects from CCSD(T) calculations. The CCSD(T) CBS limit interaction energies of the 26 crystal structures range from -3.19 to -6.77 kcal mol (-1) and average -5.01 kcal mol (-1). The CCSD(T) CBS limit interaction energies of the optimized complexes increase by roughly 1.5 kcal mol (-1) on average to -6.54 kcal mol (-1) (ranging from -5.93 to -7.05 kcal mol (-1)). Corrections for higher-order correlation effects are extremely important for both sets of structures and are responsible for the modest increase in the interaction energy after optimization. The MP2 method overbinds the crystal structures by 2.31 kcal mol (-1) on average compared to 4.50 kcal mol (-1) for the optimized structures.     

Hybrid density functional theory for pi-stacking interactions: application to benzenes, pyridines, and DNA bases

The suitability of a hybrid density functional to qualitatively reproduce geometric and energetic details of parallel pi-stacked aromatic complexes is presented. The hybrid functional includes an ad hoc mixture of half the exact (HF) exchange with half of the uniform electron gas exchange, plus Lee, Yang, and Parr's expression for correlation energy. This functional, in combination with polarized, diffuse basis sets, gives a binding energy for the parallel-displaced benzene dimer in good agreement with the best available high-level calculations reported in the literature, and qualitatively reproduces the local MP2 potential energy surface of the parallel-displaced benzene dimer. This method was further critically compared to high-level calculations recently reported in the literature for a range of pi-stacked complexes, including monosubstituted benzene-benzene dimers, along with DNA and RNA bases, and generally agrees with MP2 and/or CCSD(T) results to within +/-2 kJ mol(-1). We also show that the resulting BH&H binding energy is closely related to the electron density in the intermolecular region. The net result is that the BH&H functional, presumably due to fortuitous cancellation of errors, provides a pragmatic, computationally efficient quantum mechanical tool for the study of large pi-stacked systems such as DNA.     

Fluorobenzene-nucleobase interactions: hydrogen bonding or pi-stacking?

Studies on modified DNA oligomers and polymerase reactions have previously demonstrated that canonical nucleobases can exhibit stable and even selective pairing with shape-complementary fluorobenzene nucleotides. Because of the fluorination of the pairing edges, hydrogen bonds are believed to be absent, and the local DNA stability has been attributed to pi-stacking and shape complementarity. Using two-color resonant two-photon ionization and fluorescence emission spectroscopies, we show here that supersonically cooled complexes of the nucleobase analogue 2-pyridone with seven substituted fluorobenzenes (1-fluorobenzene, 1,2- and 1,4-difluorobenzene, 1,3,5- and 1,2,3-trifluorobenzene, 1,2,4,5- and 1,2,3,4-tetrafluorobenzene) are hydrogen-bonded and not pi-stacked. The S1 <--> S0 vibronic spectra show intermolecular vibrational frequencies that are characteristic for doubly hydrogen bonded complexes. The bands shift to the blue with increasing hydrogen-bond strength; the measured spectral blue shifts deltanu are in excellent agreement with the ab initio calculated shifts. The spectral shifts are also linearly correlated with the calculated hydrogen-bond dissociation energies D0, published in a companion paper (Frey, J. A.; Leist, R.; Leutwyler, S. J. Phys. Chem. A 2006, 110, 4188). This correlation allows us to reliably estimate the ground-state dissociation energies as D0 approximately 6 kcal/mol of the 2-pyridone.fluorobenzene complexes from the observed spectral shifts.     

Toward true DNA base-stacking energies: MP2, CCSD(T), and complete basis set calculations

Stacking energies in low-energy geometries of pyrimidine, uracil, cytosine, and guanine homodimers were determined by the MP2 and CCSD(T) calculations utilizing a wide range of split-valence, correlation-consistent, and bond-functions basis sets. Complete basis set MP2 (CBS MP2) stacking energies extrapolated using aug-cc-pVXZ (X = D, T, and for pyrimidine dimer Q) basis sets equal to -5.3, -12.3, and -11.2 kcal/mol for the first three dimers, respectively. Higher-order correlation corrections estimated as the difference between MP2 and CCSD(T) stacking energies amount to 2.0, 0.7, and 0.9 kcal/mol and lead to final estimates of the genuine stacking energies for the three dimers of -3.4, -11.6, and -10.4 kcal/mol. The CBS MP2 stacking-energy estimate for guanine dimer (-14.8 kcal/mol) was based on the 6-31G(0.25) and aug-cc-pVDZ calculations. This simplified extrapolation can be routinely used with a meaningful accuracy around 1 kcal/mol for large aromatic stacking clusters. The final estimate of the guanine stacking energy after the CCSD(T) correction amounts to -12.9 kcal/mol. The MP2/6-31G(0.25) method previously used as the standard level to calculate aromatic stacking in hundreds of geometries of nucleobase dimers systematically underestimates the base stacking by ca. 1.0-2.5 kcal/mol per stacked dimer, covering 75-90% of the intermolecular correlation stabilization. We suggest that this correction is to be considered in calibration of force fields and other cheaper computational methods. The quality of the MP2/6-31G(0.25) predictions is nevertheless considerably better than suggested on the basis of monomer polarizability calculations. Fast and very accurate estimates of the MP2 aromatic stacking energies can be achieved using the RI-MP2 method. The CBS MP2 calculations and the CCSD(T) correction, when taken together, bring only marginal changes to the relative stability of H-bonded and stacked base pairs, with a slight shift of ca. 1 kcal/mol in favor of H-bonding. We suggest that the present values are very close to ultimate predictions of the strength of aromatic base stacking of DNA and RNA bases.     

Probing the effects of heterogeneity on delocalized pi...pi interaction energies

Dimers composed of benzene (Bz), 1,3,5-triazine (Tz), cyanogen (Cy) and diacetylene (Di) are used to examine the effects of heterogeneity at the molecular level and at the cluster level on pi...pi stacking energies. The MP2 complete basis set (CBS) limits for the interaction energies (E(int)) of these model systems were determined with extrapolation techniques designed for correlation consistent basis sets. CCSD(T) calculations were used to correct for higher-order correlation effects (deltaE(CCSD)(T)(MP2)) which were as large as +2.81 kcal mol(-1). The introduction of nitrogen atoms into the parallel-slipped dimers of the aforementioned molecules causes significant changes to E(int). The CCSD(T)/CBS E(int) for Di-Cy is -2.47 kcal mol(-1) which is substantially larger than either Cy-Cy (-1.69 kcal mol(-1)) or Di-Di (-1.42 kcal mol(-1)). Similarly, the heteroaromatic Bz-Tz dimer has an E(int) of -3.75 kcal mol(-1) which is much larger than either Tz-Tz (-3.03 kcal mol(-1)) or Bz-Bz (-2.78 kcal mol(-1)). Symmetry-adapted perturbation theory calculations reveal a correlation between the electrostatic component of E(int) and the large increase in the interaction energy for the mixed dimers. However, all components (exchange, induction, dispersion) must be considered to rationalize the observed trend. Another significant conclusion of this work is that basis-set superposition error has a negligible impact on the popular deltaE(CCSD)(T)(MP2) correction, which indicates that counterpoise corrections are not necessary when computing higher-order correlation effects on E(int). Spin-component-scaled MP2 (SCS-MP2 and SCSN-MP2) calculations with a correlation-consistent triple-zeta basis set reproduce the trends in the interaction energies despite overestimating the CCSD(T)/CBS E(int) of Bz-Tz by 20-30%.     

Probing the Watson-Crick, wobble, and sugar-edge hydrogen bond sites of uracil and thymine

The nucleobases uracil (U) and thymine (T) offer three hydrogen-bonding sites for double H-bond formation via neighboring N-H and C=O groups, giving rise to the Watson-Crick, wobble and sugar-edge hydrogen bond isomers. We probe the hydrogen bond properties of all three sites by forming hydrogen bonded dimers of U, 1-methyluracil (1MU), 3-methyluracil (3MU), and T with 2-pyridone (2PY). The mass- and isomer-specific S1 <-- S0 vibronic spectra of 2PY.U, 2PY.3MU, 2PY.1MU, and 2PY.T were measured using UV laser resonant two-photon ionization (R2PI). The spectra of the Watson-Crick and wobble isomers of 2PY.1MU were separated using UV-UV spectral hole-burning. We identify the different isomers by combining three different diagnostic tools: (1) Selective methylation of the uracil N3-H group, which allows formation of the sugar-edge isomer only, and methylation of the N1-H group, which leads to formation of the Watson-Crick and wobble isomers. (2) The experimental S1 <-- S0 origins exhibit large spectral blue shifts relative to the 2PY monomer. Ab initio CIS calculations of the spectral shifts of the different hydrogen-bonded dimers show a linear correlation with experiment. This correlation allows us to identify the R2PI spectra of the weakly populated Watson-Crick and wobble isomers of both 2PY.U and 2PY.T. (3) PW91 density functional calculation of the ground-state binding and dissociation energies De and D0 are in agreement with the assignment of the dominant hydrogen bond isomers of 2PY.U, 2PY.3MU and 2PY.T as the sugar-edge form. For 2PY.U, 2PY.T and 2PY.1MU the measured wobble:Watson-Crick:sugar-edge isomer ratios are in good agreement with the calculated ratios, based on the ab initio dissociation energies and gas-phase statistical mechanics. The Watson-Crick and wobble isomers are thereby determined to be several kcal/mol less strongly bound than the sugar-edge isomers. The 36 observed intermolecular frequencies of the nine different H-bonded isomers give detailed insight into the intermolecular force field.     

Ab initio study of substituent effects in the interactions of dimethyl ether with aromatic rings

Ab initio calculations have been used to investigate the interaction energies of dimers of dimethyl ether with benzene, hexafluorobenzene, and several monosubstituted benzenes. The potential energy curves were explored at the MP2/aug-cc-pVDZ level for two basic configurations of the dimers, one in which the oxygen atom of the dimethyl ether was pointed towards the aromatic ring and the other in which the oxygen atom was pointed away from the aromatic ring. Once the optimum intermolecular distances between the dimethyl and the aromatic ring had been determined for each of the dimers in both configurations at the MP2/aug-cc-pVDZ level, single point energy calculations were performed at the MP2/aug-cc-pVTZ level. A CCSD(T) correction term to the energy was determined and this was combined with the MP2/aug-cc-pVTZ energies to estimate the CCSD(T)/aug-cc-pVTZ interaction energies of the dimers. The estimated CCSD(T)/aug-cc-pVTZ interaction energies are predicted to be attractive for all of the dimers in both configurations and dispersion interactions are found to be a large component of the stabilization of the dimers. For the dimers with the dimethyl ether oxygen pointing towards the aromatic ring, the strengths of interaction energies are found to increase as the aromatic ring becomes more electron deficient, while for the dimers with the dimethyl ether oxygen pointing away from the aromatic ring, they increase as the aromatic ring becomes more electron rich. In both cases, the trends can be explained in terms of the electrostatic potentials of the dimethyl ether and the aromatic rings.     

Interaction of dihydrogen with small and light molecules

Second-order M?ller-Plesset (MP2) calculations (using the approximate resolution of the identity, RI-MP2), explicitly correlated MP2 (MP2-R12) calculations, and coupled-cluster calculations including all single and double excitations with a perturbative estimate of triple excitations [CCSD(T)] are performed to study the interaction of molecular hydrogen with the small molecules HF, H2O, NH3, and LiOH. Different adsorption positions are studied. In the cases of H2O and NH3, the most favorable configuration places H2 in an end-on fashion on the O or N atom, respectively. In the cases of HF and LiOH, the H2 molecule takes a side-on position on the H atom of HF or the Li atom. With respect to MP2 calculations in a triple-zeta basis, both the enlargement of the basis set and the extension of the correlation treatment (CCSD(T) vs MP2) increase the interaction energy. The basis set limit CCSD(T) estimates of the interaction energy of H2 with the HF, H2O, NH3, and LiOH molecules amount to 4.40, 2.67, 3.02, and 10.74 kJ mol-1, respectively. The interaction energy for the simultaneous interaction of H2 with two LiOH molecules does not significantly exceed the value obtained for the interaction with a single LiOH molecule. Furthermore, the interaction energies (by MP2) of H2 with glycine, the glycine dimer, and imidazolium chloride amount to 2.78, 5.00, and 6.30 kJ mol-1, respectively.     

Toward more efficient CCSD(T) calculations of intermolecular interactions in model hydrogen-bonded and stacked dimers

Interaction energies of the model H-bonded complexes, the formamide and formamidine dimers, as well as the stacked formaldehyde and ethylene dimers are calculated by the coupled cluster CCSD(T) method. These systems serve as a model for H-bonded and stacking interactions, typical in molecules participating in biological systems. We use the optimized virtual orbital space (OVOS) technique, by which the dimension of the space of virtual orbitals in coupled cluster CCSD(T) calculations can be significantly reduced. We demonstrate that when the space of virtual orbitals is reduced to 50% of the full space, which means reducing computational demands by 1 order of magnitude, the interaction energies for both H-bonded and stacked dimers are affected by no more than 0.1 kcal/mol. This error is much smaller than the error when interaction energies are calculated using limited basis sets.     

Toward a less costly but accurate calculation of the CCSD(T)/CBS noncovalent interaction energy

The popular method of calculating the noncovalent interaction energies at the coupled-cluster single-, double-, and perturbative triple-excitations [CCSD(T)] theory level in the complete basis set (CBS) limit was to add a CCSD(T) correction term to the CBS second-order Møller-Plesset perturbation theory (MP2). The CCSD(T) correction term is the difference between the CCSD(T) and MP2 interaction energies evaluated in a medium basis set. However, the CCSD(T) calculations with the medium basis sets are still very expensive for systems with more than 30 atoms. Comparatively, the domain-based local pair natural orbital coupled-cluster method [DLPNO-CCSD(T)] can be applied to large systems with over 1,000 atoms. Considering both the computational accuracy and efficiency, in this work, we propose a new scheme to calculate the CCSD(T)/CBS interaction energies. In this scheme, the MP2/CBS term keeps intact and the CCSD(T) correction term is replaced by a DLPNO-CCSD(T) correction term which is the difference between the DLPNO-CCSD(T) and DLPNO-MP2 interaction energies evaluated in a medium basis set. The interaction energies of the noncovalent systems in the S22, HSG, HBC6, NBC10, and S66 databases were recalculated employing this new scheme. The consistent and tight settings of the truncation parameters for DLPNO-CCSD(T) and DLPNO-MP2 in this noncanonical CCSD(T)/CBS calculations lead to the maximum absolute deviation and root-mean-square deviation from the canonical CCSD(T)/CBS interaction energies of less than or equal to 0.28 kcal/mol and 0.09 kcal/mol, respectively. The high accuracy and low cost of this new computational scheme make it an excellent candidate for the study of large noncovalent systems.     

Energetics and mechanism of the decomposition of trifluoromethanol

The thermal instability of alpha-fluoroalcohols is generally attributed to a unimolecular 1,2-elimination of HF, but the barrier to intramolecular HF elimination from CF3OH is predicted to be 45.1 +/- 2 kcal/mol. The thermochemical parameters of trifluoromethanol were calculated using coupled-cluster theory (CCSD(T)) extrapolated to the complete basis set limit. High barriers of 42.9, 43.1, and 45.0 kcal/mol were predicted for the unimolecular decompositions of CH2FOH, CHF2OH, and CF3OH, respectively. These barriers are lowered substantially if cyclic H-bonded dimers of CF3OH with complexation energies of approximately 5 kcal/mol are involved. A six-membered ring dimer has an energy barrier of 28.7 kcal/mol and an eight-membered dimer has an energy barrier of 32.9 kcal/mol. Complexes of CF3OH with HF lead to strong H-bonded dimers, trimers and tetramers with complexation energies of approximately 6, 11, and 16 kcal/mol, respectively. The dimer, CH3OH:HF, and the trimers, CF3OH:2HF and (CH3OH)2:HF, have decomposition energy barriers of 26.7, 20.3, and 22.8 kcal/mol, respectively. The tetramer (CH3OH:HF)2 gives rise to elimination of two HF molecules with a barrier of 32.5 kcal/mol. Either CF3OH or HF can act as catalysts for HF-elimination via an H-transfer relay. Because HF is one of the decomposition products, the decomposition reactions become autocatalytic. If the energies due to complexation for the CF3OH-HF adducts are not dissipated, the effective barriers to HF elimination are lowered from approximately 20 to approximately 9 kcal/mol, which reconciles the computational results with the experimentally observed stabilities.     

Ab initio study of the isomeric equilibrium of the HCN...H2O and H2O...HCN hydrogen-bonded clusters

An ab initio study of the stability, spectroscopic properties, and isomeric equilibrium of the hydrogen-bonded HCN...H2O and H2O...HCN isomers is presented. Density functional theory and perturbative second-order MP2 and coupled-cluster CCSD(T) calculations were carried out and binding energies obtained with correlation-consistent basis sets including extrapolation to the infinity basis set level. At the best theoretical level, CCSD(T), the H2O...HCN complex is more stable than the HCN...H2O complex by ca. 6.3 kJ mol(-1). Rotational and vibrational spectra, including anharmonic corrections, are calculated. These calculated spectroscopic data are used to obtain thermochemical contributions to the thermodynamic functions and hence the Gibbs free energy. The relative free energies are used to estimate the equilibrium constant for isomerism. We find that under typical conditions of supersonic expansion experiments (T < 150 K) H2O...HCN is essentially the only isomer present. Furthermore, our calculations indicate that the hydrogen-bonded cluster becomes favorable over the separated moieties at temperatures below 200 K.     

Stabilisation energy of C(6)H(6)...C(6)X(6) (X = F, Cl, Br, I, CN) complexes: complete basis set limit calculations at MP2 and CCSD(T) levels

Stabilisation energies of stacked structures of C(6)H(6)...C(6)X(6) (X = F, Cl, Br, CN) complexes were determined at the CCSD(T) complete basis set (CBS) limit level. These energies were constructed from MP2/CBS stabilisation energies and a CCSD(T) correction term determined with a medium basis set (6-31G**). The former energies were extrapolated using the two-point formula of Helgaker et al. from aug-cc-pVDZ and aug-cc-pVTZ Hartree-Fock energies and MP2 correlation energies. The CCSD(T) correction term is systematically repulsive. The final CCSD(T)/CBS stabilisation energies are large, considerably larger than previously calculated and increase in the series as follows: hexafluorobenzene (6.3 kcal mol(-1)), hexachlorobenzene (8.8 kcal mol(-1)), hexabromobenzene (8.1 kcal mol(-1)) and hexacyanobenzene (11.0 kcal mol(-1)). MP2/SDD** relativistic calculations performed for all complexes mentioned and also for benzene[dot dot dot]hexaiodobenzene have clearly shown that due to relativistic effects the stabilisation energy of the hexaiodobenzene complex is lower than that of hexabromobenzene complex. The decomposition of the total interaction energy to physically defined energy components was made by using the symmetry adapted perturbation treatment (SAPT). The main stabilisation contribution for all complexes investigated is due to London dispersion energy, with the induction term being smaller. Electrostatic and induction terms which are attractive are compensated by their exchange counterparts. The stacked motif in the complexes studied is very stable and might thus be valuable as a supramolecular synthon.     

pi-Stacked quadruply hydrogen-bonded dimers: pi-stacking influences H-bonding

[reaction: see text] The effects of pi-stacking on the stability of multiply hydrogen-bonded systems are investigated using hybrid DFT calculations on pi-stacked quadruply H-bonded dimers of ureidopyrimidinone in its different tautomeric forms. Both the strengths of the hydrogen bonds and the relative occurrence of tautomers are influenced by pi-stacking; electrostatics and natural bond orbital analysis are used to explain these observations. Finally, these conclusions are independent of the precise nature of the multiply hydrogen-bonded systems, including the DNA base pairs.     

Substituent effects in parallel-displaced pi-pi interactions

High-quality quantum-mechanical methods are used to examine how substituents tune pi-pi interactions between monosubstituted benzene dimers in parallel-displaced geometries. The present study focuses on the effect of the substituent across entire potential energy curves. Substituent effects are examined in terms of the fundamental components of the interaction (electrostatics, exchange-repulsion, dispersion and induction) through the use of symmetry-adapted perturbation theory. Both second-order M?ller-Plesset perturbation theory (MP2) with a truncated aug-cc-pVDZ' basis and spin-component-scaled MP2 (SCS-MP2) with the aug-cc-pVTZ basis are found to mimic closely estimates of coupled-cluster with perturbative triples [CCSD(T)] in an aug-cc-pVTZ basis. Substituents can have a significant effect on the electronic structure of the pi cloud of an aromatic ring, leading to marked changes in the pi-pi interaction. Moreover, there can also be significant direct interactions between a substituent on one ring and the pi-cloud of the other ring.     

An ab initio study of van der Waals potential energy parameters for silver clusters

We employ ab initio calculations of van der Waals complexes to study the potential energy parameters (C(6) coefficients) of van der Waals interactions for modeling of the adsorption of silver clusters on the graphite surface. Electronic structure calculations of the (Ag(2))(2), Ag(2)-H(2), and Ag(2)-C(6)H(6) complexes are performed using a coupled-cluster approach that includes single, double, and perturbative triple excitations (CCSD(T)), M?ller-Plesset second-order perturbation theory (MP2), and spin-component-scaled MP2 (SCS-MP2) methods. Using the atom pair approximation, the C(6) coefficients for silver-silver, silver-hydrogen, and silver-carbon atom systems are obtained after subtracting the energies of quadrupole-quadrupole interactions from the total electronic energy.