The sudden vector projection model for reactivity: mode specificity and bond selectivity made simple

Role of non-statistical effects in deciding the fate of HO3˙ in the atmosphere

HO3˙ has been postulated as a reservoir for OH˙ in various atmospheric reactions. Under collision-free conditions, experiments indicate that the lifetime of this species should be more than one microsecond. Interestingly, the binding energy of HO3˙ is estimated to be ∼3 kcal mol-1 by recent experimental as well as theoretical works. The value of the binding energy suggests that the lifetime of HO3˙ should be in the picosecond range, and with this lifetime, HO3˙ cannot act as a reservoir for OH˙. In the present work, using on-the-fly semiclassical dynamics, we argue that if non-RRKM effects are included, the lifetime of HO3˙ may be higher than that estimated from the binding energy.

Deciphering the Factors Controlling Hydrogen and Methyl Spillover upon Methane Dissociation on Rh/Cu(111) Single-Atom Alloy

Spillover of adsorbed species from one active site to another is a key step in heterogeneous catalysis. However, the factors controlling this step, particularly the spillover of polyatomic species, have rarely been studied. Herein, we investigate the spillover dynamics of H* and CH3* species on a single-atom alloy surface (Rh/Cu(111)) upon the dissociative chemisorption of methane (CH4), using molecular dynamics that considers both surface phonons and electron-hole pairs. These dynamical calculations are made possible by a high-dimensional potential energy surface machine learned from density functional theory data. Our results provide compelling evidence that the H* and CH3* can spill over on the metal surface at experimental temperatures and reveal novel dynamical features involving an internal motion during diffusion for CH3*. Increasing surface temperature has a minor effect on promoting spillover, as geminate recombinative desorption becomes more prevalent. However, the poisoning of the active site can be mitigated by the frequent gaseous molecular collisions that occur under ambient pressure in real-world catalysis, which transfer energy to the trapped adsorbates. Interestingly, the bulky CH3* exhibits a significant spillover advantage over the light H* due to its larger size, which facilitates energy acquisition. These insights help to advance our understanding of spillover in heterogeneous catalysis.

Intersystem Crossing Control of the Nb+ + CO2 → NbO+ + CO Reaction

The transfer of an oxygen atom from carbon dioxide (CO2) to a transition metal cation in the gas phase offers atomic level insights into single-atom catalysis for CO2 activation. Given that these reactions often involve open-shell transition metals, they may proceed through intersystem crossing between different spin manifolds. However, a definitive understanding of such spin-forbidden reaction requires dynamical calculations on multiple global potential energy surfaces (PESs) coupled by spin-orbit couplings. In this work, we report global PESs and spin-orbit couplings for three low-lying spin (quintet, triplet, and singlet) states for the reaction between the niobium cation (Nb+) and CO2, which are used to investigate the nonadiabatic reaction dynamics and kinetics. Comparison with experimental data of kinetics and collision dynamics shows satisfactory agreement. This reaction is found to be very similar to that between Ta+ + CO2. Specifically, our theoretical findings suggest that the rate-limiting step in this reaction is intersystem crossing, rather than potential barriers.

Mode-Specific Dynamics Studies for the Multichannel C2H2 + OH Reaction

The C2H2 + OH reaction is a key elementary reaction in acetylene oxidation, and the products forming in different reaction channels, such as C2H and CH3 radicals, are also important for subsequent reaction processes in the combustion process. In this work, we investigated the dynamics of the C2H2 + OH reaction with specific vibrational mode excitations and analyzed the mode specificity based on quasi-classical trajectory calculations on a recently developed full-dimensional potential energy surface. It is found that exciting OH stretching mode can promote the production of H + OCCH2 and CO + CH3, while the excitation of C-H symmetric/antisymmetric stretching mode of C2H2 can facilitate the H2O + C2H channel. Based on the prediction of vibrationally adiabatic and sudden vector projection models, the mode specificity in the C2H2 + OH reaction can be attributed to the difference in the degree of coupling between the initial motion mode and the reaction coordinate of each reaction path, which ultimately leads to the changes in rate constants and the product branching ratios. These findings can offer theoretical insights to regulate the branching ratio of the multichannel C2H2 + OH reaction.

Molecular Dynamics Simulations of Soft and Reactive Landing of Proteins Desorbed by Argon Cluster Bombardment

Reactive molecular dynamics (MD) simulations were conducted to investigate the soft and reactive landing of hyperthermal velocity proteins transferred to a vacuum using large argon clusters. Experimentally, the interaction of argon cluster ion beams (Ar1000-5000+) with a target biofilm was previously used in such a manner to transfer lysozymes onto a collector with the retention of their bioactivity, paving the way to a new solvent-free method for complex biosurface nanofabrication. However, the experiments did not give access to a microscopic view of the interactions needed for their full understanding, which can be provided by the MD model. Our reactive force field simulations clarify the landing mechanisms of the lysozymes and their fragments on collectors with different natures (gold- and hydrogen-terminated graphite). The results highlight the conditions of soft and reactive landing on rigid surfaces, the effects of the protein structure, energy, and incidence angle before landing, and the adhesion forces with the collector substrate. Many of the obtained results can be generalized to other soft and reactive landing approaches used for biomolecules such as electrospray ionization and matrix-assisted laser desorption ionization.

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.

Multistate Dynamics and Kinetics of CO2 Activation by Ta+ in the Gas Phase: Insights into Single-Atom Catalysis

