Thermal decomposition of HO2NO2 (peroxynitric acid, PNA): rate coefficient and determination of the enthalpy of formation

Rate coefficients for the gas-phase thermal decomposition of HO(2)NO(2) (peroxynitric acid, PNA) are reported at temperatures between 331 and 350 K at total pressures of 25 and 50 Torr of N(2). Rate coefficients were determined by measuring the steady-state OH concentration in a mixture of known concentrations of HO(2)NO(2) and NO. The measured thermal decomposition rate coefficients k(-)(1)(T,P) are used in combination with previously published rate coefficient data for the HO(2)NO(2) formation reaction to yield a standard enthalpy for reaction 1 of Delta(r)H degrees (298K) = -24.0 +/- 0.5 kcal mol(-1) (uncertainties are 2sigma values and include estimated systematic errors). A HO(2)NO(2) standard heat of formation, Delta(f)H degrees (298K)(HO(2)NO(2)), of -12.6 +/- 1.0 kcal mol(-1) was calculated from this value. Some of the previously reported data on the thermal decomposition of HO(2)NO(2) have been reanalyzed and shown to be in good agreement with our reported value.     

Establishment of the C(2)H(5)+O(2) reaction mechanism: a combustion archetype

The celebrated C(2)H(5)+O(2) reaction is an archetype for hydrocarbon combustion, and the critical step in the process is the concerted elimination of HO(2) from the ethylperoxy intermediate (C(2)H(5)O(2)). Master equation kinetic models fitted to measured reaction rates place the concerted elimination barrier 3.0 kcal mol(-1) below the C(2)H(5)+O(2) reactants, whereas the best previous electronic structure computations yield a barrier more than 2.0 kcal mol(-1) higher. We resolve this discrepancy here by means of the most rigorous computations to date, using focal point methods to converge on the ab initio limit. Explicit computations were executed with basis sets as large as cc-pV5Z and correlation treatments as extensive as coupled cluster through full triples with a perturbative inclusion of quadruple excitations [CCSDT(Q)]. The final predicted barrier is -3.0 kcal mol(-1), bringing the concerted elimination mechanism into precise agreement with experiment. This work demonstrates that higher correlation treatments such as CCSDT(Q) are not only feasible on systems of chemical interest but are necessary to supply accuracy beyond 0.5 kcal mol(-1), which is not obtained with the "gold standard" CCSD(T) method. Finally, we compute the enthalpy of formation of C(2)H(5)O(2) to be Delta(f)H degrees (298 K)=-5.3+/-0.5 kcal mol(-1) and Delta(f)H degrees (0 K)=-1.5+/-0.5 kcal mol(-1).     

Thermodynamic and ab initio analysis of the controversial enthalpy of formation of formaldehyde.

There are two values, -26.0 and -27.7 kcal mol(-1), that are routinely reported in literature evaluations for the standard enthalpy of formation, Delta(f) H(o)(298), of formaldehyde (CH(2)=O), where error limits are less than the difference in values. In this study, we summarize the reported literature for formaldehyde enthalpy values based on evaluated measurements and on computational studies. Using experimental reaction enthalpies for a series of reactions involving formaldehyde, in conjunction with known enthalpies of formation, its enthalpy is determined to be -26.05+/-0.42 kcal mol(-1), which we believe is the most accurate enthalpy currently available. For the same reaction series, the reaction enthalpies are evaluated using six computational methods: CBS-Q, CBS-Q//B3, CBS-APNO, G2, G3, and G3B3 yield Delta(f) H(o)(298)=-25.90+/-1.17 kcal mol(-1), which is in good agreement to our experimentally derived result. Furthermore, the computational chemistry methods G3, G3MP2B3, CCSD/6-311+G(2df,p)//B3LYP/6-31G(d), CCSD(T)/6-311+G(2df,p)//B3LYP/6-31G(d), and CBS-APNO in conjunction with isodesmic and homodesmic reactions are used to determine Delta(f) H(o)(298). Results from a series of five work reactions at the higher levels of calculation are -26.30+/-0.39 kcal mol(-1) with G3, -26.45+/-0.38 kcal mol(-1) with G3MP2B3, -26.09+/-0.37 kcal mol(-1) with CBS-APNO, -26.19+/-0.48 kcal mol(-1) with CCSD, and -26.16+/-0.58 kcal mol(-1) with CCSD(T). Results from heat of atomization calculations using seven accurate ab initio methods yields an enthalpy value of -26.82+/-0.99 kcal mol(-1). The results using isodesmic reactions are found to give enthalpies more accurate than both other computational approaches and are of similar accuracy to atomization enthalpy calculations derived from computationally intensive W1 and CBS-APNO methods. Overall, our most accurate calculations provide an enthalpy of formation in the range of -26.2 to -26.7 kcal mol(-1), which is within computational error of the suggested experimental value. The relative merits of each of the three computational methods are discussed and depend upon the accuracy of experimental enthalpies of formation required in the calculations and the importance of systematic computational errors in the work reaction. Our results also calculate Delta(f) H(o)(298) for the formyl anion (HCO(-)) as 1.28+/-0.43 kcal mol(-1).     

