The Recombination of Propargyl Radicals: Solving the Master Equation

Electronic Spectrum of α-Hydrofulvenyl Radical (C6H7), and a Simple and Accurate Recipe for Predicting Adiabatic Ionization Energies of Resonance-Stabilized Hydrocarbon Radicals

Using a combination of resonant two-photon two-color ionization (R2C2PI) and laser-induced fluorescence/dispersed fluorescence spectroscopy, we have examined the A ~  2A″ ← X ~  2A″ transition of the resonance-stabilized α-hydrofulvenyl radical, produced from methylcyclopentadiene dimer in a jet-cooled discharge. Like the related 1,4-pentadienyl and cyclohexadienyl radicals, the α-hydrofulvenyl Ã-state lifetime is orders of magnitude shorter than the predicted f-value implies, indicative of rapid nonradiative decay. The transition is fully allowed by symmetry but considerably weakened by transition moment interference. Intensity borrowing among a' modes brings about static (i.e., Condon) and vibronic (i.e., Herzberg-Teller) moments of similar size, the result being a spectrum substantially less origin-dominated than is usually observed for extensively delocalized radicals. Twenty A ~ -state modes and twelve X ~ -state modes are identified with high confidence and assignments for several others are suggested. In addition, from a series of two-color appearance potential scans with the A ~ -state zero-point level serving as an intermediate, we obtain a field-free adiabatic ionization energy (AIE) of 7.012(1) eV. For a set of 21 resonance-stabilized radicals bearing 5 to 11 carbon atoms, it emerges that the field-free AIE obtained by R2C2PI methods under jet-cooled conditions lies very close to the average of B3LYP/6-311G++(d,p) (with harmonic zero-point energy) and CBS-QB3 0 K calculations, with a mean absolute deviation of only 0.010(7) eV (approximately 1 kJ/mol). On average, this represents a nearly 10-fold improvement in accuracy over CBS-QB3 predictions for the same set of radicals.

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

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

Unconventional gas-phase preparation of the prototype polycyclic aromatic hydrocarbon naphthalene (C10H8) via the reaction of benzyl (C7H7) and propargyl (C3H3) radicals coupled with hydrogen-atom assisted isomerization

Polycyclic aromatic hydrocarbons (PAHs) are ubiquitous in the interstellar medium and in meteorites such as Murchison and Allende and signify the missing link between resonantly stabilized free radicals and carbonaceous nanoparticles (soot particles, interstellar grains). However, the predicted lifetime of interstellar PAHs of some 10⁸ years imply that PAHs should not exist in extraterrestrial environments suggesting that key mechanisms of their formation are elusive. Exploiting a microchemical reactor and coupling these data with computational fluid dynamics (CFD) simulations and kinetic modeling, we reveal through an isomer selective product detection that the reaction of the resonantly stabilized benzyl and the propargyl radicals synthesizes the simplest representative of PAHs - the 10π Hückel aromatic naphthalene (C₁₀H₈) molecule - via the novel Propargyl Addition-BenzAnnulation (PABA) mechanism. The gas-phase preparation of naphthalene affords a versatile concept of the reaction of combustion and astronomically abundant propargyl radicals with aromatic radicals carrying the radical center at the methylene moiety as a previously passed over source of aromatics in high temperature environments thus bringing us closer to an understanding of the aromatic universe we live in.

Reactions of Methyl Radicals with Aniline Acting as Hydrogen Donors and as Methyl Radical Acceptors

