Reaction-path dynamics calculations of the Cl + NH(3) hydrogen abstraction reaction: the role of the intermediate complexes

Benchmark ab initio characterization of the multi-channel Cl + CH3X [X = F, Cl, Br, I] reactive potential energy surfaces

We determine benchmark geometries and relative energies for the stationary points of the Cl + CH3X [X = F, Cl, Br, I] reactions. We consider four possible reaction pathways: hydrogen abstraction, hydrogen substitution, halogen abstraction, and halogen substitution, where the substitution processes can proceed via either Walden inversion or front-side attack. We perform geometry optimizations and obtain harmonic vibrational frequencies at the explicitly-correlated UCCSD(T)-F12b/aug-cc-pVTZ level of theory, followed by UCCSD(T)-F12b/aug-cc-pVQZ single-point computations to make finite-basis-set error negligible. To reach chemical (<1 kcal mol-1), or even subchemical (<0.5 kcal mol-1) accuracy, we include core-correlation, scalar relativistic, post-(T), spin-orbit-splitting and zero-point-energy contributions, as well, in the relative energies of all the stationary points. Our benchmark 0 K reaction enthalpies are compared to available experimental results and show good agreement. The stationary-point structures and energetics are interpreted in terms of Hammond's postulate and used to make predictions related to the dynamical behavior of these reactive systems.

Theoretical Insights into the Dynamics of Gas-Phase Bimolecular Reactions with Submerged Barriers

Much attention has been paid to the dynamics of both activated gas-phase bimolecular reactions, which feature monotonically increasing integral cross sections and Arrhenius kinetics, and their barrierless capture counterparts, which manifest monotonically decreasing integral cross sections and negative temperature dependence of the rate coefficients. In this Perspective, we focus on the dynamics of gas-phase bimolecular reactions with submerged barriers, which often involve radicals or ions and are prevalent in combustion, atmospheric chemistry, astrochemistry, and plasma chemistry. The temperature dependence of the rate coefficients for such reactions is often non-Arrhenius and complex, and the corresponding dynamics may also be quite different from those with significant barriers or those completely dominated by capture. Recent experimental and theoretical studies of such reactions, particularly at relatively low temperatures or collision energies, have revealed interesting dynamical behaviors, which are discussed here. The new knowledge enriches our understanding of the dynamics of these unusual reactions.

Benchmark Ab Initio Characterization of the Abstraction and Substitution Pathways of the Cl + CH3CN Reaction

We investigate the reaction pathways of the Cl + CH₃CN system: hydrogen abstraction, methyl substitution, hydrogen substitution, and cyanide substitution, leading to HCl + CH₂CN, ClCN/CNCl + CH₃, ClCH₂CN + H, and CH₃Cl + CN, respectively. Hydrogen abstraction is exothermic and has a low barrier, whereas the other channels are endothermic with high barriers. The latter two can proceed via a Walden inversion or front-side attack mechanism, and the front-side attack barriers are always higher. The C-side methyl substitution has a lower barrier and also a lower endothermicity than the N-side reaction. The computations utilize an accurate composite ab initio approach and the explicitly correlated CCSD(T)-F12b method. The benchmark classical and vibrationally adiabatic energies of the stationary points are determined with the most accurate CCSD(T)-F12b/aug-cc-pVQZ energies adding further contributions of the post-(T) and core correlation, scalar relativistic effects, spin–orbit coupling, and zero-point energy corrections. These contributions are found to be non-negligible to reach subchemical accuracy.

Quasi-Classical Trajectory Study of the CN + NH3 Reaction Based on a Global Potential Energy Surface

Based on a combination of valence-bond and molecular mechanics functions which were fitted to high-level ab initio calculations, we constructed an analytical full-dimensional potential energy surface, named PES-2020, for the hydrogen abstraction title reaction for the first time. This surface is symmetrical with respect to the permutation of the three hydrogens in ammonia, it presents numerical gradients and it improves the description presented by previous theoretical studies. In order to analyze its quality and accuracy, stringent tests were performed, exhaustive kinetics and dynamics studies were carried out using quasi-classical trajectory calculations, and the results were compared with the available experimental evidence. Firstly, the properties (geometry, vibrational frequency and energy) of all stationary points were found to reasonably reproduce the ab initio information used as input; due to the complicated topology with deep wells in the entrance and exit channels and a “submerged” transition state, the description of the intermediate complexes was poorer, although it was adequate to reasonably simulate the kinetics and dynamics of the title reaction. Secondly, in the kinetics study, the rate constants simulated the experimental data in the wide temperature range of 25–700 K, improving the description presented by previous theoretical studies. In addition, while previous studies failed in the description of the kinetic isotope effects, our results reproduced the experimental information. Finally, in the dynamics study, we analyzed the role of the vibrational and rotational excitation of the CN(v,j) reactant and product angular scattering distribution. We found that vibrational excitation by one quantum slightly increased reactivity, thus reproducing the only experimental measurement, while rotational excitation strongly decreased reactivity. The scattering distribution presented a forward-backward shape, associated with the presence of deep wells along the reaction path. These last two findings await experimental confirmation.