The activation of carbon dioxide (CO2) by a transition-metal cation in the gas phase is a unique model system for understanding single-atom catalysis. The mechanism of such reactions is often attributed to a "two-state reactivity" model in which the high-energy barrier of a spin state correlating with ground-state reactants is avoided by intersystem crossing (ISC) to a different spin state with a lower barrier. However, such a "spin-forbidden" mechanism, along with the corresponding dynamics, has seldom been rigorously examined theoretically, due to the lack of global potential energy surfaces (PESs). In this work, we report full-dimensional PESs of the lowest-lying quintet, triplet, and singlet states of the TaCO2+ system, machine-learned from first-principles data. These PESs and the corresponding spin-orbit couplings enable us to provide an extensive theoretical characterization of the dynamics and kinetics of the reaction between the tantalum cation (Ta+) and CO2, which have recently been investigated experimentally at high collision energies using crossed beams and velocity map imaging, as well as at thermal energies using a selected-ion flow tube apparatus. The multistate quasi-classical trajectory simulations with surface hopping reproduce most of the measured product translational and angular distributions, shedding valuable light on the nonadiabatic reaction dynamics. The calculated rate coefficients from 200 to 600 K are also in good agreement with the latest experimental measurements. More importantly, these calculations revealed that the reaction is controlled by intersystem crossing, rather than potential barriers.

Molecular Dynamics Simulation of Thermal Nonequilibrium and Chemical Reaction Processes in Hydrogen Combustion

Using reactive force field (ReaxFF) and molecular dynamics simulation, we investigate the combustion process of hydrogen-oxygen systems in initial thermal nonequilibrium states with different translational and rovibrational temperatures for oxygen. The system studied in this work contains 300 oxygen molecules and 700 hydrogen molecules with a density of 7 times the air density. For this system, the characteristic relaxation times of oxygen and hydrogen vibrational energies are 0.173 and 0.249 ns, respectively. 0.6% of hydrogen undergoes a chemical reaction with oxygen during the thermal nonequilibrium relaxation stage. For the distribution of translational energy and vibrational energy of oxygen in the thermal nonequilibrium state, the maximum mean error of the statistical distribution in the simulation and the Boltzmann distribution at temperature calculated from the average kinetic energy of molecules is about 2.25 × 10-5. At the same time, it was observed in the simulation that many-body interactions play a certain role in the combustion process. Furthermore, we compare the ignition time and temperature rise behavior of different combustion mechanisms and molecular dynamics simulations starting from the thermal equilibrium state. These results will provide meaningful references for the construction of thermal nonequilibrium combustion chemical reaction mechanisms.

Ta+ and Nb+ + CO2: intersystem crossing in ion-molecule reactions

The reactions of Ta+ and Nb+ with CO2 proceed only by a highly efficient oxygen atom transfer reaction to the respective oxide at room temperature in the gas phase. Although the product spin states are not determined, thermochemistry dictates that they must be different from ground state quintet Ta+ and Nb+, implying that intersystem crossing (ISC) has occurred. Recent reactive scattering experiments found dominant indirect dynamics for the reaction with Ta+ hinting at a bottleneck along the reaction path. The question on the nature of the bottleneck, whether it involves a crossing point or a transition state, could not be finally answered because theory located both close to each other. Here, we aim at shedding further light onto the impact of intersystem crossing on the reaction dynamics and ultimately the reactivity of transition metal ion reactions in the gas phase. We employ a combination of thermal kinetics for Ta+ and Nb+ with CO2 using a selected-ion flow tube (SIFT) apparatus and differential scattering cross sections for Nb+ + CO2 from crossed-beam velocity map imaging. The reaction with niobium again shows dominant indirect dynamics and in general very similar dynamics compared to Ta+ + CO2. At thermal energies, both reactions show sub-collisional rate constants with small negative temperature dependencies. Experiments are complemented by high level quantum chemical calculations of the minimum energy pathway. Statistical modelling well-reproduces the experimental thermal rate constants, and suggests that the Nb+ reaction is rate-limited by the intersystem crossing at thermal energies.

New Full-Dimensional Reactive Potential Energy Surface for the H4 System

As the most abundant molecule in the universe, collisions involving H2 have important implications in astrochemistry. Collisions between hydrogen molecules also represent a prototype for assessing various dynamic methods for understanding fundamental few-body processes. In this work, we develop a new and highly accurate full-dimensional potential energy surface (PES) covering all reactive channels of the H2 + H2 system, which extends our previously reported H2 + H2 nonreactive PES [J. Chem. Theory Comput., 2021, 17, 6747] by adding 39,538 additional ab initio points calculated at the MRCI/AV5Z level in the reactive channels. The global PES is represented with high fidelity (RMSE = 0.6 meV for a total of 79,000 points) by a permutation invariant polynomial neural network (PIP-NN) and is suitable for studying collision-induced dissociation, single-exchange, as well as four-center exchange reactions. Preliminary quasi-classical trajectory studies on the new PIP-NN PES reveal strong vibrational enhancement of all reaction channels.

Solving the OH + Glyoxal Problem: A Complete Theoretical Description of Post-Transition-State Energy Deposition in Activated Systems

Activated chemistry in coupled reaction systems has broadened our understanding of the chemical kinetics. In the case of intermediates formed in gas phase abstraction reactions (e.g., OH + HC(O)C(O)H (glyoxal) →HC(O)CO + H2O), it is particularly crucial to understand how the reaction energy is partitioned between product species as this determines the propensity for a given product to undergo "prompt" dissociation (e.g., HC(O)CO → HCO + CO) before the excess reaction energy is removed. An example of such an activated system is the OH + glyoxal + O2 coupled reaction system. In this work, we develop a molecular dynamics pipeline, which, combined with a master equation analysis, accurately models previous experimental measurements. This new work resolves previous complexities and discrepancies from earlier master equation modeling for this reaction system. The detailed molecular dynamics approach employed here is a powerful new tool for modeling challenging activated reaction systems.

Quantum state-to-state nonadiabatic dynamics of the charge transfer reaction H+ + NO(X2Π) → H + NO+(X1Σ+): Influence of ro-vibrational excitation of NO

