The correlation consistent composite approach (ccCA): an alternative to the Gaussian-n methods

An alternative to the Gaussian-n (G1, G2, and G3) composite methods of computing molecular energies is proposed and is named the "correlation consistent composite approach" (ccCA, ccCA-CBS-1, ccCA-CBS-2). This approach uses the correlation consistent polarized valence (cc-pVXZ) basis sets. The G2-1 test set of 48 enthalpies of formation (DeltaHf), 38 adiabatic ionization potentials (IPs), 25 adiabatic electron affinities (EAs), and 8 adiabatic proton affinities (PAs) are computed using this approach, as well as the DeltaHf values of 30 more systems. Equilibrium molecular geometries and vibrational frequencies are obtained using B3LYP density functional theory. When applying the ccCA-CBS method with the cc-pVXZ series of basis sets augmented with diffuse functions, mean absolute deviations within the G2-1 test set compared to experiment are 1.33 kcal mol(-1) for DeltaHf,0.81 kcal mol(-1) for IPs, 1.02 kcal mol(-1) for EAs, and 1.51 kcal mol(-1) for PAs, without including the "high-level correction" (HLC) contained in the original Gn methods. Whereas the HLC originated in the Gaussian-1 method as an isogyric correction, it evolved into a fitted parameter that minimized the error of the composite methods, eliminating its physical meaning. Recomputing the G1 and G3 enthalpies of formation without the HLC reveals a systematic trend where most DeltaHf values are significantly higher than experimental values. By extrapolating electronic energies to the complete basis set (CBS) limit and adding G3-like corrections for the core-valence and infinite-order electron correlation effects, ccCA-CBS-2 often underestimates the experimental DeltaHf, especially for larger systems. This is desired as inclusion of relativistic and atomic spin-orbit effects subsequently improves theoretical DeltaHf values to give a 0.81 kcal mol(-1) mean absolute deviation with ccCA-CBS-2. The ccCA-CBS method is a viable "black box" method that can be used on systems with at least 10-15 heavy atoms.     

A pseudopotential-based composite method: the relativistic pseudopotential correlation consistent composite approach for molecules containing 4d transition metals (Y-Cd)

The correlation consistent composite approach (ccCA) has proven to be an effective first-principles-based composite approach for main group and first-row transition metal species. By combining relativistic pseudopotentials and ccCA, accurate energetic and thermodynamic data for heavier elements, including transition metals, is obtainable. Relativistic pseudopotential ccCA (rp-ccCA) was formulated and tested on 25 molecules from the G3∕05 set that contain 4p elements (Ga-Kr). A 32.5% time savings was obtained using rp-ccCA, relative to ccCA employing all-electron basis sets. When implementing rp-ccCA to compute dissociation energies and enthalpies of formation for molecules from the 4p block, rp-ccCA results in a mean absolute deviation of 0.89 kcal?mol(-1) from experimental data. rp-ccCA was also applied to a set of 30 4d transition metal-containing molecules, ranging from diatomics to Mo(CO)(6), and enthalpies of formation for these species were obtained with a mean absolute deviation of 2.89 kcal mol(-1) in comparison to experimental data. Based on quality of the experimentally available enthalpies of formation, where the average value of reported experimental error bars is 3.43 kcal mol(-1), rp-ccCA is within transition metal chemical accuracy for the 4d molecule set. rp-ccCA is a pseudopotential-based composite method for transition metals and is shown to yield accurate thermodynamic results for molecules containing heavy elements Ga-Kr and Y-Cd.     

Assessment of Gaussian-3 and density-functional theories on the G3/05 test set of experimental energies

