Crossed molecular beam experiments and theoretical simulations on the multichannel reaction of toluene with atomic oxygen.

Despite extensive experimental and theoretical studies on the kinetics of the O(3P) + C7H8 (toluene) reaction and a pioneering crossed molecular beam (CMB) investigation, the branching fractions (BFs) of the CH3C6H4O(methylphenoxy) + H, C6H5O(phenoxy) + CH3, and spin-forbidden C5H5CH3 (methylcyclopentadiene) + CO product channels remain an open question, which has hampered the proper inclusion of this important reaction in the chemical modelling of various chemical environments. We report a CMB study with universal soft electron-ionization mass-spectrometric detection of the reactions O(3P,1D) + toluene at the collision energy of 34.7 kJ mol-1. From CMB data we have inferred the reaction dynamics and quantified the BFs of the primary products and the role of intersystem crossing (ISC). The CH3-elimination channel dominates (BF = 0.69 ± 0.22) in the O(3P) reaction, while the H-displacement and CO-formation channels are minor (BF = 0.22 ± 0.07 and 0.09 ± 0.05, respectively), with ISC accounting for more than 50% of the reactive flux. Synergistic transition-state theory (TST)-based master equation simulations including nonadiabatic TST on ab initio coupled triplet/singlet potential energy surfaces were employed to compute the product BFs and assist in the interpretation of the CMB results. In the light of the good agreement between the theoretical predictions for the O(3P) + toluene reaction and the CMB results as well as the absolute rate constant as a function of temperature (T) (from literature), the so-validated computational methodology was used to predict channel-specific rate constants as a function of T at 1 atm.

An experimental and theoretical investigation of the N(2D) + C6H6 (benzene) reaction with implications for the photochemical models of Titan.

We report on a combined experimental and theoretical investigation of the N(2D) + C6H6 (benzene) reaction, which is of relevance in the aromatic chemistry of the atmosphere of Titan. Experimentally, the reaction was studied (i) under single-collision conditions by the crossed molecular beams (CMB) scattering method with mass spectrometric detection and time-of-flight analysis at the collision energy (Ec) of 31.8 kJ mol-1 to determine the primary products, their branching fractions (BFs), and the reaction micromechanism, and (ii) in a continuous supersonic flow reactor to determine the rate constant as a function of temperature from 50 K to 296 K. Theoretically, electronic structure calculations of the doublet C6H6N potential energy surface (PES) were performed to assist the interpretation of the experimental results and characterize the overall reaction mechanism. The reaction is found to proceed via barrierless addition of N(2D) to the aromatic ring of C6H6, followed by formation of several cyclic (five-, six-, and seven-membered ring) and linear isomeric C6H6N intermediates that can undergo unimolecular decomposition to bimolecular products. Statistical estimates of product BFs on the theoretical PES were carried out under the conditions of the CMB experiments and at the temperatures relevant for Titan's atmosphere. In all conditions the ring-contraction channel leading to C5H5 (cyclopentadienyl) + HCN is dominant, while minor contributions come from the channels leading to o-C6H5N (o-N-cycloheptatriene radical) + H, C4H4N (pyrrolyl) + C2H2 (acetylene), C5H5CN (cyano-cyclopentadiene) + H, and p-C6H5N + H. Rate constants (which are close to the gas kinetic limit at all temperatures, with the recommended value of 2.19 ± 0.30 × 10-10 cm3 s-1 over the 50-296 K range) and BFs have been used in a photochemical model of Titan's atmosphere to simulate the effect of the title reaction on the species abundances as a function of the altitude.

A combined crossed molecular beam and theorerical study of the O(3P,1D) + acrylonitrile (CH2CHCN) reactions and implications for combustion and extraterrestrial environments.

