An ab initio G3-type/statistical theory study of the formation of indene in combustion flames. II. The pathways originating from reactions of cyclic C5 species-cyclopentadiene and cyclopentadienyl radicals

Chemically accurate ab initio Gaussian-3-type calculations of various rearrangements on the C10H11 potential energy surface have been performed to investigate the indene formation mechanism originating from the reactions of two abundant cyclic C5 species, cyclopentadiene and cyclopentadienyl radicals. Using the accurate ab initio data, statistical theory calculations have been applied to obtain high-pressure-limit thermal rate constants within the 300-3000 K temperature range, followed by calculations of relative product yields. Totally, 12 reaction pathways leading to indene and several azulene precursors, 1,5-, 1,7-, 1,8a-, and 1,3a-dihydroazulene, have been mapped out, and the relative contributions of each pathway to the formation of reaction products have been estimated. At temperatures relevant to combustion, the indene has been found as the major reaction product (>50%) followed by 1,5-dihydroazulene (25-35%), whereas all other products demonstrate either minor or negligible yields. The results of the present study have been combined with our previous data for rearrangements of the 9-H-fulvalenyl radical on the C10H9 potential energy surface to draw the detailed picture of radical-promoted reaction mechanisms leading from c-C5 species to the production of indene, naphthalene, azulene, and fulvalene in combustion. The suggested mechanism and computed product yields are consistent with the experimental data obtained in the low-temperature pyrolysis of cyclopentadiene, where indene and naphthalene have been found as the major reaction products.     

Ab initio study of chain branching reactions involving second generation products in hydrocarbon combustion mechanisms

sec-Alkyl radicals are key reactive intermediates in the hydrocarbon combustion and atmospheric decomposition mechanisms that are formed by the abstraction of hydrogen from an alkane, or as a second generation product of n-alkyl H-migrations, C-C bond scissions in branched alkyl radicals, or the bimolecular reaction between olefins and n-alkyl radicals. Since alkanes and branched alkanes, which the sec-alkyl radicals are derived from, make up roughly 40-50% of traditional fuels an understanding of their chemistry is essential to improving combustion systems. The present work investigates all H-migration reactions initiated from an sec-alkyl radical that involve the movement of a secondary hydrogen, for the 2-butyl through 4-octyl radicals, using the CBS-Q, G2, and G4 composite methods. The resulting thermodynamic and kinetic parameters are compared to similar reactions in n-alkyl radicals in order to determine underlying trends. Particular attention is paid to the effect of cis/trans and 1,3-diaxial interactions on activation energies and rate coefficients. When combined with our previous work on n-alkyl radical H-migrations, a complete picture of H-migrations in unbranched alkyl radicals is obtained. This full data set suggests that the directionality of the remaining branched chains has a minimal effect on the rate coefficients for all but the largest viable transition states, which is in stark contrast to the differences predicted by the structurally similar dimethylcycloalkanes. In fact the initial location of the secondary radical site has a greater effect on the rate than does the directionality of the remaining alkyl chains. The activation energies for secondary to secondary reactions are much closer to those of the secondary to primary H-migrations. However, the rate coefficients are found to be closer to the corresponding primary to primary reaction values. A significant ramification of these results is that there will be multiple viable reaction pathways for these reactions instead of only one dominant pathway as previously believed.     

Reaction of phenyl radical with propylene as a possible source of indene and other polycyclic aromatic hydrocarbons: an ab initio/RRKM-ME study

