Comparing the dynamical effects of symmetric and antisymmetric stretch excitation of methane in the Cl+CH4 reaction

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.

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.

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.

Current Status of the X + C2H6 [X ≡ H, F(2P), Cl(2P), O(3P), OH] Hydrogen Abstraction Reactions: A Theoretical Review

This paper is a detailed review of the chemistry of medium-size reactive systems using the following hydrogen abstraction reactions with ethane, X + C₂H₆ → HX + C₂H₅; X ≡ H, F(²P), Cl(²P), O(³P) and OH, and focusing attention mainly on the theoretical developments. These bimolecular reactions range from exothermic to endothermic systems and from barrierless to high classical barriers of activation. Thus, the topography of the reactive systems changes from reaction to reaction with the presence or not of stabilized intermediate complexes in the entrance and exit channels. The review begins with some reflections on the inherent problems in the theory/experiment comparison. When one compares kinetics or dynamics theoretical results with experimental measures, one is testing both the potential energy surface describing the nuclei motion and the kinetics or dynamics method used. Discrepancies in the comparison may be due to inaccuracies of the surface, limitations of the kinetics or dynamics methods, and experimental uncertainties that also cannot be ruled out. The paper continues with a detailed review of some bimolecular reactions with ethane, beginning with the reactions with hydrogen atoms. The reactions with halogens present a challenge owing to the presence of stabilized intermediate complexes in the entrance and exit channels and the influence of the spin-orbit states on reactivity. Reactions with O(³P) atoms lead to three surfaces, which is an additional difficulty in the theoretical study. Finally, the reactions with the hydroxyl radical correspond to a reactive system with ten atoms and twenty-four degrees of freedom. Throughout this review, different strategies in the development of analytical potential energy surfaces describing these bimolecular reactions have been critically analyzed, showing their advantages and limitations. These surfaces are fitted to a large number of ab initio calculations, and we found that a huge number of calculations leads to accurate surfaces, but this information does not guarantee that the kinetics and dynamics results match the experimental measurements.

Imaging the Mode-Specificity in Cl + CH3D(v1-I, v1-II, v4 = 1; |jK⟩ = |10⟩) → CH2D(41) + HCl(v)

We report a crossed-beam imaging experiment on the title reactions at two collisional energies (E_c) of 5.3 and 10 kcal mol^-1. Both the integral cross sections relative to the ground-state reactivity and the differential cross sections were measured and compared. We found that one-quantum excitations of the CH₃-stretching vibrations of the CH₃D reagent exerted profound mode-specificity in forming the umbrella-mode-excited CH₂D(4₁) products with the vibrational efficacy of v₄ > v₁-I > v₁-II at both E_c values. The concomitantly formed HCl(v) coproducts were vibrationally cold. Interestingly, the branching ratios of (v = 1)/(v = 0) appeared invariant to the initial stretch-modes of excitation at E_c = 5.3 kcal mol^-1, yet exhibited a pronounced mode-specific dependency in the order of v₁-II > v₁-I > v₄ at E_c = 10.3 kcal mol^-1. This large and E_c-dependent disparity between the two Fermi-coupled reagents, v₁-I and v₁-II, is particularly significant and could be another facet─in addition to that in the recently reported vibrational enhancement factors─of the Fermi-phase-induced interference effect manifested in the product vibrational branching ratio. The pair-correlated angular distributions (v_CH2D, v_HCl)_s = (4₁, 0)_s in the three stretch-excited reactions were globally alike and resembled that of the ground-state reaction pair (0₀, 0)_g, suggestive of a direct abstraction mechanism of the peripheral type. This is in sharp contrast to all other vibrationally excited pairs of (1₁, 0)_s, (3₁, 0)_s, and (6₁, 0)_s previously reported in the CH₂D + HCl isotopic channel, for which both the direct abstraction and a time-delayed resonance pathway partake.

Describing Polytopal Rearrangements of Fluxional Molecules with Curvilinear Coordinates Derived from Normal Vibrational Modes: A Conceptual Extension of Cremer-Pople Puckering Coordinates

