Pressure and temperature dependence of the reaction of vinyl radical with alkenes II: Measured rates and predicted product distributions for vinyl + propene

This work reports measurements of the absolute rate coefficients and Rice-Ramsperger-Kassel-Markus (RRKM) master equation (ME) simulations of the C₂H₃ + C₃H₆ reaction. Direct kinetic studies were performed over a temperature range of 300-700 K and pressures of 15, 25, and 100 Torr. Vinyl radicals were generated by laser photolysis of vinyl iodide at 266 nm, and time-resolved absorption spectroscopy was used to probe vinyl radicals through absorption at 423.2 nm. A weighted modified Arrhenius fit to the experimental rate constant is k₁ = (1.3 ± 0.2) × 10⁻¹² cm³ molecule⁻¹ s⁻¹(T/1000)^1.6 exp[−(1510 ± 80/T)]. Fifteen stationary points and 48 transition states on the C₅H₉ potential energy surface (PES) were calculated using the G3 method in Gaussian 03. RRKM/ME simulations were performed using VariFlex on a simplified PES to predict pressure dependent rate coefficients and branching fractions for the major channels. For temperatures between 350 and 700 K, the calculated rate coefficient agrees with the experimental rate coefficient within 20%. At low temperatures, the primary products are the initial adducts 4-penten-2-yl and 2-methyl-3-buten-1-yl. At higher temperatures, the dominant products are 1,3-butadiene + methyl, allyl + ethene, and 1,3-pentadiene + H. Although C₂H₃ + C₃H₆ → allyl + ethene is thermodynamically favored, the simulations predict that it does not become the dominant product until 1700 K.     

Kinetics and thermochemistry of the addition of atomic chlorine to acetylene

Atomic Cl was generated by pulsed laser photolysis at 193 nm of CCl₄, and was monitored by time resolved resonance fluorescence in the course of reaction with excess C₂H₂, diluted in Ar bath gas at pressures from 13 to 800 mbar. At 288 K simple pseudo first order kinetics were observed. Over 365-430 K bi-exponential decays were obtained, because of equilibration between the β-chlorovinyl adduct and the reactants. The ratios of forward and reverse rate constants yield Δ_fH₂₉₈(CHCHCl) = 274.0 ± 1.0 kJ mol⁻¹ via a Third-Law analysis of the carbon-chlorine bond strength. The thermochemistry compares well with that predicted by an initio theory. The effective second-order rate constant was pressure dependent and was analyzed using Troe’s unimolecular formalism. Over the whole temperature range the low-pressure limiting value for addition, with Ar bath gas, is given by k₀ = 4.1 × 10⁻³⁰ (T/300 K)^−2.47 cm⁶ molecule⁻² s⁻¹.     

Metal ion-molecule kinetics at combustion temperatures: The reaction of Ca with O2

The high-temperature photochemistry (HTP) technique, previously used for reactions of neutral species, has been adapted to the study of atomic metal ion-molecule reactions. Ca⁺ ions were generated by 193 nm multi-photon photolysis of calcium acetyl acetonate and its pyrolysis fragments. The relative ion concentrations were monitored by laser-induced fluorescence at 393.4 nm. Ar was used as the bath gas. The data for the Ca⁺ + O₂ + M → CaO₂⁺ + M association reaction (1) are fitted by k₁(907-1425 K) = 3.5 × 10⁻³² exp(+3161 K/T) cm⁶ molecule⁻² s⁻¹. Combining with an approximate k₁(296 K) value in the literature leads to k₁(296-1425 K) = 5.8 × 10⁻²² (T/K)^−2.9 exp(−601 K/T) cm⁶ molecule⁻² s⁻¹. Over much of the observed temperature range reaction (1) has much smaller rate coefficients than the corresponding neutral Ca association reaction. Reaction (1) is shown to behave very similarly to the O₂ association reaction with neutral K atoms, with which Ca⁺ is iso-electronic. This suggests that the initial step is ion-pair complex formation of the superoxide Ca²⁺(O₂)⁻, which is also consistent with results from density functional calculations. The k₁ values are rationalized via Troe’s unimolecular formalism, which leads to good accord with the experiments.     

