Calibration study of the CCSD(T)-F12a/b methods for C2 and small hydrocarbons

Explicitly correlated CCSD(T)-F12a/b methods combined with basis sets specifically designed for this technique have been tested for their ability to reproduce standard CCSD(T) benchmark data covering 16 small molecules composed of hydrogen and carbon. The standard method calibration set was obtained with very large one-particle basis sets, including some aug-cc-pV7Z and aug-cc-pV8Z results. Whenever possible, the molecular properties (atomization energies, structures, and harmonic frequencies) were extrapolated to the complete basis set limit in order to facilitate a direct comparison of the standard and explicitly correlated approaches without ambiguities arising from the use of different basis sets. With basis sets of triple-ζ quality or better, the F12a variant was found to overshoot the presumed basis set limit, while the F12b method converged rapidly and uniformly. Extrapolation of F12b energies to the basis set limit was found to be very effective at reproducing the best standard method atomization energies. Even extrapolations based on the small cc-pVDZ-F12/cc-pVTZ-F12 combination proved capable of a mean absolute deviation of 0.20 kcal/mol. The accuracy and simultaneous cost savings of the F12b approach are such that it should enable high quality property calculations to be performed on chemical systems that are too large for standard CCSD(T).     

Sources of error in electronic structure calculations on small chemical systems

The sources of error in electronic structure calculations arising from the truncation of the one-particle and n-particle expansions are examined with very large correlation consistent basis sets, in some cases up through valence 10-zeta quality, and coupled cluster methods, up through connected quadruple excitations. A limited number of full configuration interaction corrections are also considered. For cases where full configuration interaction calculations were unavailable or prohibitively expensive, a continued fraction approximation was used. In addition, errors arising from corevalence and relativistic corrections are also probed for a number of small chemical systems. The accuracies of several formulas for estimating total energies and atomization energies in the complete basis set limit are compared in light of the present large basis set findings. In agreement with previous work, the CCSD(T) method is found to provide results that are closer to the CCSDTQ and full configuration-interaction results than the less approximate CCSDT method.     

Accurate ab initio binding energies of alkaline earth metal clusters

The effects of basis set superposition error (BSSE) and core-correlation on the electronic binding energies of alkaline earth metal clusters Y(n) (Y = Be, Mg, Ca; n = 2-4) at the Moller-Plesset second-order perturbation theory (MP2) and the single and double coupled cluster method with perturbative triples correction (CCSD(T)) levels are examined using the correlation consistent basis sets cc-pVXZ and cc-pCVXZ (X = D, T, Q, 5). It is found that, while BSSE has a negligible effect for valence-electron-only-correlated calculations for most basis sets, its magnitude becomes more pronounced for all-electron-correlated calculations, including core electrons. By utilizing the negligible effect of BSSE on the binding energies for valence-electron-only-correlated calculations, in combination with the negligible core-correlation effect at the CCSD(T) level, accurate binding energies of these clusters up to pentamers (octamers in the case of the Be clusters) are estimated via the basis set extrapolation of ab initio CCSD(T) correlation energies of the monomer and cluster with only the cc-pVDZ and cc-pVTZ sets, using the basis set and correlation-dependent extrapolation formula recently devised. A comparison between the CCSD(T) and density functional theory (DFT) binding energies is made to identify the most appropriate DFT method for the study of these clusters.     

Accurate calculations of intermolecular interaction energies using explicitly correlated wave functions

