The Recombination of Propargyl Radicals: Solving the Master Equation

Kinetics of enol formation from reaction of OH with propene

Kinetics of enol generation from propene has been predicted in an effort to understand the presence of enols in flames. A potential energy surface for reaction of OH with propene was computed by CCSD(T)/cc-pVDZ//B3LYP/cc-pVTZ calculations. Rate constants of different product channels and branching ratios were then calculated using the Master Equation formulation (J. Phys. Chem. A 2006, 110, 10528). Of the two enol products, ethenol is dominant over propenol, and its pathway is also the dominant pathway for the OH + propene addition reactions to form bimolecular products. In the temperature range considered, hydrogen abstraction dominated propene + OH consumption by a branching ratio of more than 90%. Calculated rate constants of enol formation were included in the Utah Surrogate Mechanism to model the enol profile in a cyclohexane premixed flame. The extended model shows consistency with experimental data and gives 5% contribution of ethenol formation from OH + propene reaction, the rest coming from ethene + OH.

The simulated photoelectron spectrum of 1-propynide

The negative ion photoelectron spectrum of 1-propynide is computed by employing the multimode vibronic coupling approach. A three-state quasidiabatic Hamiltonian, H(d), is reported, which accurately represents the ab initio determined equilibrium geometries and harmonic frequencies of the ground X (2)A(1) state as well as the low-lying Jahn-Teller distorted components of the A (2)E excited state. It also reproduces both the minimum energy crossing point (MECP) on the symmetry-required (2)E(x)-(2)E(y) conical intersection seam and the MECP on the same symmetry (2)A(1)-(2)E(x) conical intersection seam. H(d) includes all terms through second order in internal coordinates for both the diagonal and off-diagonal blocks. It is centered at the (2)E(x)-(2)E(y) MECP and is determined using ab initio gradients and derivative couplings near both the (2)E(x)-(2)E(y) MECP and the X (2)A(1) equilibrium geometry. This construction is enabled by a recently reported normal equation based algorithm. The C(3v) symmetry of the system is used to significantly reduce the computational cost of the ab initio treatment. This H(d) is then expressed in a vibronic basis that is chosen for its ability to reduce the dimension of the vibronic expansion. The vibronic Hamiltonian matrix is diagonalized to obtain a negative ion photoelectron spectrum for 1-propynide-h(3). The determined spectrum compares favorably with previous spectroscopic results. In particular, the lines attributable to the (2)E state are found to be much weaker than those corresponding to the (2)A(1) state of 1-propynyl. This diminution of the (2)E state is attributable principally to the (2)E(x)-(2)A(1) conical intersection rather than an intrinsically small electronic transition moment for the production of the (2)E state.

On the Combination Reactions of Hydrogen Atoms with Resonance-Stabilized Hydrocarbon Radicals

Procedures for accurately predicting the kinetics of H atom associations with resonance stabilized hydrocarbon radicals are described and applied to a series of reactions. The approach is based on direct CASPT2/cc-pvdz evaluations of the orientation dependent interaction energies within variable reaction coordinate transition state theory. One-dimensional corrections to the interaction energies are estimated from a CASPT2/aug-cc-pvdz minimum energy path (MEP) on the specific reaction of interest and a CASPT2/aug-cc-pvtz MEP for the H + CH3 reaction. A dynamical correction factor of 0.9 is also applied. For the H + propargyl, allyl, cyclopentadienyl, and benzyl reactions, where the experimental values appear to be quite well determined, theory and experiment agree to within their error bars. Predictions are also made for the combinations with triplet propargylene, CH2CCCH, CH3CCCH2, CH2CHCCH2, CH3CHCCH, cyclic-C4H5, CH2CCCCH, and CHCCHCCH.

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

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

Photodissociation of benzene under collision-free conditions: an ab initio/Rice-Ramsperger-Kassel-Marcus study