Ab initio G3(MP2,CC)//B3LYP/6-311G** calculations have been performed to investigate the potential energy surface (PES) and mechanism of the reaction of phenyl radical with propylene followed by kinetic RRKM-ME calculations of rate constants and product branching ratios at various temperatures and pressures. The reaction can proceed either by direct hydrogen abstraction producing benzene and three C(3)H(5) radicals [1-propenyl (CH(3)CHCH), 2-propenyl (CH(3)CCH(2)), and allyl (CH(2)CHCH(2))] or by addition of phenyl to the CH or CH(2) units of propylene followed by rearrangements on the C(9)H(11) PES producing nine different products after H or CH(3) losses. The H abstraction channels are found to be kinetically preferable at temperatures relevant to combustion and to contribute 55-75% to the total product yield in the 1000-2000 K temperature range, with the allyl radical being the major product (~45%). The relative contributions of phenyl addition channels are calculated to be ~35% at 1000 K, decreasing to ~15% at 2000 K, with styrene + CH(3) and 3-phenylpropene + H being the major products. Collisional stabilization of C(6)H(5) + C(3)H(6) addition complexes is computed to be significant only at temperatures up to 1000-1200 K, depending on the pressure, and maximizes at low temperatures of 300-700 K reaching up to 90% of the total product yield. At T > 1200 K collisional stabilization becomes negligible, whereas the dissociation products, styrene plus methyl and 3-phenylpropene + H, account for up to 45% of the total product yield. The production of bicyclic aromatic species including indane C(9)H(10) is found to be negligible at all studied conditions indicating that the phenyl addition to propylene cannot be a source of polycyclic aromatic hydrocarbons (PAH) on the C(9)H(11) PES. Alternatively, the formation of a PAH molecule, indene C(9)H(8), can be accomplished through secondary reactions after activation of a major product of the C(6)H(5) + C(3)H(6) addition reaction, 3-phenylpropene, by direct hydrogen abstraction by small radicals, such as H, OH, CH(3), etc. It is shown that at typical combustion temperatures 77-90% of C(9)H(9) radicals formed by H-abstraction from 3-phenylpropene undergo a closure of a cyclopentene ring via low barriers and then lose a hydrogen atom producing indene. This results in 7.0-14.5% yield of indene relative to the initial C(6)H(5) + C(3)H(6) reactants within the 1000-2000 K temperature range.     

Ab initio study of the thermodynamic and kinetic stability of the ammonium radical

An increasing interest in the possible existence of the NH₄ radical has emerged in recent years. In this paper we report an ab initio UHF CI study of the ammonium radical, an investigation of parts of the energy surface around NH₄ and a theoretical prediction of the kinetic parameters of the radical formation and dissociation reactions within the framework of the TST theory. The ground state of the ammonium radical appears to be of the Rydberg type. Its ionization potential is found to be 4.29 eV. The NH₄ formation reaction from NH₂ + H₂ is very slightly exothermic whereas the reaction from NH₃ + H is slightly endothermic. We find a transition state of C_3v symmetry for the dissociation of NH₄ into NH₃ + H. The insertion of H₂ into NH₂ occurs according to a two-step mechanism whose determining step corresponds to the crossing of a saddle point with C_s symmetry previously obtained in the study of the reaction NH₂ + H₂ → NH₃ + H. Finally, we predict for NH₄ and ND₄ lifetimes of 0.1 and 1.4 μs respectively.     

Ab initio study of water clustering in the presence of a methyl radical

The geometrical structure and binding energy of small clusters of methyl radical and water molecules (up to five water molecules) in gas phase and water media have been investigated at the MP2 level of theory using 6-311++G(2df,2p) basis set. The complexes characterized contain OH···O, CH···O, and OH···C attractive interactions with stabilization energies in the range 6–143 kJ mol^?1. The solvent has an enhancing influence on the stabilities of studied clusters. The atoms in molecules theory were also applied to explain the nature of the complexes. The interaction energies have been partitioned with the natural energy decomposition analysis showing that the most important attractive term corresponds to the charge transfer one.     

Ab initio scf and Ci study of the electronic structure of the trichlorine radical

Ab initio SCF and CI calculations using a double-zeta plus polarization basis set have been carried out on the trichlorine radical Cl₃ to determine its electronic structure. The minimum in energy is determined for a bent structure at a bond angle of 146° and bond lengths of 2.18 Å (SCF ) or 2.22 Å (CI ). At linear geometry a ²Π_u state is found to be lowest, approximately 7 kcal above the bent minimum, followed by a ²∑_g⁺ state, which is around 4 kcal higher. This situation suggests that already for low quantum numbers a complex vibrational pattern in the Cl₃ infrared spectrum is to be expected due to spin-orbit coupling as well as coupling of electronic, vibrational, and rotational motion.     

