Branko Ruscic

W4 theory for computational thermochemistry: In pursuit of confident sub-kJ/mol predictions

In an attempt to improve on our earlier W3 theory [A. D. Boese et al., J. Chem. Phys. 120, 4129 (2004)] we consider such refinements as more accurate estimates for the contribution of connected quadruple excitations (T4), inclusion of connected quintuple excitations (T5), diagonal Born-Oppenheimer corrections (DBOC), and improved basis set extrapolation procedures. Revised experimental data for validation purposes were obtained from the latest version of the Active Thermochemical Tables thermochemical network. The recent CCSDT(Q) method offers a cost-effective way of estimating T4, but is insufficient by itself if the molecule exhibits some nondynamical correlation. The latter considerably slows down basis set convergence for T4, and anomalous basis set convergence in highly polar systems makes two-point extrapolation procedures unusable. However, we found that the CCSDTQ-CCSDT(Q) difference converges quite rapidly with the basis set, and that the formula 1.10[CCSDT(Q)cc-pVTZ+CCSDTQcc-pVDZ-CCSDT(Q)cc-pVDZ] offers a very reliable as well as fairly cost-effective estimate of the basis set limit T4 contribution. The T5 contribution converges very rapidly with the basis set, and even a simple double-zeta basis set appears to be adequate. The largest T5 contribution found in the present work is on the order of 0.5 kcal/mol (for ozone). DBOCs are significant at the 0.1 kcal/mol level in hydride systems. Post-CCSD(T) contributions to the core-valence correlation energy are only significant at that level in systems with severe nondynamical correlation effects. Based on the accumulated experience, a new computational thermochemistry protocol for first- and second-row main-group systems, to be known as W4 theory, is proposed. Its computational cost is not insurmountably higher than that of the earlier W3 theory, while performance is markedly superior. Our W4 atomization energies for a number of key species are in excellent agreement (better than 0.1 kcal/mol on average, 95% confidence intervals narrower than 1 kJ/mol) with the latest experimental data obtained from Active Thermochemical Tables. Lower-cost variants are proposed: the sequence W1-->W2.2-->W3.2-->W4lite-->W4 is proposed as a converging hierarchy of computational thermochemistry methods. A simple a priori estimate for the importance of post-CCSD(T) correlation contributions (and hence a pessimistic estimate for the error in a W2-type calculation) is proposed.

High-accuracy extrapolated ab initio thermochemistry. III. Additional improvements and overview

Effects of increased basis-set size as well as a correlated treatment of the diagonal Born-Oppenheimer approximation are studied within the context of the high-accuracy extrapolated ab initio thermochemistry (HEAT) theoretical model chemistry. It is found that the addition of these ostensible improvements does little to increase the overall accuracy of HEAT for the determination of molecular atomization energies. Fortuitous cancellation of high-level effects is shown to give the overall HEAT strategy an accuracy that is, in fact, higher than most of its individual components. In addition, the issue of core-valence electron correlation separation is explored; it is found that approximate additive treatments of the two effects have limitations that are significant in the realm of <1 kJ mol(-1) theoretical thermochemistry.

IUPAC Critical Evaluation of Thermochemical Properties of Selected Radicals. Part I

This is the first part of a series of articles reporting critically evaluated thermochemical properties of selected free radicals. The present article contains datasheets for 11 radicals: CH, CH2(triplet), CH2(singlet), CH3, CH2OH, CH3O, CH3CO, C2H5O, C6H5CH2, OH, and NH2. The thermochemical properties discussed are the enthalpy of formation, as well as the heat capacity, integrated heat capacity, and entropy of the radicals. One distinguishing feature of the present evaluation is the systematic utilization of available kinetic, spectroscopic and ion thermochemical data as well as high-level theoretical results.

Vacuum ultraviolet photoionization mass spectrometric study of C60

Gaseous C60 has been studied by photoionization mass spectrometry between the ionization threshold and 40.8 eV. An adiabatic threshold of 7.57±0.01 eV is observed, which may be slightly low due to hot bands. The energy derivative of the photoion yield curve is in rough agreement with the He i photoelectron spectrum on the positions of some peaks, but others are weak or absent. The discrepancy is not attributed to autoionization, but rather to selection rules governing the ejection of low energy electrons into high angular momentum waves. C++60 is observed at higher energies, and becomes ∼0.6 as intense as C+60 at 40.8 eV. The photoion yield curve of C++60, approximately linear well above threshold, appears to exhibit curvature near threshold, thwarting an attempt to make a distinction between two alternative values of the second ionization potential. Fragmentation to form C+58 is only observed at the highest energy, 40.8 eV. The unimolecular decay is modelled by quasiequilibrium theory. In this model, the...

