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.     

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.     

The correlation-consistent composite approach: application to the G3/99 test set

The correlation-consistent composite approach (ccCA), an ab initio composite technique for computing atomic and molecular energies, recently has been shown to successfully reproduce experimental data for a number of systems. The ccCA is applied to the G3/99 test set, which includes 223 enthalpies of formation, 88 adiabatic ionization potentials, 58 adiabatic electron affinities, and 8 adiabatic proton affinities. Improvements on the original ccCA formalism include replacing the small basis set quadratic configuration interaction computation with a coupled cluster computation, employing a correction for scalar relativistic effects, utilizing the tight-d forms of the second-row correlation-consistent basis sets, and revisiting the basis set chosen for geometry optimization. With two types of complete basis set extrapolation of MP2 energies, ccCA results in an almost zero mean deviation for the G3/99 set (with a best value of -0.10 kcal mol(-1)), and a 0.96 kcal mol(-1) mean absolute deviation, which is equivalent to the accuracy of the G3X model chemistry. There are no optimized or empirical parameters included in the computation of ccCA energies. Except for a few systems to be discussed, ccCA performs as well as or better than Gn methods for most systems containing first-row atoms, while for systems containing second-row atoms, ccCA is an improvement over Gn model chemistries.     

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.     

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.     

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.     

Prediction of hydrocarbon enthalpies of formation by various thermochemical schemes

The correlation consistent composite approach (ccCA) has been used to compute the enthalpies of formation (ΔH_f′s) for 60 closed‐shell, neutral hydrocarbon molecules selected from an established set (Cioslowski et al., J. Chem. Phys. 2000 , 113, 9377). This set of thermodynamic values includes ΔH_f's for hydrocarbons that span a range of molecular sizes, degrees of aromaticity, and geometrical configurations, and, as such, provides a rigorous assessment of ccCA's applicability to a variety of hydrocarbons. The ΔH_f's were calculated via atomization energies, isodesmic reactions, and hypohomodesmotic reactions. In addition, for 12 of the aromatic molecules in the set that are larger than benzene, the energies of ring‐conserved isodesmic reactions were used to calculate the ΔH_f′s. Using an atomization energy approach to determine the ΔH_f′s, the lowest mean absolute deviation (MAD) from experiment achieved by ccCA for the 60 hydrocarbons was 1.10 kcal mol^?1. The use of the mixed Gaussian/inverse exponential complete basis set extrapolation scheme (ccCA‐P) in conjunction with hypohomodesmotic reaction energies resulted in a MAD of 0.87 kcal mol^?1. This value is compared with MADs of 1.17, 1.18, and 1.28 kcal mol^?1 obtained via the Gaussian‐4 (G4), Gaussian‐3 (G3), and Gaussian‐3(MP2) [G3(MP2)] methods, respectively (using the hypohomodesmotic reactions). © 2012 Wiley Periodicals, Inc.     

Domain-based local pair natural orbital methods within the correlation consistent composite approach

Ab initio composite approaches have been utilized to model and predict main group thermochemistry within 1 kcal mol⁻¹, on average, from well-established reliable experiments, primarily for molecules with less than 30 atoms. For molecules of increasing size and complexity, such as biomolecular complexes, composite methodologies have been limited in their application. Therefore, the domain-based local pair natural orbital (DLPNO) methods have been implemented within the correlation consistent composite approach (ccCA) framework, namely DLPNO-ccCA, to reduce the computational cost (disk space, CPU (central processing unit) time, memory) and predict energetic properties such as enthalpies of formation, noncovalent interactions, and conformation energies for organic biomolecular complexes including one of the largest molecules examined via composite strategies, within 1 kcal mol⁻¹, after calibration with 119 molecules and a set of linear alkanes. © 2019 Wiley Periodicals, Inc.     

Toward accurate theoretical thermochemistry of first row transition metal complexes

The recently developed correlation consistent Composite Approach for transition metals (ccCA-TM) was utilized to compute the thermochemical properties for a collection of 225 inorganic molecules containing first row (3d) transition metals, ranging from the monohydrides to larger organometallics such as Sc(C(5)H(5))(3) and clusters such as (CrO(3))(3). Ostentatiously large deviations of ccCA-TM predictions stem mainly from aging and unreliable experimental data. For a subset of 70 molecules with reported experimental uncertainties less than or equal to 2.0 kcal mol(-1), regardless of the presence of moderate multireference character in some molecules, ccCA-TM achieves transition metal chemical accuracy of ±3.0 kcal mol(-1) as defined in our earlier work [J. Phys. Chem. A2007, 111, 11269-11277] by giving a mean absolute deviation of 2.90 kcal mol(-1) and a root-mean-square deviation of 3.91 kcal mol(-1). As subsets are constructed with decreasing upper limits of reported experimental uncertainties (5.0, 4.0, 3.0, 2.0, and 1.0 kcal mol(-1)), the ccCA-TM mean absolute deviations were observed to monotonically drop off from 4.35 to 2.37 kcal mol(-1). In contrast, such a trend is missing for DFT methods as exemplified by B3LYP and M06 with mean absolute deviations in the range 12.9-14.1 and 10.5-11.0 kcal mol(-1), respectively. Salient multireference character, as demonstrated by the T(1)/D(1) diagnostics and the weights (C(0)(2)) of leading electron configuration in the complete active self-consistent field wave function, was found in a significant amount of molecules, which can still be accurately described by the single reference ccCA-TM. The ccCA-TM algorithm has been demonstrated as an accurate, robust, and widely applicable model chemistry for 3d transition metal-containing species with versatile bonding features.     