In this work a new curvilinear coordinate system is presented for the comprehensive description of polytopal rearrangements of N-coordinate compounds (N = 4-7) and systems containing an N-coordinate subunit. It is based on normal vibrational modes and a natural extension of the Cremer-Pople puckering coordinates ( J. Am. Chem. Soc. 1975, 97, 1354) together with the Zou-Izotov-Cremer deformation coordinates ( J. Phys. Chem. A 2011, 115, 8731) for ring structures to N-coordinate systems. We demonstrate that the new curvilinear coordinates are ideal reaction coordinates describing fluxional rearrangement pathways by revisiting the Berry pseudorotation and the lever mechanism in sulfur tetrafluoride, the Berry pseudorotation and two Muetterties' mechanisms in pentavalent compounds, the chimeric pseudorotation in iodine pentafluoride, Bailar and Ray-Dutt twists in hexacoordinate tris-chelates as well as the Bartell mechanism in iodine heptafluoride. The results of our study reveal that this dedicated curvilinear coordinate system can be applied to most coordination compounds opening new ways for the systematic modeling of fluxional processes.

Quantum dynamics of polyatomic dissociative chemisorption on transition metal surfaces: mode specificity and bond selectivity

Recent advances in quantum dynamical characterization of polyatomic dissociative chemisorption on accurate global potential energy surfaces are critically reviewed.

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

CONSPECTUS: Mode specificity is defined by the differences in reactivity due to excitations in various reactant modes, while bond selectivity refers to selective bond breaking in a reaction. These phenomena not only shed light on reaction dynamics but also open the door for laser control of reactions. The existence of mode specificity and bond selectivity in a reaction indicates that not all forms of energy are equivalent in promoting the reactivity, thus defying a statistical treatment. They also allow the enhancement of reactivity and control product branching ratio. As a result, they are of central importance in chemistry. This Account discusses recent advances in our understanding of these nonstatistical phenomena. In particular, the newly proposed sudden vector projection (SVP) model and its applications are reviewed. The SVP model is based on the premise that the collision in many direct reactions is much faster than intramolecular vibrational energy redistribution in the reactants. In such a sudden limit, the coupling of a reactant mode with the reaction coordinate at the transition state, which dictates its ability to promote the reaction, is approximately quantified by the projection of the former onto the latter. The SVP model can be considered as a generalization of the venerable Polanyi's rules, which are based on the location of the barrier. The SVP model is instead based on properties of the saddle point and as a result capable of treating the translational, rotational, and multiple vibrational modes in reactions involving polyatomic reactants. In case of surface reactions, the involvement of surface atoms can also be examined. Taking advantage of microscopic reversibility, the SVP model has also been used to predict product energy disposal in reactions. This simple yet powerful rule of thumb has been successfully demonstrated in many reactions including uni- and bimolecular reactions in the gas phase and gas-surface reactions. The success of the SVP model underscores the importance of the transition state in controlling mode-specific and bond-selective chemistry.

Relative efficacy of vibrational vs. translational excitation in promoting atom-diatom reactivity: rigorous examination of Polanyi's rules and proposition of sudden vector projection (SVP) model

To provide a systematic and rigorous re-examination of the well-known Polanyi's rules, excitation functions of several A + BC(v = 0, 1) reactions are determined using the Chebyshev real wave packet method on accurate potential energy surfaces. Reactions with early (F + H2 and F + HCl), late (Cl + H2), and central (H∕D∕Mu + H2, where Mu is a short-lived light isotope of H) barriers are represented. Although Polanyi's rules are in general consistent with the quantum dynamical results, their predictions are strictly valid only in certain energy ranges divided by a cross-over point. In particular, vibrational excitation of the diatomic reactant typically enhances reactivity more effectively than translational excitation at high energies, while reverse is true at low energies. This feature persists irrespective of the barrier location. A sudden vector projection model is proposed as an alternative to Polanyi's rules. It is found to give similar, but more quantitative, predictions about mode selectivity in these reactions, and has the advantage to be extendible to reactions involving polyatomic molecules.

Vibrational Enhancement Factor of the Cl + CHD3(v1 = 1) Reaction: Rotational-Probe Effects

The vibrational enhancement factor in the Cl + CHD3(v1 = 1) reaction is revisited over the collisional energy range of 2-5.9 kcal mol(-1). Contrary to the previous results obtained by probing the low-|N, K⟩ states of CD3(v = 0) products, CH stretching excitation becomes more efficacious than the same amount of translational energy in promoting the HCl(v) + CD3(v = 0) product pairs when all-|N, K⟩ states are probed. Whereas the new vibrational enhancement factors, which are three to four times larger than the previous report, agree reasonably well with a recent reduced-dimensionality quantum dynamics calculation, a cautious note is made on the different initial |J,K⟩ rotational selections of the CHD3 reactants in the present theory-experiment comparison.