Solvation properties of N-substituted cis and trans amides are not identical: significant enthalpy and entropy changes are revealed by the use of variable temperature 1H NMR in aqueous and chloroform solutions and ab initio calculations

The cis/trans conformational equilibrium of N-methyl formamide (NMF) and the sterically hindered tert-butylformamide (TBF) was investigated by the use of variable temperature gradient 1H NMR in aqueous solution and in the low dielectric constant and solvation ability solvent CDCl3 and various levels of first principles calculations. The trans isomer of NMF in aqueous solution is enthalpically favored relative to the cis (deltaH(o) = -5.79 +/- 0.18 kJ mol(-1)) with entropy differences at 298 K (298 x deltaS(o) = -0.23 +/- 0.17 kJ mol(-1)) playing a minor role. The experimental value of the enthalpy difference strongly decreases (deltaH(o) = -1.72 +/- 0.06 kJ mol(-1)), and the contribution of entropy at 298 K (298 x deltaS(o) = -1.87 +/- 0.06 kJ mol(-1)) increases in the case of the sterically hindered tert-butylformamide. The trans isomer of NMF in CDCl3 solution is enthalpically favored relative to the cis (deltaH(o) = -3.71 +/- 0.17 kJ mol(-1)) with entropy differences at 298 K (298 x deltaS(o) = 1.02 +/- 0.19 kJ mol(-1)) playing a minor role. In the sterically hindered tert-butylformamide, the trans isomer is enthalpically disfavored (deltaH(o) = 1.60 +/- 0.09 kJ mol(-1)) but is entropically favored (298 x deltaS(o) = 1.71 +/- 0.10 kJ mol(-1)). The results are compared with literature data of model peptides. It is concluded that, in amide bonds at 298 K and in the absence of strongly stabilizing sequence-specific inter-residue interactions involving side chains, the free energy difference of the cis/trans isomers and both the enthalpy and entropy contributions are strongly dependent on the N-alkyl substitution and the solvent. The significant decreasing enthalpic benefit of the trans isomer in CDCl3 compared to that in H2O, in the case of NMF and TBF, is partially offset by an adverse entropy contribution. This is in agreement with the general phenomenon of enthalpy versus entropy compensation. B3LY/6-311++G** and MP2/6-311++G** quantum chemical calculations confirm the stability orders of isomers and the deltaG decrease in going from water to CHCl3 as solvent. However, the absolute calculated values, especially for TBF, deviate significantly from the experimental values. Consideration of the solvent effects via the PCM approach on NMF x H2O and TBF x H2O supermolecules improves the agreement with the experimental results for TBF isomers, but not for NMF.     

Theoretical prediction of the heats of formation of C2H5O* radicals derived from ethanol and of the kinetics of beta-C-C scission in the ethoxy radical

Thermochemical parameters of three C(2)H(5)O* radicals derived from ethanol were reevaluated using coupled-cluster theory CCSD(T) calculations, with the aug-cc-pVnZ (n = D, T, Q) basis sets, that allow the CC energies to be extrapolated at the CBS limit. Theoretical results obtained for methanol and two CH(3)O* radicals were found to agree within +/-0.5 kcal/mol with the experiment values. A set of consistent values was determined for ethanol and its radicals: (a) heats of formation (298 K) DeltaHf(C(2)H(5)OH) = -56.4 +/- 0.8 kcal/mol (exptl: -56.21 +/- 0.12 kcal/mol), DeltaHf(CH(3)C*HOH) = -13.1 +/- 0.8 kcal/mol, DeltaHf(C*H(2)CH(2)OH) = -6.2 +/- 0.8 kcal/mol, and DeltaHf(CH(3)CH(2)O*) = -2.7 +/- 0.8 kcal/mol; (b) bond dissociation energies (BDEs) of ethanol (0 K) BDE(CH(3)CHOH-H) = 93.9 +/- 0.8 kcal/mol, BDE(CH(2)CH(2)OH-H) = 100.6 +/- 0.8 kcal/mol, and BDE(CH(3)CH(2)O-H) = 104.5 +/- 0.8 kcal/mol. The present results support the experimental ionization energies and electron affinities of the radicals, and appearance energy of (CH(3)CHOH+) cation. Beta-C-C bond scission in the ethoxy radical, CH(3)CH2O*, leading to the formation of C*H3 and CH(2)=O, is characterized by a C-C bond energy of 9.6 kcal/mol at 0 K, a zero-point-corrected energy barrier of E0++ = 17.2 kcal/mol, an activation energy of Ea = 18.0 kcal/mol and a high-pressure thermal rate coefficient of k(infinity)(298 K) = 3.9 s(-1), including a tunneling correction. The latter value is in excellent agreement with the value of 5.2 s(-1) from the most recent experimental kinetic data. Using RRKM theory, we obtain a general rate expression of k(T,p) = 1.26 x 10(9)p(0.793) exp(-15.5/RT) s(-1) in the temperature range (T) from 198 to 1998 K and pressure range (p) from 0.1 to 8360.1 Torr with N2 as the collision partners, where k(298 K, 760 Torr) = 2.7 s(-1), without tunneling and k = 3.2 s(-1) with the tunneling correction. Evidence is provided that heavy atom tunneling can play a role in the rate constant for beta-C-C bond scission in alkoxy radicals.     

