Direct ab initio molecular dynamics study on a microsolvated SN2 reaction of OH-(H2O) with CH3Cl

Reaction dynamics for a microsolvated SN2 reaction OH-(H2O)+CH3Cl have been investigated by means of the direct ab initio molecular dynamics method. The relative center-of-mass collision energies were chosen as 10, 15, and 25 kcal/mol. Three reaction channels were found as products. These are (1) a channel leading to complete dissociation (the products are CH3OH+Cl- +H2O: denoted by channel I), (2) a solvation channel (the products are Cl-(H2O)+CH3OH: channel II), and (3) a complex formation channel (the products are CH3OH...H2O+Cl-: channel III). The branching ratios for the three channels were drastically changed as a function of center-of-mass collision energy. The ratio of complete dissociation channel (channel I) increased with increasing collision energy, whereas that of channel III decreased. The solvation channel (channel II) was minor at all collision energies. The selectivity of the reaction channels and the mechanism are discussed on the basis of the theoretical results.     

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.     

Reaction pathways and free energy barriers for alkaline hydrolysis of insecticide 2-trimethylammonioethyl methylphosphonofluoridate and related organophosphorus compounds: electrostatic and steric effects

Reaction pathways and free energy barriers for alkaline hydrolysis of the highly neurotoxic insecticide 2-trimethylammonioethyl methylphosphonofluoridate and related organophosphorus compounds were studied by performing first-principles electronic structure calculations on representative methylphosphonofluoridates, (RO)CH3P(O)F, in which R = CH2CH2N+(CH3)3, CH3, CH2CH2C(CH3)3, CH2CH2CH(CH3)2, CH(CH3)CH2N+(CH3)3, and CH(CH3)CH2N(CH3)2. The dominant reaction pathway was found to be associated with a transition state in which the attacking nucleophile OH- and the leaving group F- are positioned on opposite sides of the plane formed by the three remaining atoms attached to the phosphorus in order to minimize the electrostatic repulsion between these two groups. The free energy barriers calculated for the rate-determining step of the dominant pathway are 12.5 kcal/mol when R = CH2CH2N+(CH3)3, 15.5 kcal/mol when R = CH3, 17.9 kcal/mol when R = CH2CH2C(CH3)3, 16.5 kcal/mol when R = CH2CH2CH(CH3)2, 13.4 kcal/mol when R = CH(CH3)CH2N+(CH3)3, and 18.7 kcal/mol when R = CH(CH(3))CH(2)N(CH(3))(2). The calculated free energy barriers are in good agreement with available experimentally derived activation free energies, i.e. 14.7 kcal/mol when R = CH(3), 13.4 kcal/mol when R = CH2CH2N+(CH3)3, and 13.9 kcal/mol when R = CH(CH3)CH2N+(CH3)3. A detailed analysis of the calculated energetic results and available experimental data suggests that the net charge of the molecule (M) being hydrolyzed is a prominent factor affecting the free energy barrier (DeltaG) for the alkaline hydrolysis of phosphodiesters, phosphonofluoridates, and related organophosphorus compounds. The electrostatic interactions between the attacking nucleophile OH- and the molecule M being hydrolyzed favor such an order of the free energy barrier: DeltaG(M(+)+OH-) < DeltaG(M0+OH-) < DeltaG(M(-)+OH-), where M+, M0, and M- represent the cationic, neutral, and anionic molecules, respectively. The change of the substituent R in (RO)CH(3)P(O)F from CH3 to CH2CH2N+(CH3)3 is associated with both the electrostatic and steric effects on the free energy barrier, but the electrostatic effect dominates the substituent shift of the free energy barrier. This helps to better understand why the alkaline hydrolysis of (RO)CH3P(O)F with R = CH2CH2N+(CH3)3 and CH(CH3)CH2N+(CH3)3 is significantly faster than that with R = CH3. The effect of electrostatic interaction also helps to understand why the rate constants for the alkaline hydrolysis of phosphodiesters, such as intramolecular second messenger adenosine 3',5'-phosphate (cAMP), are generally smaller than those for the alkaline hydrolysis of the phosphonofluoridates and related phosphotriesters.     

Insight into selected reactions in low-temperature dimethyl ether combustion from Born-Oppenheimer molecular dynamics

