Formation of benzene in the interstellar medium

Reaction dynamics of S(3P) with 1,3-butadiene and isoprene: crossed-beam scattering, low-temperature flow experiments, and high-level electronic structure calculations

Sulfur atoms serve as key players in diverse chemical processes, from astrochemistry at very low temperature to combustion at high temperature. Building upon our prior findings, showing cyclization to thiophenes following the reaction of ground-state sulfur atoms with dienes, we here extend this investigation to include many additional reaction products, guided by detailed theoretical predictions. The outcomes highlight the complex formation of products during intersystem crossing (ISC) to the singlet surfaces. Here, we employed crossed-beam velocity map imaging and high-level ab initio methods to explore the reaction of S(3P) with 1,3-butadiene and isoprene under single-collision conditions and in low-temperature flows. For the butadiene reaction, our experimental results show the formation of thiophene via H2 loss, a 2H-thiophenyl radical through H loss, and thioketene through ethene loss at a slightly higher collision energy compared to previous observations. Complementary Chirped-Pulse Fourier-Transform mmWave spectroscopy (CP-FTmmW) measurements in a uniform flow confirmed the formation of thioketene in the reaction at 20 K. For the isoprene reaction, we observed analogous products along with the 2H-thiophenyl radical arising from methyl loss and C3H4S (loss of ethene or H2 + acetylene). CP-FTmmW detected the formation of thioformaldehyde via loss of 1,3-butadiene, again in the 20 K flow. Coupled-cluster calculations on the pathways found by the automated kinetic workflow code KinBot support these findings and indicate ISC to the singlet surface, leading to the generation of various long-lived intermediates, including 5-membered heterocycles.

Noncovalent Interactions Steer the Formation of Polycyclic Aromatic Hydrocarbons

Aromatic molecules play an important role in the chemistry of astronomical environments such as the cold interstellar medium (ISM) and (exo)planetary atmospheres. The observed abundances of (polycyclic) aromatic hydrocarbons such as benzonitrile and cyanonaphthalenes are, however, highly underestimated by astrochemical models. This demonstrates the need for more experimentally verified reaction pathways. The low-temperature ion-molecule reaction of benzonitrile•+ with acetylene is studied here using a multifaceted approach involving kinetics and spectroscopic probing of the reaction products. A fast radiative association reaction via an in situ experimentally observed prereactive complex shows the importance of noncovalent interactions in steering the pathway during cold ion-molecule reactions. Product structures of subsequent reactions are unambiguously identified using infrared action spectroscopy and reveal the formation of nitrogen-containing, linked bicyclic structures such as phenylpyridine•+ and benzo-N-pentalene+ structures. The results, contradicting earlier assumptions on the product structure, demonstrate the importance of spectroscopic probing of reaction products and emphasize the possible formation of linked bicyclic molecules and benzo-N-pentalene+ structures in astronomical environments.

Rapid heating rates define the volatile emission and regolith composition of (3200) Phaethon

Asteroid (3200) Phaethon experiences extreme solar radiant heating ( ~ 750 °C) during perihelion (0.14 au), leading to comet-like activity. The regolith composition and mechanism of volatile emission are unknown but key to understanding JAXA's DESTINY+ mission data (fly-by in 2029) and the fate of near-Sun asteroids more generally. By subjecting CM chondrite fragments to fast, open system, cyclic heating (2-20 °C/min), simulating conditions on Phaethon we demonstrate that rapid heating rates combine with the low permeability, resulting in reactions between volatile gases and decomposing minerals. The retention of S-bearing gas limits the thermal decomposition of Fe-sulphides, allowing these minerals to survive repeated heating cycles. Slow escape of S-bearing gases provides a mechanism for repeated gas release from a thermally processed surface and, therefore the comet-like activity without requiring surface renewal to expose fresh material each perihelion cycle. We predict Phaethon regolith is composed of olivine, Fe-sulphides, Ca-sulphates and hematite.

The isomer distribution of C6H6 products from the propargyl radical gas-phase recombination investigated by threshold-photoelectron spectroscopy

The resonance-stabilization of the propargyl radical (C3H3) makes it among the most important reactive intermediates in extreme environments and grants it a long enough lifetime to recombine in both terrestrial combustion media and cold molecular clouds in space. This makes the propargyl self-reaction a pivotal step in the formation of benzene, the first aromatic ring, to eventually lead to polycyclic aromatic hydrocarbons in a variety of environments. In this work, by producing propargyl radicals in a flow tube where propyne reacted with F atoms and probing the reaction products by mass-selected threshold-photoelectron spectroscopy (TPES), we identified eight C6H6 products in total, including benzene. On top of providing the first comprehensive measurements of the branching ratios of the eight identified C6H6 isomers in the propargyl self reaction products (4 mbar, 298 K conditions), this study also highlights the advantages and disadvantages of using isomer-selective TPES to identify and quantify reaction products.

Unconventional gas-phase synthesis of biphenyl and its atropisomeric methyl-substituted derivatives

The biphenyl molecule (C12H10) acts as a fundamental molecular backbone in the stereoselective synthesis of organic materials due to its inherent twist angle causing atropisomerism in substituted derivatives and in molecular mass growth processes in circumstellar environments and combustion systems. Here, we reveal an unconventional low-temperature phenylethynyl addition-cyclization-aromatization mechanism for the gas-phase preparation of biphenyl (C12H10) along with ortho-, meta-, and para-substituted methylbiphenyl (C13H12) derivatives through crossed molecular beams and computational studies providing compelling evidence on their formation via bimolecular gas-phase reactions of phenylethynyl radicals (C6H5CC, X2A1) with 1,3-butadiene-d6 (C4D6), isoprene (CH2C(CH3)CHCH2), and 1,3-pentadiene (CH2CHCHCHCH3). The dynamics involve de-facto barrierless phenylethynyl radical additions via submerged barriers followed by facile cyclization and hydrogen shift prior to hydrogen atom emission and aromatization to racemic mixtures (ortho, meta) of biphenyls in overall exoergic reactions. These findings not only challenge our current perception of biphenyls as high temperature markers in combustion systems and astrophysical environments, but also identify biphenyls as fundamental building blocks of complex polycyclic aromatic hydrocarbons (PAHs) such as coronene (C24H12) eventually leading to carbonaceous nanoparticles (soot, grains) in combustion systems and in deep space thus affording critical insight into the low-temperature hydrocarbon chemistry in our universe.

