Accurate ab initio based potential energy surface and kinetics of the Cl + NH3 → HCl + NH2 reaction

The gas-phase reaction Cl + NH3 → HCl + NH2 is a prototypical hydrogen abstraction reaction, whose minimum energy path involves several intermediate complexes. In this work, a full-dimensional, spin-orbit corrected potential energy surface (SOC PES) is constructed for the ground electronic state of the Cl + NH3 reaction. About 52 000 energy points are sampled and calculated at the UCCSD(T)-F12a/aug-cc-pVTZ level, in which the data points located in the entrance channel are spin-orbit corrected. The spin-orbit corrections are predicted by a fitted three-dimensional energy surface from about 7520 energy points in the entrance channel at the level of CASSCF (15e, 11o)/aug-cc-pVTZ. The fundamental-invariant neural network method is utilized to fit the SOC PES, resulting in a total root mean square error of 0.12 kcal mol-1. The calculated thermal rate constants of the Cl + NH3 → HCl + NH2 reaction on the SOC PES with the soft-zero-point energy constraint agree reasonably well with the available experimental values.

Role of the Vibrational and Translational Energies in the CN(v)+C2H6(ν1, ν2, ν5 and ν9) Reactions. A Theoretical QCT Study

Quasi-classical trajectory (QCT) calculations were conducted on the newly developed full-dimensional potential energy surface, PES-2023, to analyse two critical aspects: the influence of vibrational versus translational energy in promoting reactivity, and the impact of vibrational excitation within similar vibrational modes. The former relates to Polanyi's rules, while the latter concerns mode selectivity. Initially, the investigation revealed that independent vibrational excitation by a single quantum of ethane's symmetric and asymmetric stretching modes (differing by only 15 cm-1) yielded comparable dynamics, reaction cross-sections, HCN(v) vibrational product distributions, and scattering distributions. This observation dismisses any significant mode selectivity. Moreover, an equivalent amount of energy provided as translational energy (at total energies of 9.6 and 20.0 kcal mol-1) gave rise to slightly lower reactivity compared to the same amount of energy provided as vibrational energy. This effect is more evident at low energies, presenting a counterintuitive scenario in an 'early transition state' reaction. These findings challenge the straightforward application of Polanyi's rules in polyatomic systems. Regarding CN(v) vibrational excitation, our calculations reveal that the reaction cross-section remains practically unaffected by this vibrational excitation, suggesting that the CN stretching mode is a spectator mode. The results were rationalized by considering several factors: the strong coupling between different vibrational modes, and between vibrational modes and the reaction coordinate; and a significant vibrational energy redistribution within the ethane reactant before collision. This redistribution creates an unphysical energy flow, resulting in loss of adiabaticity and vibrational memory before the reactants' collision. These theoretical findings require future confirmation through experimental or theoretical quantum mechanical studies, which are currently unavailable.

Theoretical study of the O(3P) + C2H6 reaction based on a new ab initio-based global potential energy surface

A new analytical potential energy surface was developed for the first time for the nine-body O(3P) + C2H6 hydrogen abstraction reaction, named PES-2020, which was fitted to explicitly-correlated high-level electronic structure calculations. This surface simulates the topography of the reactive system, from reactants to products, OH(v,j) + C2H5. The adiabatic energy of reaction, ΔHr(0 K) = -2.33 kcal mol-1, reproduces the experimental evidence, and the barrier height, 10.70 kcal mol-1, agrees with the ab initio calculations used as input. In addition, an intermediate complex in the exit channel is observed, which is stabilized with respect to the products of the reaction. Based on PES-2020 a dynamics study was carried out, where quasi-classical trajectory calculations were performed for collision energies in the range of 7.0-60.0 kcal mol-1, which covers high collision energy regions. The reaction cross section increases with collision energy; the largest fraction of available energy is deposited as translational energy (44-66%), and the scattering distribution evolves from backward to forward with collision energy. These findings reproduce previous theoretical calculations using electronic structure calculations of lower levels. However, where these previous studies failed, viz. in rotational and vibrational OH(v,j) distributions, PES-2020 reproduces practically quantitatively the experimental evidence, i.e., cold vibration and rotation, the rotational distribution peaking at j = 1-3 depending on the collision energy. In sum, this behaviour is typical of gas-phase hydrogen abstraction reactions with direct mechanism and high reaction barrier.

Quasi-Classical Trajectory Dynamics Study of the Cl(2P) + C2H6 → HCl(v,j) + C2H5 Reaction. Comparison with Experiment

To understand and simulate the dynamics behavior of the title reaction, QCT calculations were performed on a recently developed global analytical potential energy surface, PES-2017. These calculations combine the classical description of the dynamics with pseudoquantization in the reactants and products to perform a theoretical/experimental comparison on the same footing. Thus, in the products a series of constraints are included to analyze the HCl(v = 0,j) product, which is experimentally detected. At collision energies of 5.5 and 6.7 kcal mol^-1 the largest fraction of available energy is deposited as translation, 67%, while the ethyl radical shows significant internal energy, 27%, and so it does not act as a spectator of the reaction, thus reproducing recent experimental evidence. The HCl(v=0, j) rotational distribution is cold, peaking at j = 2, only one unit hotter than experiment, which represents an error of 0.12 kcal mol^-1. At a collision energy of 5.5 kcal mol^-1 product translational distribution is slightly hotter than experiment, but at 6.7 kcal mol^-1 agreement with recent experiments is practically quantitative, suggesting that the first experiments should be revised. In addition, we observe that the HCl(v=0, j) scattering distribution shifts from isotropic at low values of j to backward at high values of j, which is in agreement with experimental data. Finally, no evidence was found for the "chattering" mechanism suggested to explain the low translational energy of the HCl product in the backward scattering region. In sum, agreement with experiments of a series of sensible dynamic properties permits us to be optimistic on the quality and accuracy of the theoretical tools used in the present work, QCT and PES-2017.

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.