Classical trajectory and statistical adiabatic channel study of the dynamics of capture and unimolecular bond fission. IV. Valence interactions between atoms and linear rotors

Influence of Molecular Parameters on Rate Constants of Thermal Dissociation/Recombination Reactions: The Reaction System CF4 ⇄ CF3 + F

The possibilities to extract incompletely characterized molecular parameters from experimental thermal rate constants for dissociation and recombination reactions are explored. The reaction system CF₄ (+M) ⇄ CF₃ + F (+M) is chosen as a representative example. A set of falloff curves is constructed and compared with the available experimental database. Agreement is achieved by minor (unfortunately not separable) adjustments of reaction enthalpy and collisional energy transfer parameters.

Statistical theory for the reaction N + OH → NO + H: thermal low-temperature rate constants

The reaction N + OH → NO + H involves the intermediate formation of NOH adducts which in part rearrange to HNO conformers. A statistical treatment of the process is developed in which an initial adiabatic channel capture of the reactants is accompanied by partial primary redissociation of the N⋯OH collision pairs. A criterion for the extent of this primary redissociation in competition to the formation of randomized, long-lived, complex of NOH is proposed. The NOH adducts then may decompose to NO + H, rearrange in a unimolecular process to HNO, or undergo secondary redissociation back to the reactants N + OH, while HNO may also decompose to NO + H. As the reactants N(⁴S) + OH(²Π) have open electronic shells, non-Born-Oppenheimer effects have to be considered. Their influence on thermal rate constants of the reaction at low temperatures is illustrated and compared with such effects in other reactions such as C(³P) + OH(²Π).

Thermal Decomposition of 3-Bromopropene. A Theoretical Kinetic Investigation

A detailed kinetic study of the gas-phase thermal decomposition of 3-bromopropene over wide temperature and pressure ranges was performed. Quantum chemical calculations employing the density functional theory methods B3LYP, BMK, and M06-2X and the CBS-QB3 and G4 ab initio composite models provide the relevant part of the potential energy surfaces and the molecular properties of the species involved in the CH2═CH-CH2Br → CH2═C═CH2 + HBr (1) and CH2═CH-CH2Br → CH2═CH-CH2 + Br (2) reaction channels. Transition-state theory and unimolecular reaction rate theory calculations show that the simple bond fission reaction ( 2 ) is the predominant decomposition channel and that all reported experimental studies are very close to the high-pressure limit of this process. Over the 500-1400 K range a rate constant for the primary dissociation of k2,∞ = 4.8 × 10(14) exp(-55.0 kcal mol(-1)/RT) s(-1) is predicted at the G4 level. The calculated k1,∞ values lie between 50 to 260 times smaller. A value of 10.6 ± 1.5 kcal mol(-1) for the standard enthalpy of formation of 3-bromopropene at 298 K was estimated from G4 thermochemical calculations.

Vibrational excitation influence on the H + BrO → HBr + O reaction: a quasi-classical trajectory investigation

Quasi-classical trajectory calculations are employed to investigate the vibrational excitation effect on the scalar and vector properties of the H + BrO → HBr + O reaction using a X1A′ state ab initio potential energy surface (J. Chem. Phys. 2000, 113, 4598). The reaction probability, cross section, and rate constant are carried out with the effect of the collision energy (Ecol = 0.1–6 kcal/mol) and vibrational levels (v = 0–3). A significant vibrational dependency has been observed in the reaction probability and cross section at a relatively low collision energy area and has also been found in a low-temperature (T < 150 K) region of the rate constant. In addition, two product angular distributions, P(θr) and P(ϕr), and two generalized polarization-dependent differential cross sections, PDDCS00 and PDDCS20, are calculated as well. All of these scalar and vector properties have shown sensitive behaviors to the vibrational levels.

Electronic nonadiabatic effects in low temperature radical-radical reactions. I. C(3P) + OH(2Π)

The formation of collision complexes, as a first step towards reaction, in collisions between two open-electronic shell radicals is treated within an adiabatic channel approach. Adiabatic channel potentials are constructed on the basis of asymptotic electrostatic, induction, dispersion, and exchange interactions, accounting for spin-orbit coupling within the multitude of electronic states arising from the separated reactants. Suitable coupling schemes (such as rotational + electronic) are designed to secure maximum adiabaticity of the channels. The reaction between C((3)P) and OH((2)Π) is treated as a representative example. The results show that the low temperature association rate coefficients in general cannot be represented by results obtained with a single (generally the lowest) potential energy surface of the adduct, asymptotically reaching the lowest fine-structure states of the reactants, and a factor accounting for the thermal population of the latter states. Instead, the influence of non-Born-Oppenheimer couplings within the multitude of electronic states arising during the encounter markedly increases the capture rates. This effect extends up to temperatures of several hundred K.