Product branching ratios in photodissociation of phenyl radical: a theoretical ab initio/Rice-Ramsperger-Kassel-Marcus study

Ab initio CCSD(T)/CBS//B3LYP/6-311G** calculations of the potential energy surface for possible dissociation channels of the phenyl radical are combined with microcanonical Rice-Ramsperger-Kassel-Marcus calculations of reaction rate constants in order to predict statistical product branching ratios in photodissociation of c-C(6)H(5) at various wavelengths. The results indicate that at 248 nm the photodissociation process is dominated by the production of ortho-benzyne via direct elimination of a hydrogen atom from the phenyl radical. At 193 nm, the statistical branching ratios are computed to be 63.4%, 21.1%, and 14.4% for the o-C(6)H(4) + H, l-C(6)H(4) ((Z)-hexa-3-ene-1,5-diyne) + H, and n-C(4)H(3) + C(2)H(2) products, respectively, in a contradiction with recent experimental measurements, which showed C(4)H(3) + C(2)H(2) as the major product. Although two lower energy pathways to the i-C(4)H(3) + C(2)H(2) products are identified, they appeared to be kinetically unfavorable and the computed statistical branching ratio of i-C(4)H(3) + C(2)H(2) does not exceed 1%. To explain the disagreement with experiment, we optimized conical intersections between the ground and the first excited electronic states of C(6)H(5) and, based on their structures and energies, suggested the following photodissociation mechanism at 193 nm: c-C(6)H(5) 1 → absorption of a photon → electronically excited 1 → internal conversion to the lowest excited state → conversion to the ground electronic state via conical intersections at CI-2 or CI-3 → non-statistical decay of the vibrationally excited radical favoring the formation of the n-C(4)H(3) + C(2)H(2) products. This scenario can be attained if the intramolecular vibrational redistribution in the CI-2 or CI-3 structures in the ground electronic state is slower than their dissociation to n-C(4)H(3) + C(2)H(2) driven by the dynamical preference.     

Ab initio MRD CI calculations for the electron spectrum of the C3 radical

Calculations using the MRD CI method are reported for the ground and low lying excited states of C₃. Transitions from the 3σ_u, 4σ_g and 1π_u MO's into 1π_g are considered, as well as the 1π_u → 3s Rydberg species and the corresponding ionization, and good agreement with experimental data is obtained where comparison is possible. Potential curves calculated for the ground and (1π_u → 1π_g) ¹Σ⁺_u excited state are discussed.     

Ab initio study of the inversion barrier in NF3+

Ab initio Hartree–Fock, Møller–Plesset perturbation theory (MP 2), and quadratic configuration interaction, using single and double substitutions (QCISD ), calculations were carried out for the NF₃⁺ ion. Optimized structures were examined at the various levels of theory. Calculation of the inversion barrier height shows the importance of optimizing the geometry at the post-Hartree–Fock level and the inclusion of polarization functions. The best calculated inversion barrier was 13.3 kcal/mol, compared to an experimental value of 17.3 kcal/mol. The dissociation transition state was computed to determine the well depth of the NF₃⁺ ion and its stability toward dissociation. The computed well depth was 28 and 48 kcal/mol at the SCF and MP 2 levels, respectively. © 1994 John Wiley & Sons, Inc.     

Ab initio scf and CI study of the ground and excited states of the HN2 radical

Non-empirical SCF and CI calculations are reported for the HN₂, free radical in various low-lying electronic states. The nature of the angular and N-N and N-H stretching potential curves of each of these species is investigated, including a study of the dissociative behavior of such states. The ground state is found to be only very slightly bound with respect to NH stretch, in contrast to what is observed for isoelectronic HCO, The vertical electronic spectrum of HN₂, appears to be marked by a single long wavelength transition (1.95 eV) from the bent (124°) ‘A’ ground state to the linear ²Π excited species, but at least four other intra-valence and an additional n → 3s Rydberg species are indicated in the 5.5–8.0 eV absorbing region.     