Explicitly correlated second-order M?ller-Plesset (MP2-F12) calculations of intermolecular interaction energies for the S22 benchmark set of Jurecka, Sponer, Cerny, and Hobza (Chem. Phys. Phys. Chem. 2006, 8, 1985) are presented and compared with standard MP2 results. The MP2 complete basis set limits are estimated using basis set extrapolation and augmented quadruple-zeta and quintuple-zeta basis sets. Already with augmented double-zeta basis sets the MP2-F12 interaction energies are found to be closer to the complete basis set limits than standard MP2 calculations with augmented quintuple-zeta basis sets. Various possible approximations in the MP2-F12 method are systematically tested. Best results are obtained with localized orbitals and the diagonal MP2-F12/C(D) ansatz. Hybrid approximations, in which some contributions of the auxiliary basis set are neglected and which considerably reduce the computational cost, have a negligible effect on the interaction energies. Also the orbital-invariant fixed-amplitude approximation of Ten-no leads to only slightly less accurate results. Preliminary results for the neon and benzene dimers, obtained with the recently proposed CCSD(T)-F12a approximation, indicate that the CCSD(T) basis set limits can also be very closely approached using augmented triple-zeta basis sets.     

Simple coupled-cluster singles and doubles method with perturbative inclusion of triples and explicitly correlated geminals: The CCSD(T)R12 model

To approach the complete basis set limit of the "gold-standard" coupled-cluster singles and doubles plus perturbative triples [CCSD(T)] method, we extend the recently proposed perturbative explicitly correlated coupled-cluster singles and doubles method, CCSD(2)(R12) [E. F. Valeev, Phys. Chem. Chem. Phys. 8, 106 (2008)], to account for the effect of connected three-electron correlations. The natural choice of the zeroth-order Hamiltonian produces a perturbation expansion with rigorously separable second-order energy corrections due to the explicitly correlated geminals and conventional triple and higher excitations. The resulting CCSD(T)(R12) energy is defined as a sum of the standard CCSD(T) energy and an amplitude-dependent geminal correction. The method is technically very simple: Its implementation requires no modification of the standard CCSD(T) program and the formal cost of the geminal correction is small. We investigate the performance of the open-shell version of the CCSD(T)(R12) method as a possible replacement of the standard complete-basis-set CCSD(T) energies in the high accuracy extrapolated ab initio thermochemistry model of Stanton et al. [J. Chem. Phys. 121, 11599 (2004)]. Correlation contributions to the heat of formation computed with the new method in an aug-cc-pCVXZ basis set have mean absolute basis set errors of 2.8 and 1.0 kJmol when X is T and Q, respectively. The corresponding errors of the standard CCSD(T) method are 9.1, 4.0, and 2.1 kJmol when X=T, Q, and 5. Simple two-point basis set extrapolations of standard CCSD(T) energies perform better than the explicitly correlated method for absolute correlation energies and atomization energies, but no such advantage found when computing heats of formation. A simple Schwenke-type two-point extrapolation of the CCSD(T)(R12)aug-cc-pCVXZ energies with X=T,Q yields the most accurate heats of formation found in this work, in error on average by 0.5 kJmol and at most by 1.7 kJmol.     

Theoretical reference values for the AE6 and BH6 test sets from explicitly correlated coupled-cluster theory

We propose a new computational protocol to obtain highly accurate theoretical reference data. This protocol employs the explicitly correlated coupled-cluster method with iterative single and double excitations as well as perturbative triple excitations, CCSD(T)(F12), using quadruple-z\zeta basis sets. Higher excitations are accounted for by conventional CCSDT(Q) calculations using double-z\zeta basis sets, while core/core-valence correlation effects are estimated by conventional CCSD(T) calculations using quadruple-z\zeta basis sets. Finally, scalar-relativistic effects are accounted for by conventional CCSD(T) calculations using triple-z\zeta basis sets. In the present article, this protocol is applied to the popular test sets AE6 and BH6. An error analysis shows that the new reference values obtained by our computational protocol have an uncertainty of less than 1 kcal/mol (chemical accuracy). Furthermore, concerning the atomization energies, a cancellation of the basis set incompleteness error in the CCSD(T)(F12) perturbative triples contribution with the corresponding error in the contribution from higher excitations is observed. This error cancellation is diminished by the CCSD(T*)(F12) method. Thus, we recommend the use of the CCSD(T*)(F12) method only for small- and medium-sized basis sets, while the CCSD(T)(F12) approach is preferred for high-accuracy calculations in large basis sets.     

Truncation of the correlation consistent basis sets: extension to third-row (Ga-Kr) molecules