Thermochemistry of methyl and ethyl nitro, RNO2, and nitrite, RONO, organic compounds

Computational quantum theory is employed to determine the thermochemical properties of n-alkyl nitro and nitrite compounds: methyl and ethyl nitrites, CH3ONO and C2H5ONO, plus nitromethane and nitroethane, CH3NO2 and C2H5NO2, at 298.15 K using multilevel G3, CBS-QB3, and CBS-APNO composite methods employing both atomization and isodesmic reaction analysis. Structures and enthalpies of the corresponding aci-tautomers are also determined. The enthalpies of formation for the most stable conformers of methyl and ethyl nitrites at 298 K are determined to be -15.64 +/- 0.10 kcal mol-1 (-65.44 +/- 0.42 kJ mol-1) and -23.58 +/- 0.12 kcal mol-1 (-98.32 +/- 0.58 kJ mol-1), respectively. DeltafHo(298 K) of nitroalkanes are correspondingly evaluated at -17.67 +/- 0.27 kcal mol-1 (-74.1 +/- 1.12 kJ mol-1) and -25.06 +/- 0.07 kcal mol-1 (-121.2 +/- 0.29 kJ mol-1) for CH3NO2 and C2H5NO2. Enthalpies of formation for the aci-tautomers are calculated as -3.45 +/- 0.44 kcal mol-1 (-14.43 +/- 0.11 kJ mol-1) for aci-nitromethane and -14.25 +/- 0.44 kcal mol-1 (-59.95 +/- 1.84 kJ mol-1) for the aci-nitroethane isomers, respectively. Data are evaluated against experimental and computational values in the literature with recommendations. A set of thermal correction parameters to atomic (H, C, N, O) enthalpies at 0 K is developed, to enable a direct calculation of species enthalpy of formation at 298.15 K, using atomization reaction and computation outputs.     

Experimental and ab initio study of the HO2.CH3OH complex: thermodynamics and kinetics of formation

Near-infrared spectroscopy was used to monitor HO2 formed by pulsed laser photolysis of Cl2-O2-CH3OH-N2 mixtures. On the microsecond time scale, [HO2] exhibited a time dependence consistent with a mechanism in which [HO2] approached equilibrium via HO2 + HO2.CH3OH (3, -3). The equilibrium constant for reaction 3, K(p), was measured between 231 and 261 K at 50 and 100 Torr, leading to standard reaction enthalpy and entropy values (1 sigma) of delta(r) = -37.4 +/- 4.8 kJ mol(-1) and delta(r) = -100 +/- 19 J mol(-1) K(-1). The effective bimolecular rate constant, k3, for formation of the HO2.CH3OH complex is .10(-15).exp[(1800 +/- 500)/T] cm3 molecule(-1) s(-1) at 100 Torr (1 sigma). Ab initio calculations of the optimized structure and energetics of the HO2.CH3OH complex were performed at the CCSD(T)/6-311++G(3df,3pd)//MP2(full)/6-311++G(2df,2pd) level. The complex was found to have a strong hydrogen bond (D(e) = 43.9 kJ mol(-1)) with the hydrogen in HO2 binding to the oxygen in CH3OH. The calculated enthalpy for association is delta(r) = -36.8 kJ mol(-1). The potentials for the torsion about the O2-H bond and for the hydrogen-bond stretch were computed and 1D vibrational levels determined. After explicitly accounting for these degrees of freedom, the calculated Third Law entropy of association is delta(r) = -106 J mol(-1) K(-1). Both the calculated enthalpy and entropy of association are in reasonably good agreement with experiment. When combined with results from our previous study (Christensen et al. Geophys. Res. Lett. 2002, 29; doi:10.1029/2001GL014525), the rate coefficient for the reaction of HO2 with the complex, HO2 + HO2.CH3OH, is determined to be (2.1 +/- 0.7) x 10(-11) cm3 molecule(-1) s(-1). The results of the present work argue for a reinterpretation of the recent measurement of the HO2 self-reaction rate constant by Stone and Rowley (Phys. Chem. Chem. Phys. 2005, 7, 2156). Significant complex concentrations are present at the high methanol concentrations used in that work and lead to a nonlinear methanol dependence of the apparent rate constant. This nonlinearity introduces substantial uncertainty in the extrapolation to zero methanol.     

