Chemical dynamics of the formation of the 1,3-butadiynyl radical (C4H(X2Sigma+)) and its isotopomers

The reaction of dicarbon molecules in their electronic ground, C2(X1Sigma(g)+), and first excited state, C2(a3Pi(u)), with acetylene, C2H2(X1Sigma(g)+), to synthesize the 1,3-butadiynyl radical, C4H(X2Sigma+), plus a hydrogen atom was investigated at six different collision energies between 10.6 and 47.5 kJ mol(-1) under single collision conditions. These studies were contemplated by crossed molecular beam experiments of dicarbon with three acetylene isotopomers C2D2(X1Sigma(g)+), C2HD (X1Sigma+), and 13C2H2(X1Sigma(g)+) to elucidate the role of intersystem crossing (ISC) and of the symmetry of the reaction intermediate(s) on the center-of-mass functions. On the singlet surface, dicarbon was found to react with acetylene through an indirect reaction mechanism involving a diacetylene intermediate. The latter fragmented via a loose exit transition state via an emission of a hydrogen atom to form the 1,3-butadiynyl radical C4H(X2Sigma+). The D(infinity)(h) symmetry of the decomposing diacetylene intermediate results in collision-energy invariant, isotropic (flat) center-of-mass angular distributions of this microchannel. Isotopic substitution experiments suggested that at least at a collision energy of 29 kJ mol(-1), the diacetylene isotopomers are long-lived with respect to their rotational periods. On the triplet surface, the reaction involved three feasible addition complexes located in shallower potential energy wells as compared to singlet diacetylene. The involvement of the triplet surface accounted for the asymmetry of the center-of-mass angular distributions. The detection of the 1,3-butadiynyl radical, C4H(X2Sigma+), in the crossed beam reaction of dicarbon molecules with acetylene presents compelling evidence that the 1,3-butadiynyl radical can be formed via bimolecular reactions involving carbon clusters in extreme environments such as circumstellar envelopes of dying carbon stars and combustion flames.     

Formation of the 2,4-pentadiynyl-1 radical (H2CCCCCH, X2B1) in the crossed beams reaction of dicarbon molecules with methylacetylene

The chemical dynamics to synthesize the 2,4-pentadiynyl-1 radical, HCCCCCH(2)(X(2)B(1)), via the neutral-neutral reaction of dicarbon with methylacetylene, was examined in a crossed molecular beams experiment at a collision energy of 37.6 kJ mol(-1). The laboratory angular distribution and time-of-flight spectra of the 2,4-pentadiynyl-1 radical and its fragmentation patterns were recorded at m/z = 63-60 and m/z = 51-48. Our findings suggest that the reaction dynamics are indirect and dictated by an initial attack of the dicarbon molecule to the pi electron density of the methylacetylene molecule to form cyclic collision complexes. The latter ultimately rearranged via ring opening to methyldiacetylene, CH(3)-C triple bond C-C triple bond C-H. This structure decomposed via atomic hydrogen emission to the 2,4-pentadiynyl-1 radical; here, the hydrogen atom was found to be emitted almost parallel to the total angular momentum as suggested by the experimentally observed sideways scattering. The overall reaction was strongly exoergic by 182 +/- 10 kJ mol(-1). The identification of the resonance-stabilized free 2,4-pentadiynyl-1 radical represents a solid background for the title reaction to be included into more refined reaction networks modeling the chemistry of circumstellar envelopes and also of sooting combustion flames.     

A crossed molecular beams and ab initio study on the formation of C6H3 radicals. an interface between resonantly stabilized and aromatic radicals

