A restricted-open-shell complete-basis-set model chemistry

A restricted-open-shell model chemistry based on the complete basis set-quadratic Becke3 (CBS-QB3) model is formulated and denoted ROCBS-QB3. As the name implies, this method uses spin-restricted wave functions, both for the direct calculations of the various components of the electronic energy and for extrapolating the correlation energy to the complete-basis-set limit. These modifications eliminate the need for empirical corrections that are incorporated in standard CBS-QB3 to compensate for spin contamination when spin-unrestricted wave functions are used. We employ an initial test set of 19 severely spin-contaminated species including doublet radicals and both singlet and triplet biradicals. The mean absolute deviation (MAD) from experiment for the new ROCBS-QB3 model (3.6+/-1.5 kJ mol(-1)) is slightly smaller than that of the standard unrestricted CBS-QB3 version (4.8+/-1.5 kJ mol(-1)) and substantially smaller than the MAD for the unrestricted CBS-QB3 before inclusion of the spin correction (16.1+/-1.5 kJ mol(-1)). However, when applied to calculate the heats of formation at 298 K for the moderately spin-contaminated radicals in the G2/97 test set, ROCBS-QB3 does not perform quite as well as the standard unrestricted CBS-QB3, with a MAD from experiment of 3.8+/-1.6 kJ mol(-1) (compared with 2.9+/-1.6 kJ mol(-1) for standard CBS-QB3). ROCBS-QB3 performs marginally better than standard CBS-QB3 for the G2/97 set of ionization energies with a MAD of 4.1+/-0.1 kJ mol(-1) (compared with 4.4+/-0.1 kJ mol(-1)) and electron affinities with a MAD of 3.9+/-0.2 kJ mol(-1) (compared with 4.3+/-0.2 kJ mol(-1)), but the differences in MAD values are comparable to the experimental uncertainties. Our overall conclusion is that ROCBS-QB3 eliminates the spin correction in standard CBS-QB3 with no loss in accuracy.     

Conformational and thermodynamic properties of gaseous levulinic acid

Molecular modeling is used to determine low-energy conformational structures and thermodynamic properties of levulinic acid in the gas phase. Structure and IR vibrational frequencies are obtained using density functional and M?ller-Plesset perturbation theories. Electronic energies are computed using G3//B3LYP and CBS-QB3 model chemistries. Computed anharmonic frequencies are consistent with reported experimental data. Population analysis shows a boat- and a chainlike structure to be most abundant at 298 K, with increasing proportions of two other conformers at higher temperatures. Population mean distribution values for thermodynamic quantities are derived. At 298 K and 1 atm, the enthalpy of formation, entropy, and heat capacity are -613.1 ± 1.0 kJ·mol(-1), 407.4 J·mol(-1)·K(-1), and 132.3 J·mol(-1)·K(-1), respectively.     

Accurate prediction of thermodynamic properties of alkyl peroxides by combining density functional theory calculation with least-square calibration

Owing to the significance in kinetic modeling of the oxidation and combustion mechanisms of hydrocarbons, a fast and relatively accurate method was developed for the prediction of Delta(f)H(298)(o) of alkyl peroxides. By this method, a raw Delta(f)H(298)(o) value was calculated from the optimized geometry and vibration frequencies at B3LYP/6-31G(d,p) level and then an accurate Delta(f)H(298)(o) value was obtained by a least-square procedure. The least-square procedure is a six-parameter linear equation and is validated by a leave-one out technique, giving a cross-validation squared correlation coefficient q(2) of 0.97 and a squared correlation coefficient of 0.98 for the final model. Calculated results demonstrated that the least-square calibration leads to a remarkable reduction of error and to the accurate Delta(f)H(298)(o) values within the chemical accuracy of 8 kJ mol(-1) except (CH(3))(2)CHCH(2)CH(2)CH(2)OOH which has an error of 8.69 kJ mol(-1). Comparison of the results by CBS-Q, CBS-QB3, G2, and G3 revealed that B3LYP/6-31G(d,p) in combination with a least-square calibration is reliable in the accurate prediction of the standard enthalpies of formation for alkyl peroxides. Standard entropies at 298 K and heat capacities in the temperature range of 300-1500 K for alkyl peroxides were also calculated using the rigid rotor-harmonic oscillator approximation.     