Thermochemistry of the fluoroformyloxyl radical: a computational study based on coupled cluster theory

The standard enthalpy of formation of FCO(2) (X (2)B(2)) was determined by a computational approach based on coupled cluster theory [CCSD(T)] with energies extrapolated to the basis-set limit, with additional corrections accounting for core-valence correlation, scalar relativity, spin-orbit coupling, and zero-point vibrational motions. Utilizing a variety of independent reaction schemes, our best estimate is Delta(f)H(o)(0)(FCO(2)) = -86.0 +/- 0.6 kcal mol(-1) [Delta(f)H(o)(298) )(FCO(2)) = -86.7 +/- 0.6 kcal mol(-1)], which is shown to be more accurate than previous theoretical and experimental values. The chosen computational procedure was also applied to HCO (X (2)A'), where we find excellent agreement with experiment, and to FCO (X (2)A'), where we recommend an improved value of Delta(f)H(o)(0)(FCO) = -42.1 +/- 0.5 kcal mol(-1) [ Delta(f)H(o)(298)(FCO) = -42.0 +/- 0.5 kcal mol(-1)]. Further theoretical results concern the C-F bond dissociation energy, electron affinity, ionization energy, first and second excitation energies in FCO(2), fluoride ion affinity of CO(2), and equilibrium geometries of the molecules treated presently. For FCO (X (2)A') we propose an improved equilibrium structure: r(e)(CF) = 132.5(2) pm, r(e)(CO) = 116.7(2) pm, and theta(e)(FCO) = 127.8(2)(o).     

Formation of nitric acid in the gas-phase HO2 + NO reaction: effects of temperature and water vapor

A high-pressure turbulent flow reactor coupled with a chemical ionization mass spectrometer was used to investigate the minor channel (1b) producing nitric acid, HNO3, in the HO2 + NO reaction for which only one channel (1a) is known so far: HO2 + NO --> OH + NO2 (1a), HO2 + NO --> HNO3 (1b). The reaction has been investigated in the temperature range 223-298 K at a pressure of 200 Torr of N2 carrier gas. The influence of water vapor has been studied at 298 K. The branching ratio, k1b/k1a, was found to increase from (0.18(+0.04/-0.06))% at 298 K to (0.87(+0.05/-0.08))% at 223 K, corresponding to k1b = (1.6 +/- 0.5) x 10(-14) and (10.4 +/- 1.7) x 10(-14) cm3 molecule(-1) s(-1), respectively at 298 and 223 K. The data could be fitted by the Arrhenius expression k1b = 6.4 x 10(-17) exp((1644 +/- 76)/T) cm3 molecule(-1) s(-1) at T = 223-298 K. The yield of HNO3 was found to increase in the presence of water vapor (by 90% at about 3 Torr of H2O). Implications of the obtained results for atmospheric radicals chemistry and chemical amplifiers used to measure peroxy radicals are discussed. The results show in particular that reaction 1b can be a significant loss process for the HO(x) (OH, HO2) radicals in the upper troposphere.     

Accurate computational determination of the binding energy of the SO3 x H2O complex

Reliable thermochemical data for the reaction SO3 + H2O<-->SO3 x H2O (1a) are of crucial importance for an adequate modeling of the homogeneous H2SO4 formation in the atmosphere. We report on high-level quantum chemical calculations to predict the binding energy of the SO3 x H2O complex. The electronic binding energy is accurately computed to De = 40.9+/-1.0 kJ/mol = 9.8+/-0.2 kcal/mol. By using harmonic frequencies from density functional theory calculations (B3LYP/cc-pVTZ and TPSS/def2-TZVP), zero-point and thermal energies were calculated. From these data, we estimate D0 = -Delta H(1a)0(0 K) = 7.7+/-0.5 kcal/mol and Delta H(1a)0(298 K) = -8.3+/-1.0 kcal/mol.     

Thermochemical parameters of CHFO and CF2O

Thermochemical properties of CHFO and CF 2O and their derivatives were calculated by using coupled-cluster theory (U)CCSD(T) calculations with the aug-cc-pV nZ ( n = D, T, Q, 5) basis sets extrapolated to the complete basis set limit with additional corrections. The predicted properties include the following. Enthalpies of formation (298 K, kcal/mol): Delta H f (CF 2O) = -144.7, Delta H f(CHFO) = -91.1, Delta H f (CFO (*)) = -41.6. Bond dissociation energy (0 K, kcal/mol): BDE(CFO-F) = 120.7, BDE(CHO-F) = 119.1, BDE(CFO-H) = 100.2. Ionization potential (eV): IP 1(CF 2O) = 13.04, IP 2(CF 2O) = 14.09, IP 1(CHFO) = 12.41, IP 2(CHFO) = 13.99, IP 1(CFO (*)) = 9.34. Proton affinity (298 K, kcal/mol), PA O(CF 2O) = 148.8, PA O(CHFO) = 156.7, PA F(CHFO) = 154.5 kcal/mol. Electron affinity: EA(CFO (*)) = 2.38 eV. Triplet-singlet separation gap (eV): Delta E T1-S0(CF 2O) = 4.47, Delta E T1-S0(CHFO) = 4.36. Triplet-triplet transition energy (eV): Delta E T2-T1(CF 2O) = 0.44. The new calculated values contribute to solving some persistent discrepancies in the literature. The effects of F-atoms on thermochemical parameters are not linearly additive, and the changes are largely dominated by the first F-substitution. On the basis of the calculated proton affinities of CF 2O and CF 3OH, the nucleophilicities of the oxygen atoms are, within computational errors, the same in both compounds.     

Heats of formation of triplet ethylene, ethylidene, and acetylene

Heats of formation of the lowest triplet state of ethylene and the ground triplet state of ethylidene have been predicted by high level electronic structure calculations. Total atomization energies obtained from coupled-cluster CCSD(T) energies extrapolated to the complete basis set limit using correlation consistent basis sets (CBS), plus additional corrections predict the following heats of formation in kcal/mol: DeltaH0r(C2H4,3A1) = 80.1 at 0 K and 78.5 at 298 K, and DeltaH0t(CH3CH,3A' ') = 86.8 at 0 K and 85.1 at 298 K, with an error of less than +/-1.0 kcal/mol. The vertical and adiabatic singlet-triplet separation energies of ethylene were calculated as DeltaES-T,vert = 104.1 and DeltaES-T,adia = 65.8 kcal/mol. These results are in excellent agreement with recent quantum Monte Carlo (DMC) values of 103.5 +/- 0.3 and 66.4 +/- 0.3 kcal/mol. Both sets of computational values differ from the experimental estimate of 58 +/- 3 kcal/mol for the adiabatic splitting. The computed singlet-triplet gap at 0 K for acetylene is DeltaES-T,adia(C2H2) = 90.5 kcal/mol, which is in notable disagreement with the experimental value of 82.6 kcal/mol. The heat of formation of the triplet is DeltaH0tC2H2,3B2) = 145.3 kcal/mol. There is a systematic underestimation of the singlet-triplet gaps in recent photodecomposition experiments by approximately 7 to 8 kcal/mol. For vinylidene, we predict DeltaH0t(H2CC,1A1) = 98.8 kcal/mol at 298 K (exptl. 100.3 +/- 4.0), DeltaH0t(H2CC,3B2) = 146.2 at 298 K, and an energy gap DeltaES-T-adia(H2CC) = 47.7 kcal/mol.     

