An encompassed representation of timescale hierarchies in first-order reaction network.

Complex networks are pervasive in various fields such as chemistry, biology, and sociology. In chemistry, first-order reaction networks are represented by a set of first-order differential equations, which can be constructed from the underlying energy landscape. However, as the number of nodes increases, it becomes more challenging to understand complex kinetics across different timescales. Hence, how to construct an interpretable, coarse-graining scheme that preserves the underlying timescales of overall reactions is of crucial importance. Here, we develop a scheme to capture the underlying hierarchical subsets of nodes, and a series of coarse-grained (reduced-dimensional) rate equations between the subsets as a function of time resolution from the original reaction network. Each of the coarse-grained representations guarantees to preserve the underlying slow characteristic timescales in the original network. The crux is the construction of a lumping scheme incorporating a similarity measure in deciphering the underlying timescale hierarchy, which does not rely on the assumption of equilibrium. As an illustrative example, we apply the scheme to four-state Markovian models and Claisen rearrangement of allyl vinyl ether (AVE), and demonstrate that the reduced-dimensional representation accurately reproduces not only the slowest but also the faster timescales of overall reactions although other reduction schemes based on equilibrium assumption well reproduce the slowest timescale but fail to reproduce the second-to-fourth slowest timescales with the same accuracy. Our scheme can be applied not only to the reaction networks but also to networks in other fields, which helps us encompass their hierarchical structures of the complex kinetics over timescales.

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.

Deciphering Time Scale Hierarchy in Reaction Networks.

Markovian dynamics on complex reaction networks are one of the most intriguing subjects in a wide range of research fields including chemical reactions, biological physics, and ecology. To represent the global kinetics from one node (corresponding to a basin on an energy landscape) to another requires information on multiple pathways that directly or indirectly connect these two nodes through the entire network. In this paper we present a scheme to extract a hierarchical set of global transition states (TSs) over a discrete-time Markov chain derived from first-order rate equations. The TSs can naturally take into account the multiple pathways connecting any pair of nodes. We also propose a new type of disconnectivity graph (DG) to capture the hierarchical organization of different time scales of reactions that can capture changes in the network due to changes in the time scale of observation. The crux is the introduction of the minimum conductance cut (MCC) in graph clustering, corresponding to the dividing surface across the network having the "smallest" transition probability between two disjoint subnetworks (superbasins on the energy landscape) in the network. We present a new combinatorial search algorithm for finding this MCC. We apply our method to a reaction network of Claisen rearrangement of allyl vinyl ether that consists of 23 nodes and 66 links (saddles on the energy landscape) connecting them. We compare the kinetic properties of our DG to those of the transition matrix of the rate equations and show that our graph can properly reveal the hierarchical organization of time scales in a network.

Kinetic Analysis for the Multistep Profiles of Organic Reactions: Significance of the Conformational Entropy on the Rate Constants of the Claisen Rearrangement.

The significance of kinetic analysis as a tool for understanding the reactivity and selectivity of organic reactions has recently been recognized. However, conventional simulation approaches that solve rate equations numerically are not amenable to multistep reaction profiles consisting of fast and slow elementary steps. Herein, we present an efficient and robust approach for evaluating the overall rate constants of multistep reactions via the recursive contraction of the rate equations to give the overall rate constants for the products and byproducts. This new method was applied to the Claisen rearrangement of allyl vinyl ether, as well as a substituted allyl vinyl ether. Notably, the profiles of these reactions contained 23 and 84 local minima, and 66 and 278 transition states, respectively. The overall rate constant for the Claisen rearrangement of allyl vinyl ether was consistent with the experimental value. The selectivity of the Claisen rearrangement reaction has also been assessed using a substituted allyl vinyl ether. The results of this study showed that the conformational entropy in these flexible chain molecules had a substantial impact on the overall rate constants. This new method could therefore be used to estimate the overall rate constants of various other organic reactions involving flexible molecules.

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.

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.