Product Branching Ratios of the C(3P) + C2H3(2A‘) and CH(2Π) + C2H2(1Σg+) Reactions and Photodissociation of H2CC⋮CH(2B1) at 193 and 242 nm: an ab Initio/RRKM Study

Interstellar hide and go seek: C3H4O. There and back (again)

The molecular species C₃H₄O represents a striking example of an astrochemical conundrum. With more than 60 structural isomers theoretically possible, to date only acrolein (CH₂CHCHO) has been identified in the Sgr B2(N) region of the interstellar medium (ISM). The topography of the singlet potential energy surface is complicated, with three low-lying minima predicted to be almost isoenergetic: cis and trans-acrolein, and methylketene (CH₃CHCO). Our CCSD(T)/cc-pVTZ calculations confirm that methylketene is energetically lower than cis-acrolein, lying only 1.9 kJ mol^-1 above the trans-isomer, which is the global minimum. In this respect, methylketene is a promising candidate for interstellar observation. Unfortunately, however, despite several searches its astronomical detection has been unsuccessful. To this end, the key question is whether in fact methylketene exists as a discrete chemical entity in the ISM at all? In this paper, we present a detailed examination of the C₃H₄O potential energy surface, with specific focus on formation pathways. CCSD(T)/cc-pVTZ calculations enable a more elaborate interpretation of reaction mechanisms than was published hitherto. Our results show that gauche-propargyl alcohol and syn and anti-allenol emerge as interesting new targets for observational astronomers in TMC-1: given the recent discovery of the propargyl radical in this region, barrierless product channels involving OH˙ lend support to their candidacy as possible interstellar species. Finally, this work provides accurate spectral data of these three potential molecules, to be used for searches in interstellar space.

Theoretical Calculations of the 242 nm Absorption of Propargyl Radical

The propargyl radical, the most stable isomer of neutral C₃H₃, is important in combustion reactions, and a number of spectroscopic and reaction dynamics studies have been performed over the years. However, theoretical calculations have never been able to find a state that can generate strong absorption around 242 nm as seen in experiments. In this study, we calculated the low-lying electronic energy levels of the propargyl radical using the highly accurate multireference configuration interaction singles and doubles method with triples and quadruples treated perturbatively [denoted as MRCISD(TQ)]. Calculations indicate that this absorption can be attributed to a Franck-Condon-allowed electronic transition from the ground ²B₁ state to the Rydberg-like excited state 1²A₁. Further insight into the behavior of the multireference perturbative theory methods, GVVPT2 and GVVPT3, on a very challenging system are also obtained.

Photoinduced C–H bond fission in prototypical organic molecules and radicals

We survey and assess current knowledge regarding the primary photochemistry of hydrocarbon molecules and radicals.

Isomer-specific detection in the UV photodissociation of the propargyl radical by chirped-pulse mm-wave spectroscopy in a pulsed quasi-uniform flow

Isomer-specific detection and product branching fractions in the UV photodissociation of the propargyl radical is achieved through the use of chirped-pulse Fourier-transform mm-wave spectroscopy in a pulsed quasi-uniform flow (CPUF). Propargyl radicals are produced in the 193 nm photodissociation of 1,2-butadiene. Absorption of a second photon leads to H atom elimination giving three possible C₃H₂ isomers: singlets cyclopropenylidene (c-C₃H₂) and propadienylidene (l-C₃H₂), and triplet propargylene (³HCCCH). The singlet products and their appearance kinetics in the flow are directly determined by rotational spectroscopy, but due to the negligible dipole moment of propargylene, it is not directly monitored. However, we exploit the time-dependent kinetics of H-atom catalyzed isomerization to infer the branching to propargylene as well. We obtain the overall branching among H loss channels to be 2.9% (+1.1/-0.5) l-C₃H₂ + H, 16.8% (+3.2/-1.3) c-C₃H₂ + H, and 80.2 (+1.8/-4.2) ³HCCCH + H. Our findings are qualitatively consistent with earlier RRKM calculations in that the major channel in the photodissociation of the propargyl radical at 193 nm is to ³HCCCH + H; however, a greater contribution to the energetically most favorable isomer, c-C₃H₂ + H is observed in this work. We do not detect the predicted HCCC + H₂ channel, but place an upper bound on its yield of 1%.