Heats of formation and singlet-triplet separations of hydroxymethylene and 1-hydroxyethylidene

Thermochemical parameters of hydroxymethylene (HC:OH) and 1-hydroxyethylidene (CH3C:OH) were evaluated by using coupled-cluster, CCSD(T), theory, in conjunction with the augmented correlation consistent, aug-cc-pVnZ, basis sets, with n = D, T, Q, and 5, extrapolated to the complete basis set limit. The predicted value at 298 K for Delta Hf(CH2O) is -26.0 +/- 1 kcal/mol, as compared to an experimental value of -25.98 +/- 0.01 kcal/mol, and for Delta Hf(CH:OH) it is 26.1 +/- 1 kcal/mol. The hydroxymethylene-formaldehyde energy gap is 52.1 +/- 0.5 kcal/mol, the singlet-triplet separation of hydroxymethylene is Delta E(ST)(HC:OH) = 25.3 +/- 0.5 kcal/mol, the proton affinity is PA(HC:OH) = 222.5 +/- 0.5 kcal/mol, and the ionization energy is IEa(HC:OH) = 8.91 +/- 0.04 eV. The predicted value at 298 K for Delta Hf(CH3CHO) is -39.1 +/- 1 kcal/mol as compared to an experimental value of -40.80 +/- 0.35 kcal/mol, and for Delta Hf(CH3C:OH) it is 11.2 +/- 1 kcal/mol. The hydroxyethylidene-acetaldehyde energy gap is 50.6 +/- 0.5 kcal/mol, the singlet-triplet separation of 1-hydroxyethylidene is Delta E(ST)(CH3C:OH) = 30.5 +/- 0.5 kcal/mol, the proton affinity is PA(CH3C:OH) = 234.7 +/- 0.5 kcal/mol, and the ionization energy is IEa(CH3C:OH) = 8.18 +/- 0.04 eV. The calculated energy differences between the carbene and aldehyde isomers, and, thus, the heats of formation of the carbenes, differ from the experimental values by 2.5 kcal/mol.     

