A complete basis set model chemistry. V. Extensions to six or more heavy atoms

Abstract

The major source of error in most ab initio calculations of molecular energies is the truncation of the one-electron basis set. Extrapolation to the complete basis set second-order (CBS2) limit using the N−1 asymptotic convergence of N-configuration pair natural orbital (PNO) expansions can be combined with the use of relatively small basis sets for the higher-order (i.e., MP3, MP4, and QCI) correlation energy to develop cost effective computational models. Following this strategy, three new computational models denoted CBS-4, CBS-q, and CBS-Q, are introduced. The mean absolute deviations (MAD) from experiment for the 125 energies of the G2 test set are 2.0, 1.7, and 1.0 kcal/mol, respectively. These results compare favorably with the MAD for the more costly G2(MP2), G2, and CBS-QCI/APNO models (1.6, 1.2, and 0.5 kcal/mol, respectively). The error distributions over the G2 test set are indistinguishable from Gaussian distribution functions for all six models, indicating that the rms errors can be interpreted in the same way that experimental uncertainties are used to assess reliability. However, a broader range of examples reveals special difficulties presented by spin contamination, high molecular symmetry, and localization problems in molecules with multiple lone pairs on the same atom. These characteristics can occasionally result in errors several times the size expected from the Gaussian distributions. Each of the CBS models has a range of molecular size for which it is the most accurate computational model currently available. The largest calculations reported for these models include: The CBS-4 heat of formation of tetranitrohydrazine (91.5±5 kcal/mol), the CBS-4 and CBS-q isomerization energies for the conversion of azulene to naphthalene (ΔHcalc=−35.2±1.0 kcal/mol, ΔHexp=−35.3±2.2 kcal/mol), and the CBS-Q heat of formation of SF6 (ΔHcalc=−286.6±1.3 kcal/mol, ΔHexp=−288.3±0.2 kcal/mol). The CBS-Q value for the dissociation energy of a C–H bond in benzene (113.1±1.3 kcal/mol) is also in agreement with the most recent experimental result (112.0±0.6 kcal/mol). The CBS-QCI/APNO model is applicable to the prediction of the C–H bond dissociation energies for the primary (100.7±0.7 kcal/mol calc.) and secondary (97.7±0.7 kcal/mol calc., 97.1±0.4 kcal/mol exp.) hydrogens of propane.

An ignition delay and kinetic modeling study of methane, dimethyl ether, and their mixtures at high pressures

“…We have therefore re-evaluated k(T, p) for this reaction via quantum chemistry and TST/RRKM/ME computations. Geometry optimisation, frequency and singlepoint energy calculations were carried out using the CBS-QB3 [69], CBS-APNO [70] and G3 [71] compound methods. Relaxed potential energy surface scans of dihedral angles were carried out in 10 • increments using the B3LYP method [72,73] with a CBSB7 basis set, and the computed rotational barriers were used in 1-D hindered internal rotation treatments.…”

Kinetic studies of atmospherically relevant silicon chemistry : Part I: Silicon atom reactions

“…Accurate energies (±7 kJ mol -1 ) were then calculated at the CBS-Q level of theory. 39,40 As a benchmark, the heats of formation for SiO, SiO 2 , SiH, SiH 2 ( 1 A 1 ) and SiC calculated at this level of theory (Table 3) who reported different behaviour of this species according to the nature of the buffer gas.…”

In situ generation of HCN for mass spectrometric studies

“…Geometries, absolute energies, vibrational frequencies, and enthalpies for isomers of HC 2 N 2 Ϫ were calculated using density functional theory (B3LYP [20], perturbation theory (MP2 [21,22] and Complete Basis Set methods (CBS [23]). …”