Pressure and Temperature Effects on Product Channels of the Propargyl (HC⋮CCH2) Combination Reaction and the Formation of the “First Ring”

Radical-Radical Reactions in Molecular Weight Growth: The Phenyl + Propargyl Reaction

The mechanism for hydrocarbon ring growth in sooting environments is still the subject of considerable debate. The reaction of phenyl radical (C₆H₅) with propargyl radical (H₂CCCH) provides an important prototype for radical-radical ring-growth pathways. We studied this reaction experimentally over the temperature range of 300-1000 K and pressure range of 4-10 Torr using time-resolved multiplexed photoionization mass spectrometry. We detect both the C₉H₈ and C₉H₇ + H product channels and report experimental isomer-resolved product branching fractions for the C₉H₈ product. We compare these experiments to theoretical kinetics predictions from a recently published study augmented by new calculations. These ab initio transition state theory-based master equation calculations employ high-quality potential energy surfaces, conventional transition state theory for the tight transition states, and direct CASPT2-based variable reaction coordinate transition state theory (VRC-TST) for the barrierless channels. At 300 K only the direct adducts from radical-radical addition are observed, with good agreement between experimental and theoretical branching fractions, supporting the VRC-TST calculations of the barrierless entrance channel. As the temperature is increased to 1000 K we observe two additional isomers, including indene, a two-ring polycyclic aromatic hydrocarbon, and a small amount of bimolecular products C₉H₇ + H. Our calculated branching fractions for the phenyl + propargyl reaction predict significantly less indene than observed experimentally. We present further calculations and experimental evidence that the most likely cause of this discrepancy is the contribution of H atom reactions, both H + indenyl (C₉H₇) recombination to indene and H-assisted isomerization that converts less stable C₉H₈ isomers into indene. Especially at low pressures typical of laboratory investigations, H-atom-assisted isomerization needs to be considered. Regardless, the experimental observation of indene demonstrates that the title reaction leads, either directly or indirectly, to the formation of the second ring in polycyclic aromatic hydrocarbons.

Excited state photochemically driven surface formation of benzene from acetylene ices on Pluto and in the outer solar system

NASA's New Horizons mission unveiled a diverse landscape of Pluto's surface with massive regions being neutral in color, while others like Cthulhu Macula range from golden-yellow to reddish comprising up to half of Pluto's carbon budget. Here, we demonstrate in laboratory experiments merged with electronic structure calculations that the photolysis of solid acetylene - the most abundant precipitate on Pluto's surface - by low energy ultraviolet photons efficiently synthesizes benzene and polycyclic aromatic hydrocarbons via excited state photochemistry thus providing critical molecular building blocks for the colored surface material. Since low energy photons deliver doses to Pluto's surface exceeding those from cosmic rays by six orders of magnitude, these processes may significantly contribute to the coloration of Pluto's surface and of hydrocarbon-covered surfaces of Solar System bodies such as Triton in general. This discovery critically enhances our perception of the distribution of aromatic molecules and carbon throughout our Solar System.

Product Detection of the CH(X2Π) Radical Reaction with Cyclopentadiene: A Novel Route to Benzene

The reaction of the methylidyne radical (CH(X²Π)) with cyclopentadiene (c-C₅H₆) is studied in the gas phase at 4 Torr and 373 K using a multiplexed photoionization mass spectrometer. Under multiple collision conditions, the dominant product channel observed is the formation of C₆H₆ + H. Fitting the photoionization spectrum using reference spectra allows for isomeric resolution of C₆H₆ isomers, where benzene is the largest contributor with a relative branching fraction of 90 (±5)%. Several other C₆H₆ isomers are found to have smaller contributions, including fulvene with a branching fraction of 8 (±5)%. Master Equation calculations for four different entrance channels on the C₆H₇ potential energy surface are performed to explore the competition between CH cycloaddition to a C═C bond vs CH insertion into C-H bonds of cyclopentadiene. Previous studies on CH addition to unsaturated hydrocarbons show little evidence for the C-H insertion pathway. The present computed branching fractions support benzene as the sole cyclic product from CH cycloaddition, whereas fulvene is the dominant product from two of the three pathways for CH insertion into the C-H bonds of cyclopentadiene. The combination of experiment with Master Equation calculations implies that insertion must account for ∼10 (±5)% of the overall CH + cyclopentadiene mechanism.