Quantum state-to-state nonadiabatic dynamics of the charge transfer reaction H+ + NO(X2Π, vi = 1, 3, ji = 0, 1) → H + NO+(X1Σ+) has been studied based on the recently constructed diabatic potential energy matrix. It was found that the vibrational excitation of reactant NO inhibits the reactivity, while the rotational excitation of reactant NO has little effect on the reaction probability. These attributes were also observed in the semi-classical trajectory calculations employed in the adiabatic representation. Such an inhibitory effect of the vibrational excitation of reactant NO was owing to lower accessibility of the conical intersection and avoided crossing regions, which are located in the wells with respect to the Π diabat, as evidenced by the analysis of the population of the time-independent wave functions. Calculated vibrational state distributions of the product show that the decrease of the reaction mainly leads to the less formation of low vibrational states (vf < 6), and the product vibrational state distributions are more evenly populated for vi = 1 and 3, suggesting a non-statistical behavior. However, the overall shapes of the product rotational distributions remain unchanged, indicating that the redistribution of energy into the rotation of product NO is sufficient in the charge transfer process between H+ and NO. While the reaction is dominated by the forward and backward scattering in differential cross sections (DCSs), consistent with the complex-forming mechanism, a clear forward bias in the DCSs appears, indicating that the occurrence of the reaction is not sufficiently long to undergo the whole phase space of the interaction configurations.

High vibrational excitation of the reagent transforms the late-barrier H + HOD reaction into an early-barrier reaction

Polanyi's rules predict that a late-barrier reaction yields vibrationally cold products; however, experimental studies showed that the H2 product from the late-barrier H + H2O(|04⟩-) and H + HOD(vOH = 4) reactions is vibrationally hot. Here, we report a potential-averaged five-dimensional state-to-state quantum dynamics study for the H + HOD(vOH = 0-4) → H2 + OD reactions on a highly accurate potential energy surface with the total angular momentum J = 0. It is found that with the HOD vibration excitation increasing from vOH = 1 to 4, the product H2 becomes increasingly vibrationally excited and manifests a typical characteristic of an early barrier reaction for vOH = 3 to 4. Analysis of the scattering wave functions revealed that vibrational excitation in the breaking OH bond moves the location of dynamical saddle point from product side to reactant side, transforming the reaction into an early barrier reaction. Interestingly, pronounced oscillatory structures in the total and product vibrational-state-resolved reaction probabilities were observed for the H + HOD(vOH = 3, 4) reactions, in particular at low collision energies, which originate from the Feshbach resonance states trapped in the bending/torsion excited vibrational adiabatic potential wells in the entrance region due to van der Waals interactions.

Combined Quantum Mechanical and Quasi-Classical State-to-State Dynamical Study on the Isotopic Effect in H/D + LiH+/LiD+ → H2/HD/D2 + Li+ Reactions

Coriolis-coupled quantum mechanical (QM-CC) and quasi-classical trajectory (QCT) calculations are carried out to investigate the dynamics of the H(D) + LiH+(v = 0, j = 0) → H2(HD) (v', j') + Li+ reactions on the ground electronic state potential energy surface reported by Martinazzo et al. (Martinazzo et al., J. Chem. Phys. 2003, 119, 11241). The QM-CC and QCT results at the initial state-selected and state-to-state levels are used to investigate the validity and accuracy of the QCT method for these exoergic barrierless reactions. Furthermore, the QCT method is used to understand the isotopic effects on reaction observables like total and state-to-state integral cross section, differential cross section, product energy disposal, and rate constants of H(D) + LiH+(v = 0, j = 0) → H2(HD) (v', j') + Li+ and H(D) + LiD+(v = 0, j = 0) → HD(D2) (v', j') + Li+ reactions. Attempts are also made to understand the impact of the isotopic substitution on the reaction mechanism. It is observed that QM-CC and QCT results closely follow each other at the initial state-selected and state-to-state levels. Noticeable kinematic effects of reagents on the reactivity and mechanism of the reactions are also observed.

State-to-State Dynamics in Mode-Selective Polyatomic Reactions

Chemical reactions are intrinsically quantum mechanical transformations of reactants to products. Recent experimental and theoretical advances have enabled the exploration of reaction dynamics with a quantum state resolution for both reactants and products. To this end, reactions involving more than three atoms are of particular interest, because they exhibit rich dynamics concerning the role of different reactant modes in controlling reactivity and product energy disposal. A clear understanding of the state-to-state dynamics requires new paradigms. In this Perspective, we examine some new concepts that have emerged from recent state-to-state studies of polyatomic reactions and illustrate the key role played by the transition state.

Differential Cross Sections for the H + H2S → H2 + SH Reaction: A Full-Dimensional State-to-State Quantum Dynamics Study

We utilized the time-dependent wave packet approach to compute the first full-dimensional (6D) state-to-state differential cross sections (DCSs) for the title reaction with the initial nonrotating H2S in the ground and the (100) and (001) vibrational excited states. It is found that the fundamental symmetric and asymmetric stretching excitations of H2S exhibit almost the same influence on the DCS, but unlike the H + H2O → H2 + OH reaction, they greatly increase the vibrational excitation of both the H2 and SH products. The hot vibrational state distributions of H2 are consistent with the prediction of product energy disposal by Polanyi's rules for an early barrier reaction. Because the incident H atom reacts strongly with both the ground and excited S-H states, the large populations of product SH(v2 = 1), which are very close to the relative reactivity of the initial S-H(v = 0) state, can still be explained by the local mode picture for H2S and the nonreacting SH bond's spectator nature.

Formation of N(2D) from Hyperthermal Collisions between O(3P) and NO(X2Π)

Hyperthermal collisions between O(3P) and NO(X2Π) could lead to the formation of the first electronically excited atomic nitrogen (N(2D)), which plays a key role in plasma formation in shock-heated air. This process is facilitated mainly by four doublet states, and to a much lesser extent by two quartet states. In this work, we report quasi-classical trajectory studies of this reactive process using the four analytical adiabatic potential energy surfaces for the doublet states developed previously from fitting high-level ab initio data. The reactions were found to be strongly enhanced by vibrational excitation of the NO reactant, consistent with the existence of potential energy barriers in the exit channel. Despite the large endothermicity of the reaction, the rate coefficient is significant at high temperatures, suggesting a possible role of this reaction in the hyperthermal kinetics in the shock layer of a hypersonic (re)entry vehicle.