The G3/99 test set [L. A. Curtiss, K. Raghavachari, P. C. Redfern, and J. A. Pople, J. Chem. Phys. 112, 7374 (2000)] of thermochemical data for validation of quantum chemical methods is expanded to include 78 additional energies including 14 enthalpies of formation of the first- and second-row nonhydrogen molecules, 58 energies of molecules containing the third-row elements K, Ca, and Ga-Kr, and 6 hydrogen-bonded complexes. The criterion used for selecting the additional systems is the same as before, i.e., experimental uncertainties less than +/- 1 kcal/mol. This new set, referred to as the G3/05 test set, has a total of 454 energies. The G3 and G3X theories are found to have mean absolute deviations of 1.13 and 1.01 kcal/mol, respectively, when applied to the G3/05 test set. Both methods have larger errors for the nonhydrogen subset of 79 species for which they have mean absolute deviations of 2.10 and 1.64 kcal/mol, respectively. On all of the other types of energies the G3 and G3X methods are very reliable. The G3/05 test set is also used to assess density-functional methods including a series of new functionals. The most accurate functional for the G3/05 test set is B98 with a mean absolute deviation of 3.33 kcal/mol, compared to 4.14 kcal/mol for B3LYP. The latter functional has especially large errors for larger molecules with a mean absolute deviation of 9 kcal/mol for molecules having 28 or more valence electrons. For smaller molecules B3LYP does as well or better than B98 and the other functionals. It is found that many of the density-functional methods have significant errors for the larger molecules in the test set.     

Computational s-block thermochemistry with the correlation consistent composite approach

The correlation consistent composite approach (ccCA) is a model chemistry that has been shown to accurately compute gas-phase enthalpies of formation for alkali and alkaline earth metal oxides and hydroxides (Ho, D. S.; DeYonker, N. J.; Wilson, A. K.; Cundari, T. R. J. Phys. Chem. A 2006, 110, 9767).The ccCA results contrast to more widely used model chemistries where calculated enthalpies of formation for such species can be in error by up to 90 kcal mol-1. In this study, we have applied ccCA to a more general set of 42 s-block molecules and compared the ccCA DeltaHf values to values obtained using the G3 and G3B model chemistries. Included in this training set are water complexes such as Na(H2O)n+ where n = 1 - 4, dimers and trimers of ionic compounds such as (LiCl)2 and (LiCl)3, and the largest ccCA computation to date: Be(acac)2, BeC10H14O4. Problems with the G3 model chemistries seem to be isolated to metal-oxygen bonded systems and Be-containing systems, as G3 and G3B still perform quite well with a 2.7 and 2.6 kcal mol-1 mean absolute deviation (MAD), respectively, for gas-phase enthalpies of formation. The MAD of the ccCA is only 2.2 kcal mol-1 for enthalpies of formation (DeltaHf) for all compounds studied herein. While this MAD is roughly double that found for a ccCA study of >350 main group (i.e., p-block) compounds, it is commensurate with typical experimental uncertainties for s-block complexes. Some molecules where G3/G3B and ccCA computed DeltaHf values deviate significantly from experiment, such as (LiCl)3, NaCN, and MgF, are inviting candidates for new experimental and high-level theoretical studies.     

Hydrogen atom and hydride anion addition to adenine: structures and energetics.

The radicals and anions derived from the 9H tautomer of adenine by adding a hydrogen atom to one of the four double bonds of the adenine framework have been studied. Computations were carried out using a carefully calibrated density functional (B3LYP) method and basis set (DZP++). Optimized geometries, energies, and vibrational frequencies are predicted for eight radicals and anions. The radicals are found to lie in a range of 22 kcal mol(-1), with the radical derived by addition to the C(8) carbon atom being the lowest lying energetically. The anions are predicted to be bound species in the gas phase with an energetic range of 43 kcal mol(-1). Anions produced by addition of a hydride ion to adenine carbon atoms are found to be the most favorable. Six of the anions are predicted to be stable species with respect to electron detachment. The adiabatic electron affinities, vertical electron affinities, and vertical detachment energies are computed for the first time. Electron affinities for these radicals range from 0.0 to 2.0 eV. Radicals produced by addition to a nitrogen atom have near-zero adiabatic electron affinities, while radicals produced by addition at carbon atoms have considerably higher electron affinities.     

The C2H3 + O2 reaction revisited: is multireference treatment of the wave function really critical?