Acrylonitrile (CH2CHCN) is ubiquitous in space (molecular clouds, solar-type star forming regions, and circumstellar envelopes) and is also abundant in the upper atmosphere of Titan. The reaction O(3P) + CH2CHCN can be of relevance in the chemistry of the interstellar medium because of the abundance of atomic oxygen. The oxidation of acrylonitrile is also important in combustion as the thermal decomposition of pyrrolic and pyridinic structures present in fuel-bound nitrogen generates many nitrogen-bearing compounds, including acrylonitrile. Despite its relevance, limited information exists on this reaction. We report a combined experimental and theoretical investigation of the reactions of acrylonitrile with both ground 3P and excited 1D atomic oxygen. From product angular and time-of-flight distributions in crossed molecular beam experiments with mass spectrometric detection at a collision energy, Ec, of 31.4 kJ mol-1, we have identified the primary reaction products and determined their branching fractions (BFs). Theoretical calculations of the relevant triplet and singlet potential energy surfaces (PESs) were performed to interpret the experimental results and elucidate the reaction mechanism. Adiabatic statistical calculations of product BFs for the decomposition of the main triplet and singlet intermediates have been carried out. Combining the experimental and theoretical results, we conclude that the O(3P) reaction leads to two main product channels: (i) CH2CNH (ketenimine) + CO (dominant with a BF of 0.87 ± 0.05), formed via efficient intersystem crossing from the entrance triplet PES to the underlying singlet PES, and (ii) HCOCHCN + H (minor, with a BF of 0.13 ± 0.05), occurring adiabatically on the triplet PES. Our study suggests the inclusion of this reaction as a possible destruction pathway of CH2CHCN and a possible formation route of CH2CNH in the interstellar medium. The O(1D) + CH2CHCN reaction mainly leads to the formation of CH2CNH + CO adiabatically on the singlet PES. This result can improve models related to the chemistry of interstellar ice and cometary comas, where O(1D) reactions can play a role. Overall, our results are expected to be useful for improving the models of combustion and extraterrestrial environments.

Reactions O(3P, 1D) + HCCCN(X1Σ+) (Cyanoacetylene): Crossed-Beam and Theoretical Studies and Implications for the Chemistry of Extraterrestrial Environments.

Cyanoacetylene (HCCCN), the first member of the cyanopolyyne family (HCnN, where n = 3, 5, 7, ...), is of particular interest in astrochemistry being ubiquitous in space (molecular clouds, solar-type protostars, protoplanetary disks, circumstellar envelopes, and external galaxies) and also relatively abundant. It is also abundant in the upper atmosphere of Titan and comets. Since oxygen is the third most abundant element in space, after hydrogen and helium, the reaction O + HCCCN can be of relevance in the chemistry of extraterrestrial environments. Despite that, scarce information exists not only on the reactions of oxygen atoms with cyanoacetylene but with nitriles in general. Here, we report on a combined experimental and theoretical investigation of the reactions of cyanoacetylene with both ground 3P and excited 1D atomic oxygen and provide detailed information on the primary reaction products, their branching fractions (BFs), and the overall reaction mechanisms. More specifically, the reactions of O(3P, 1D) with HCCCN(X1Σ+) have been investigated under single-collision conditions by the crossed molecular beams scattering method with mass spectrometric detection and time-of-flight analysis at the collision energy, Ec, of 31.1 kJ/mol. From product angular and time-of-flight distributions, we have identified the primary reaction products and determined their branching fractions (BFs). Theoretical calculations of the relevant triplet and singlet potential energy surfaces (PESs) were performed to assist the interpretation of the experimental results and clarify the reaction mechanism. Adiabatic statistical calculations of product BFs for the decomposition of the main triplet and singlet intermediates have also been carried out. Merging together the experimental and theoretical results, we conclude that the O(3P) reaction is characterized by a minor adiabatic channel leading to OCCCN (cyanoketyl) + H (experimental BF = 0.10 ± 0.05), while the dominant channel (BF = 0.90 ± 0.05) occurs via intersystem crossing to the underlying singlet PES and leads to formation of 1HCCN (cyanomethylene) + CO. The O(1D) reaction is characterized by the same two channels, with the relative CO/H yield being slightly larger. Considering the recorded reactive signal and the calculated entrance barrier, we estimate that the rate coefficient for reaction O(3P) + HC3N at 300 K is in the 10-12 cm3 molec-1 s-1 range. Our results are expected to be useful to improve astrochemical and photochemical models. In addition, they are also relevant in combustion chemistry, because the thermal decomposition of pyrrolic and pyridinic structures present in fuel-bound nitrogen generates many nitrogen-bearing compounds, including cyanoacetylene.

OH(2Π) + C2H4 Reaction: A Combined Crossed Molecular Beam and Theoretical Study.

The reaction between the ground-state hydroxyl radical, OH(2Π), and ethylene, C2H4, has been investigated under single-collision conditions by the crossed molecular beam scattering technique with mass-spectrometric detection and time-of-flight analysis at the collision energy of 50.4 kJ/mol. Electronic structure calculations of the underlying potential energy surface (PES) and statistical Rice-Ramsperger-Kassel-Marcus (RRKM) calculations of product branching fractions on the derived PES for the addition pathway have been performed. The theoretical results indicate a temperature-dependent competition between the anti-/syn-CH2CHOH (vinyl alcohol) + H, CH3CHO (acetaldehyde) + H, and H2CO (formaldehyde) + CH3 product channels. The yield of the H-abstraction channel could not be quantified with the employed methods. The RRKM results predict that under our experimental conditions, the anti- and syn-CH2CHOH + H product channels account for 38% (in similar amounts) of the addition mechanism yield, the H2CO + CH3 channel for ∼58%, while the CH3CHO + H channel is formed in negligible amount (<4%). The implications for combustion and astrochemical environments are discussed.

