Variational transition state theory with least‐action tunneling calculations for the kinetic isotope effects in the Cl+H2 reaction: Tests of extended‐LEPS, information‐theoretic, and diatomics‐in‐molecules potential energy surfaces

Development of the Time-Independent Methods for the Cl + H2/F + HD Reaction Using Hyper-Spherical Coordinates Including (Full) Spin-Orbit Characteristics

Recently, a combined study of high-resolution molecular crossed beam experiment and accurate full-dimensional time-dependent theory, including full spin-orbit characteristics on the effect of electronic spin and orbital angular momenta in the F + HD reaction, was reported by some of us, focusing on the partial wave resonance phenomenon (Science 2021, 371, 936-940). It revealed that the time-dependent theory could explain all of the details observed in the high-resolution experiment. Here, we develop two time-independent close-coupling methods using hyperspherical coordinates, including the two-state model, where only a part of the spin-orbit characteristics is considered, and the six-state model, where the full spin-orbit characteristics is considered. With these two newly developed theoretical models and the adiabatic theoretical model, the detailed reaction dynamics of the F + HD (v = 0, j = 0) reaction and the Cl + H2 (v = 0, j = 0) reaction are investigated and compared. Some of the results are compared with the time-dependent quantum wave packet theory and the experimental observations, and good agreements have been obtained, which suggests the validity of the pure-procession approximation in the six-state model using different theoretical methods. This work demonstrates the ability of the reactive scattering theory including full spin-orbit characteristics for describing the reactions of a halogen atom plus hydrogen molecule and its isotopologues.

Significant Isotope Effects from the Nonadiabatic Couplings in the Cl(2P) + HD(v = 0, j = 0) Reaction

The impact of non-Born-Oppenheimer couplings on the isotopic effects in the reaction of the Cl(2P) atom with the HD (v = 0, j = 0) molecule is investigated with our recently developed nonadiabatic time-independent quantum scattering methods, where the full open-shell characteristics are included in the six-state model, and also with the recently developed two-state model solving by time-independent methods, where part of the open-shell characteristic is included. The same reaction is also calculated with the simple adiabatic model using the lowest adiabatic potential energy surface. Compared with the results from different models, it is found that the reactivity of the Cl + HD → HCl + D channel is significantly overestimated in the adiabatic model. In contrast, the reactivity of the other channel agrees well with the nonadiabatic models. This is due to the van der Waals well in the reactant channel being changed a lot by including the nonadiabatic couplings. These quantum dynamics calculations suggest that sometimes the adiabatic model should be used with caution; otherwise, it may result in significant deviations for some reactions.

Experimental and Theoretical Study of the Vibrationally Excited Reaction Cl + D2 (v = 1, j = 0) → DCl + D

Vibrationally excited reaction of Cl + D₂ (v = 1, j = 0) → DCl + D was investigated by a high-resolution crossed beam experiment, with D₂ molecules in the vibrationally excited state prepared by the scheme of stimulated Raman pumping. Differential cross sections (DCSs) were obtained at three collision energies of 4.03, 4.93, and 5.68 kcal/mol. Backward scattering is dominant for both DCl (v' = 0) and DCl (v' = 1) products, and no forward scattering signal was observed at these three collision energies. Collision-energy-dependent DCS in the backward scattering direction was measured at collision energies between 3.62 and 5.97 kcal/mol. Comparing with the DCSs from the vibrational ground state, it is found that the vibrational excitation of D₂ molecules significantly enhances the reactivity because of the later barrier nature of the reaction. No obvious oscillatory structure was found in the collision-energy-dependent DCS in the backward scattering direction, suggesting that the title reaction proceeds via a direct abstraction mechanism.

Degradation of Carbonyl Hydroperoxides in the Atmosphere and in Combustion

Oxygenates with carbonyl and hydroperoxy functional groups are important intermediates that are generated during the autoxidation of organic compounds in the atmosphere and during the autoignition of transport fuels. In the troposphere, the degradation of carbonyl hydroperoxides leads to low-vapor-pressure polyfunctional species that may precipitate in clouds and fog droplets or to the formation of secondary organic aerosols (SOAs). In combustion, the fate of carbonyl hydroperoxides is important for the performance of advanced combustion engines, especially for autoignition. A key fate of the carbonyl hydroperoxides is reaction with OH radicals, for which kinetics data are experimentally unavailable. Here, we study 4-hydroperoxy-2-pentanone (CH₃C(═O)CH₂CH(OOH)CH₃) as a model compound to clarify the kinetics of OH reactions with carbonyl hydroperoxides, in particular H atom abstraction and OH addition reactions. With a combination of electronic structure calculations, we determine previously missing thermochemical data, and with multipath variational transition state theory (MP-VTST), a multidimensional tunneling (MT) approximation, multiple-structure anharmonicity, and torsional potential anharmonicity, we obtained much more accurate rate constants than the ones that can computed by conventional single-structure harmonic transition state theory (TST) and than the empirically estimated rate constants that are currently used in atmospheric and combustion modeling. The roles of various factors in determining the rates are elucidated. The pressure-dependent rate constants for the addition reaction are computed using system-specific quantum RRK theory. The calculated temperature range is 298-2400 K, and the pressure range is 0.01-100 atm. The accurate thermodynamic and kinetics data determined in this work are indispensable in the global modeling of SOAs in atmospheric science and in the detailed understanding and prediction of ignition properties of hydrocarbons and alternative fuels.