The systematic reduction of the commonly used correlation consistent basis sets [cc-pVnZ where n=D(2), T(3), Q(4), and 5] as a means to reduce computational cost has been extended to hydrogen-containing third-row (Ga-Kr) molecules of the G2 test suite. Coupled cluster with singles, doubles, and quasiperturbative triple excitations [CCSD(T)] calculations were performed using both the full correlation consistent basis sets and a series of truncated basis sets in order to assess the impact of basis set reduction upon the structures and energies of the species. The impact that truncation of the basis sets for hydrogen has upon extrapolation of energies to the complete basis set limit also has been examined, and the cost savings that can be achieved are discussed. Overall, basis set reduction can be accomplished which preserves the systematic convergence behavior of the full correlation consistent basis sets.     

A simple and efficient CCSD(T)-F12 approximation

A new explicitly correlated CCSD(T)-F12 approximation is presented and tested for 23 molecules and 15 chemical reactions. The F12 correction strongly improves the basis set convergence of correlation and reaction energies. Errors of the Hartree-Fock contributions are effectively removed by including MP2 single excitations into the auxiliary basis set. Using aug-cc-pVTZ basis sets the CCSD(T)-F12 calculations are more accurate and two orders of magnitude faster than standard CCSD(T)/aug-cc-pV5Z calculations.     

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.     

Ab initio study of hydrazinyl radical: toward a DMBE potential energy surface

A series of stationary structures of the hydrazinyl radical have been characterized by optimization at the CCSD(T)/cc-pVTZ level of theory. CCSD(T)/aug-cc-pVXZ single-point calculations have also been carried out at the optimized geometries with basis sets of different cardinal numbers (X = T, Q), which were used to obtain accurate energies via extrapolation to the complete basis set limit. A discussion on the analytical modeling of the potential energy surface of hydrazinyl is also presented.     

Extrapolation of electron correlation energies to finite and complete basis set targets

The electron correlation energy of two-electron atoms is known to converge asymptotically as approximately (L+1)(-3) to the complete basis set limit, where L is the maximum angular momentum quantum number included in the basis set. Numerical evidence has established a similar asymptotic convergence approximately X(-3) with the cardinal number X of correlation-consistent basis sets cc-pVXZ for coupled cluster singles and doubles (CCSD) and second order perturbation theory (MP2) calculations of molecules. The main focus of this article is to probe for deviations from asymptotic convergence behavior for practical values of X by defining a trial function X(-beta) that for an effective exponent beta=beta(eff)(X,X+1,X+N) provides the correct energy E(X+N), when extrapolating from results for two smaller basis sets, E(X) and E(X+1). This analysis is first applied to "model" expansions available from analytical theory, and then to a large body of finite basis set results (X=D,T,Q,5,6) for 105 molecules containing H, C, N, O, and F, complemented by a smaller set of 14 molecules for which accurate complete basis set limits are available from MP2-R12 and CCSD-R12 calculations. beta(eff) is generally found to vary monotonically with the target of extrapolation, X+N, making results for large but finite basis sets a useful addition to the limited number of cases where complete basis set limits are available. Significant differences in effective convergence behavior are observed between MP2 and CCSD (valence) correlation energies, between hydrogen-rich and hydrogen-free molecules, and, for He, between partial-wave expansions and correlation-consistent basis sets. Deviations from asymptotic convergence behavior tend to get smaller as X increases, but not always monotonically, and are still quite noticeable even for X=5. Finally, correlation contributions to atomization energies (rather than total energies) exhibit a much larger variation of effective convergence behavior, and extrapolations from small basis sets are found to be particularly erratic for molecules containing several electronegative atoms. Observed effects are discussed in the light of results known from analytical theory. A carefully calibrated protocol for extrapolations to the complete basis set limit is presented, based on a single "optimal" exponent beta(opt)(X,X+1,infinity) for the entire set of molecules, and compared to similar approaches reported in the literature.     

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

The accuracy of diffusion quantum Monte Carlo simulations in the determination of molecular equilibrium structures