The present investigation theoretically reports the comprehensive kinetic mechanism of the reaction between aniline and the methyl radical over a wide range of temperatures (300-2000 K) and pressures (76-76,000 Torr). The potential energy surface of the C₆H₅NH₂ + CH₃ reaction has been established at the CCSD(T)//M06-2X/6-311++G(3df,2p) level of theory. The conventional transition-state theory (TST) was utilized to calculate rate constants for the elementary reaction channels, while the stochastic RRKM-based master equation framework was applied for the T- and P-dependent rate-coefficient calculation of multiwell reaction paths. Hindered internal rotation and Eckart tunneling treatments were included. The H-abstraction from the -NH₂ group of aniline (to form P1 (C₆H₅NH + CH₄)) has been found to compete with the CH₃-addition on the C atom at the ortho site of aniline (to form IS2) with the atmospheric rate expressions (in cm³ molecule^-1 s^-1) as k_a1 = 7.5 × 10^-23 T^3.04 exp[(-40.63 ± 0.29 kJ·mol^-1)/RT] and k_b2 = 2.29 × 10^-3 T^-3.19 exp[(-56.94 ± 1.17 kJ·mol^-1)/RT] for T = 300-2000 K and P = 760 Torr. Even though rate constants of several reaction channels decrease with increasing pressures, the total rate constant k_total = 7.71 × 10^-17 T^1.20 exp[(-40.96 ± 2.18 kJ·mol^-1)/RT] of the title reaction still increases as the pressure increases in the range of 76-76,000 Torr. The calculated enthalpy changes for some species are in good agreement with the available experimental data within their uncertainties (the maximum deviation between theory and experiment is ∼11 kJ·mol^-1). The T1 diagnostic and spin contamination analysis for all species involved have also been observed. This work provides sound quality rate coefficients for the title reaction, which will be valuable for the development of detailed combustion reaction mechanisms for hydrocarbon fuels.

Automated Reaction Kinetics of Gas-Phase Organic Species over Multiwell Potential Energy Surfaces

Automation of rate-coefficient calculations for gas-phase organic species became possible in recent years and has transformed how we explore these complicated systems computationally. Kinetics workflow tools bring rigor and speed and eliminate a large fraction of manual labor and related error sources. In this paper we give an overview of this quickly evolving field and illustrate, through five detailed examples, the capabilities of our own automated tool, KinBot. We bring examples from combustion and atmospheric chemistry of C-, H-, O-, and N-atom-containing species that are relevant to molecular weight growth and autoxidation processes. The examples shed light on the capabilities of automation and also highlight particular challenges associated with the various chemical systems that need to be addressed in future work.

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

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

Excited state photochemically driven surface formation of benzene from acetylene ices on Pluto and in the outer solar system

NASA's New Horizons mission unveiled a diverse landscape of Pluto's surface with massive regions being neutral in color, while others like Cthulhu Macula range from golden-yellow to reddish comprising up to half of Pluto's carbon budget. Here, we demonstrate in laboratory experiments merged with electronic structure calculations that the photolysis of solid acetylene - the most abundant precipitate on Pluto's surface - by low energy ultraviolet photons efficiently synthesizes benzene and polycyclic aromatic hydrocarbons via excited state photochemistry thus providing critical molecular building blocks for the colored surface material. Since low energy photons deliver doses to Pluto's surface exceeding those from cosmic rays by six orders of magnitude, these processes may significantly contribute to the coloration of Pluto's surface and of hydrocarbon-covered surfaces of Solar System bodies such as Triton in general. This discovery critically enhances our perception of the distribution of aromatic molecules and carbon throughout our Solar System.

Mechanism, Thermochemistry, and Kinetics of the Reversible Reactions: C2H3 + H2 ⇌ C2H4 + H ⇌ C2H5

High-level coupled cluster theory, in conjunction with Active Thermochemical Tables (ATcT) and E,J-resolved master equation calculations were used in a study of the title reactions, which play an important role...

Concluding Remarks on Faraday Discussion on Unimolecular Reactions

This Faraday Discussion, marking the centenary of Lindemann’s explanation of the pressure-dependence of unimolecular reactions, presented recent advances in measuring and computing collisional energy transfer efficiencies, microcanonical rate coefficients, pressure-dependent...

Product Detection of the CH(X2Π) Radical Reaction with Cyclopentadiene: A Novel Route to Benzene