Variational transition state theory: theoretical framework and recent developments

This article reviews the fundamentals of variational transition state theory (VTST), its recent theoretical development, and some modern applications.

A time-dependent wave packet study of the vibronic and spin-orbit interactions in the dynamics of Cl(2P)+H2→HCl(X̃ 1Σg+)+H(2S) reaction

We investigate the vibronic and spin-orbit (SO) coupling effects in the state-selected dynamics of the title reaction with the aid of a time-dependent wave packet approach. The ab initio potential energy surfaces of Capecchi and Werner [Science 296, 715 (2002)] have been employed for this purpose. Collinear approach of the Cl((2)P) atom to the H(2) molecule splits the degeneracy of the (2)P state and gives rise to (2)Sigma and (2)Pi electronic states. These two surfaces form a conical intersection at this geometry. These states transform as 1 (2)A('), 1 (2)A("), and 2 (2)A('), respectively, at the nonlinear configurations of the nuclei. In addition, the SO interaction due to Cl atom further splits these states into (2)Sigma(1/2), (2)Pi(3/2), and (2)Pi(1/2) components at the linear geometry. The ground-state reagent Cl((2)P(3/2))+H(2) correlates with (2)Sigma(1/2) and (2)Pi(3/2), where as the SO excited reagent Cl(*)((2)P(1/2))+H(2) correlates with (2)Pi(1/2) at the linear geometry. In order to elucidate the impact of the vibronic and SO coupling effects on the initial state-selected reactivity of these electronic states we carry out quantum scattering calculations based on a flux operator formalism and a time-dependent wave packet approach. In this work, total reaction probabilities and the time dependence of electronic population of the system by initiating the reaction on each of the above electronic states are presented. The role of conical intersection alone on the reaction dynamics is investigated with a coupled two-state model and for the total angular momentum J=0 (neglecting the electronic orbital angular momentum) both in a diabatic as well as in the adiabatic electronic representation. The SO interaction is then included and the dynamics is studied with a coupled three-state model comprising six diabatic surfaces for the total angular momentum J=0.5 neglecting the Coriolis Coupling terms of the Hamiltonian. Companion calculations are carried out for the uncoupled adiabatic and diabatic surfaces in order to explicitly reveal the impact of two different surface coupling mechanisms in the dynamics of this prototypical reaction.

Determining transition states from kinetic isotope effects

BEBOVA-based TS determination has been very successful in elucidating enzyme mechanisms at a level of detail that would be otherwise inaccessible. The resulting TS structures have been used successfully as the basis for designing TS mimics as enzyme inhibitors with dissociation constants to 10(-11) M. The structure interpolation approach has systematized the process of finding a TS, increasing both the speed and the accuracy of TS determination. The combination of information from several TSs into a unified model increases the accuracy of the process significantly and results in an extremely sensitive probe of changes in TS with varying reaction conditions (i.e., enzymatic vs nonenzymatic reactions, different enzymes, or different nucleophiles). The TS determination process is summarized in Fig. 15.

Kinetic Isotope Effects and Transition State Geometries. A Theoretical Investigation of E2 Model Systems

Ab initio calculations at the MP2/6-31+G level have been performed on E2 model systems to investigate whether differences in kinetic isotope effects correlate with changes in transition state geometries. By combining various nucleophiles (NH(2)(-), OH(-), F(-), PH(2)(-), SH(-), Cl(-)) and leaving groups (NH(3), Br(-), Cl(-), F(-), SH(-)) for reactions of the type Nu(-) + CH(3)CH(2)X, a large diversity of transition structures from reactant-like to product-like are generated. For each reaction one primary and two different alpha-secondary kinetic isotope effects are calculated. The primary kinetic isotope effects depend strongly on the nucleophilic placement in the periodic system, which mainly is due to differences in equilibrium isotope effects. When this effect is subtracted, the primary kinetic isotope effects display the expected maximum for symmetric transition structures, although the maximum is broad. The secondary kinetic isotope effects associated with the leaving group provide a qualitative correlation with the hybridization at the carbon, but the corresponding effects at the carbon where the hydrogen abstraction takes place is uncorrelated with the transition state geometry.