For a test set of 17 first-row small molecules, the equilibrium structures are calculated with Ornstein-Uhlenbeck diffusion quantum Monte Carlo simulations guiding by trial wave functions constructed from floating spherical Gaussian orbitals and spherical Gaussian geminals. To measure performance of the Monte Carlo calculations, the mean deviation, the mean absolute deviation, the maximum absolute deviation, and the standard deviation of Monte Carlo calculated equilibrium structures with respect to empirical equilibrium structures are given. This approach is found to yield results having a uniformly high quality, being consistent with empirical equilibrium structures and surpassing calculated values from the coupled cluster model with single, double, and noniterative triple excitations [CCSD(T)] with the basis sets of cc-pCVQZ and cc-pVQZ. The nonrelativistic equilibrium atomization energies are also presented to assess performance of the calculated methods. The mean absolute deviations regarding experimental atomization energy are 0.16 and 0.21 kcal/mol for the Monte Carlo and CCSD(T)/cc-pCV(56)Z calculations, respectively.     

Benchmark database of accurate (MP2 and CCSD(T) complete basis set limit) interaction energies of small model complexes, DNA base pairs, and amino acid pairs

MP2 and CCSD(T) complete basis set (CBS) limit interaction energies and geometries for more than 100 DNA base pairs, amino acid pairs and model complexes are for the first time presented together. Extrapolation to the CBS limit is done by using two-point extrapolation methods and different basis sets (aug-cc-pVDZ - aug-cc-pVTZ, aug-cc-pVTZ - aug-cc-pVQZ, cc-pVTZ - cc-pVQZ) are utilized. The CCSD(T) correction term, determined as a difference between CCSD(T) and MP2 interaction energies, is evaluated with smaller basis sets (6-31G** and cc-pVDZ). Two sets of complex geometries were used, optimized or experimental ones. The JSCH-2005 benchmark set, which is now available to the chemical community, can be used for testing lower-level computational methods. For the first screening the smaller training set (S22) containing 22 model complexes can be recommended. In this case larger basis sets were used for extrapolation to the CBS limit and also CCSD(T) and counterpoise-corrected MP2 optimized geometries were sometimes adopted.     

From CCSD(T)/aug-cc-pVTZ-J to CCSD(T) complete basis set limit isotropic nuclear magnetic shieldings via affordable DFT/CBS calculations

It is shown that a linear correlation exists between nuclear shielding constants for nine small inorganic and organic molecules (N(2), CO, CO(2), NH(3), CH(4), C(2)H(2), C(2)H(4), C(2)H(6) and C(6)H(6)) calculated with 47 methods (42 DFT methods, RHF, MP2, SOPPA, SOPPA(CCSD), CCSD(T)) and the aug-cc-pVTZ-J basis set and corresponding complete basis set results, estimated from calculations with the family of polarization-consistent pcS-n basis sets. This implies that the remaining basis set error of the aug-cc-pVTZ-J basis set is very similar in DFT and CCSD(T) calculations. As the aug-cc-pVTZ-J basis set is significantly smaller, CCSD(T)/aug-cc-pVTZ-J calculations allow in combination with affordable DFT/pcS-n complete basis set calculations the prediction of nuclear shieldings at the CCSD(T) level of nearly similar accuracy as those, obtained by fitting results obtained from computationally demanding pcS-n calculations at the CCSD(T) limit. A significant saving of computational efforts can thus be achieved by scaling inexpensive CCSD(T)/aug-cc-pVTZ-J calculations of nuclear isotropic shieldings with affordable DFT complete basis set limit corrections.     

Ab initio calculations on halogen-bonded complexes and comparison with density functional methods