Unsaturated Dinitriles Formation Routes in Extraterrestrial Environments: A Combined Experimental and Theoretical Investigation of the Reaction between Cyano Radicals and Cyanoethene (C2H3CN).

The reaction between cyano radicals (CN, X2Σ+) and cyanoethene (C2H3CN) has been investigated by a combined approach coupling crossed molecular beam (CMB) experiments with mass spectrometric detection and time-of-flight analysis at a collision energy of 44.6 kJ mol-1 and electronic structure calculations to determine the relevant potential energy surface. The experimental results can be interpreted by assuming the occurrence of a dominant reaction pathway leading to the two but-2-enedinitrile (1,2-dicyanothene) isomers (E- and Z-NC-CH═CH-CN) in a H-displacement channel and, to a much minor extent, to 1,1-dicyanoethene, CH2C(CN)2. In order to derive the product branching ratios under the conditions of the CMB experiments and at colder temperatures, including those relevant to Titan and to cold interstellar clouds, we have carried out RRKM statistical calculations using the relevant potential energy surface of the investigated reaction. We have also estimated the rate coefficient at very low temperatures by employing a semiempirical method for the treatment of long-range interactions. The reaction has been found to be barrierless and fast also under the low temperature conditions of cold interstellar clouds and the atmosphere of Titan. Astrophysical implications and comparison with literature data are also presented. On the basis of the present work, 1,2-dicyanothene and 1,1-dicyanothene are excellent candidates for the search of dinitriles in the interstellar medium.

The N(2D) + CH2CHCN (Vinyl Cyanide) Reaction: A Combined Crossed Molecular Beam and Theoretical Study and Implications for the Atmosphere of Titan.

The reaction of electronically excited nitrogen atoms, N(2D), with vinyl cyanide, CH2CHCN, has been investigated under single-collision conditions by the crossed molecular beam (CMB) scattering method with mass spectrometric detection and time-of-flight (TOF) analysis at the collision energy, Ec, of 31.4 kJ/mol. Synergistic electronic structure calculations of the doublet potential energy surface (PES) have been performed to assist in the interpretation of the experimental results and characterize the overall reaction micromechanism. Statistical (Rice-Ramsperger-Kassel-Marcus, RRKM) calculations of product branching fractions (BFs) on the theoretical PES have been carried out at different values of temperature, including the one corresponding to the temperature (175 K) of Titan's stratosphere and at a total energy corresponding to the Ec of the CMB experiment. According to our theoretical calculations, the reaction is found to proceed via barrierless addition of N(2D) to the carbon-carbon double bond of CH2═CH-CN, followed by the formation of cyclic and linear intermediates that can undergo H, CN, and HCN elimination. In competition, the N(2D) addition to the CN group is also possible via a submerged barrier, leading ultimately to N2 + C3H3 formation, the most exothermic of all possible channels. Product angular and TOF distributions have been recorded for the H-displacement channels leading to the formation of a variety of possible C3H2N2 isomeric products. Experimentally, no evidence of CN, HCN, and N2 forming channels was observed. These findings were corroborated by the theory, which predicts a variety of competing product channels, following N(2D) addition to the double bond, with the main ones, at Ec = 31.4 kJ/mol, being six isomeric H forming channels: c-CH(N)CHCN + H (BF = 35.0%), c-CHNCHCN + H (BF = 28.1%), CH2NCCN + H (BF = 26.3%), c-CH2(N)CCN(cyano-azirine) + H (BF = 7.4%), trans-HNCCHCN + H (BF = 1.6%), and cis-HNCCHCN + H (BF = 1.3%), while C-C bond breaking channels leading to c-CH2(N)CH(2H-azirine) + CN and c-CH2(N)C + HCN are predicted to be negligible (0.02% and 0.2%, respectively). The highly exothermic N2 + CH2CCH (propargyl) channel is also predicted to be negligible because of the very high isomerization barrier from the initial addition intermediate to the precursor intermediate able to lead to products. The predicted product BFs are found to have, in general, a very weak energy dependence. The above cyclic and linear products containing an additional C-N bond could be potential precursors of more complex, N-rich organic molecules that contribute to the formation of the aerosols on Titan's upper atmosphere. Overall, the results are expected to have a significant impact on the gas-phase chemistry of Titan's atmosphere and should be properly included in the photochemical models.

Tracing the Primordial Chemical Life of Glycine: A Review from Quantum Chemical Simulations.