Full-dimensional automated potential energy surface development and detailed dynamics for the CH2OO + NH3 reaction

With the help of the ROBOSURFER program package, a global full-dimensional potential energy surface (PES) for the reaction of the Criegee intermediate, CH2OO, with the NH3 molecule is developed iteratively using different ab initio methods and the monomial symmetrization fitting approach. The final permutationally-invariant analytical PES is constructed based on 23447 geometries and the corresponding ManyHF-based CCSD(T)-F12b/cc-pVTZ-F12 energies. The accuracy of the PES is confirmed by the excellent agreement of its stationary-point properties and one-dimensional potential energy curves compared with the corresponding ab initio data. The reaction probabilities and integral cross sections are calculated for the ground-state and several vibrationally excited-state reactions by quasi-classical trajectory simulations. Remarkable is that the maximum impact parameter b where reactivity vanishes is almost independent of collision energy ranging from 1 to 40 kcal mol-1, and the reaction probability increases with increasing collision energy for this negative-barrier reaction. At the same time, a slight mode-specificity effect is observed. In addition, the deuterium effect is investigated and the sudden vector projection is discussed.

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.

Vibrational Mode-Specific Dynamics of the OH + C2H6 Reaction

We investigate the effects of the initial vibrational excitations on the dynamics of the OH + C2H6 → H2O + C2H5 reaction using the quasi-classical trajectory method and a full-dimensional analytical ab initio potential energy surface. Excitation of the initial CH, CC, and OH stretching modes enhances, slightly inhibits, and does not affect the reactivity, respectively. Translational energy activates the early-barrier title reaction more efficiently than OH and CC stretching excitations, in accord with the Polanyi rules whereas CH stretching modes have similar or higher efficacy than translation, showing that these rules are not always valid in polyatomic processes. Scattering angle, initial attack angle, and product translational energy distributions show the dominance of direct stripping with increasing collision energy, side-on OH and isotropic C2H6 attack preferences, and substantial reactant-product translational energy transfer without any significant mode specificity. The reactant vibrational excitation energy of OH and C2H6 flows into the H2O and C2H5 product vibrations, respectively, whereas product rotations are not affected. The computed mode-specific H2O vibrational distributions show that initial OH excitation appears in the asymmetric stretching vibration of the H2O product and allow comparison with experiments.

Theoretical studies on the kinetics and dynamics of the BeH+ + H2O reaction: comparison with the experiment

The reaction of BeH+ with background gaseous H2O may play a role in qubit loss for quantum information processing with Be+ as trapped ions, and yet its reaction mechanism has not been well understood until now. In this work, a globally accurate, full-dimensional ground-state potential energy surface (PES) for the BeH+ + H2O reaction was constructed by fitting a total of 170 438 ab initio energy points at the level of RCCSD(T)-F12/aug-cc-pVTZ using the fundamental invariant-neural network method. The total root-mean-square error of the final PES was 0.178 kcal mol-1. For comparison, quasi-classical trajectory calculations were carried out on the PES at an experimental temperature of 150 K. The obtained thermal rate constant and product branching ratio of the BeD+ + H2O reaction agreed quite well with experimental results. In addition, the vibrational state distributions and energy disposals of the products were calculated and rationalized using the sudden vector projection model.

Enhancing reactivity of SiO+ ions by controlled excitation to extreme rotational states

Optical pumping of molecules provides unique opportunities for control of chemical reactions at a wide range of rotational energies. This work reports a chemical reaction with extreme rotational excitation of a reactant and its kinetic characterization. We investigate the chemical reactivity for the hydrogen abstraction reaction SiO⁺ + H₂ → SiOH⁺ + H in an ion trap. The SiO⁺ cations are prepared in a narrow rotational state distribution, including super-rotor states with rotational quantum number (j) as high as 170, using a broad-band optical pumping method. We show that the super-rotor states of SiO⁺ substantially enhance the reaction rate, a trend reproduced by complementary theoretical studies. We reveal the mechanism for the rotational enhancement of the reactivity to be a strong coupling of the SiO⁺ rotational mode with the reaction coordinate at the transition state on the dominant dynamical pathway.

State-to-State Dynamics in Mode-Specific Reactions of Cl + CH3D(v1-I, v1-II, and v4 = 1; |10⟩): Loss of Memory or Not

Several decades of the study of reaction dynamics culminate in the concept of mode specificity and bond selectivity in polyatomic systems. Until very recently, the main concern of those studies has been total reactivity and little attention has been paid to the mode-specific effects on the more detailed product-state and angular distributions. Conventional wisdom would anticipate that the fine detail should reveal a more pronounced mode dependency. However, a few recent studies showed that the product distributions could appear to be surprisingly insensitive to the modes of internal excitation of reagents. This counterintuitive finding led to a concept of loss of memory. Here, we present detailed experimental results in the reactions of the Cl atom with three distinct stretching-excited CH₃D(v_CH3 = 1) reagents. In conjunction with the previous reports on various aspects of this reaction, such a comprehensive set of data enables us to perform an in-depth examination of the validity of this new concept.

Stereodynamical control of the H + HD → H2 + D reaction through HD reagent alignment