The crossed molecular beams reaction of dicarbon molecules, C(2)(X(1)Σ(g)(+)/a(3)Π(u)) with vinylacetylene was studied under single collision conditions at a collision energy of 31.0 kJ mol(-1) and combined with electronic structure calculations on the singlet and triplet C(6)H(4) potential energy surfaces. The investigations indicate that both reactions on the triplet and singlet surfaces are dictated by a barrierless addition of the dicarbon unit to the vinylacetylene molecule and hence indirect scattering dynamics via long-lived C(6)H(4) complexes. On the singlet surface, ethynylbutatriene and vinyldiacetylene were found to decompose via atomic hydrogen loss involving loose exit transition states to form exclusively the resonantly stabilized 1-hexene-3,4-diynyl-2 radical (C(6)H(3); H(2)CCCCCCH; C(2v)). On the triplet surface, ethynylbutatriene emitted a hydrogen atom through a tight exit transition state located about 20 kJ mol(-1) above the separated stabilized 1-hexene-3,4-diynyl-2 radical plus atomic hydrogen product; to a minor amount (<5%) theory predicts that the aromatic 1,2,3-tridehydrobenzene molecule is formed. Compared to previous crossed beams and theoretical investigations on the formation of aromatic C(6)H(x) (x = 6, 5, 4) molecules benzene, phenyl, and o-benzyne, the decreasing energy difference from benzene via phenyl and o-benzyne between the aromatic and acyclic reaction products, i.e., 253, 218, and 58 kJ mol(-1), is narrowed down to only ～7 kJ mol(-1) for the C(6)H(3) system (aromatic 1,2,3-tridehydrobenzene versus the resonantly stabilized free radical 1-hexene-3,4-diynyl-2). Therefore, the C(6)H(3) system can be seen as a "transition" stage among the C(6)H(x) (x = 6-1) systems, in which the energy gap between the aromatic isomer (x = 6, 5, 4) is reduced compared to the acyclic isomer as the carbon-to-hydrogen ratio increases and the acyclic isomer becomes more stable (x = 1, 2).     

Investigating the chemical dynamics of the reaction of ground-state carbon atoms with acetylene and its isotopomers

We investigated the multichannel reaction of ground-state carbon atoms with acetylene, C2H2 (X1Sigmag+), to form the linear and cyclic C3H isomers (atomic hydrogen elimination pathway) as well as tricarbon plus molecular hydrogen. The experiments were conducted under single-collision conditions at three different collision energies between 8.0 kJ mol-1 and 31.0 kJ mol-1. Our studies were complemented by crossed molecular beam experiments of carbon with three isotopomers C2D2(X1Sigmag+), C2HD (X1Sigma+), and 13C2H2 (X1Sigmag+) to clarify a potential intersystem crossing (ISC), the effect of the symmetry of the reaction intermediates on the center-of-mass angular distributions, the collision energy-dependent branching ratios of the atomic versus molecular hydrogen elimination pathways, and deuterium-enrichment processes. The results are discussed in light of recent electronic structure and dynamics calculations.     

Crossed molecular beam study on the reaction of boron atoms, B(2Pj), with allene, H2CCCH2(X1A1)

The reaction of ground state boron atoms, 11B(2Pj), with allene, H2CCCH2(X1A1), was studied under single collision conditions at a collision energy of 21.5 kJ mol(-1) utilizing the crossed molecular beam technique; the experimental data were combined with electronic structure calculations on the 11BC3H4 potential energy surface. The chemical dynamics were found to be indirect and initiated by an addition of the boron atom to the pi-electron density of the allene molecule leading ultimately to a cyclic reaction intermediate. The latter underwent ring-opening to yield an acyclic intermediate H2CCBCH2. As derived from the center-of-mass functions, this structure was long-lived with respect to its rotational period and decomposed via an atomic hydrogen loss through a tight exit transition state to form the closed shell, C2v symmetric H-C is equivalent C-B=CH2 molecule. A brief comparison of the product isomers formed in the reaction of boron atoms with methylacetylene is also presented.     