Glycine (Gly), NH2CH2COOH, is the simplest amino acid. Although it has not been directly detected in the interstellar gas-phase medium, it has been identified in comets and meteorites, and its synthesis in these environments has been simulated in terrestrial laboratory experiments. Likewise, condensation of Gly to form peptides in scenarios resembling those present in a primordial Earth has been demonstrated experimentally. Thus, Gly is a paradigmatic system for biomolecular building blocks to investigate how they can be synthesized in astrophysical environments, transported and delivered by fragments of asteroids (meteorites, once they land on Earth) and comets (interplanetary dust particles that land on Earth) to the primitive Earth, and there react to form biopolymers as a step towards the emergence of life. Quantum chemical investigations addressing these Gly-related events have been performed, providing fundamental atomic-scale information and quantitative energetic data. However, they are spread in the literature and difficult to harmonize in a consistent way due to different computational chemistry methodologies and model systems. This review aims to collect the work done so far to characterize, at a quantum mechanical level, the chemical life of Gly, i.e., from its synthesis in the interstellar medium up to its polymerization on Earth.

The Reaction N(2D) + CH3CCH (Methylacetylene): A Combined Crossed Molecular Beams and Theoretical Investigation and Implications for the Atmosphere of Titan.

The reaction of excited nitrogen atoms N(2D) with CH3CCH (methylacetylene) was investigated under single-collision conditions by the crossed molecular beams (CMB) scattering method with mass spectrometric detection and time-of-flight analysis at the collision energy (Ec) of 31.0 kJ/mol. Synergistic electronic structure calculations of the doublet potential energy surface (PES) were performed to assist the interpretation of the experimental results and characterize the overall reaction micromechanism. Theoretically, the reaction is found to proceed via a barrierless addition of N(2D) to the carbon-carbon triple bond of CH3CCH and an insertion of N(2D) into the CH bond of the methyl group, followed by the formation of cyclic and linear intermediates that can undergo H, CH3, and C2H elimination or isomerize to other intermediates before unimolecularly decaying to a variety of products. Kinetic calculations for addition and insertion mechanisms and statistical (Rice-Ramsperger-Kassel-Marcus) computations of product branching fractions (BFs) on the theoretical PES were performed at different values of total energy, including the one corresponding to the temperature (175 K) of Titan's stratosphere and that of the CMB experiment. Up to 14 competing product channels were statistically predicted, with the main ones, at Ec = 31.0 kJ/mol, being the formation of CH2NH (methanimine) + C2H (ethylidyne) (BF = 0.41), c-C(N)CH + CH3 (BF = 0.32), CH2CHCN (acrylonitrile) + H (BF = 0.12), and c-CH2C(N)CH + H (BF = 0.04). Of the 14 possible channels, seven correspond to H displacement channels of different exothermicity, for a total H channel BF of ∼0.25 at Ec = 31.0 kJ/mol. Experimentally, dynamical information could only be obtained about the overall H channels. In particular, the experiment corroborates the formation of acrylonitrile + H, which is the most exothermic of all 14 reaction channels and is theoretically calculated to be the dominant H-forming channel (BF = 0.12). The products containing a novel C-N bond could be potential precursors to form other nitriles (C2N2, C3N) or more complex organic species containing N atoms in planetary atmospheres, such as those of Titan and Pluto. Overall, the results are expected to have a potentially significant impact on the understanding of the gas-phase chemistry of Titan's atmosphere and the modeling of that atmosphere.

Crossed-Beam and Theoretical Studies of the O(3P, 1D) + Benzene Reactions: Primary Products, Branching Fractions, and Role of Intersystem Crossing.

