Computational Study on the Water Corrosion Process at Schreibersite (Fe2NiP) Surfaces: from Phosphide to Phosphates

Phosphorus (P) is a fundamental element for whatever form of life, in the same way as the other biogenic macroelements (SONCH). The prebiotic origin of P is still a matter of debate, as the phosphates present on earth are trapped in almost insoluble solid matrixes (apatites) and, therefore, hardly available for inclusion in living systems in the prebiotic era. The most accepted theories regard a possible exogenous origin during the Archean Era, through the meteoritic bombardment, when tons of reactive P in the form of phosphide ((Fe,Ni)3P, schreibersite mineral) reached the primordial earth, reacting with water and providing oxygenated phosphorus compounds (including phosphates). In the last 20 years, laboratory experiments demonstrated that the corrosion process of schreibersite by water indeed leads to reactive phosphates that, in turn, react with other biological building blocks (nucleosides and simple sugars) to form more complex molecules (nucleotides and complex sugars). In the present paper, we study the water corrosion of different crystalline surfaces of schreibersite by means of periodic DFT (density functional theory) simulations. Our results show that water adsorbs molecularly on the most stable (110) surface but dissociates on the less stable (001) one, giving rise to further reactivity. Indeed, subsequent water adsorptions, up to the water monolayer coverage, show that, on the (001) surface, iron and nickel atoms are the first species undergoing the corrosion process and, in a second stage, the phosphorus atoms also get involved. When adsorbing up to three and four water molecules per unit cell, the most stable structures found are the phosphite and phosphate forms of phosphorus, respectively. Simulation of the vibrational spectra of the considered reaction products revealed that the experimental band at 2423 cm-1 attributed to the P-H stretching frequency is indeed predicted for a phosphite moiety attached to the schreibersite (001) surface upon chemisorption of up to three water molecules.

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 (CH₂CHCN) 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(³P) + CH₂CHCN 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 ³P and excited ¹D atomic oxygen. From product angular and time-of-flight distributions in crossed molecular beam experiments with mass spectrometric detection at a collision energy, E_c, 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(³P) reaction leads to two main product channels: (i) CH₂CNH (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 CH₂CHCN and a possible formation route of CH₂CNH in the interstellar medium. The O(¹D) + CH₂CHCN reaction mainly leads to the formation of CH₂CNH + CO adiabatically on the singlet PES. This result can improve models related to the chemistry of interstellar ice and cometary comas, where O(¹D) reactions can play a role. Overall, our results are expected to be useful for improving the models of combustion and extraterrestrial environments.

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

The reaction between the ground-state hydroxyl radical, OH(²Π), and ethylene, C₂H₄, 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-CH₂CHOH (vinyl alcohol) + H, CH₃CHO (acetaldehyde) + H, and H₂CO (formaldehyde) + CH₃ 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-CH₂CHOH + H product channels account for 38% (in similar amounts) of the addition mechanism yield, the H₂CO + CH₃ channel for ∼58%, while the CH₃CHO + H channel is formed in negligible amount (<4%). The implications for combustion and astrochemical environments are discussed.

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 (HC_nN, 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 ³P and excited ¹D 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(³P, ¹D) with HCCCN(X¹Σ⁺) 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, E_c, 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(³P) 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 ¹HCCN (cyanomethylene) + CO. The O(¹D) 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(³P) + HC₃N at 300 K is in the 10^-12 cm³ 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.

Intersystem crossing in the entrance channel of the reaction of O(3P) with pyridine

Two quantum effects can enable reactions to take place at energies below the barrier separating reactants from products: tunnelling and intersystem crossing between coupled potential energy surfaces. Here we show that intersystem crossing in the region between the pre-reactive complex and the reaction barrier can control the rate of bimolecular reactions for weakly coupled potential energy surfaces, even in the absence of heavy atoms. For O(3P) plus pyridine, a reaction relevant to combustion, astrochemistry and biochemistry, crossed-beam experiments indicate that the dominant products are pyrrole and CO, obtained through a spin-forbidden ring-contraction mechanism. The experimental findings are interpreted-by high-level quantum-chemical calculations and statistical non-adiabatic computations of branching fractions-in terms of an efficient intersystem crossing occurring before the high entrance barrier for O-atom addition to the N-atom lone pair. At low to moderate temperatures, the computed reaction rates prove to be dominated by intersystem crossing.

Reaction N(2D) + CH2CCH2 (Allene): An Experimental and Theoretical Investigation and Implications for the Photochemical Models of Titan

We report on a combined experimental and theoretical investigation of the N(²D) + CH₂CCH₂ (allene) reaction of relevance in the atmospheric chemistry of Titan. Experimentally, the reaction was investigated (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 (E _c) of 33 kJ/mol to determine the primary products and the reaction micromechanism and (ii) in a continuous supersonic flow reactor to determine the rate constant as a function of temperature from 50 to 296 K. Theoretically, electronic structure calculations of the doublet C₃H₄N 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(²D) to one of the two equivalent carbon-carbon double bonds of CH₂CCH₂, followed by the formation of several cyclic and linear isomeric C₃H₄N intermediates that can undergo unimolecular decomposition to bimolecular products with elimination of H, CH₃, HCN, HNC, and CN. The kinetic experiments confirm the barrierless nature of the reaction through the measurement of rate constants close to the gas-kinetic rate at all temperatures. Statistical estimates of product branching fractions (BFs) on the theoretical PES were carried out under the conditions of the CMB experiments at room temperature and at temperatures (94 and 175 K) relevant for Titan. Up to 14 competing product channels were statistically predicted with the main ones at E _c = 33 kJ/mol being formation of cyclic-CH₂C(N)CH + H (BF = 87.0%) followed by CHCCHNH + H (BF = 10.5%) and CH₂CCNH + H (BF = 1.4%) the other 11 possible channels being negligible (BFs ranging from 0 to 0.5%). BFs under the other conditions are essentially unchanged. Experimental dynamical information could only be obtained on the overall H-displacement channel, while other possible channels could not be confirmed within the sensitivity of the method. This is also in line with theoretical predictions as the other possible channels are predicted to be negligible, including the HCN/HNC + C₂H₃ (vinyl) channels (overall BF < 1%). The dynamics and product distributions are dramatically different with respect to those observed in the isomeric reaction N(²D) + CH₃CCH (propyne), where at a similar E _c the main product channels are CH₂NH (methanimine) + C₂H (BF = 41%), c-C(N)CH + CH₃ (BF = 32%), and CH₂CHCN (vinyl cyanide) + H (BF = 12%). Rate coefficients (the recommended value is 1.7 (±0.2) × 10^-10 cm³ s^-1 over the 50-300 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 abundance (including any new products formed) as a function of the altitude.

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(²D), with vinyl cyanide, CH₂CHCN, 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, E_c, 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 E_c of the CMB experiment. According to our theoretical calculations, the reaction is found to proceed via barrierless addition of N(²D) to the carbon–carbon double bond of CH₂=CH–CN, followed by the formation of cyclic and linear intermediates that can undergo H, CN, and HCN elimination. In competition, the N(²D) addition to the CN group is also possible via a submerged barrier, leading ultimately to N₂ + C₃H₃ 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 C₃H₂N₂ isomeric products. Experimentally, no evidence of CN, HCN, and N₂ forming channels was observed. These findings were corroborated by the theory, which predicts a variety of competing product channels, following N(²D) addition to the double bond, with the main ones, at E_c = 31.4 kJ/mol, being six isomeric H forming channels: c-CH(N)CHCN + H (BF = 35.0%), c-CHNCHCN + H (BF = 28.1%), CH₂NCCN + H (BF = 26.3%), c-CH₂(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-CH₂(N)CH(2H-azirine) + CN and c-CH₂(N)C + HCN are predicted to be negligible (0.02% and 0.2%, respectively). The highly exothermic N₂ + CH₂CCH (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.

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, X²Σ⁺) and cyanoethene (C₂H₃CN) 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, CH₂C(CN)₂. 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.

Non-energetic Formation of Ethanol via CCH Reaction with Interstellar H2O Ices. A Computational Chemistry Study

Ethanol (CH₃CH₂OH) is a relatively common molecule, often found in star-forming regions. Recent studies suggest that it could be a parent molecule of several so-called interstellar complex organic molecules (iCOMs). However, the formation route of this species remains under debate. In the present work, we study the formation of ethanol through the reaction of CCH with one H₂O molecule belonging to the ice as a test case to investigate the viability of chemical reactions based on a "radical + ice component" scheme as an alternative mechanism for the synthesis of iCOMs, beyond the usual radical-radical coupling. This has been done by means of DFT calculations adopting two clusters of 18 and 33 water molecules as ice models. Results indicate that CH₃CH₂OH can potentially be formed by this proposed reaction mechanism. The reaction of CCH with H₂O on the water ice clusters can be barrierless (because of the help of boundary icy water molecules acting as proton-transfer assistants), leading to the formation of vinyl alcohol precursors (H₂CCOH and CHCHOH). Subsequent hydrogenation of vinyl alcohol yielding ethanol is the only step presenting a low activation energy barrier. We finally discuss the astrophysical implications of these findings.

Water Interaction with Fe2NiP Schreibersite (110) Surface: a Quantum Mechanical Atomistic Perspective

Phosphorus is an element of primary importance for all living creatures, being present in many biological activities in the form of phosphate (PO₄ ^3-). However, there are still open questions about the origin of this specific element and on the transformation that allowed it to be incorporated in biological systems. The most probable source of prebiotic phosphorus is the intense meteoritic bombardment during the Archean era, a few million years after the solar system formation, which brought tons of iron-phosphide materials (schreibersite) on the early Earth crust. It was recently demonstrated that by simple wetting/corrosion processes from this material, various oxygenated phosphorus compounds are produced. In the present work, the wetting process of schreibersite (Fe₂NiP) was studied by computer simulations using density functional theory, with the PBE functional supplemented with dispersive interactions through a posteriori empirical correction. To start disentangling the complexity of the system, only the most stable (110) surface of Fe₂NiP was used simulating different water coverages, from which structures, water binding energies, and vibrational spectra have been predicted. The computed (ana-)harmonic infrared spectra have been compared with the experimental ones, thus, confirming the validity of the adopted methodology and models.

Crossed-beam and theoretical studies of multichannel nonadiabatic reactions: branching fractions and role of intersystem crossing for O(3P)+1,3-butadiene

Atomic oxygen reactions can contribute significantly to the oxidation of unsaturated aliphatic and aromatic hydrocarbons. The reaction mechanism is started by electrophilic O atom addition to the unsaturated bond(s) to...

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(²D) with CH₃CCH (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 (E_c) 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(²D) to the carbon–carbon triple bond of CH₃CCH and an insertion of N(²D) into the CH bond of the methyl group, followed by the formation of cyclic and linear intermediates that can undergo H, CH₃, and C₂H 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 E_c = 31.0 kJ/mol, being the formation of CH₂NH (methanimine) + C₂H (ethylidyne) (BF = 0.41), c-C(N)CH + CH₃ (BF = 0.32), CH₂CHCN (acrylonitrile) + H (BF = 0.12), and c-CH₂C(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 E_c = 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 (C₂N₂, C₃N) 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(³P). During the past decade, we have performed in-depth investigations of the reactions of O(³P) 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 (C₆H₆) molecule. Despite extensive experimental and theoretical studies on the kinetics and dynamics of the O(³P) + C₆H₆ reaction, the relative importance of the C₆H₅O (phenoxy) + H open-shell products and of the spin-forbidden C₅H₆ (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(³P) + benzene at a collision energy (E_c) of 8.2 kcal/mol, focusing on the occurrence of the phenoxy + H and spin-forbidden C₅H₆ + 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(¹D) + benzene at the same E_c. 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(³P) + benzene reaction, the so-validated computational methodology has been used to predict (i) the BF between the C₆H₅O + H and C₅H₆ + 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 (C₅H₆ + 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.

Theoretical Study of the Extent of Intersystem Crossing in the O(3P) + C6H6 Reaction with Experimental Validation

The extent of intersystem crossing in the O(³P) + C₆H₆ reaction, a prototypical system for spin-forbidden reactions in oxygenated aromatic molecules, is theoretically evaluated for the first time. Calculations are performed using nonadiabatic transition-state theory coupled with stochastic master equation simulations and Landau-Zener theory. It is found that the dominant intersystem crossing pathways connect the T2 and S0 potential energy surfaces through at least two distinct minimum-energy crossing points. The calculated channel-specific rate constants and intersystem crossing branching fractions differ from previous literature estimates and provide valuable kinetic data for the investigation of benzene and polycyclic aromatic hydrocarbons oxidation in interstellar, atmospheric, and combustion conditions. The theoretical results are supported by crossed molecular beam experiments with electron ionization mass-spectrometric detection and time-of-flight analysis at 8.2 kcal/mol collision energy. This system is a suitable benchmark for theoretical and experimental studies of intersystem crossing in aromatic species.

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 SiO₂ samples, selected as a prototype of silicate minerals. The experiments were carried out after reaction of formamide at 160 °C on amorphous SiO₂ (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.

Angular Distribution of Ion Products in the Double Photoionization of Propylene Oxide

A photoelectron-photoion-photoion coincidence technique, using an ion imaging detector and tunable synchrotron radiation in the 18.0–37.0 eV photon energy range, inducing the ejection of molecular valence electrons, has been applied to study the double ionization of the propylene oxide, a simple prototype chiral molecule. The experiment performed at the Elettra Synchrotron Facility (Trieste, Italy) allowed to determine angular distributions for ions produced by the two-body dissociation reactions following the Coulomb explosion of the intermediate (C₃H₆O)²⁺ molecular dication. The analysis of the coincidence spectra recorded at different photon energies was done in order to determine the dependence of the β anisotropy parameter on the photon energy for the investigated two-body fragmentation channels. In particular, the reaction leading to CH3+ + C₂H₃O⁺ appears to be characterized by an increase of β, from β ≈ 0.00 up to β = 0.59, as the photon energy increases from 29.7 to 37.0 eV, respectively. This new observation confirms that the dissociation channel producing CH3+ and C₂H₃O⁺ final ions can occur with two different microscopic mechanisms as already indicated by the bimodality obtained in the kinetic energy released (KER) distributions as a function of the photon energy in a recent study. Energetic considerations suggest that experimental data are compatible with the formation of two different stable isomers of C₂H₃O⁺: acetyl and oxiranyl cations. These new experimental data are inherently relevant and are mandatory information for further experimental and theoretical investigations involving oriented chiral molecules and linearly or circularly polarized radiation. This work is in progress in our laboratory.

An Experimental and Theoretical Investigation of 1-Butanol Pyrolysis

Bioalcohols are a promising family of biofuels. Among them, 1-butanol has a strong potential as a substitute for petrol. In this manuscript, we report on a theoretical and experimental characterization of 1-butanol thermal decomposition, a very important process in the 1-butanol combustion at high temperatures. Advantage has been taken of a flash pyrolysis experimental set-up with mass spectrometric detection, in which the brief residence time of the pyrolyzing mixture inside a short, resistively heated SiC tube allows the identification of the primary products of the decomposing species, limiting secondary processes. Dedicated electronic structure calculations of the relevant potential energy surface have also been performed and RRKM estimates of the rate coefficients and product branching ratios up to 2,000 K are provided. Both electronic structure and RRKM calculations are in line with previous determinations. According to the present study, the H₂O elimination channel leading to 1-butene is more important than previously believed. In addition to that, we provide experimental evidence that butanal formation by H₂ elimination is not a primary decomposition route. Finally, we have experimental evidence of a small yield of the CH₃ elimination channel.

Combined Experimental-Theoretical Study of the OH + CO → H + CO2 Reaction Dynamics

A combined experimental-theoretical study is performed to advance our understanding of the dynamics of the prototypical tetra-atom, complex-forming reaction OH + CO → H + CO₂, which is also of great practical relevance in combustion, Earth's atmosphere, and, potentially, Mars's atmosphere and interstellar chemistry. New crossed molecular beam experiments with mass spectrometric detection are analyzed together with the results from previous experiments and compared with quasi-classical trajectory (QCT) calculations on a new, full-dimensional potential energy surface (PES). Comparisons between experiment and theory are carried out both in the center-of-mass and laboratory frames. Good agreement is found between experiment and theory, both for product angular and translational energy distributions, leading to the conclusion that the new PES is the most accurate at present in elucidating the dynamics of this fundamental reaction. Yet, small deviations between experiment and theory remain and are presumably attributable to the QCT treatment of the scattering dynamics.

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, NH₂CN, another interstellar species. The second most important channel is the one leading to methyl cyanamide, CH₃NHCN (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 cm³ molec^-1 s^-1, β = -1.80, γ = 7.79 K for the channel leading to NH₂CN + CH₃; α = 5.05 × 10^-13 cm³ molec^-1 s^-1, β = -1.82, γ = 7.93 K for the channel leading to CH₃NHCN + H.

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

Isomer-Specific Chemistry in the Propyne and Allene Reactions with Oxygen Atoms: CH3CH + CO versus CH2CH2 + CO Products

We report direct experimental and theoretical evidence that, under single-collision conditions, the dominant product channels of the O((3)P) + propyne and O((3)P) + allene isomeric reactions lead in both cases to CO formation, but the coproducts are singlet ethylidene ((1)CH3CH) and singlet ethylene (CH2CH2), respectively. These data, which settle a long-standing issue on whether ethylidene is actually formed in the O((3)P) + propyne reaction, suggest that formation of CO + alkylidene biradicals may be a common mechanism in O((3)P) + alkyne reactions, in contrast to formation of CO + alkene molecular products in the corresponding isomeric O((3)P) + diene reactions, either in combustion or other gaseous environments. These findings are of fundamental relevance and may have implications for improved combustion models. Moreover, we predict that the so far neglected (1)CH3CH + CO channel is among the main reaction routes also when the C3H4O singlet potential energy surface is accessed from the OH + C3H3 (propargyl) entrance channel, which are radical species playing a key role in many combustion systems.

Relevance of the Channel Leading to Formaldehyde + Triplet Ethylidene in the O((3)P) + Propene Reaction under Combustion Conditions

Comprehension of the detailed mechanism of O((3)P) + unsaturated hydrocarbon reactions is complicated by the existence of many possible channels and intersystem crossing (ISC) between triplet and singlet potential energy surfaces (PESs). We report synergic experimental/theoretical studies of the O((3)P) + propene reaction by combining crossed molecular beams experiments using mass spectrometric detection at 9.3 kcal/mol collision energy (Ec) with high-level ab initio electronic structure calculations of the triplet PES and RRKM/master equation computations of branching ratios (BRs) including ISC. At high Ec's and temperatures higher than 1000 K, main products are found to be formaldehyde (H2CO) and triplet ethylidene ((3)CH3CH) formed in a reaction channel that has never been identified or considered significant in previous kinetics studies at 300 K and that, as such, is not included in combustion kinetics models. Global and channel-specific rate constants were computed and are reported as a function of temperature and pressure. This study shows that BRs of multichannel reactions useful for combustion modeling cannot be extrapolated from room-temperature kinetics studies.

Quasiclassical Trajectory Calculations of the N((2)D) + H2O Reaction Elucidating the Formation Mechanism of HNO and HON Seen in Molecular Beam Experiments

The N((2)D) + H2O is a reaction with competitive product channels, passing through several intermediates. Dynamics of this reaction had been investigated by two of the present authors at two collision energies, Ec, using the crossed molecular beams mass spectrometric method ( Faraday Discuss. 2001 , 119 , 27 - 49 ). The complicated mechanism of this reaction and puzzling results encouraged us to investigate the reaction in a joint experimental/theoretical study. Quasiclassical trajectory (QCT) calculations on an ab initio potential energy surface describing all channels of the title reaction are done with a focus on the N/H exchange channels. Interesting results of QCT calculations, in very good agreement with experimental data, reveal subtle details of the reaction dynamics of the title reaction to HNO/HON + H exit channels by disentangling the different routes to formation of the two possible HNO/HON isomers and therefore assisting in a critical manner the derivation of the reaction mechanism. Results of the present study show that the nonstatistical HNOH intermediate governs exit channels; therefore, the HON channel is as important as that of HNO. The study also confirms that the H2 + NO molecular channel is negligible even though the barrier to its formation is calculated to be well below the reactant asymptote.

The fast C(3P) + CH3OH reaction as an efficient loss process for gas-phase interstellar methanol

The reaction between ground state atomic carbon and methanol is shown to be an efficient destruction mechanism for interstellar methanol.

Gas-phase synthesis of phenyl oxoborane (C6H5BO) via the reaction of boron monoxide with benzene

Organyl oxoboranes (RBO) are valuable reagents in organic synthesis due to their role in Suzuki coupling reactions. However, organyl oxoboranes (RBO) are only found in trimeric forms (RBO3) commonly known as boronic acids or boroxins; obtaining their monomers has proved a complex endeavor. Here, we demonstrate an oligomerization-free formation of organyl oxoborane (RBO) monomers in the gas phase by a radical substitution reaction under single-collision conditions in the gas phase. Using the cross molecular beams technique, phenyl oxoborane (C6H5BO) is formed through the reaction of boronyl radicals (BO) with benzene (C6H6). The reaction is indirect, initially forming a van der Waals complex that isomerizes below the energy of the reactants and eventually forming phenyl oxoborane by hydrogen emission in an overall exoergic radical-hydrogen atom exchange mechanism.

A crossed beam and ab initio investigation on the formation of boronyldiacetylene (HCCCC11BO; X1Σ+) via the reaction of the boron monoxide radical (11BO; X2Σ+) with diacetylene (C4H2; X1Σg(+))

The reaction dynamics of the boron monoxide radical ((11)BO; X(2)Σ(+)) with diacetylene (C4H2; X(1)Σg(+)) were investigated at a nominal collision energy of 17.5 kJ mol(-1) employing the crossed molecular beam technique and supported by ab initio and statistical (RRKM) calculations. The reaction is governed by indirect (complex forming) scattering dynamics with the boron monoxide radical adding with its boron atom to the carbon-carbon triple bond of the diacetylene molecule at one of the terminal carbon atoms without entrance barrier. This leads to a doublet radical intermediate (C4H2(11)BO), which undergoes unimolecular decomposition through hydrogen atom emission from the C1 carbon atom via a tight exit transition state located about 18 kJ mol(-1) above the separated products. This process forms the hitherto elusive boronyldiacetylene molecule (HCCCC(11)BO; X(1)Σ(+)) in a bimolecular gas phase reaction under single collision conditions. The overall reaction was determined to be exoergic by 62 kJ mol(-1). The reaction dynamics are compared to the isoelectronic diacetylene (C4H2; X(1)Σg(+))-cyano radical (CN; X(2)Σ(+)) system studied previously in our group. The characteristics of boronyl-diacetylene and the boronyldiacetylene molecule (HCCCC(11)BO; X(1)Σ(+)) as well as numerous intermediates are reported for the first time.