Dimethyl ether is under consideration as an alternative diesel fuel. Its combustion chemistry is as yet ill-characterized. Here we use Born-Oppenheimer molecular dynamics (BOMD) based on DFT-B3LYP forces to investigate the short-time dynamics of selected features of the low-temperature dimethyl ether (DME) oxidation potential energy surface. Along the chain propagation pathway, we run BOMD simulations from the transition state involving the decomposition of (*)CH(2)OCH(2)OOH to two CH(2)=O and an (*)OH radical. We predict that formaldehyde C-O stretch overtones are excited, consistent with laser photolysis experiments. We also predict that O-H overtones are excited for the (*)OH formed from (*)CH(2)OCH(2)OOH dissociation. We also investigate short-time dynamics involved in chain branching. First, we examine the isomerization transition state of (*)OOCH(2)OCH(2)OOH --> HOOCH(2)OCHOOH. The latter species is predicted to be a short-lived metastable radical that decomposes within 500 fs to hydroperoxymethyl formate (HPMF; HOOCH(2)OC(=O)H) and the first (*)OH of chain branching. The dissociation of HOOCH(2)OCHOOH exhibits non-RRKM behavior in its lifetime profile, which may be due to conformational constraints or slow intramolecular vibrational energy transfer (IVR) from the nascent H-O bond to the opposite end of the radical, where O-O scission occurs to form HPMF and (*)OH. In a few trajectories, we see HOOCH(2)OCHOOH recross back to (*)OOCH(2)OCH(2)OOH because the isomerization is endothermic, with only an 8 kcal/mol barrier to recrossing. Therefore, some inhibition of chain-branching may be due to recrossing. Second, trajectories run from the transition state leading to the direct decomposition of HPMF (an important source of the second (*)OH radical in chain branching) to HCO, (*)OH, and HC(=O)OH show that these products can recombine to form many other possible products. These products include CH(2)OO + HC(=O)OH, H(2)O + CO + HC(=O)OH, HC(=O)OH + HC(=O)OH, and HC(=O)C(=O)H + H(2)O, which (save CH(2)OO + HC(=O)OH) are all more thermodynamically stable than the original HCO + (*)OH + HC(=O)OH products. Moreover, the multitude of extra products suggest that standard statistical rate theories cannot completely describe the reaction kinetics of significantly oxygenated compounds such as HPMF. These secondary products consume the second (*)OH required for explosive combustion, suggesting an inhibition of DME fuel combustion is likely.     

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.     

Unimolecular dissociation and thermochemistry of CH3OOH

The unimolecular dissociation of CH3OOH is investigated by exciting the molecule in the region of its 5nu(OH) band and probing the resulting OH fragments using laser-induced fluorescence. The measured OH fragment rotational and translational energies are used to determine the CH3O-OH bond dissociation energy, which we estimate to be approximately 42.6+/-1 kcal/mol. Combining this value with the known heats of formation of the fragments also gives an estimate for the heat of formation of CH3OOH which at 0 K we determine to be deltaH(f)0=-27+/-1 kcal/mol. This experimental value is in good agreement with the results of ab initio calculations carried out at the CCSD(T)/complete basis set limit which finds the heat of formation of CH3OOH at 0 K to be deltaH(f)0=-27.3 kcal/mol.     

Base-induced decomposition of alkyl hydroperoxides in the gas phase. Part 3. Kinetics and dynamics in HO- + CH3OOH, C2H5OOH, and tert-C4H9OOH reactions

The E(CO)2 elimination reactions of alkyl hydroperoxides proceed via abstraction of an alpha-hydrogen by a base: X(-) + R(1)R(2)HCOOH --> HX + R(1)R(2)C=O + HO(-). Efficiencies and product distributions for the reactions of the hydroxide anion with methyl, ethyl, and tert-butyl hydroperoxides are studied in the gas phase. On the basis of experiments using three isotopic analogues, HO(-) + CH3OOH, HO(-) + CD3OOH, and H(18)O(-) + CH3OOH, the overall intrinsic reaction efficiency is determined to be 80% or greater. The E(CO)2 decomposition is facile for these methylperoxide reactions, and predominates over competing proton transfer at the hydroperoxide moiety. The CH3CH2OOH reaction displays a similar E(CO)2 reactivity, whereas proton transfer and the formation of HOO(-) are the exclusive pathways observed for (CH3)3COOH, which has no alpha-hydrogen. All results are consistent with the E(CO)2 mechanism, transition state structure, and reaction energy diagrams calculated using the hybrid density functional B3LYP approach. Isotope labeling for HO(-) + CH3OOH also reveals some interaction between H2O and HO(-) within the E(CO)2 product complex [H2O...CH2=O...HO(-)]. There is little evidence, however, for the formation of the most exothermic products H2O + CH2(OH)O(-), which would arise from nucleophilic condensation of CH2=O and HO(-). The results suggest that the product dynamics are not totally statistical but are rather direct after the E(CO)2 transition state. The larger HO(-) + CH3CH2OOH system displays more statistical behavior during complex dissociation.     

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).     