The ab initio/Rice-Ramsperger-Kassel-Marcus (RRKM) approach has been applied to investigate the photodissociation mechanism of benzene at various wavelengths upon absorption of one or two UV photons followed by internal conversion into the ground electronic state. Reaction pathways leading to various decomposition products have been mapped out at the G2M level and then the RRKM and microcanonical variational transition state theories have been applied to compute rate constants for individual reaction steps. Relative product yields (branching ratios) for C(6)H(5)+H, C(6)H(4)+H(2), C(4)H(4)+C(2)H(2), C(4)H(2)+C(2)H(4), C(3)H(3)+C(3)H(3), C(5)H(3)+CH(3), and C(4)H(3)+C(2)H(3) have been calculated subsequently using both numerical integration of kinetic master equations and the steady-state approach. The results show that upon absorption of a 248 nm photon dissociation is too slow to be observable in molecular beam experiments. In photodissociation at 193 nm, the dominant dissociation channel is H atom elimination (99.6%) and the minor reaction channel is H(2) elimination, with the branching ratio of only 0.4%. The calculated lifetime of benzene at 193 nm is about 11 micros, in excellent agreement with the experimental value of 10 micros. At 157 nm, the H loss remains the dominant channel but its branching ratio decreases to 97.5%, while that for H(2) elimination increases to 2.1%. The other channels leading to C(3)H(3)+C(3)H(3), C(5)H(3)+CH(3), C(4)H(4)+C(2)H(2), and C(4)H(3)+C(2)H(3) play insignificant role but might be observed. For photodissociation upon absorption of two UV photons occurring through the neutral "hot" benzene mechanism excluding dissociative ionization, we predict that the C(6)H(5)+H channel should be less dominant, while the contribution of C(6)H(4)+H(2) and the C(3)H(3)+C(3)H(3), CH(3)+C(5)H(3), and C(4)H(3)+C(2)H(3) radical channels should significantly increase.

Crossed-beam radical-radical reaction dynamics of O(3P)+C3H3-->H(2S)+C3H2O

The radical-radical oxidation reaction, O(3P)+C3H3 (propargyl)-->H(2S)+C3H2O (propynal), was investigated using vacuum-ultraviolet laser-induced fluorescence spectroscopy in a crossed-beam configuration, together with ab initio and statistical calculations. The barrierless addition of O(3P) to C3H3 is calculated to form energy-rich addition complexes on the lowest doublet potential energy surface, which subsequently undergo direct decomposition steps leading to the major reaction products, H+C3H(2)O (propynal). According to the nascent H-atom Doppler-profile analysis, the average translational energy of the products and the fraction of the average transitional energy to the total available energy were determined to be 5.09+/-0.36 kcal/mol and 0.077, respectively. On the basis of a comparison with statistical prior calculations, the reaction mechanism and the significant internal excitation of the polyatomic propynal product can be rationalized in terms of the formation of highly activated, short-lived addition-complex intermediates and the adiabaticity of the excess available energy along the reaction coordinate.

Conformation-specific spectroscopy of 3-benzyl-1,5-hexadiyne

3-Benzyl-1,5-hexadiyne (BHD) was studied by a combination of methods, including resonance-enhanced-two-photon ionization, UV-UV hole-burning spectroscopy, resonant ion-dip infrared spectroscopy, and rotational band contour analysis. There are five conformations of BHD observed in the expansion with their 1<-- S0 origins occurring at 37520, 37565, 37599, 37605, and 37631 cm(-1). DFT calculations predict six low energy conformations. Conformational assignments have been made by comparison of the experimental infrared spectra in the alkyl and acetylenic CH stretch region to DFT vibrational frequency and infrared intensity calculations. Rotational band contours provided further confirmation of these assignments. The electronic origin shifts of BHD compare favorably to the electronic origin shifts of 5-phenyl-1-pentyne with the exception of one conformation. This conformation is unique in that it is the only structure with both acetylenic groups in the gauche position over the ring. This gauche-gauche conformation exhibits a perpendicular (b-type) transition and produces extensive vibronic coupling reminiscent of symmetric monosubstituted benzenes.

Isomeric Product Distributions from the Self-Reaction of Propargyl Radicals

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

Isomer-Specific Spectroscopy and Conformational Isomerization Energetics ofo-,m-, andp-Ethynylstyrenes

