High-level ab initio computations of structures and interaction energies of naphthalene dimers: origin of attraction and its directionality

The intermolecular interaction energies of naphthalene dimers have been calculated by using an aromatic intermolecular interaction model (a model chemistry for the evaluation of intermolecular interactions between aromatic molecules). The CCSD(T) (coupled cluster calculations with single and double substitutions with noniterative triple excitations) interaction energy at the basis set limit has been estimated from the second-order M?ller-Plesset perturbation interaction energy near saturation and the CCSD(T) correction term obtained using a medium-size basis set. The estimated interaction energies of the set of geometries explored in this work show that two structures emerge as being the lowest energy, and may effectively be considered as isoenergetic on the basis of the errors inherent in out extrapolation procedure. These structures are the slipped-parallel (Ci) structure (-5.73 kcal/mol) and the cross (D2d) structure (-5.28 kcal/mol). The T-shaped (C2v) and sandwich (D2h) dimers are substantially less stable (-4.34 and -3.78 kcal/mol, respectively). The dispersion interaction is found to be the major source of attraction in the naphthalene dimer. The electrostatic interaction is substantially smaller than the dispersion interaction. The large dispersion interaction is the cause of the large binding energies of the cross and slipped-parallel dimers.     

Intermolecular interactions of nitrobenzene-benzene complex and nitrobenzene dimer: significant stabilization of slipped-parallel orientation by dispersion interaction

The CCSD(T) level interaction energies of eight orientations of nitrobenzene-benzene complexes and nine orientations of nitrobenzene dimers at the basis set limit have been estimated. The calculated interaction energy of the most stable slipped-parallel (C(s)) nitrobenzene-benzene complex was -4.51 kcal/mol. That of the most stable slipped-parallel (antiparallel) (C(2h)) nitrobenzene dimer was -6.81 kcal/mol. The interaction energies of these complexes are significantly larger than that of the benzene dimer. The T-shaped complexes are substantially less stable. Although nitrobenzene has a polar nitro group, electrostatic interaction is always considerably weaker than the dispersion interaction. The dispersion interaction in these complexes is larger than that in the benzene dimer, which is the cause of the preference of the slipped-parallel orientation in these complexes.     

Intermolecular interaction between hexafluorobenzene and benzene: Ab initio calculations including CCSD(T) level electron correlation correction

The intermolecular interaction energy of hexafluorobenzene-benzene has been calculated with the ARS-E model (a model chemistry for the evaluation of the intermolecular interaction energy between aromatic systems using extrapolation), which was formerly called the AIMI model. The CCSD(T) interaction energy at the basis-set limit has been estimated from the MP2 interaction energy at the basis-set limit and the CCSD(T) correction term obtained using a medium-sized basis set. The slipped-parallel (Cs) complex has the largest (most negative) interaction energy (-5.38 kcal/mol). The sandwich (C6v) complex is slightly less stable (-5.07 kcal/mol). The interaction energies of two T-shaped (C2v) complexes are very small (-1.74 and -0.88 kcal/mol). The calculated interaction energy of the slipped-parallel complex is about twice as large as that of the benzene dimer. The dispersion interaction is found to be the major source of attraction in the complex, although electrostatic interaction also contributes to the attraction. The dispersion interaction increases the relative stability of the slipped-parallel benzene dimer and the hexafluorobenzene-benzene complex compared to T-shaped ones. The electrostatic interaction is repulsive in the slipped-parallel benzene dimer, whereas it stabilizes the slipped-parallel hexafluorobenzene-benzene complex. Both electrostatic and dispersion interactions stabilize the slipped-parallel hexafluorobenzene-benzene complex, which is the cause of the preference of the slipped-parallel orientation and the larger interaction energy of the complex compared to the benzene dimer.     

Origin of attraction and directionality of the pi/pi interaction: model chemistry calculations of benzene dimer interaction.

A model chemistry for the evaluation of intermolecular interaction between aromatic molecules (AIMI Model) has been developed. The CCSD(T) interaction energy at the basis set limit has been estimated from the MP2 interaction energy near the basis set limit and the CCSD(T) correction term obtained by using a medium size basis set. The calculated interaction energies of the parallel, T-shaped,and slipped-parallel benzene dimers are -1.48, -2.46, and -2.48 kcal/mol, respectively. The substantial attractive interaction in benzene dimer, even where the molecules are well separated, shows that the major source of attraction is not short-range interactions such as charge-transfer but long-range interactions such as electrostatic and dispersion. The inclusion of electron correlation increases attraction significantly. The dispersion interaction is found to be the major source of attraction in the benzene dimer. The orientation dependence of the dimer interaction is mainly controlled by long-range interactions. Although electrostatic interaction is considerably weaker than dispersion interaction, it is highly orientation dependent. Dispersion and electrostatic interactions are both important for the directionality of the benzene dimer interaction.     