Active Thermochemical Tables: thermochemistry for the 21st century

Active Thermochemical Tables (ATcT) are a good example of a significant breakthrough in chemical science that is directly enabled by the US DOE SciDAC initiative. ATcT is a new paradigm of how to obtain accurate, reliable, and internally consistent thermochemistry and overcome the limitations that are intrinsic to the traditional sequential approach to thermochemistry. The availability of high-quality consistent thermochemical values is critical in many areas of chemistry, including the development of realistic predictive models of complex chemical environments such as combustion or the atmosphere, or development and improvement of sophisticated high-fidelity electronic structure computational treatments. As opposed to the traditional sequential evolution of thermochemical values for the chemical species of interest, ATcT utilizes the Thermochemical Network (TN) approach. This approach explicitly exposes the maze of inherent interdependencies normally ignored by the conventional treatment, and allows, inter alia, a statistical analysis of the individual measurements that define the TN. The end result is the extraction of the best possible thermochemistry, based on optimal use of all the currently available knowledge, hence making conventional tabulations of thermochemical values obsolete. Moreover, ATcT offer a number of additional features that are neither present nor possible in the traditional approach. With ATcT, new knowledge can be painlessly propagated through all affected thermochemical values. ATcT also allows hypothesis testing and evaluation, as well as discovery of weak links in the TN. The latter provides pointers to new experimental or theoretical determinations that can most efficiently improve the underlying thermochemical body of knowledge.

Photochemical formation and transport of ozone in Athens, Greece

Abstract An evaluation of the diurnal variation of the hourly ozone concentrations measured at five sites in Greater Athens from June until early September 1984 indicates that photosmog episodes in Greater Athens are associated with the sea breeze circulation. Due to local air circulation in the Athens basin, precursors of O3 are transported to and accumulated in the Saronikos Bay during the morning hours while the land breeze is blowing. At noon, when the sea breeze sets in, the O3 formed over the sea is brought back to the coast and to central Athens where it increases the local O3 concentration by a factor of 3–5 within a few hours. The O3 levels often remain high throughout the night. During the photochemical smog episodes, all of them accompanied by well-developed sea breezes, the U.S. Air Quality Standard of 120 ppb O3 was exceeded for 4–7 h day−1. Peak O3 concentrations up to nearly 200 ppb were recorded in the smog episodes. Relatively high O3 concentrations were measured on the island of Aegina. They tend to remain high during the night and can be attributed only to primary pollutant transport from Greater Athens advected by the land breeze. The O3 values obtained at Mount Immitos (1000 m above MSL) suggest that, first, the sea breeze inhibits the influence of vertical thermal convection up to heights above 600 m, and second, no O3 is noticeable from above the tropopause to ground level or from long-range transport.

Rate Constants for the Thermal Decomposition of Ethanol and Its Bimolecular Reactions with OH and D: Reflected Shock Tube and Theoretical Studies