Polyfunctional methodology for improved DFT thermochemical predictions

Statistical error distributions for enthalpies of formation as predicted by 18 different density functionals have been analyzed using a test set of 675 molecules. Systematic errors, dependent on the number of valence electrons, have been identified for some functionals. A simple empirical correction makes a significant improvement in the prediction error for these single functionals. Linear combinations of enthalpy estimates from different density functionals are identified that exploit the error correlations among the functionals and allow for further improvements in the accuracy of thermodynamic predictions. A good compromise between accuracy and computational efforts is achieved by the BLUE (best linear unbiased estimator) combination of three functionals, B3LYP, BLYP, and VSXC (polyfunctional 3 or PF3). The PF3 method has a mean absolute deviation (MAD) from experiment of 2.4 kcal/mol on the G3 set of 271 molecules. This can be compared to the MAD of 4.9 kcal/mol for B3LYP and 1.2 kcal/mol for the more costly G3 method. On the larger set of 675 molecules, the MAD for PF3 is 3.0 kcal/mol. Opportunities for further improvements in the accuracy of this method are discussed.     

Accurate thermochemistry from corrected Hartree–Fock results: rapid estimation of nearly experimental quality total energy using the small 6-31G(d) basis set

Gaussian-3 ground-state total electronic energies have been approximated using single point 6-31G(d) basis set Harteee–Fock self-consistent-field (HF-SCF) total energies and partial charges based on our earlier rapid estimation of correlation energy from partial charges method. Sixty-five closed-shell neutral molecules (composed of H, C, N, O, and F atoms) of the G2/97 thermochemistry database were selected for the present study. The main feature in this work is that the␣basis set error has been treated by the least squares fit of rapid estimation of basis set error and correlation energy from partial charges (REBECEP) atomic parameters. With these parameters a rather accurate closed-shell ground-state electronic total energy can be obtained from a small basis set HF-SCF calculation in the vicinity of stationary points. The average absolute deviation of the best REBECEP enthalpies of formation from the experimental enthalpies of formation is 1.39 kcal/mol for the test set of 65 enthalpies of neutral molecules. Received: 11 December 2000 / Accepted: 6 February 2001/Published online: 11 October 2001     

Improved modification for the density-functional theory calculation of thermodynamic properties for C-H-O composite compounds

A three-parametric modification equation and the least-squares approach are adopted to calibrating hybrid density-functional theory energies of C(1)-C(10) straight-chain aldehydes, alcohols, and alkoxides to accurate enthalpies of formation DeltaH(f) and Gibbs free energies of formation DeltaG(f), respectively. All calculated energies of the C-H-O composite compounds were obtained based on B3LYP6-311++G(3df,2pd) single-point energies and the related thermal corrections of B3LYP6-31G(d,p) optimized geometries. This investigation revealed that all compounds had 0.05% average absolute relative error (ARE) for the atomization energies, with mean value of absolute error (MAE) of just 2.1 kJ/mol (0.5 kcal/mol) for the DeltaH(f) and 2.4 kJ/mol (0.6 kcal/mol) for the DeltaG(f) of formation.     

Enthalpy of the gas-phase CO2 + Mg reaction from ab initio total energies

