Experimental and reduced dimensionality quantum rate coefficients for H2(D2)+CN→H(D)CN+H(D)

Energetics and kinetics of various cyano radical hydrogen abstractions

The cyano radical (CN) is an abundant, open-shell molecule found in a variety of environments, including the atmosphere, the interstellar medium and combustion processes. In these environments, it often reacts with small, closed-shell molecules via hydrogen abstraction. Both carbon and nitrogen atoms of the cyano radical are reactive sites, however the carbon is more reactive with reaction barrier heights generally between 2-15 kcal mol^-1 lower than those of the analogous nitrogen. The CN + HX → HCN/HNC + X, with X = H, CH₃, NH₂, OH, F, SiH₃, PH₂, SH, Cl, C₂H, CN reactions have been studied at a high-level of theory, including CCSD(T)-F12a. Furthermore, kinetics were obtained over the 100-1000 K temperature range, showing excellent agreement with those rate constants that have been determined experimentally.

A general formulation of the quasiclassical trajectory method for reduced-dimensionality reaction dynamics calculations

Dimension reduction by freezing the unimportant coordinates is widely used in intramolecular and reaction dynamics calculations when the solution of the accurate full-dimensional nuclear Schrödinger equation is not feasible. In this paper we report on a novel form of the exact classical internal-coordinate Hamiltonian for full and reduced-dimensional vibrational motion of polyatomic molecules with the purpose of using it in quasiclassical trajectory (QCT) calculations. The derivation is based on the internal to body-fixed frame transformation, as in the t-vector formalism, however it does not require the introduction of rotational variables to allow cancellation of non-physical rotations within the body-fixed frame. The formulas needed for QCT calculations: normal mode analysis and state sampling as well as for following the dynamics and normal-mode quantum number assignment at instantaneous states are presented. The procedure is demonstrated on the CH4, CD4, CH3D and CHD3 isotopologs of methane using three reduced-dimensional models, which were previously used in quantum reactive scattering studies of the CH4 + X → CH3 + HX type reactions. The reduced-dimensional QCT methodology formulated this way combined with full-dimensional QCT calculations makes possible the classical validation of reduced-dimensional models that are used in the quantum mechanical description of the nuclear dynamics in reactive systems [A. Vikár et al., J. Phys. Chem. A, 2016, 120, 5083-5093].

Kinetics study of the CN + CH4 hydrogen abstraction reaction based on a new ab initio analytical full-dimensional potential energy surface

We have developed an analytical full-dimensional potential energy surface, named PES-2017, for the gas-phase hydrogen abstraction reaction between the cyano radical and methane.

Deformed transition-state theory: Deviation from Arrhenius behavior and application to bimolecular hydrogen transfer reaction rates in the tunneling regime

A formulation is presented for the application of tools from quantum chemistry and transition-state theory to phenomenologically cover cases where reaction rates deviate from Arrhenius law at low temperatures. A parameter d is introduced to describe the deviation for the systems from reaching the thermodynamic limit and is identified as the linearizing coefficient in the dependence of the inverse activation energy with inverse temperature. Its physical meaning is given and when deviation can be ascribed to quantum mechanical tunneling its value is calculated explicitly. Here, a new derivation is given of the previously established relationship of the parameter d with features of the barrier in the potential energy surface. The proposed variant of transition state theory permits comparison with experiments and tests against alternative formulations. Prescriptions are provided and implemented to three hydrogen transfer reactions: CH₄  + OH → CH₃  + H₂ O, CH₃ Cl + OH → CH₂ Cl + H₂ O and H₂  + CN → H + HCN, widely investigated both experimentally and theoretically. © 2016 Wiley Periodicals, Inc.

A THEORETICAL STUDY FOR H2 + CN ↔ HCN + H REACTION AND ITS KINETIC ISOTOPE EFFECTS WITH VARIATIONAL TRANSITION STATE THEORY

We present variational transition state theory (VTST) calculations for the H 2 + CN → HCN + H (R1) and D 2 + CN → DCN + D (R2) reactions and their reverses based on a global many-body expansion potential energy surface (PES) for ground-state H 2 CN (ter Horst MA, Schatz GC, Harding LB, J Chem Phys105:558, 1996). It is found that the tunneling effects are negligible over the 200–2000 K temperature range and non-negligible over 100–200 K for R1 and R2 reactions. The C–N bond acts almost as a spectator for both reactions. The present VTST rate constants are in good agreement with the available experimental results and the previous theoretical predictions for R1 and R2 reactions except for the overestimation of rate constants by VTST at lower temperatures that may be caused by recrossing effect. Additionally, the kinetic isotope effects are important for the forward R1 and R2 reactions, but not for the reverses of R1 and R2.

Dynamics of CN+alkane reactions by crossed-beam dc slice imaging

The hydrogen atom abstraction reactions of CN (X (2)Sigma(+)) with alkanes have been studied using the crossed molecular beam technique with dc slice ion imaging at collision energies of 7.5 and 10.8 kcalmol. The product alkyl radical images were obtained via single photon ionization at 157 nm for the reactions of CN (X (2)Sigma(+)) with n-butane, n-pentane, n-hexane, and cyclohexane. From analysis of the images, we obtained the center-of-mass frame product angular distributions and translational energy distributions directly. The results indicate that the products are largely backscattered and that most of the available energy ( approximately 80%-85%) goes to the internal energy of the products. The reaction dynamics is discussed in light of recent kinetics data, theoretical calculations, and results for related halogen and oxygen atom reactions.

