Time- and Isomer-Resolved Measurements of Sequential Addition of Acetylene to the Propargyl Radical

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.

Comprehensive Kinetics on the C7H7 Potential Energy Surface under Combustion Conditions

The automated kinetics workflow code, KinBot, was used to explore and characterize the regions of the C₇H₇ potential energy surface that are relevant to combustion environments and especially soot inception. We first explored the lowest-energy region, which includes the benzyl, fulvenallene + H, and cyclopentadienyl + acetylene entry points. We then expanded the model to include two higher-energy entry points, vinylpropargyl + acetylene and vinylacetylene + propargyl. The automated search was able to uncover the pathways from the literature. In addition, three important new routes were discovered: a lower-energy pathway connecting benzyl with vinylcyclopentadienyl, a decomposition mechanism from benzyl that results in side-chain hydrogen atom loss to produce fulvenallene + H, and shorter and lower energy routes to the dimethylene-cyclopentenyl intermediates. We systematically reduced the extended model to a chemically relevant domain composed of 63 wells, 10 bimolecular products, 87 barriers, and 1 barrierless channel and constructed a master equation using the CCSD(T)-F12a/cc-pVTZ//ωB97X-D/6-311++G(d,p) level of theory to provide rate coefficients for chemical modeling. Our calculated rate coefficients show excellent agreement with measured ones. We also simulated concentration profiles and calculated branching fractions from the important entry points to provide an interpretation of this important chemical landscape.

Automated Reaction Kinetics of Gas-Phase Organic Species over Multiwell Potential Energy Surfaces

Automation of rate-coefficient calculations for gas-phase organic species became possible in recent years and has transformed how we explore these complicated systems computationally. Kinetics workflow tools bring rigor and speed and eliminate a large fraction of manual labor and related error sources. In this paper we give an overview of this quickly evolving field and illustrate, through five detailed examples, the capabilities of our own automated tool, KinBot. We bring examples from combustion and atmospheric chemistry of C-, H-, O-, and N-atom-containing species that are relevant to molecular weight growth and autoxidation processes. The examples shed light on the capabilities of automation and also highlight particular challenges associated with the various chemical systems that need to be addressed in future work.

A Conical Intersection Influences the Ground State Rearrangement of Fulvene to Benzene

The rearrangement of fulvene to benzene is believed to play an important role in the formation of soot during hydrocarbon combustion. Previous work has identified two possible mechanisms for the rearrangement─a unimolecular path and a hydrogen-atom-assisted, bimolecular path. Computational results to date have suggested that the unimolecular mechanism faces a barrier of about 74 kcal/mol, which makes it unable to compete with the bimolecular mechanism under typical combustion conditions. This computed barrier is about 10 kcal/mol higher than the experimental value, which is an unusually large discrepancy for modern electronic structure theory. In the present work, we have reinvestigated the unimolecular mechanism computationally, and we have found a second transition state that is approximately 10 kcal/mol lower in energy than the previously identified one and, therefore, in excellent agreement with the experimental value. The existence of two transition states for the same rearrangement arises because there is a conical intersection between the two lowest singlet states which occurs in the vicinity of the reaction coordinates. The two possible paths around the cone on the lower adiabatic surface give rise to the two distinct saddle points. The lower barrier for the unimolecular mechanism now makes it competitive with the bimolecular one, according to our calculations. In support of this conclusion, we have reanalyzed some previous experimental results on anisole pyrolysis, which leads to benzene as a significant product and have shown that the unimolecular and bimolecular mechanisms for fulvene → benzene must be occurring competitively in that system. Finally, we have identified that similar conical intersections arise during the isomerizations of benzofulvene and isobenzofulvene to naphthalene.

A molecular beam and computational study on the barrierless gas phase formation of (iso)quinoline in low temperature extraterrestrial environments

Despite remarkable progress toward the understanding of the formation pathways leading to polycyclic aromatic hydrocarbons (PAHs) in combustion systems and in deep space, the complex reaction pathways leading to nitrogen-substituted PAHs (NPAHs) at low temperatures of molecular clouds and hydrocarbon-rich, nitrogen-containing atmospheres of planets and their moons like Titan have remained largely obscure. Here, we demonstrate through laboratory experiments and computations that the simplest prototype of NPAHs - quinoline and isoquinoline (C₉H₇N) - can be synthesized via rapid and de-facto barrier-less reactions involving o-, m- and p-pyridinyl radicals (C₅H₄N˙) with vinylacetylene (C₄H₄) under low-temperature conditions.

Continuous Butadiyne Addition to Propargyl: A Radical-Efficient Pathway for Polycyclic Aromatic Hydrocarbons