Valence shell threshold photoelectron spectroscopy of C3Hx (x = 0–3)

We present the photoelectron spectra of C₃H_x (x = 0–3) formed in a microwave discharge flow-tube reactor by consecutive H abstractions from C₃H₄ (C₃H_x + F → C₃H_x−1 + HF (x = 1–4)), but also from F + CH₄ schemes by secondary reactions.

The interstellar chemistry of C3H and C3H2 isomers

We report the detection of linear and cyclic isomers of C3H and C3H2 towards various starless cores and review the corresponding chemical pathways involving neutral (C3Hx with x=1,2) and ionic (C3Hx+ with x = 1,2,3) isomers. We highlight the role of the branching ratio of electronic Dissociative Recombination (DR) reactions of C3H2+ and C3H3+ isomers showing that the statistical treatment of the relaxation of C3H* and C3H2* produced in these DR reactions may explain the relative c,l-C3H and c,l-C3H2 abundances. We have also introduced in the model the third isomer of C3H2 (HCCCH). The observed cyclic-to-linear C3H2 ratio vary from 110 ± 30 for molecular clouds with a total density around 1×104 molecules.cm-3 to 30 ± 10 for molecular clouds with a total density around 4×105 molecules.cm-3, a trend well reproduced with our updated model. The higher ratio for low molecular cloud densities is mainly determined by the importance of the H + l-C3H2 → H + c-C3H2 and H + t-C3H2 → H + c-C3H2 isomerization reactions.

Theoretical Study and Design of Phosphorescent Cyclometalated (C∧C*)PtII(acac) Complexes: The Substituent Effect Controls the Radiative and Nonradiative Decay Processes

Density functional theory (DFT) and time-dependent DFT calculations were performed to evaluate the influence of substituent effect of (1) R = 4-Me, (2) R = 4-OMe, and (3) R = 2,3-OC₆H₄ on the phenyl ring of (C^∧C*)Pt^II(acac) (C^∧C* = phenylimidazole, acac = acetylacetone), respectively, on absorption and phosphorescent spectra properties, as well as the radiative and nonradiative processes. We found that emissions of complexes 2 and 3 originate from the Kasha-like T₁ state, whereas that of complex 1 originates from non-Kasha T₂ state. Compared with the emission of complex 1, the emission peaks of 2 and 3 are red-shifted, which is attributed to p-π and π-π conjugation effects resulting from the electron-donating groups -OCH₃ and -OC₆H₄ with ligand C^∧C*, respectively. The radiative rate constants (κ_r) of 2 and 3 are larger than that of 1, namely, κ_r(1) < κ_r(2) < κ_r(3), indicating that κ_r can be efficiently increased by enlarging π-conjugation at the main ligand of (C^∧C*)Pt^II(acac), which can cause the increase of spin-orbit coupling (SOC) matrix elements. At the same time, the activation energy barriers for the rate-limiting step can be largely raised accompanied by enlarging the ability of electron-donation of the substituent group at the main ligand of (C^∧C*)Pt^II(acac), which can cause the decrease of the nonradiative rate constant (κ_nr), namely, κ_nr(1) > κ_nr(2) > κ_nr(3). According to Φ_P = κ_r/(κ_r + κ_nr), the quantum yields should have the sequence Φ_P(1) < Φ_P(2) < Φ_P(3), which is in accordance with the experiment. In addition, to guide experimental synthesis of highly efficient (C^∧C*)Pt^II(acac), a new complex 4 through extending the π-conjugation in the C^∧C* ligand of (C^∧C*)Pt^II(acac) was theoretically designed, which has a larger quantum yield than 1-3.

