07.04.2013 Views

Essentials of Computational Chemistry

Essentials of Computational Chemistry

Essentials of Computational Chemistry

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

382 10 THERMODYNAMIC PROPERTIES<br />

Table 10.4 Predicted 298 K enthalpies <strong>of</strong> formation (kcal mol −1 ) for hydroxylamine<br />

From<br />

Level <strong>of</strong> theory Eq. (10.61) Eq. (10.62) Eq. (10.63)<br />

AM1 −32.34 −31.31<br />

HF/cc-pVDZ −12.14 −12.02<br />

HF/cc-pVQZ −8.83 −13.06<br />

B3LYP/aug-cc-pVTZ −12.18 −9.69<br />

MP2/cc-pVQZ −8.61 −12.09<br />

CCSD(T)/cc-pVQZ// −11.56 −10.61<br />

CCSD(T)/cc-pVDZ<br />

BAC-MP4 −12.98 −11.09<br />

G2 −10.60 −11.78 −11.53<br />

G2MP2 −11.46 −11.69 −11.67<br />

G3B3 −10.02 −11.51 −11.53<br />

G3 −9.36 −11.15 −11.28<br />

CBS-Q −10.03 −12.18 −11.16<br />

Statistical a −10.29 ± 0.70 −11.66 ± 0.34 −11.43 ± 0.19<br />

a Average ± standard deviation from G2, G2MP2, G3B3, G3, and CBS-Q.<br />

calculation far exceeds every other entry in the table. Since the G2, G2MP2, G3B3, G3,<br />

and CBS-Q models (all discussed in Chapter 7) are cheaper than CCSD(T)/cc-pVQZ and<br />

moreover designed specifically for the purpose <strong>of</strong> computing enthalpies <strong>of</strong> formation, there<br />

is ample reason to focus more closely on their performance.<br />

Rather than attempting to rationalize why any one <strong>of</strong> these composite levels might<br />

be more or less good than another, let us examine their joint performance. The final<br />

row <strong>of</strong> Table 10.4 provides the means and standard deviations <strong>of</strong> the predicted H o<br />

f,298<br />

(H2NOH) values from these levels for Eqs. (10.61) to (10.63). The largest standard<br />

deviation is associated with Eq. (10.61), the next largest with Eq. (10.62), and the smallest,<br />

only 0.19 kcal mol −1 , with Eq. (10.63). This trend is entirely consistent with the above<br />

discussion <strong>of</strong> the relative quality <strong>of</strong> the three isodesmic equations, and provides some<br />

quantitative feel for how difficult the accurate computation <strong>of</strong> an atomization energy really<br />

is. Given this analysis, it appears reasonable to take the average value from the last five<br />

methods and Eq. (10.63) as a best estimate: −11.4 kcal mol −1 . Further support for this<br />

choice comes from considering a different reaction, namely<br />

H2 + H2O2 ⇀↽ 2H2O (10.64)<br />

Note that this is the analog to Eq. (10.62) with hydrogen peroxide replacing hydroxylamine.<br />

In this case, all enthalpies <strong>of</strong> formation are known experimentally to high accuracy, so the<br />

performance <strong>of</strong> the various theoretical models may be directly assessed. Applying the same<br />

averaging procedure, one finds that the models predict an enthalpy <strong>of</strong> formation for H2O2<br />

that is too negative by 0.3 kcal mol −1 . Note that if one assumes that this correction may be<br />

applied to the results from Eq. (10.62) for hydroxylamine, one predicts −11.4 kcal mol −1 ,<br />

in perfect agreement with the uncorrected results from Eq. (10.63).<br />

Note that the average atomization energy prediction differs from −11.4 kcal mol −1 by<br />

only 1.1 kcal mol −1 , which is about the range <strong>of</strong> accuracy typically quoted for the models

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!