Global ab Initio Potential Energy Surfaces for Low-Lying Doublet States of NO3

Application of Reaction Path Search Calculations to Potential Energy Surface Fits

There has been significant progress in recent years in the use of machine learning techniques to model high-dimensional reactive potential energy surfaces using large-scale data obtained from ab initio electronic structure calculations. In these methods, the strategy used to gather data becomes a key issue as the molecular size increases. In this work, we examine the applicability of the reaction path search algorithm implemented in the Global Reaction Route Mapping (GRRM) code as a data-gathering approach. The electronic energies and gradients sampled by using the GRRM calculation are directly used in potential energy surface fitting to a permutationally invariant polynomial function. This simple approach was applied to the HNS and HCNO reaction systems, and we found that the fitted potential energy surfaces reasonably reproduce the features of the electronic structure calculations used in the GRRM calculations. This suggests that the GRRM sampling scheme can be used to construct an initial potential energy surface.

Roaming Reactions and Dynamics in the van der Waals Region

Roaming reactions were first clearly identified in photodissociation of formaldehyde 15 years ago, and roaming dynamics are now recognized as a universal aspect of chemical reactivity. These reactions typically involve frustrated near-dissociation of a quasibound system to radical fragments, followed by reorientation at long range and intramolecular abstraction. The consequences can be unexpected formation of molecular products, depletion of the radical pool in chemical systems, and formation of products with unusual internal state distributions. In this review, I examine some current aspects of roaming reactions with an emphasis on experimental results, focusing on possible quantum effects in roaming and roaming dynamics in bimolecular systems. These considerations lead to a more inclusive definition of roaming reactions as those for which key dynamics take place at long range.

A generalized unimolecular impulsive model for curved reaction path

This work aims to introduce a generalized impulsive model for unimolecular dissociation processes. This model allows us to take into account the curvature of the reaction path explicitly. It is a generalization of the previously developed multi-center impulsive model [P.-Y. Tsai and K.-C. Lin, J. Phys. Chem. A 119, 29 (2015)]. Several limitations of conventional impulsive models are eliminated by this study: (1) Unlike conventional impulsive models, in which a single molecular geometry is responsible for the impulse determination, the gradients on the whole dissociation path are taken into account. The model can treat dissociation pathways with large curvatures and loose saddle points. (2) The method can describe the vibrational excitation of polyatomic fragments due to the bond formation by multi-center impulse. (3) The available energy in conventional impulsive models is separated into uncoupled statistical and impulsive energy reservoirs, while the interplay between these reservoirs is allowed in the new model. (4) The quantum state correlation between fragments can be preserved in analysis. Dissociations of several molecular systems including the roaming pathways of formaldehyde, nitrate radical, acetaldehyde, and glyoxal are chosen as benchmarks. The predicted photofragment energy and vector distributions are consistent with the experimental results reported previously. In these examples, the capability of the new model to treat the curved dissociation path, loose saddle points, polyatomic fragments, and multiple-body dissociation is verified. As a cheaper computational tool with respect to ab initio on-the-fly direct dynamic simulations, this model can provide detailed information on the energy disposal, quantum state correlation, and stereodynamics in unimolecular dissociation processes.

Quasiclassical trajectory studies of 18O(3P) + NO2 isotope exchange and reaction to O2 + NO on D0 and D1 potentials

We report quasiclassical trajectory calculations for the bimolecular reaction (18)O((3)P) + NO2 on the recent potential energy surfaces of the ground (D0) and first excited (D1) states of NO3 [B. Fu, J. M. Bowman, H. Xiao, S. Maeda, and K. Morokuma, J. Chem. Theory. Comput. 9, 893 (2013)]. The branching ratio of isotope exchange versus O2 + NO formation, as well as the product angular distributions and energy and rovibrational state distributions are presented. The calculations are done at the collision energy of relevance to recent crossed beam experiments [K. A. Mar, A. L. Van Wyngarden, C.-W. Liang, Y. T. Lee, J. J. Lin, and K. A. Boering, J. Chem. Phys. 137, 044302 (2012)]. Very good agreement is achieved between the current calculations and these experiments for the branching ratio and final translational energy and angular distributions of isotope exchange products (16)O((3)P) + NO2 and O2 + NO formation products. The reactant (18)O atom results in (18)O(16)O but not N(18)O for the O2 + NO formation product channel, consistent with the experiment. In addition, the detailed vibrational and rotational state information of diatomic molecules calculated currently for the (34)O2 + NO formation channel on D0 and D1 states are in qualitative agreement with the previous experimental and theoretical results of the photodissociation of NO3 and are consistent with older thermal bimolecular kinetics measurements.