Insights into gas-phase reaction mechanisms of small carbon radicals using isomer-resolved product detection

Gas-phase radical reactions of CN and CH with small hydrocarbons are overviewed with emphasis on isomer-resolved product detection.

Formation of C3H2, C5H2, C7H2, and C9H2 from reactions of CH, C3H, C5H, and C7H radicals with C2H2

A three-dimensional velocity distribution contour of C_2n+1H₂ produced from the reaction of C_2n−1H (n = 1–4) with C₂H₂ in crossed molecular beams.

Kinetics of Propargyl Radical Dissociation

Due to the prominent role of the propargyl radical for hydrocarbon growth within combustion environments, it is important to understand the kinetics of its formation and loss. The ab initio transition state theory-based master equation method is used to obtain theoretical kinetic predictions for the temperature and pressure dependence of the thermal decomposition of propargyl, which may be its primary loss channel under some conditions. The potential energy surface for the decomposition of propargyl is first mapped at a high level of theory with a combination of coupled cluster and multireference perturbation calculations. Variational transition state theory is then used to predict the microcanonical rate coefficients, which are subsequently implemented within the multiple-well multiple-channel master equation. A variety of energy transfer parameters are considered, and the sensitivity of the thermal rate predictions to these parameters is explored. The predictions for the thermal decomposition rate coefficient are found to be in good agreement with the limited experimental data. Modified Arrhenius representations of the rate constants are reported for utility in combustion modeling.

Gas-Phase Synthesis of Boronylallene (H2CCCH(BO)) under Single Collision Conditions: A Crossed Molecular Beams and Computational Study

The gas phase reaction between the boron monoxide radical (¹¹BO; X²Σ⁺) and allene (H₂CCCH₂; X¹A₁) was investigated experimentally under single collision conditions using the crossed molecular beam technique and theoretically exploiting ab initio electronic structure and statistical (RRKM) calculations. The reaction was found to follow indirect (complex forming) scattering dynamics and proceeded via the formation of a van der Waals complex (¹¹BOC₃H₄). This complex isomerized via addition of the boron monoxide radical (¹¹BO; X²Σ⁺) with the radical center located at the boron atom to the terminal carbon atom of the allene molecule forming a H₂CCCH₂¹¹BO intermediate on the doublet surface. The chemically activated H₂CCCH₂¹¹BO intermediate underwent unimolecular decomposition via atomic hydrogen elimination from the terminal carbon atom holding the boronyl group through a tight exit transition state to synthesize the boronylallene product (H₂CCCH¹¹BO) in a slightly exoergic reaction (55 ± 11 kJ mol^-1). Statistical (RRKM) calculations suggest that minor reaction channels lead to the products 3-propynyloxoborane (CH₂(¹¹BO)CCH) and 1-propynyloxoborane (CH₃CC¹¹BO) with fractions of 1.5% and 0.2%, respectively. The title reaction was also compared with the cyano (CN; X²Σ⁺)-allene and boronyl-methylacetylene reactions to probe similarities, but also differences of these isoelectronic systems. Our investigation presents a novel gas phase synthesis and characterization of a hitherto elusive organyloxoborane (RBO) monomer-boronylallene-which is inherently tricky to isolate in the condensed phase except in matrix studies; our work further demonstrates that the crossed molecular beams approach presents a useful tool in investigating the chemistry and synthesis of highly reactive organyloxoboranes.

A crossed molecular beam and ab-initio investigation of the reaction of boron monoxide (BO; X2Σ+) with methylacetylene (CH3CCH; X1A1): competing atomic hydrogen and methyl loss pathways