This letter revisits critical intermediates and transition states of the C2H3 + O2 reaction. To obtain their accurate relative energies, ab initio calculations are performed using sophisticated single and multireference theoretical methods with various basis sets. The energy difference between two crucial transition states, for ring opening in dioxiranylmethyl radical and its isomerization to C2H3OO, is calculated as approximately 2 kcal/mol both at multireference MRCI and at single-reference CCSD(T) levels extrapolated to the complete basis set limit. The deviation from the earlier G2M(RCC,MP2) value (approximately 7 kcal/mol) is caused by a deficiency of the 6-311+G(3df,2p) basis set as compared to correlation-consistent Dunning's basis sets.     

Gaussian-4 theory

The Gaussian-4 theory (G4 theory) for the calculation of energies of compounds containing first- (Li-F), second- (Na-Cl), and third-row main group (K, Ca, and Ga-Kr) atoms is presented. This theoretical procedure is the fourth in the Gaussian-n series of quantum chemical methods based on a sequence of single point energy calculations. The G4 theory modifies the Gaussian-3 (G3) theory in five ways. First, an extrapolation procedure is used to obtain the Hartree-Fock limit for inclusion in the total energy calculation. Second, the d-polarization sets are increased to 3d on the first-row atoms and to 4d on the second-row atoms, with reoptimization of the exponents for the latter. Third, the QCISD(T) method is replaced by the CCSD(T) method for the highest level of correlation treatment. Fourth, optimized geometries and zero-point energies are obtained with the B3LYP density functional. Fifth, two new higher level corrections are added to account for deficiencies in the energy calculations. The new method is assessed on the 454 experimental energies in the G305 test set [L. A. Curtiss, P. C. Redfern, and K. Raghavachari, J. Chem. Phys. 123, 124107 (2005)], and the average absolute deviation from experiment shows significant improvement from 1.13 kcal/mol (G3 theory) to 0.83 kcal/mol (G4 theory). The largest improvement is found for 79 nonhydrogen systems (2.10 kcal/mol for G3 versus 1.13 kcal/mol for G4). The contributions of the new features to this improvement are analyzed and the performance on different types of energies is discussed.     

Performance of the correlation consistent composite approach for transition states: a comparison to G3B theory

The correlation consistent composite approach (ccCA) was applied to the prediction of reaction barrier heights (i.e., transition state energy relative to reactants and products) for a standard benchmark set of reactions comprised of both hydrogen transfer reactions and nonhydrogen transfer reactions (i.e., heavy-atom transfer, SN2, and unimolecular reactions). The ccCA method was compared against G3B for the same set of reactions. Error metrics indicate that ccCA achieves "chemical accuracy" with a mean unsigned error (MUE) of 0.89 kcal/mol with respect to the benchmark data for barrier heights; G3B has a mean unsigned error of 1.94 kcal/mol. Further, the greater accuracy of ccCA for predicted reaction barriers is compared to other benchmarked literature methods, including density functional (BB1K, MUE=1.16 kcal/mol) and wavefunction-based [QCISD(T), MUE=1.10 kcal/mol] methods.     

Quantitative computational thermochemistry of transition metal species

The correlation consistent Composite Approach (ccCA), which has been shown to achieve chemical accuracy (+/-1 kcal mol-1) for a large benchmark set of main group and s-block metal compounds, is used to compute enthalpies of formation for a set of 17 3d transition metal species. The training set includes a variety of metals, ligands, and bonding types. Using the correlation consistent basis sets for the 3d transition metals, we find that gas-phase enthalpies of formation can be efficiently calculated for inorganic and organometallic molecules with ccCA. However, until the reliability of gas-phase transition metal thermochemistry is improved, both experimentally and theoretically, a large experimental training set where uncertainties are near +/-1 kcal mol-1 (akin to commonly used main group benchmarking sets) remains an ambitious goal. For now, an average deviation of +/-3 kcal mol-1 appears to be the initial goal of "chemical accuracy" for ab initio transition metal model chemistries. The ccCA is also compared to a more robust but relatively expensive composite approach primarily utilizing large basis set coupled cluster computations. For a smaller training set of eight molecules, ccCA has a mean absolute deviation (MAD) of 3.4 kcal mol-1 versus the large basis set coupled-cluster-based model chemistry, which has a MAD of 3.1 kcal mol-1. However, the agreement for transition metal complexes is more system dependent than observed in previous benchmark studies of composite methods and main group compounds.     