Gas-phase synthesis of benzene via the propargyl radical self-reaction

Polycyclic aromatic hydrocarbons (PAHs) have been invoked in fundamental molecular mass growth processes in our galaxy. We provide compelling evidence of the formation of the very first ringed aromatic and building block of PAHs-benzene-via the self-recombination of two resonantly stabilized propargyl (C₃H₃) radicals in dilute environments using isomer-selective synchrotron-based mass spectrometry coupled to theoretical calculations. Along with benzene, three other structural isomers (1,5-hexadiyne, fulvene, and 2-ethynyl-1,3-butadiene) and o-benzyne are detected, and their branching ratios are quantified experimentally and verified with the aid of computational fluid dynamics and kinetic simulations. These results uncover molecular growth pathways not only in interstellar, circumstellar, and solar systems environments but also in combustion systems, which help us gain a better understanding of the hydrocarbon chemistry of our universe.

Are the three hydroxyphenyl radical isomers created equal?--The role of the phenoxy radical

We have investigated the thermal decomposition of the three hydroxyphenyl radicals (˙C6H4OH) in a heated microtubular reactor. Intermediates and products were identified isomer-selectively applying photoion mass-selected threshold photoelectron spectroscopy with vacuum ultraviolet synchrotron radiation. Similarly to the phenoxy radical (C6H5-O˙), hydroxyphenyl decomposition yields cyclopentadienyl (c-C5H5) radicals in a decarbonylation reaction at elevated temperatures. This finding suggests that all hydroxyphenyl isomers first rearrange to form phenoxy species, which subsequently decarbonylate, a mechanism which we also investigate computationally. Meta- and para-radicals were selectively produced and spectroscopically detectable, whereas the ortho isomer could not be traced due to its fast rethermalization and rapid decomposition in the reactor. A smaller barrier to isomerization to phenoxy was found to be the reason for this observation. Since hydroxyphenyl species may be present under typical sooting conditions in flames, the resonantly stabilized cyclopentadienyl radical adds to the hydrocarbon pool and can contribute to the formation of polycyclic aromatic hydrocarbons, which are precursors in soot formation.

Experimental and theoretical studies into the formation of C4-C6 products in partially chlorinated hydrocarbon pyrolysis systems: a probabilistic approach to congener-specific yield predictions

This work presents a study of the pyrolytic formation of vinylacetylene and benzene congeners formed from chlorinated hydrocarbon precursors, a complex, multipath polymerization system formed in a monomer-rich environment. (Co-)pyrolyses of dichloro- and trichloroethylene yield a rich array of products, and assuming a single dominant underlying growth mechanism, this (on comparing expected and observed products) allows a number of potentially competing channels to C4 and C6 products to be ruled out. Poor congener/isomer descriptions rule out even-carbon radical routes, and the absence of C3 and C5 products rule out odd-carbon processes. Vinylidenes appear unable to describe the increased reactivity of acetylenes with chlorination noted in our experiments, leaving molecular acetylene dimerization processes and, in C6 systems, the closely related Diels-Alder cyclization as the likely reaction mechanism. The feasibility of these routes is further supported by ab initio calculations. However, some of the most persuasive evidence is provided by congener-specific yield predictions enabled by the construction of a probability tree analogue of kinetic modeling. This approach is relatively quick to construct, provides surprisingly accurate predictions, and may be a very useful tool in screening for important reaction channels in poorly understood congener- or isomer-rich reaction systems.

Molecular mechanisms in the pyrolysis of unsaturated chlorinated hydrocarbons: formation of benzene rings. 2. Experimental and kinetic modeling studies