A systematic theoretical investigation on a series of dimeric complexes formed between some halocarbon molecules and electron donors has been carried out by employing both ab initio and density functional methods. Full geometry optimizations are performed at the Moller-Plesset second-order perturbation (MP2) level of theory with the Dunning's correlation-consistent basis set, aug-cc-pVDZ. Binding energies are extrapolated to the complete basis set (CBS) limit by means of two most commonly used extrapolation methods and the aug-cc-pVXZ (X = D, T, Q) basis sets series. The coupled cluster with single, double, and noniterative triple excitations [CCSD(T)] correction term, determined as a difference between CCSD(T) and MP2 binding energies, is estimated with the aug-cc-pVDZ basis set. In general, the inclusion of higher-order electron correlation effects leads to a repulsive correction with respect to those predicted at the MP2 level. The calculations described herein have shown that the CCSD(T) CBS limits yield binding energies with a range of -0.89 to -4.38 kcal/mol for the halogen-bonded complexes under study. The performance of several density functional theory (DFT) methods has been evaluated comparing the results with those obtained from MP2 and CCSD(T). It is shown that PBEKCIS, B97-1, and MPWLYP functionals provide accuracies close to the computationally very expensive ab initio methods.     

Accuracy of basis-set extrapolation schemes for DFT-RPA correlation energies in molecular calculations

We construct a reference benchmark set for atomic and molecular random phase approximation (RPA) correlation energies in a density functional theory framework at the complete basis-set limit. This set is used to evaluate the accuracy of some popular extrapolation schemes for RPA all-electron molecular calculations. The results indicate that for absolute energies, accurate results, clearly outperforming raw data, are achievable with two-point extrapolation schemes based on quintuple- and sextuple-zeta basis sets. Moreover, we show that results in good agreement with the benchmark can also be obtained by using a semiempirical extrapolation procedure based on quadruple- and quintuple-zeta basis sets. Finally, we analyze the performance of different extrapolation schemes for atomization energies.     

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.     

Relativistic effects determined using the Douglas-Kroll contracted basis sets and correlation consistent basis sets with small-core relativistic pseudopotentials

The coupled cluster approximation with single, double, and quasiperturbative triple excitations [CCSD(T)] was used in combination with the Douglas-Kroll contracted correlation consistent basis sets [cc-pVnZ-DK, where n = D(2), T(3), Q(4), and 5] and small-core relativistic pseudopotentials (PP) with correlation consistent polarized valence basis sets (cc-pVnZ-PP and aug-cc-pVnZ-PP) to investigate the impact of scalar relativistic corrections on energetic and structural properties of small molecules containing third-row (Ga-Kr) atoms. These molecules were taken from the Gaussian-2 extended test set for third-row atoms. Atomization energies, ionization energies, electron affinities, and proton affinities for molecules in the test set were determined and compared with nonrelativistic results which were obtained in a recent study in which the standard and augmented correlation consistent basis sets were used in combination with CCSD(T). Several schemes were used to extrapolate the energies to the complete basis set limit.     

Infinite Basis set extrapolation for Double Hybrid Density Functional Theory 1: effect of applying various extrapolation functions

We applied the Infinite Basis (IB) set extrapolation and Double Hybrid Density Functional Theory (DHDF) to calculate the databases of atomization energies, ionization energies, electron affinities, reaction barrier heights, proton affinities, alkyl bond dissociation energies, and noncovalent interactions. The Complete Basis Set (CBS) limit is estimated by extrapolating the hybrid density functional theory and PT2 energies using extrapolation functions including exponential, inverse power, modified exponential, and the combination of the these functions. We found that the combination of B2KPLYP/cc-pV[D|T]Z (which is the extrapolation based on the energies calculated in cc-pVDZ and cc-pVTZ) gives results in quadruple-ζ quality. However, if we want to reach the ～2 kcal/mol chemical accuracy limit, the cc-pV[T|Q]Z is required. Similar results with various extrapolation functions obtained, because the IB parameters were determined by minimizing the averaged mean unsigned error of the calculated databases. We generalized the IB set extrapolation to include more than two basis sets, but we found that extrapolation with two basis sets is satisfactory to give reasonable results. The largest error occurred in the databases of the electron affinities and the weak interactions between the noble gas and the nonpolar molecules. We expect that performing the DHDF-IB scheme with the basis sets augmented by diffuse basis functions will further improve the results.