Reliable modeling of hydrocarbon oxidation relies critically on knowledge of the branching fractions (BFs) as a function of temperature (T) and pressure (p) for the products of the reaction of the hydrocarbon with atomic oxygen in its ground state, O(3P). During the past decade, we have performed in-depth investigations of the reactions of O(3P) with a variety of small unsaturated hydrocarbons using the crossed molecular beam (CMB) technique with universal mass spectrometric (MS) detection and time-of-flight (TOF) analysis, combined with synergistic theoretical calculations of the relevant potential energy surfaces (PESs) and statistical computations of product BFs, including intersystem crossing (ISC). This has allowed us to determine the primary products, their BFs, and extent of ISC to ultimately provide theoretical channel-specific rate constants as a function of T and p. In this work, we have extended this approach to the oxidation of one of the most important species involved in the combustion of aromatics: the benzene (C6H6) molecule. Despite extensive experimental and theoretical studies on the kinetics and dynamics of the O(3P) + C6H6 reaction, the relative importance of the C6H5O (phenoxy) + H open-shell products and of the spin-forbidden C5H6 (cyclopentadiene) + CO and phenol adduct closed-shell products are still open issues, which have hampered the development of reliable benzene combustion models. With the CMB technique, we have investigated the reaction dynamics of O(3P) + benzene at a collision energy (Ec) of 8.2 kcal/mol, focusing on the occurrence of the phenoxy + H and spin-forbidden C5H6 + CO and phenol channels in order to shed further light on the dynamics of this complex and important reaction, including the role of ISC. Concurrently, we have also investigated the reaction dynamics of O(1D) + benzene at the same Ec. Synergistic high-level electronic structure calculations of the underlying triplet/singlet PESs, including nonadiabatic couplings, have been performed to complement and assist the interpretation of the experimental results. Statistical (RRKM)/master equation (ME) computations of the product distribution and BFs on these PESs, with inclusion of ISC, have been performed and compared to experiment. In light of the reasonable agreement between the CMB experiment, literature kinetic experimental results, and theoretical predictions for the O(3P) + benzene reaction, the so-validated computational methodology has been used to predict (i) the BF between the C6H5O + H and C5H6 + CO channels as a function of collision energy and temperature (at 0.1 and 1 bar), showing that their increase progressively favors radical (phenoxy + H)-forming over molecule (C5H6 + CO and phenol stabilization)-forming channels, and (ii) channel-specific rate constants as a function of T and p, which are expected to be useful for improved combustion models.

Monitoring the Reactivity of Formamide on Amorphous SiO2 by In-Situ UV-Raman Spectroscopy and DFT Modeling.

Formamide has been recognized in the literature as a key species in the formation of the complex molecules of life, such as nucleobases. Furthermore, several studies reported the impact of mineral phases as catalysts for its decomposition/polymerization processes, increasing the conversion and also favoring the formation of specific products. Despite the progresses in the field, in situ studies on these mineral-catalyzed processes are missing. In this work, we present an in situ UV-Raman characterization of the chemical evolution of formamide over amorphous SiO2 samples, selected as a prototype of silicate minerals. The experiments were carried out after reaction of formamide at 160 °C on amorphous SiO2 (Aerosil OX50) either pristine or pre-calcined at 450 °C, to remove a large fraction of surface silanol groups. Our measurements, interpreted on the basis of density functional B3LYP-D3 calculations, allow to assign the spectra bands in terms of specific complex organic molecules, namely, diaminomaleonitrile (DAMN), 5-aminoimidazole (AI), and purine, showing the role of the mineral surface on the formation of relevant prebiotic molecules.

Combined Experimental and Theoretical Studies of the O(3P) + 1-Butene Reaction Dynamics: Primary Products, Branching Fractions, and Role of Intersystem Crossing.

Information on the detailed mechanism and dynamics (primary products, branching fractions (BFs), and channel specific rate constants as a function of temperature) for many important combustion reactions of O(3P) with unsaturated hydrocarbons is still lacking. We report synergistic experimental/theoretical studies on the mechanism and dynamics of the O(3P) + 1-C4H8 (1-butene) reaction by combining crossed molecular beam (CMB) experiments with soft electron ionization mass-spectrometric detection and time-of-flight analysis at 10.5 kcal/mol collision energy (Ec) to high-level ab initio electronic structure calculations of the underlying triplet and singlet potential energy surfaces (PESs) and statistical Rice-Ramsperger-Kassel-Marcus/Master Equation (RRKM/ME) computations of BFs including intersystem crossing (ISC). The reactive interaction of O(3P) with 1-butene is found to mainly break apart the 4-carbon atom chain, leading to the radical product channels ethyl + vinoxy (BF = 0.34 ± 0.11), methyl + C3H5O (BF = 0.28 ± 0.09), formyl + propyl (BF = 0.17 ± 0.05), ethyl + acetyl (BF = 0.014 ± 0.007), and butanal radical (ethylvinoxy) + H (BF = 0.013 ± 0.004), and molecular product channels formaldehyde + propenylidene/propene (BF = 0.15 ± 0.05) and butenone (ethyl ketene) + H2 (BF = 0.037 ± 0.015). As some of these products can only be formed via ISC from triplet to singlet PESs, from BFs an extent of ISC of 50% is inferred. This value is significantly larger than that recently observed for O(3P) + propene (22%) at similar Ec, underlying the question of how important it is to consider nonadiabatic effects for these and similar combustion reactions. Comparison of the derived BFs with those of statistical (RRKM/ME) simulations on the ab initio coupled triplet/singlet PESs shows good agreement, warranting the use of the RRKM/ME approach to provide information on the variation of BFs with temperature and to derive channel specific rate constants as a function of temperature (T) and pressure (p). Notably, ISC is predicted to decrease strongly with increasing temperature (from about 70% at 300 K to 46% at Ec = 10.5 kcal/mol, and about 1% at 2000 K). The present results lead to a detailed understanding of the complex reaction mechanism of O(3P) + 1-butene and, by providing channel specific rate constants as a function of T and p, should facilitate the improvement of current fossil-fuel (1-butene) as well as biofuel (1-butanol) combustion models.