A crossed beams study of the reaction of carbon atoms, C(3Pj), with vinyl cyanide, C2H3CN(X 1A')--investigating the formation of cyano propargyl radicals

The chemical dynamics of the reaction of ground state carbon atoms, C(3Pj), with vinyl cyanide, C2H3CN(X 1A'), were examined under single collision conditions at collision energies of 29.9 and 43.9 kJ mol(-1) using the crossed molecular beams approach. The experimental studies were combined with electronic structure calculations on the triplet C4H3N potential energy surface (H. F. Su, R. I. Kaiser, A. H. H. Chang, J. Chem. Phys., 2005, 122, 074320). Our investigations suggest that the reaction follows indirect scattering dynamics via addition of the carbon atom to the carbon-carbon double bond of the vinyl cyanide molecule yielding a cyano cyclopropylidene collision complex. The latter undergoes ring opening to form cis/trans triplet cyano allene which fragments predominantly to the 1-cyano propargyl radical via tight exit transition states; the 3-cyano propargyl isomer was inferred to be formed at least a factor of two less; also, no molecular hydrogen elimination channel was observed experimentally. These results are in agreement with the computational studies predicting solely the existence of a carbon versus hydrogen atom exchange pathway and the dominance of the 1-cyano propargyl radical product. The discovery of the cyano propargyl radical in the reaction of atomic carbon with vinyl cyanide under single collision conditions implies that this molecule can be an important reaction intermediate in combustion flames and also in extraterrestrial environments (cold molecular clouds, circumstellar envelopes of carbon stars) which could lead to the formation of cyano benzene (C6H5CN) upon reaction with a propargyl radical.     

Unraveling the dynamics of the C(3P,1D) + C2H2 reactions by the crossed molecular beam scattering technique

A detailed investigation of the dynamics of the reactions of ground- and excited-state carbon atoms, C(3P) and C(1D), with acetylene is reported over a wide collision energy range (3.6-49.1 kJ mol-1) using the crossed molecular beam (CMB) scattering technique with electron ionization mass spectrometric detection and time-of-flight (TOF) analysis. We have exploited the capability of (a) generating continuous intense supersonic beams of C(3P, 1D), (b) crossing the two reactant beams at different intersection angles (45, 90, and 135 degrees ) to attain a wide range of collision energies, and (c) tuning the energy of the ionizing electrons to low values (soft ionization) to suppress interferences from dissociative ionization processes. From angular and TOF distribution measurements of products at m/z=37 and 36, the primary reaction products of the C(3P) and C(1D) reactions with C2H2 have been identified to be cyclic (c)-C3H + H, linear (l)-C3H + H, and C3 + H2. From the data analysis, product angular and translational energy distributions in the center-of-mass (CM) system for both the linear and cyclic C3H isomers as well as the C3 product from C(3P) and for l/c-C3H and C3 from C(1D) have been derived as a function of collision energy from 3.6 to 49.1 kJ mol-1. The cyclic/linear C3H ratio and the C3/(C3 + c/l-C3H) branching ratios for the C(3P) reaction have been determined as a function of collision energy. The present findings have been compared with those from previous CMB studies using pulsed beams; here, a marked contrast is noted in the CM angular distributions for both C3H- and C3-forming channels from C(3P) and their trend with collision energy. Consequently, the interpretation of the reaction dynamics derived in the present work contradicts that previously proposed from the pulsed CMB studies. The results have been discussed in the light of the available theoretical information on the relevant triplet and singlet C3H2 ab initio potential energy surfaces (PESs). In particular, the branching ratios for the C(3P) + C2H2 reaction have been compared with the available theoretical predictions (approximate quantum scattering calculations and quasiclassical trajectory calculations on ab initio triplet PESs and, very recent, statistical calculations on ab initio triplet PESs as well as on ab initio triplet/singlet PESs including nonadiabatic effects, that is, intersystem crossing). While the experimental branching ratios have been corroborated by the statistical predictions, strong disagreement has been found with the results of the dynamical calculations. The astrophysical implications of the present results have been noted.     

A crossed molecular beam study on the formation of hexenediynyl radicals (H(2)CCCCCCH; C(6)H(3) (X(2)A')) via reactions of tricarbon molecules, C(3)(X(1)Sigma(g)(+)), with allene (H(2)CCCH(2); X(1)A(1)) and methylacetylene (CH(3)CCH; X(1)A(1))

Crossed molecular beams experiments have been utilized to investigate the reaction dynamics between two closed shell species, i.e. the reactions of tricarbon molecules, C(3)(X(1)Sigma(g)(+)), with allene (H(2)CCCH(2); X(1)A(1)), and with methylacetylene (CH(3)CCH; X(1)A(1)). Our investigations indicated that both these reactions featured characteristic threshold energies of 40-50 kJ mol(-1). The reaction dynamics are indirect and suggested the reactions proceeded via an initial addition of the tricarbon molecule to the unsaturated hydrocarbon molecules forming initially cyclic reaction intermediates of the generic formula C(6)H(4). The cyclic intermediates isomerize to yield eventually the acyclic isomers CH(3)CCCCCH (methylacetylene reaction) and H(2)CCCCCCH(2) (allene reaction). Both structures decompose via atomic hydrogen elimination to form the 1-hexene-3,4-diynyl-2 radical (C(6)H(3); H(2)CCCCCCH). Future flame studies utilizing the Advanced Light Source should therefore investigate the existence of 1-hexene-3,4-diynyl-2 radicals in high temperature methylacetylene and allene flames. Since the corresponding C(3)H(3), C(4)H(3), and C(5)H(3) radicals have been identified via their ionization potentials in combustion flames, the existence of the C(6)H(3) isomer 1-hexene-3,4-diynyl-2 can be predicted as well.     

Trajectory surface hopping study of the O(3P) + ethylene reaction dynamics

Spin-orbit coupling (SOC) induced intersystem crossing (ISC) has long been believed to play a crucial role in determining the product distributions in the O(3P) + C2H4 reaction. In this paper, we present the first nonadiabatic dynamics study of the title reaction at two center-of-mass collision energies: 0.56 eV, which is barely above the H-atom abstraction barrier on the triplet surface, and 3.0 eV, which is in the hyperthermal regime. The calculations were performed using a quasiclassical trajectory surface hopping (TSH) method with the potential energy surface generated on the fly at the unrestricted B3LYP/6-31G(d,p) level of theory. To simplify our calculations, nonadiabatic transitions were only considered when the singlet surface intersects the triplet surface. At the crossing points, Landau-Zener transition probabilities were computed assuming a fixed spin-orbit coupling parameter, which was taken to be 70 cm-1 in most calculations. Comparison with a recent crossed molecular beam experiment at 0.56 eV collision energy shows qualitative agreement as to the primary product branching ratios, with the CH3 + CHO and H + CH2CHO channels accounting for over 70% of total product formation. However, our direct dynamics TSH calculations overestimate ISC so that the total triplet/singlet ratio is 25:75, compared to the observed 43:57. Smaller values of SOC reduce ISC, resulting in better agreement with the experimental product relative yields; we demonstrate that these smaller SOC values are close to being consistent with estimates based on CASSCF calculations. As the collision energy increases, ISC becomes much less important and at 3.0 eV, the triplet to singlet branching ratio is 71:29. As a result, the triplet products CH2 + CH2O, H + CH2CHO and OH + C2H3 dominate over the singlet products CH3 + CHO, H2 + CH2CO, etc.     

Unraveling the formation of HCPH(X2A') molecules in extraterrestrial environments: crossed molecular beam study of the reaction of carbon atoms, C(3Pj), with phosphine, PH3(X1A1)

The reaction between ground state carbon atoms, C(3P(j)), and phosphine, PH3(X(1)A1), was investigated at two collision energies of 21.1 and 42.5 kJ mol(-1) using the crossed molecular beam technique. The chemical dynamics extracted from the time-of-flight spectra and laboratory angular distributions combined with ab initio calculations propose that the reaction proceeds on the triplet surface via an addition of atomic carbon to the phosphorus atom. This leads to a triplet CPH3 complex. A successive hydrogen shift forms an HCPH2 intermediate. The latter was found to decompose through atomic hydrogen emission leading to the cis/trans-HCPH(X(2)A') reaction products. The identification of cis/trans-HCPH(X(2)A') molecules under single collision conditions presents a potential pathway to form the very first carbon-phosphorus bond in extraterrestrial environments like molecular clouds and circumstellar envelopes, and even in the postplume chemistry of the collision of comet Shoemaker-Levy 9 with Jupiter.     

Chemical dynamics of the CH(X2Π) + C2H4(X1A1g), CH(X2Π) + C2D4(X1A1g), and CD(X2Π) + C2H4(X1A1g) reactions studied under single collision conditions

The crossed beam reactions of the methylidyne radical with ethylene (CH(X(2)Π) + C(2)H(4)(X(1)A(1g))), methylidyne with D4-ethylene (CH(X(2)Π) + C(2)D(4)(X(1)A(1g))), and D1-methylidyne with ethylene (CD(X(2)Π) + C(2)H(4)(X(1)A(1g))) were conducted at nominal collision energies of 17-18 kJ mol(-1) to untangle the chemical dynamics involved in the formation of distinct C(3)H(4) isomers methylacetylene (CH(3)CCH), allene (H(2)CCCH(2)), and cyclopropene (c-C(3)H(4)) via C(3)H(5) intermediates. By tracing the atomic hydrogen and deuterium loss pathways, our experimental data suggest indirect scattering dynamics and an initial addition of the (D1)-methylidyne radical to the carbon-carbon double bond of the (D4)-ethylene reactant forming a cyclopropyl radical intermediate (c-C(3)H(5)/c-C(3)D(4)H/c-C(3)H(4)D). The latter was found to ring-open to the allyl radical (H(2)CCHCH(2)/D(2)CCHCD(2)/H(2)CCDCH(2)). This intermediate was found to be long lived with life times of at least five times its rotational period and decomposed via atomic hydrogen/deuterium loss from the central carbon atom (C2) to form allene via a rather loose exit transition state in an overall strongly exoergic reaction. Based on the experiments with partially deuterated reactants, no compelling evidence could be provided to support the formation of the cyclopropene and methylacetylene isomers under single collision conditions. Likewise, hydrogen/deuterium shifts in the allyl radical intermediates or an initial insertion of the (D1)-methylidyne radical into the carbon-hydrogen/deuterium bond of the (D4)-ethylene reactant were found to be-if at all-of minor importance. Our experiments propose that in hydrocarbon-rich atmospheres of planets and their moons such as Saturn's satellite Titan, the reaction of methylidyne radicals should lead predominantly to the hitherto elusive allene molecule in these reducing environments.     

Singlet and triplet potential surfaces for the O2+C2H4 reaction

Electronic structure calculations at the CASSCF and UB3LYP levels of theory with the aug-cc-pVDZ basis set were used to characterize structures, vibrational frequencies, and energies for stationary points on the ground state triplet and singlet O(2)+C(2)H(4) potential energy surfaces (PESs). Spin-orbit couplings between the PESs were calculated using state averaged CASSCF wave functions. More accurate energies were obtained for the CASSCF structures with the MRMP2/aug-cc-pVDZ method. An important and necessary aspect of the calculations was the need to use different CASSCF active spaces for the different reaction paths on the investigated PESs. The CASSCF calculations focused on O(2)+C(2)H(4) addition to form the C(2)H(4)O(2) biradical on the triplet and singlet surfaces, and isomerization reaction paths ensuing from this biradical. The triplet and singlet C(2)H(4)O(2) biradicals are very similar in structure, primarily differing in their C-C-O-O dihedral angles. The MRMP2 values for the O(2)+C(2)H(4)→C(2)H(4)O(2) barrier to form the biradical are 33.8 and 6.1 kcal/mol, respectively, for the triplet and singlet surfaces. On the singlet surface, C(2)H(4)O(2) isomerizes to dioxetane and ethane-peroxide with MRMP2 barriers of 7.8 and 21.3 kcal/mol. A more exhaustive search of reaction paths was made for the singlet surface using the UB3LYP/aug-cc-pVDZ theory. The triplet and singlet surfaces cross between the structures for the O(2)+C(2)H(4) addition transition states and the biradical intermediates. Trapping in the triplet biradical intermediate, following (3)O(2)+C(2)H(4) addition, is expected to enhance triplet→singlet intersystem crossing.     

Reaction dynamics of phenyl radicals (C6H5, X2A') with methylacetylene (CH3CCH(X1A1)), allene (H2CCCH2(X1A1)), and their D4-isotopomers

Crossed molecular beam experiments were utilized to untangle the reaction dynamics to form 1-phenylmethylacetylene [CH(3)CCC(6)H(5)] and 1-phenylallene [C(6)H(5)HCCCH(2)] in the reactions of phenyl radicals with methylacetylene and allene, respectively, over a range of collision energies from 91.4 to 161.1 kJ mol(-1). Both reactions proceed via indirect scattering dynamics and are initiated by an addition of the phenyl radical to the terminal carbon atom of the methylacetylene and allene reactants to form short-lived doublet C(9)H(9) collision complexes CH(3)CCHC(6)H(5) and C(6)H(5)H(2)CCCH(2). Studies with isotopically labeled reactants and the information on the energetics of the reactions depict that the energy randomization in the decomposing intermediates is incomplete. The collision complexes undergo atomic hydrogen losses via tight exit transition states leading to 1-phenylmethylacetylene [CH(3)CCC(6)H(5)] and 1-phenylallene [C(6)H(5)HCCCH(2)]. The possible role of both C(9)H(8) isomers as precursors to PAHs in combustion flames and in the chemistry of circumstellar envelopes of dying carbon stars is discussed.     

Crossed molecular beam studies of the reactions of allyl radicals, C3H5(X2A2), with methylacetylene (CH3CCH(X1A1)), allene (H2CCCH2(X1A1)), and their isotopomers

The chemical dynamics of the reaction of allyl radicals, C(3)H(5)(X(2)A(2)), with two C(3)H(4) isomers, methylacetylene (CH(3)CCH(X(1)A(1))) and allene (H(2)CCCH(2)(X(1)A(1))) together with their (partially) deuterated counterparts, were unraveled under single-collision conditions at collision energies of about 125 kJ mol(-1) utilizing a crossed molecular beam setup. The experiments indicate that the reactions are indirect via complex formation and proceed via an addition of the allyl radical with its terminal carbon atom to the terminal carbon atom of the allene and of methylacetylene (alpha-carbon atom) to form the intermediates H(2)CCHCH(2)CH(2)CCH(2) and H(2)CCHCH(2)CHCCH(3), respectively. The lifetimes of these intermediates are similar to their rotational periods but too short for a complete energy randomization to occur. Experiments with D4-allene and D4-methylacetylene verify explicitly that the allyl group stays intact: no hydrogen emission was observed but only the release of deuterium atoms from the perdeuterated reactants. Further isotopic substitution experiments with D3-methylacetylene combined with the nonstatistical nature of the reaction suggest that the intermediates decompose via hydrogen atom elimination to 1,3,5-hexatriene, H(2)CCHCH(2)CHCCH(2), and 1-hexen-4-yne, H(2)CCHCH(2)CCCH(3), respectively, via tight exit transition states located about 10-15 kJ mol(-1) above the separated products. The overall reactions were found to be endoergic by 98 +/- 4 kJ mol(-1) and have characteristic threshold energies to reaction between 105 and 110 kJ mol(-1). Implications of these findings to combustion and interstellar chemistry are discussed.     

Experimental and theoretical investigations of the inelastic and reactive scattering dynamics of O(3p) + D2

This paper presents a combined experimental and theoretical study of the dynamics of O((3)P) + D(2) collisions, with emphasis on a center-of-mass (c.m.) collision energy of 25 kcal mol(-1). The experiments were conducted with a crossed-molecular-beams apparatus, employing a laser detonation source to produce hyperthermal atomic oxygen and mass spectrometric detection to measure the product angular and time-of-flight distributions. The novel beam source, which enabled these experiments to be conducted, contributed unique challenges to the experiments and to the analysis, so the experimental methods and approach to the analysis are discussed in detail. Three different levels of theory were used: (1) quasiclassical trajectories (QCT), (2) time-independent quantum scattering calculations based on high-quality potential surfaces for the two lower-energy triplet states, and (3) trajectory-surface-hopping (TSH) studies that couple the triplet surfaces with the lowest singlet surface using a spin-orbit Hamiltonian derived from ab-initio calculations. The latter calculations explore the importance of intersystem crossing in the dynamics. Both experiment and theory show that inelastically scattered O atoms scatter almost exclusively in the forward direction, with little or no loss of translational energy. For the reaction, O((3)P) + D(2) --> OD + D, the experiment shows that, on average, approximately 50% of the available energy goes into product translation and that the OD product angular distributions are largely backward-peaked. These results may be interpreted in light of the QCT and TSH calculations, leading to the conclusion that the reaction occurs mainly on triplet potential energy surfaces with, at most, minor intersystem crossing to a singlet surface. Reaction on either of the two low-lying reactive triplet surfaces proceeds through a rebound mechanism in which the angular distributions are backward-peaked and the product OD is both vibrationally and rotationally excited. The quantum scattering results are in good agreement with QCT calculations, indicating that quantum effects are relatively small for this reaction at a collision energy of 25 kcal mol(-1).     

A theoretical study for the reaction of vinyl cyanide C2H3CN(X1A') with the ground state carbon atom C(3P) in cold molecular clouds

The reaction of the ground state atomic carbon, C(3P), with simple unsaturated nitrile, C2H3CN(X1A' (vinyl cyanide), is investigated theoretically to explore the probable routes for the formation of carbon-nitrogen-bearing species in extraterrestrial environments particularly of ultralow temperature. Five collision complexes without entrance barrier as a result of the carbon atom addition to the pi systems of C2H3CN are characterized. The B3YLP/6-311G(d,p) level of theory is utilized in obtaining the optimized geometries, harmonic frequencies, and energies of the intermediates, transition states, and products along the isomerization and dissociation pathways of each collision complex. Subsequently, with the facilitation of computed RRKM rate constants at collision energy of 0-10 kcal/mol, the most probable paths for each collision complexes are determined, of which the CCSD(T)/6-311G(d,p) energies are calculated. The major products predicted are exclusively due to the hydrogen atom dissociations, while the products of H2, CN, and CH2 decompositions are found negligible. Among many possible H-elimination products, cyano propargyl (p4) and 3-cyano propargyl (p5) are the most probable, in which p5 can be formed via two intermediates, cyano allene (i8) and cyano vinylmethylene (i6), while p4 is yielded from i8. The study suggests this class of reaction is an important route to the synthesis of unsaturated nitriles at the temperature as low as 10 K, and the results are valuable for future chemical models of interstellar clouds.     

Prediction of product branching ratios in the C(3P)+C2H2-->l-C3H+H/c-C3H+H/C3+H2 reaction using ab initio coupled clusters calculations extrapolated to the complete basis set combined with Rice-Ramsperger-Kassel-Marcus and radiationless transition theories

Ab initio CCSD(T) calculations of intermediates and transition states on the singlet and triplet C3H2 potential energy surfaces extrapolated to the complete basis set limit are combined with statistical computations of energy-dependent rate constants of the C(3P)+C2H2 reaction under crossed molecular beam conditions. Rice-Ramsperger-Kassel-Marcus theory is applied for isomerization and dissociation steps within the same multiplicity and radiationless transition and nonadiabatic transition state theories are used for singlet-triplet intersystem crossing rates. The calculated rate constants are utilized to predict product branching ratios. The results demonstrate that, in qualitative agreement with available experimental data, c-C3H+H and C3+H2 are the most probable products at low collision energies, whereas l-C3H+H becomes dominant at higher Ec above approximately 25 kJ/mol.     

Ab initio/Rice-Ramsperger-Kassel-Marcus study of the singlet C4H4 potential energy surface and of the reactions of C2(X1 Sigmag+) with C4H4(X1A1g) and C(1D) with C3H4 (allene and methylacetylene)

Ab initio modified Gaussian-2 G2M(RCC,MP2) calculations have been performed for various isomers and transition states on the singlet C4H4 potential energy surface. The computed relative energies and molecular parameters have then been used to calculate energy-dependent rate constants for different isomerization and dissociation processes in the C4H4 system employing Rice-Ramsperger-Kassel-Marcus theory and to predict branching ratios of possible products of the C2(1Sigmag+)+C2H4, C(1D)+H2CCCH2, and C(1D)+H3CCCH reactions under single-collision conditions. The results show that C2 adds to the double C=C bond of ethylene without a barrier to form carbenecyclopropane, which then isomerizes to butatriene by a formal C2 "insertion" into the C-C bond of the C2H4 fragment. Butatriene can rearrange to the other isomers of C4H4, including allenylcarbene, methylenecyclopropene, vinylacetylene, methylpropargylene, cyclobutadiene, tetrahedrane, methylcyclopropenylidene, and bicyclobutene. The major decomposition products of the chemically activated C4H4 molecule formed in the C2(1Sigmag+)+C2H4 reaction are calculated to be acetylene+vinylidene (48.6% at Ecol = 0) and 1-buten-3-yne-2-yl radical [i-C4H3(X2A'), H2C=C=C=CH*]+H (41.3%). As the collision energy increases from 0 to 10 kcal/mol, the relative yield of i-C4H3+H grows to 52.6% and that of C2H2+CCH2 decreases to 35.5%. For the C(1D)+allene reaction, the most important products are also i-C4H3+H (55.2%) and C2H2+CCH2 (30.1%), but for C(1D)+methylacetylene, which accesses a different region of the C4H4 singlet potential energy surface, the calculated product branching ratios differ significantly: 65%-69% for i-C4H3+H, 18%-14% for C2H2+CCH2, and approximately 8% for diacetylene+H2.     

Potential energy surface and product branching ratios for the reaction of dicarbon, C(2)(X(1)Sigma(g)(+)), with methylacetylene, CH(3)CCH(X(1)A(1)): 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.     

Gas‐Phase Synthesis of the Benzyl Radical (C6H5CH2)

Dicarbon (C₂), the simplest bare carbon molecule, is ubiquitous in the interstellar medium and in combustion flames. A gas‐phase synthesis is presented of the benzyl radical (C₆H₅CH₂) by the crossed molecular beam reaction of dicarbon, C₂(X¹Σ_g⁺, a³Π_u), with 2‐methyl‐1,3‐butadiene (isoprene; C₅H₈; X¹A′) accessing the triplet and singlet C₇H₈ potential energy surfaces (PESs) under single collision conditions. The experimental data combined with ab initio and statistical calculations reveal the underlying reaction mechanism and chemical dynamics. On the singlet and triplet surfaces, the reactions involve indirect scattering dynamics and are initiated by the barrierless addition of dicarbon to the carbon–carbon double bond of the 2‐methyl‐1,3‐butadiene molecule. These initial addition complexes rearrange via multiple isomerization steps, leading eventually to the formation of C₇H₇ radical species through atomic hydrogen elimination. The benzyl radical (C₆H₅CH₂), the thermodynamically most stable C₇H₇ isomer, is determined as the major product.