Theoretical kinetic study of the formic acid catalyzed Criegee intermediate isomerization: multistructural anharmonicity and atmospheric implications

We performed a theoretical study on the double hydrogen shift isomerization reaction of a six carbon atom Criegee intermediate (C6-CI), catalyzed by formic acid (HCOOH), to produce vinylhydroperoxide (VHP), C6-CI + HCOOH → VHP + HCOOH. This Criegee intermediate can serve as a surrogate for larger CIs derived from important volatile organic compounds like monoterpenes, whose reactivity is not well understood and which are difficult to handle computationally. The reactant HCOOH exerts a pronounced catalytic effect on the studied reaction by lowering the barrier height, but the kinetic enhancement is hindered by the multistructural anharmonicity. First, the rigid ring-structure adopted by the saddle point to facilitate simultaneous transfer of two atoms does not allow the formation of as many conformers as those formed by the reactant C6-CI. And second, the flexible carbon chain of C6-CI facilitates the formation of stabilizing intramolecular C-HO hydrogen bonds; this stabilizing effect is less pronounced in the saddle point structure due to its tightness and steric effects. Thus, the contribution of the reactant C6-CI conformers to the multistructural partition function is larger than that of the saddle point conformers. The resulting low multistructural anharmonicity factor partially cancels out the catalytic effect of the carboxylic acid, yielding in a moderately large rate coefficient, k(298 K) = 4.9 × 10-13 cm3 molecule-1 s-1. We show that carboxylic acids may promote the conversion of stabilized Criegee intermediates into vinylhydroperoxides in the atmosphere, which generates OH radicals and leads to secondary organic aerosols, thereby affecting the oxidative capacity of the atmosphere and ultimately the climate.

Ab initio and transition state theory study of the OH + HO2 → H2O + O2(3Σg-)/O2(1Δg) reactions: yield and role of O2(1Δg) in H2O2 decomposition and in combustion of H2

Reactions of hydroxyl (OH) and hydroperoxyl (HO₂) are important for governing the reactivity of combustion systems. We performed post-CCSD(T) ab initio calculations at the W3X-L//CCSD = FC/cc-pVTZ level to explore the triplet ground-state and singlet excited-state potential energy surfaces of the OH + HO₂ → H₂O + O₂(³Σ_g^-)/O₂(¹Δ_g) reactions. Using microcanonical and multistructural canonical transition state theories, we calculated the rate constant for the triplet and singlet channels over the temperature range 200-2500 K, represented by k(T) = 3.08 × 10¹²T^0.07 exp(1151/RT) + 8.00 × 10¹²T^0.32 exp(-6896/RT) and k(T) = 2.14 × 10⁶T^1.65 exp(-2180/RT) in cm³ mol^-1 s^-1, respectively. The branching ratios show that the yield of singlet excited oxygen is small (<0.5% below 1000 K). To ascertain the importance of singlet oxygen channel, our new kinetic information was implemented into the kinetic model for hydrogen combustion recently updated by Konnov (Combust. Flame, 2015, 162, 3755-3772). The updated kinetic model was used to perform H₂O₂ thermal decomposition simulations for comparison against shock tube experiments performed by Hong et al. (Proc. Combust. Inst., 2013, 34, 565-571), and to estimate flame speeds and ignition delay times in H₂ mixtures. The simulation predicted a larger amount of O₂(¹Δ_g) in H₂O₂ decomposition than that predicted by Konnov's original model. These differences in the O₂(¹Δ_g) yield are due to the use of a higher ab initio level and a more sophisticated methodology to compute the rate constant than those used in previous studies, thereby predicting a significantly larger rate constant. No effect was observed on the rate of the H₂O₂ decomposition and on the flame speeds and ignition delay times of different H₂-oxidizer mixtures. However, if the oxidizer is seeded with O₃, small differences appear in the flame speed. Given that O₂(¹Δ_g) is much more reactive than O₂(³Σ_g^-), we do not preclude an effect of the singlet channel of the titled reaction in other combustion systems, especially in systems where excited oxygen plays an important role.

Dynamics study of the OH + NH3 hydrogen abstraction reaction using QCT calculations based on an analytical potential energy surface