The reaction of the methylidyne radical (CH(X²Π)) with cyclopentadiene (c-C₅H₆) is studied in the gas phase at 4 Torr and 373 K using a multiplexed photoionization mass spectrometer. Under multiple collision conditions, the dominant product channel observed is the formation of C₆H₆ + H. Fitting the photoionization spectrum using reference spectra allows for isomeric resolution of C₆H₆ isomers, where benzene is the largest contributor with a relative branching fraction of 90 (±5)%. Several other C₆H₆ isomers are found to have smaller contributions, including fulvene with a branching fraction of 8 (±5)%. Master Equation calculations for four different entrance channels on the C₆H₇ potential energy surface are performed to explore the competition between CH cycloaddition to a C═C bond vs CH insertion into C-H bonds of cyclopentadiene. Previous studies on CH addition to unsaturated hydrocarbons show little evidence for the C-H insertion pathway. The present computed branching fractions support benzene as the sole cyclic product from CH cycloaddition, whereas fulvene is the dominant product from two of the three pathways for CH insertion into the C-H bonds of cyclopentadiene. The combination of experiment with Master Equation calculations implies that insertion must account for ∼10 (±5)% of the overall CH + cyclopentadiene mechanism.

Gas-phase synthesis of benzene via the propargyl radical self-reaction

Polycyclic aromatic hydrocarbons (PAHs) have been invoked in fundamental molecular mass growth processes in our galaxy. We provide compelling evidence of the formation of the very first ringed aromatic and building block of PAHs-benzene-via the self-recombination of two resonantly stabilized propargyl (C₃H₃) radicals in dilute environments using isomer-selective synchrotron-based mass spectrometry coupled to theoretical calculations. Along with benzene, three other structural isomers (1,5-hexadiyne, fulvene, and 2-ethynyl-1,3-butadiene) and o-benzyne are detected, and their branching ratios are quantified experimentally and verified with the aid of computational fluid dynamics and kinetic simulations. These results uncover molecular growth pathways not only in interstellar, circumstellar, and solar systems environments but also in combustion systems, which help us gain a better understanding of the hydrocarbon chemistry of our universe.

Discovery of the propargyl radical (CH2CCH) in TMC-1: one of the most abundant radicals ever found and a key species for cyclization to benzene in cold dark clouds

We present the first identification in interstellar space of the propargyl radical (CH2CCH). This species was observed in the cold dark cloud TMC-1 using the Yebes 40m telescope. The six strongest hyperfine components of the 20,2-10,1 rotational transition, lying at 37.46 GHz, were detected with signal-to-noise ratios in the range 4.6-12.3 σ. We derive a column density of 8.7 × 1013 cm-2 for CH2CCH, which translates to a fractional abundance relative to H2 of 8.7 × 10-9. This radical has a similar abundance to methyl acetylene, with an abundance ratio CH2CCH/CH3CCH close to one. The propargyl radical is thus one of the most abundant radicals detected in TMC-1, and it is probably the most abundant organic radical with a certain chemical complexity ever found in a cold dark cloud. We constructed a gas-phase chemical model and find calculated abundances that agree with, or fall two orders of magnitude below, the observed value depending on the poorly constrained low-temperature reactivity of CH2CCH with neutral atoms. According to the chemical model, the propargyl radical is essentially formed by the C + C2H4 reaction and by the dissociative recombination of C3Hn + ions with n = 4-6. The propargyl radical is believed to control the synthesis of the first aromatic ring in combustion processes, and it probably plays a key role in the synthesis of large organic molecules and cyclization processes to benzene in cold dark clouds.

Direct Measurement of Radical-Catalyzed C6H6 Formation from Acetylene and Validation of Theoretical Rate Coefficients for C2H3 + C2H2 and C4H5 + C2H2 Reactions