Ab initio study of the formation of C3H3+ from the reaction of CH3+ with acetylene

Ab initio molecular orbital theory has been used to study the mechanism of the formation of C₃H₃⁺ from the reaction of CH₃⁺ with acetylene. The highest level geometry optimizations and frequencies were computed at MP2-FC/6-31G**; single point energies of all the critical structures were computed to the MP4-FC/6-31G**//MP2-FC/6-31G** theory level. One of the three alternative transition structures leading to the formation of C₃H₃⁺ gives the cyclopropenyl cation and the other two the propargyl cation. The proportions of C₃H₂D⁺ and C₃HD₂⁺ obtained when CD₃⁺ reacts with acetylene, and the composite nature of the metastable peak observed for the [C₃H₅]⁺→[C₃H₃]⁺ + H₂ fragmentation are explained by assuming a different degree of deuterium scrambling depending on the energy of the system. © 1996 by John Wiley & Sons, Inc.     

Ab initio scf and ci study of the electronic structure of BH3CO

Results of ab initio SCF and CI calculations employing a Gaussian basis set of double-zeta quality are reported for BH₃CO. The heat of formation for the gas-phase reaction, BH₃ + CO → BH₃CO, is calculated as ?10.98 kcal mol^?1 within the SCF approximation, and as ?14.56 kcal mol^?1 if the CI treatment is included. This is in good agreement with the estimated experimental value of ?16.6 kcal mol^?1. The energy of rearrangement of the BH₃ fragment from D_3h to C_3v symmetry in BH₃CO is calculated as 15.97 kcal mol^?1. Molecular properties have been studied in terms of the calculated electron populations, the dipole moment and the electric Field gradient of ¹¹B in BH₃CO.     

Ab initio calculation of the heats of formation of CF3OH and CF2O

The heats of formation of the trifluoromethanol (CF₃OH) and carbonic difluoride (CF₂O) molecules are calculated with theoretical methods of demonstrated high accuracy. The results are used to assess the accuracy of the reported value of the heat of formation of CF₂O and to provide an accurate estimate for CF₃OH. The implications of these results for atmospheric chemistry are discussed.     

Ab initio MO Calculations of benzene + TCNE and naphthalene + TCNE complexes with STO-3G π-split basis set

Ab initio SCF MO calculations have been carried out on benzene + TCNE (tetracyanoethylene) and naphthalene + TCNE complexes with the STO -3G, STO -3G π-split (STO -3G for π orbitals and a split basis for π orbitals), and 4–31G basis sets. The interaction energy, gross charges, dipole moment, and the electron density in the middle plane of the complexes have also been evaluated. The STO -3G π-split basis set is appropriate for the calculation of large π–π stacking complexes from two points of view, production of reliable results and ease of computations. The approximation scheme based on the semiorthogonalized orgitals is revealed to be very efficient to save CPU time and storage in such calculations. The stable conformation and the charge-transfer interaction of the two complexes are discussed on the basis of the calculated quantities.     

Ab initio/density functional theory and multichannel RRKM study for the ClO + CH2O reaction

The potential energy surface for the CH(2)O + ClO reaction was calculated at the QCISD(T)/6-311G(2d,2p)//B3LYP/6-311G(d,p) level of theory. The rate constants for the lower barrier reaction channels producing HOCl + HCO, H atom, OCH(2)OCl, cis-HC(O)OCl and trans-HC(O)OCl have been calculated by TST and multichannel RRKM theory. Over the temperature range of 200-2000 K, the overall rate constants were k(200-2000K) = 1.19 x 10(-13)T(0.79) exp(-3000.00/T). At 250 K, the calculated overall rate constant was 5.80 x 10(-17) cm(3) molecule(-1) s(-1), which was in good agreement with the experimental upper limit data. The calculated results demonstrated that the formation of HOCl + HCO was the dominant reaction channel and was exothermic by 9.7 kcal/mol with a barrier of 5.0 kcal/mol. When it retrograded to the reactants CH(2)O + ClO, an energy barrier of 14.7 kcal/mol is required. Furthermore, when HOCl decomposed into H + ClO, the energy required was 93.3 kcal/mol. These results suggest that the decomposition in both the forward and backward directions for HOCl would be difficult in the ground electronic state.     