Ab initio chemical kinetics for the reactions of HNCN with O(P) and O2

The kinetics and mechanisms of the reactions of cyanomidyl radical (HNCN) with oxygen atoms and molecules have been investigated by ab initio calculations with rate constant prediction. The doublet and quartet state potential energy surfaces (PESs) of the two reactions have been calculated by single-point calculations at the CCSD(T)/6-311+G(3df, 2p) level based on geometries optimized at the CCSD/6-311++G(d, p) level. The rate constants for various product channels of the two reactions in the temperature range of 300-3000 K are predicted by variational transition state and RRKM theories. The predicted total rate constants of the O(³P) + HNCN reaction at 760 Torr Ar pressure can be represented by the expressions k_total (O + HNCN) = 3.12 × 10⁻¹⁰ × T^−0.05 exp (−37/T) cm³ molecule⁻¹ s⁻¹ at T = 300-3000 K. The branching ratios of primary channels of the O(³P) + HNCN are predicted: k₁ for producing the NO + CNH accounts for 0.72-0.64, k₂ + k₉ for producing the ³NH + NCO accounts for 0.27-0.32, and k₆ for producing the CN + HNO accounts for 0.01-0.07 in the temperature range studied. Meanwhile, the predicted total rate constants of the O₂ + HNCN reaction at 760 Torr Ar pressure can be represented by the expression, k_total(O₂ + HNCN) = 2.10 × 10⁻¹⁶ × T^1.28exp (−12200/T) cm³ molecule⁻¹ s⁻¹ at T = 300-3000 K. The predicted branching ratio for k₁₁ + k₁₃ producing HO₂ + ³NCN as the primary products accounts for 0.98-1.00 in the temperature range studied.     

Kinetics for the reaction of phenyl radical with phenylacetylene and styrene

The kinetics for the reactions of C₆H₅ with phenylacetylene and styrene have been measured by CRDS in the temperature range 297-409 K under an Ar pressure of 3.6 Torr. The total rate constants can be given by the following Arrhenius expressions (in units of cm³ mol⁻¹ s⁻¹): k₁(C₆H₅ + C₆H₅C₂H) = 10^13.0±0.1exp [−(2430 ± 150)/RT] and k₂(C₆H₅ + C₆H₅C₂H₃) = 10^13.3±0.1 exp [−(2570 ± 180)/T]. Additional DFT and MP2 calculations have been carried out to assist our interpretation of the measured kinetic data. The addition of C₆H₅ to the terminal CH_x (x = 1 or 2) sites is predicted to be the dominant channel for both reactions. The calculated bimolecular rate constants are in reasonable agreement with experimental values for the temperature range studied.     

Theoretical study on the kinetics for OH reactions with CH3OH and C2H5OH

Kinetics and mechanisms for reactions of OH with methanol and ethanol have been investigated at the CCSD(T)/6-311 + G(3df, 2p)//MP2/6-311 + G(3df, 2p) level of theory. The total and individual rate constants, and product branching ratios for the reactions have been computed in the temperature range 200-3000 K with variational transition state theory by including the effects of multiple reflections above the wells of their pre-reaction complexes, quantum-mechanical tunneling and hindered internal rotations. The predicted results can be represented by the expressions k₁ = 4.65 × 10⁻²⁰ × T^2.68 exp(414/T) and k₂ = 9.11 × 10⁻²⁰ × T^2.58 exp(748/T) cm³ molecule⁻¹ s⁻¹ for the CH₃OH and C₂H₅OH reactions, respectively. These results are in reasonable agreements with available experimental data except that of OH + C₂H₅OH in the high temperature range. The former reaction produces 96-89% of the H₂O + CH₂OH products, whereas the latter process produces 98-70% of H₂O + CH₃CHOH and 2-21% of the H₂O + CH₂CH₂OH products in the temperature range computed (200-3000 K).     