Global potential energy surface, vibrational spectrum, and reaction dynamics of the first excited (Ã (2)A(')) state of HO(2)

The authors report extensive high-level ab initio studies of the first excited (Ã (2)A(')) state of HO(2). A global potential energy surface (PES) was developed by spline-fitting 17 000 ab initio points at the internal contracted multireference configuration interaction (icMRCI) level with the AVQZ basis set. To ascertain the spectroscopic accuracy of the PES, the near-equilibrium region of the molecule was also investigated using three interpolating moving least-squares-based PESs employing dynamically weighted icMRCI methods in the complete basis set limit. Vibrational energy levels on all four surfaces agree well with each other and a new assignment of some vibrational features is proposed. In addition, the dynamics of both the forward and reverse directions of the H+O(2)(ã (1)Δ(g))↔OH+O reaction (J=0) were studied using an exact wave packet method. The reactions are found to be dominated by sharp resonances.

Contribution of the radical-complex mechanism to the rate of the reaction CH(3) + O(2) (+ M) → CH(3)O(2) (+ M) at high pressures

Earlier experimental studies of the falloff curves of the reaction CH(3) + O(2) (+ M) → CH(3)O(2) (+ M) in the bath gases M = Ar and N(2) (Fernandes et al., J. Phys. Chem. A 2006, 110, 4442), in addition to the usual behavior of the energy-transfer (ET) mechanism, showed first evidence for a participation of the radical-complex (RC) mechanism in the reaction at pressures above about 300 bar and at temperatures below 400 K. By extending these measurements to the bath gas M = CO(2), more pronounced deviations from the ET mechanism were now observed. This unambiguously confirms the presence of the RC mechanism at high pressures in a medium-sized molecular system, analogous to earlier observations for larger systems such as the dimerization of benzyl radicals (Luther et al., Phys. Chem. Chem. Phys. 2004, 6, 4133).

Classical Trajectory and Statistical Adiabatic Channel Study of the Dynamics of Capture and Unimolecular Bond Fission. VII. Thermal Capture and Specific Rate Constants k(E,J) for the Dissociation of Molecular Ions

Specific rate constants, k(E,J), and thermal capture rate constants, k cap(T), are determined by statistical adiabatic channel model/classical trajectory (SACM/CT) calculations for unimolecular dissociation and the reverse association reactions of representative polyatomic molecular ions. Simple short-range valence/long-range ion-induced dipole model potentials without reverse barriers have been employed, using the reactions C8H10 + ⇔ C7H7 + + CH3 and C9H12 + ⇔ C7H7 + + C2H5 as illustrative examples. Simplified representations of k(E) and k cap(T) from rigid activated complex Rice–Ramsperger–Kassel–Marcus (RRKM) theory are compared with the SACM/CT treatment and with experimental results. The Massey parameters of the transitional mode dynamics, for the systems considered, are smaller than unity such that their dynamics is nonadiabatic while the dynamics of the conserved modes is adiabatic. Because of the long-range/short-range switching character of the potential, simple rigid activated complex RRKM theory cannot be used without modifications. The effects of a shifting of the effective bottle-neck of the dynamics with increasing energy towards smaller interfragment distances in the present cases are amplified by a shift into a range of increasing anisotropy of the potential. As a consequence, the thermal capture rate constants markedly decrease with increasing temperature.

Rate constant for OH(2 Pi)+O(3P)-->H(2S)+O2(3 Sigma g-) reaction on an improved ab initio potential energy surface and implications for the interstellar oxygen problem

The authors report a global potential energy surface for the ground electronic state of HO(2)(X (2)A(")), which improves upon the XXZLG potential [Xu and et al., J. Chem. Phys. 122, 244305 (2005)] with additional high-level ab initio points for the long-range interaction potential in the O+OH channel. Exact J=0 quantum mechanical reaction probabilities were calculated on the new potential and the rate constant for the title reaction was obtained using a J-shifting method. The calculated rate constant is in good agreement with available experimental values and our results predict a significantly lower rate at temperature range below 30 K, offering a possible explanation for the "interstellar oxygen problem."

Pressure and temperature dependence of the recombination of p-fluorobenzyl radicals

The rate constants of the recombination reaction of p-fluorobenzyl radicals, p-F-C6H4CH2 + p-F-C6H4CH2 (+M) --> C14H12F2 (+M), have been measured over the pressure range 0.2-800 bar and the temperature range 255-420 K. Helium, argon, and CO2 were employed as bath gases (M). At pressures below 0.9 bar in Ar and CO2, and 40 bar in He, the rate constant k1 showed no dependence on the pressure and the nature of the bath gas, clearly indicating that it had reached the limiting high-pressure value of the energy-transfer (ET) mechanism (k(1,infinity)ET). A value of k(1,infinity)ET(T) = (4.3 +/- 0.5) x 10(-11) (T/300 K)(-0.2) cm3 molecule(-1) s(-1) was determined. At pressures above about 5 bar, the k1 values in Ar and CO2 were found to gradually increase in a pressure range where according to energy-transfer mechanism, they should remain at the constant value k(1,infinity)ET. The enhancement of the recombination rate constant beyond the value k(1,infinity)ET increased in the order He < Ar < CO2, and it became more pronounced with decreasing temperature. The dependences of k1 on pressure, temperature, and the bath gas were similar to previous observations in the recombination of benzyl radicals. The effect of fluorine-substitution of the benzyl ring on k1 values is discussed. The present results confirm the significant role of radical complexes in the recombination kinetics of benzyl-type radicals in the gas-liquid transition range. The observations on a rate enhancement beyond the experimental value of k(1,infinity)ET at elevated densities up to the onset of diffusion-control are consistently explained by the kinetic contribution of a "radical-complex" mechanism which is solely based on standard van der Waals interaction between radicals and bath gases.

Falloff Curves for the Recombination Reaction Cl + FC(O)O + M → FC(O)OCl + M

The pressure dependence of the recombination reaction Cl + FC(O)O + M --> FC(O)OCl + M has been investigated at 296 K. FC(O)O radicals and Cl atoms were generated by laser flash photodissociation of FC(O)OO(O)CF at 193 nm in mixtures with Cl2 and He or SF6 over the total pressure range 8-645 Torr. The measured FC(O)O radical and F atom yields in the photolysis are 0.33 +/- 0.06 and 0.67 +/- 0.06. The reaction lies in the falloff range approaching the high-pressure limit. The extrapolations toward the limiting low- and high-pressure ranges were carried out using a reduced falloff curves formalism, which includes a recent implementation for the strong-collision broadening factors. The resulting values for the low-pressure rate coefficients are (2.2 +/- 0.4) x 10(-28)[He], (4.9 +/- 0.9) x 10(-28)[SF6], (1.9 +/- 0.3) x 10(-28)[Cl2] and (5.9 +/- 1.1) x 10(-28)[FC(O)OO(O)CF] cm3 molecule(-1) s(-1). The derived high-pressure rate coefficient is (4.4 +/- 0.8) x 10(-11) cm3 molecule(-1) s(-1). For the reaction Cl + FC(O)OCl --> Cl2 + FC(O)O a rate coefficient of (1.6 +/- 0.3) x 10(-11) cm3 molecule(-1) s(-1) was determined. The high-pressure rate coefficient was theoretically interpreted using SACM/CT calculations on an ab initio electronic potential computed at the G3S level of theory. Standard heat of formation values of -99.9 and -102.5 kcal mol(-1) were computed at the G3//B3LYP/6-311++G(3df,3pd) level of theory for cis-FC(O)OCl and trans-FC(O)OCl, respectively. The computed electronic barrier for the conversion between the trans and cis conformers is 8.9 kcal mol(-1). On the basis of the present results, the above reactions are expected to have a negligible impact on stratospheric ozone levels.

Reaction kinetics of the addition of atomic sulfur to nitric oxide

The reaction of S((3)P(J)) with NO ((2)Pi) in an Ar bath gas has been studied by the laser photolysis-resonance fluorescence technique over 300-810 K at pressures from 60 to 800 mbar. The observed second-order rate constants are close to the low-pressure limit. Fitting of Troe's formalism to experiment, with an estimated F(cent) = 0.78 exp(-T/7445) and k(infinity) given subsequently, yields k(0) = (6.2+/-0.6) x 10(-33) exp(+ (940+/-40)/T) cm(6) molecule(-2) s(-1). Error limits are +/-25%. A theoretical analysis of this value suggests that the average energy transferred during collisions between Ar and the excited intermediate is DeltaE = -360(-160) (+90) cm(-1). Over 300-800 K, the high-pressure limit is predicted to be k(infinity) = 2.2 x 10(-10) (T/300)(0.24) cm(3) molecule(-1) s(-1). Doublet and quartet adducts between S and NO were characterized via CBS-QB3 theory. The kinetic data can be rationalized with SNO ((2)A(')) as the major product, and an ab initio estimate of Delta(f)H(298) for SNO is 176+/-8 kJ mol(-1).

Correlation of Adiabatic States between Perturbed Rotor and Renner-Teller Limits for a Closed-Shell Ion + Open Shell Diatom System

The adiabatic channel interaction potentials for a system charge + dipole-open-shell-quadrupole linear rotor are calculated in the perturbed rotor and slightly anharmonic Renner-Teller limits. In both cases, first order charge-dipole and charge-quadrupole interactions are not additive; this leads to novel expressions for the energy levels. The result can be used in the construction of the adiabatic correlation diagrams between free rotor states with account for the