Reaction of methacrolein with the hydroxyl radical in air: incorporation of secondary O2 addition into the MACR + OH master equation

Methacrolein (MACR) plays an important role in atmospheric chemistry within the planetary boundary layer, as it is one of the major oxidation products of isoprene and has a short lifetime toward the hydroxyl radical (OH). In this study, quantum chemical techniques and statistical reaction rate theory have been used to simulate the addition of OH to MACR at conditions representative of the troposphere. In this chemically activated reaction, the time scales for product formation versus collisional deactivation of the vibrationally excited adduct are explicitly considered. Furthermore, the subsequent addition of O(2) is also incorporated within a single master equation, so as to investigate doubly activated peroxyl radical formation. The major reaction product of OH addition to MACR is the HOCH(2)C(?)(CH(3))CHO radical formed via addition to the outer (β) carbon. This radical is predominantly in the Z isomer although around a third of the population is quenched as the higher-energy E isomer. Calculated rate constants agree well with experiment when using M06-2X/aug-cc-pVTZ barrier heights, but are somewhat overpredicted using G3SX energies. The overall rate constant is controlled by competition between dissociation of the MACR···OH van der Waals complex back to reactants and isomerization on to MACR-OH adducts, which takes place on a time scale of several nanoseconds, but collisional deactivation of the MACR-OH adducts occurs on a time scale that is around an order of magnitude longer. When O(2) addition is included in the master equation, we observe that the MACR-OH adducts are removed by reaction with O(2) on a similar time scale to collisional deactivation. Around 50% of the subsequent peroxyl radical population is formed with some identifiable excess vibrational energy above singly activated [MACR-OH-O(2)]*, with around 20% provided with an additional 20 kcal mol(-1) (>40 kcal mol(-1) relative to quenched MACR-OH-O(2)) that can go into further unimolecular reaction. This double activation process is expected to lead to some prompt unimolecular decomposition of excited [MACR-OH-O(2)]** peroxyl radicals to yield products including hydroxyacetone and methylglyoxal, regenerating the initiating OH radical in the process.     

Global potential energy surfaces for O((3)P) + H2O((1)A1) collisions

Global analytic potential energy surfaces for O((3)P) + H(2)O((1)A(1)) collisions, including the OH + OH hydrogen abstraction and H + OOH hydrogen elimination channels, are presented. Ab initio electronic structure calculations were performed at the CASSCF + MP2 level with an O(4s3p2d1f)/H(3s2p) one electron basis set. Approximately 10(5) geometries were used to fit the three lowest triplet adiabatic states corresponding to the triply degenerate O((3)P) + H(2)O((1)A(1)) reactants. Transition state theory rate constant and total cross section calculations using classical trajectories to collision energies up to 120?kcal mol(-1) (～11?km s(-1) collision velocity) were performed and show good agreement with experimental data. Flux-velocity contour maps are presented at selected energies for H(2)O collisional excitation, OH + OH, and H + OOH channels to further investigate the dynamics, especially the competition and distinct dynamics of the two reactive channels. There are large differences in the contributions of each of the triplet surfaces to the reactive channels, especially at higher energies. The present surfaces should support quantitative modeling of O((3)P) + H(2)O((1)A(1)) collision processes up to ～150?kcal mol(-1).     

Thermochemical and kinetic analysis on the reactions of O2 with products from OH addition to isobutene, 2-hydroxy-1,1-dimethylethyl, and 2-hydroxy-2-methylpropyl radicals: HO2 formation from oxidation of neopentane, Part II