The gas-phase reaction of boron monoxide ((11)BO; X(2)Σ(+)) with methylacetylene (CH3CCH; X(1)A1) was investigated experimentally using crossed molecular beam technique at a collision energy of 22.7 kJ mol(-1) and theoretically using state of the art electronic structure calculation, for the first time. The scattering dynamics were found to be indirect (complex forming reaction) and the reaction proceeded through the barrier-less formation of a van-der-Waals complex ((11)BOC3H4) followed by isomerization via the addition of (11)BO(X(2)Σ(+)) to the C1 and/or C2 carbon atom of methylacetylene through submerged barriers. The resulting (11)BOC3H4 doublet radical intermediates underwent unimolecular decomposition involving three competing reaction mechanisms via two distinct atomic hydrogen losses and a methyl group elimination. Utilizing partially deuterated methylacetylene reactants (CD3CCH; CH3CCD), we revealed that the initial addition of (11)BO(X(2)Σ(+)) to the C1 carbon atom of methylacetylene was followed by hydrogen loss from the acetylenic carbon atom (C1) and from the methyl group (C3) leading to 1-propynyl boron monoxide (CH3CC(11)BO) and propadienyl boron monoxide (CH2CCH(11)BO), respectively. Addition of (11)BO(X(2)Σ(+)) to the C1 of methylacetylene followed by the migration of the boronyl group to the C2 carbon atom and/or an initial addition of (11)BO(X(2)Σ(+)) to the sterically less accessible C2 carbon atom of methylacetylene was followed by loss of a methyl group leading to the ethynyl boron monoxide product (HCC(11)BO) in an overall exoergic reaction (78 ± 23 kJ mol(-1)). The branching ratios of these channels forming CH2CCH(11)BO, CH3CC(11)BO, and HCC(11)BO were derived to be 4 ± 3%, 40 ± 5%, and 56 ± 15%, respectively; these data are in excellent agreement with the calculated branching ratios using statistical RRKM theory yielding 1%, 38%, and 61%, respectively.