The Composition and Chemistry of Titan's Atmosphere

In this review I summarize the current state of knowledge about the composition of Titan's atmosphere and our current understanding of the suggested chemistry that leads to that observed composition. I begin with our present knowledge of the atmospheric composition, garnered from a variety of measurements including Cassini-Huygens, the Atacama Large Millimeter/submillimeter Array, and other ground- and space-based telescopes. This review focuses on the typical vertical profiles of gases at low latitudes rather than global and temporal variations. The main body of the review presents a chemical description of how complex molecules are believed to arise from simpler species, considering all known "stable" molecules-those that have been uniquely identified in the neutral atmosphere. The last section of the review is devoted to the gaps in our present knowledge of Titan's chemical composition and how further work may fill those gaps.

Elucidating the chemical dynamics of the elementary reactions of the 1-propynyl radical (CH3CC; X2A1) with 2-methylpropene ((CH3)2CCH2; X1A1)

Exploiting the crossed molecular beam technique, we studied the reaction of the 1-propynyl radical (CH3CC; X2A1) with 2-methylpropene (isobutylene; (CH3)2CCH2; X1A1) at a collision energy of 38 ± 3 kJ mol-1. The experimental results along with ab initio and statistical calculations revealed that the reaction has no entrance barrier and proceeds via indirect scattering dynamics involving C7H11 intermediates with lifetimes longer than their rotation period(s). The reaction is initiated by the addition of the 1-propynyl radical with its radical center to the π-electron density at the C1 and/or C2 position in 2-methylpropene. Further, the C7H11 intermediate formed from the C1 addition either emits atomic hydrogen or undergoes isomerization via [1,2-H] shift from the CH3 or CH2 group prior to atomic hydrogen loss preferentially leading to 1,2,4-trimethylvinylacetylene (2-methylhex-2-en-4-yne) as the dominant product. The molecular structures of the collisional complexes promote hydrogen atom loss channels. RRKM results show that hydrogen elimination channels dominate in this reaction, with a branching ratio exceeding 70%. Since the reaction of the 1-propynyl radical with 2-methylpropene has no entrance barrier, is exoergic, and all transition states involved are located below the energy of the separated reactants, bimolecular collisions are feasible to form trimethylsubstituted 1,3-enyne (p1) via a single collision event even at temperatures as low as 10 K prevailing in cold molecular clouds such as G+0.693. The formation of trimethylsubstituted vinylacetylene could serve as the starting point of fundamental molecular mass growth processes leading to di- and trimethylsubstituted naphthalenes via the HAVA mechanism.

Electron impact electronic excitation of benzene: Theory and experiment

We report experimental differential cross sections (DCSs) for electron impact excitation of bands I to V of benzene at incident energies of 10, 12.5, 15, and 20 eV. They are compared to calculations using the Schwinger multichannel method while accounting for up to 437 open channels. For intermediate scattering angles, the calculations reveal that the most intense band (V) emerges from surprisingly similar contributions from all its underlying states (despite some preference for the dipole-allowed transitions). They further shed light on intricate multichannel couplings between the states of bands I to V and higher-lying Rydberg states. In turn, the measurements support a vibronic coupling mechanism for excitation of bands II and IV and also show an unexpected forward peak in the spin-forbidden transition accounting for band III. Overall, there is decent agreement between theory and experiment at intermediate angles and at lower energies and in terms of the relative DCSs of the five bands. Discrepancies between the present and previous experiment regarding bands IV and V draw attention to the need of additional experimental investigations. We also report measured DCSs for vibrational excitation of combined C-H stretching modes.

Possible Roles of Transition Metal Cations in the Formation of Interstellar Benzene via Catalytic Acetylene Cyclotrimerization

Polycyclic aromatic hydrocarbons (PAHs) are ubiquitous interstellar molecules. However, the formation mechanisms of PAHs and even the simplest cyclic aromatic hydrocarbon, benzene, are not yet fully understood. Recently, we reported the statistical and dynamical properties in the reaction mechanism of Fe+-catalyzed acetylene cyclotrimerization, whereby three acetylene molecules are directly converted to benzene. In this study, we extended our previous work and explored the possible role of the complex of other 3d transition metal cations, TM+ (TM = Sc, Ti, Mn, Co, and Ni), as a catalyst in acetylene cyclotrimerization. Potential energy profiles for bare TM+-catalyst (TM = Sc and Ti), for TM+NC--catalyst (TM = Sc, Ti, Mn, Co, and Ni), and for TM+-(H2O)8-catalyst (TM = Sc and Ti) systems were obtained using quantum chemistry calculations, including the density functional theory levels. The calculation results show that the scandium and titanium cations act as efficient catalysts in acetylene cyclotrimerization and that reactants, which contain an isolated acetylene and (C2H2)2 bound to a bare (ligated) TM cation (TM = Sc and Ti), can be converted into a benzene-metal-cation product complex without an entrance barrier. We found that the number of electrons in the 3d orbitals of the transition metal cation significantly contributes to the catalytic efficiency in the acetylene cyclotrimerization process. On-the-fly Born-Oppenheimer molecular dynamics (BOMD) simulations of the Ti+-NC- and Ti+-(H2O)8 complexes were also performed to comprehensively understand the nuclear dynamics of the reactions. The computational results suggest that interstellar benzene can be produced via acetylene cyclotrimerization reactions catalyzed by transition metal cation complexes.

Prediction and Rationalization of Abundant C-N-H Molecules in Different Environments