Reactions of halogenated hydroperoxides and peroxyl and alkoxyl radicals from isoflurane in aqueous solution

Model systems, based on aqueous solutions containing isoflurane (CHF(2)OCHClCF(3)) as an example, have been studied in the presence and absence of methionine (MetS) to evaluate reactive fates of halogenated hydroperoxides and peroxyl and alkoxyl radicals. Primary peroxyl radicals, CHF(2)OCH(OO*)CF(3), generated upon 1-e-reduction of isoflurane react quantitatively with MetS via an overall two-electron oxidation mechanism to the corresponding sulfoxide (MetSO). This reaction is accompanied by the formation of oxyl radicals CHF(2)OCH(O*)CF(3) that quantitatively rearrange by a 1,2-hydrogen shift to CHF(2)OC*(OH)CF(3). According to quantum-chemical calculations, this reaction is exothermic (DeltaH = -5.1 kcal/mol) in contrast to other potentially possible pathways. These rearranged CHF(2)OC*(OH)CF(3) radicals react further via either of two pathways: (i) direct addition of oxygen or (ii) deprotonation followed by fluoride elimination resulting in CHF(2)OC(O)CF(2)*. Route i yields the corresponding CHF(2)OC(OO*)(OH)CF(3) peroxyl radicals, which eliminate H+/O(2)*-. The resulting ester, CHF(2)OC(O)CF(3), hydrolyzes further, accounting for the formation of HF, trifluoroacetic acid, and formic acid with a contribution of 45% and 80% in air- and oxygen-saturated solutions, respectively. A competitive pathway (ii) involves the reactions of the secondary peroxyl radicals, CHF(2)OC(O)CF(2)OO*. The two more stable of the three above mentioned peroxyl radicals can be distinguished through their reaction with MetS. Although the primary CHF(2)OCH(OO*)CF(3) oxidizes MetS to MetSO in a 2-e step, the majority of the secondarily formed CHF(2)OC(O)CF(2)OO* reacts with MetS via a 1-e transfer mechanism, yielding CHF(2)OC(O)CF(2)OO-, which eventually suffers a total breakup into CHF(2)O- + CO(2) + CF(2)O. Quantum-chemical calculations show that this reaction is highly exothermic (DeltaH = -81 kcal/mol). In air-saturated solution this pathway accounts for about 35% of the overall isoflurane degradation. Minor products (10% each), namely, oxalic acid and carbon monoxide originate from oxyl radicals, CHF(2)OC(O)CF(2)O* and CHF(2)OCH(O*)CF(3). An isoflurane-derived hydroperoxide CHF(2)OCH(OOH)CF(3) in high yield was generated in radiolysis of air-saturated solutions containing isoflurane and formate either via a H-atom abstraction from formate by the isoflurane-derived peroxyl radicals or by their cross-termination reaction with superoxide O(2)*-. CHF(2)OCH(OOH)CF(3), is an unstable intermediate whose multistep hydrolysis is giving H(2)O(2) + 2HF + HC(O)OH + CF(3)CH(OH)(2). In the absence of MetS, about 55% of CHF(2)OCH(OO*)CF(3) undergo termination via the Russell mechanism and 27% are involved in cross-termination with superoxide (O(2)*-) and peroxyl radicals derived from t-BuOH (used to scavenge *OH radicals). The remaining 18% of the primary peroxyl radicals undergo termination via formation of alkoxyl radicals, CHF(2)OCH(O*)CF(3).     

Theoretical study of the antioxidant properties of pyridoxine

Molecules acting as antioxidants capable of scavenging reactive oxygen species (ROS) are of utmost importance in the living cell. The antioxidative properties of pyridoxine (vitamin B6) have recently been discovered. In this study, we have analyzed the reactivity of pyridoxine toward the ROS (.-)OH, (.-)OOH, and (.-)O(2)- at the density functional theory level (functionals B3LYP and MPW1B95). Two reaction types have been studied as follows: addition to the aromatic ring atoms and hydrogen/proton abstraction. Our results show that (.-)OH is the most reactive species, while (.-)OOH displays low reactivity and (.-)O2(-) does not react at all with pyridoxine. The most exergonic reactions are those where (.-)H is removed from the CH(2)OH groups or the ring-bound OH group and range from -33 to -39 kcal/mol. The most exergonic addition reactions occur by attacking the carbon atoms bonded to nitrogen but with an energy gain of only 6 kcal/mol.     