Prealigning nonpolar reacting molecules leads to large stereodynamical effects because of their weak steering interaction en route to the reaction barrier. However, experimental limitations in preparing aligned molecules efficiently have hindered the investigation of steric effects in bimolecular reactions involving hydrogen. Here, we report a high-resolution crossed-beam study of the reaction H + HD(v = 1, j = 2) → H₂(v', j') + D at collision energies of 0.50, 1.20, and 2.07 electron volts in which the vibrationally excited hydrogen deuteride (HD) molecules were prepared in two collision configurations, with their bond preferentially aligned parallel and perpendicular to the relative velocity of collision partners. Notable stereodynamical effects in differential cross sections were observed. Quantum dynamics calculations revealed that strong constructive interference in the perpendicular configuration plays an important role in the stereodynamical effects observed.

Imaging the Complex-Forming Reaction Dynamics in Al + CO2 → AlO + CO

For indirect reactions involving more than one intermediate complex from reactant valley to product valley, the reaction dynamics becomes very complicated for researchers due to competition between pathways. In order to resolve the large discrepancy between theoretical and experimental studies on the linear or bent structures of complexes involved in the title endothermic reaction, we performed a crossed-beam experiment at a large collision energy (E_c) range with mapping of the velocity distributions of Al reactants and AlO products. We found that the reaction proceeds through different complex-forming mechanisms with the increase of E_c: at low E_c near the reaction threshold, only low impact-parameter collisions contribute through a collinear Al-OCO short-lived complex, and at high E_c, the bent-structure complexes, formed by either isomerization or noncollinear collisions, come into play.

Dynamically Hidden Reaction Paths in the Reaction of CF3+ + CO

Reaction paths on a potential energy surface are widely used in quantum chemical studies of chemical reactions. The recently developed global reaction route mapping (GRRM) strategy automatically constructs a reaction route map, which provides a complete picture of the reaction. Here, we thoroughly investigate the correspondence between the reaction route map and the actual chemical reaction dynamics for the CF₃⁺ + CO reaction studied by guided ion beam tandem mass spectrometry (GIBMS). In our experiments, FCO⁺, CF₂⁺, and CF⁺ product ions were observed, whereas if the collision partner is N₂, only CF₂⁺ is observed. Interestingly, for reaction with CO, GRRM-predicted reaction paths leading to the CF⁺ + F₂CO product channel are found at a barrier height of about 2.5 eV, whereas the experimentally obtained threshold for CF⁺ formation was 7.48 ± 0.15 eV. In other words, the ion was not obviously observed in the GIBMS experiment, unless a much higher collision energy than the requisite energy threshold was provided. On-the-fly molecular dynamics simulations revealed a mechanism that hides these reaction paths, in which a non-statistical energy distribution at the first collisionally reached transition state prevents the reaction from proceeding along some reaction paths. Our results highlight the existence of dynamically hidden reaction paths that may be inaccessible in experiments at specific energies and hence the importance of reaction dynamics in controlling the destinations of chemical reactions.

Stretching-mode specificity in the Cl + CH3D(v1-I, v1-II, and v4 = 1; |jK〉) reactions: dependency on the initial |jK〉 selectivity

The title reactions were studied at a collisional energy of 5.4 kcal mol^-1 in a crossed-beam product-imaging experiment. The dynamics attributes of the dominant ground-state CH₂D(0₀) and the accompanied C-D bend-excited CH₂D(6₁) products were imaged in reactions with totally 16 ro-vibrationally selected states of the CH₃D(v_i, |jK〉) reagents. We found that all three vibrational excitations yielded marked |jK〉-dependent rate-enhancements in forming the (0₀, 0/1)_s product pairs. Furthermore, for a given rotational |jK〉-mode, a vibrational-mode propensity of v₄ > v₁-I > v₁-II in rate promotion and a clear manifestation of the Fermi-phase-induced interference effect of the latter two were observed. Compared to the reactivity of the rotationless state |jK〉 = |00〉, a minute rotational-excitation of all three stretch-excited CH₃D(v_i = 1) reagents could yield significantly higher reaction rates for the product pair (0₀, 0)_s, but not so for (0₀, 1)_s. The signals in forming the (6₁, 0)_s pair were clearly notable but smaller than that of the ground-state reaction product pair, (0₀, 0)_g. An opposite propensity of v₁-II ≈ v₁-I > v₄ with a milder dependency on the initial |jK〉-states was observed. The angular distributions of the (0₀, 0)_s pairs were nearly identical for all ro-vibrationally excited reagents, displaying the typical trait for a direct abstraction of the rebound mechanism. Similar distributions were found for the (6₁, 0)_s pairs; yet, both pairs deviated substantially from the peripheral feature of the ground-state reaction pair of (0₀, 0)_g. Those of the (0₀, 1)_s pairs in reactions with v₄-excitation featured a prominent forward-peaking distribution-suggestive of a time-delayed, resonance-mediated pathway, again with little dependency on the initial |jK〉-states. As for the reactions with the two Fermi-dyads, v₁-I and v₁-II, albeit showing globally similar distributions to that for v₄, a substantial variation with the initial rotational-mode excitation could be discerned in the forward-peaking features. To unravel the impact of the Fermi-phase on the |jK〉-dependent attributes, we adopted a comparative approach by contrasting the observations in reactions with the Fermi-dyad reagents (the superposition states) to those with the pure-state reagents. Remarkable distinctions are unveiled and elucidated with the unexplained results explicitly pointed out, which call for future theoretical investigations for deeper understanding.

Nonadiabatic Reactive Quenching of OH(A2Σ+) by H2: Origin of High Vibrational Excitation in the H2O Product

The nonadiabatic dynamics of the reactive quenching channel of the OH(A²Σ⁺) + H₂/D₂ collisions is investigated with a semiclassical surface hopping method, using a recently developed four-state diabatic potential energy matrix (DPEM). In agreement with experimental observations, the H₂O/HOD products are found to have significant vibrational excitation. Using a Gaussian binning method, the H₂O vibrational state distribution is determined. The preferential energy disposal into the product vibrational modes is rationalized by an extended Sudden Vector Projection model, in which the h and g vectors associated with the conical intersection are found to have large projections with the product normal modes. However, our calculations did not find significant insertion trajectories, suggesting the need for further improvement of the DPEM.