Structures and heats of formation of the neutral and ionic PNO, NOP, and NPO systems from electronic structure calculations

High level ab initio electronic structure calculations using the coupled cluster CCSD(T) method with augmented correlation-consistent basis sets extrapolated to the complete basis set limit have been performed on the PNO, NOP, and NPO isomers and their corresponding anions and cations. Geometries for all species were optimized up through the aug-cc-pV(Q+d)Z level and vibrational frequencies were calculated with the aug-cc-pV(T+d)Z basis set. The most stable of the three isomers is NPO and it is predicted to have a heat of formation of 23.3 kcal/mol. PNO is predicted to be only 1.7 kcal/mol higher in energy. The calculated adiabatic ionization potential of NPO is 12.07 eV and the calculated adiabatic electron affinity is 2.34 eV. The calculated adiabatic ionization potential of PNO is 10.27 eV and the calculated adiabatic electron affinity is only 0.24 eV. NOP is predicted to be much higher in energy by 29.9 kcal/mol. The calculated rotational constants for PNO and NPO should allow for these species to be spectroscopically distinguished. The adiabatic bond dissociation energies for the P[Single Bond]N, P[Single Bond]O, and N[Single Bond]O bonds in NPO and PNO are the same within approximately 10 kcal/mol and fall in the range of 72-83 kcal/mol.     

Gaussian-3 and related methods for accurate thermochemistry

In this paper, we present an overview of Gaussian-3 (G3) theory, a composite technique that employs a sequence of ab initio molecular orbital calculations to derive a total energy of a given molecular species. This method provides accurate energies of molecular systems for the calculation of enthalpies of formation, ionization potentials, electron affinities, proton affinities, etc. Also covered in this review are several variants of G3 theory including one based on scale factors (G3S) and an extended version (G3X) that uses improved geometries and larger Hartree-Fock basis sets. Finally, the G3/99 test set of accurate experimental data that is used for critical assessment is described. Overall, G3 theory has a mean absolute deviation from experiment of 1.07 kcal mol⁻¹ for the G3/99 test set and G3S theory has a similar accuracy of 1.08 kcal mol⁻¹. G3X theory is significantly more accurate with the mean absolute deviation from experiment decreasing from 1.07 kcal mol⁻¹ (G3) to 0.95 kcal mol⁻¹ (G3X). The scaled version of G3X theory shows a similar improvement. Received: 23 January 2002 / Accepted: 7 April 2002 / Published online: 4 July 2002     

Accurate energetics of small molecules containing third-row atoms Ga-Kr: a comparison of advanced ab initio and density functional theory

Advanced ab initio [coupled cluster theory through quasiperturbative triple excitations (CCSD(T))] and density functional (B3LYP) computational chemistry approaches were used in combination with the standard and augmented correlation consistent polarized valence basis sets [cc-pVnZ and aug-cc-pVnZ, where n=D(2), T(3), Q(4), and 5] to investigate the energetic and structural properties of small molecules containing third-row (Ga-Kr) atoms. These molecules were taken from the Gaussian-2 (G2) extended test set for third-row atoms. Several different schemes were used to extrapolate the calculated energies to the complete basis set (CBS) limit for CCSD(T) and the Kohn-Sham (KS) limit for B3LYP. Zero point energy and spin orbital corrections were included in the results. Overall, CCSD(T) atomization energies, ionization energies, proton affinities, and electron affinities are in good agreement with experiment, within 1.1 kcal/mol when the CBS limit has been determined using a series of two basis sets of at least triple zeta quality. For B3LYP, the overall mean absolute deviation from experiment for the three properties and the series of molecules is more significant at the KS limit, within 2.3 and 2.6 kcal/mol for the cc-pVnZ and aug-cc-pVnZ basis set series, respectively.     