Kinetics of the H + NCO reaction

Ab initio transition state theory (TST) based master equation simulations are used to predict the temperature and pressure dependence of the H + NCO reaction rate and product branching. The barrierless entrance channels to form singlet HNCO and NCOH are studied with variable reaction coordinate TST employing a potential energy surface obtained from multi-reference configuration interaction ab initio calculations. The remaining channels, including reactions on the triplet surface, are studied with standard TST methods employing high level electronic structure results. The energy transfer parameters for the master equation simulations arise from a fit to the experimentally observed HNCO dissociation rate. The lowest energy threshold to formation of bimolecular products, ³NH + CO, lies well below the reactants. The bottleneck for intersystem crossing, which precedes the formation of ³NH + CO from the singlet adducts, becomes the dominant bottleneck for that channel at quite low energies relative to reactants. The effect of this bottleneck is studied with model calculations designed to reproduce detailed experimental observations of photolysis branching ratios. This bottleneck greatly reduces the flux from H + NCO to ³NH + CO via the singlet adducts. As a result, stabilization and reaction on solely the triplet surface are significant components of the overall rate. The present predictions for the high pressure and collisionless limit rate coefficients are accurately reproduced over the 200-2500 K range by the expressions, 1.53 × 10⁻⁵T^−1.86exp(−399/T) + 1.07 × 10³T^−3.15exp(−15219/T) and 5.62 × 10⁻¹²T^0.493exp(148/T) cm³ molecule⁻¹ s⁻¹, respectively, where T is in K. These predictions are in reasonably satisfactory agreement with the somewhat discordant experimental rate measurements.     

Temperature and pressure effects on formation and decomposition of phenylvinylperoxy radicals in the C6H5C2H2 + O2 reaction

The effects of temperature and pressure on the formation and decomposition of C₆H₅C₂H₂O₂ in the C₆H₅C₂H₂ + O₂ reaction have been investigated at temperatures from 298 to 378 K by directly monitoring the C₆H₅C₂H₂O₂ radical in the visible region by cavity ringdown spectrometry (CRDS). The rate constant for the C₆H₅C₂H₂ + O₂ association and that for fragmentation of C₆H₅C₂H₂O₂ were found to be k₁ (C₆H₅C₂H₂ + O₂ → C₆H₅C₂H₂O₂) = (3.20 ± 1.19) × 10¹¹ exp(+760/T) cm³ mol⁻¹ s⁻¹ and k₂ (C₆H₅C₂H₂ O₂ → C₆H₅CHO + HCO) = (1.68 ± 0.13) × 10⁴ s⁻¹, respectively. Additional kinetic measurements by pulsed laser photolysis/mass spectrometry show that C₆H₅CHO was produced in the C₆H₅C₂H₂ + O₂ reaction as predicted and the formation of C₆H₅CHO from the decomposition of C₆H₅C₂H₂O₂ is temperature-independent, consistent with the CRDS experimental data.     

Specific mass shift of potassium 5P1/2 state

We have measured the isotope shift between ⁴¹K and ³⁹K in the 4s_1/2 → 5p_1/2 transition at 405 nm using saturation spectroscopy. Our measured isotope shift is 456.1 ± 0.8 MHz, implying a residual isotope shift (sum of specific mass shift and field shift) of −52.7 ± 0.8 MHz. We deduce a specific mass shift of −40 ± 5 MHz, which would imply that the 5p_1/2 state has a considerably larger specific mass shift than the 4p_1/2 state. We have in addition measured the 5p_1/2 hyperfine splitting for ⁴¹K.     

Synthesis and magnetic properties of the size-controlled Mn-Zn ferrite nanoparticles by oxidation method