The ozonolysis of acetylene--a quantum chemical investigation.

The ozonolysis of acetylene was investigated using CCSD(T), CASPT2, and B3LYP-DFT in connection with a 6-311+G(2d,2p) basis set. The reaction is initiated by the formation of a van der Waals complex followed by a [4pi + 2pi] cycloaddition between ozone and acetylene (activation enthalpy DeltaH(a)(298) = 9.6 kcal/mol; experiment, 10.2 kcal/mol), yielding 1,2,3-trioxolene, which rapidly opens to alpha-ketocarbonyl oxide 5. Alternatively, an O atom can be transferred from ozone to acetylene (DeltaH(a)(298) = 15.6 kcal/mol), thus leading to formyl carbene, which can rearrange to oxirene or ketene. The key compound in the ozonolysis of acetylene is 5 because it is the starting point for the isomerization to the corresponding dioxirane 19 (DeltaH(a)(298) = 16.9 kcal/mol), for the cyclization to trioxabicyclo[2.1.0]pentane 10 (DeltaH(a)(298) = 19.5 kcal/mol), for the formation of hydroperoxy ketene 15 (DeltaH(a)(298) = 20.6 kcal/mol), and for the rearrangement to dioxetanone 9 (DeltaH(a)(298) = 23.6 kcal/mol). Compounds 19, 10, 15, and 9 rearrange or decompose with barriers between 13 and 16 kcal/mol to yield as major products formanhydride, glyoxal, formaldehyde, formic acid, and (to a minor extent) glyoxylic acid. Hence, the ozonolysis of acetylene possesses a very complicated reaction mechanism that deserves intensive experimental studies.     

A photoelectron and TPEPICO investigation of the acetone radical cation

The valence shell photoelectron spectrum, threshold photoelectron spectrum, and threshold photoelectron photoion coincidence (TPEPICO) mass spectra of acetone have been measured using synchrotron radiation. New vibrational progressions have been observed and assigned in the X 2B2 state photoelectron bands of acetone-h6 and acetone-d6, and the influence of resonant autoionization on the threshold electron yield has been investigated. The dissociation thresholds for fragment ions up to 31 eV have been measured and compared to previous values. In addition, kinetic modeling of the threshold region for CH3* and CH4 loss leads to new values of 78 +/- 2 kJ mol(-1) and 75 +/- 2 kJ mol(-1), respectively, for the 0 K activation energies for these two processes. The result for the methyl loss channel is in reasonable agreement with, but slightly lower than, that of 83 +/- 1 kJ mol(-1) derived in a recent TPEPICO study by Fogleman et al. The modeling accounts for both low-energy dissociation channels at two different ion residence times in the mass spectrometer. Moreover, the effects of the ro-vibrational population distribution, the electron transmission efficiency, and the monochromator band-pass are included. The present activation energies yield a Delta(f)H298 for CH3CO+ of 655 +/- 3 kJ mol(-1), which is 4 kJ mol(-1) lower than that reported by Fogleman et al. The present Delta(f)H298 for CH3CO+ can be combined with the Delta(f)H298 for CH2CO (-47.5 +/- 1.6 kJ mol(-1)) and H+ (1530 kJ mol(-1)) to yield a 298 K proton affinity for ketene of 828 +/- 4 kJ mol(-1), in good agreement with the value (825 kJ mol(-1)) calculated at the G2 level of theory. The measured activation energy for CH4 loss leads to a Delta(f)H298 (CH2CO+*) of 873 +/- 3 kJ mol(-1).     

