Effects of Roaming Trajectories on the Transition State Theory Rates of a Reduced-Dimensional Model of Ketene Isomerization

Dynamics and decay rates of a time-dependent two-saddle system

The framework of transition state theory (TST) provides a powerful way for analyzing the dynamics of physical and chemical reactions. While TST has already been successfully used to obtain reaction rates for systems with a single time-dependent saddle point, multiple driven saddles have proven challenging because of their fractal-like phase space structure. This paper presents the construction of an approximately recrossing-free dividing surface based on the normally hyperbolic invariant manifold in a time-dependent two-saddle model system. Based on this, multiple methods for obtaining instantaneous (time-resolved) decay rates of the underlying activated complex are presented and their results discussed.

Identifying reaction pathways in phase space via asymptotic trajectories

In this paper, we revisit the concepts of the reactivity map and the reactivity bands as an alternative to the use of perturbation theory for the determination of the phase space geometry of chemical reactions. We introduce a reformulated metric, called the asymptotic trajectory indicator, and an efficient algorithm to obtain reactivity boundaries. We demonstrate that this method has sufficient accuracy to reproduce phase space structures such as turnstiles for a 1D model of the isomerization of ketene in an external field. The asymptotic trajectory indicator can be applied to higher dimensional systems coupled to Langevin baths as we demonstrate for a 3D model of the isomerization of ketene.

Roaming Reactions and Dynamics in the van der Waals Region

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

Theories and simulations of roaming

The phenomenon of roaming in chemical reactions has now become both commonly observed in experiment and extensively supported by theory and simulations. Roaming occurs in highly-excited molecules when the trajectories of atomic motion often bypass the minimum energy pathway and produce reaction in unexpected ways from unlikely geometries. The prototypical example is the unimolecular dissociation of formaldehyde (H₂CO), in which the "normal" reaction proceeds through a tight transition state to yield H₂ + CO but for which a high fraction of dissociations take place via a "roaming" mechanism in which one H atom moves far from the HCO, almost to dissociation, and then returns to abstract the second H atom. We review below the theories and simulations that have recently been developed to address and understand this new reaction phenomenon.

Toward Understanding the Roaming Mechanism in H + MgH → Mg + HH Reaction

The roaming mechanism in the reaction H + MgH →Mg + HH is investigated by classical and quantum dynamics employing an accurate ab initio three-dimensional ground electronic state potential energy surface. The reaction dynamics are explored by running trajectories initialized on a four-dimensional dividing surface anchored on three-dimensional normally hyperbolic invariant manifold associated with a family of unstable orbiting periodic orbits in the entrance channel of the reaction (H + MgH). By locating periodic orbits localized in the HMgH well or involving H orbiting around the MgH diatom, and following their continuation with the total energy, regions in phase space where reactive or nonreactive trajectories may be trapped are found. In this way roaming reaction pathways are deduced in phase space. Patterns similar to periodic orbits projected into configuration space are found for the quantum bound and resonance eigenstates. Roaming is attributed to the capture of the trajectories in the neighborhood of certain periodic orbits. The complex forming trajectories in the HMgH well can either return to the radical channel or "roam" to the MgHH minimum from where the molecule may react.

Deconstructing field-induced ketene isomerization through Lagrangian descriptors

Phase space contours (shown in color) constructed using the method of Lagrangian descriptors resolve the separatrices governing state transitions on the reaction-path potential energy surface (shown in white) for field-induced ketene isomerization.

Phase Space Structures Explain Hydrogen Atom Roaming in Formaldehyde Decomposition

We re-examine the prototypical roaming reaction--hydrogen atom roaming in formaldehyde decomposition--from a phase space perspective. Specifically, we address the question "why do trajectories roam, rather than dissociate through the radical channel?" We describe and compute the phase space structures that define and control all possible reactive events for this reaction, as well as provide a dynamically exact description of the roaming region in phase space. Using these phase space constructs, we show that in the roaming region, there is an unstable periodic orbit whose stable and unstable manifolds define a conduit that both encompasses all roaming trajectories exiting the formaldehyde well and shepherds them toward the H2···CO well.

Roaming dynamics in ion-molecule reactions: phase space reaction pathways and geometrical interpretation

A model Hamiltonian for the reaction CH4(+) -> CH3(+) + H, parametrized to exhibit either early or late inner transition states, is employed to investigate the dynamical characteristics of the roaming mechanism. Tight/loose transition states and conventional/roaming reaction pathways are identified in terms of time-invariant objects in phase space. These are dividing surfaces associated with normally hyperbolic invariant manifolds (NHIMs). For systems with two degrees of freedom NHIMS are unstable periodic orbits which, in conjunction with their stable and unstable manifolds, unambiguously define the (locally) non-recrossing dividing surfaces assumed in statistical theories of reaction rates. By constructing periodic orbit continuation/bifurcation diagrams for two values of the potential function parameter corresponding to late and early transition states, respectively, and using the total energy as another parameter, we dynamically assign different regions of phase space to reactants and products as well as to conventional and roaming reaction pathways. The classical dynamics of the system are investigated by uniformly sampling trajectory initial conditions on the dividing surfaces. Trajectories are classified into four different categories: direct reactive and non-reactive trajectories, which lead to the formation of molecular and radical products respectively, and roaming reactive and non-reactive orbiting trajectories, which represent alternative pathways to form molecular and radical products. By analysing gap time distributions at several energies, we demonstrate that the phase space structure of the roaming region, which is strongly influenced by nonlinear resonances between the two degrees of freedom, results in nonexponential (nonstatistical) decay.

Dynamic effects dictate the mechanism and selectivity of dehydration-rearrangement reactions of protonated alcohols [Me2 (R)CCH(OH2 )Me](+) (R=Me, Et, iPr) in the gas phase

The gas-phase dehydration-rearrangement (DR) reactions of protonated alcohols [Me2 (R)CCH(OH2 )Me](+) [R=Me (ME), Et (ET), and iPr (I-PR)] were studied by using static approaches (intrinsic reaction coordinate (IRC), Rice-Ramsperger-Kassel-Marcus theory) and dynamics (quasiclassical trajectory) simulations at the B3LYP/6-31G(d) level of theory. The concerted mechanism involves simultaneous water dissociation and alkyl migration, whereas in the stepwise reaction pathway the dehydration step leads to a secondary carbocation intermediate followed by alkyl migration. Internal rotation (IR) can change the relative position of the migrating alkyl group and the leaving group (water), so distinct products may be obtained: [Me(R)CCH(Me)Me⋅⋅⋅OH2 ](+) and [Me(Me)CCH(R)Me⋅⋅⋅OH2 ](+) . The static approach predicts that these reactions are concerted, with the selectivity towards these different products determined by the proportion of the conformers of the initial protonated alcohols. These selectivities are explained by the DR processes being much faster than IR. These results are in direct contradiction with the dynamics simulations, which indicate a predominantly stepwise mechanism and selectivities that depend on the alkyl groups and dynamics effects. Indeed, despite the lifetimes of the secondary carbocations being short (<0.5 ps), IR can take place and thus provide a rich selectivity. These different selectivities, particularly for ET and I-PR, are amenable to experimental observation and provide evidence for the minor role played by potential-energy surface and the relevance of the dynamics effects (non-IRC pathways, IR) in determining the reaction mechanisms and product distribution (selectivity).