Kinetics, mechanism, and thermochemistry of the gas phase reaction of atomic chlorine with dimethyl sulfoxide

A laser flash photolysis-resonance fluorescence technique has been employed to study the kinetics of the reaction of chlorine atoms with dimethyl sulfoxide (CH3S(O)CH3; DMSO) as a function of temperature (270-571 K) and pressure (5-500 Torr) in nitrogen bath gas. At T = 296 K and P > or = 5 Torr, measured rate coefficients increase with increasing pressure. Combining our data with literature values for low-pressure rate coefficients (0.5-3 Torr He) leads to a rate coefficient for the pressure independent H-transfer channel of k1a = 1.45 x 10(-11) cm3 molecule(-1) s(-1) and the following falloff parameters for the pressure-dependent addition channel in N2 bath gas: k(1b,0) = 2.53 x 10(-28) cm6 molecule(-2) s(-1); k(1b,infinity) = 1.17 x 10(-10) cm3 molecule(-1) s(-1), F(c) = 0.503. At the 95% confidence level, both k1a and k1b(P) have estimated accuracies of +/-30%. At T > 430 K, where adduct decomposition is fast enough that only the H-transfer pathway is important, measured rate coefficients are independent of pressure (30-100 Torr N2) and increase with increasing temperature. The following Arrhenius expression adequately describes the temperature dependence of the rate coefficients measured at over the range 438-571 K: k1a = (4.6 +/- 0.4) x 10(-11) exp[-(472 +/- 40)/T) cm3 molecule(-1) s(-1) (uncertainties are 2sigma, precision only). When our data at T > 430 K are combined with values for k1a at temperatures of 273-335 K that are obtained by correcting reported low-pressure rate coefficients from discharge flow studies to remove the contribution from the pressure-dependent channel, the following modified Arrhenius expression best describes the derived temperature dependence: k1a = 1.34 x 10(-15)T(1.40) exp(+383/T) cm3 molecule(-1) s(-1) (273 K < or = T < or = 571 K). At temperatures around 330 K, reversible addition is observed, thus allowing equilibrium constants for Cl-DMSO formation and dissociation to be determined. A third-law analysis of the equilibrium data using structural information obtained from electronic structure calculations leads to the following thermochemical parameters for the association reaction: delta(r)H(o)298 = -72.8 +/- 2.9 kJ mol(-1), deltaH(o)0 = -71.5 +/- 3.3 kJ mol(-1), and delta(r)S(o)298 = -110.6 +/- 4.0 J K(-1) mol(-1). In conjunction with standard enthalpies of formation of Cl and DMSO taken from the literature, the above values for delta(r)H(o) lead to the following values for the standard enthalpy of formation of Cl-DMSO: delta(f)H(o)298 = -102.7 +/- 4.9 kJ mol(-1) and delta(r)H(o)0 = -84.4 +/- 5.8 kJ mol(-1). Uncertainties in the above thermochemical parameters represent estimated accuracy at the 95% confidence level. In agreement with one published theoretical study, electronic structure calculations using density functional theory and G3B3 theory reproduce the experimental adduct bond strength quite well.     

A molecular dynamics study of the thermodynamic properties of calcium apatites. 1. Hexagonal phases

Structural and thermodynamic properties of crystal hexagonal calcium apatites, Ca10(PO4)6(X)2 (X = OH, F, Cl, Br), were investigated using an all-atom Born-Huggins-Mayer potential by a molecular dynamics technique. The accuracy of the model at room temperature and atmospheric pressure was checked against crystal structural data, with maximum deviations of ca. 4% for the haloapatites and 8% for hydroxyapatite. The standard molar lattice enthalpy, delta(lat)H298(o), of the apatites was calculated and compared with previously published experimental results, the agreement being better than 2%. The molar heat capacity at constant pressure, C(p,m), in the range 298-1298 K, was estimated from the plot of the molar enthalpy of the crystal as a function of temperature, H(m) = (H(m,298) - 298C(p,m)) + C(p,m)T, yielding C(p,m) = 694 +/- 68 J x mol(-1) x K(-1), C(p,m) = 646 +/- 26 J x mol(-1) x K(-1), C(p,m) = 530 +/- 34 J x mol(-1) x K(-1), and C(p,m) = 811 +/- 42 J x mol(-1) x K(-1) for hydroxy-, fluor-, chlor-, and bromapatite, respectively. High-pressure simulation runs, in the range 0.5-75 kbar, were performed in order to estimate the isothermal compressibility coefficient, kappaT, of those compounds. The deformation of the compressed solids is always elastically anisotropic, with BrAp exhibiting a markedly different behavior from those displayed by HOAp and ClAp. High-pressure p-V data were fitted to the Parsafar-Mason equation of state with an accuracy better than 1%.     

