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.     

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.     

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

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.     

CH/pi interactions in methane clusters with polycyclic aromatic hydrocarbons

Geometries and interaction energies for methane clusters with naphthalene and pyrene were studied. Estimated CCSD(T) interaction energies for the clusters at the basis set limit were -1.92 and -2.50 kcal mol(-1), respectively. Dispersion is mainly responsible for the attraction. Electrostatic interaction is very small. Although the benzene-methane cluster prefers a monodentate structure, in which a C-H bond of the methane points toward the benzene, the methane clusters with the polycyclic aromatic hydrocarbons do not prefer monodentate structures. In the benzene-methane cluster, the weak electrostatic interaction stabilizes the monodentate structure. On the other hand the dispersion interaction controls the orientation of methane in the naphthalene and pyrene clusters. The dispersion interactions in these clusters are significantly larger than those in the benzene-methane cluster. The methane prefers the orientation which is suitable for stabilization by dispersion. Hydrogen atoms of the methane locate above the centers of hexagonal rings of the polycyclic aromatic hydrocarbons in the stable structures. The structures have a small steric repulsion and this positions them only a short distance from the aromatic plane. The large dispersion contribution to the attraction shows that interactions between carbon atoms are mainly responsible for the attraction, and that hydrogen atoms are not important for the attraction. This shows that the interactions between the methane and polycyclic aromatic hydrocarbons are not pi-hydrogen bonds.     

Approximations to complete basis set-extrapolated, highly correlated non-covalent interaction energies

The first-principles calculation of non-covalent (particularly dispersion) interactions between molecules is a considerable challenge. In this work we studied the binding energies for ten small non-covalently bonded dimers with several combinations of correlation methods (MP2, coupled-cluster single double, coupled-cluster single double (triple) (CCSD(T))), correlation-consistent basis sets (aug-cc-pVXZ, X = D, T, Q), two-point complete basis set energy extrapolations, and counterpoise corrections. For this work, complete basis set results were estimated from averaged counterpoise and non-counterpoise-corrected CCSD(T) binding energies obtained from extrapolations with aug-cc-pVQZ and aug-cc-pVTZ basis sets. It is demonstrated that, in almost all cases, binding energies converge more rapidly to the basis set limit by averaging the counterpoise and non-counterpoise corrected values than by using either counterpoise or non-counterpoise methods alone. Examination of the effect of basis set size and electron correlation shows that the triples contribution to the CCSD(T) binding energies is fairly constant with the basis set size, with a slight underestimation with CCSD(T)∕aug-cc-pVDZ compared to the value at the (estimated) complete basis set limit, and that contributions to the binding energies obtained by MP2 generally overestimate the analogous CCSD(T) contributions. Taking these factors together, we conclude that the binding energies for non-covalently bonded systems can be accurately determined using a composite method that combines CCSD(T)∕aug-cc-pVDZ with energy corrections obtained using basis set extrapolated MP2 (utilizing aug-cc-pVQZ and aug-cc-pVTZ basis sets), if all of the components are obtained by averaging the counterpoise and non-counterpoise energies. With such an approach, binding energies for the set of ten dimers are predicted with a mean absolute deviation of 0.02 kcal/mol, a maximum absolute deviation of 0.05 kcal/mol, and a mean percent absolute deviation of only 1.7%, relative to the (estimated) complete basis set CCSD(T) results. Use of this composite approach to an additional set of eight dimers gave binding energies to within 1% of previously published high-level data. It is also shown that binding within parallel and parallel-crossed conformations of naphthalene dimer is predicted by the composite approach to be 9% greater than that previously reported in the literature. The ability of some recently developed dispersion-corrected density-functional theory methods to predict the binding energies of the set of ten small dimers was also examined.     

Origin of attraction, magnitude, and directionality of interactions in benzene complexes with pyridinium cations