Model chemistry calculations of thiophene dimer interactions: origin of pi-stacking

The intermolecular interaction energies of thiophene dimers have been calculated by using an aromatic intermolecular interaction (AIMI) model (a model chemistry for the evaluation of intermolecular interactions between aromatic molecules). The CCSD(T) interaction energy at the basis set limit has been estimated from the MP2 interaction energy near the basis set limit and the CCSD(T) correction term obtained by using a medium-size basis set. The calculated interaction energies of the parallel and perpendicular thiophene dimers are -1.71 and -3.12 kcal/mol, respectively. The substantial attractive interaction in the thiophene dimer, even where the molecules are well separated, shows that the major source of attraction is not short-range interactions such as charge transfer but rather long-range interactions such as electrostatic and dispersion. The inclusion of electron correlation increases the attraction significantly. The dispersion interaction is found to be the major source of attraction in the thiophene dimer. The calculated total interaction energy of the thiophene dimer is highly orientation dependent. Although electrostatic interaction is substantially weaker than dispersion interaction, it is highly orientation dependent, and therefore electrostatic interaction play an important role in the orientation dependence of the total interaction energy. The large attractive interaction in the perpendicular dimer is the cause of the preference for the herringbone structure in the crystals of nonsubstituted oligothiophenes (alpha-terthienyls), and the steric repulsion between the beta-substituents is the cause of the pi-stacked structure in the crystals of some beta-substituted oligothiophenes.     

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%.     

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.     

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.     

Reliable structures and energetics for two new delocalized pi...pi prototypes: cyanogen dimer and diacetylene dimer

Two new prototype delocalized pi[dot dot dot]pi complexes are introduced: the dimers of cyanogen, (N[triple bond]C-C[triple bond]N)(2), and diacetylene, (HC[triple bond]C-C[triple bond]CH)(2). These dimers have properties similar to larger delocalized pi...pi systems such as benzene dimer but are small enough that they can be probed in far greater detail with high accuracy electronic structure methods. Parallel-slipped and T-shaped structures of both cyanogen dimer and diacetylene dimer have been optimized with 15 different procedures. The effects of basis set size, theoretical method, counterpoise correction, and the rigid monomer approximation on the structure and energetics of each dimer have been examined. MP2 and CCSD(T) optimized geometries for all four dimer structures are reported, as well as estimates of the CCSD(T) complete basis set (CBS) interaction energy for every optimized geometry. The data reported here suggest that future optimizations of delocalized pi[dot dot dot]pi clusters should be carried out with basis sets of triple-zeta quality. Larger basis sets and the expensive counterpoise correction to the molecular geometry are not necessary. The rigid monomer approximation has very little effect on structure and energetics of these dimers and may be used without consequence. Due to a consistent cancellation of errors, optimization with the MP2 method leads to CCSD(T)/CBS interaction energies that are within 0.2 kcal mol(-1) of those for structures optimized with the CCSD(T) method. Future studies that aim to resolve structures separated by a few tenths of a kcal mol(-1) should consider the effects of optimization with the CCSD(T) method.     

Benzene-pyridine interactions predicted by the effective fragment potential method

The accurate representation of nitrogen-containing heterocycles is essential for modeling biological systems. In this study, the general effective fragment potential (EFP2) method is used to model dimers of benzene and pyridine, complexes for which high-level theoretical data -including large basis spin-component-scaled second-order perturbation theory (SCS-MP2), symmetry-adapted perturbation theory (SAPT), and coupled cluster with singles, doubles, and perturbative triples (CCSD(T))-are available. An extensive comparison of potential energy curves and components of the interaction energy is presented for sandwich, T-shaped, parallel displaced, and hydrogen-bonded structures of these dimers. EFP2 and CCSD(T) potential energy curves for the sandwich, T-shaped, and hydrogen-bonded dimers have an average root-mean-square deviation (RMSD) of 0.49 kcal/mol; EFP2 and SCS-MP2 curves for the parallel displaced dimers have an average RMSD of 0.52 kcal/mol. Additionally, results are presented from an EFP2 Monte Carlo/simulated annealing (MC/SA) computation to sample the potential energy surface of the benzene-pyridine and pyridine dimers.     

Performance of spin-component-scaled Møller-Plesset theory (SCS-MP2) for potential energy curves of noncovalent interactions