Implementation of pseudopotential in the G3 theory for molecules containing first-, second-, and non-transition third-row atoms

Compact effective pseudopotential (CEP) is adapted in the G3 theory providing a theoretical alternative referred to as G3CEP for calculations involving the first-, second-, and non-transition third-row elements. These modifications tried to preserve as much as possible the original characteristics of G3. G3CEP was used in the study of 247 enthalpies of formation, 22 atomization energies, 104 ionization potentials, 63 electron affinities, and 10 proton affinities, resulting in the calculation of 446 species for the first-, second-, and third-row atoms. The final average total absolute deviation was of 1.29 kcal mol(-1) against 1.16 kcal mol(-1) from all-electron G3 for the same calculations. The CPU time has been reduced by 7% to 56%, depending on the size of the molecules and the type of atoms considered.     

The 6-31B(d) basis set and the BMC-QCISD and BMC-CCSD multicoefficient correlation methods

Three new multicoefficient correlation methods (MCCMs) called BMC-QCISD, BMC-CCSD, and BMC-CCSD-C are optimized against 274 data that include atomization energies, electron affinities, ionization potentials, and reaction barrier heights. A new basis set called 6-31B(d) is developed and used as part of the new methods. BMC-QCISD has mean unsigned errors in calculating atomization energies per bond and barrier heights of 0.49 and 0.80 kcal/mol, respectively. BMC-CCSD has mean unsigned errors of 0.42 and 0.71 kcal/mol for the same two quantities. BMC-CCSD-C is an equally effective variant of BMC-CCSD that employs Cartesian rather than spherical harmonic basis sets. The mean unsigned error of BMC-CCSD or BMC-CCSD-C for atomization energies, barrier heights, ionization potentials, and electron affinities is 22% lower than G3SX(MP2) at an order of magnitude less cost for gradients for molecules with 9-13 atoms, and it scales better (N6 vs N,7 where N is the number of atoms) when the size of the molecule is increased.     

Proton-bound homodimers involving second-row atoms