Mode Selectivity for a "Central" Barrier Reaction: Eight-Dimensional Quantum Studies of the O((3)P) + CH4 → OH + CH3 Reaction on an Ab Initio Potential Energy Surface

The dynamics of a combustion reaction, namely, O((3)P) + CH4 → OH + CH3, is investigated with an eight-dimensional quantum model that includes representatives of all vibrational modes of CH4 and with a full-dimensional quasi-classical trajectory (QCT) method. The calculated excitation functions for the ground vibrational state CH4 agree well with experiment. Both quantum and QCT results suggest that excitation of the stretching modes of CH4 enhances the reaction, while the bending and umbrella modes have a smaller impact on reactivity, again consistent with experimental findings. However, none of the vibrational excitations has comparable efficiency in promoting the reaction as translational energy.

STEREODYNAMICS STUDY OF THE REACTION Cl(2p3/2) + C2D6 (v = 0, j = 0) → DCl + C2D5

The product angular momentum polarization of the Cl + C2D6 → DCl + C2D5 reaction is calculated via the quasiclassical trajectory (QCT) method at the collision energy of 0.25 eV. A new London–Eyring–Polanyi–Sato (LEPS) potential energy surface (PES) is used in this reaction. There is a "late" barrier and a "deep" well on this new LEPS PES. The four polarization-dependent "generalized" differential cross sections (PDDCSs) are presented in the center-of-mass frame. In the meantime, the distributions of P(ϕr), P(θr), and P(θr, ϕr) are calculated. The calculations are in good agreement with the experimental data. In addition, the rotational alignment factors [Formula: see text], [Formula: see text], and [Formula: see text] in the stationary-target frame (STF) are also calculated.

False estimates of stimulated Raman pumping efficiency caused by the optical Stark effect

One technique for measuring the fraction of molecules pumped to the excited state in stimulated Raman pumping (SRP) is to record the depletion of molecules in the lower state by resonance enhanced multiphoton ionization (REMPI). The presence of electric fields on the order of 10(7) V/cm arising from the pulsed SRP laser beams is sufficient to shift the line position of the REMPI transition to such an extent that the estimate of the pumping efficiency is overestimated unless this shift is accounted for.

Quasiclassical trajectory calculations of correlated product distributions for the F + CHD3(v1 = 0, 1) reactions using an ab initio potential energy surface

We report quasiclassical trajectory (QCT) calculations of the correlated product distributions and branching ratios of the reactions F+CHD(3)(v(1)=0,1)-->HF(v)+CD(3)(v) and DF(v)+CHD(2)(v) using a recently published ab initio-based full-dimensional global potential energy surface [G. Czako et al., J. Chem. Phys. 130, 084301 (2009)]. Harmonic normal mode analysis is done for the methyl products to determine the classical actions of each normal mode and then standard histogram binning and Gaussian binning (GB) methods are employed to obtain quantum state-resolved probabilities of the products. QCT calculations have been performed for both the vibrationally ground state and the CH stretching excited F+CHD(3)(v(1)=0,1) reactions at eight different collision energies in the 0.5-7.0 kcal/mol range. HF and DF vibrationally state-resolved rotational and angular distributions, CD(3) and CHD(2) mode-specific vibrational distributions, and correlated vibrationally state-specific distributions for the product pairs have been obtained and the correlated results were compared to the experiment. We find that the use of GB can be advantageous especially in the threshold regions. The CH stretching excitation in the reactant does not change the CD(3) vibrational distributions significantly, whereas the HF molecules become vibrationally and rotationally hotter. On the other hand in the case of the DF+CHD(2) channel the initially excited CH stretch appears mainly "intact" in the CHD(2) product and the DF distributions are virtually the same as formed from the ground state CHD(3) reaction. The computed results qualitatively agree with recent crossed molecular beam experiment [W. Zhang et al., Science 325, 303 (2009)] that (a) CHD(2)(v(1)=1) is the most populated product state of the F+CHD(3)(v(1)=1) reaction and this reaction produces much less CHD(2)(v=0) compared to the reaction F+CHD(3)(v=0); (b) the cross section ratio of CHD(2)(v(1)=1):CHD(2)(v=0) formed from the reactions F+CHD(3)(v(1)=1):F+CHD(3)(v=0) is less than 1 and shows little collision energy dependency; (c) the reactant CH stretch excitation increases the DF:HF ratio at low collision energies; (d) the correlated vibrational and angular distributions for DF(v)+CHD(2)(v(1)=0,1) from the ground state and stretch-excited reactions, respectively, are almost identical.