The infrared and ultraviolet spectroscopy of o-, m-, and p-ethynylstyrene isomers (oES, mES, and pES) were studied by a combination of methods, including resonance-enhanced two-photon ionization (R2PI), UV-UV hole-burning spectroscopy (UVHB), resonant ion-dip infrared spectroscopy (RIDIRS), and rotationally resolved fluorescence excitation spectroscopy. In addition, the newly developed method of stimulated emission pumping-population transfer spectroscopy (SEP-PTS) was used to determine the energy threshold to conformational isomerization in m-ethynylstyrene. The S(1) <-- S(0) origin transitions of oES and pES occur at 32 369 and 33 407 cm(-1), respectively. In mES, the cis and trans conformations are calculated to be close in energy. In the R2PI spectrum of mES, the two most prominent peaks (32672 and 32926 cm(-1)) were confirmed by UVHB spectroscopy to be S(1) <-- S(0) origins of these two conformers. The red-shifted conformer was identified as the cis structure by least-squares fitting of the rotationally resolved fluorescence excitation spectrum of the origin band. There are also two possible conformations in oES, but transitions due to only one were observed experimentally, as confirmed by UVHB spectroscopy. Density functional theory calculations (B3LYP/6-31+G) predict that the cis-ortho conformer, in which the substituents point toward each other, is about 8 kJ/mol higher in energy than the trans-ortho isomer, and should only be about 5% of the room temperature population of oES. Ground-state infrared spectra in the C-H stretch region (3000-3300 cm(-1)) of each isomer were obtained with RIDIRS. In all three structural isomers, the acetylenic C-H stretch fundamental was split by Fermi resonance. Infrared spectra were also recorded in the excited electronic state, using a UV-IR-UV version of RIDIR spectroscopy. In all three isomers the acetylenic C-H stretch fundamental was unshifted from the ground state, but no Fermi resonance was seen. The first observed and last unobserved transitions in the SEP-PT spectrum were used to place lower and upper bounds on the barrier to cis --> trans isomerization in m-ethynylstyrene of 990-1070 cm(-1). Arguments are given for the lack of a kinetic shift in the measurement. The analogous trans --> cis barrier is in the same range (989-1065 cm(-1)), indicating that the relative energies of the zero-point levels of the two isomers are (E(ZPL)(cis) - E(ZPL)(trans))= -75 to +81 cm(-1). Both the barrier heights and relative energies of the minima are close to those determined by DFT (Becke3LYP/6-31+G) calculations.

Ab Initio and Direct Dynamics Studies of the Reaction of Singlet Methylene with Acetylene and the Lifetime of the Cyclopropene Complex

The energetics of the (1)CH(2) + C(2)H(2) --> H + C(3)H(3) reaction are accurately calculated using an extrapolated coupled-cluster/complete basis set (CBS) method based on the cc-pVDZ, cc-pVTZ, and cc-pVQZ basis sets. The reaction enthalpy (0 K) is predicted to be -20.33 kcal/mol. This reaction has no classical barrier in either the entrance or exit channel. However, there are several stable intermediates-cyclopropene (c-C(3)H(4)), allene (CH(2)CCH(2)), and propyne (CH(3)CCH)-along the minimum energy path. These intermediates with zero-point energy corrections lie below the reactants by 87.11 (c-C(3)H(4)), 109.69 (CH(2)CCH(2)), and 110.78 kcal/mol (CH(3)CCH). The vibrationally adiabatic ground-state (VAG) barrier height for c-C(3)H(4) isomerization to allene is obtained as 45.2 kcal/mol, and to propyne as 37.2 kcal/mol. In addition, the (1)CH(2) + C(2)H(2) reaction is investigated utilizing the dual-level "scaling all correlation" (SAC) ab initio method of Truhlar et al., i.e., the UCCSD(SAC)/cc-pVDZ theory. Results show that the reaction occurs via long-lived complexes. The lifetime of the cyclopropene intermediate is obtained as 3.2 +/- 0.4 ps. It is found that the intermediate propyne can be formed directly from reactants through the insertion of (1)CH(2) into a C-H bond of C(2)H(2). However, compared to the major mechanism in which the propyne is produced through a ring-opening of the cyclopropene complex, this reaction pathway is much less favorable. Finally, the theoretical thermal rate constant exhibits a negative temperature dependence, which is in excellent agreement with the previous results. The temperature dependence is consistent with the earlier RRKM results but weaker than the experimental observations at high temperatures.

Thermochemistry of disputed soot formation intermediates C4H3 and C4H5