Exploring the Gas Phase Synthesis of the Elusive Class of Boronyls and the Mechanism of Boronyl Radical Reactions under Single Collision Conditions.

Until recently, the chemistry of boronyl (BO), a diatomic radical isolectronic with the cyano (CN) species, has remained unknown. The boronyl group is characterized by a boron-oxygen multiple bond, and because of the inherent electron deficiency of the boron atom, boronyls (RBO) are highly reactive and typically only exist in their cyclotrimeric form (RBO)3. Due to their invaluable role as reactants, the isolation of the monomers in gas phase experiments has been extensively sought after by the organic synthesis and physical organic chemistry communities but never achieved. Besides the interests from a physical organic and synthetic point of view, boronyls also play a role as reaction intermediates in boron-assisted rocket propulsion systems. In this Account, we review recent experimental work in which gas phase organo boronyl monomers (RBO) are formed via bimolecular reactions of the boronyl radical (BO) with C2-C6 unsaturated hydrocarbons. The investigated hydrocarbons are widely exploited as fuels, and their reactions with boronyl radicals under single collision conditions lead to the formation of organo boronyls. Our studies also elucidate the mechanisms of their formation reactions thus furnishing a comprehension at the molecular level of this reaction class. The variety of the employed hydrocarbon substrates has allowed us to systematically classify the chemical behavior of the boronyl radicals. With the exception of the case of the dimethylacetylene reaction, the boron monoxide radical versus atomic hydrogen exchange mechanisms were always open leading to the formation of highly unsaturated organo boronyl monomers (RBO), which could be easily identified because they cannot trimerize under single collision conditions. Besides the hydrogen displacement pathway, methylacetylene, dimethylacetylene, and propylene, carrying one or two methyl groups, were also found to eliminate a methyl group. In all systems, the reactions were barrierless, indirect, and initiated by addition of the boron monoxide radical to the π electron density of the hydrocarbon molecule, with the radical center located at the boron atom of the BO radical, thus leading to doublet radical intermediates. These intermediates either decompose via hydrogen or methyl loss or isomerize prior to their decomposition via atomic hydrogen or migration of the BO moiety. A consistent trend suggests that all exit transition states are rather tight with those involved in the hydrogen atom loss depicting exit barriers of typically 25 to 35 kJ mol-1, whereas the methyl loss pathways are associated with tighter exit transition states located about 30-50 kJ mol-1 above the separated products. Further, the overall energetics suggest that those bimolecular reactions are exoergic by 40-90 kJ mol-1. These findings confirm that this reaction class leads to the formation of highly unsaturated organo boronyl molecules.

Low temperature kinetics and theoretical studies of the reaction CN + CH3NH2: a potential source of cyanamide and methyl cyanamide in the interstellar medium.

The reaction between cyano radicals (which are ubiquitous in interstellar clouds) and methylamine (a molecule detected in various interstellar sources) has been investigated in a synergistic experimental and theoretical study. The reaction has been found to be very fast in the entire range of temperatures investigated (23-297 K) by using a CRESU apparatus coupled to pulsed laser photolysis - laser induced fluorescence. The global experimental rate coefficient is given by In addition, dedicated electronic structure calculations of the underlying potential energy surface have been performed, together with capture theory and RRKM calculations. The experimental data have been interpreted in the light of the theoretical calculations and the product branching ratio has been established. According to the present study, in the range of temperatures investigated the title reaction is an efficient interstellar route of formation of cyanamide, NH2CN, another interstellar species. The second most important channel is the one leading to methyl cyanamide, CH3NHCN (an isomer of aminoacetonitrile), via a CN/H exchange mechanism with a yield of 12% of the global reaction in the entire range of temperatures explored. For a possible inclusion in future astrochemical models we suggest, by referring to the usual expression the following values: α = 3.68 × 10-12 cm3 molec-1 s-1, ß = -1.80, γ = 7.79 K for the channel leading to NH2CN + CH3; α = 5.05 × 10-13 cm3 molec-1 s-1, ß = -1.82, γ = 7.93 K for the channel leading to CH3NHCN + H.