Various highly accurate ab initio composite methods of Gaussian-n (G1, G2, G3), their variations (G2(MP2), G3(MP2), G3//B3LYP, G3(MP2)//B3LYP), and complete basis set (CBS-Q, CBS-Q//B3LYP) series of models were applied to compute reaction enthalpies of the ground-state reaction of CO2 with Mg. All model chemistries predict highly endothermic reactions, with DeltaH(298) = 63.6-69.7 kcal x mol(-1). The difference between the calculated reaction enthalpies and the experimental value, evaluated with recommended experimental standard enthalpies of formation for products and reactants, is more than 20 kcal x mol(-1) for all methods. This difference originates in the incorrect experimental enthalpy of formation of gaseous MgO given in thermochemical databases. When the theoretical formation enthalpy for MgO calculated by a particular method is used, the deviation is reduced to 1.3 kcal x mol(-1). The performance of the methodologies used to calculate the heat of this particular reaction and the enthalpy of formation of MgO are discussed.     

Assessment of Gaussian-4 theory for the computation of enthalpies of formation of large organic molecules

Gas-phase enthalpies of formation of 122 relatively large organic molecules with up to 15 non-hydrogen atoms have been calculated at the Gaussian-4 (G4) level of theory using the atomization reaction procedure. The calculated values were compared with experimental data published mainly last years. Particular attention has been given to nitro compounds and nitrogen, oxygen, and sulfur containing heterocyclic compounds. The expected accumulation of systematic errors as the molecular size increases was not observed with increasing the number of non-hydrogen atoms from 6 to 15. The largest mean absolute deviation between experimental and G4 enthalpies of formation, 10.7 kJ/mol, was revealed for nitro compounds. All theoretical values for nitro compounds were underestimated by 5–15 kJ/mol. The best agreement with experiment with mean absolute deviation of 4.5 kJ/mol was observed for compounds which types were widely presented in the original test set of G4 method. The mean absolute deviations for nitrogen heterocycles (6.8 kJ/mol) and oxygen and sulfur heterocycles (9.1 kJ/mol) are noticeably larger. Experimental enthalpies of formation of four compounds (N,N-dinitromethanamine, 2,3,5,6-tetrachloronitrobenzene, 2-methyl-2H-tetrazole, and proline) were suggested to be unreliable from comparison with the G4 values calculated from atomization energies and isodesmic reactions.     

Investigation of the use of B3LYP zero-point energies and geometries in the calculation of enthalpies of formation

The use of B3LYP/6–31G^* zero-point energies and geometries in the calculation of enthalpies of formation has been investigated for the enlarged G2 test set of 148 molecules [J. Chem. Phys. 106 (1997) 1063]. A scale factor of 0.96 for the B3LYP zero-point energies gives an average absolute deviation nearly the same as scaled HF/6–31G^* zero-point energies for G2, G2(MP2), and B3LYP/6–311 + G(3df,2p) enthalpies. A scale factor of 0.98, which has been recommended in some studies, increases the average absolute deviation by about 0.2 kcal/mol. Geometries from B3LYP/6–31G^* are found to do as well as MP2/6–31G^* geometries in the calculation of the enthalpies of formation.     

Theoretical thermochemistry: Enthalpies of formation of a set of nitrogen‐containing compounds

Gas‐phase standard state formation enthalpies of 63 nitrogen‐containing compounds (NCCs) were studied by computational methods. Gaussian‐n and complete basis set composite methods were applied. After the calculation of the set of NCCs, the results were analyzed in several ways. All the seven selected methods depicted various precisions in this work. According to the calculations and further data processing, G4 was proven to be appropriate and capable of formation enthalpy calculations on NCCs with a mean absolute deviation of 0.63 kcal/mol. Thus, G4 calculation may help us on choosing experimental values and conducting predictions. Isodesmic reactions using G4 methods are also conducted and further correction was given. The high accuracy of G4 method makes it reliable for calculating gas‐phase enthalpies of formations and allows them to serve as a valuable check on the accuracy of reported experimental data. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

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.     

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.     

Quantum monte carlo study of heats of formation and bond dissociation energies of small hydrocarbons

A quantum Monte Carlo (QMC) benchmark study of heats of formation at 298 K and bond dissociation energies (BDEs) of 22 small hydrocarbons is reported. Diffusion Monte Carlo (DMC) results, obtained using a simple product trial wavefunctions consisting of a single determinant and correlation function, are compared to experiment and to other theory including a version of complete basis set theory (CBS‐Q) and density functional theory (DFT) with the B3LYP functional. For heats of formation, the findings are a mean absolute deviation from experiment of 1.2 kcal/mol for CBS‐Q, 2.0 kcal/mol for B3LYP, and 2.2 kcal/mol for DMC. The mean absolute deviation of 31 BDEs is 2.0 kcal/mol for CBS‐Q, 4.2 kcal/mol for B3LYP, and 2.5 kcal/mol for DMC. These findings are for 17 BDEs of closed‐shell molecules that have mean absolute deviations from experiment of 1.7 kcal/mol (CBS‐Q), 4.0 kcal/mol (B3LYP), and 2.2 kcal/mol (DMC). The corresponding results for the 14 BDEs of open‐shell molecules studied are 2.4 kcal/mol (CBS‐Q), 4.3 kcal/mol (B3LYP), and 2.9 kcal/mol (DMC). The DMC results provide a baseline from which improvement using multideterminant trial functions can be measured. © 2005 Wiley Periodicals, Inc. Int J Chem Kinet 37: 583–592, 2005