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.     

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.     

Ab initio calculations of structures and interaction energies of toluene dimers including CCSD(T) level electron correlation correction

The intermolecular interaction energy of the toluene dimer has been calculated with the ARS-F model (a model chemistry for the evaluation of intermolecular interaction energy between ARomatic Systems using Feller's method), which was formerly called as the AIMI model III. 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 Moller-Plesset perturbation interaction energy at the basis set limit obtained by Feller's method and the CCSD(T) correction term obtained using a medium-size basis set. The cross (C(2)) dimer has the largest (most negative) interaction energy (-4.08 kcal/mol). The antiparallel (C(2h)) and parallel (C(S)) dimers (-3.77 and -3.41 kcal/mol, respectively) are slightly less stable. The dispersion interaction is found to be the major source of attraction in the toluene dimer. The dispersion interaction mainly determines the relative stability of the stacked three dimers. The electrostatic interaction of the stacked three dimers is repulsive. Although the T-shaped and slipped-parallel benzene dimers are nearly isoenergetic, the stacked toluene dimers are substantially more stable than the T-shaped toluene dimer (-2.62 kcal/mol). The large dispersion interaction in the stacked toluene dimers is the cause of their enhanced stability.     

Competition between H···π and H···O interactions in furan heterodimers

Here the interactions of furan with HZ (Z = CCH, CCF, CN, Cl, and F) are studied using a variety of electron correlation methods (MP2, CCSD(T), DFT-SAPT) and correlation-consistent triple- and quadruple-ζ basis sets including complete basis set (CBS) extrapolation. For Fu-HF all methods agree that a n-type structure with a hydrogen bridge between the oxygen lone-pair of furan and the hydrogen atom of HF is the global minimum structure. It is found to be significantly more stable than a π-type structure where the hydrogen atom of HF points toward the π system of furan. For the other four dimers MP2 and DFT-SAPT predict the π-type structure to be somewhat more stable, while CCSD(T) favors the n-type structure as the global minimum for Fu-HCl and predicts both structures as nearly isoenergetic for Fu-HCCH and Fu-HCCF. From a geometrical point of view, the Fu-HCN dimer structures are more related to those of the Fu-HCl complex than to Fu-HCCH. The different behavior of HCCF and HF upon complexation with furan evidence the effect of the presence of a π system in the aggregation of fluorine derivatives. It is shown that aggregates of furan cannot be understood by means of dipole-dipole and electrostatic analysis only. Yet, through a combined and detailed analysis of DFT-SAPT energy contributions and resonance effects on the molecular charge distributions a consistent explanation of the aggregation of furan with both π electron rich molecules and halogen hydrides is provided.     

Structural variations and bonding in gold halides: a quantum chemical study of monomeric and dimeric gold monohalide and gold trihalide molecules, AuX, Au2X2, AuX3, and Au2X6 (X = F, Cl, Br, I)

The molecular structures of all gold mono- and trihalides and of their dimers have been calculated at the B3LYP, MP2, and CCSD(T) levels of theory by using relativistic pseudopotentials for all atoms except fluorine. Our computations support the experimental observation that the relative stability of the monohalides increases from the fluoride toward the iodide, while the stability trend of the trihalides is the opposite. The potential energy surface (PES) of all gold trihalides has been investigated. These molecules are typical Jahn-Teller systems; the trigonal planar D3h-symmetry geometry does not correspond to the minimum energy structure for any of them. At the same time, the amount and character of their Jahn-Teller distortion changes gradually from AuF3 to AuI3. The minimum energy geometry is a T-shaped structure for AuF3 and AuCl3, with a Y-shaped transition-state structure. For AuI3, the Y-shaped structure lies lower than the T-shaped structure on the PES. For AuBr3 and AuI3, neither of them is the global minimum but instead an L-shaped structure, which lies outside the Jahn-Teller PES. This structure can be considered to be a donor-acceptor system, or a closed-shell interaction, with I2 acting as donor and AuI as acceptor. The dimers of gold monohalides have very short gold-gold distances and demonstrate the aurophilic interaction. The dimers of the trihalides are planar molecules with two bridging halogen atoms.     

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.     

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.     

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.     

Estimating the basis set superposition error in the CCSD(T)(F12) explicitly correlated method using the example of a water dimer

Changes in the basis set superposition errors upon transitioning from conventional CCSD(T) to the CCSD(T)(F12) explicitly correlated method is studied using the example of a water dimer. A comparison of the compensation errors for CCSD(T) and CCSD(T)(F12) reveals a substantial reduction in the superposition error upon use of the latter. Numerical experiments with water dimers show it is possible theoretically predict an equilibrium distance between oxygen atoms that is similar to the experimental data (2.946 Å), as is the predicted energy of dissociation of a dimer (5.4 ± 0.7 kcal/mol). It is found that the structural and energy parameters of hydrogen bonds in water dimers can be calculated precisely even with two-exponential correlation-consistent basis sets if we use the explicitly correlated approach and subsequently correct the basis set superposition error.     

Halogen bonded complexes between volatile anaesthetics (chloroform, halothane, enflurane, isoflurane) and formaldehyde: a theoretical study

The structures and intermolecular interactions in the halogen bonded complexes of anaesthetics (chloroform, halothane, enflurane and isoflurane) with formaldehyde were studied by ab initio MP2 and CCSD(T) methods. The CCSD(T)/CBS calculated binding energies of these complexes are between -2.83 and -4.21 kcal mol(-1). The largest stabilization energy has been found for the C-Br···O bonded halothane···OCH(2) complex. In all complexes the C-X bond length (where X = Cl, Br) is slightly shortened, in comparison to a free compound, and an increase of the C-X stretching frequency is observed. The electrostatic interaction was excluded as being responsible for the C-X bond contraction. It is suggested that contraction of the C-X bond length can be explained in terms of the Pauli repulsion (the exchange overlap) between the electron pairs of oxygen and halogen atoms in the investigated complexes. This is supported by the DFT-SAPT results, which indicate that the repulsive exchange energy overcompensates the electrostatic one. Moreover, the dispersion and electrostatic contributions cover about 95% of the total attraction forces, in these complexes.     

How accurate is the density functional theory combined with symmetry-adapted perturbation theory approach for CH-pi and pi-pi interactions? A comparison to supermolecular calculations for the acetylene-benzene dimer

Five different orientations of the acetylene-benzene dimer including the T-shaped global minimum structure are used to assess the accuracy of the density functional theory combined with symmetry adapted perturbation theory (DFT-SAPT) approach in its density-fitting implementation (DF-DFT-SAPT) for the study of CH-pi and pi-pi interactions. The results are compared with the outcome of counterpoise corrected supermolecular calculations employing second-order M?ller-Plesset (MP2), spin-component scaled MP2 (SCS-MP2) and single and double excitation coupled cluster theory including perturbative triple excitations (CCSD(T)). For all considered orientations MP2 predicts much deeper potential energy curves with considerably shifted minima compared to CCSD(T) and DFT-SAPT. In spite of being an improvement over the results of MP2, SCS-MP2 tends to underestimate the well depth while DFT-SAPT, employing an asymptotically corrected hybrid exchange-correlation potential in conjunction with the adiabatic local density approximation for the exchange-correlation kernel, is found to be in excellent agreement with CCSD(T). Furthermore, DFT-SAPT provides a detailed understanding of the importance of the electrostatic, induction and dispersion contributions to the total interaction energy and their repulsive exchange corrections.     

On the nature of the stabilization of benzene···dihalogen and benzene···dinitrogen complexes: CCSD(T)/CBS and DFT-SAPT calculations

The structure and stabilization energies of benzene (and methylated benzenes)···X(2) (X=F, Cl, Br, N) complexes were investigated by performing CCSD(T)/complete basis set limit and density functional theory/symmetry-adapted perturbation theory (DFT-SAPT) calculations. The global minimum of the benzene···dihalogen complexes corresponds to the T-shaped structure, whereas that of benzene···dinitrogen corresponds to the sandwich one. The different binding motifs of these complexes arise from the different quadrupole moments of dihalogens and dinitrogen. The different sign of the quadrupole moments of these diatomics is explained based on the electrostatic potential (ESP). Whereas all dihalogens, including difluorine, possess a positive σ hole, such a positive area of the ESP is completely missing in the case of dinitrogen. Moreover, benzene···X(2) (X=Br, Cl) complexes are stronger than benzene···X(2) (X=F, N) complexes. When analyzing DFT-SAPT electrostatic, dispersion, induction, and δ(Hartree-Fock) energies, we recapitulate that the former complexes are stabilized mainly by dispersion energy, followed by electrostatic energy, whereas the latter complexes are stabilized mostly by the dispersion interaction. The charge-transfer energy of benzene···dibromine complexes, and surprisingly, also of methylated benzenes···dibromine complexes is only moderate, and thus, not responsible for their stabilization. Benzene···dichlorine and benzene···dibromine complexes can thus be characterized merely as complexes with a halogen bond rather than as charge-transfer complexes.     

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.     

Investigation of the Benzene–Naphthalene and Naphthalene–Naphthalene Potential Energy Surfaces: DFT/CCSD(T) Correction Scheme

The potential energy surfaces of the naphthalene dimer and benzene–naphthalene complexes are investigated using the recently developed DFT/CCSD(T) correction scheme [J. Chem. Phys. 2008 , 128, 114 102]. One and three minima are located on the PES of the benzene–naphthalene and the naphthalene dimer complexes, respectively, all of which are of the parallel‐displaced type. The stabilities of benzene–naphthalene and the naphthalene dimer are ?4.2 and ?6.2 kcal mol^?1, respectively. Unlike the benzene dimer, where the T‐shaped complex is the global minimum, the lowest‐energy T‐shaped structure is about 0.2 and 1.6 kcal mol^?1 above the global minimum on the benzene–naphthalene and the naphthalene dimer potential energy surfaces, respectively.     

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.     

Characterization of the potential energy surfaces of two small but challenging noncovalent dimers: (P2)2 and (PCCP)2

This work characterizes eight stationary points of the P₂ dimer and six stationary points of the PCCP dimer, including a newly identified minimum on both potential energy surfaces. Full geometry optimizations and corresponding harmonic vibrational frequencies were computed with the second‐order Møller–Plesset (MP2) electronic structure method and six different basis sets: aug‐cc‐pVXZ, aug‐cc‐pV(X+d)Z, and aug‐cc‐pCVXZ where X = T, Q. A new L‐shaped structure with C₂ symmetry is the only minimum for the P₂ dimer at the MP2 level of theory with these basis sets. The previously reported parallel‐slipped structure with C_2h symmetry and a newly identified cross configuration with D₂ symmetry are the only minima for the PCCP dimer. Single point energies were also computed using the canonical MP2 and CCSD(T) methods as well as the explicitly correlated MP2‐F12 and CCSD(T)‐F12 methods and the aug‐cc‐pVXZ (X = D, T, Q, 5) basis sets. The energetics obtained with the explicitly correlated methods were very similar to the canonical results for the larger basis sets. Extrapolations were performed to estimate the complete basis set (CBS) limit MP2 and CCSD(T) binding energies. MP2 and MP2‐F12 significantly overbind the P₂ and PCCP dimers relative to the CCSD(T) and CCSD(T)‐F12 binding energies by as much as 1.5 kcal mol^?1 for the former and 5.0 kcal mol^?1 for the latter at the CBS limit. The dominant attractive component of the interaction energy for each dimer configuration was dispersion according to several symmetry‐adapted perturbation theory analyses. © 2014 Wiley Periodicals, Inc.     

Explicit correlation and basis set superposition error: the structure and energy of carbon dioxide dimer

We have investigated the slipped parallel and t-shaped structures of carbon dioxide dimer [(CO(2))(2)] using both conventional and explicitly correlated coupled cluster methods, inclusive and exclusive of counterpoise (CP) correction. We have determined the geometry of both structures with conventional coupled cluster singles doubles and perturbative triples theory [CCSD(T)] and explicitly correlated cluster singles doubles and perturbative triples theory [CCSD(T)-F12b] at the complete basis set (CBS) limits using custom optimization routines. Consistent with previous investigations, we find that the slipped parallel structure corresponds to the global minimum and is 1.09 kJ mol(-1) lower in energy. For a given cardinal number, the optimized geometries and interaction energies of (CO(2))(2) obtained with the explicitly correlated CCSD(T)-F12b method are closer to the CBS limit than the corresponding conventional CCSD(T) results. Furthermore, the magnitude of basis set superposition error (BSSE) in the CCSD(T)-F12b optimized geometries and interaction energies is appreciably smaller than the magnitude of BSSE in the conventional CCSD(T) results. We decompose the CCSD(T) and CCSD(T)-F12b interaction energies into the constituent HF or HF CABS, CCSD or CCSD-F12b, and (T) contributions. We find that the complementary auxiliary basis set (CABS) singles correction and the F12b approximation significantly reduce the magnitude of BSSE at the HF and CCSD levels of theory, respectively. For a given cardinal number, we find that non-CP corrected, unscaled triples CCSD(T)-F12b/VXZ-F12 interaction energies are in overall best agreement with the CBS limit.     

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.     

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