Reactivity boundaries to separate the fate of a chemical reaction associated with an index-two saddle

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.

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.

Can reactions follow non-traditional second-order saddle pathways avoiding transition states?

We report here an ab initio (CASSCF/6-31+G*) trajectory simulation study on the mechanisms of the denitrogenation of 1-pyrazoline and its subsituted analogue that reveals reaction pathways via a high energy second-order saddle (SOS) region. This mechanism involves the molecule adopting a five-membered planar structure contrary to the traditional boat-like transition state. The SOS offers a trifurcation point where a pathway branches into three, different from the single pathway associated with a transitions state. We observe that the molecules following the SOS path exhibit distinctive dynamical features and form products with high translational energies and low rotational energies compared to those following the traditional pathways. In addition, the SOS pathway provides an alternative mechanism for the formation of stereo-selective products. Interestingly, although the reaction proceeds via a trimethylene diradical intermediate, the simulations show that the product cyclopropane is formed with a major single inversion of the configuration consistent with experimental observations. They also reveal mechanisms that do not follow the minimum energy paths and exhibit non-statistical dissociation dynamics.

Reactivity boundaries for chemical reactions associated with higher-index and multiple saddles

Reactivity boundaries that divide the origin and destination of trajectories are of crucial importance to reveal the mechanism of reactions, which was recently found to exist robustly even at high energies for index 1 saddles [Phys. Rev. Lett. 105, 048304 (2010)]. Here we revisit the concept of the reactivity boundary and propose a more general definition that can involve a single reaction associated with a bottleneck composed of higher-index saddles and/or several saddle points with different indices, where the normal form theory, based on expansion around a single stationary point, does not work. We numerically demonstrate the reactivity boundary by using a reduced model system of the H(5)(+) cation where the proton exchange reaction takes place through a bottleneck composed of two index 2 saddle points and two index 1 saddle points. The cross section of the reactivity boundary in the reactant region of the phase space reveals which initial conditions are effective in making the reaction happen and thus sheds light on the reaction mechanism.