The thermal decomposition of ethanol and its reactions with OH and D have been studied with both shock tube experiments and ab initio transition state theory-based master equation calculations. Dissociation rate constants for ethanol have been measured at high T in reflected shock waves using OH optical absorption and high-sensitivity H-atom ARAS detection. The three dissociation processes that are dominant at high T are C2H5OH--> C2H4+H2O (A) -->CH3+CH2OH (B) -->C2H5+OH (C).The rate coefficient for reaction C was measured directly with high sensitivity at 308 nm using a multipass optical White cell. Meanwhile, H-atom ARAS measurements yield the overall rate coefficient and that for the sum of reactions B and C , since H-atoms are instantaneously formed from the decompositions of CH(2)OH and C(2)H(5) into CH(2)O + H and C(2)H(4) + H, respectively. By difference, rate constants for reaction 1 could be obtained. One potential complication is the scavenging of OH by unreacted ethanol in the OH experiments, and therefore, rate constants for OH+C2H5OH-->products (D)were measured using tert-butyl hydroperoxide (tBH) as the thermal source for OH. The present experiments can be represented by the Arrhenius expression k=(2.5+/-0.43) x 10(-11) exp(-911+/-191 K/T) cm3 molecule(-1) s(-1) over the T range 857-1297 K. For completeness, we have also measured the rate coefficient for the reaction of D atoms with ethanol D+C2H5OH-->products (E) whose H analogue is another key reaction in the combustion of ethanol. Over the T range 1054-1359 K, the rate constants from the present experiments can be represented by the Arrhenius expression, k=(3.98+/-0.76) x10(-10) exp(-4494+/-235 K/T) cm3 molecule(-1) s(-1). The high-pressure rate coefficients for reactions B and C were studied with variable reaction coordinate transition state theory employing directly determined CASPT2/cc-pvdz interaction energies. Reactions A , D , and E were studied with conventional transition state theory employing QCISD(T)/CBS energies. For the saddle point in reaction A , additional high-level corrections are evaluated. The predicted reaction exo- and endothermicities are in good agreement with the current Active Thermochemical Tables values. The transition state theory predictions for the microcanonical rate coefficients in ethanol decomposition are incorporated in master equation calculations to yield predictions for the temperature and pressure dependences of reactions A - C . With modest adjustments (<1 kcal/mol) to a few key barrier heights, the present experimental and adjusted theoretical results yield a consistent description of both the decomposition (1-3) and abstraction kinetics (4 and 5). The present results are compared with earlier experimental and theoretical work.

The ethyl radical: photoionization and theoretical studies

The ethyl radical has been produced by the F+C2H6 reaction, and studied by photoionization mass spectrometry. The adiabatic ionization potential (IP) is found to be 8.117±0.008 eV. The radical and cation have been studied by the G1 level of molecular orbital theory. The calculated adiabatic ionization potential is 8.06 eV. On the basis of both experiment and theory, it is concluded that the cation in its ground state has a bridged structure. The adiabatic value resolves a previous inconsistency involving ΔH○f(C2H5), ΔH○f(C2H+5) and IP (C2H5). Broad autoionizing structures are observed at higher energy. These features can be rationalized as Rydberg states converging to the first excited state 3A‘ calculated to be at 12.05 eV.

Photoionization studies of GeHn (n=2-4)

The adiabatic ionization potential of GeH+4 (GeH4) is measured by photoionization mass spectrometry to be ≤10.53 eV and perhaps as low as 10.44 eV. This is about 0.8 eV (∼9 vibrational quanta) lower than the value reported by photoelectron spectroscopy. This result, analogous to that found for SiH+4 (SiH4), implies a marked Jahn–Teller distortion of GeH+4. The appearance potentials of GeH+2 and GeH+3 from GeH4 are ≤10.772±0.009 eV and 56.4 (59); D0 (HGe–H) 53.7 (63). These specific bond energies, when appropriately normalized, display the same pattern as the SiHn (but not the CHn) bond energies and provid...

Heats of formation of C6H5?, C6H5+, and C6H5NO by threshold photoelectron photoion coincidence and active thermochemical tables analysis

Threshold photoelectron photoion coincidence has been used to prepare selected internal energy distributions of nitrosobenzene ions [C(6)H(5)NO(+)]. Dissociation to C(6)H(5)(+) + NO products was measured over a range of internal energies and rate constants from 10(3) to 10(7) s(-1) and fitted with the statistical theory of unimolecular decay. A 0 K dissociative photoionization onset energy of 10.607 ± 0.020 eV was derived by using the simplified statistical adiabatic channel model. The thermochemical network of Active Thermochemical Tables (ATcT) was expanded to include phenyl and phenylium, as well as nitrosobenzene. The current ATcT heats of formation of these three species at 0 K (298.15 K) are 350.6 (337.3) ± 0.6, 1148.7 (1136.8) ± 1.0, and 215.6 (198.6) ± 1.5 kJ mol(-1), respectively. The resulting adiabatic ionization energy of phenyl is 8.272 ± 0.010 eV. The new ATcT thermochemistry for phenyl entails a 0 K (298.15 K) C-H bond dissociation enthalpy of benzene of 465.9 (472.1) ± 0.6 kJ mol(-1). Several related thermochemical quantities from ATcT, including the current enthalpies of formation of benzene, monohalobenzenes, and their ions, as well as interim ATcT values for the constituent atoms, are also given.