Polycyclic aromatic hydrocarbons (PAHs) play a crucial role in soot inception, interstellar evolution, and nanomaterial synthesis. Although several mechanisms, such as hydrogen-abstraction acetylene/vinylacetylene addition, have previously been proposed, PAH formation and growth are not yet fully understood. We propose an alternate PAH growth mechanism wherein propargyl radical reacts with butadiyne to form larger radicals containing newly fused aromatic rings. Butadiyne is an important intermediate in hydrocarbon oxidation and carbon rich stars, while propargyl is one of the most important resonantly stabilized radicals that persists for long times. Our proposed mechanism is validated by quantum chemical calculations, elementary reaction experiments, laminar flame analysis, and kinetic modeling. Our findings challenge the conventional wisdom that radical site regeneration, being central to PAH growth, requires sequential hydrogen elimination and/or abstraction. In our proposed mechanism, PAH growth does not depend on abundant free radical consumption, and could, therefore, help explain carbonaceous nanoparticle coalescence in radical-deficient reaction environments.

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.

Heteroatom (N, P, As, Sb, Bi) Effects on the Resonance-Stabilized 2-, 3-, and 4-Picolyl Radicals

Important recent experimental studies have allowed the isomer-selective identification of the 2-, 3-, and 4-picolyl radicals. The picolyl radicals and their valence isoelectronic P, As, Sb, and Bi congeners are investigated here. For the three observed parent radicals, the theoretical ionization potentials agree with experiment to within 0.02 eV. Two rules are proposed for predicting vertical ionization potentials (E_VIE) and relative energies. The E_VIE values for these radicals will be higher when large percentages of the SOMO orbitals are distributed on the atoms with greater electronegativities. The cations of these systems were also studied along with the closed-shell methylpyridines and their P, As, Sb, and Bi analogs. The energies for the cationic species will lie lower when high percentages of π natural localized molecular orbitals occur on the more electronegative atoms. The structures of the 2- and 4-isomers strongly depend upon the heteroatoms, with the C-C linkages adopting a single-double alternating bond manner when the heteroatoms become heavier. The 3-isomers adopt roughly equal C-C bond distances with small changes from N to Bi.