Theoretical dynamics studies of the CH3 + HBr → CH4 + Br reaction: integral cross sections, rate constants and microscopic mechanism

Quantum and quasi-classical dynamics calculations have been performed for the reaction of HBr with CH₃. The accurate ab initio-based potential energy surface function developed earlier for this reaction displays a potential well corresponding to a reactant complex and a submerged potential barrier. The integral cross sections were calculated on this potential energy surface using both a six-degree-of-freedom reduced dimensional quantum dynamics and the quasi-classical trajectory method and very good agreement was found between the two approaches. The cross sections were found to diverge when the collision energy decreases, indicating that the reactant attraction is responsible for the dynamics at low collision energy. The quantum mechanical and the quasi-classical rate constants also agree very well and almost exactly reproduce the experimental results at low temperatures up to 540 K. The negative activation energy observed experimentally is confirmed by the calculations and is a consequence of the long-range attraction between the reactants. From the classical trajectories mechanistic details have been extracted. It is found that at very low collision energy, the reacting system crosses the potential barrier because the forces within the complex guide them, although some 30% is reflected from the product side of the barrier. When the collision energy increases, the system does not follow the most favorable path and the reactants are, with increasing probability, reflected from the repulsive walls of the nonreactive parts of the reactants, providing a picture beyond the decreasing excitation function.

Fermi-phase-induced interference in the reaction between Cl and vibrationally excited CH3D

Mode selectivity is a well-established concept in chemical dynamics. A polyatomic molecule possesses multiple vibrational modes and the mechanical couplings between them can result in complicated anharmonic motions that defy a simple oscillatory description. A prototypical example of this is Fermi-coupled vibration, in which an energy-split eigenstate executes coherent nuclear motion that is comprised of the constituent normal modes with distinctive phases. Will this vibrational phase affect chemical reactivity? How can this phase effect be disentangled from more classical amplitude effects? Here, to address these questions, we study the reaction of Cl with a pair of Fermi states of CH₃D(v₁-I and v₁-II). We find that the reactivity ratio of (v₁-I)/(v₁-II) in forming the CH₂D(v = 0) + HCl(v) products deviates significantly from that permitted by the conventional reactivity-borrowing framework. Based on a proposed metric, this discrepancy can only be explained when the scattering interferences mediated by the CH₃D vibrational phases are explicitly considered, which expands the concept of vibrational control of chemical reactivity into the quantum regime.

Reaction Pathway Control via Reactant Vibrational Excitation and Impact on Product Vibrational Distributions: The O + HO2 → OH + O2 Atmospheric Reaction

Chemical reactions often have multiple pathways, the control of which is of fundamental and practical importance. In this Letter, we examine the dynamics of the O + HO₂ → OH + O₂ reaction, which plays an important role in atmospheric chemistry, using quasi-classical trajectories on a recently developed full-dimensional potential energy surface (PES). This reaction has two pathways leading to the same products: the H abstraction pathway (O_a + HO_bO_c → O_aH + O_bO_c) and the O abstraction pathway (O_a + HO_bO_c → O_bH + O_aO_c). Under thermal conditions, the reaction is dominated by the latter channel, which is barrierless, leading to vibrational excitation of the O₂ product. However, we demonstrate that excitation of the HO₂ reactant in its O-H (v₁) vibrational mode results in dramatic switching of the reaction pathway to the activated H abstraction channel, which leads to a highly excited OH product vibrational state distribution. The implications of such dynamical effects in the atmospheric chemistry are discussed.

Toward a Comprehensive Understanding of Mode-Specific Dynamics of Polyatomic Reactions: A Full-Dimensional Quantum Dynamics Study of the H + NH3 Reaction

Mode specificity not only sheds light on reaction dynamics but also opens the door for chemical reaction control. This work reports a state-of-the-art full-dimensional quantum dynamics study on the prototypical hydrogen abstraction reaction of hydrogen with ammonia, which serves as a benchmark for advancing our fundamental understanding of polyatomic reaction dynamics. By taking advantage of the (3 + 1) Radau-Jacobi coordinates, the bond-specific probabilities are resolved with the reactant NH₃ initiated from either a non-degenerate or degenerate stretching vibrational state. The observed different atom-specific abstraction probabilities from individual states of the degenerate pair are rationalized in the local mode representation according to the different vibrational energy deposited in each N-H bond. It is verified that the three H atoms in NH₃ have the same abstraction probability as that from the degenerate pair and the linear combination of the degenerate pair gives the same reaction probability. In addition, the symmetric and asymmetric stretching modes of the reactant NH₃ have comparable efficacies on driving the reaction.

Mode Selective Chemistry for the Dissociation of Methane on Efficient Ni/Pt-Bimetallic Alloy Catalysts

The mode selectivity of methane dissociation is studied on three different Ni/Pt-bimetallic alloy surfaces using a fully quantum approach based on reaction path Hamiltonian. Dissociative sticking probability depends on the...

Spiers Memorial Lecture: theory of unimolecular reactions

One hundred years ago, at an earlier Faraday Discussion meeting, Lindemann presented a mechanism that provides the foundation for contemplating the pressure dependence of unimolecular reactions. Since that time, our...

Time-dependent quantum dynamics study of the F + C2H6 → HF + C2H5 reaction

The reaction probabilities, integral cross sections, energy efficiency and rate constants are investigated for the F + C₂H₆ reaction using the quantum reaction dynamics, wave packet method. The ground-state integral cross section calculated using a six-degree-of-freedom approach is in very good agreement with the quasi-classical trajectory results. We find that the H-CH₂CH₃ stretching motion has the largest enhancement to reactivity, followed by the H-CH₂-CH₃ bending motion. However, the stretching motion between CH₂ and CH₃ slightly hinders the reactivity. The energy-form efficacy based on an equal amount of total energy shows that translational energy is more effective in enhancing the reactivity than vibrational energy of the H-CH₂CH₃ stretching motion at a relatively lower translational energy, while the reverse is true at a relatively high translational energy. An energy-shifting method is employed to calculate the full-dimensional rate constants. The quantum rate constants agree well with one of the two main experimental measurements, and the activation energy has an excellent agreement with the one calculated using canonical variational transition-state theory.