The mechanism of formation of benzene rings during the pyrolysis of dichloro- and trichloroethylenes has been investigated by the method of laser powered homogeneous pyrolysis coupled with product analysis by gas chromatography. Additionally, selected (co)pyrolyses between the chlorinated ethylenes, CH2Cl2, C4Cl4, C4Cl6, and C2H2 have been performed to explicitly probe the roles of 2C3 and C4/C2 reaction pairs in aromatic growth. The presence of odd-carbon products in neat C4Cl6 pyrolyses indicates that 2C3 processes are operative in these systems; however, comparison with product yields from C2HCl3 suggests that C4/C2 processes dominate most other systems. This is further evidenced by an absence of C3 and other odd-carbon species in (co)pyrolyses with dichloromethane which should seed C3-based growth. The reactions of perchlorinated C4 species C4Cl5, C4Cl3, and C4Cl4 with C2Cl2 were subsequently explored through extensive kinetic simulations of the possible reaction pathways based on previous kinetic models and the exhaustive quantum chemical investigations of our preceding work. The experimental and theoretical results strongly suggest that, at moderate temperatures, aromatic ring formation from chlorinated ethylenes normally follows a Diels-Alder coupling of C4 and C2 molecular units followed by internal shifts; the one exception is the C4Cl4 + C2Cl2 system, where steric factors lead to the formation of nonaromatic products. There is little evidence for radical-based routes in these systems.

Product branching ratios in photodissociation of phenyl radical: a theoretical ab initio/Rice-Ramsperger-Kassel-Marcus study

Ab initio CCSD(T)/CBS//B3LYP/6-311G** calculations of the potential energy surface for possible dissociation channels of the phenyl radical are combined with microcanonical Rice-Ramsperger-Kassel-Marcus calculations of reaction rate constants in order to predict statistical product branching ratios in photodissociation of c-C(6)H(5) at various wavelengths. The results indicate that at 248 nm the photodissociation process is dominated by the production of ortho-benzyne via direct elimination of a hydrogen atom from the phenyl radical. At 193 nm, the statistical branching ratios are computed to be 63.4%, 21.1%, and 14.4% for the o-C(6)H(4) + H, l-C(6)H(4) ((Z)-hexa-3-ene-1,5-diyne) + H, and n-C(4)H(3) + C(2)H(2) products, respectively, in a contradiction with recent experimental measurements, which showed C(4)H(3) + C(2)H(2) as the major product. Although two lower energy pathways to the i-C(4)H(3) + C(2)H(2) products are identified, they appeared to be kinetically unfavorable and the computed statistical branching ratio of i-C(4)H(3) + C(2)H(2) does not exceed 1%. To explain the disagreement with experiment, we optimized conical intersections between the ground and the first excited electronic states of C(6)H(5) and, based on their structures and energies, suggested the following photodissociation mechanism at 193 nm: c-C(6)H(5) 1 → absorption of a photon → electronically excited 1 → internal conversion to the lowest excited state → conversion to the ground electronic state via conical intersections at CI-2 or CI-3 → non-statistical decay of the vibrationally excited radical favoring the formation of the n-C(4)H(3) + C(2)H(2) products. This scenario can be attained if the intramolecular vibrational redistribution in the CI-2 or CI-3 structures in the ground electronic state is slower than their dissociation to n-C(4)H(3) + C(2)H(2) driven by the dynamical preference.

Pressure and temperature dependence of the reaction of vinyl radical with ethylene

This work reports measurements of absolute rate coefficients and Rice-Ramsperger-Kassel-Marcus (RRKM) master equation simulations of the C2H3+C2H4 reaction. Direct kinetic studies were performed over a temperature range of 300-700 K and pressures of 20 and 133 mbar. Vinyl radicals (H2C=CH) were generated by laser photolysis of vinyl iodide (C2H3I) at 266 nm, and time-resolved absorption spectroscopy was used to probe vinyl radicals through absorption at 423.2 nm. Measurements at 20 mbar are in good agreement with previous determinations at higher temperature. A weighted three-parameter Arrhenius fit to the experimental rate constant at 133 mbar, with the temperature exponent fixed, gives k=(7+/-1)x10(-14) cm3 molecule(-1) s(-1) (T/298 K)2 exp[-(1430+/-70) K/T]. RRKM master equation simulations, based on G3 calculations of stationary points on the C4H7 potential energy surface, were carried out to predict rate coefficients and product branching fractions. The predicted branching to 1-methylallyl product is relatively small under the conditions of the present experiments but increases as the pressure is lowered. Analysis of end products of 248 nm photolysis of vinyl iodide/ethylene mixtures at total pressures between 27 and 933 mbar provides no direct evidence for participation of 1-methylallyl.