The addition of vinylic radicals to acetylene is an important step contributing to the formation of polycyclic aromatic hydrocarbons in combustion. The overall reaction 3C₂H₂ → C₆H₆ could result in large benzene yields, but without accurate rate parameters validated by experiment, the extent of aromatic ring formation from this pathway is uncertain. The addition of vinyl radicals to acetylene was investigated using time-resolved photoionization time-of-flight mass spectrometry at 500 and 700 K and 5–50 Torr. The formation of C₆H₆ was observed at all conditions, attributed to sequential addition to acetylene followed by cyclization. Vinylacetylene (C₄H₄) was observed with increasing yield from 500 to 700 K, attributed to the β-scission of the thermalized 1,3-butadien-1-yl radical and the chemically activated reaction C₂H₃ + C₂H₂ → C₄H₄ + H. The measured kinetics and product distributions are consistent with a kinetic model constructed using pressure- and temperature-dependent reaction rate coefficients computed from previously reported ab initio calculations. The experiments provide direct measurements of the hypothesized C₄H₅ intermediates and validate predictions of pressure-dependent addition reactions of vinylic radicals to C₂H₂, which are thought to play a key role in soot formation.

Theoretical Study of the Kinetics of the Gas-Phase Reaction between Phenyl and Amino Radicals

The potential energy surface (PES) of the C₆H₅ + NH₂ reaction has been investigated by using ab initio CCSD(T)//B3LYP/6-311++G(3df,2p) calculations. The conventional transition-state theory (TST) and the variable reaction coordinate-TST (VRC-TST) have been used to predict the rate constants for the channels possessing tight and barrierless transition states, respectively. The Rice-Ramsperger-Kassel-Marcus/Master equation (RRKM/ME) theory has been utilized to determine the pressure-dependent rate constants for these reactions. The PES shows that the reaction begins with an exothermic barrierless addition of NH₂ to C₆H₅ producing the vital intermediate state, namely, aniline (C₆H₅NH₂, IS1). Once IS1 is generated, it can further isomerize to various intermediate states, which can give rise to different products, including PR4 (4,5,6-trihydro-1-amino phenyl + H₂), PR5 (3,4,5,6-tetrahydro phenyl + NH₃), PR6 (2,3,5,6-tetrahydro-1-imidogen phenyl + H₂), PR9 (3,4,5,6-tetrahydro-1-imidogen phenyl + H₂), and PR10 (2,5,6-trihydro-1-amino phenyl + H₂), of which the most stable product, PR5, was formed by the most favorable channel going through the two advantageous transition states T1/11 (-28.9 kcal/mol) and T11P5 (-21.5 kcal/mol). The calculated rate constants for the low-energy channel, 1.37 × 10^-9 and 2.16 × 10^-11 cm³ molecule^-1 s^-1 at T = 300, P = 1 Torr and T = 2000 K, P = 760 Torr, respectively, show that the title reaction is almost pressure- and temperature-dependent. The negative temperature-dependent rate coefficients can be expressed in the modified Arrhenius form of k ₁ = 8.54 × 10¹³ T ^-7.20 exp (-7.07 kcal·mol^-1/RT) and k ₂ = 2.42 × 10¹⁵ T ^-7.61 exp (-7.75 kcal·mol^-1/RT) at 1 and 10 Torr, respectively, and in the temperature range of 300-2000 K. The forward and reverse rate coefficients as well as the high-pressure equilibrium constants of the C₆H₅ + NH₂ ↔ IS1 process were also predicted; their values revealed that its kinetics do not depend on pressure at low temperature but strongly depend on pressure at high temperature. Moreover, the predicted formation enthalpies of reactants and the enthalpy changes of some channels are in good agreement with the experimental results.

Rate Coefficients of the Reaction of OH with Allene and Propyne at High Temperatures

Allene (H₂C═C═CH₂; a-C₃H₄) and propyne (CH₃C≡CH; p-C₃H₄) are important species in various chemical environments. In combustion processes, the reactions of hydroxyl radicals with a-C₃H₄ and p-C₃H₄ are critical in the overall fuel oxidation system. In this work, rate coefficients of OH radicals with allene (OH + H₂C═C═CH₂ → products) and propyne (OH + CH₃C≡CH → products) were measured behind reflected shock waves over the temperature range of 843-1352 K and pressures near 1.5 atm. Hydroxyl radicals were generated by rapid thermal decomposition of tert-butyl hydroperoxide ((CH₃)₃-CO-OH), and monitored by narrow line width laser absorption of the well-characterized R₁(5) electronic transition of the OH A-X (0,0) electronic system near 306.7 nm. Results show that allene reacts faster with OH radicals than propyne over the temperature range of this study. Measured rate coefficients can be expressed in Arrhenius form as follows: k_allene+OH(T) = 8.51(±0.03) × 10^-22T^3.05 exp(2215(±3)/T), T = 843-1352 K; k_propyne+OH(T) = 1.30(±0.07) × 10^-21T^3.01 exp(1140(±6)/T), T = 846-1335 K.