Crossed beam study of the atom-radical reaction of ground state carbon atoms (C(3P)) with the vinyl radical (C2H3(X2A'))

The atom-radical reaction of ground state carbon atoms (C((3)P)) with the vinyl radical (C(2)H(3)(X(2)A')) was conducted under single collision conditions at a collision energy of 32.3 ± 2.9 kJ mol(-1). The reaction dynamics were found to involve a complex forming reaction mechanism, which is initiated by the barrier-less addition of atomic carbon to the carbon-carbon-double bond of the vinyl radical forming a cyclic C(3)H(3) radical intermediate. The latter has a lifetime of at least 1.5 times its rotational period and decomposes via a tight exit transition state located about 45 kJ mol(-1) above the separated products through atomic hydrogen loss to the cyclopropenylidene isomer (c-C(3)H(2)) as detected toward cold molecular clouds and in star forming regions.

Crossed-beam and theoretical studies of the S(1D) + C2H2 reaction

The reaction dynamics of excited sulfur atoms, S((1)D), with acetylene has been investigated by the crossed-beam scattering technique with mass spectrometric detection and time-of-flight (TOF) analysis at the collision energy of 35.6 kJ mol(-1). These studies have been made possible by the development of intense continuous supersonic beams of S((3)P,(1)D) atoms. From product angular and TOF distributions, center-of-mass product angular and translational energy distributions are derived. The S((1)D) + C(2)H(2) reaction is found to lead to formation of HCCS (thioketenyl) + H, while the only other energetically allowed channels, those leading to CCS((3)Sigma(-), (1)Delta) + H(2), are not observed to occur to an appreciable extent. The dynamics of the H-elimination channel is discussed and elucidated. The interpretation of the scattering results is assisted by synergic high-level ab initio electronic structure calculations of stationary points and product energetics for the C(2)H(2)S ground-state singlet potential energy surface. In addition, by exploiting the novel capability of performing product detection by means of a tunable electron-impact ionizer, we have obtained the first experimental information on the ionization energy of thioketenyl radical, HCCS, as synthesized in the reactive scattering experiment. This has been complemented by ab initio calculations of the adiabatic and vertical ionization energies for the ground-state radical. The theoretically derived value of 9.1 eV confirms very recent, accurate calculations and is corroborated by the experimentally determined ionization threshold of 8.9 +/- 0.3 eV for the internally warm HCCS produced from the title reaction.

Quantum Chemical Evaluation of the Astrochemical Significance of Reactions between S Atom and Acetylene or Ethylene

Addition-elimination reactions of S atom in its 3P ground state with acetylene (C2H2) and ethylene (C2H4) were characterized with both molecular orbital and density functional theory calculations employing correlation consistent basis sets in order to assess the likelihood that either reaction might play a general role in astrochemistry or a specific role in the formation of S2 (X3Sigmag-) via a mechanism proposed by Saxena, P. P.; Misra, A. Mon. Not. R. Astron. Soc. 1995, 272, 89. The acetylene and ethylene reactions proceed through C2H2S (3A' ') and C2H4S (3A' ') intermediates, respectively, to yield HCCS (2Pi) and C2H3S (2A'). Substantial barriers were found in the exit channels for every combination of method and basis set considered in this work, which effectively precludes hydrogen elimination pathways for both S + C2H2 and S + C2H4 in the ultracold interstellar medium where only very modest barriers can be surmounted and processes without barriers tend to predominate. However, if one or both intermediates are formed and stabilized efficiently under cometary or dense interstellar cloud conditions, they could serve as temporary reservoirs for the S atom and participate in reactions such as S + C2H2S --> S2 + C2H2 or S + C2H4S --> S2 + C2H4. For formation and stabilization to be efficient, the reaction must possess a barrier height small enough to be surmountable at low temperatures yet large enough to prevent redissociation to reactants. Barrier heights computed with B3LYP and large basis sets are very low, but more rigorous QCISD(T) and RCCSD(T) results indicate that the barrier heights are closer to 3-4 kcal/mol. The calculations therefore indicate that S + C2H2 or S + C2H4 could contribute to the formation of S2 in comets and may serve as a means to gauge coma temperature. The energetics of the ethylene reaction are more favorable.

Unimolecular dissociation of the propargyl radical intermediate of the CH+C2H2 and C+C2H3 reactions

This paper examines the unimolecular dissociation of propargyl (HCCCH2) radicals over a range of internal energies to probe the CH+HCCH and C+C2H3 bimolecular reactions from the radical intermediate to products. The propargyl radical was produced by 157 nm photolysis of propargyl chloride in crossed laser-molecular beam scattering experiments. The H-loss and H2 elimination channels of the nascent propargyl radicals were observed. Detection of stable propargyl radicals gave an experimental determination of 71.5 (+5-10) kcal/mol as the lowest barrier to dissociation of the radical. This barrier is significantly lower than predictions for the lowest barrier to the radical's dissociation and also lower than calculated overall reaction enthalpies. Products from both H2+HCCC and H+C3H2 channels were detected at energies lower than what has been theoretically predicted. An HCl elimination channel and a minor C-H fission channel were also observed in the photolysis of propargyl chloride.

The temperature-dependence of rapid low temperature reactions: experiment, understanding and prediction

Despite the success of the CRESU (Cinétique de Réaction en Ecoulement Supersonique Uniforme) method in measuring rate coefficients for neutral-neutral reactions of radicals down close to the very low temperatures prevalent in dense interstellar clouds (ISCs), there are still many reactions of potential importance in the chemistry of these objects for which there have been no measurements of low temperature rate coefficients. One important class of reactions is that between atomic and molecular free radicals and unsaturated hydrocarbons; that is, alkynes and alkenes. Based on semi-empirical arguments and correlations of 'room temperature' rate coefficients, k(298 K), for reactions of this type with the difference between the ionisation energy of the alkyne/alkene and the electron affinity of the radical, we suggest which reactions between the radicals, C(3P), O(3P), N(4S), CH, C2H and CN, and carbon chain molecules (Cn) and cyanopolyynes (HC2nCN and NCC2nCN) are likely to be fast at the temperature of dense ISCs. These reactions and rate coefficients have been incorporated into a purely gas-phase model (osu2005) of ISC chemistry. The results of these calculations are presented and discussed.

Potential Energy Surface and Product Branching Ratios for the Reaction of Dicarbon, C2(XΣg+), with Methylacetylene, CH3CCH(X1A1): An Ab Initio/RRKM Study

Ab initio calculations of the potential energy surface for the C(2)(X(1)Sigma(g)(+)) + CH(3)CCH(X(1)A(1)) reaction have been carried at the G2M level of theory. The calculations show that the dicarbon molecule in the ground singlet electronic state can add to methylacetylene without a barrier producing a three-member or a four-member ring intermediate, which can rapidly rearrange to the most stable H(3)CCCCCH isomer on the C(5)H(4) singlet surface. This isomer can then lose a hydrogen atom (H) or molecular hydrogen (H(2)) from the CH(3) group with the formation of H(2)CCCCCH and HCCCCCH, respectively. Alternatively, H atom migrations and three-member-ring closure/opening rearrangements followed by H and H(2) losses can lead to other isomers of the C(5)H(3) and C(5)H(2) species. According to the calculated energetics, the C(2)(X(1)Sigma(g)(+)) + CH(3)CCH reaction is likely to be a major source of the C(5)H(3) radicals (in particular, the most stable H(2)CCCCCH and HCCCHCCH isomers, which are relevant to the formation of benzene through the reactions with CH(3)). Among heavy-fragment product channels, only C(3)H(3) + C(2)H and c-C(3)H(2) + C(2)H(2) might compete with C(5)H(3) + H and C(5)H(2) + H(2). RRKM calculations of reaction rate constants and product branching ratios depending on the reactive collision energy showed that the major reaction products are expected to be H(2)CCCCCH + H (64-66%) and HCCCHCCH + H (34-30%), with minor contributions from HCCCCCH + H(2) (1-2%), HCCCHCC + H(2) (up to 1%), C(3)H(3) + C(2)H (up to 1%), and c-C(3)H(2) + C(2)H(2) (up to 0.1%) if the energy randomization is complete. The calculations also indicate that the C(2)(X(1)Sigma(g)(+)) + CH(3)CCH(X(1)A(1)) reaction can proceed by direct H-abstraction of a methyl hydrogen to form C(3)H(3) + C(2)H almost without a barrier.

Kinetic measurements on methylidyne radical reactions with several hydrocarbons at low temperatures

The temperature dependences of the methylidyne radical reactions with methane, allene, methylacetylene and propene were studied. This work was carried out in a supersonic flow reactor coupled with pulsed laser photolysis (PLP) and laser-induced fluorescence (LIF) techniques. Three Laval nozzles were designed to provide uniform supersonic expansions of nitrogen at Mach 2 and of argon at Mach 2 and 3 to reach low temperatures, e.g. 170, 128 and 77 K, respectively. CH radicals were produced by PLP of CHBr3 at 266 nm and probed by LIF. The exponential decays of the CH fluorescence were acquired, hydrocarbons being introduced in excess. The rate constants for the CH+CH4 reaction are in good agreement with the temperature dependence proposed by Canosa et al. (A. Canosa, I. R. Sims, D. Travers, I. W. M. Smith and B. R. Rowe, Astron. Astrophys., 1997, 323, 644-651, ) i.e. 3.96x10(-8)(T/K)(-1.04) exp(-36.1 K/T) in the range 23-298 K. The rate constants of the CH+C3H4(allene), CH+C3H4(methylacetylene) and CH+C3H6(propene) reactions exhibit a small temperature dependence between 77 and 170 K, with a maximum rate around 100 K close to (4.3-4.6)x10(-10) cm3 molecule-1 s-1.