H Atom Yields from the Reactions of CN Radicals with C2H2, C2H4, C3H6, trans-2-C4H8, and iso-C4H8

The kinetics and H atom channel yield at both 298 and 195 K have been determined for reactions of CN radicals with C2H2 (1.00+/-0.21, 0.97+/-0.20), C2H4 (0.96+/-0.032, 1.04+/-0.042), C3H6 (pressure dependent), iso-C4H8 (pressure dependent), and trans-2-C4H8 (0.039+/-0.019, 0.029+/-0.047) where the first figure in each bracket is the H atom yield at 298 K and the second is that at 195 K. The kinetics of all reactions were studied by monitoring both CN decay and H atom growth by laser-induced fluorescence at 357.7 and 121.6 nm, respectively. The results are in good agreement with previous studies where available. The rate coefficients for the reaction of CN with trans-2-butene and iso-butene have been measured at 298 and 195 K for the first time, and the rate coefficients are as follows: k298K=(2.93+/-0.23)x10(-10) cm3 molecule(-1) s(-1), k195K=(3.58+/-0.43)x10(-10) cm3 molecule(-1) s(-1) and k298K=(3.17+/-0.10)x10(-10) cm3 molecule(-1) s(-1), k195K=(4.32+/-0.35)x10(-10) cm3 molecule(-1) s(-1), respectively, where the errors represent a combination of statistical uncertainty (2sigma) and an estimate of possible systematic errors. A potential energy surface for the CN+C3H6 reaction has been constructed using G3X//UB3LYP electronic structure calculations identifying a number of reaction channels leading to either H, CH3, or HCN elimination following the formation of initial addition complexes. Results from the potential energy surface calculations have been used to run master equation calculations with the ratio of primary:secondary addition, the average amount of downward energy transferred in a collision DeltaEd, and the difference in barrier heights between H atom elimination and an H atom 1, 2 migration as variable parameters. Excellent agreement is obtained with the experimental 298 K H atom yields with the following parameter values: secondary addition complex formation equal to 80%, DeltaEd=145 cm(-1), and the barrier height for H atom elimination set 5 kJ mol(-1) lower than the barrier for migration. Finally, very low temperature master equation simulations using the best fit parameters have been carried out in an increased precision environment utilizing quad-double and double-double arithmetic to predict H and CH3 yields for the CN+C3H6 reaction at temperatures and pressures relevant to Titan. The H and CH3 yields predicted by the master equation have been parametrized in a simple equation for use in modeling.

Interpolating moving least-squares methods for fitting potential energy surfaces: An application to the H2CN unimolecular reaction

Classical trajectories have been used to compute rates for the unimolecular reaction H2CN-->H+HCN on a fitted ab initio potential energy surface (PES). The ab initio energies were obtained from CCSD(T)/aug-cc-pvtz electronic structure calculations. The ab initio energies were fitted by the interpolating moving least-squares (IMLS) method. This work continues the development of the IMLS method for producing ab initio PESs for use in molecular dynamics simulations of many-atom systems. A dual-level scheme was used in which the preliminary selection of data points was done using a low-level theory and the points used for fitting the final PES were obtained at the desired higher level of theory. Classical trajectories were used on various low-level IMLS fits to tune the fit to the unimolecular reaction under study. Procedures for efficiently picking data points, selecting basis functions, and defining cutoff limits to exclude distant points were investigated. The accuracy of the fitted PES was assessed by comparing interpolated values of quantities to the corresponding ab initio values. With as little as 330 ab initio points classical trajectory rate constants were converged to 5%-10% and the rms error over the six-dimensional region sampled by the trajectories was a few tenths of a kcal/mol.

Infrared laser spectroscopy of the CH3–HCN radical complex stabilized in helium nanodroplets

The CH3-HCN and CD3-HCN radical complexes have been formed in helium nanodroplets by sequential pickup of a CH3 (CD3) radical and a HCN molecule and have been studied by high-resolution infrared laser spectroscopy. The complexes have a hydrogen-bonded structure with C3v symmetry, as inferred from the analysis of their rotationally resolved nu = 1 <-- 0 H-CN vibrational bands. The A rotational constants of the complexes are found to change significantly upon vibrational excitation of the C-H stretch of HCN within the complex, DeltaA = A'-A" = -0.04 cm(-1) (for CH3-HCN), whereas the B rotational constants are found to be 2.9 times smaller than that predicted by theory. The reduction in B can be attributed to the effects of helium solvation, whereas the large DeltaA is found to be a sensitive probe of the vibrational averaging dynamics of such weakly bound systems. The complex has a permanent electric dipole moment of 3.1 +/- 0.2 D, as measured by Stark spectroscopy. A vibration-vibration resonance is observed to couple the excited C-H stretching vibration of HCN within the complex to the lower-frequency C-H stretches of the methyl radical. Deuteration of the methyl radical was used to detune these levels from resonance, increasing the lifetime of the complex by a factor of 2. Ab initio calculations for the energies and molecular parameters of the stationary points on the CN+CH4 --> HCN+CH3 potential-energy surface are also presented.