Geometries and interaction energies of benzene complexes with pyridine, pyridinium, N-methylpyridinium were studied by ab initio molecular orbital calculations. Estimated CCSD(T) interaction energies of the complexes at the basis set limit were -3.04, -14.77, and -9.36 kcal/mol, respectively. The interactions in the pyridinium and N-methylpyridinium complexes should be categorized into a cation/pi interaction, because the electrostatic and induction interactions greatly contribute to the attraction. On the other hand, the interaction in the pyridine complex is a pi/pi interaction. The dispersion interaction is mainly responsible for the attraction in the benzene-pyridine complex. Short-range interactions including charge-transfer interactions are not important for the attraction in the three complexes. The most stable pyridinium complex has a T-shaped structure, in which the N-H bond points toward the benzene, while the N-methylpyridinium complex prefers a slipped-parallel structure. The benzene-pyridine complex has two nearly isoenergetic (Slipped-parallel and T-shaped) structures.     

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.     

Correction to DFT interaction energies by an empirical dispersion term valid for a range of intermolecular distances

The computation of intermolecular interaction energies via commonly used density functionals is hindered by their inaccurate inclusion of medium and long range dispersion interactions. Practical computation of inter- and intra-macrobiomolecule interaction energies, in particular, requires a fairly accurate yet not overly expensive methodology. It is also desirable to compute intermolecular energies not only at their equilibrium (lowest energy) configurations but also over a range of biophysically relevant distances. We present a method to compute intermolecular interaction energies by including an empirical correction for dispersion which is valid over a range of intermolecular distances. This is achieved by optimizing parameters that moderate the empirical correction by explicit comparison of density functional (B3LYP) energies with distance-dependent (DD) reference values obtained at the CCSD(T)/CBS limit. The resulting method, hereafter referred to as B3LYP-DD, yields interaction energies with an accuracy generally better than 1 kcal mol(-1) for different types of noncovalent complexes, over a range of intermolecular distances and interaction strengths, relative to the expensive CCSD(T)/CBS standard. For a training set of dispersion interacting complexes, B3LYP-DD interaction energies in combination with diffuse functions display absolute errors equal to or smaller than 0.68 kcal mol(-1). The empirical correction does not significantly increase the computational cost as compared to standard density functional calculations. Applications relevant to biomolecular energy and structure, such as prediction of DNA base-pair interactions, are also presented.     

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.     

Amide–π interactions between formamide and benzene

High‐level ab initio calculations have been carried out using a formamide–benzene model system to evaluate amide–π interactions. The interaction energies were estimated as a sum of the CCSD(T) correlation contribution and the HF energy at the complete basis set limit, for the geometries of the model structures at the energy minimum obtained by potential energy surface (PES) scans. NH/π geometry in a face‐on configuration was found to be the most attractive among the various geometries considered, with interaction energy of ?3.75 kcal/mol. An interaction energy of ?2.08 kcal/mol was calculated for the stacked N/Center type geometry, where the nitrogen atom of formamide points directly toward the center of the aromatic ring. The weakest C?O/π geometry, where a carbonyl oxygen atom points toward the plane of the aromatic ring, was found to have energy minimum at an intermolecular distance of 3.67 Å from the PES, with a repulsive interaction energy less than 1 kcal/mol. However, if there are simultaneous attractive interactions with other parts of the molecule besides the amide group, the weak repulsion could be easily overcome, to give a C?O/π geometry interaction. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2009     

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.     

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.     

Theoretical Study of the Benzene Dimer by the Density‐Functional‐Theory Formalism Based on Electron‐Density Partitioning

The density‐functional approach based on the partition into subsystems was applied to study the benzene dimer. For several structures, the calculated interaction energy and intermolecular distance were compared with the previous theoretical results. A good agreement with high level ab initio correlated methods was found. For instance, the interaction energies obtained in this work and the CCSD(T) method agree within 0.1 – 0.6 kcal/mol depending on the structure of the dimer. The structure with the largest interaction energy is T‐shaped, in agreement with CCSD(T) results. The T‐shaped structure of benzene dimer was suggested by several experimental measurements. The calculated interaction energy of 2.09 kcal/mol agrees also well with experimental estimates based on the dissociation energy which ranges from 1.6±0.2 to 2.4±0.4 kcal/mol and the estimated zero‐point vibration energy of 0.3 – 0.5 kcal/mol.