Ab initio study on the structures of fluorinated formates and hydrogen abstraction reaction with OH radical

The conformational potential energy surfaces for mono- and difluoromethyl formate have been determined by using a modified G2(MP2) level of calculations. The structures and vibrational frequencies for the conformers of mono- and difluoromethyl formate have been reported. The hydrogen abstraction reaction channels between these two formates and OH radicals have been studied at the same level of theory. Using the standard transition state theory and taking into account the effect of tunneling across the reaction barrier, we have estimated the rate constant for hydrogen abstraction by OH radical. The effect of successive fluorine substitution for methyl hydrogen on the conformational stability and on the hydrogen abstraction rate has been analyzed.     

Ab initio UHF and SWX α study of the GeF3 and SiF3 radicals

Results of ab initio UHF and SWX _α calculations concerning the geometrical and electronic properties of the GeF₃ and SiF₃ radicals are presented. Comparison of the theoretical estimates of inversion barriers, ionization potentials and electron affinities obtained for three XF₃ radicals (X=Ge, Si, C) is made.     

Theoretical study of the dynamics of Cl + O3 reaction I. Ab initio potential energy surface and quasiclassical trajectory results

We present a global full dimensional potential energy surface (PES) for the Cl + O(3)→ ClO + O(2) reaction, which is an elementary step in a catalytic cycle that leads to the destruction of ozone in the stratosphere. The PES is constructed by interpolation of quantum chemistry data using the method developed by Collins and co-workers. Ab initio data points (energy, gradients and Hessian matrix elements) have been calculated at the UQCISD/aug-cc-pVDZ (unrestricted quadratic configuration interaction with single and double excitations) level of theory. The ab initio calculations predict a markedly non-coplanar (dihedral angle of 80°) transition state for the reaction, located very early in the reactant valley and slightly below the energy of the reactants as long as the spin-orbit splitting is neglected. Quasiclassical trajectory (QCT) calculations have been carried out at several collision energies to investigate the reaction dynamics. The QCT excitation function shows no threshold, displays a minimum at a collision energy of 2.5 kcal mol(-1), and then increases monotonically at larger collision energies. This behaviour is consistent with a barrierless reaction dominated by an oxygen-abstraction mechanism. The calculated product vibrational distributions (strongly inverted for ClO) and rate constants are compared with experimental determinations. Differential cross sections (DCS) summed over all final states are found to be in fairly good agreement with those derived from crossed molecular beam experiments.     

Ab initio study of some peroxides. HOOH, CH3OOH and, CH3OOCH3

The geometries of HOOH, CH₃OOH, and CH₃OOCH₃, were optimized with different basis sets (3-21G, 6-31G^*(*) and D95^**) at different levels of theory (HF, MP2, MP4, and CI). HF/3-21G optimizations result in planar trans conformations for all three peroxides. HF/6-31G^** calculations predict skew conformations for HOOH and CH₃OOH, but a planar trans struture for CH₃OOCH₃. For the larger basis set the calculated bond lengths, especially the O-O bonds, are too short. Optimizations for HOOH including electron correlation at the MP2, MP3, MP4, CI, and CCD level improve the agreement for bond lengths and the OOH angle, but result in dihedral angles Which are too large by 3– 8°. In the case of CH₃OOCH₃, similar calculations at the MP2 and CI level predict planar trans structures instead of the experimentally observed skew conformation. On the other hand, MP4 single point calculations at MP2 optimized parameters result in a correct skew structure. For all three peroxides a computationally “economic” method, i.e., single point calculations at MP2 or MP4 level with HF/3-21G optimized parameters, result in close agreement between calculated and experimental structures.