Searching for resonances in the reaction Cl+CH4→HCl+CH3: Quantum versus quasiclassical dynamics and comparison with experiments

A quantum-mechanical (QM) and quasiclassical trajectory (QCT) study was performed on the title reaction, using a pseudotriatomic ab initio based surface. Probabilities and integral cross sections present some clear peaks versus the collision energy E(col), which we assign to Feshbach resonances of the transition state, where the light H atom oscillates between the heavy Cl and CH(3) groups. For ground-state reactants, reactivity is essentially of quantum origin (QCT observables and oscillations are smaller, or much smaller, than QM ones), and the calculated integral cross section and product distributions are in reasonable agreement with the experiment. The reaction occurs through an abstraction mechanism, following both a direct and an indirect mechanism. The quasiclassical trajectory calculations show the participation of a short-lived collision complex in the microscopic reaction mechanism. Finally, QCT differential cross sections of Cl+CH(4)-->HCl (nu(')=0 and 1)+CH(3) oscillate versus E(col), whereas experimentally this only occurs for HCl (nu(')=1). This theoretical result and other oscillating properties found here could, however, be related to the existence of a Feshbach resonance for the production of HCl (nu(')=1), as suggested by experimentalists.

Role of the C-H stretch mode excitation in the dynamics of the Cl + CHD3 reaction: a quasi-classical trajectory calculation

To analyze the effect of the C-H stretch mode excitation on the dynamics of the Cl + CHD3 gas-phase abstraction reaction, an exhaustive state-to-state dynamics study was performed. This reaction can evolve along two channels: H-abstraction, CD3 + ClH, and D-abstraction, CHD2 + ClD. On an analytical potential energy surface constructed previously by our group, named PES-2005, quasi-classical trajectory calculations were performed at a collision energy of 0.18 eV, including corrections to avoid zero-point energy leakage along the trajectories. First, strong coupling between different vibrational modes in the entry valley was observed; i.e., the reaction is vibrationally nonadiabatic. Second, for the ground-state CHD3(nu=0) reaction, the diatomic fragments appeared in their ground states, and the H- and D-abstraction reactions showed similar reactivities. However, when the reactivity per atom is considered, the H is three times more reactive than the D atom. Third, when the C-H stretch mode is excited by one quantum, CHD3(nu1=1), the H-abstraction is strongly favored, and the C-H stretch excitation is maintained in the product CHD2(nu1=1) + ClD channel; i.e., the reaction shows mode selectivity, reproducing the experimental evidence, and also the reactivity of the vibrational ground state is increased, in agreement with experiment. Fourth, the state-to-state angular distributions of the CD3 and CHD2 products showed the products to be practically sideways for the reactant ground state, while the C-H excitation yielded a more forward scattering, reproducing the experimental data. The role of the zero-point energy correction was also analyzed, and we find that the dynamics results are very sensitive on how the ZPE issue is treated. Finally, a comparison is made with the similar H + CHD3(nu1=0,1) and Cl + CH4(nu1=0,1) reactions.

Quasi-classical trajectory calculations analyzing the dynamics of the C-H stretch mode excitation in the H+CHD3 reaction

A state-to-state dynamics study was performed at a collision energy of 1.53 eV to analyze the effect of the C-H stretch mode excitation on the dynamics of the gas-phase H+CHD3 reaction, which can evolve along two channels, H-abstraction, CD3+H2, and D-abstraction, CHD2+HD. Quasi-classical trajectory calculations were performed on an analytical potential energy surface constructed previously by our group. First, strong coupling between different vibrational modes in the entry channel was observed; i.e., the reaction is non-adiabatic. Second, we found that the C-H stretch mode excitation has little influence on the product rotational distributions for both channels, and on the vibrational distribution for the CD3+H2 channel. However, it has significant influence on the product vibrational distribution for the CHD2+HD channel, where the C-H stretch excitation is maintained in the products, i.e., the reaction shows mode selectivity, reproducing the experimental evidence. Third, the C-H stretch excitation by one quantum increases the reactivity of the vibrational ground-state, in agreement with experiment. Fourth, the state-to-state angular distributions of the CD3 and CHD2 products are reported, finding that for the reactant ground-state the products are practically sideways, whereas the C-H excitation yields a more forward scattering.