Reaction mechanism between carbonyl oxide and hydroxyl radical: a theoretical study

The reaction mechanism of carbonyl oxide with hydroxyl radical was investigated by using CASSCF, B3LYP, QCISD, CASPT2, and CCSD(T) theoretical approaches with the 6-311+G(d,p), 6-311+G(2df, 2p), and aug-cc-pVTZ basis sets. This reaction involves the formation of H2CO + HO2 radical in a process that is computed to be exothermic by 57 kcal/mol. However, the reaction mechanism is very complex and begins with the formation of a pre-reactive hydrogen-bonded complex and follows by the addition of HO radical to the carbon atom of H2COO, forming the intermediate peroxy-radical H2C(OO)OH before producing formaldehyde and hydroperoxy radical. Our calculations predict that both the pre-reactive hydrogen-bonded complex and the transition state of the addition process lie energetically below the enthalpy of the separate reactants (DeltaH(298K) = -6.1 and -2.5 kcal/mol, respectively) and the formation of the H2C(OO)OH adduct is exothermic by about 74 kcal/mol. Beyond this addition process, further reaction mechanisms have also been investigated, which involve the abstraction of a hydrogen of carbonyl oxide by HO radical, but the computed activation barriers suggest that they will not contribute to the gas-phase reaction of H2COO + HO.     

Thermochemistry of the hypobromous and hypochlorous acids, HOBr and HOCl

The enthalpies of formation of HOBr and HOCl have been estimated by employing coupled cluster theory in conjunction with the correlation consistent basis sets and corrections for core-valence, relativistic, and anharmonic effects. We have employed three different reactions to estimate the DeltaH(o)(f,298)(HOBr), namely, the atomization reaction and two homodesmic reactions. Our best estimation is DeltaH(o)(f,298) (HOBr) = -15.3 +/- 0.6 kcal/mol and is very likely to lie toward the more negative values. The present value is 1.4 kcal/mol lower than the widely used experimental determination of Ruscic and Berkowitz (J. Chem. Phys. 1994, 101, 7795), DeltaH(o)(f,298)(HOBr) > -13.93 +/- 0.42 kcal/mol. However, it is closer to the more recent measurement of Lock et al. (J. Phys. Chem. 1996, 100, 7972), DeltaH(o)(f,298)(HOBr) = -14.8 +/- 1 kcal/mol. In the case of HOCl we have determined DeltaH(o)(f,298)(HOCl) = -18.1 +/- 0.3 kcal/mol, just in the middle of the two experimental values proposed, -17.8 +/- 0.5 kcal/mol (JANAF), obtained from equilibrium constant measurements, and -18.36 +/- 0.03 kcal/mol (Joens, J. A. J. Phys. Chem. A 2001, 105, 11041), determined from the measurements of the Cl-OH bond energy. If our conclusions are correct, several enthalpies of formation that have been determined by experimental chemists, Orlando and Burholder (J. Phys. Chem. 1995, 99, 1143), and theoretical chemists, Lee (J. Phys. Chem. 1995, 99, 15074), need to be revised, since a larger value was used for DeltaH(o)(f,298)(HOBr). Employing the results obtained by Orlando and Burkholder for Br(2)O we propose DeltaH(o)(f,298)(Br(2)O) = 24.9 +/- 0.6 kcal/mol, and employing Lee's enthalpies of reaction we propose the following DeltaH(o)(f,298): for BrBrO, HBrO, ClOBr, ClBrO, BrClO, BrCN, BrNC, BrNO, BrON, FOBr, and FBrO, 39.5 +/- 1, 41.0 +/- 1, 22.7 +/- 1.5, 34.2 +/- 1.5, 40.9 +/- 1.5, 43.7 +/- 1.5, 80.1 +/- 1.5, 22.3 +/- 1, 46.2 +/- 1, 17.3 +/- 1.5, and 6.3 +/- 1.5 kcal/mol, respectively. We expect that this work will stimulate new experimental measurements of the thermodynamic properties of HOBr and HOCl.     

Thermochemistry of the gaseous vanadium chlorides VCl, VCl2, VCl3, and VCl4