Gas-Phase Formation of C5H6 Isomers via the Crossed Molecular Beam Reaction of the Methylidyne Radical (CH; X2Π) with 1,2-Butadiene (CH3CHCCH2; X1A')

The bimolecular gas-phase reaction of the methylidyne radical (CH; X²Π) with 1,2-butadiene (CH₂CCHCH₃; X¹A') was investigated at a collision energy of 20.6 kJ mol^-1 under single collision conditions. Combining our laboratory data with high-level electronic structure calculations, we reveal that this bimolecular reaction proceeds through the barrierless addition of the methylidyne radical to the carbon-carbon double bonds of 1,2-butadiene leading to doublet C₅H₇ intermediates. These collision adducts undergo a nonstatistical unimolecular decomposition through atomic hydrogen elimination to at least the cyclic 1-vinyl-cyclopropene (p5/p26), 1-methyl-3-methylenecyclopropene (p28), and 1,2-bis(methylene)cyclopropane (p29) in overall exoergic reactions. The barrierless nature of this bimolecular reaction suggests that these cyclic C₅H₆ isomers might be viable targets to be searched for in cold molecular clouds like TMC-1.

Combustion in the future: The importance of chemistry

Combustion involves chemical reactions that are often highly exothermic. Combustion systems utilize the energy of chemical compounds released during this reactive process for transportation, to generate electric power, or to provide heat for various applications. Chemistry and combustion are interlinked in several ways. The outcome of a combustion process in terms of its energy and material balance, regarding the delivery of useful work as well as the generation of harmful emissions, depends sensitively on the molecular nature of the respective fuel. The design of efficient, low-emission combustion processes in compliance with air quality and climate goals suggests a closer inspection of the molecular properties and reactions of conventional, bio-derived, and synthetic fuels. Information about flammability, reaction intensity, and potentially hazardous combustion by-products is important also for safety considerations. Moreover, some of the compounds that serve as fuels can assume important roles in chemical energy storage and conversion. Combustion processes can furthermore be used to synthesize materials with attractive properties. A systematic understanding of the combustion behavior thus demands chemical knowledge. Desirable information includes properties of the thermodynamic states before and after the combustion reactions and relevant details about the dynamic processes that occur during the reactive transformations from the fuel and oxidizer to the products under the given boundary conditions. Combustion systems can be described, tailored, and improved by taking chemical knowledge into account. Combining theory, experiment, model development, simulation, and a systematic analysis of uncertainties enables qualitative or even quantitative predictions for many combustion situations of practical relevance. This article can highlight only a few of the numerous investigations on chemical processes for combustion and combustion-related science and applications, with a main focus on gas-phase reaction systems. It attempts to provide a snapshot of recent progress and a guide to exciting opportunities that drive such research beyond fossil combustion.

On the mechanism of soot nucleation

New ideas and theoretical results offer a solution to soot particle inception following critical examination of prior proposals.

Role of ring-enlargement reactions in the formation of aromatic hydrocarbons

Ring-enlargement reactions can provide a fast route towards the formation of six-membered single-ring or polycyclic aromatic hydrocarbons (PAHs).

Multiphoton dissociation dynamics of the indenyl radical at 248 nm and 193 nm

Photofragment translational spectroscopy is used to investigate the unimolecular photodissociation of the indenyl radical (C9H7). C9H7 radicals are generated by photodetachment of C9H7 - anions and are dissociated at 248 nm (5.00 eV) and 193 nm (6.42 eV). The following product channels are definitively observed at both wavelengths: C2H2 + C7H5, C2H2 + C3H3 + C4H2, and C2H2 + C2H2 + C5H3. The three-body product channels are energetically inaccessible from single photon excitation at either dissociation wavelength. This observation, in combination with calculated dissociation rates and laser power studies, implies that all dissociation seen in this experiment occurs exclusively through multiphoton processes in which the initial C9H7 radical absorbs two photons sequentially prior to dissociation to two or three fragments. The corresponding translational energy distributions for each product channel peak well below the maximum available energy for two photons and exhibit similar behavior regardless of dissociation wavelength. These results suggest that all products are formed by internal conversion to the ground electronic state, followed by dissociation.

A combined experimental and computational study on the reaction dynamics of the 1-propynyl radical (CH3CC; X2A1) with ethylene (H2CCH2; X1A1g) and the formation of 1-penten-3-yne (CH2CHCCCH3; X1A')

The crossed molecular beam reactions of the 1-propynyl radical (CH₃CC; X²A₁) with ethylene (H₂CCH₂; X¹A_1g) and ethylene-d₄ (D₂CCD₂; X¹A_1g) were performed at collision energies of 31 kJ mol^-1 under single collision conditions. Combining our laboratory data with ab initio electronic structure and statistical Rice-Ramsperger-Kassel-Marcus (RRKM) calculations, we reveal that the reaction is initiated by the barrierless addition of the 1-propynyl radical to the π-electron density of the unsaturated hydrocarbon of ethylene leading to a doublet C₅H₇ intermediate(s) with a life time(s) longer than the rotation period(s). The reaction eventually produces 1-penten-3-yne (p1) plus a hydrogen atom with an overall reaction exoergicity of 111 ± 16 kJ mol^-1. About 35% of p1 originates from the initial collision complex followed by C-H bond rupture via a tight exit transition state located 22 kJ mol^-1 above the separated products. The collision complex (i1) can also undergo a [1,2] hydrogen atom shift to the CH₃CHCCCH₃ intermediate (i2) prior to a hydrogen atom release; RRKM calculations suggest that this pathway contributes to about 65% of p1. In higher density environments such as in combustion flames and circumstellar envelopes of carbon stars close to the central star, 1-penten-3-yne (p1) may eventually form the cyclopentadiene (c-C₅H₆) isomer via hydrogen atom assisted isomerization followed by hydrogen abstraction to the cyclopentadienyl radical (c-C₅H₅) as an important pathway to key precursors to polycyclic aromatic hydrocarbons (PAHs) and to carbonaceous nanoparticles.

Kinetics of the OH Radical Reaction with Fulvenallene from 298 to 450 K

Self-recombination and cross-reactions of large resonant stabilized hydrocarbon radicals such as fulvenallenyl (C7H5) are predicted to form polycyclic aromatic hydrocarbons in combustion and the interstellar medium. Although fulvenallenyl is likely to be present in these environments, large uncertainties remain about its formation mechanisms. We have investigated the formation of fulvenallenyl by reacting the OH radical with fulvenallene (C7H6) over the 298 to 450 K temperature range and at a pressure of 5 Torr (667 Pa). The reaction rate coefficient is found to be 8.8(±1.7) × 10(-12) cm(3) s(-1) at room temperature with a negative temperature dependence that can be fit from 298 to 450 K to k(T) = 8.8(±1.7) × 10(-12) (T/298 K)(-6.6(±1.1)) exp[-(8.72(±3.03) kJ mol(-1))/(R((1/T) - (1/298 K)))] cm(3) s(-1). The comparison of the experimental data with calculated abstraction rate coefficients suggests that over the experimental temperature range, association of the OH radical to fulvenallene plays a significant role likely leading to a low fulvenallenyl branching fraction.

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.