Size-controlled Mn_0.67Zn_0.33Fe₂O₄ nanoparticles in the wide range from 80 to 20 nm have been synthesized, for the first time, using the oxidation method. It has been demonstrated that the particle size can be tailor-made by varying the concentration of the oxidant. The magnetization of the 80 nm particles was 49 A m² kg⁻¹ compared to 34 A m² kg⁻¹ for the 20 nm particles. The Curie temperatures for all the samples are found to be within 630±5 K suggesting that there is no size-dependent cation distribution. The critical particle size for the superparamagnetic limit is found to be about 25 nm. The effective magnetic anisotropy constant is experimentally determined to be 7.78 kJ m⁻³ for the 25 nm particles, which is about an order of magnitude higher than that of the bulk ferrite.     

Absorption cross-sections off HFC-23 at atmospheric conditions

Spectral absorption cross-sections, k_ν (cm⁻¹ atm⁻¹), were measured in the 8.64 μm band system of HFC-23 using a high-resolution Fourier transform spectrometer. Temperature and total (N₂-broadening) pressure were selected to represent atmospheric layers according to various atmospheric models. The measured absolute intensity (or integrated cross-section) was obtained to be 0.693±0.029×10⁻¹⁶ cm molecule⁻¹.     

Thermal isomerization and cyclization of 1-naphthylacetylene: Experimental results, quantum chemical and transition state theory calculations

The isomerization of 1-naphthylacetylene diluted in argon was studied behind reflected shock waves in a 2 in i.d. single pulse shock tube over the temperature range 1000-1250 K and overall densities of ∼3 × 10⁻⁵ mol/cm³. The only reaction product found in the post shock mixtures was acenaphthylene. The first order rate constant of the isomerization was found to be k = 2.08 × 10¹² exp(−54.2 × 10³/RT) s⁻¹, where R is expressed in units of cal/K mol. Potential energy surfaces of the cyclization reaction 1-naphthylacetylene → acenaphthylene and the isomerization 1-naphthylacetylene → 2-naphthylacetylene were calculated using the Becke three-parameter hybrid method with Lee-Yang-Parr correlation functional approximation (B3LYP). Structure, energy and frequency calculations were carried out with the Dunning correlation consistent polarized double ζ (cc-pVDZ) basis set. The rate constant (k_∞) for the 1-naphthylacetylene → acenaphthylene cyclization was calculated using transition state theory, the value obtained is k_∞ = 3.52 × 10¹² exp(−55.9 × 10³/RT) s⁻¹, where R is expressed in units of cal/K mol. The agreement between the experiment and the calculations is very good. RRKM calculations were done to transfer k_∞ to the pressure of the single pulse shock tube experiments. In view of high temperature and the large molecule involved the deviation from k_∞ is very small. The isomerization 1-naphthylacetylene → 2-naphthylacetylene proceeds via the formation of an unstable intermediate 1,2-naphthalenocyclobutene and has a high barrier of ∼83.5 kcal/mol. In view of this high barrier, the isomerization cannot compete with the cyclization that proceeds with a barrier of ∼56 kcal/mol.     

Thermal decomposition of toluene: Overall rate and branching ratio

The two-channel thermal decomposition of toluene, C₆H₅CH₃ → C₆H₅CH₂ + H (1) and C₆H₅CH₃ → C₆H₅ + CH₃ (2), was investigated in shock tube experiments over the temperature range of 1400-1780 K at a pressure of 1.5 (±0.1) bar. Rate coefficients for reactions (1) and (2) were determined by monitoring benzyl radical (C₆H₅CH₂) absorption at 266 nm during the decomposition of toluene diluted in argon and modeling the temporal behavior of the benzyl concentration with a kinetic model. The first-order rate coefficients determined at a pressure of 1.5 bar are expressed by k₁(T) = 2.09 × 10¹⁵ exp (−87510 [cal/mol]/RT) [s⁻¹] and k₂(T) = 2.66 × 10¹⁶ exp (−97880 [cal/mol]/RT) [s⁻¹]. The resulting branching ratio, k₁/(k₁ + k₂), ranges from 0.8 at 1350 K to 0.6 at 1800 K.     