Theoretical determinations of the ambient conformational distribution and unimolecular decomposition of n-propylperoxy radical

The conformational distribution and unimolecular decomposition pathways for the n-propylperoxy radical have been generated at the CBS-QB3, B3LYP/6-31+G and mPW1K/6-31+G levels of theory. At each of the theoretical levels, the 298 K Boltzmann distributions and rotational profiles indicate that all five unique rotamers of the n-propylperoxy radical can be expected to be present in significant concentrations at thermal equilibrium. At the CBS-QB3 level, the 298 K distribution of rotamers is predicted to be 28.1, 26.4, 19.6, 14.0, and 11.9% for the gG, tG, gT, gG', and tT conformations, respectively. The CBS-QB3 C-OO bond dissociation energy (DeltaH298 K) for the n-propylperoxy radical has been calculated to be 36.1 kcal/mol. The detailed CBS-QB3 potential energy surface for the unimolecular decomposition of the n-propylperoxy radical indicates that important bimolecular products could be derived from two 1,4-H transfer mechanisms available at T < 500 K, primarily via an activated n-propylperoxy adduct.     

大气中一元酸催化亚硫酸分解反应的从头算及动力学研究

采用从头算的理论计算方法,在CCSD（T）/aug-cc-pVDZ//MP2/aug-cc-pVDZ水平下,对甲酸、乙酸、丙酸和硝酸催化亚硫酸分解成二氧化硫和水的微观反应机理进行了理论研究.结果表明,这4种一元酸均可使亚硫酸分解反应能垒显著降低,降低幅度由大到小的顺序为丙酸 > 乙酸 > 甲酸 > 硝酸.其中,丙酸甚至使反应能垒从裸反应时的99.84 kJ/mol降至27.24 kJ/mol.在此基础上,计算了200~320 K温度范围内4种一元酸催化亚硫酸分解反应的速率常数,并结合它们在大气中的实际浓度计算了有效速率常数.结果表明,在实际大气环境中,乙酸对亚硫酸分解反应的催化效果最为显著.当乙酸存在时,亚硫酸在室温下的大气寿命仅为0.02 s.     

Activation energies of pericyclic reactions: performance of DFT, MP2, and CBS-QB3 methods for the prediction of activation barriers and reaction energetics of 1,3-dipolar cycloadditions, and revised activation enthalpies for a standard set of hydrocarbon pericyclic reactions

Activation barriers and reaction energetics for the three main classes of 1,3-dipolar cycloadditions, including nine different reactions, were evaluated with the MPW1K and B3LYP density functional methods, MP2, and the multicomponent CBS-QB3 method. The CBS-QB3 values were used as standards for 1,3-dipolar cycloaddition activation barriers and reaction energetics, and the density functional theory (DFT) and MP2 methods were benchmarked against these values. The MPW1K/6-31G* method and basis set performs best for activation barriers, with a mean absolute deviation (MAD) value of 1.1 kcal/mol. The B3LYP/6-31G* method and basis set performs best for reaction enthalpies, with a MAD value of 2.4 kcal/mol, while the MPW1K method shows large errors for reaction energetics. The MP2 method gives the expected systematic underestimation of barriers. Concerted and nearly synchronous transition structures are predicted by all DFT and MP2 methods. Also reported are revised estimated 0 K experimental activation enthalpies for a standard set of hydrocarbon pericyclic reactions and updated comparisons to experiment for DFT, ab initio, and multicomponent methods. B3LYP and MPW1K methods with MAD values of 1.5 and 2.1 kcal/mol, respectively, fortuitously outperform the multicomponent CBS-QB3 method, which has a MAD value of 2.3. The MAD value of the O3LYP functional improves to 2.4 kcal/mol from the previously reported 3.0 kcal/mol.     

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.     