Unimolecular dissociation of a neopentyl radical to isobutene and methyl radical is competitive with the neopentyl association with O2 ((3)Sigma(g)-) in thermal oxidative systems. Furthermore, both isobutene and the OH radical are important primary products from the reactions of neopentyl with O2. Consequently, the reactions of O2 with the 2-hydroxy-1,1-dimethylethyl and 2-hydroxy-2-methylpropyl radicals resulting from the OH addition to isobutene are important to understanding the oxidation of neopentane and other branched hydrocarbons. Reactions that correspond to the association of radical adducts with O2((3)Sigma(g)-) involve chemically activated peroxy intermediates, which can isomerize and react to form one of several products before stabilization. The above reaction systems were analyzed with ab initio and density functional calculations to evaluate the thermochemistry, reaction paths, and kinetics that are important in neopentyl radical oxidation. The stationary points of potential energy surfaces were analyzed based on the enthalpies calculated at the CBS-Q level. The entropies, S(degrees)298, and heat capacities, C(p)(T), (0 相似文献   

Mechanism of the hydration of carbon dioxide: direct participation of H2O versus microsolvation

Thermochemical parameters of carbonic acid and the stationary points on the neutral hydration pathways of carbon dioxide, CO 2 + nH 2O --> H 2CO 3 + ( n - 1)H 2O, with n = 1, 2, 3, and 4, were calculated using geometries optimized at the MP2/aug-cc-pVTZ level. Coupled-cluster theory (CCSD(T)) energies were extrapolated to the complete basis set limit in most cases and then used to evaluate heats of formation. A high energy barrier of approximately 50 kcal/mol was predicted for the addition of one water molecule to CO 2 ( n = 1). This barrier is lowered in cyclic H-bonded systems of CO 2 with water dimer and water trimer in which preassociation complexes are formed with binding energies of approximately 7 and 15 kcal/mol, respectively. For n = 2, a trimeric six-member cyclic transition state has an energy barrier of approximately 33 (gas phase) and a free energy barrier of approximately 31 (in a continuum solvent model of water at 298 K) kcal/mol, relative to the precomplex. For n = 3, two reactive pathways are possible with the first having all three water molecules involved in hydrogen transfer via an eight-member cycle, and in the second, the third water molecule is not directly involved in the hydrogen transfer but solvates the n = 2 transition state. In the gas phase, the two transition states have comparable energies of approximately 15 kcal/mol relative to separated reactants. The first path is favored over in aqueous solution by approximately 5 kcal/mol in free energy due to the formation of a structure resembling a (HCO 3 (-)/H 3OH 2O (+)) ion pair. Bulk solvation reduces the free energy barrier of the first path by approximately 10 kcal/mol for a free energy barrier of approximately 22 kcal/mol for the (CO 2 + 3H 2O) aq reaction. For n = 4, the transition state, in which a three-water chain takes part in the hydrogen transfer while the fourth water microsolvates the cluster, is energetically more favored than transition states incorporating two or four active water molecules. An energy barrier of approximately 20 (gas phase) and a free energy barrier of approximately 19 (in water) kcal/mol were derived for the CO 2 + 4H 2O reaction, and again formation of an ion pair is important. The calculated results confirm the crucial role of direct participation of three water molecules ( n = 3) in the eight-member cyclic TS for the CO 2 hydration reaction. Carbonic acid and its water complexes are consistently higher in energy (by approximately 6-7 kcal/mol) than the corresponding CO 2 complexes and can undergo more facile water-assisted dehydration processes.     

Model of uptake of OH radicals on nonreactive solids

A model of adsorption and recombination of OH radicals was developed for nonreactive solid surfaces of atmospheric interest. A parametrization of this heterogeneous mechanism was carried out to determine the role of the catalytic properties of these solid surfaces, taking into account the adsorption energy, defects, surface diffusion, and chemical reactions in the gas-solid interface. The uptake process was simulated for diffusion-controlled chemical reactions on the surface on the basis of Langmuir-Hinshelwood and Eley-Rideal mechanisms. Using an analytical approach and the Monte Carlo technique, we show the dependencies of the uptake probability of the heterogeneous reactions on the OH concentration and adsorption energy. The model is employed in the analysis of the empirically derived uptake coefficient for water ice, Al(2)O(3), NaCl, NH(4)NO(3), NH(4)HSO(4), and (NH(4))(2)SO(4). We found the following values for the free energy of adsorption of OH radicals: E(ice) = 7.3-7.6 kcal/mol, E(Al)(2)(O)(3) = 11-11.7 kcal/mol, E(NH)(4)(NO)(3) = 10.2 kcal/mol, E(NaCl) = 10.2 kcal/mol, E(NH)(4)(HSO)(4) = 9.8 kcal/mol, and E((NH)(4))(2)(SO)(4) = 9.8 kcal/mol. The atmospheric implications of the catalytic reactions of OH with adsorbed reactive molecules are discussed. The results of the modeling of the uptake process showed that the heterogeneous decay rate can exceed the corresponding gas-phase reaction rate under atmospheric conditions.     

Formation of a Criegee intermediate in the low-temperature oxidation of dimethyl sulfoxide

Dimethyl sulfoxide (DMSO) is the major sulfur-containing constituent of the Marine Boundary Layer. It is a significant source of H2SO4 aerosol/particles and methane sulfonic acid via atmospheric oxidation processes, where the mechanism is not established. In this study, several new, low-temperature pathways are revealed in the oxidation of DMSO using CBS-QB3 and G3MP2 multilevel and B3LYP hybrid density functional quantum chemical methods. Unlike analogous hydrocarbon peroxy radicals the chemically activated DMSO peroxy radical, [CH3S(=O)CH2OO*]*, predominantly undergoes simple dissociation to a methylsulfinyl radical CH3S*(=O) and a Criegee intermediate, CH2OO, with the barrier to dissociation 11.3 kcal mol(-1) below the energy of the CH3S(=O)CH2* + O2 reactants. The well depth for addition of O2 to the CH3S(=O)CH2 precursor radical is 29.6 kcal mol(-1) at the CBS-QB3 level of theory. We believe that this reaction may serve an important role in atmospheric photochemical and irradiated biological (oxygen-rich) media where formation of initial radicals is facilitated even at lower temperatures. The Criegee intermediate (carbonyl oxide, peroxymethylene) and sulfinyl radical can further decompose, resulting in additional chain branching. A second reaction channel important for oxidation processes includes formation (via intramolecular H atom transfer) and further decomposition of hydroperoxide methylsulfoxide radical, *CH2S(=O)CH2OOH over a low barrier of activation. The initial H-transfer reaction is similar and common in analogous hydrocarbon radical + O2 reactions; but the subsequent very low (3-6 kcal mol(-1)) barrier (14 kcal mol(-1) below the initial reagents) to beta-scission products is not common in HC systems. The low energy reaction of the hydroperoxide radical is a beta-scission elimination of *CH2S(=O)CH2OOH into the CH2=S=O + CH2O + *OH product set. This beta-scission barrier is low, because of the delocalization of the *CH2 radical center through the -S(=O) group, to the -CH2OOH fragment in the transition state structure. The hydroperoxide methylsulfoxide radical can also decompose via a second reaction channel of intramolecular OH migration, yielding formaldehyde and a sulfur-centered hydroxymethylsulfinyl radical HOCH2S*(=O). The barrier of activation relative to initial reagents is 4.2 kcal mol(-1). Heats of formation for DMSO, DMSO carbon-centered radical and Criegee intermediate are evaluated at 298 K as -35.97 +/- 0.05, 13.0 +/- 0.2 and 25.3 +/- 0.7 kcal mol(-1) respectively using isodesmic reaction analysis. The [CH3S*(=O) + CH2OO] product set is shown to form a van der Waals complex that results in O-atom transfer reaction and the formation of new products CH3SO2* radical and CH2O. Proper orientation of the Criegee intermediate and methylsulfinyl radical, as a pre-stabilized pre-reaction complex, assist the process. The DMSO radical reaction is also compared to that of acetonyl radical.     

Mechanisms of the reactions of W AND W+ with H2O: computational studies

The mechanism of the reactions of W and W(+) with the water molecule have been studied for several lower-lying electronic states of tungsten centers at the CCSD(T)/6-311G(d,p)+SDD and B3LYP/6-31G(d,p)+SDD levels of theory. It is shown that these reactions are essentially multistate processes, during which lower-lying electronic states of the systems cross several times. They start with the formation of initial prereaction M(H(2)O) complexes with M-H(2)O bonding energies of 9.6 and 48.2 kcal/mol for M = W and W(+), followed by insertion of the metal center into an O-H bond with 20.0 and 53.3 kcal/mol barriers for neutral and cationic systems, respectively. The overall process of M + H(2)O --> t-HM(OH) is calculated to be highly exothermic, 48.4 and 48.8 kcal/mol for M = W and W(+). From the HM(OH) intermediate the reaction may proceed via several different channels, among which the stepwise HM(OH) --> HMO + H --> (H)(2)MO and concerted HM(OH) --> (H)(2)MO pathways are more favorable and can compete (energetically) with each other. For the neutral system (M = W), the concerted process is the most favorable, whereas for the charged system (M = W(+)), the stepwise pathway is slightly more favorable. From the energetically most favorable intermediate (H)(2)MO the reactions proceed via H(2)-molecule formation with a 53.1 kcal/mol activation barrier for the neutral system. For the cationic system, H-H formation and dissociation is an almost barrierless process. The overall reaction of W and W(+) with the water molecule leading to H(2) + MO formation is found to be exothermic by 48.2 and 39.8 kcal/mol, respectively. In the gas phase with the collision-less conditions the reactions W((7)S) + H(2)O --> H(2) + WO((3)Sigma(+)), and W(+)((6)D) + H(2)O --> H(2) + WO(+)((4)Sigma(+)) are expected to proceed via a 10.4 and 5.1 kcal/mol overall energy barrier corresponding to the first O-H dissociation at the TS1. On the basis of these PESs, we predict kinetic rate constants for the reactions of W and W(+) with H(2)O.     

Reaction mechanism of naphthyl radicals with molecular oxygen. 1. Theoretical study of the potential energy surface

Potential energy surfaces (PESs) of the reactions of 1- and 2-naphthyl radicals with molecular oxygen have been investigated at the G3(MP2,CC)//B3LYP/6-311G** level of theory. Both reactions are shown to be initiated by barrierless addition of O(2) to the respective radical sites of C(10)H(7). The end-on O(2) addition leading to 1- and 2-naphthylperoxy radicals exothermic by 45-46 kcal/mol is found to be more preferable thermodynamically than the side-on addition. At the subsequent reaction step, the chemically activated 1- and 2-C(10)H(7)OO adducts can eliminate an oxygen atom leading to the formation of 1- and 2-naphthoxy radical products, respectively, which in turn can undergo unimolecular decomposition producing indenyl radical + CO via the barriers of 57.8 and 48.3 kcal/mol and with total reaction endothermicities of 14.5 and 10.2 kcal/mol, respectively. Alternatively, the initial reaction adducts can feature an oxygen atom insertion into the attacked C(6) ring leading to bicyclic intermediates a10 and a10' (from 1-naphthyl + O(2)) or b10 and b10' (from 2-naphthyl + O(2)) composed from two fused six-member C(6) and seven-member C(6)O rings. Next, a10 and a10' are predicted to decompose to C(9)H(7) (indenyl) + CO(2), 1,2-C(10)H(6)O(2) (1,2-naphthoquinone) + H, and 1-C(9)H(7)O (1-benzopyranyl) + CO, whereas b10 and b10' would dissociate to C(9)H(7) (indenyl) + CO(2), 2-C(9)H(7)O (2-benzopyranyl) + CO, and 1,2-C(10)H(6)O(2) (1,2-naphthoquinone) + H. On the basis of this, the 1-naphthyl + O(2) reaction is concluded to form the following products (with the overall reaction energies given in parentheses): 1-naphthoxy + O (-15.5 kcal/mol), indenyl + CO(2) (-123.9 kcal/mol), 1-benzopyranyl + CO (-97.2 kcal/mol), and 1,2-naphthoquinone + H (-63.5 kcal/mol). The 2-naphthyl + O(2) reaction is predicted to produce 2-naphthoxy + O (-10.9 kcal/mol), indenyl + CO(2) (-123.7 kcal/mol), 2-benzopyranyl + CO (-90.7 kcal/mol), and 1,2-naphthoquinone + H (-63.2 kcal/mol). Simplified kinetic calculations using transition-state theory computed rate constants at the high-pressure limit indicate that the C(10)H(7)O + O product channels are favored at high temperatures, while the irreversible oxygen atom insertion first leading to the a10 and a10' or b10 and b10' intermediates and then to their various decomposition products is preferable at lower temperatures. Among the decomposition products, indenyl + CO(2) are always most favorable at lower temperatures, but the others, 1,2-C(10)H(6)O(2) (1,2-naphthoquinone) + H (from a10 and b10'), 1-C(9)H(7)O (1-benzopyranyl) + CO (from a10'), and 2-C(10)H(7)O (2-benzopyranyl) + O (from b10 and minor from b10'), may notably contribute or even become major products at higher temperatures.     

On the mechanism of (PCP)Ir-catalyzed acceptorless dehydrogenation of alkanes: a combined computational and experimental study

Pincer complexes of the type ((R)PCP)IrH(2), where ((R)PCP)Ir is [eta(3)-2,6-(R(2)PCH(2))(2)C(6)H(3)]Ir, are the most effective catalysts reported to date for the "acceptorless" dehydrogenation of alkanes to yield alkenes and free H(2). We calculate (DFT/B3LYP) that associative (A) reactions of ((Me)PCP)IrH(2) with model linear (propane, n-PrH) and cyclic (cyclohexane, CyH) alkanes may proceed via classical Ir(V) and nonclassical Ir(III)(eta(2)-H(2)) intermediates. A dissociative (D) pathway proceeds via initial loss of H(2), followed by C-H addition to ((Me)PCP)Ir. Although a slightly higher energy barrier (DeltaE(+ +)) is computed for the D pathway, the calculated free-energy barrier (DeltaG(+ +)) for the D pathway is significantly lower than that of the A pathway. Under standard thermodynamic conditions (STP), C-H addition via the D pathway has DeltaG(o)(+ +) = 36.3 kcal/mol for CyH (35.1 kcal/mol for n-PrH). However, acceptorless dehydrogenation of alkanes is thermodynamically impossible at STP. At conditions under which acceptorless dehydrogenation is thermodynamically possible (for example, T = 150 degrees C and P(H)2 = 1.0 x 10(-7) atm), DeltaG(+ +) for C-H addition to ((Me)PCP)Ir (plus a molecule of free H(2)) is very low (17.5 kcal/mol for CyH, 16.7 kcal/mol for n-PrH). Under these conditions, the rate-determining step for the D pathway is the loss of H(2) from ((Me)PCP)IrH(2) with DeltaG(D)(+ +) approximately DeltaH(D)(+ +) = 27.2 kcal/mol. For CyH, the calculated DeltaG(o)(+ +) for C-H addition to ((Me)PCP)IrH(2) on the A pathway is 35.2 kcal/mol (32.7 kcal/mol for n-PrH). At catalytic conditions, the calculated free energies of C-H addition are 31.3 and 33.7 kcal/mol for CyH and n-PrH addition, respectively. Elimination of H(2) from the resulting "seven-coordinate" Ir-species must proceed with an activation enthalpy at least as large as the enthalpy change of the elimination step itself (DeltaH approximately 11-13 kcal/mol), and with a small entropy of activation. The free energy of activation for H(2) elimination (DeltaG(A)(+ +)) is hence found to be greater than ca. 36 kcal/mol for both CyH and n-PrH under catalytic conditions. The overall free-energy barrier of the A pathway is calculated to be higher than that of the D pathway by ca. 9 kcal/mol. Reversible C-H(D) addition to ((R)PCP)IrH(2) is predicted to lead to H/D exchange, because the barriers for hydride scrambling are extremely low in the "seven-coordinate" polyhydrides. In agreement with calculation, H/D exchange is observed experimentally for several deuteriohydrocarbons with the following order of rates: C(6)D(6) > mesitylene-d(12) > n-decane-d(22) > cyclohexane-d(12). Because H/D exchange in cyclohexane-d(12) solution is not observed even after 1 week at 180 degrees C, we estimate that the experimental barrier to cyclohexane C-D addition is greater than 36.4 kcal/mol. This value is considerably greater than the experimental barrier for the full catalytic dehydrogenation cycle for cycloalkanes (ca. 31 kcal/mol). Thus, the experimental evidence, in agreement with calculation, strongly indicates that the A pathway is not kinetically viable as a segment of the "acceptorless" dehydrogenation cycle.     

Ab initio investigation of the electronic and geometric structure of magnesium diboride, MgB2

Employing multireference variational (MRCI) and coupled cluster (CC) methods combined with quadruple-zeta quality correlation-consistent basis set, we have studied 36 states of the magnesium diboride (MgB(2)) molecule as well as 17 states of the experimentally unknown diatomic MgB. For all states of MgB(2), we report geometries, atomization energies, and dipole moments, while for the first 5 states, potential energy profiles have been also constructed. The ground state is formally of (1)A(1) V-shaped symmetry with an atomization energy of 108.1(109) kcal/mol at the MRCI(MRCI + Davidson correction) level. The first excited state ((3)B(1)) is less than 1 kcal/mol above the X(1)A(1) state, with the next state of linear Mg-B-B geometry (b(3)Sigma(-)) located 10 kcal/mol higher. In all states, bent or linear, the bonding is complicated and unconventional because of the extraordinary bonding agility of the boron atom(s).     

Ab initio study on the decomposition of first excited state HOOO radicals

The HOOO radical plays a crucial role in atmospheric processes involving the OH radical and O(2) molecule. We present an ab initio molecular orbital theory study on the decomposition reaction of the first excited state HOOO((2)A') with respect to OH and O(2). The geometries and harmonic vibrational frequencies of all stationary points are calculated at the CASSCF and MRCI levels of theory in conjunction with the 6-31+G(d,p) basis set. The potential energy profile of the decomposition reaction is studied at the CASSCF/6-31+G(d,p) level of theory, in which the complete valence orbitals and electrons are included in the active space. The energies of the potential energy profile are further refined at the CASPT2 and MRCI levels of the theory. Additionally, we have determined the interesting reaction process: the HOOO((2)A') radical with C(s) symmetry does not dissociate to OH((2)Pi) and O(2)((3)Sigma(-)(g)) directly as this is forbidden by orbital symmetry, but dissociates to OH((2)Pi) and O(2)((3)Sigma(-)(g)) via the change in symmetry from C(s) to C(infinity v) symmetry with a low barrier.     

Reaction of the i-C4H5 (CH2CCHCH2) radical with O2

The resonantly stabilized radical i-C(4)H(5) (CH(2)CCHCH(2)) is an important intermediate in the combustion of unsaturated hydrocarbons and is thought to be involved in the formation of polycyclic aromatic hydrocarbons through its reaction with acetylene (C(2)H(2)) to form benzene + H. This study uses quantum chemistry and statistical reaction rate theory to investigate the mechanism and kinetics of the i-C(4)H(5) + O(2) reaction as a function of temperature and pressure, and unlike most resonantly stabilized radicals we show that i-C(4)H(5) is consumed relatively rapidly by its reaction with molecular oxygen. O(2) addition occurs at the vinylic and allenic radical sites in i-C(4)H(5), with respective barriers of 0.9 and 4.9 kcal mol(-1). Addition to the allenic radical form produces an allenemethylperoxy radical adduct with only around 20 kcal mol(-1) excess vibrational energy. This adduct can isomerize to the ca. 14 kcal mol(-1) more stable 1,3-divinyl-2-peroxy radical via concerted and stepwise processes, both steps with barriers around 10 kcal mol(-1) below the entrance channel energy. Addition of O(2) to the vinylic radical site in i-C(4)H(5) directly forms the 1,3-divinyl-2-peroxy radical with a small barrier and around 36.8 kcal mol(-1) of excess energy. The 1,3-divinyl-2-peroxy radical isomerizes via ipso addition of the O(2) moiety followed by O atom insertion into the adjacent C-C bond. This process forms an unstable intermediate that ultimately dissociates to give the vinyl radical, formaldehyde, and CO. At higher temperatures formation of vinylacetylene + HO(2), the vinoxyl radical + ketene, and the 1,3-divinyl-2-oxyl radical + O paths have some importance. Because of the adiabatic transition states for O(2) addition, and significant reverse dissociation channels in the peroxy radical adducts, the i-C(4)H(5) + O(2) reaction proceeds to new products with rate constant of around 10(11) cm(3) mol(-1) s(-1) at typical combustion temperatures (1000-2000 K). For fuel-rich flames we show that the reaction of i-C(4)H(5) with O(2) is likely to be faster than that with C(2)H(2), bringing into question the importance of the i-C(4)H(5) + C(2)H(2) reaction in initiating ring formation in sooting flames.     

Comparative DFT study of the spin trapping of methyl, mercapto, hydroperoxy, superoxide, and nitric oxide radicals by various substituted cyclic nitrones

The thermodynamics of the spin trapping of various cyclic nitrones with biologically relevant radicals such as methyl, mercapto, hydroperoxy, superoxide anion, and nitric oxide was investigated using computational methods. A density functional theory (DFT) approach was employed in this study at the B3LYP/6-31+G(d,p)//B3LYP/6-31G(d) level. The order of increasing favorability for Delta G(rxn) (kcal/mol) of the radical reaction with various nitrones, in general, follows a trend similar to their respective experimental reduction potentials as well as their experimental second-order rate constants in aqueous solution: NO (14.57) < O2*- (-7.51) < *O2H (-13.92) < *SH (-16.55) < *CH3 (-32.17) < *OH (-43.66). The same qualitative trend is predicted upon considering the effect of solvation using the polarizable continuum model (PCM): i.e., NO (14.12) < O2*- (9.95) < *O2H (-6.95) < *SH (-13.57) < *CH3 (-32.88) < *OH (-38.91). All radical reactions with these nitrones are exoergic, except for NO (and O2*- in the aqueous phase), which is endoergic, and the free energy of activation (Delta G) for the NO additions ranges from 17.7 to 20.3 kcal/mol. This study also predicts the favorable formation of certain adducts that exhibit intramolecular H-bonding interactions, nucleophilic addition, or H-atom transfer reactions. The spin density on the nitronyl N of the superoxide adducts reveals conformational dependences. The failure of nitrones to trap NO at normal conditions was theoretically rationalized due to the endoergic reaction parameters.     

Characterization of singlet ground and low-lying electronic excited states of phosphaethyne and isophosphaethyne

The singlet ground ((approximate)X(1)Sigma1+) and excited (1Sigma-,1Delta) states of HCP and HPC have been systematically investigated using ab initio molecular electronic structure theory. For the ground state, geometries of the two linear stationary points have been optimized and physical properties have been predicted utilizing restricted self-consistent field theory, coupled cluster theory with single and double excitations (CCSD), CCSD with perturbative triple corrections [CCSD(T)], and CCSD with partial iterative triple excitations (CCSDT-3 and CC3). Physical properties computed for the global minimum ((approximate)X(1)Sigma+HCP) include harmonic vibrational frequencies with the cc-pV5Z CCSD(T) method of omega1=3344 cm(-1), omega2=689 cm(-1), and omega3=1298 cm(-1). Linear HPC, a stationary point of Hessian index 2, is predicted to lie 75.2 kcal mol(-1) above the global minimum HCP. The dissociation energy D0[HCP((approximate)X(1)Sigma+)-->H(2S)+CP(X2Sigma+)] of HCP is predicted to be 119.0 kcal mol(-1), which is very close to the experimental lower limit of 119.1 kcal mol(-1). Eight singlet excited states were examined and their physical properties were determined employing three equation-of-motion coupled cluster methods (EOM-CCSD, EOM-CCSDT-3, and EOM-CC3). Four stationary points were located on the lowest-lying excited state potential energy surface, 1Sigma- -->1A", with excitation energies Te of 101.4 kcal mol(-1) (1A"HCP), 104.6 kcal mol(-1)(1Sigma-HCP), 122.3 kcal mol(-1)(1A" HPC), and 171.6 kcal mol(-1)(1Sigma-HPC) at the cc-pVQZ EOM-CCSDT-3 level of theory. The physical properties of the 1A" state with a predicted bond angle of 129.5 degrees compare well with the experimentally reported first singlet state ((approximate)A1A"). The excitation energy predicted for this excitation is T0=99.4 kcal mol(-1) (34 800 cm(-1),4.31 eV), in essentially perfect agreement with the experimental value of T0=99.3 kcal mol(-1)(34 746 cm(-1),4.308 eV). For the second lowest-lying excited singlet surface, 1Delta-->1A', four stationary points were found with Te values of 111.2 kcal mol(-1) (2(1)A' HCP), 112.4 kcal mol(-1) (1Delta HPC), 125.6 kcal mol(-1)(2(1)A' HCP), and 177.8 kcal mol(-1)(1Delta HPC). The predicted CP bond length and frequencies of the 2(1)A' state with a bond angle of 89.8 degrees (1.707 A, 666 and 979 cm(-1)) compare reasonably well with those for the experimentally reported (approximate)C(1)A' state (1.69 A, 615 and 969 cm(-1)). However, the excitation energy and bond angle do not agree well: theoretical values of 108.7 kcal mol(-1) and 89.8 degrees versus experimental values of 115.1 kcal mol(-1) and 113 degrees. of 115.1 kcal mol(-1) and 113 degrees.     

Energetics and mechanism of the decomposition of trifluoromethanol

The thermal instability of alpha-fluoroalcohols is generally attributed to a unimolecular 1,2-elimination of HF, but the barrier to intramolecular HF elimination from CF3OH is predicted to be 45.1 +/- 2 kcal/mol. The thermochemical parameters of trifluoromethanol were calculated using coupled-cluster theory (CCSD(T)) extrapolated to the complete basis set limit. High barriers of 42.9, 43.1, and 45.0 kcal/mol were predicted for the unimolecular decompositions of CH2FOH, CHF2OH, and CF3OH, respectively. These barriers are lowered substantially if cyclic H-bonded dimers of CF3OH with complexation energies of approximately 5 kcal/mol are involved. A six-membered ring dimer has an energy barrier of 28.7 kcal/mol and an eight-membered dimer has an energy barrier of 32.9 kcal/mol. Complexes of CF3OH with HF lead to strong H-bonded dimers, trimers and tetramers with complexation energies of approximately 6, 11, and 16 kcal/mol, respectively. The dimer, CH3OH:HF, and the trimers, CF3OH:2HF and (CH3OH)2:HF, have decomposition energy barriers of 26.7, 20.3, and 22.8 kcal/mol, respectively. The tetramer (CH3OH:HF)2 gives rise to elimination of two HF molecules with a barrier of 32.5 kcal/mol. Either CF3OH or HF can act as catalysts for HF-elimination via an H-transfer relay. Because HF is one of the decomposition products, the decomposition reactions become autocatalytic. If the energies due to complexation for the CF3OH-HF adducts are not dissipated, the effective barriers to HF elimination are lowered from approximately 20 to approximately 9 kcal/mol, which reconciles the computational results with the experimentally observed stabilities.     

Kinetics of the OH + ClOOCl and OH + Cl2O reactions: experiment and theory

The rate coefficients for the reactions OH + ClOOCl --> HOCl + ClOO (eq 5) and OH + Cl2O --> HOCl + ClO (eq 6) were measured using a fast flow reactor coupled with molecular beam quadrupole mass spectrometry. OH was detected using resonance fluorescence at 309 nm. The measured Arrhenius expressions for these reactions are k5 = (6.0 +/- 3.5) x 10(-13) exp((670 +/- 230)/T) cm(3) molecule(-1) s(-1) and k6 = (5.1 +/- 1.5) x 10(-12) exp((100 +/- 92)/T) cm(3) molecule(-1) s(-1), respectively, where the uncertainties are reported at the 2sigma level. Investigation of the OH + ClOOCl potential energy surface using high level ab initio calculations indicates that the reaction occurs via a chlorine abstraction from ClOOCl by the OH radical. The lowest energy pathway is calculated to proceed through a weak ClOOCl-OH prereactive complex that is bound by 2.6 kcal mol(-1) and leads to ClOO and HOCl products. The transition state to product formation is calculated to be 0.59 kcal mol(-1) above the reactant energy level. Inclusion of the OH + ClOOCl rate data into an atmospheric model indicates that this reaction contributes more than 15% to ClOOCl loss during twilight conditions in the Arctic stratosphere. Reducing the rate of ClOOCl photolysis, as indicated by a recent re-examination of the ClOOCl UV absorption spectrum, increases the contribution of the OH + ClOOCl reaction to polar stratospheric loss of ClOOCl.     

Dissociative addition of water to a main group 13 metal cluster: computational study of the reaction of gallium dimer with H2O

We have investigated the lowest triplet and singlet potential energy surfaces (PESs) for the reaction of Ga(2) dimer with water. Under thermal conditions, we predict formation of the triplet ground state addition complex Ga(2)···OH(2)((3)B(1)) involving Ga···O···Ga bridge interaction. At the coupled cluster CCSD(T)/AE (CCSD(T)/ECP) computational levels, Ga(2)···OH(2)((3)B(1)) is bound by 5.5 (5.7) kcal/mol with respect to the ground state reactants Ga(2)((3)Π(u)) + H(2)O. Identification of the addition complex is in agreement with the experimental evidence from matrix isolation infrared (IR) spectroscopy reported recently by Macrae and Downs. The located minimum energy crossing points (MECPs) between the triplet and singlet energy surfaces on the entrance channel of Ga(2) + H(2)O are not expected to be energetically accessible under the matrix conditions, consistent with the lack of occurrence of Ga(2) insertion into the O-H bond under such conditions. The computed energies and harmonic and anharmonic vibrational frequencies for the triplet and singlet Ga(2)(H)(OH) insertion isomers indicate the singlet double-bridged Ga(μ-H)(μ-OH)Ga isomer to be the most stable and support the experimental IR identification of this species. The energy barrier for elimination of H(2) from the second most stable singlet HGa(μ-OH)Ga insertion isomer found to be 13.9 (12.9) kcal/mol is also consistent with the available experimental data.     

Measurement of the rate coefficient for collisional removal of O2(X3Sigma- g, upsilon=1) by O(3P)

We report a laboratory measurement of the rate coefficient for the collisional removal of O(2)(X(3)Sigma(g) (-),upsilon=1) by O((3)P) atoms. In the experiments, 266-nm laser light photodissociates ozone in a mixture of molecular oxygen and ozone. The photolysis step produces vibrationally excited O(2)(a(1)Delta(g)) that is rapidly converted to O(2)(X(3)Sigma(g) (-),upsilon=1-3) in a near-resonant electronic energy-transfer process with ground-state O(2). In parallel, a large amount of O((1)D) atoms is generated that promptly relaxes to O((3)P). Under the conditions of the experiments, only collisions with the photolytically produced O((3)P) atoms control the lifetime of O(2)(X(3)Sigma(g) (-),upsilon=1), because its removal by molecular oxygen at room temperature is extremely slow. Tunable 193-nm laser light monitors the temporal evolution of the O(2)(X(3)Sigma(g) (-),upsilon=1) population by detection of laser-induced fluorescence near 360 nm. The removal rate coefficient for O(2)(X(3)Sigma(g) (-),upsilon=1) by O((3)P) atoms is (3.2+/-1.0)x10(-12) cm(3) s(-1) (2sigma) at a temperature of 315+/-15 K (2sigma). This result is essential for the analysis and correct interpretation of the 6.3-mum H(2)O(nu(2)) band emission in the Earth's mesosphere and indicates that the deactivation of O(2)(X (3)Sigma(g) (-),upsilon=1) by O((3)P) atoms is significantly faster than the nominal values recently used in atmospheric models.