High-level ab initio quantum chemical calculations (G4(MP2)//MP2/6-311+G(2df,p)) have been used to examine homodimers of second-row bases, and to compare the results with those obtained previously for the first-row analogs. The relationship between the binding energies of the dimers and the proton affinities (PAs) of the bases follows the same pattern as that for the first-row systems, with the binding energies initially increasing with increasing proton affinity but subsequently decreasing. This may be attributed to the opposing effects of increased PA on the hydrogen-bond donor and hydrogen-bond acceptor. The binding energies are generally smaller for the second-row dimers than for the corresponding first-row dimers. There is an increased tendency for asymmetrical hydrogen bonds in homodimers of the second-row compared with first-row dimers. This may be attributed to the lower electronegativities of second-row atoms relative to their first-row counterparts, and to the longer internuclear separation between the hydrogen-bonded second-row atoms.     

Approaching chemical accuracy with quantum Monte Carlo

A quantum Monte Carlo study of the atomization energies for the G2 set of molecules is presented. Basis size dependence of diffusion Monte Carlo atomization energies is studied with a single determinant Slater-Jastrow trial wavefunction formed from Hartree-Fock orbitals. With the largest basis set, the mean absolute deviation from experimental atomization energies for the G2 set is 3.0 kcal/mol. Optimizing the orbitals within variational Monte Carlo improves the agreement between diffusion Monte Carlo and experiment, reducing the mean absolute deviation to 2.1 kcal/mol. Moving beyond a single determinant Slater-Jastrow trial wavefunction, diffusion Monte Carlo with a small complete active space Slater-Jastrow trial wavefunction results in near chemical accuracy. In this case, the mean absolute deviation from experimental atomization energies is 1.2 kcal/mol. It is shown from calculations on systems containing phosphorus that the accuracy can be further improved by employing a larger active space.     

Thermochemical properties of selenium fluorides, oxides, and oxofluorides

The bond dissociation energies (BDEs), fluoride and fluorocation affinities, and electron affinities of SeF(n) (n = 1-6), SeOF(n) (n = 0-4), and SeO(2)F(n) (n = 0-2) have been predicted with coupled cluster CCSD(T) theory extrapolated to the complete basis set limit. To achieve near chemical accuracy, additional corrections were added to the complete basis set binding energies based on frozen core coupled cluster theory energies. These included corrections for core-valence effects, scalar relativistic effects, for first-order atomic spin-orbit effects, and vibrational zero point energies. The adiabatic BDEs contain contributions from product reorganization energies and, therefore, can be much smaller than the diabatic BDEs and can vary over a wide range. For thermochemical calculations, the adiabatic values must be used, whereas for bond strength and kinetic considerations, the diabatic values should be used when only small displacements of the atoms without change of the geometry of the molecule are involved. The adiabatic Se-F BDEs of SeF(n) (n = 1-6) are SeF(6) = 90, SeF(5) = 27, SeF(4) = 93, SeF(3) = 61, SeF(2) = 86, and SeF = 76 kcal/mol, and the corresponding diabatic values are SeF(6) = 90, SeF(5) = 88, SeF(4) = 93, SeF(3) = 74, SeF(2) = 86, and SeF = 76 kcal/mol. The adiabatic Se-O BDEs of SeO(n) (n = 1-3), SeOF(n) (n = 1-4), and SeO(2)F(n) (n = 1,2) range from 23 to 107 kcal/mol, whereas the diabatic ones range from 62 to 154 kcal/mol. The adiabatic Se-F BDEs of SeOF(n) (n = 1-4) and SeO(2)F(n) (n = 1,2) range from 20 to 88 kcal/mol, whereas the diabatic ones range from 73 to 112 kcal/mol. The fluoride affinities of SeF(n), (n = 1-6), SeO(n), (n = 1-3), SeOF(n), (n = 1-4), and SeO(2)F(n) (n = 1,2) range from 15 to 121 kcal/mol, demonstrating that the Lewis acidity of these species covers the spectrum from very weak (SeF(6)) to very strong (SeO(3)) acids. The electron affinities which are a measure of the oxidizing power of a species, span a wide range from 1.56 eV in SeF(4) to 5.16 eV in SeF(5) and for the free radicals are much higher than for the neutral molecules. Another interesting feature of these molecules and ions stems from the fact that many of them possess both a Se free valence electron pair and a free unpaired valence electron, raising the questions of their preferred location and their influence on the Se-F and Se═O bond strengths.     

Theoretical study of the spectroscopically relevant parameters for the detection of HNP(q) and HPN(q) (q = 0, +1, -1) in the gas phase

High level ab initio electronic structure calculations at different levels of theory have been performed on HNP and HPN neutrals, anions, and cations. This includes standard coupled cluster CCSD(T) level with augmented correlation-consistent basis sets, internally contacted multi-reference configuration interaction, and the newly developed CCSD(T)-F12 methods in connection with the explicitly correlated basis sets. Core-valence correction and scalar relativistic effects were examined. We present optimized equilibrium geometries, harmonic vibrational frequencies, rotational constants, adiabatic ionization energies, electron affinities, vertical detachment energies, and relative energies. In addition, the three-dimensional potential energy surfaces of HNP(-1,0,+1) and of HPN(-1,0,+1) were generated at the (R)CCSD(T)-F12b∕cc-pVTZ-F12 level. The anharmonic terms and fundamentals were derived using second order perturbation theory. For HNP, our best estimate for the adiabatic ionization energy is 7.31 eV, for the adiabatic electron affinity is 0.47 eV. The higher energy isomer, HPN, is 23.23 kcal∕mol above HNP. HPN possesses a rather large adiabatic electron affinity of 1.62 eV. The intramolecular isomerization pathways were computed. Our calculations show that HNP(-) to HPN(-) reaction is subject to electron detachment.     

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