Vibrational mode-specific dynamics of the F(2P3/2) + C2H6 → HF + C2H5 reaction

We investigate the competing effect of vibrational and translational excitation and the validity of the Polanyi rules in the early- and negative-barrier F(²P_3/2) + C₂H₆ → HF + C₂H₅ reaction by performing quasi-classical dynamics simulations on a recently developed full-dimensional multi-reference analytical potential energy surface. The effect of five normal-mode excitations of ethane on the reactivity, the mechanism, and the post-reaction energy flow is followed through a wide range of collision energies. Promoting effects of vibrational excitations and interaction time, related to the slightly submerged barrier, are found to be suppressed by the early-barrier-induced translational enhancement, in contrast to the slightly late-barrier Cl + C₂H₆ reaction. The excess vibrational energy mostly converts into ethyl internal excitation while collision energy is transformed into product separation. The substantial reaction energy excites the HF vibration, which tends to show mode-specificity and translational energy dependence as well. With increasing collision energy, direct stripping becomes dominant over the direct rebound and indirect mechanisms, being basically independent of reactant excitation.

Dissociative photodetachment of H3O2-: a full-dimensional quantum dynamics study

The transition state is a central concept of chemistry. Photoelectron-photofragment coincidence (PPC) spectroscopy has been proven as an attractive method to study the transition state dynamics. Within a state-of-the-art full-dimensional quantum mechanical model, the dissociative photodetachment dynamics of H₃O₂^- is investigated on accurate anion and neutral potential energy surfaces. The calculated PPC spectrum of H₃O₂^- agrees well with the experimental measurement. The dissociative product OH is exclusively populated on the ground vibrational state, implying the character of the spectator bond. In contrast, the product H₂O is predominantly populated in the ground and fundamental states of the symmetric and antisymmetric stretching modes, which is caused by the strong coupling between the antisymmetric motion of the transferred H atom in the transient intermediate [H₃O₂]* and both stretching modes of the product H₂O.

The Role of Zero-Point Vibration and Reactant Attraction in Exothermic Bimolecular Reactions with Submerged Potential Barriers: Theoretical Studies of the R + HBr → RH + Br (R = CH3, HO) Systems

The dynamics of the reactions CH₃ + HBr → CH₄ + Br and HO + HBr → H₂O + Br have been studied using the quasiclassical trajectory method to explore the interplay of the vibrational excitation of the breaking bond and the potential energy surface characterized by a prereaction van der Waals well and a submerged barrier to reaction. The attraction between the reactants is favorable for the reaction, because it brings together the reactants without any energy investment. The reaction can be thought to be controlled by capture. The trajectory calculations indeed provide excitation functions typical to capture: the reaction cross sections diverge when the collision energy is reduced toward zero. Excitation of reactant vibration accelerates both reactions. The barrier on the potential surface is so early that the coupling between the degrees of freedom at the saddle point geometry is negligible. However, the trajectory calculations show that when the breaking bond is stretched at the time of the encounter, an attractive force arises, as if the radical approached a HBr molecule whose bond is partially broken. As a result, the dynamics of the reaction are controlled more by the temporary "dynamical", vibrationally induced than by the "static" van der Waals attraction even when the reactants are in vibrational ground state. The cross sections are shown to drop to very small values when the amplitude of the breaking bond's vibration is artificially reduced, which provides an estimate of the reactivity due to the "static" attraction. Without zero-point vibration these reactions would be very slow, which is a manifestation of a unique quantum effect. Reactions where the reactivity is determined by dynamical factors such as the vibrationally enhanced attraction are found to be beyond the range of applicability of Polanyi's rules.

Mode-specific dynamics in multichannel reaction NH+ + H2

The vibrational- and rotational-mode specificity of the multichannel NH⁺ + H₂ reaction is studied on a recently constructed ab initio-based global potential energy surface using an initial state selected quasi-classical trajectory method, and the trajectories are analyzed using an isometric feature mapping and k-means approach. All excitation modes promote two reactions (R1: NH'⁺ + H₂ → NH⁺ + HH' and R4: NH'⁺ + H₂ → NH₂⁺ + H') where both NH and HH bonds are broken, but reduce the reactivity of the proton-transfer reaction R2 (NH'⁺ + H₂ → N + H'H₂⁺) at low collision energies. For the hydrogen-transfer reaction R3 (NH'⁺ + H₂ → HNH'⁺ + H), the rotational excitation of NH⁺ enhances the reactivity remarkably, while its vibrational excitation has an inhibiting effect on the reaction. The trajectory analyses show that the vibrational and rotational excitations of NH⁺ make R3 tend to go over a submerged saddle point instead of extracting hydrogen atoms directly. On the other hand, the motions of the H₂ reactant facilitate the enhancement of the reactivity but they do not affect the mechanism of R3. In addition, the results suggest that the coupling of the isometric feature mapping and the k-means approach in the trajectory analysis is an appropriate tool for reaction-dynamics studies.

Vibrational mode-specificity in the dynamics of the Cl + C2H6 → HCl + C2H5 reaction