Nitrogen-14 NMR Study on Solvent Exchange of the Octahedral Cobalt(II) Ion in Neat 1,3-Propanediamine and n-Propylamine at Various Temperatures and Pressures. Tetrahedral-Octahedral Equilibrium of the Solvated Cobalt(II) Ion in n-Propylamine As Studied by EXAFS and Electronic Absorption Spectroscopy

Solvated cobalt(II) ions in neat 1,3-propanediamine (tn) and n-propylamine (pa) have been characterized by electronic absorption spectroscopy and extended X-ray absorption fine structure (EXAFS) spectroscopy. The equilibrium between tetrahedral and octahedral geometry for cobalt(II) ion has been observed in a neat pa solution, but not in neat diamine solutions such as tn and ethylenediamine (en). The thermodynamic parameters and equilibrium constant at 298 K for the geometrical equilibrium in pa were determined to be DeltaH degrees = -36.1 +/- 2.3 kJ mol(-1), DeltaS degrees = -163 +/- 8 J mol(-1) K(-1), and K(298) = 6.0 x 10(-3) M(-2), where K = [Co(pa)(6)(2+)]/{[Co(pa)(4)(2+)][pa](2)}. The equilibrium is caused by the large entropy gain in formation of the tetrahedral cobalt(II) species. The solvent exchange of cobalt(II) ion with octahedral geometry in tn and pa solutions has been studied by the (14)N NMR line-broadening method. The activation parameters and rate constants at 298 K for the solvent exchange reactions are as follows: DeltaH() = 49.3 +/- 0.9 kJ mol(-1), DeltaS() = 25 +/- 3 J mol(-1) K(-1), DeltaV() = 6.6 +/- 0.3 cm(3) mol(-1) at 302.1 K, and k(298) = 2.9 x 10(5) s(-1) for the tn exchange, and DeltaH() = 36.2 +/- 1.2 kJ mol(-1), DeltaS() = 35 +/- 6 J mol(-1) K(-1), and k(298) = 2.0 x 10(8) s(-1) for the pa exchange. By comparison of the activation parameters with those for the en exchange of cobalt(II) ion, it has been confirmed that the kinetic chelate strain effect is attributed to the large activation enthalpy for the bidentate chelate opening and that the enthalpic effect is smaller in the case of the six-membered tn chelate compared with the five-membered en chelate.     

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.     

Ab initio group contribution method for activation energies of hydrogen abstraction reactions.

The group contribution method for activation energies is applied to hydrogen abstraction reactions. To this end an ab initio database was constructed, which consisted of activation energies calculated with the ab initio CBS-QB3 method for a limited set of well-chosen homologous reactions. CBS-QB3 is shown to predict reaction rate coefficients within a factor of 2-4 and Arrhenius activation energies within 3-5 kJ mol(-1) of experimental data. Activation energies in the set of homologous reactions vary over 156 kJ mol(-1) with the structure of the abstracting radical and over 94 kJ mol(-1) with the structure of the abstracted hydrocarbon. The parameters required for the group contribution method, the so-called standard activation group additivity values, were determined from this database. To test the accuracy of the group contribution method, a large set of 88 additional activation energies were calculated from first principles and compared with the predictions from the group contribution method. It was found that the group contribution method yields accurate activation energies for hydrogen-transfer reactions between hydrogen molecules, alkylic hydrocarbons, and vinylic hydrocarbons, with the largest deviations being less than 6 kJ mol(-1). For reactions between allylic and propargylic hydrocarbons, the transition state is believed to be stabilized by resonance effects, thus requiring the introduction of an appropriate correction term to obtain a reliable prediction of the activation energy for this subclass of hydrogen abstraction reactions.     

ADSORPTION OF MACROPOROUS PHOSPHONIC ACID RESIN FOR INDIUM