An examination of the performance of density-fitted, spin-component-scaled, second-order M?ller-Plesset theory (SCS-MP2), SCS-MP2 with parameters optimized for nucleic acids (SCSN-MP2), and their local-correlation variants, SCS-LMP2 and SCSN-LMP2, is presented for the sandwich and T-shaped benzene dimers, the methane-benzene and H(2)S-benzene complexes, and the methane dimer over entire potential energy curves. These are compared to benchmark-quality estimates of the complete-basis-set limit for coupled-cluster theory through perturbative triple excitations, CCSD(T)/CBS. With the exception of the methane dimer, SCSN-LMP2/CBS tends to outperform SCS-LMP2/CBS with maximum relative errors of 6 and 18%, respectively, at the optimal CCSD(T)/CBS intermolecular distances. For the methane dimer, errors for SCS(N)-(L)MP2/CBS remain in the 0.2-0.3 kcal mol(-1) range, corresponding to a larger relative error of 40-50%. Although the local MP2 methods perform very similarly to their conventional counterparts when aug-cc-pVTZ or larger basis sets are used, in the absence of counterpoise correction the local approximation becomes significantly worse for the aug-cc-pVDZ basis set. The changes due to local correlation approximations for the aug-cc-pVDZ basis are reduced when diffuse functions are neglected for hydrogen atoms.     

True stabilization energies for the optimal planar hydrogen-bonded and stacked structures of guanine...cytosine, adenine...thymine, and their 9- and 1-methyl derivatives: complete basis set calculations at the MP2 and CCSD(T) levels and comparison with experiment

Planar H-bonded and stacked structures of guanine...cytosine (G.C), adenine...thymine (A...T), 9-methylguanine...1-methylcytosine (mG...mC), and 9-methyladenine...1-methylthymine (mA...mT) were optimized at the RI-MP2 level using the TZVPP ([5s3p2d1f/3s2p1d]) basis set. Planar H-bonded structures of G...C, mG...mC, and A...T correspond to the Watson-Crick (WC) arrangement, in contrast to mA...mT for which the Hoogsteen (H) structure is found. Stabilization energies for all structures were determined as the sum of the complete basis set limit of MP2 energies and a (DeltaE(CCSD(T)) - DeltaE(MP2)) correction term evaluated with the cc-pVDZ(0.25,0.15) basis set. The complete basis set limit of MP2 energies was determined by two-point extrapolation using the aug-cc-pVXZ basis sets for X = D and T and X = T and Q. This procedure is required since the convergency of the MP2 interaction energy for the present complexes is rather slow, and it is thus important to include the extrapolation to the complete basis set limit. For the MP2/aug-cc-pVQZ level of theory, stabilization energies for all complexes studied are already very close to the complete basis set limit. The much cheaper D-->T extrapolation provided a complete basis set limit close (by less than 0.7 kcal/mol) to the more accurate T-->Q term, and the D-->T extrapolation can be recommended for evaluation of complete basis set limits of more extended complexes (e.g. larger motifs of DNA). The convergency of the (DeltaE(CCSD(T)) - DeltaE(MP2)) term is known to be faster than that of the MP2 or CCSD(T) correlation energy itself, and the cc-pVDZ(0.25,0.15) basis set provides reasonable values for planar H-bonded as well as stacked structures. Inclusion of the CCSD(T) correction is essential for obtaining reliable relative values for planar H-bonding and stacking interactions; neglecting the CCSD(T) correction results in very considerable errors between 2.5 and 3.4 kcal/mol. Final stabilization energies (kcal/mol) for the base pairs studied are very substantial (A...T WC, 15.4; mA...mT H, 16.3; A...T stacked, 11.6; mA...mT stacked, 13.1; G...C WC, 28.8; mG...mC WC, 28.5; G...C stacked, 16.9; mG...mC stacked, 18.0), much larger than published previously. On the basis of comparison with experimental data, we conclude that our values represent the lower boundary of the true stabilization energies. On the basis of error analysis, we expect the present H-bonding energies to be fairly close to the true values, while stacked energies are still expected to be about 10% too low. The stacking energy for the mG...mC pair is considerably lower than the respective H-bonding energy, but it is larger than the mA...mT H-bonding energy. This conclusion could significantly change the present view on the importance of specific H-bonding interactions and nonspecific stacking interactions in nature, for instance, in DNA. Present stabilization energies for H-bonding and stacking energies represent the most accurate and reliable values and can be considered as new reference data.     

Estimated MP2 and CCSD(T) interaction energies of n-alkane dimers at the basis set limit: comparison of the methods of Helgaker et al. and Feller