Reaction Dynamics of O((3)P) + Propyne: II. Primary Products, Branching Ratios, and Role of Intersystem Crossing from Ab Initio Coupled Triplet/Singlet Potential Energy Surfaces and Statistical Calculations.

The mechanism of the O((3)P) + CH3CCH reaction was investigated using a combined experimental/theoretical approach. Experimentally the reaction dynamics was studied using crossed molecular beams (CMB) with mass-spectrometric detection and time-of-flight analysis at 9.2 kcal/mol collision energy. Theoretically master equation (ME) simulations were performed on a potential energy surface (PES) determined using high-level ab initio electronic structure calculations. In this paper (II) the theoretical results are described and compared with experiments, while in paper (I) are reported and discussed the results of the experimental study. The PES was investigated by determining structures and vibrational frequencies of wells and transition states at the CASPT2/aug-cc-pVTZ level using a minimal active space. Energies were then determined at the CASPT2 level increasing systematically the active space and at the CCSD(T) level extrapolating to the complete basis set limit. Two separate portions of the triplet PES were investigated, as O((3)P) can add either on the terminal or the central carbon of the unsaturated propyne bond. Minimum energy crossing points (MECPs) between the triplet and singlet PESs were searched at the CASPT2 level. The calculated spin-orbit coupling constants between the T1 and S0 electronic surfaces were ∼25 cm(-1) for both PESs. The portions of the singlet PES that can be accessed from the MECPs were investigated at the same level of theory. The system reactivity was predicted integrating stochastically the one-dimensional ME using Rice-Ramsperger-Kassel-Marcus theory to determine rate constants on the full T1/S0 PESs, accounting explicitly for intersystem crossing (ISC) using the Landau-Zener model. The computational results are compared both with the branching ratios (BRs) determined experimentally in the companion paper (I) and with those estimated in a recent kinetic study at 298 K. The ME results allow to interpret the main system reactivity: CH3CCO + H and CH3 + HCCO are the major channels active on the triplet PES and are formed from the wells accessed after O addition to the terminal and central C, respectively; (1)CH3CH + CO and C2H3 + HCO are the major channels active on the singlet PES and are formed from the methylketene and acrolein wells after ISC. However, also a large number of minor channels (∼15) are active, so that the system reactivity is quite complicated. The comparison between computational and experimental BRs is quite good for the kinetic study, while some discrepancy with the CMB estimations suggests that dynamic non-ergodic effects may influence the system reactivity. Channel specific rate constants are calculated in the 300-2250 K and 1-30 bar temperature and pressure ranges. It is found that as the temperature increases the H abstraction reaction, whose contribution is negligible in the experimental conditions, increases in relevance, and the extent of ISC decreases from ∼80% at 300 K to less than 2% at 2250 K.

Reaction Dynamics of O((3)P) + Propyne: I. Primary Products, Branching Ratios, and Role of Intersystem Crossing from Crossed Molecular Beam Experiments.

We performed synergic experimental/theoretical studies on the mechanism of the O((3)P) + propyne reaction by combining crossed molecular beams experiments with mass-spectrometric detection and time-of-flight analysis at 9.2 kcal/mol collision energy (Ec) with ab initio electronic structure calculations at a high level of theory of the relevant triplet and singlet potential energy surfaces (PESs) and statistical calculations of branching ratios (BRs) taking into account intersystem crossing (ISC). In this paper (I) we report the results of the experimental investigation, while the accompanying paper (II) shows results of the theoretical investigation with comparison to experimental results. By exploiting soft electron ionization detection to suppress/mitigate the effects of the dissociative ionization of reactants, products, and background gases, product angular and velocity distributions at different charge-to-mass ratios were measured. From the laboratory data angular and translational energy distributions in the center-of-mass system were obtained for the five competing most important product channels, and product BRs were derived. The reactive interaction of O((3)P) with propyne under single collision conditions is mainly leading to the rupture of the three-carbon atom chain, with production of the radical products methylketenyl + atomic hydrogen (BR = 0.04), methyl + ketenyl (BR = 0.10), and vinyl + formyl (BR = 0.11) and the molecular products ethylidene/ethylene + carbon monoxide (BR = 0.74) and propandienal + molecular hydrogen (BR = 0.01). Because some of the products can only be formed via ISC from the entrance triplet to the low-lying singlet PES, we infer from their BRs an amount of ISC larger than 80%. This value is dramatically large when compared to the negligible ISC reported for the O((3)P) reaction with the simplest alkyne, acetylene. At the same time, it is much larger than that (∼20%) recently observed in the related reaction of the three-carbon atom alkene, O((3)P) + propene at a comparable Ec. This poses the question of how important it is to consider nonadiabatic effects and their dependence on molecular structure for this kind of combustion reactions. The prevalence of the CH3 over the H displacement channels is not explained by invoking a preference for the addition on the methyl-substituted acetylenic carbon atom, but rather it is believed to be due to the different tendencies of the two addition triplet intermediates CH3CCHO (preferentially leading to H elimination) and CH3COCH (preferentially leading to CH3 elimination) to undergo ISC to the underlying singlet PES. It is concluded that the main coproduct of the CO forming channel is singlet ethylidene ((1)CH3CH) rather than ground-state ethylene. By comparing the derived BRs with those very recently derived from kinetics studies at room temperature and 4 Torr we obtained information on how BRs vary with collision energy. The extent of ISC is estimated to remain essentially constant (∼85%) with increasing Ec from ∼1 to ∼10 kcal/mol. The present experimental results shed light on the mechanism of the title reaction at energies comparable to those involved in combustion and, when compared with the results from the statistical calculations on the ab initio coupled PESs (see accompanying paper II), lead to an in-depth understanding of the rather complex reaction mechanism of O + propyne. The overall results are expected to contribute to the development of more refined models of hydrocarbon combustion.