Are the three hydroxyphenyl radical isomers created equal?--The role of the phenoxy radical

We have investigated the thermal decomposition of the three hydroxyphenyl radicals (˙C6H4OH) in a heated microtubular reactor. Intermediates and products were identified isomer-selectively applying photoion mass-selected threshold photoelectron spectroscopy with vacuum ultraviolet synchrotron radiation. Similarly to the phenoxy radical (C6H5-O˙), hydroxyphenyl decomposition yields cyclopentadienyl (c-C5H5) radicals in a decarbonylation reaction at elevated temperatures. This finding suggests that all hydroxyphenyl isomers first rearrange to form phenoxy species, which subsequently decarbonylate, a mechanism which we also investigate computationally. Meta- and para-radicals were selectively produced and spectroscopically detectable, whereas the ortho isomer could not be traced due to its fast rethermalization and rapid decomposition in the reactor. A smaller barrier to isomerization to phenoxy was found to be the reason for this observation. Since hydroxyphenyl species may be present under typical sooting conditions in flames, the resonantly stabilized cyclopentadienyl radical adds to the hydrocarbon pool and can contribute to the formation of polycyclic aromatic hydrocarbons, which are precursors in soot formation.

Ab initio pressure-dependent reaction kinetics of methyl propanoate radicals

The unimolecular dissociation and isomerization kinetics of the three methyl propanoate (MP) radicals, CH₃CH₂C(O)OĊH₂ (MP-m), CH₃ĊHC(O)OCH₃ (MP-α), and ĊH₂CH₂C(O)OCH₃ (MP-β), are theoretically investigated using high-level ab initio methods and the Rice–Ramsperger–Kassel–Marcus (RRKM)/master equation (ME) theory.

Experimental and theoretical studies into the formation of C4-C6 products in partially chlorinated hydrocarbon pyrolysis systems: a probabilistic approach to congener-specific yield predictions

This work presents a study of the pyrolytic formation of vinylacetylene and benzene congeners formed from chlorinated hydrocarbon precursors, a complex, multipath polymerization system formed in a monomer-rich environment. (Co-)pyrolyses of dichloro- and trichloroethylene yield a rich array of products, and assuming a single dominant underlying growth mechanism, this (on comparing expected and observed products) allows a number of potentially competing channels to C4 and C6 products to be ruled out. Poor congener/isomer descriptions rule out even-carbon radical routes, and the absence of C3 and C5 products rule out odd-carbon processes. Vinylidenes appear unable to describe the increased reactivity of acetylenes with chlorination noted in our experiments, leaving molecular acetylene dimerization processes and, in C6 systems, the closely related Diels-Alder cyclization as the likely reaction mechanism. The feasibility of these routes is further supported by ab initio calculations. However, some of the most persuasive evidence is provided by congener-specific yield predictions enabled by the construction of a probability tree analogue of kinetic modeling. This approach is relatively quick to construct, provides surprisingly accurate predictions, and may be a very useful tool in screening for important reaction channels in poorly understood congener- or isomer-rich reaction systems.

C8H6 Thermal Chemistry. 7-Methylenecyclohepta-1,3,5-dienyne (Heptafulvyne) by Flash Vacuum Thermolysis–Matrix Isolation. Chemical Activation in the Rearrangements of Phenylenedicarbenes and of Benzocyclobutadiene to Phenylacetylene