To understand the reactivity and mechanism of the OH + NH3 → H2O + NH2 gas-phase reaction, which evolves through wells in the entrance and exit channels, a detailed dynamics study was carried out using quasi-classical trajectory calculations. The calculations were performed on an analytical potential energy surface (PES) recently developed by our group, PES-2012 [Monge-Palacios et al. J. Chem. Phys. 138, 084305 (2013)]. Most of the available energy appeared as H2O product vibrational energy (54%), reproducing the only experimental evidence, while only the 21% of this energy appeared as NH2 co-product vibrational energy. Both products appeared with cold and broad rotational distributions. The excitation function (constant collision energy in the range 1.0-14.0 kcal mol(-1)) increases smoothly with energy, contrasting with the only theoretical information (reduced-dimensional quantum scattering calculations based on a simplified PES), which presented a peak at low collision energies, related to quantized states. Analysis of the individual reactive trajectories showed that different mechanisms operate depending on the collision energy. Thus, while at high energies (E(coll) ≥ 6 kcal mol(-1)) all trajectories are direct, at low energies about 20%-30% of trajectories are indirect, i.e., with the mediation of a trapping complex, mainly in the product well. Finally, the effect of the zero-point energy constraint on the dynamics properties was analyzed.

Ab initio based potential energy surface and kinetics study of the OH + NH3 hydrogen abstraction reaction

A full-dimensional analytical potential energy surface (PES) for the OH + NH3 → H2O + NH2 gas-phase reaction was developed based exclusively on high-level ab initio calculations. This reaction presents a very complicated shape with wells along the reaction path. Using a wide spectrum of properties of the reactive system (equilibrium geometries, vibrational frequencies, and relative energies of the stationary points, topology of the reaction path, and points on the reaction swath) as reference, the resulting analytical PES reproduces reasonably well the input ab initio information obtained at the coupled-cluster single double triple (CCSD(T)) = FULL/aug-cc-pVTZ//CCSD(T) = FC/cc-pVTZ single point level, which represents a severe test of the new surface. As a first application, on this analytical PES we perform an extensive kinetics study using variational transition-state theory with semiclassical transmission coefficients over a wide temperature range, 200-2000 K. The forward rate constants reproduce the experimental measurements, while the reverse ones are slightly underestimated. However, the detailed analysis of the experimental equilibrium constants (from which the reverse rate constants are obtained) permits us to conclude that the experimental reverse rate constants must be re-evaluated. Another severe test of the new surface is the analysis of the kinetic isotope effects (KIEs), which were not included in the fitting procedure. The KIEs reproduce the values obtained from ab initio calculations in the common temperature range, although unfortunately no experimental information is available for comparison.

Constructing Potential Energy Surfaces for Polyatomic Systems: Recent Progress and New Problems

Different methods of constructing potential energy surfaces in polyatomic systems are reviewed, with the emphasis put on fitting, interpolation, and analytical (defined by functional forms) approaches, based on quantum chemistry electronic structure calculations. The different approaches are reviewed first, followed by a comparison using the benchmark H + CH4 and the H + NH3 gas-phase hydrogen abstraction reactions. Different kinetics and dynamics properties are analyzed for these reactions and compared with the available experimental data, which permits one to estimate the advantages and disadvantages of each method. Finally, we analyze different problems with increasing difficulty in the potential energy construction: spin-orbit coupling, molecular size, and more complicated reactions with several maxima and minima, which test the soundness and general applicability of each method. We conclude that, although the field of small systems, typically atom-diatom, is mature, there still remains much work to be done in the field of polyatomic systems.

QCT and QM calculations of the Cl(2P) + NH3 reaction: influence of the reactant well on the dynamics

A detailed dynamics study, using both quasi-classical trajectory (QCT) and reduced-dimensional quantum mechanical (QM) calculations, was carried out to understand the reactivity and mechanism of the Cl((2)P) + NH(3)→ HCl + NH(2) gas-phase reaction, which evolves through deep wells in the entry and exit channels. The calculations were performed on an analytical potential energy surface recently developed by our group, PES-2010 [M. Monge-Palacios, C. Rangel, J. C. Corchado and J. Espinosa-Garcia, Int. J. Quantum. Chem., 2011], together with a simplified model surface, mod-PES, in which the reactant well is removed to analyze its influence. The main finding was that the QCT and QM methods show a change of the reaction probability with collision energy, suggesting a change of the atomic-level mechanism of reaction with energy. This change disappeared when the mod-PES was used, showing that the behaviour at low energies is a direct consequence of the existence of the reactant well. Analysis of the trajectories showed that different mechanisms operate depending on the collision energy. Thus, while at high energies (E(coll) > 5 kcal mol(-1)) practically all trajectories are direct, at low energies (E(coll) < 3 kcal mol(-1)) the trajectories are indirect, i.e., with the mediation of a trapping complex in the entry and/or the exit wells. The reactant complex allows repeated encounters between the reactants, increasing the reaction probability at low energies. The differential cross section results reinforce this change of mechanism, showing also the influence of the reactant well on this reaction. Thus, the PES-2010 surface yields a forward-backward symmetry in the scattering, while when the reactant well is removed with the mod-PES the shape is more isotropic.