Accurate isomeric energy differences and standard enthalpies of formation for disputed intermediates in soot formation, C(4)H(3) and C(4)H(5), have been determined through systematic extrapolations of ab initio energies. Electron correlation has been included through second-order Z-averaged perturbation theory (ZAPT2), and spin-restricted, open-shell coupled-cluster methods through triple excitations [ROCCSD, ROCCSD(T), and ROCCSDT] utilizing the correlation-consistent hierarchy of basis sets, cc-pVXZ (X = D, T, Q, 5, and 6), followed by extrapolations to the complete basis set limit via the focal point method of Allen and co-workers. Reference geometries were fully optimized at the ROCCSD(T) level with a TZ(2d1f,2p1d) basis set. Our analysis finds that the resonance-stabilized i-C(4)H(3) and i-C(4)H(5) isomers lie 11.8 and 10.7 kcal mol(-1) below E-n-C(4)H(3) and E-n-C(4)H(5), respectively, several kcal mol(-1) (more, less) than reported in recent (diffusion Monte Carlo, B3LYP density-functional) studies. Moreover, in these systems Gaussian-3 (G3) theory suffers from large spin contamination in electronic wave functions, poor reference geometries, and anomalous vibrational frequencies, but fortuitous cancellation of these sizable errors leads to isomerization energies apparently accurate to 1 kcal mol(-1). Using focal-point extrapolations for isodesmic reactions, we determine the enthalpies of formation (delta(f)H(0) (composite function)) for i-C(4)H(3), Z-n-C(4)H(3), E-n-C(4)H(3), i-C(4)H(5), Z-n-C(4)H(5), and E-n-C(4)H(5) to be 119.0, 130.8, 130.8, 78.4, 89.7, and 89.1 kcal mol(-1), respectively. These definitive values remove any remaining uncertainty surrounding the thermochemistry of these isomers in combustion models, allowing for better assessment of whether even-carbon pathways contribute to soot formation.

Crossed beams study of the reaction 1CH2+C2H2→C3H3+H

The reaction of electronically excited singlet methylene (1CH2) with acetylene (C2H2) was studied using the method of crossed molecular beams at a mean collision energy of 3.0 kcal/mol. The angular and velocity distributions of the propargyl radical (C3H3) products were measured using single photon ionization (9.6 eV) at the advanced light source. The measured distributions indicate that the mechanism involves formation of a long-lived C3H4 complex followed by simple C-H bond fission producing C3H3+H. This work, which is the first crossed beams study of a reaction involving an electronically excited polyatomic molecule, demonstrates the feasibility of crossed molecular beam studies of reactions involving 1CH2.

A combined crossed beam and theoretical investigation of O(3P)+C3H3→C3H2+OH

The radical-radical reaction dynamics of ground-state atomic oxygen [O(3P)] with propargyl radicals (C3H3) has first been investigated in a crossed beam configuration. The radical reactants O(3P) and C3H3 were produced by the photodissociation of NO2 and the supersonic flash pyrolysis of precursor propargyl bromide, respectively. A new exothermic channel of O(3P) + C3H3 --> C3H2 + OH was identified and the nascent distributions of the product OH in the ground vibrational state (X 2Pi:nu" = 0) showed bimodal rotational excitations composed of the low- and high-N" components without spin-orbit propensities. The averaged ratios of Pi(A')/Pi(A") were determined to be 0.60 +/- 0.28. With the aid of ab initio theory it is predicted that on the lowest doublet potential energy surface, the reaction proceeds via the addition complexes formed through the barrierless addition of O(3P) to C3H3. The common direct abstraction pathway through a collinear geometry does not occur due to the high entrance barrier in our low collision energy regime. In addition, the major reaction channel is calculated to be the formation of propynal (CHCCHO) + H, and the counterpart C3H2 of the probed OH product in the title reaction is cyclopropenylidene (1c-C3H2) after considering the factors of barrier height, reaction enthalpy and structural features of the intermediates formed along the reaction coordinate. On the basis of the statistical prior and rotational surprisal analyses, the ratio of population partitioning for the low- and high-N" is found to be about 1:2, and the reaction is described in terms of two competing addition-complex mechanisms: a major short-lived dynamic complex and a minor long-lived statistical complex. The observed unusual reaction mechanism stands in sharp contrast with the reaction of O(3P) with allyl radical (C3H5), a second significant conjugated hydrocarbon radical, which shows totally dynamic processes [J. Chem. Phys. 117, 2017 (2002)], and should be understood based upon the characteristic electronic structures and reactivity of the intermediates on the potential energy surface.