The adsorption kinetics and mechanism of a nrovel chelate resin, macroporous phosphonic acid resin (PAR) for In(Ill) were investigated Tile statically saturated adsorption capacity is 216mg&#183;g^-1resin at 298K in HAc-NaAc medium. Tire apparent adsorption rate constant is k298=4.84&#215;10^-5 s^-1. Tile adsorption behavior of PAR for In(Ill) obeys the Freundlich isotherm. The thermodynamic adsorption parameters, enthalpy change △H, free energy change △G and entropychange △S of PAR for In(Ⅲ) are 11.5kmol, -12.6kJmol and 80.8Jmol.K, respectively. The apparent activation energy is Ea=3.5k.l/mol. Tire molar coordination ratio of the functional group of PAR to In(Ⅲ) is about 3:1.     

Thinking out of the black box: accurate barrier heights of 1,3-dipolar cycloadditions of ozone with acetylene and ethylene

Accurate barriers for the 1,3-dipolar cycloadditions of ozone with acetylene and ethylene have been determined via the systematic extrapolation of ab initio energies within the focal point approach of Allen and co-workers. Electron correlation has been accounted for primarily via coupled cluster theory, including single, double, and triple excitations, as well as a perturbative treatment of connected quadruple excitations [CCSD, CCSD(T), CCSDT, and CCSDT(Q)]. For the concerted [4 + 2] cycloadditions, the final recommended barriers are DeltaH(0K) = 9.4 +/- 0.2 and 5.3 +/- 0.2 kcal mol(-1) for ozone adding to acetylene and ethylene, respectively. These agree with recent results of Cremer et al. and Anglada et al., respectively. The reaction energy for O3 + C2H2 exhibits a protracted convergence with respect to inclusion of electron correlation, with the CCSDT/cc-pVDZ and CCSDT(Q)/cc-pVDZ values differing by 2.3 kcal mol-1. Recommended enthalpies of formation (298 K) for cycloadducts 1,2,3-trioxole and 1,2,3-trioxolane are +32.8 and -1.6 kcal mol(-1), respectively. Popular composite ab initio approaches [CBS-QB3, CBS-APNO, G3, G3B3, G3(MP2)B3, G4, G4(MP3), and G4(MP2)] predict a range of barrier heights for these systems. The CBS-QB3 computed barrier for ozone and acetylene, DeltaH(0K) = 4.4 kcal mol(-1), deviates by 5 kcal mol(-1) from the focal point value. CBS-QB3 similarly underestimates the barrier for the reaction of ozone and ethylene, yielding a prediction of only 0.7 kcal mol(-1). The errors in the CBS-QB3 results are significantly larger than mean errors observed in application to the G2 test set. The problem is traced to the nontransferability of MP2 basis set effects in the case of these reaction barriers. The recently published G4 and G4(MP2) approaches perform substantially better for O3 + C2H2, predicting enthalpy barriers of 9.0 and 8.4 kcal mol(-1), respectively. For the prediction of these reaction barriers, the additive corrections applied in the majority of the composite approaches considered lead to worse agreement with the reference focal point values than would be obtained relying only on single point energies evaluated at the highest level of theory utilized within each composite method.     

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.     

Thermochemical properties, DeltafH degrees (298), S degrees (298), and Cp degrees (T), for n-butyl and n-pentyl hydroperoxides and the alkyl and peroxy radicals, transition states, and kinetics for intramolecular hydrogen shift reactions of the peroxy radicals