Unimolecular decomposition of ethyl hydroperoxide: ab initio/Rice-Ramsperger-Kassel-Marcus theoretical prediction of rate constants

Alkyl hydroperoxides are found to be important intermediates in the combustion and oxidation processes of hydrocarbons. However, studies of ethyl hydroperoxide (CH(3)CH(2)OOH) are limited. In this work, kinetics and mechanisms for unimolecular decomposition of CH(3)CH(2)OOH have been investigated. The potential energy surface of decomposition reactions have first been predicted at the CCSD(T)/6-311+G(3df,2p)//B3LYP/6-311G(d,p) level. The results show that the formation of CH(3)CH(2)O + OH via O-O direct bond dissociation is dominant, the branching ratio of which is over 99% in the whole temperature range from 300 to 1000 K, and its rate constant can be expressed as k1 = 9.26 × 10(52)T(-11.91)exp(-26879/T) s(-1) at 1 atm. The rate constants of the reaction CH(3)CH(2)OOH → CH(3)CH(2)O + OH at different temperatures and pressures have been calculated, which can help us to comprehend the reactions of CH(3)CH(2)OOH at experimental conditions.     

Explanation of the unusual temperature dependence of the atmospherically important OH + H(2)S --> H(2)O + HS reaction and prediction of the rate constant at combustion temperatures

Rate constants for the OH + H2S --> H2O + HS reaction, which is important for both atmospheric chemistry and combustion, are calculated by direct dynamics with the M06-2X density functional using the MG3S basis set. Energetics are compared to high-level MCG3/3//MC-QCISD/3 wave function theory and to results obtained by other density functionals. We employ canonical variational transition-state theory with multidimensional tunneling contributions and scaled generalized normal-mode frequencies evaluated in redundant curvilinear coordinates with anharmonicity included in the torsion. The transition state has a quantum mechanically distinguishable, nonsuperimposable mirror image that corresponds to a separate classical reaction path; the effect of the multiple paths is examined through use of a symmetry number and by torsional methods. Calculations with the reference-potential Pitzer-Gwinn treatment of the torsional mode agree with experiment, within experimental scatter, and predict a striking temperature dependence of the activation energy, increasing from -0.1 kcal/mol at 200 K to 0.2, 1.0, 3.4, and 9.8 kcal/mol at 300, 500, 1000, and 2400 K. The unusual temperature dependence arises from a dynamical bottleneck at an energy below reactants, following an addition complex on the reaction path with a classical binding energy of 4.4 kcal/mol. As a way to check the mechanism, kinetic isotope effects of the OH + D2S and OD + D2S reactions have been predicted.     

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.     

Theoretical mechanism study of UF6 hydrolysis in the gas phase

The mechanism of the gas-phase reaction UF 6 + H 2O --> UOF 4 + 2HF is explored using relativistic density functional theory calculations. Initially, H 2O coordinates with UF 6 to form a 1:1 complex UF 6.H 2O. Over an activation energy barrier of about 19 kcal/mol, H 2O transfers a H atom to a nearby ligand F, resulting in UF 5OH + HF. The eliminated HF or another H 2O molecule may form a hydrogen bond with UF 5OH. Starting from UF 5OH, the second HF elimination results in UOF 4. If UF 5OH is in the isolated form, UF 5OH --> UOF 4 + HF takes place over a barrier of 24 kcal/mol. If UF 5OH is hydrogen-bonded with H 2O or HF, the conversion barrier is less than 10 kcal/mol. Once formed, the unstable UOF 4 tends to associate with additional ligands and hydrogen-bonding donors. The calculated binding energies indicate the significance of such interactions, which may have profound impact on further HF eliminating reactions. The IR spectra features can be used to indicate the formation and interaction type of the intermediates and products.     

Proton affinity of methyl peroxynitrate

