Crossing the dividing surface of transition state theory. II. Recrossing times for the atom-diatom interaction

Toward an Accurate Black-Box Tool for the Kinetics of Gas-Phase Reactions Involving Barrier-less Elementary Steps

An enhanced computational protocol has been devised for the accurate characterization of gas-phase barrier-less reactions in the framework of the reaction-path (RP) and variable reaction coordinate variational transition-state theory. In particular, the synergistic combination of density functional theory and Monte Carlo sampling to optimize reactive fluxes led to a reliable yet effective computational workflow. A black-box strategy has been developed for selecting the most suited density functional with reference to a high-level one-dimensional reference potential. At the same time, different descriptions of hindered rotations are automatically selected, depending on the corresponding harmonic frequencies along the RP. The performance of the new tool is investigated by means of two prototypical reactions involving different degrees of static and dynamic correlation, namely, H2S + Cl and CH3 + CH3. The remarkable agreement of the computed kinetic parameters with the available experimental data confirms the accuracy and robustness of the proposed approach. Together with their intrinsic interest, these results also pave the way toward systematic investigations of gas-phase reactions involving barrier-less elementary steps by a reliable, user-friendly tool, which can be confidently used also by nonspecialists.

Phase space geometry of isolated to condensed chemical reactions

The complexity of gas and condensed phase chemical reactions has generally been uncovered either approximately through transition state theories or exactly through (analytic or computational) integration of trajectories. These approaches can be improved by recognizing that the dynamics and associated geometric structures exist in phase space, ensuring that the propagator is symplectic as in velocity-Verlet integrators and by extending the space of dividing surfaces to optimize the rate variationally, respectively. The dividing surface can be analytically or variationally optimized in phase space, not just over configuration space, to obtain more accurate rates. Thus, a phase space perspective is of primary importance in creating a deeper understanding of the geometric structure of chemical reactions. A key contribution from dynamical systems theory is the generalization of the transition state (TS) in terms of the normally hyperbolic invariant manifold (NHIM) whose geometric phase-space structure persists under perturbation. The NHIM can be regarded as an anchor of a dividing surface in phase space and it gives rise to an exact non-recrossing TS theory rate in reactions that are dominated by a single bottleneck. Here, we review recent advances of phase space geometrical structures of particular relevance to chemical reactions in the condensed phase. We also provide conjectures on the promise of these techniques toward the design and control of chemical reactions.

Crossing the dividing surface of transition state theory. IV. Dynamical regularity and dimensionality reduction as key features of reactive trajectories

The atom-diatom interaction is studied by classical mechanics using Jacobi coordinates (R, r, θ). Reactivity criteria that go beyond the simple requirement of transition state theory (i.e., P_R* > 0) are derived in terms of specific initial conditions. Trajectories that exactly fulfill these conditions cross the conventional dividing surface used in transition state theory (i.e., the plane in configuration space passing through a saddle point of the potential energy surface and perpendicular to the reaction coordinate) only once. Furthermore, they are observed to be strikingly similar and to form a tightly packed bundle of perfectly collimated trajectories in the two-dimensional (R, r) configuration space, although their angular motion is highly specific for each one. Particular attention is paid to symmetrical transition states (i.e., either collinear or T-shaped with C_2v symmetry) for which decoupling between angular and radial coordinates is observed, as a result of selection rules that reduce to zero Coriolis couplings between modes that belong to different irreducible representations. Liapunov exponents are equal to zero and Hamilton's characteristic function is planar in that part of configuration space that is visited by reactive trajectories. Detailed consideration is given to the concept of average reactive trajectory, which starts right from the saddle point and which is shown to be free of curvature-induced Coriolis coupling. The reaction path Hamiltonian model, together with a symmetry-based separation of the angular degree of freedom, provides an appropriate framework that leads to the formulation of an effective two-dimensional Hamiltonian. The success of the adiabatic approximation in this model is due to the symmetry of the transition state, not to a separation of time scales. Adjacent trajectories, i.e., those that do not exactly fulfill the reactivity conditions have similar characteristics, but the quality of the approximation is lower. At higher energies, these characteristics persist, but to a lesser degree. Recrossings of the dividing surface then become much more frequent and the phase space volumes of initial conditions that generate recrossing-free trajectories decrease. Altogether, one ends up with an additional illustration of the concept of reactive cylinder (or conduit) in phase space that reactive trajectories must follow. Reactivity is associated with dynamical regularity and dimensionality reduction, whatever the shape of the potential energy surface, no matter how strong its anharmonicity, and whatever the curvature of its reaction path. Both simplifying features persist during the entire reactive process, up to complete separation of fragments. The ergodicity assumption commonly assumed in statistical theories is inappropriate for reactive trajectories.

Crossing the dividing surface of transition state theory. III. Once and only once. Selecting reactive trajectories

The purpose of the present work is to determine initial conditions that generate reacting, recrossing-free trajectories that cross the conventional dividing surface of transition state theory (i.e., the plane in configuration space passing through a saddle point of the potential energy surface and perpendicular to the reaction coordinate) without ever returning to it. Local analytical equations of motion valid in the neighborhood of this planar surface have been derived as an expansion in Poisson brackets. We show that the mere presence of a saddle point implies that reactivity criteria can be quite simply formulated in terms of elements of this series, irrespective of the shape of the potential energy function. Some of these elements are demonstrated to be equal to a sum of squares and thus to be necessarily positive, which has a profound impact on the dynamics. The method is then applied to a three-dimensional model describing an atom-diatom interaction. A particular relation between initial conditions is shown to generate a bundle of reactive trajectories that form reactive cylinders (or conduits) in phase space. This relation considerably reduces the phase space volume of initial conditions that generate recrossing-free trajectories. Loci in phase space of reactive initial conditions are presented. Reactivity is influenced by symmetry, as shown by a comparative study of collinear and bent transition states. Finally, it is argued that the rules that have been derived to generate reactive trajectories in classical mechanics are also useful to build up a reactive wave packet.

Crossing the dividing surface of transition state theory. I. Underlying symmetries and motion coordination in multidimensional systems

The objective of the present paper is to show the existence of motion coordination among a bundle of trajectories crossing a saddle point region in the forward direction. For zero total angular momentum, no matter how complicated the anharmonic part of the potential energy function, classical dynamics in the vicinity of a transition state is constrained by symmetry properties. Trajectories that all cross the plane R = R* at time t = 0 (where R* denotes the position of the saddle point) with the same positive translational momentum P(R*) can be partitioned into two sets, denoted "gerade" and "ungerade," which coordinate their motions. Both sets have very close average equations of motion. This coordination improves tremendously rapidly as the number of degrees of freedom increases. This property can be traced back to the existence of time-dependent constants of the motion.