The effect of duration of diffusion on Ag diffusion coefficients in YBa2Cu3O7

The Ag diffusion in superconducting YBa₂Cu₃O₇ (YBaCuO) ceramic has been studied over the duration of the diffusion range 5-24 h in the temperature range 700-850 °C by the energy-dispersive X-ray fluorescence (EDXRF) technique. For the excitation of silver atoms, an annular Am-241 radioisotope source (50 mCi) emitting 59.543 keV photons was used. The temperature dependences of silver diffusion coefficients in grains (D₁) and over the grain boundaries in the range 700-850 °C (D₂) are described by the relations D₁=1.4×10⁻² exp[−(1.18±0.10)/kT] and D₂=3.1×10⁻⁴ exp[−(0.87±0.10)/kT].     

The conductance and capacitance-frequency characteristics of Au/pyronine-B/p-type Si/Al contacts

The rectifying junction characteristics of the organic compound pyronine-B (PYR-B) film on a p-type Si substrate have been studied. The PYR-B has been evaporated onto the top of p-Si surface. The barrier height and ideality factor values of 0.67 ± 0.02 eV and 2.02 ± 0.03 for this structure have been obtained from the forward bias current-voltage (I-V) characteristics. The energy distribution of the interface states and their relaxation time have been determined from the forward bias capacitance-frequency and conductance-frequency characteristics in the energy range of ((0.42 ± 0.02) − E_v)-((0.66 ± 0.02) − E_v) eV. The interface state density values ranges from (4.21 ± 0.14) × 10¹³ to (3.82 ± 0.24) × 10¹³ cm⁻² eV⁻¹. Furthermore, the relaxation time ranges from (1.65 ± 0.23) × 10⁻⁵ to (8.12 ± 0.21) × 10⁻⁴ s and shows an exponential rise with bias from the top of the valance band towards the midgap.     

High-pressure laminar flame speeds and kinetic modeling of carbon monoxide/hydrogen combustion

Laminar flame speeds were accurately measured for CO/H₂/air and CO/H₂/O₂/helium mixtures at different equivalence ratios and mixing ratios by the constant-pressure spherical flame technique for pressures up to 40 atmospheres. A kinetic mechanism based on recently published reaction rate constants is presented to model these measured laminar flame speeds as well as a limited set of other experimental data. The reaction rate constant of CO + HO₂ → CO₂ + OH was determined to be k = 1.15 × 10⁵T^2.278 exp(−17.55 kcal/RT) cm³ mol⁻¹ s⁻¹ at 300-2500 K by ab initio calculations. The kinetic model accurately predicts our measured flame speeds and the non-premixed counterflow ignition temperatures determined in our previous study, as well as homogeneous system data from literature, such as concentration profiles from flow reactor and ignition delay time from shock tube experiments.     

Vibronic and cation spectroscopy of m-chloroaniline

We applied the resonant two-photon ionization and mass-analyzed threshold ionization spectroscopic techniques to record the vibronic and cation spectra of m-chloroaniline. The band origin of the first electronic transition was found to be 33 658 ± 2 cm⁻¹, whereas the adiabatic ionization energy was determined to be 63 958 ± 5 cm⁻¹. Within our experimental detection limit, these measured values are the same for both of the ³⁵Cl and ³⁷Cl isotopomers. The observed active modes of this molecule in the electronically excited S₁ and cationic ground D₀ states mainly involve the in-plane ring deformation and substituent-sensitive bending vibrations.     

Quantitative measurement of integrated band intensities of benzene vapor in the mid-infrared at 278, 298, and 323 K