Ab initio G3-type/statistical theory study of the formation of indene in combustion flames. I. Pathways involving benzene and phenyl radical

Ab initio G3(MP2,CC)//B3LYP calculations of the potential energy surface (PES) for the formation of indene involving hydrocarbon species abundant in combustion, including benzene, phenyl, propargyl, and methyl radicals, and acetylene, have been performed to investigate the build-up of an additional cyclopenta moiety over the existing six-member aromatic ring. They were followed by statistical calculations of high-pressure-limit thermal rate constants in the temperature range of 300-3000 K for all reaction steps utilizing conventional Rice-Ramsperger-Kassel-Marcus (RRKM) and transition-state (TST) theories. The hydrogen abstraction acetylene addition (HACA) type mechanism, which involves the formation of benzyl radical followed by addition of acetylene, is shown to have low barriers (12-16 kcal/mol) and to be a viable candidate to account for indene formation in combustion flames, such as the 1,3-butadiene flame, where this mechanism was earlier suggested as the major indene formation route (Granata et al. Combust. Flame 2002, 131, 273). The mechanism of indene formation involving the addition of propargyl radical to benzene and rearrangements on the C9H9 PES is demonstrated to have higher barriers for all reaction steps as compared to an alternative pathway, which starts from the recombination of phenyl and propargyl radicals and then proceeds by activation of the C9H8 adducts by H abstraction or elimination followed by five-member ring closure in C9H7 and H addition to the 2-indenyl radical. The suggested pathways represent potentially important contributors to the formation of indene in combustion flames, and the computed rate constants can be utilized in kinetic simulations of the reaction mechanisms leading to indene and to higher cyclopentafused polycyclic aromatic hydrocarbons (CP-PAH).

Association rate constants for reactions between resonance-stabilized radicals: C3H3 + C3H3, C3H3 + C3H5, and C3H5 + C3H5

Reactions between resonance-stabilized radicals play an important role in combustion chemistry. The theoretical prediction of rate coefficients and product distributions for such reactions is complicated by the fact that the initial complex-formation steps and some dissociation steps are barrierless. In this paper direct variable reaction coordinate transition state theory (VRC-TST) is used to predict accurately the association rate constants for the self and cross reactions of propargyl and allyl radicals. For each reaction, a set of multifaceted dividing surfaces is used to account for the multiple possible addition channels. Because of their resonant nature the geometric relaxation of the radicals is important. Here, the effect of this relaxation is explicitly calculated with the UB3LYP/cc-pvdz method for each mutual orientation encountered in the configurational integrals over the transition state dividing surfaces. The final energies are obtained from CASPT2/cc-pvdz calculations with all pi-orbitals in the active space. Evaluations along the minimum energy path suggest that basis set corrections are negligible. The VRC-TST approach was also used to calculate the association rate constant and the corresponding number of states for the C(6)H(5) + H --> C(6)H(6) exit channel of the C(3)H(3) + C(3)H(3) reaction, which is also barrierless. For this reaction, the interaction energies were evaluated with the CASPT2(2e,2o)/cc-pvdz method and a 1-D correction is included on the basis of CAS+1+2+QC/aug-cc-pvtz calculations for the CH(3) + H reference system. For the C(3)H(3) + C(3)H(3) reaction, the VRC-TST results for the energy and angular momentum resolved numbers of states in the entrance channels and in the C(6)H(5) + H exit channel are incorporated in a master equation simulation to determine the temperature and pressure dependence of the phenomenological rate coefficients. The rate constants for the C(3)H(3) + C(3)H(3) and C(3)H(5) + C(3)H(5) self-reactions compare favorably with the available experimental data. To our knowledge there are no experimental rate data for the C(3)H(3) + C(3)H(5) reaction.

Shock tube pyrolysis of 1,2,4,5-hexatetraene

