Reactivity Boundaries to Separate the Fate of a Chemical Reaction Associated with an Index-two saddle

Reactivity boundaries that divide the destination and the origin of trajectories are of crucial importance to reveal the mechanism of reactions. We investigate whether such reactivity boundaries can be extracted for higher index saddles in terms of a nonlinear canonical transformation successful for index-one saddles by using a model system with an index-two saddle. It is found that the true reactivity boundaries do not coincide with those extracted by the transformation taking into account a nonlinearity in the region of the saddle even for small perturbations, and the discrepancy is more pronounced for the less repulsive direction of the index-two saddle system. The present result indicates an importance of the global properties of the phase space to identify the reactivity boundaries, relevant to the question of what reactant and product are in phase space, for saddles with index more than one

Reactivity Boundaries to Separate the Fate of a Chemical Reaction Associated with Multiple Saddles

Reactivity boundaries that divide the origin and destination of trajectories are crucial of importance to reveal the mechanism of reactions, which was recently found to exist robustly even at high energies for index-one 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 made up 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 made up of two index-two saddle points and two index-one 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.Comment: 12 pages, 7 figure

Transition State Theory For Laser-Driven Reactions

Recent developments in transition state theory brought about by dynamical systems theory are extended to time-dependent systems such as laser-driven reactions. Using time-dependent normal form theory, the authors construct a reaction coordinate with regular dynamics inside the transition region. The conservation of the associated action enables one to extract time-dependent invariant manifolds that act as separatrices between reactive and nonreactive trajectories and thus make it possible to predict the ultimate fate of a trajectory. They illustrate the power of our approach on a driven Henon-Heiles system, which serves as a simple example of a reactive system with several open channels. The present generalization of transition state theory to driven systems will allow one to study processes such as the control of chemical reactions through laser pulses

Transition state theory for laser-driven reactions

Recent developments in Transition State Theory brought about by dynamical systems theory are extended to time-dependent systems such as laser-driven reactions. Using time-dependent normal form theory, we construct a reaction coordinate with regular dynamics inside the transition region. The conservation of the associated action enables one to extract time-dependent invariant manifolds that act as separatrices between reactive and non-reactive trajectories and thus make it possible to predict the ultimate fate of a trajectory. We illustrate the power of our approach on a driven H´enon-Heiles system, which serves as a simple example of a reactive system with several open channels. The present generalization of Transition State Theory to driven systems will allow one to study processes such as the control of chemical reactions through laser pulses

Transition state theory for laser-driven reactions

This is a pre-print.Recent developments in Transition State Theory brought about by dynamical systems theory are extended to time-dependent systems such as laser-driven reactions. Using time-dependent normal form theory, we construct a reaction coordinate with regular dynamics inside the transition region. The conservation of the associated action enables one to extract time-dependent invariant manifolds that act as separatrices between reactive and non-reactive trajectories and thus make it possible to predict the ultimate fate of a trajectory. We illustrate the power of our approach on a driven H´enon-Heiles system, which serves as a simple example of a reactive system with several open channels. The present generalization of Transition State Theory to driven systems will allow one to study processes such as the control of chemical reactions through laser pulses

Why and how do systems react in thermally fluctuating environments?

Many chemical reactions, including those of biological importance, take place in thermally fluctuating environments. Compared to isolated systems, there arise markedly different features due to the effects of energy dissipation through friction and stochastic driving by random forces reflecting the fluctuation of the environment. Investigation of how robustly the system reacts under the influence of thermal fluctuation, and elucidating the role of thermal fluctuation in the reaction are significant subjects in the study of chemical reactions. In this article, we start with overviewing the generalized Langevin equation (GLE), which has long been used and continues to be a powerful tool to describe a system surrounded by a thermal environment. It has been also generalized further to treat a nonstationary environment, in which the conventional fluctuation-dissipation theorem no longer holds. Then, within the framework of the Langevin equation we present a method recently developed to extract a new reaction coordinate that is decoupled from all the other coordinates in the region of a rank-one saddle linking the reactant and the product. The reaction coordinate is buried in nonlinear couplings among the original coordinates under the influence of stochastic random force. It was ensured that the sign of this new reaction coordinate (= a nonlinear functional of the original coordinates, velocities, friction, and random force) at any instant is sufficient to determine in which region, the reactant or the product, the system finally arrives. We also discuss how one can extend the method to extract such a coordinate from the GLE framework in stationary and nonstationary environments, where memory effects exist in dynamics of the reaction

Phase space geometry of dynamics passing through saddle coupled with spatial rotation

Nonlinear reaction dynamics through a rank-one saddle is investigated for many-particle system with spatial rotation. Based on the recently developed theories of the phase space geometry in the saddle region, we present a theoretical framework to incorporate the spatial rotation which is dynamically coupled with the internal vibrational motions through centrifugal and Coriolis interactions. As an illustrative simple example, we apply it to isomerization reaction of HCN with some nonzero total angular momenta. It is found that no-return transition state (TS) and a set of impenetrable reaction boundaries to separate the "past" and "future" of trajectories can be identified analytically under rovibrational couplings. The three components of the angular momentum are found to have distinct effects on the migration of the "anchor" of the TS and the reaction boundaries through rovibrational couplings and anharmonicities in vibrational degrees of freedom. This method provides new insights in understanding the origin of a wide class of reactions with nonzero angular momentum

Dynamic reaction coordinate in thermally fluctuating environment in the framework of the multidimensional generalized Langevin equations

A framework recently developed for the extraction of a dynamic reaction coordinate to mediate reactions buried in multidimensional Langevin equation is extended to the generalized Langevin equations without a priori assumption of the forms of the potential (in general, nonlinearly coupled systems) and the friction kernel. The equation of motion with memory effect can be transformed into an equation without memory at the cost of an increase in the dimensionality of the system, and hence the theoretical framework developed for the (nonlinear) Langevin formulation can be generalized to the non-Markovian process with colored noise. It is found that the increased dimension can be physically interpreted as effective modes of the fluctuating environment. As an illustrative example, we apply this theory to a multidimensional generalized Langevin equation for motion on the Müller-Brown potential surface with an exponential friction kernel. Numerical simulations find a boundary between the highly reactive region and the less reactive region in the space of initial conditions. The location of the boundary is found to depend significantly on both the memory kernel and the nonlinear couplings. The theory extracts a reaction coordinate whose sign determines the fate of the reaction taking into account the thermally fluctuating environments, the memory effect, and the nonlinearities. It is found that the location of the boundary of reactivity is satisfactorily reproduced as the zero of the statistical average of the new reaction coordinate, which is an analytical functional of both the original position coordinates and velocities of the system, and of the properties of the environment

Robust Existence of a Reaction Boundary to Separate the Fate of a Chemical Reaction

Nonlinear dynamics around a rank-one saddle is investigated in a high energy regime above the reaction threshold. The transition state (TS) is considered as a surface of a "point of no return" through which all reactive trajectories pass only once in the process of climbing over the saddle before being captured in the product state. A no-return TS ceases to exist above a certain high energy regime. However, even at high energies where the no-return TS can no longer exist, it is shown that "an impenetrable barrier" in the phase space robustly persists, which acts as a boundary between reactive and nonreactive trajectories. This implies that we can yet predict the fate of reactions even when the no-return TS may not exist. As an example, we show the analysis of dynamical systems theory for a hydrogen atom in crossed electric and magnetic fields