The MP2 (the second-order M?ller-Plesset calculation) and CCSD(T) (coupled cluster calculation with single and double substitutions with noniterative triple excitations) interaction energies of all-trans n-alkane dimers were calculated using Dunning's [J. Chem. Phys. 90, 1007 (1989)] correlation consistent basis sets. The estimated MP2 interaction energies of methane, ethane, and propane dimers at the basis set limit [EMP2(limit)] by the method of Helgaker et al. [J. Chem. Phys. 106, 9639 (1997)] from the MP2/aug-cc-pVXZ (X=D and T) level interaction energies are very close to those estimated from the MP2/aug-cc-pVXZ (X=T and Q) level interaction energies. The estimated EMP2(limit) values of n-butane to n-heptane dimers from the MP2/cc-pVXZ (X=D and T) level interaction energies are very close to those from the MP2/aug-cc-pVXZ (X=D and T) ones. The EMP2(limit) values estimated by Feller's [J. Chem. Phys. 96, 6104 (1992)] method from the MP2/cc-pVXZ (X=D, T, and Q) level interaction energies are close to those estimated by the method of Helgaker et al. from the MP2/cc-pVXZ (X=T and Q) ones. The estimated EMP2(limit) values by the method of Helgaker et al. using the aug-cc-pVXZ (X=D and T) are close to these values. The estimated EMP2(limit) of the methane, ethane, propane, n-butane, n-pentane, n-hexane, n-heptane, n-octane, n-nonane, and n-decane dimers by the method of Helgaker et al. are -0.48, -1.35, -2.08, -2.97, -3.92, -4.91, -5.96, -6.68, -7.75, and -8.75 kcal/mol, respectively. Effects of electron correlation beyond MP2 are not large. The estimated CCSD(T) interaction energies of the methane, ethane, propane, and n-butane dimers at the basis set limit by the method of Helgaker et al. (-0.41, -1.22, -1.87, and -2.74 kcal/mol, respectively) from the CCSD(T)/cc-pVXZ (X=D and T) level interaction energies are close to the EMP2(limit) obtained using the same basis sets. The estimated EMP2(limit) values of the ten dimers were fitted to the form m0+m1X (X is 1 for methane, 2 for ethane, etc.). The obtained m0 and m1 (0.595 and -0.926 kcal/mol) show that the interactions between long n-alkane chains are significant. Analysis of basis set effects shows that cc-pVXZ (X=T, Q, or 5), aug-cc-pVXZ (X=D, T, Q, or 5) basis set, or 6-311G** basis set augmented with diffuse polarization function is necessary for quantitative evaluation of the interaction energies between n-alkane chains.     

Differences in structure, energy, and spectrum between neutral, protonated, and deprotonated phenol dimers: comparison of various density functionals with ab initio theory

We have carried out extensive calculations for neutral, cationic protonated, anionic deprotonated phenol dimers. The structures and energetics of this system are determined by the delicate competition between H-bonding, H-π interaction and π-π interaction. Thus, the structures, binding energies and frequencies of the dimers are studied by using a variety of functionals of density functional theory (DFT) and M?ller-Plesset second order perturbation theory (MP2) with medium and extended basis sets. The binding energies are compared with those of highly reliable coupled cluster theory with single, double, and perturbative triple excitations (CCSD(T)) at the complete basis set (CBS) limit. The neutral phenol dimer is unique in the sense that its experimental rotational constants have been measured. The geometry of the neutral phenol dimer is governed by the hydrogen bond formed by two hydroxyl groups and the H-π interaction between two aromatic rings, while the structure of the protonated/deprotonated phenol dimers is additionally governed by the electrostatic and induction effects due to the short strong hydrogen bond (SSHB) and the charges populated in the aromatic rings in the ionic systems. Our salient finding is the substantial differences in structure between neutral, protonated, and deprotonated phenol dimers. This is because the neutral dimer involves in both H(π)···O and H(π)···π interactions, the protonated dimer involves in H(π)···π interactions, and the deprotonated dimer involves in a strong H(π)···O interaction. It is important to compare the reliability of diverse computational approaches employed in quantum chemistry on the basis of the calculational results of this system. MP2 calculations using a small cc-pVDZ basis set give reasonable structures, but those using extended basis sets predict wrong π-stacked structures due to the overestimation of the dispersion energies of the π-π interactions. A few new DFT functionals with the empirical dispersion give reliable results consistent with the CCSD(T)/CBS results. The binding energies of the neutral, cationic protonated, and anionic deprotonated phenol dimers are estimated to be more than 28.5, 118.2, and 118.3 kJ mol(-1), respectively. The energy components of the intermolecular interactions for the neutral, protonated and deprotonated dimers are analyzed.     

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.     