Alkyl radicals in atmospheric and combustion environments undergo a rapid association with molecular oxygen (3O2) to form an alkyl peroxy radical (ROO*). One important reaction of these peroxy radicals is the intramolecular H-shift (intramolecular abstraction) to form a hydroperoxide alkyl radical (R'*COOH), where the hydroperoxide alkyl radical may undergo chemical activation reaction with O2 and result in chain branching at moderate to low temperatures. The thermochemistry and trends in kinetic parameters for the hydrogen shift reactions from each carbon (4-8-member-ring TST's) in n-butyl and n-pentyl peroxy radicals (CCCCOO* and CCCCCOO*) are analyzed using density functional and ab initio calculation methods. Thermochemical properties, DeltafH degrees (298 K), C-H bond energies, S degrees (298 K), and Cp degrees (T) of saturated linear C4 and C5 aliphatic peroxides (ROOH), as well as the corresponding hydroperoxide alkyl radicals (R'*COOH), are determined. DeltafH degrees (298 K) are obtained from isodesmic reactions and the total energies of the CBS-QB3 and B3LYP computational methods. Contributions to the entropy and the heat capacity from translation, vibration, and external rotation are calculated using the rigid-rotor-harmonic-oscillator approximation based on the CBS-QB3 frequencies and structures. The results indicate that pre-exponential factors, A(T), decrease with the increase of the ring size (4-8-member-ring TS, H-atom included). The DeltaH for 4-, 5-, 6-, and 7-member rings in n-butyl (and n-pentyl) peroxy are 40.8 (40.8), 31.4 (31.5), 20.5 (20.0), 22.6-p (19.4) kcal mol(-1), respectively. The DeltaH for the 8-member ring in n-pentylperoxy is 23.8-p kcal mol(-1), All abstractions are from secondary (-CH2-) groups except those marked (-p), which are from primary sites. Enthalpy and barrier values from the B3LYP/6-311++G(2d,p) and BHandHLYP/6-311G(d,p) methods are compared with CBS-QB3 results. The B3LYP results show good agreement with the higher level CBS-QB3 calculation method; the BHandH barriers for the intramolecular peroxy H-shifts are not acceptable.     

Mechanism of CO exchange on cis-[M(CO)2X2]- complexes (M=Rh,Ir;X=Cl, Br, or I)

The CO exchange on cis-[M(CO)2X2]- with M = Ir (X = Cl, la; X = Br, 1b; X = I, 1c) and M = Rh (X = Cl, 2a; X = Br, 2b; X = I, 2c) was studied in dichloromethane. The exchange reaction [cis-[M(CO)2X2]- + 2*CO is in equilibrium cis-[M(*CO)2X2]- + 2CO (exchange rate constant: kobs)] was followed as a function of temperature and carbon monoxide concentration (up to 6 MPa) using homemade high gas pressure NMR sapphire tubes. The reaction is first order for both CO and cis-[M(CO)2X2]- concentrations. The second-order rate constant, k2(298) (=kobs)[CO]), the enthalpy, deltaH*, and the entropy of activation, deltaS*, obtained for the six complexes are respectively as follows: la, (1.08 +/- 0.01) x 10(3) L mol(-1) s(-1), 15.37 +/- 0.3 kJ mol(-1), -135.3 +/- 1 J mol(-1) K(-1); 1b, (12.7 +/- 0.2) x 10(3) L mol(-1) s(-1), 13.26 +/- 0.5 kJ mol(-1), -121.9 +/- 2 J mol(-1) K(-1); 1c, (98.9 +/- 1.4) x 10(3) L mol(-1) s(-1), 12.50 +/- 0.6 kJ mol(-1), -107.4 +/- 2 J mol(-1) K(-1); 2a, (1.62 +/- 0.02) x 10(3) L mol(-1) s(-1), 17.47 +/- 0.4 kJ mol(-1), -124.9 +/- 1 J mol(-1) K(-1); 2b, (24.8 +/- 0.2) x 10(3) L mol(-1) s(-1), 11.35 +/- 0.4 kJ mol(-1), -122.7 +/- 1 J mol(-1) K(-1); 2c, (850 +/- 120) x 10(3) L mol(-1), s(-1), 9.87 +/- 0.8 kJ mol(-1), -98.3 +/- 4 J mol(-1) K(-1). For complexes la and 2a, the volumes of activation were measured and are -20.9 +/- 1.2 cm3 mol(-1) (332.0 K) and -17.2 +/- 1.0 cm3 mol(-1) (330.8 K), respectively. The second-order kinetics and the large negative values of the entropies and volumes of activation point to a limiting associative, A, exchange mechanism. The reactivity of CO exchange follows the increasing trans effect of the halogens (Cl < Br < I), and this is observed on both metal centers. For the same halogen, the rhodium complex is more reactive than the iridium complex. This reactivity difference between rhodium and iridium is less marked for chloride (1.5: 1) than for iodide (8.6:1) at 298 K.     