Methylenecycloheptadienyne 11 (heptafulvyne) is obtained very cleanly by flash vacuum thermolysis (FVT) of the diazobenzocyclobutene precursor 8 at 400°C followed by isolation as a neat solid at 77 K or in an Ar matrix at 7–10 K. Compound 11 is a yellow solid, stable till ~–100°C in the neat state. The diazo compound itself (2) is observable by IR spectroscopy following mild decomposition of the tosylhydrazone salt 1 at 115°C. FVT of 8 at 200°C also generates diazo compound 2 as observed by IR spectroscopy and on-line mass spectrometry. FVT of 8 at 600–800°C causes rearrangement of 11 to phenylacetylene 12 and benzocyclobutadiene 13. Mechanisms for the rearrangements are proposed. Facile rearrangement of benzocyclobutadiene to phenylacetylene is ascribed to chemical activation, which is also seen to be involved in the rearrangement of p-, m-, and o-phenylenebiscarbenes 25–27 to phenylacetylene 12.

Molecular mechanisms in the pyrolysis of unsaturated chlorinated hydrocarbons: formation of benzene rings. 2. Experimental and kinetic modeling studies

The mechanism of formation of benzene rings during the pyrolysis of dichloro- and trichloroethylenes has been investigated by the method of laser powered homogeneous pyrolysis coupled with product analysis by gas chromatography. Additionally, selected (co)pyrolyses between the chlorinated ethylenes, CH2Cl2, C4Cl4, C4Cl6, and C2H2 have been performed to explicitly probe the roles of 2C3 and C4/C2 reaction pairs in aromatic growth. The presence of odd-carbon products in neat C4Cl6 pyrolyses indicates that 2C3 processes are operative in these systems; however, comparison with product yields from C2HCl3 suggests that C4/C2 processes dominate most other systems. This is further evidenced by an absence of C3 and other odd-carbon species in (co)pyrolyses with dichloromethane which should seed C3-based growth. The reactions of perchlorinated C4 species C4Cl5, C4Cl3, and C4Cl4 with C2Cl2 were subsequently explored through extensive kinetic simulations of the possible reaction pathways based on previous kinetic models and the exhaustive quantum chemical investigations of our preceding work. The experimental and theoretical results strongly suggest that, at moderate temperatures, aromatic ring formation from chlorinated ethylenes normally follows a Diels-Alder coupling of C4 and C2 molecular units followed by internal shifts; the one exception is the C4Cl4 + C2Cl2 system, where steric factors lead to the formation of nonaromatic products. There is little evidence for radical-based routes in these systems.

Ab initio kinetics for the decomposition of hydroxybutyl and butoxy radicals of n-butanol

The decomposition kinetics of the hydroxybutyl and butoxy radicals (C4H9O) arising via H abstraction from n-butanol were studied theoretically with ab initio transition-state-theory-based master equation analyses. Stationary points on the C4H9O potential energy surface were calculated at either the RQCISD(T)/CBS//B3LYP/6-311++G(d,p) level or the RQCISD(T)/CBS//CASPT2/aug-cc-pVDZ level. Unimolecular pressure- and temperature-dependent rate coefficients were calculated over broad ranges of temperature (300-2500 K) and pressure (1.3 × 10(-3) to 10(2) atm) by solving the time-dependent multiple-well master equation. The "well merging" phenomenon was observed and analyzed for its influence on the branching ratios and rate coefficients. The theoretical predictions were compared with the available experimental and theoretical data and any discrepancies were analyzed. The predicted rate coefficients are represented with forms that may readily be used in combustion modeling of n-butanol.

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.

Selecting Methods to Solve Multi-Well Master Equations

There are several competing methods commonly used to solve energy grained master equations describing gas-phase reactive systems. When it comes to selecting an appropriate method for any particular problem, there is little guidance in the literature. In this paper we directly compare several variants of spectral and numerical integration methods from the point of view of computer time required to calculate the solution and the range of temperature and pressure conditions under which the methods are successful. The test case used in the comparison is an important reaction in combustion chemistry and incorporates reversible and irreversible bimolecular reaction steps as well as isomerizations between multiple unimolecular species. While the numerical integration of the ODE with a stiff ODE integrator is not the fastest method overall, it is the fastest method applicable to all conditions.