Stacking of the mutagenic DNA base analog 5-bromouracil

The potential energy surface of the stacked 5-bromouracil/uracil (BrU/U) dimer has been investigated in the gas phase and in solution (water and 1,4-dioxane), modeled by a continuum solvent using the polarizable continuum model. Minima and transition states were optimized using DFT (the M06-2X density functional and the 6-31+G(d) basis set). Six stacked gas-phase BrU/U minima were located: four in the face-to-back orientation and two face-to-face. The global minimum in the gas phase is a face-to-face structure with a twist angle of 60° and a zero-point energy-corrected interaction energy of ?10.7 kcal/mol. The BrU/U potential energy surface is geometrically and energetically similar to that of U/U (Hunter and Van Mourik in J Comput Chem 33:2161, 2012). Energy calculations were also performed on experimental geometries of stacked dimers (47 containing BrU stacking with either adenine, cytosine, guanine or thymine and 51 containing thymine also stacking with one of those four bases) taken from DNA structures in the Protein Data Bank. Single-point interaction energies were computed at different levels of theory including MP2, CCSD(T) and DFT using the mPW2PLYP-D double-hybrid functional augmented with an empirical dispersion term, using basis sets ranging from aug-cc-pVDZ to aug-cc-pVQZ. No strong evidence was found for the suggestion that the mutagenicity of BrU is due to enhanced stacking of BrU compared to the corresponding stacked dimers involving thymine.     

Model of peptide bond-aromatic ring interaction: correlated ab initio quantum chemical study

Aromatic ring-peptide bond interactions (modeled as benzene and formamide, N-methylformamide and N-methylacetamide) are studied by means of advanced computational chemistry methods: second-order M?ller-Plesset (MP2), coupled-cluster single and double excitation model [CCSD(T)], and density functional theory with dispersion (DFT-D). The geometrical preferences of these interactions as well as their interaction energy content, in both parallel and T-shaped arrangements, are investigated. The stabilization energy reaches a value of over 5 kcal/mol for the N-methylformamide-benzene complex at the CCSD(T)/complete basis set (CBS) level. Decomposition of interaction energy by the DFT-symmetry-adapted perturbation treatment (SAPT) technique shows that the parallel and T-shaped arrangements, although similar in their total interaction energies, differ significantly in the proportion of electrostatic and dispersion terms.     

Estimates of the ab initio limit for sulfur-pi interactions: the H2S-benzene dimer

The interaction between aromatic rings and sulfur atoms in the side chains of amino acids is a factor in the formation and stabilization of alpha-helices in proteins. We studied the H(2)S-benzene dimer as the simplest possible prototype of sulfur-pi interactions. High-quality potential energy curves were obtained using coupled-cluster theory with single, double, and perturbative triple substitutions (CCSD(T)) and a large, augmented quadruple-zeta basis set (aug-cc-pVQZ). The equilibrium intermonomer distance for the hydrogens-down C(2)(v) configuration is 3.8 A with an interaction energy of -2.74 kcal mol(-1). Extrapolating the binding energy to the complete basis set limit gives -2.81 kcal mol(-1). This binding energy is comparable to that of H(2)O-benzene or of the benzene dimer, and the equilibrium distance is in close agreement with experiment. Other orientations of the dimer were also considered at less complete levels of theory. A considerable reduction in binding for the sulfur-down configuration, together with an energy decomposition analysis, indicates that the attraction in H(2)S-benzene is best thought of as arising from a favorable electrostatic interaction between partially positive hydrogens in H(2)S with the negatively charged pi-cloud of the benzene.     

Origins of the stability of imidazole-imidazole, benzene-imidazole, and benzene-indole dimers: CCSD(T)/CBS and SAPT calculations

The respective structures and stabilities of imidazole-imidazole, benzene-imidazole, and benzene-indole dimers have been investigated using different DFT-D functional, MP2, CCSD(T), and SAPT levels of theory with a medium basis set. Comparative analysis of binding energies and structural parameters of the dimers points to a preference for stacking contact or hydrogen bond in an imidazole-imidazole dimer. In contrast, a T-shaped configuration with H-π interaction is maximally advantageous for benzene-imidazole and benzene-indole dimers. High-level ab initio calculations at the CCSD(T)/CBS and DFT-SAPT levels show that classical hydrogen-bonded tilted imidazole-imidazole dimer is a global minimum structure and that it has high electrostatic energy. However, for benzene-imidazole and benzene-indole dimers, the global minimum (N-H···π) structure has high electrostatic energy as well as dispersion energy.