We report a detailed dynamics study on the mode-specificity of the Cl + C₂H₆ → HCl + C₂H₅ H-abstraction reaction. We perform quasi-classical trajectory simulations using a recently developed high-level ab initio full-dimensional potential energy surface by exciting five different vibrational modes of ethane at four collision energies. We find that all the studied vibrational excitations, except that of the CC-stretching mode, clearly promote the title reaction, and the vibrational enhancements are consistent with the predictions of the Sudden Vector Projection (SVP) model, with the largest effect caused by the CH-stretching excitations. Intramolecular vibrational redistribution is also monitored for the differently excited ethane molecule. Our results indicate that the mechanism of the reaction changes with increasing collision energy, with no mode-specificity at high energies. The initial translational energy mostly converts into product recoil, while a significant part of the excess vibrational energy remains in the ethyl radical. An interesting competition between translational and vibrational energies is observed for the HCl vibrational distribution: the effect of exciting the low-frequency ethane modes, having small SVP values, is suppressed by translational excitation, whereas a part of the excess vibrational energy pumped into the CH-stretching modes (larger SVP values) efficiently flows into the HCl vibration.

Kinetic and Dynamic Studies of the F(2P) + ND3 → DF + ND2 Reaction

The fast F reaction with NH₃ poses a big challenge to experimental studies because of secondary chemical and collisional reactions. The quasi-classical trajectory method is utilized to investigate the mode specificity, product energy disposal, and temperature dependence of the thermal rate coefficient of F + ND₃ → DF + ND₂ on a recently developed potential energy surface. The effect of isotopic substitution is explored by comparing the F + ND₃ reaction with the F + NH₃ reaction. The computed results permit a better understanding of the F + ammonia reaction. The DF vibrational state has a Λ-type distribution, in accordance with the experimental measurement by the fast flow reactor technique. The product ND₂ is dominantly populated in the ground state, and a considerable amount of ND₂ is produced in the fundamental states of the bending mode. The similar vibrational state distributions of HF and NH₂ in the F + NH₃ reaction indicate a weak isotopic substitution effect on the product energy disposal. Exciting the umbrella mode of ND₃ suppresses the reaction at low energies below 5 kcal mol^-1, in sharp contrast to the observation in the F + NH₃ reaction. These dynamical behaviors can be partially explained by the sudden vector projection model. In addition, the thermal rate coefficient of F + ND₃ shows no temperature dependence in the range between 150 and 2000 K. There exists an inverse kinetic isotope effect at temperatures from 150 to 1500 K.

Pair-Correlated Imaging of Cl + CH3D(v4, v1-I, v1-II = 1, |jK⟩) → CH2D(vi) + HCl(v)

The title reactions were studied at a collisional energy of 10.0 kcal mol^-1 in a crossed-beam, product-imaging experiment. In terms of integral cross sections, all three CH₃-stretching excited CH₃D(v_CH3 = 1) reagents promote the reactivity in forming the predominant product pair of (v_CH2D, v_HCl)_s = (0₀, 0/1)_s with a prominent mode-propensity of v₄ > v₁-I > v₁-II, where v₄ denotes the degenerate mode of CH₃ asymmetric stretch and v₁-I and v₁-II are a pair of Fermi-coupled, symmetric-stretch states. The vibrationally excited CH₂D product pairs of (6₁, 0)_s, (1₁, 0)_s, and (3₁, 0)_s appear to be minor channels and display a reverse propensity of v₄ < v₁-I ≈ v₁-II for (6₁, 0)_s, while v₄ > v₁-I for (1₁, 0)_s. Based on the observed angular distributions, we conjecture that, irrespective of the initial mode of excitation, the (0₀, 0)_s product pair proceeds by a direct abstraction of the peripheral type, whereas the (0₀,1)_s pair is mediated via a resonance pathway. Intriguingly, the angular distributions of the excited product pairs-(6₁, 0)_s, (1₁, 0)_s, and (3₁, 0)_s-are remarkably similar and comprise the traits of both the peripheral mechanism and resonance pathway. Possible interpretation and implication are suggested. In addition, due to the spectral overlap of the REMPI bands and heavily congested image features, a robust data analysis method is developed, which enables us to extract the dynamics attributes of a weak feature buried in the proximate, more intense ones with high fidelity.

Fifty years of nucleophilic substitution in the gas phase

Bimolecular nucleophilic substitution ( S N 2 ) reactions have become a model system for the investigation of structure-reactivity relationships, stereochemistry, solvent influences, and detailed atomistic dynamics. In this review, the progress during five decades of experimental and theoretical research on gas phase S N 2 reactions is discussed. Many advancements of the employed methods have led to a tremendous increase in our understanding of the properties and the dynamics of these reactions. For reactions involving six atoms a quantitative agreement of the differential reactive scattering cross sections has already been achieved, in the future it is expected that even larger polyatomic reactions systems become tractable. Furthermore, studies with higher precision, improved reactant control, and a more accurate theoretical treatment of quantum effects are envisioned.

Theoretical Study of the Energy Disposal Mechanism and the State-Resolved Quantum Dynamics of the H + LiH+ → H2 + Li+ Reaction

Despite several studies in the literature, the detailed quantum state-to-state level mechanism of the astrophysically important exoergic barrierless H + LiH⁺ → H₂ + Li⁺ reaction is yet to be understood. In this work, we have investigated the energy disposal mechanism of the reaction in terms of integral reaction cross section, product internal state distributions, differential cross section, and rate constant. Fully converged and Coriolis coupled quantum mechanical calculations based on a time-dependent wave packet method have been performed at the state-to-state level on the ab initio electronic ground state potential energy surface (PES) constructed by Martinazzo et al. (J. Chem. Phys. 2003, 119, 11241-11248). The agreement between the present quantum mechanical and previous quasi-classical results is found even at very low relative translational energies of reagents. A non-statistical inverse Boltzmann vibrational distribution for the product is found. This is attributed to the "attractive" nature of the underlying PES, which facilitates the excess energy release mostly as product vibration (60-80%). The energy disposal in products is found to be unaffected by the rovibrational excitation of the reagent diatom due to the weak coupling between the vibrational modes of the reagent and the product. The mild effect of collision energy on the product energy disposal is ascribed to the effective coupling between the translational modes of the reagent and the product. It is found that the collisions lead to the formation of the product H₂ in its rovibrationally excited levels.