Pressure broadened (1 atm. N₂) laboratory spectra of benzene vapor (in natural abundance) were recorded at 278, 298, and 323 K, covering 600-6500 cm⁻¹. The spectra were recorded at a resolution of 0.112 cm⁻¹ using a commercial Fourier transform spectrometer. The pressure of each benzene vapor sample was measured using high-precision capacitance manometers, and a minimum of nine sample pressures were recorded for each temperature. The samples were introduced into a temperature-stabilized static cell (19.94(1) cm pathlength) that was hard-mounted into the spectrometer. From these data a fit composite spectrum was calculated for each temperature. The number density for the three composite spectra was normalized to 296 K. The spectra give the absorption coefficient (cm² molecule⁻¹, naperian units) as a function of wavenumber. From these spectra integrated band intensities (cm molecule⁻¹ and atm⁻¹ cm⁻²) for intervals corresponding to the stronger benzene bands were calculated and were compared with previously reported values. We discuss and quantify error sources and estimate our systematic (NIST Type-B) errors to be 3% for the stronger bands. The measured absorption coefficients and integrated band intensities are useful for remote sensing applications such as measurements of planetary atmospheres and assessment of the environmental impact of terrestrial oil fire emissions.     

A shock-tube and theory study of the dissociation of acetone and subsequent recombination of methyl radicals

The dissociation of acetone: CH₃COCH₃ → CH₃CO + CH₃, quickly followed by CH₃CO → CH₃ + CO, has been examined with Laser-Schlieren measurements in incident shock waves over 32-717 Torr and 1429-1936 K using 5% acetone dilute in krypton. A few very low pressure experiments (∼10 Torr) were used in a marginal effort to resolve the extremely fast vibrational relaxation of this molecule. This effort was partly motivated as a test for molecular, “roaming methyl” reactions, and also as a source of methyl radicals to test the application of a recent high-temperature mechanism for ethane decomposition [J.H. Kiefer, S. Santhanam, N.K. Srinivasan, R.S. Tranter, S.J. Klippenstein, M.A. Oehlschlaeger, Proc. Combust. Inst. 30 (2005) 1129-1135] on the reverse methyl combination. The gradient profiles show strong initial positive gradients and following negative values fully consistent with methyl radical formation and its following recombination. Thus C-C fission is certainly a large part of the process and molecular channels cannot be responsible for more than 30% of the dissociation. Rates obtained for the C-C fission show strong falloff well fit by variable reaction coordinate transition state theory when combined with a master equation. The calculated barrier is 82.8 kcal/mol, the fitted 〈ΔE〉_down = 400 (T/298) cm⁻¹, similar to what was found in a recent study of C-C fission in acetaldehyde, and the extrapolated k_∞ = 10^25.86 T^−2.72 exp(−87.7 (kcal/mol)/RT), which agrees with the literature rate for CH₃ + CH₃CO. Large negative (exothermic) gradients appearing late from methyl combination are accurately fit in both time of onset and magnitude by the earlier ethane dissociation mechanism. The measured dissociation rates are in close accord with one earlier shock-tube study [K. Sato, Y. Hidaka, Combust. Flame 122 (2000) 291-311], but show much less falloff than the high pressure experiments of Ernst et al. [J. Ernst, K. Spindler, H.Gg. Wagner, Ber. Bunsenges. Phys. Chem. 80 (1976) 645-650].     

Multi-step laser excitation of the highly excited states of zinc

We report experimental data on the highly excited states of zinc in the energy range 74,625-75,740 cm⁻¹ using two-step laser excitation scheme in conjunction with a thermionic diode ion detector. The 4s4p ³P₁ inter-combination level at 32501.399 cm⁻¹ was populated using a frequency doubled dye laser. The 4s5s ³S₁ level at 53672.28 cm⁻¹ gets populated from the ASE (amplified spontaneous emission) of the second step dye laser. The Rydberg series 4snp ³P₂ (12 ? n ? 60), 4snp ¹P₁ (16 ? n ? 30) and parity forbidden transitions 4sns ³S₁ (19 ? n ? 44) have been observed. A two parameter fit to excitation energies of the observed series yields the binding energy of the 4s5s ³S₁ level as 22097.03 ± 0.03 cm⁻¹ and consequently, the first ionization potential of zinc is determined as 75769.31 ± 0.05 cm⁻¹, that is in excellent agreement with the earlier work.