Theoretical study of the thermodynamics and kinetics of hydrogen abstractions from hydrocarbons

Thermochemical and kinetic data were calculated at four cost-effective levels of theory for a set consisting of five hydrogen abstraction reactions between hydrocarbons for which experimental data are available. The selection of a reliable, yet cost-effective method to study this type of reactions for a broad range of applications was done on the basis of comparison with experimental data or with results obtained from computationally demanding high level of theory calculations. For this benchmark study two composite methods (CBS-QB3 and G3B3) and two density functional theory (DFT) methods, MPW1PW91/6-311G(2d,d,p) and BMK/6-311G(2d,d,p), were selected. All four methods succeeded well in describing the thermochemical properties of the five studied hydrogen abstraction reactions. High-level Weizmann-1 (W1) calculations indicated that CBS-QB3 succeeds in predicting the most accurate reaction barrier for the hydrogen abstraction of methane by methyl but tends to underestimate the reaction barriers for reactions where spin contamination is observed in the transition state. Experimental rate coefficients were most accurately predicted with CBS-QB3. Therefore, CBS-QB3 was selected to investigate the influence of both the 1D hindered internal rotor treatment about the forming bond (1D-HR) and tunneling on the rate coefficients for a set of 21 hydrogen abstraction reactions. Three zero curvature tunneling (ZCT) methods were evaluated (Wigner, Skodje & Truhlar, Eckart). As the computationally more demanding centrifugal dominant small curvature semiclassical (CD-SCS) tunneling method did not yield significantly better agreement with experiment compared to the ZCT methods, CD-SCS tunneling contributions were only assessed for the hydrogen abstractions by methyl from methane and ethane. The best agreement with experimental rate coefficients was found when Eckart tunneling and 1D-HR corrections were applied. A mean deviation of a factor 6 on the rate coefficients is found for the complete set of 21 reactions at temperatures ranging from 298 to 1000 K. Tunneling corrections play a critical role in obtaining accurate rate coefficients, especially at lower temperatures, whereas the hindered rotor treatment only improves the agreement with experiment in the high-temperature range.     

Ab initio thermochemistry and kinetics for carbon-centered radical addition and beta-scission reactions

A quantitative comparison of ab initio calculated rate coefficients using five computational methods and five different approaches of treating hindered internal rotation and tunneling with experimental values of rate coefficients for nine carbon-centered radical additions/beta scissions at 300, 600, and 1000 K is performed. The high-accuracy compound methods, CBS-QB3 and G3B3, and the density functionals, MPW1PW91, BB1K, and BMK, have been evaluated using the following approaches: (i) the harmonic oscillator approximation; (ii) the hindered internal rotor approximation for the internal rotation about the forming/breaking bond in the transition state and product; and the hindered internal rotation approximation combined with (iii) Wigner, (iv) Skodje and Truhlar, and (v) Eckart zero-curvature tunneling corrections. The density functional theory (DFT) based values for beta-scission rate coefficients deviate significantly from the experimental ones at 300 K, and the DFT methods do not accurately predict the equilibrium coefficient. The hindered rotor approximation offers a significant improvement in the agreement with experimental rate coefficients as compared to the harmonic oscillator treatment, especially at higher temperatures. Tunneling correction factors are smaller than 1.40 at 300 K and 1.03 at 1000 K. For both the CBS-QB3 method, including the hindered rotor treatment but excluding tunneling corrections, and the G3B3 method, including hindered rotor and Eckart tunneling corrections, a mean factor of deviation with experimentally observed values of 3 is found.