Gaseous equilibria in the V-Ag-Cl system were studied at elevated temperatures by effusion-beam mass spectrometry, where the pertinent species were generated by reaction of Cl 2(g) with V + Ag granules in the effusion cell source. Reaction enthalpies were derived from the equilibrium data, and the standard enthalpies of formation at 298 K of gaseous VCl, VCl2, and VCl3 were found to be +49.7, -34.8, and -85.6 kcal mol(-1), respectively. The corresponding bond dissociation energies at 298 K are D(V-Cl) = 102.9 kcal, D(ClV-Cl) = 113.5 kcal, D(Cl2V-Cl) = 79.8 kcal, and D(Cl3V-Cl) = 69.5 kcal. From these data, the dissociation energy D degrees 0(VCl) = 101.9 kcal mol(-1) or 4.42 eV is obtained. An alternate value, Delta(f)H(o)298(VCl 3,g) = -87.0 kcal mol (-1) was derived from third-law analysis of literature sublimation data for VCl3(s). In addition, literature thermochemical data on VCl4(g) were re-evaluated, leading to Delta(f)H(o)298 = -126.1 kcal mol (-1). The results are compared with various estimates in the literature.     

Structure and reactivity of aquacarbonylruthenium(II) complexes. An X-ray and oxygen-17 NMR study

The water exchange on [Ru(CO)(H2O-eq)4(H2O-ax)](tos)2 (1), [Ru(CO)2(H2O-eq)2(H2O-ax)2](tos)2 (2), and [Ru(CO)3(H2O)3](ClO4)2 (3), the 17O exchange between the bulk water and the carbonyl oxygens have been studied by 17O NMR spectroscopy, and the X-ray crystallographic structures of 1 and 2 have been determined. The water exchange of equatorially and axially coordinated water molecules on 1 and 2 follow an Id mechanism and are characterized by keq298 (s-1), delta H++ (kJ/mol), and delta S++ (J/(mol K)) of (2.54 +/- 0.05) x 10(-6), 111.6 +/- 0.4, and 22.4 +/- 1 (1-eq); (3.54 +/- 0.02) x 10(-2) and 81 (1-ax); (1.58 +/- 0.14) x 10(-7), 120.3 +/- 2, and 28.4 +/- 4 (2-eq); and (4.53 +/- 0.08) x 10(-4), 97.9 +/- 1, and 19.3 +/- 3 (2-ax). The observed reactivities correlate with the strength of the Ru-OH2 bonds, as expressed by their length obtained by X-ray studies: 2.079 (1-eq), 2.140 (1-ax), 2.073 (2-eq), and 2.110 (2-ax) A. 3 is strongly acidic witha pKa of -0.14 at 262 K. Therefore, the acid-dependent water exchange can take place through 3 or Ru(CO)3(H2O)3OH+ with an estimated keq298 of 10(-4)/10(-3) s-1 and kOH262 of 0.053 +/- 0.006 s-1. The 17O exchange rate between the bulk water and the carbonyl oxygens increases from 1 to 2 to 3. For 1 an upper limit of 10(-8) s-1 was estimated. For 2, no acid dependence of kRuCO between 0.1 and 1 m Htos was observed. At 312.6 K, in 0.1 and 1 m Htos, kRuCO = (1.18 +/- 0.03) x 10(-4). For the tricarbonyl complex, the exchange can proceed through 3 or Ru(CO)3(H2O)2OH+ with kRuCO and kRuOHCO of, respectively, 0.003 +/- 0.002 and 0.024 +/- 0.003 s-1, with a ruthenacarboxylic acid intermediate.     

Gas-phase acidities and O-H bond dissociation enthalpies of phenol, 3-methylphenol, 2,4,6-trimethylphenol, and ethanoic acid

Energy-resolved, competitive threshold collision-induced dissociation (TCID) methods are used to measure the gas-phase acidities of phenol, 3-methylphenol, 2,4,6-trimethylphenol, and ethanoic acid relative to hydrogen cyanide, hydrogen sulfide, and the hydroperoxyl radical using guided ion beam tandem mass spectrometry. The gas-phase acidities of Delta(acid)H298(C6H5OH) = 1456 +/- 4 kJ/mol, Delta(acid)H298(3-CH3C6H4OH) = 1457 +/- 5 kJ/mol, Delta(acid)H298(2,4,6-(CH3)3C6H2OH) = 1456 +/- 4 kJ/mol, and Delta(acid)H298(CH3COOH) = 1457 +/- 6 kJ/mol are determined. The O-H bond dissociation enthalpy of D298(C6H5O-H) = 361 +/- 4 kJ/mol is derived using the previously published experimental electron affinity for C6H5O, and thermochemical values for the other species are reported. A comparison of the new TCID values with both experimental and theoretical values from the literature is presented.