Classical versus quantum vibrational state distributions for the benchmark polyatomic reaction OH+D2: Checking the validity of the QCT method

First-principles mode-specific reaction dynamics

Controlling the outcome of chemical reactions by exciting specific vibrational and/or rotational modes of the reactants is one of the major goals of modern reaction dynamics studies. In the present Perspective, we focus on first-principles vibrational and rotational mode-specific dynamics computations on reactions of neutral and anionic systems beyond six atoms such as X + C2H6 [X = F, Cl, OH], HX + C2H5 [X = Br, I], OH- + CH3I, and F- + CH3CH2Cl. The dynamics simulations utilize high-level ab initio analytical potential energy surfaces and the quasi-classical trajectory method. Besides initial state specificity and the validity of the Polanyi rules, mode-specific vibrational-state assignment for polyatomic product species using normal-mode analysis and Gaussian binning is also discussed and compared with experiment.

The CN(X 2Σ+) + C2H6 reaction: Dynamics study based on an analytical full-dimensional potential energy surface

The hydrogen abstraction reaction of the cyano radical with molecules of ethane presents some interesting points in the chemistry from ultra-cold to combustion environments especially with regard to HCN(v) product vibrational distribution. In order to understand its dynamics, a new analytical full-dimensional potential energy surface was developed, named PES-2023. It uses a combination of valence bond and mechanic molecular terms as the functional form, fitted to high-level ab initio calculations at the explicitly correlated CCSD(T)-F12/aug-cc-pVTZ level on a reduced and selected number of points describing the reactive process. The new surface showed a continuous and smooth behavior, describing reasonably the topology of the reaction: high exothermicity, low barrier, and presence of intermediate complexes in the entrance and exit channels. Using quasi-classical trajectory calculations (QCT) on the new PES-2023, a dynamics study was performed at room temperature with special emphasis on the HCN(v1,v2,v3) product stretching and bending vibrational excitations, and the results were compared with the experimental evidence, which presented discrepancies in the bending excitation. The available energy was mostly deposited as HCN(v) vibrational energy with the vibrational population inverted in the CH stretching mode and not inverted in the CN stretching and bending modes, thus simulating the experimental evidence. Other dynamics properties at room temperature were also analyzed; cold rotational energy distribution was found, associated with a linear and soft transition state, and backward scattering distribution was found, associated with a rebound mechanism.

Dynamics and Novel Mechanisms of SN2 Reactions on ab Initio Analytical Potential Energy Surfaces

We describe a novel theoretical approach to the bimolecular nucleophilic substitution (S_N2) reactions that is based on analytical potential energy surfaces (PESs) obtained by fitting a few tens of thousands high-level ab initio energy points. These PESs allow computing millions of quasi-classical trajectories thereby providing unprecedented statistical accuracy for S_N2 reactions, as well as performing high-dimensional quantum dynamics computations. We developed full-dimensional ab initio PESs for the F^- + CH₃Y [Y = F, Cl, I] systems, which describe the direct and indirect, complex-forming Walden-inversion, the frontside attack, and the new double-inversion pathways as well as the proton-transfer channels. Reaction dynamics simulations on the new PESs revealed (a) a novel double-inversion S_N2 mechanism, (b) frontside complex formation,

Gaussian binning of the vibrational distributions for the Cl + CH4(v(4/2) = 0, 1) → H + CH3Cl(n(1)n(2)n(3)n(4)n(5)n(6)) reactions

We test several binning techniques to obtain mode-specific final-state distributions for polyatomic reactions. Normal mode analysis is done after an exact transformation to the Eckart frame. Standard histogram binning (HB) and three different variants of the energy-based Gaussian binning (1GB) are employed to obtain the probabilities of the vibrational states. We consider the two major issues of the polyatomic quasiclassical product analysis, i.e., (1) rounding the classical action to the nearest integer can result in unphysical states and (2) the normal-mode analysis can break down for highly distorted geometries. We show that 1GB can handle issue 1 when the total vibrational energy is evaluated in the normal mode space using the harmonic approximation and both issues 1 and 2 can be solved when the total vibrational energy is calculated exactly in the Cartesian space. We found that anharmonicity in the quantized energy levels does not have a significant effect on the final-state distributions. Quasiclassical trajectory calculations are performed for the reactant ground-state and bending-excited Cl((2)P(3/2)) + CH(4)(v(4/2) = 0, 1) → H + CH(3)Cl reactions using an ab initio potential energy surface. The product analysis techniques are successfully applied to the CH(3)Cl product molecules and some qualitative features of the results are discussed.