1,2,4,5-Hexatetraene (1245HT) is, according to theory, a key intermediate to benzene from propargyl radicals in a variety of flames; however, it has only been experimentally observed once in previous studies of the C3H3 + C3H3 reaction. To determine if it is indeed an intermediate to benzene formation, 1245HT was synthesized, via a Grignard reaction, and pyrolysized in a single-pulse shock tube at two nominal pressures of 22 and 40 bar over a temperature range from 540 to 1180 K. At temperatures T < 700 K, 1245HT converts efficiently to 3,4-dimethylenecyclobutene (34DMCB) with a rate constant of k = 10(10.16) x exp(-23.4 kcal/RT), which is in good agreement with the one calculated by Miller and Klippenstein. At higher temperatures, various C6H6 isomers were generated, which is consistent with theory and earlier experimental studies. Thus, the current work strongly supports the theory that 1245HT plays a bridging role in forming benzene from propargyl radicals. RRKM modeling of the current data set has also been carried out with the Miller-Klippenstein potential. It was found that the theory gives reasonably good predictions of the experimental observations of 1245HT, 1,5-hexadiyne (15HD), and 34DMCB in the current study and in our earlier studies of 15HD pyrolysis and propargyl recombination; however, there is considerable discrepancy between experiment and theory for the isomerization route of 1,2-hexadien-5-yne (12HD5Y) --> 2-ethynyl-1,3-butadiene (2E13BD) --> fulvene.

An Optimized Semidetailed Submechanism of Benzene Formation from Propargyl Recombination

The self-reaction of propargyl (C3H3) radicals has been widely suggested as one of the key routes forming benzene in a variety of aliphatic flames. Currently, in the majority of aromatic models, the C3H3 + C3H3 submechanism often contains one or two C6H6 isomers and a few global reaction steps, which do not adequately represent the actual recombination chemistry. Recent experimental and theoretical studies on the direct propargyl recombination and subsequent C6H6 isomerization have provided sufficient information to revisit and revise the C3H3 + C3H3 reaction submechanism. In the present work, a semidetailed kinetic model consisting of seven isomeric C6H6 species and 14 reaction steps was constructed based on the most recent potential energy surface for this system. The trial model was subjected to systemic optimization by use of a recently developed physically bounded Gauss-Newton (PGN) method against detailed species profiles of direct propargyl recombination and 1,5-hexadiyne (15HD) isomerization obtained from experiments at high temperatures in a shock tube and at low temperatures in a flow reactor, which were all measured at very high pressure (shock tube) or atmospheric (flow reactor) conditions. Predictions of the optimized model were in excellent agreement with all experimental measurements. The optimized C3H3 + C3H3 reaction subset was also tested for flame modeling. Two different aromatic chemistry models that incorporate benzene formation from propargyl radicals as a single step reaction were modified to include the complete submechanism for propargyl recombination. The updated models predict significant percentages of three isomeric species [2-ethynyl-1,3-butadiene (2E13BD), fulvene, and benzene] in premixed fuel-rich acetylene and ethylene flames, reflecting the observed flame structures.

Isomeric Product Distributions from the Self-Reaction of Propargyl Radicals

We have investigated the isomeric C6H6 product distributions of the self-reaction of propargyl (C3H3) radicals at two nominal pressures of 25 and 50 bar over the temperature range 720-1350 K. Experiments were performed using propargyl iodide as the radical precursor in a high-pressure single-pulse shock tube with a residence time of 1.6-2.0 ms. The relative yields of the C6H6 products are strongly temperature dependent, and the main products are 1,5-hexadiyne (15HD), 1,2-hexadiene-5-yne (12HD5Y), 3,4-dimethylenecyclobutene (34DMCB), 2-ethynyl-1,3-butadiene (2E13BD), fulvene, and benzene, with the minor products being cis- and trans-1,3-hexadiene-5-yne (13HD5Y). 1,2,4,5-Hexatetraene (1245HT) was observed below 750 K but the concentrations were too low to be quantified. The experimentally determined entry branching ratios are: 44% 15HD, 38% 12HD5Y, and 18% 1245HT, which is efficiently converted to 34DMCB. Following the initial recombination step, various C6H6 isomers are formed by thermal rearrangement. The experimentally observed concentrations for the C6H6 species are in good agreement with earlier experiments on 15HD thermal rearrangement.