The equilibrium structures, harmonic vibrational frequencies of methyl peroxynitrate, and structures of protonated methyl peroxynitrate have been investigated using ab initio methods. The methods include the single- and double-excitation quadratic configuration (QCISD) methods and the QCISD(T) method, which incorporates a perturbational estimate of the effects of corrected triple excitation. The lowest-energy gas-phase form of protonated methyl peroxynitrate is a complex between CH3OOH and NO2+. The CH3OOH.NO2+ complex is bound by 22 +/- 2 kcal/mol. The estimated proton affinity of methyl peroxynitrate is 178.8 +/- 3 kcal/mol. A general trend for the proton affinity of ROO-NO2 (peroxynitrates) compounds is discussed.     

Gas-surface scattering dynamics of CO2, NO2, and O3 in collisions with model organic surfaces

High-energy (70 kJ/mol) molecular beams of CO(2), NO(2), and O(3) were scattered from long-chain methyl (CH(3)-), hydroxyl (OH-), and perfluoro (CF(3)(CF(2))(8)-, or F-) ω-functionalized alkanethiol self-assembled monolayers (SAMs) on gold to study the dynamics of energy exchange and thermal accommodation of atmospherically important triatomic molecules on model organic surfaces. Overall, the extent of energy transfer in gas collisions with all of the surfaces studied was substantial. Specifically, the triatomics scatter from each surface only after dissipating greater than 80% of their incident energy. Furthermore, although the OH-SAM is a more rigid surface, the extent of energy transfer and accommodation of these molecules to the CH(3)- and OH-SAMs were approximately the same. The similar scattering dynamics are likely due to significant gas-surface attractive forces between the triatomics and the OH terminal groups, which compensate for the rigidity of this monolayer. In contrast to the OH- and CH(3)-SAMs, the dominant pathway in collisions of the gases with the F-SAM was impulsive scattering. The portion of molecules that accommodated (<40%) to the F-SAM was about half of the amount that accommodated (～70%) to the CH(3)- and OH-SAMs. Although differences in the surface properties had a significant effect on the dynamics, variances in the chemical and physical properties of the three gases, CO(2), NO(2), and O(3), were found to have little effect on the extent of energy transfer and accommodation for collisions with any one surface.     

Vibrational overtone induced elimination reactions within hydrogen-bonded molecular clusters: the dynamics of water catalyzed reactions in CH2FOH.H2On

The dynamics of the light initiated OH-overtone induced elimination reactions CH(2)FOH.(H(2)O)(n) + hnu--> HF + CH(2)O + n(H(2)O), n = 1-3, are studied using classical trajectory simulations where the ab initio potential energy surface is computed "on-the-fly". Hydrogen bonding to the water is found to lower the barrier to reaction by over 20 kcal mol(-1) and modifies the mechanism to a concerted multiple H-atom transfer process. The reaction process is found to occur on a rapid timescale, <100 fs, and involves the hydronium ion as an intermediate. An essential aspect of dynamics is the successful competition of reaction with energy dissipation through water evaporation from the cluster.     

The Stability of α‐Hydroperoxyalkyl Radicals

High‐level ab initio and Born–Oppenheimer molecular dynamic calculations have been carried out on a series of hydroperoxyalkyl (α‐QOOH) radicals with the aim of investigating the stability and unimolecular decomposition mechanism into QO+OH of these species. Dissociation was shown to take place through rotation of the C?O(OH) bond rather than through elongation of the CO?OH bond. Through the C?O(OH) rotation, the unpaired electron of the radical overlaps with the electron density on the O?OH bond, and from this overlap the C=O π bond forms and the O?OH bond breaks spontaneously. The CH₂OOH, CH(CH₃)OOH, CH(OH)OOH, and α‐hydroperoxycycloheptadienyl radical were found to decompose spontaneously, but the CH(CHO)OOH has a decomposition energy barrier of 5.95 kcal mol^?1 owing to its steric and electronic features. The systems studied in this work provide the first insights into how structural and electronic effects govern the stabilizing influence on elusive α‐QOOH radicals.     

Proton affinity of methyl nitrite and methyl peroxynitrite: implications for measuring branching ratios of alkyl nitrates and nitrites

Geometry optimizations for methyl nitrite and methyl peroxynitrite, along with various protonated isomers for each, have been investigated using ab initio and density functional methods. The lowest energy structure for protonated methyl nitrite is a complex between CH3OH and NO(+). For methyl peroxynitrite, the lowest energy protonated structure is a complex between CH3OOH and NO(+). Their respective proton affinities are estimated to be 195.2 and 195.8 kcal/mol at the QCISD(T)/6-311++G(3df,3pd) level of theory. The results, compared with past studies, suggest an alternative method for directly measuring branching ratios for production of alkyl nitrates and nitrites.