The extreme chemical diversity of CmNnHk molecules is at the same time very important (central in organic chemistry) and difficult to rationalize in the sense that some molecules are abundant and easy to synthesize, while others are rare and difficult to make. Using the recently developed criteria of molecular "magicity", combined with evolutionary structure prediction and quantum-chemical calculations, we study these molecules in a wide range of compositions (0 ≤ m ≤ 13, 0 ≤ n ≤ 4, and 0 ≤ k ≤ 14). "Magic" molecules are defined as those that are lower in energy than any isochemical mixture of molecules with the nearest compositions. The predicted "magic" molecules are in good agreement with compounds found in versatile environments (interstellar and circumstellar media, Titan's lower atmosphere, and crude oil fractions) and in experimental chemistry. This work shows the predictive power of our approach, capable of predicting and explaining stable molecules in complex systems.

Gas-phase formation of the resonantly stabilized 1-indenyl (C9H7•) radical in the interstellar medium

The 1-indenyl (C9H7•) radical, a prototype aromatic and resonantly stabilized free radical carrying a six- and a five-membered ring, has emerged as a fundamental molecular building block of nonplanar polycyclic aromatic hydrocarbons (PAHs) and carbonaceous nanostructures in deep space and combustion systems. However, the underlying formation mechanisms have remained elusive. Here, we reveal an unconventional low-temperature gas-phase formation of 1-indenyl via barrierless ring annulation involving reactions of atomic carbon [C(3P)] with styrene (C6H5C2H3) and propargyl (C3H3•) with phenyl (C6H5•). Macroscopic environments like molecular clouds act as natural low-temperature laboratories, where rapid molecular mass growth to 1-indenyl and subsequently complex PAHs involving vinyl side-chained aromatics and aryl radicals can occur. These reactions may account for the formation of PAHs and their derivatives in the interstellar medium and carbonaceous chondrites and could close the gap of timescales of their production and destruction in our carbonaceous universe.

Evaluated electron scattering cross section dataset for gaseous benzene in the energy range 0.1-1000 eV

In this study, a complete and self-consistent cross section dataset for electron transport simulations through gaseous benzene in the energy range 0.1-1000 eV has been critically compiled. Its reliability has been evaluated through a joint experimental and computational procedure. To accomplish this, the compiled dataset has been used as input for event-by-event Monte Carlo simulations of the magnetically confined electron transport through gaseous benzene, and the simulated transmitted intensity has been compared with the experimental one for different incident energies and benzene gas pressures.

Interstellar Benzene Formation Mechanisms via Acetylene Cyclotrimerization Catalyzed by Fe+ Attached to Water Ice Clusters: Quantum Chemistry Calculation Study

Benzene is the simplest building block of polycyclic aromatic hydrocarbons and has previously been found in the interstellar medium. Several barrierless reaction mechanisms for interstellar benzene formation that may operate under low-temperature and low-pressure conditions in the gas phase have been proposed. In this work, we studied different mechanisms for interstellar benzene formation based on acetylene cyclotrimerization catalyzed by Fe⁺ bound to solid water clusters through quantum chemistry calculations. We found that benzene is formed via a single-step process with one transition state from the three acetylene molecules on the Fe⁺(H₂O)_n (n = 1, 8, 10, 12 and 18) cluster surface. Moreover, the obtained mechanisms differed from those of single-atom catalysis, in which benzene is sequentially formed via multiple steps.

Graph-Driven Reaction Discovery: Progress, Challenges, and Future Opportunities

Graph-based descriptors, such as bond-order matrices and adjacency matrices, offer a simple and compact way of categorizing molecular structures; furthermore, such descriptors can be readily used to catalog chemical reactions (i.e., bond-making and -breaking). As such, a number of graph-based methodologies have been developed with the goal of automating the process of generating chemical reaction network models describing the possible mechanistic chemistry in a given set of reactant species. Here, we outline the evolution of these graph-based reaction discovery schemes, with particular emphasis on more recent methods incorporating graph-based methods with semiempirical and ab initio electronic structure calculations, minimum-energy path refinements, and transition state searches. Using representative examples from homogeneous catalysis and interstellar chemistry, we highlight how these schemes increasingly act as "virtual reaction vessels" for interrogating mechanistic questions. Finally, we highlight where challenges remain, including issues of chemical accuracy and calculation speeds, as well as the inherent challenge of dealing with the vast size of accessible chemical reaction space.

A crossed molecular beams and computational study on the unusual reactivity of banana bonds of cyclopropane (c-C3H6; ) through insertion by ground state carbon atoms (C(3Pj))

The mechanism and chemical dynamics of the reaction of ground electronic state atomic carbon C(³P_j) with cyclopropane c-C₃H₆ have been explored by combining crossed molecular beams experiments with electronic structure calculations of the pertinent triplet C₄H₆ potential energy surface and statistical computations of product branching ratios under single-collision conditions. The experimental findings suggest that the reaction proceeds via indirect scattering dynamics through triplet C₄H₆ reaction intermediate(s) leading to C₄H₅ product(s) plus atomic hydrogen via a tight exit transition state, with the overall reaction exoergicity evaluated as 231 ± 52 kJ mol^-1. The calculations indicate that C(³P_j) can easily insert into one of the three equivalent C-C 'banana' bonds of cyclopropane overcoming a low barrier of only 2 kJ mol^-1 following the formation of a van der Waals reactant complex stabilized by 15 kJ mol^-1. The carbon atom insertion into one of the six C-H bonds is also feasible via a slightly higher barrier of 5 kJ mol^-1. These results highlight an unusual reactivity of cyclopropane's banana C-C bonds, which behave more like unsaturated C-C bonds with a π-character than saturated σ C-C bonds, which are known to be generally unreactive toward the ground electronic state atomic carbon such as in ethane (C₂H₆). The statistical theory predicts the overall product branching ratios at the experimental collision energy as 50% for 1-butyn-4-yl, 33% for 1,3-butadien-2-yl, i-C₄H₅, and 11% for 1,3-butadien-1-yl, n-C₄H₅, with i-C₄H₅ (230 kJ mol^-1 below the reactants) favored by the C-C insertion providing the best match with the experimentally observed reaction exoergicity. The C(³P_j) + c-C₃H₆ reaction is predicted to be a source of C₄H₅ radicals under the conditions where its low entrance barriers can be overcome, such as in planetary atmospheres or in circumstellar envelopes but not in cold molecular clouds. Both i- and n-C₄H₅ can further react with acetylene eventually producing the first aromatic ring and hence, the reaction of the atomic carbon with c-C₃H₆ can be considered as an initial step toward the formation of benzene.

Triple ionization and fragmentation of benzene trimers following ultrafast intermolecular Coulombic decay

Intermolecular interactions involving aromatic rings are ubiquitous in biochemistry and they govern the properties of many organic materials. Nevertheless, our understanding of the structures and dynamics of aromatic clusters remains incomplete, in particular for systems beyond the dimers, despite their high presence in many macromolecular systems such as DNA and proteins. Here, we study the fragmentation dynamics of benzene trimer that represents a prototype of higher-order aromatic clusters. The trimers are initially ionized by electron-collision with the creation of a deep-lying carbon 2s⁻¹ state or one outer-valence and one inner-valence vacancies at two separate molecules. The system can thus relax via ultrafast intermolecular decay mechanisms, leading to the formation of C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${}_{6}{{{{{{{{{\rm{H}}}}}}}}}_{6}}^{+}\cdot$$\end{document}6H6+⋅C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${}_{6}{{{{{{{{{\rm{H}}}}}}}}}_{6}}^{+}\cdot$$\end{document}6H6+⋅C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${}_{6}{{{{{{{{{\rm{H}}}}}}}}}_{6}}^{+}$$\end{document}6H6+ trications and followed by a concerted three-body Coulomb explosion. Triple-coincidence ion momentum spectroscopy, accompanied by ab-initio calculations and further supported by strong-field laser experiments, allows us to elucidate the details on the fragmentation dynamics of benzene trimers.

Higher-order aromatic clusters are prevalent in biochemical systems, but a full understanding of their structural and dynamical properties is lacking. Here, the authors demonstrate that inner-valence ionization can induce ultrafast relaxation and further fragmentation mechanisms in benzene trimers.

Gas-phase synthesis of racemic helicenes and their potential role in the enantiomeric enrichment of sugars and amino acids in meteorites

The molecular origins of homochirality on Earth is not understood well, particularly how enantiomerically enriched molecules of astrobiological significance like sugars and amino acids might have been synthesized on icy grains in space preceding their delivery to Earth. Polycyclic aromatic hydrocarbons (PAHs) identified in carbonaceous chondrites could have been processed in molecular clouds by circularly polarized light prior to the depletion of enantiomerically enriched helicenes onto carbonaceous grains resulting in chiral islands. However, the fundamental low temperature reaction mechanisms leading to racemic helicenes are still unknown. Here, by exploiting synchrotron based molecular beam photoionization mass spectrometry combined with electronic structure calculations, we provide compelling testimony on barrierless, low temperature pathways leading to racemates of [5] and [6]helicene. Astrochemical modeling advocates that gas-phase reactions in molecular clouds lead to racemates of helicenes suggesting a pathway for future astronomical observation and providing a fundamental understanding for the origin of homochirality on early Earth.

Deuterium Isotope Fractionation of Polycyclic Aromatic Hydrocarbons in Meteorites as an Indicator of Interstellar/Protosolar Processing History

The stable isotope composition of soluble and insoluble organic compounds in carbonaceous chondrites can be used to determine the provenance of organic molecules in space. Deuterium enrichment in meteoritic organics could be a residual signal of synthetic reactions occurring in the cold interstellar medium or an indicator of hydrothermal parent-body reactions. δD values have been measured in grains and bulk samples for a wide range of meteorites; however, these reservoirs are highly variable and may have experienced fractionation during thermal and/or aqueous alteration. Among the plethora of organic compounds in meteorites are polycyclic aromatic hydrocarbons (PAHs), which are stable and abundant in carbonaceous chondrites, and their δD ratio may preserve evidence about their formation environment as well as the influence of parent-body processes. This study tests hypotheses about the potential links between PAHs-deuteration concentrations and their formation conditions by examining the δD ratio of PAHs in three CM carbonaceous chondrites representing an aqueous alteration gradient. We use deuterium enrichments in soluble 2–5-ring PAHs as an indicator of either photon-driven deuteration due to unimolecular photodissociation in warm regions of space, gas-phase ion–molecule reactions in cold interstellar regions of space, or UV photolysis in ices. We also test hypothesized reaction pathways during parent-body processing that differ between partially and fully aromatized PAHs. New methodological approaches were developed to extract small, volatile PAHs without fractionation. Our results suggest that meteoritic PAHs could have formed through reactions in cold regions, with possible overprinting of deuterium enrichment during aqueous parent-body alteration, but the data could not rule out PAH alteration in icy mantles as well.

"Magic" Molecules and a New Look at Chemical Diversity of Hydrocarbons

We address the question why among the multitude of imaginable C_nH_m compositions some are easily synthesizable and abundant in nature, while others are not. To shed light on this problem we borrow approaches from nanocluster study, where stability with respect to neighboring compositions is used as a criterion of "magic" (particularly stable) clusters. By merging this criterion with predictions of lowest-energy structures of all C_nH_m molecules in a wide range of compositions (n ≤ 20, m ≤ 42) we provide guidelines for predicting the presence or absence of certain hydrocarbon molecules in various environments, their relative abundance and reactivity/inertness. The resulting maps of stability show the increased stability of C_2nH₂ compounds, polyaromatic hydrocarbons, and diamondoids, which is supported by experimental studies of the interstellar medium, flames, and petroleum fractions. This approach can be applied to any other molecular system, rationalizing the diversity of known compounds and predicting new potentially synthesizable molecules.

Formation of the Elusive Silylenemethyl Radical (HCSiH2; X2B2) via the Unimolecular Decomposition of Triplet Silaethylene (H2CSiH2; a3A″)

We investigated the formation of small organosilicon molecules─potential precursors to silicon-carbide dust grains ejected by dying carbon-rich asymptotic giant branch stars─in the gas phase via the reaction of atomic carbon (C) in its ³P electronic ground state with silane (SiH₄; X¹A₁) using the crossed molecular beams technique. The reactants collided under single collision conditions at a collision energy of 13.0 ± 0.2 kJ mol^-1, leading to the formation of the silylenemethyl radical (HCSiH₂; X²B₂) via the unimolecular decomposition of triplet silaethylene (H₂CSiH₂; a³A″). The silaethylene radical was formed via hydrogen migration of the triplet silylmethylene (HCSiH₃; X³A″) radical, which in turn was identified as the initial collision complex accessed via the barrierless insertion of atomic carbon into the silicon-hydrogen bond of silane. Our results mark the first observation of the silylenemethyl radical, where previously only its thermodynamically more stable methylsilylidyne (CH₃Si; X²A″) and methylenesilyl (CH₂SiH; X²A') isomers were observed in low-temperature matrices. Considering the abundance of silane and the availability of atomic carbon in carbon-rich circumstellar environments, our results suggest that future astrochemical models should be updated to include contributions from small saturated organosilicon molecules as potential precursors to pure gaseous silicon-carbides and ultimately to silicon-carbide dust.

Hexachlorobenzene-negative ion formation in electron attachment experiments

In this contribution, we report a comprehensive study on the formation of hexachlorobenzene (C₆Cl₆) negative ions probed by low-energy electron interactions from 0 up to 12 eV in a gas-phase crossed beam experiment. The anionic yields as a function of the electron energy reveal a rich fragmentation pattern of the dissociative electron attachment process, yet the most intense ion has been assigned to the non-dissociated parent anion that survives long enough within the detection time window. Other less intense fragment anions have been assigned as Cl^-, Cl₂^-, C₆Cl₄^-, and C₆Cl₅^-. The experimental results are accompanied by quantum chemical calculations at various levels of accuracy, providing an insight into the electronic structure, thermochemical thresholds, electron affinities and structures of neutral and anionic molecular species. The electron attachment process induces a considerable geometry change in the temporary-negative ion relative to the neutral molecule, where the most intense fragment anion assigned to Cl^- can be formed solely through a curve crossing involving a π*/σ* coupling. The yield of chlorine anions shows a signature of vibrational excitation reminiscent of a Jahn-Teller distortion.

Rotational and Vibrational Spectra of the Pyridyl Radicals: A Coupled-Cluster Study

Pyridyl is a prototypical nitrogen-containing aromatic radical that may be a key intermediate in the formation of nitrogen-containing aromatic molecules under astrophysical conditions. On meteorites, a variety of complex molecules with nitrogen-containing rings have been detected with nonterrestrial isotopic abundances, and larger nitrogen-containing polycyclic aromatic hydrocarbons (PANHs) have been proposed to be responsible for certain unidentified infrared emission bands in the interstellar medium. In this work, the three isomers of pyridyl (2-, 3-, and 4-pyridyl) have been investigated with coupled cluster methods. For each species, structures were optimized at the CCSD(T)/cc-pwCVTZ level of theory and force fields were calculated at the CCSD(T)/ANO0 level of theory. Second-order vibrational perturbation theory (VPT2) was used to derive anharmonic vibrational frequencies and vibrationally corrected rotational constants, and resonances among vibrational states below 3500 cm^-1 were treated variationally with the VPT2+K method. The results yield a complete set of spectroscopic parameters needed to simulate the pure rotational spectrum of each isomer, including electron-spin, spin-spin, and nuclear hyperfine interactions, and the calculated hyperfine parameters agree well with the limited available data from electron paramagnetic resonance spectroscopy. For the handful of experimentally measured vibrational frequencies determined from photoelectron spectroscopy and matrix isolation spectroscopy, the typical agreement is comparable to experimental uncertainty. The predicted parameters for rotational spectroscopy reported here can guide new experimental investigations into the yet-unobserved rotational spectra of these radicals.

Gas-Phase Study of the Elementary Reaction of the D1-Ethynyl Radical (C2D; X2Σ+) with Propylene (C3H6; X1A') under Single-Collision Conditions

The bimolecular gas-phase reactions of the D1-ethynyl radical (C₂D; X²Σ⁺) with propylene (C₃H₆; X¹A') and partially substituted D3-3,3,3-propylene (C₂H₃CD₃; X¹A') were studied under single collision conditions utilizing the crossed molecular beams technique. Combining our laboratory data with electronic structure and statistical calculations, the D1-ethynyl radical is found to add without barrier to the C1 and C2 carbons of the propylene reactant, resulting in doublet C₅H₆D intermediate(s) with lifetime(s) longer than their rotational period(s). These intermediates undergo isomerization and unimolecular decomposition via atomic hydrogen loss through tight exit transition states forming predominantly cis/trans-3-penten-1-yne ((HCC)CH═CH(CH₃)) and, to a minor amount, 3-methyl-3-buten-1-yne ((HCC)C(CH₃)═CH₂) via overall exoergic reactions. Although the title reaction does not lead to the cyclopentadiene molecule (c-C₅H₆, X¹A₁), high-temperature environments can convert the identified acyclic C₅H₆ isomers through hydrogen atom assisted isomerization to cyclopentadiene (c-C₅H₆, X¹A₁). Since both the ethynyl radical and propylene reactants have been observed in cold interstellar environments such as TMC-1 and the reaction is exoergic and all barriers lie below the energy of the separated reactants, these C₅H₆ product isomers are predicted to form in those low-temperature regions.

Formation and growth mechanisms of polycyclic aromatic hydrocarbons: A mini-review

Polycyclic aromatic hydrocarbons (PAHs) are mostly formed during the incomplete combustion of organic materials, but their importance and presence in materials science, and astrochemistry has also been proven. These carcinogenic persistent organic pollutants are essential in the formation of combustion generated particles as well. Due to their significant impact on the environment and human health, to understand the formation and growth of PAHs is essential. Therefore, the most important growth mechanisms are reviewed, and presented here from the past four decades (1981-2021) to initiate discussions from a new perspective. Although, the collected and analyzed observations are derived from both experimental, and computational studies, it is neither a systematic nor a comprehensive review. Nevertheless, the mechanisms were divided into three main categories, acetylene additions (e.g. HACA), vinylacetylene additions (HAVA), and radical reactions, and discussed accordingly.

Single photon double and triple ionization of allene

Double and triple ionization of allene are investigated using electron-electron, ion-ion, electron-electron-ion and electron-electron-ion-ion (ee, ii, eei, eeii) coincidence spectroscopies at selected photon energies. The results provide supporting evidence for a previously proposed roaming mechanism in H₃⁺ formation by double ionization. The lowest vertical double ionization energy is found to be 27.9 eV, while adiabatic double ionization is not accessed by vertical ionization at the neutral geometry. The triple ionization energy is found to be close to 50 eV in agreement with theoretical predictions. The doubly charged parent ion is stable up to about 2 eV above the threshold, after which dissociations by charge separation and by double charge retention occur with comparable intensities. Fragmentation to H⁺ + C₃H₃⁺ starts immediately above the threshold as a slow (metastable) decay with 130.5 ± 9.9 ns mean lifetime.

Halide–propene complexes: validated DSD-PBEP86-D3BJ calculations and photoelectron spectroscopy

Anion photoelectron spectroscopy has been used to determine the electron binding energies of the X−⋯C3H6 (X = Cl, Br, I) complexes.

Directed gas phase preparation of ethynylallene (H2CCCHCCH; X1A′) via the crossed molecular beam reaction of the methylidyne radical (CH; X2Π) with vinylacetylene (H2CCHCCH; X1A′)

The elementary reaction of the methylidyne radical with vinylacetylene leading to the predominant formation of ethynylallene and atomic hydrogen via indirect scattering dynamics.

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.

The Effect of Cluster Size on the Intra-Cluster Ionic Polymerization Process

Polyaromatic hydrocarbons (PAHs) are widespread in the interstellar medium (ISM). The abundance and relevance of PAHs call for a clear understanding of their formation mechanisms, which, to date, have not been completely deciphered. Of particular interest is the formation of benzene, the basic building block of PAHs. It has been shown that the ionization of neutral clusters can lead to an intra-cluster ionic polymerization process that results in molecular growth. Ab-initio molecular dynamics (AIMD) studies in clusters consisting of 3-6 units of acetylene modeling ionization events under ISM conditions have shown maximum aggregation of three acetylene molecules forming bonded C₆H₆⁺ species; the larger the number of acetylene molecules, the higher the production of C₆H₆⁺. These results lead to the question of whether clusters larger than those studied thus far promote aggregation beyond three acetylene units and whether larger clusters can result in higher C₆H₆⁺ production. In this study, we report results from AIMD simulations modeling the ionization of 10 and 20 acetylene clusters. The simulations show aggregation of up to four acetylene units producing bonded C₈H₈⁺. Interestingly, C₈H₈⁺ bicyclic species were identified, setting a precedent for their astrochemical identification. Comparable reactivity rates were shown with 10 and 20 acetylene clusters.

Pure hydrocarbon cycles in TMC-1: Discovery of ethynyl cyclopropenylidene, cyclopentadiene and indene

We report the detection for the first time in space of three new pure hydrocarbon cycles in TMC-1: c-C3HCCH (ethynyl cyclopropenylidene), c-C5H6 (cyclopentadiene) and c-C9H8 (indene). We derive a column density of 3.1 × 1011 cm-2 for the former cycle and similar values, in the range (1-2) × 1013 cm-2, for the two latter molecules. This means that cyclopentadiene and indene, in spite of their large size, are exceptionally abundant, only a factor of five less abundant than the ubiquitous cyclic hydrocarbon c-C3H2. The high abundance found for these two hydrocarbon cycles, together with the high abundance previously found for the propargyl radical (CH2CCH) and other hydrocarbons like vinyl and allenyl acetylene (Agúndez et al. 2021; Cernicharo et al. 2021a,b), start to allow us to quantify the abundant content of hydrocarbon rings in cold dark clouds and to identify the intermediate species that are probably behind the in situ bottom-up synthesis of aromatic cycles in these environments. While c-C3HCCH is most likely formed through the reaction between the radical CCH and c-C3H2, the high observed abundances of cyclopentadiene and indene are difficult to explain through currently proposed chemical mechanisms. Further studies are needed to identify how are five- and six-membered rings formed under the cold conditions of clouds like TMC-1.

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.

Detection of two interstellar polycyclic aromatic hydrocarbons via spectral matched filtering

Unidentified infrared emission bands are ubiquitous in many astronomical sources. These bands are widely, if not unanimously, attributed to collective emissions from polycyclic aromatic hydrocarbon (PAH) molecules, yet no single species of this class has been identified in space. Using spectral matched filtering of radio data from the Green Bank Telescope, we detected two nitrile-group-functionalized PAHs, 1- and 2-cyanonaphthalene, in the interstellar medium. Both bicyclic ring molecules were observed in the TMC-1 molecular cloud. In this paper, we discuss potential in situ gas-phase PAH formation pathways from smaller organic precursor molecules.

Reaction mechanisms of C(3PJ) and C+(2PJ) with benzene in the interstellar medium from quantum mechanical molecular dynamics simulations

While spectroscopic data on small hydrocarbons in interstellar media in combination with crossed molecular beam (CMB) experiments have provided a wealth of information on astrochemically relevant species, much of the underlying mechanistic pathways of their formation remain elusive. Therefore, in this work, the chemical reaction mechanisms of C(3PJ) + C6H6 and C+(2P) + C6H6 systems using the quantum mechanical molecular dynamics (QMMD) technique at the PBE0-D3(BJ) level of theory is investigated, mimicking a CMB experiment. Both the dynamics of the reactions as well as the electronic structure for the purpose of the reaction network are evaluated. The method is validated for the first reaction by comparison to the available experimental data. The reaction scheme for the C(3PJ) + C6H6 system covers the literature data, e.g. the major products are the 1,2-didehydrocycloheptatrienyl radical (C7H5) and benzocyclopropenyl radical (C6H5-CH), and it reveals the existence of less common pathways for the first time. The chemistry of the C+(2PJ) + C6H6 system is found to be much richer, and we have found that this is because of more exothermic reactions in this system in comparison to those in the C(3PJ) + C6H6 system. Moreover, using the QMMD simulation, a number of reaction paths have been revealed that produce three distinct classes of reaction products with different ring sizes. All in all, at all the collision energies and orientations, the major product is the heptagon molecular ion for the ionic system. It is also revealed that the collision orientation has a dominant effect on the reaction products in both systems, while the collision energy mostly affects the charged system. These simulations both prove the applicability of this approach to simulate crossed molecular beams, and provide fundamental information on reactions relevant for the interstellar medium.

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.

Low-temperature gas-phase formation of indene in the interstellar medium

Polycyclic aromatic hydrocarbons (PAHs) are fundamental molecular building blocks of fullerenes and carbonaceous nanostructures in the interstellar medium and in combustion systems. However, an understanding of the formation of aromatic molecules carrying five-membered rings-the essential building block of nonplanar PAHs-is still in its infancy. Exploiting crossed molecular beam experiments augmented by electronic structure calculations and astrochemical modeling, we reveal an unusual pathway leading to the formation of indene (C₉H₈)-the prototype aromatic molecule with a five-membered ring-via a barrierless bimolecular reaction involving the simplest organic radical-methylidyne (CH)-and styrene (C₆H₅C₂H₃) through the hitherto elusive methylidyne addition-cyclization-aromatization (MACA) mechanism. Through extensive structural reorganization of the carbon backbone, the incorporation of a five-membered ring may eventually lead to three-dimensional PAHs such as corannulene (C₂₀H₁₀) along with fullerenes (C₆₀, C₇₀), thus offering a new concept on the low-temperature chemistry of carbon in our galaxy.

Gas-Phase Formation of 1-Methylcyclopropene and 3-Methylcyclopropene via the Reaction of the Methylidyne Radical (CH; X2Π) with Propylene (CH3CHCH2; X1A')

The crossed molecular beam reactions of the methylidyne radical (CH; X²Π) with propylene (CH₃CHCH₂; X¹A') along with (partially) substituted reactants were conducted at collision energies of 19.3 kJ mol^-1. Combining our experimental data with ab initio electronic structure and statistical calculations, the methylidyne radical is revealed to add barrierlessly to the carbon-carbon double bond of propylene reactant resulting in a cyclic doublet C₄H₇ intermediate with a lifetime longer than its rotation period. These adducts undergo a nonstatistical unimolecular decomposition via atomic hydrogen loss through tight exit transition states forming the cyclic products 1-methylcyclopropene and 3-methylcyclopropene with overall reaction exoergicities of 168 ± 25 kJ mol^-1. These C₄H₆ isomers are predicted to exist even in low-temperature environments such as cold molecular clouds like TMC-1, since the reaction is barrierless and exoergic, all transition states are below the energy of the separated reactants, and both the methylidyne radical (CH; X²Π) and propylene reactant were detected in cold molecular clouds such as TMC-1.

A Barrierless Pathway Accessing the C9H9 and C9H8 Potential Energy Surfaces via the Elementary Reaction of Benzene with 1-Propynyl

The crossed molecular beams reactions of the 1-propynyl radical (CH₃CC; X²A₁) with benzene (C₆H₆; X¹A_1g) and D6-benzene (C₆D₆; X¹A_1g) were conducted to explore the formation of C₉H₈ isomers under single-collision conditions. The underlying reaction mechanisms were unravelled through the combination of the experimental data with electronic structure and statistical RRKM calculations. These data suggest the formation of 1-phenyl-1-propyne (C₆H₅CCCH₃) via the barrierless addition of 1-propynyl to benzene forming a low-lying doublet C₉H₉ intermediate that dissociates by hydrogen atom emission via a tight transition state. In accordance with our experiments, RRKM calculations predict that the thermodynamically most stable isomer – the polycyclic aromatic hydrocarbon (PAH) indene – is not formed via this reaction. With all barriers lying below the energy of the reactants, this reaction is viable in the cold interstellar medium where several methyl-substituted molecules have been detected. Its underlying mechanism therefore advances our understanding of how methyl-substituted hydrocarbons can be formed under extreme conditions such as those found in the molecular cloud TMC-1. Implications for the chemistry of the 1-propynyl radical in astrophysical environments are also discussed.

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.

Influence of the E-Cigarette Emission Profile by the Ratio of Glycerol to Propylene Glycol in E-Liquid Composition

The use of electronic cigarettes (E-cig) is popular because of the perception that they are less addictive and safer compared to the traditional cigarettes. Nevertheless, there are still harmful effects associated with chemicals emitted from E-cig. Identifying the sources of the emitted compounds can be challenging because of the differences in the design of E-cig devices and the variability in the composition of E-cig liquids used in these devices. In this study, the emission profiles from impurity-free E-liquids containing only propylene glycol and glycerol in various percentage ratios along with two commercially available E-liquids were evaluated using gas chromatography-mass spectrometry (GC-MS). This study approach allows the elucidation of the transformation pathways of the major emitted compounds without the confounding effects of existing impurities or flavor ingredients added to E-liquids. Analysis of the vapor phases of E-cig emissions detected toxicants such as acetaldehyde, acrolein, benzaldehyde, as well as benzene, toluene, ethylbenzene, and xylene (BTEX) compounds. The amount of glycerol in the E-liquids has a major effect on the concentration of these hazardous compounds emitted because the concentration of these chemicals increased with increasing glycerol percentage in the E-liquid. Acetaldehyde and acrolein increased by 175-fold and 28-fold, respectively, when the glycerol composition was increased from 0 to 80%. Benzaldehyde, naphthalene, diphenyl ether, and glycerol along with menthol and nicotine that were present in the commercial E-liquids were also detected in the aerosol condensates. The cascade impactor data on the distribution of the nicotine and menthol in different size fractions from >2.5 to <2.5 μm allow the estimates of the extent of toxicant deposition in different parts of the pulmonary system including the oropharynx region, the trachea as well as inside the alveoli and bronchioles. In summary, users of E-cig are exposed to harmful chemicals even if the E-liquids contain only propylene glycol and glycerol without flavorings, nicotine, or impurities. Furthermore, this study shows that E-liquids containing higher percentages of glycerol will produce higher levels of toxicants compared to E-liquids with similar percentages of propylene glycol. This finding has important implications to E-cigarette vendors and manufacturers, consumers, and regulatory agencies.

The 130-370 GHz rotational spectrum of phenyl isocyanide (C6H5NC)

The analysis of phenyl isocyanide (C₆H₅NC, μ_a = 4.0 D) in its ground vibrational state and two lowest-energy excited vibrational states, ν₂₂ (141 cm^-1) and ν₃₃ (155 cm^-1), in the 130-370 GHz frequency region has been completed. Over 4500 new rotational transitions have been measured in the ground vibrational state for the most abundant isotopologue, resulting in the determination of the spectroscopic constants for a partial octic Hamiltonian with low error. The Coriolis-coupled ν₂₂-ν₃₃ dyad reported herein, containing over 3500 new transitions for each vibrational state, has been analyzed for the first time. The coupled-state least-squares fit utilizes seven coupling terms (G_a, G_a ^J, G_a ^K, G_a ^JJ, G_a ^JK, F_bc, and F_bc ^K) to address perturbation between the two vibrational states, including resonances and several nominal interstate transitions. This work results in precise determination of the energy separation between the two states, ΔE_22,33 = 9.682 248(3) cm^-1, and the Coriolis coupling coefficient, |ζ_22,33 ^a| = 0.858(9). The precise rotational and distortion constants determined in this work provide the foundation for an astronomical search for phenyl isocyanide across the radio band.

Combined Experimental and Computational Study on the Reaction Dynamics of the 1-Propynyl (CH3CC)-1,3-Butadiene (CH2CHCHCH2) System and the Formation of Toluene under Single Collision Conditions

The crossed beams reactions of the 1-propynyl radical (CH₃CC; X²A₁) with 1,3-butadiene (CH₂CHCHCH₂; X¹A_g), 1,3-butadiene- d₆ (CD₂CDCDCD₂; X¹A_g), 1,3-butadiene- d₄ (CD₂CHCHCD₂; X¹A_g), and 1,3-butadiene- d₂ (CH₂CDCDCH₂; X¹A_g) were performed under single collision conditions at collision energies of about 40 kJ mol^-1. The underlying reaction mechanisms were unraveled through the combination of the experimental data with electronic structure calculations at the CCSD(T)-F12/cc-pVTZ-f12//B3LYP/6-311G(d,p) + ZPE(B3LYP/6-311G(d,p) level of theory along with statistical Rice-Ramsperger-Kassel-Marcus (RRKM) calculations. Together, these data suggest the formation of the thermodynamically most stable C₇H₈ isomer-toluene (C₆H₅CH₃)-via the barrierless addition of 1-propynyl to the 1,3-butadiene terminal carbon atom, forming a low-lying C₇H₉ intermediate that undergoes multiple isomerization steps resulting in cyclization and ultimately aromatization following hydrogen atom elimination. RRKM calculations predict that the thermodynamically less stable isomers 1,3-heptadien-5-yne, 5-methylene-1,3-cyclohexadiene, and 3-methylene-1-hexen-4-yne are also synthesized. Since the 1-propynyl radical may be present in cold molecular clouds such as TMC-1, this pathway could potentially serve as a carrier of the methyl group incorporating itself into methyl-substituted (poly)acetylenes or aromatic systems such as toluene via overall exoergic reaction mechanisms that are uninhibited by an entrance barrier. Such pathways are a necessary alternative to existing high energy reactions leading to toluene that are formally closed in the cold regions of space and are an important step toward understanding the synthesis of polycyclic aromatic hydrocarbons (PAHs) in space's harsh extremes.

On the formation of phosphorous polycyclic aromatics hydrocarbons (PAPHs) in astrophysical environments

The formation of phosphorous-containing polycyclic aromatic hydrocarbons (PAPHs) in astrophysical contexts is proposed and analyzed by means of computational methods [B3LYP-D3BJ/ma-def2-TZVPP, MP2-F12, CCSD-F12b and CCSD(T)-F12b levels of theory]. A "bottom-up" approach based on a radical-neutral reaction scheme between acetylene (C₂H₂) and the CP radical was used investigating: (a) the synthesis of the first PAPH (C₅H₅P) "phosphinine"; (b) PAPH growth by addition of C₂H₂ to the C₅H₄P radical; (c) PAPH synthesis by addition reactions of one CP radical and nC₂H₂ to a neutral PAH. Results show: (I) the formation of the phosphinine radical has a strong thermodynamic tendency (-133.3 kcal mol^-1) and kinetic barriers ≤5.4 kcal mol^-1; (II) PAPH growth by nC₂H₂ addition on the radical phosphinine easily and exothermically produces radicals (1a- or 1-phospha-naphtalenes with kinetic barriers ≤7.1 kcal mol^-1 and reaction free energies ≤-102.5 kcal mol^-1); (III) the addition of a single CP + nC₂H₂ to a neutral benzene generates a complex chemistry where the main product is 2-phospha-naphtalene; (IV) because of the CP radical character, its barrierless addition to a PAH produces a resonant stabilized PAPH, becoming excellent candidates for addition reactions with neutral or radical hydrocarbons and PAHs; (V) the same energy trend between all four levels of theory continues a well-calibrated computational protocol to analyze complex organic reactions with astrochemical interest using electronic structure theory.