Crossed Molecular Beams and Quasiclassical Trajectory Surface Hopping Studies of the Multichannel Nonadiabatic O((3)P) + Ethylene Reaction at High Collision Energy.

The combustion relevant O((3)P) + C2H4 reaction stands out as a prototypical multichannel nonadiabatic reaction involving both triplet and singlet potential energy surfaces (PESs), which are strongly coupled. Crossed molecular beam (CMB) scattering experiments with universal soft electron ionization mass spectrometric detection have been used to characterize the dynamics of this reaction at the relatively high collision energy Ec of 13.7 kcal/mol, attained by crossing the reactant beams at an angle of 135°. This work is a full report of the data at the highest Ec investigated for this reaction. From laboratory product angular and velocity distribution measurements, angular and translational energy distributions in the center-of-mass system have been obtained for the five observed exothermic competing reaction channels leading to H + CH2CHO, H + CH3CO, CH3 + HCO, CH2 + H2CO, and H2 + CH2CO. The product branching ratios (BRs) have been derived. The elucidation of the reaction dynamics is assisted by synergic full-dimensional quasiclassical trajectory surface-hopping calculations of the reactive differential cross sections on coupled ab initio triplet/singlet PESs. This joint experimental/theoretical study extends and complements our previous combined CMB and theoretical work at the lower collision energy of 8.4 kcal/mol. The theoretically derived BRs and extent of intersystem crossing (ISC) are compared with experimental results. In particular, the predictions of the QCT results for the three main channels (those leading to vinoxy + H, methyl + HCO and methylene + H2CO formation) are compared directly with the experimental data in the laboratory frame. Good overall agreement is noted between theory and experiment, although some small, yet significant shortcomings of the theoretical differential cross section are noted. Both experiment and theory find almost an equal contribution from the triplet and singlet surfaces to the reaction, with a clear tendency of the degree of ISC to decrease with increasing Ec and with theory slightly overestimating the extent of ISC.

Intersystem crossing and dynamics in O(3P) + C2H4 multichannel reaction: experiment validates theory.

The O((3)P) + C(2)H(4) reaction, of importance in combustion and atmospheric chemistry, stands out as a paradigm reaction involving triplet- and singlet-state potential energy surfaces (PESs) interconnected by intersystem crossing (ISC). This reaction poses challenges for theory and experiments owing to the ruggedness and high dimensionality of these potentials, as well as the long lifetimes of the collision complexes. Primary products from five competing channels (H + CH(2)CHO, H + CH(3)CO, H(2) + CH(2)CO, CH(3) + HCO, CH(2) + CH(2)O) and branching ratios (BRs) are determined in crossed molecular beam experiments with soft electron-ionization mass-spectrometric detection at a collision energy of 8.4 kcal/mol. As some of the observed products can only be formed via ISC from triplet to singlet PESs, from the product BRs the extent of ISC is inferred. A new full-dimensional PES for the triplet state as well as spin-orbit coupling to the singlet PES are reported, and roughly half a million surface hopping trajectories are run on the coupled singlet-triplet PESs to compare with the experimental BRs and differential cross-sections. Both theory and experiment find almost equal contributions from the two PESs to the reaction, posing the question of how important is it to consider the ISC as one of the nonadiabatic effects for this and similar systems involved in combustion chemistry. Detailed comparisons at the level of angular and translational energy distributions between theory and experiment are presented for the two primary channel products, CH(3) + HCO and H + CH(2)CHO. The agreement between experimental and theoretical functions is excellent, implying that theory has reached the capability of describing complex multichannel nonadiabatic reactions.