Uncovering the Geometry of Barrierless Reactions Using Lagrangian Descriptors

Quantifying chaos using Lagrangian descriptors

We present and validate simple and efficient methods to estimate the chaoticity of orbits in low-dimensional conservative dynamical systems, namely, autonomous Hamiltonian systems and area-preserving symplectic maps, from computations of Lagrangian descriptors (LDs) on short time scales. Two quantities are proposed for determining the chaotic or regular nature of orbits in a system's phase space, which are based on the values of the LDs of these orbits and of nearby ones: The difference and ratio of neighboring orbits' LDs. Using as generic test models the prototypical two degree of freedom Hénon-Heiles system and the two-dimensional standard map, we find that these indicators are able to correctly characterize the chaotic or regular nature of orbits to better than 90% agreement with results obtained by implementing the Smaller Alignment Index (SALI) method, which is a well-established chaos detection technique. Further investigating the performance of the two introduced quantities, we discuss the effects of the total integration time and of the spacing between the used neighboring orbits on the accuracy of the methods, finding that even typical short time, coarse-grid LD computations are sufficient to provide reliable quantification of the systems' chaotic component, using less CPU time than the SALI. In addition to quantifying chaos, the introduced indicators have the ability to reveal details about the systems' local and global chaotic phase space structure. Our findings clearly suggest that LDs can also be used to quantify and investigate chaos in continuous and discrete low-dimensional conservative dynamical systems.

Locating Transition Zone in Phase Space

Understanding the reaction mechanism is required for better control of chemical reactions and is usually achieved by locating transition states (TSs) along a proper one-dimensional coordinate called reaction coordinate (RC). The identification of RC can be very difficult for high-dimensional realistic systems. A number of methods have been proposed to tackle this problem. A machine learning method is developed here to incorporate the influence of velocity on the reaction process. The method is also free of the unbalanced label problem resulting from the rather low fraction of configurations near the TS and can be easily extended to large systems. It locates the transition zone in the phase space and defines the dividing surface with a high transmission coefficient. Moreover, considering that the reaction environment can not only change the reaction path but also activate the reactive mode through energy transfer, we devise two measures to quantify the influence of these two factors on the reaction process and find that solvents can assist the reaction by directly doing work along the reactive mode. Not surprisingly, there is a positive correlation between the efficiency of energy transfer into the reactive mode and the reaction rate.

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.

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.

Detecting reactive islands in a system-bath model of isomerization

In this article, we study the conformational isomerization in a solvent using a system-bath model where the phase space structures relevant for the reaction dynamics are revealed. These phase space structures are an integral part of understanding the reaction mechanism, that is the pathways that reactive trajectories undertake, in the presence of a solvent. Our approach involves detecting the analogs of the reactive islands first discussed in the works by Davis, Marston, De Leon, Berne and coauthors in the system-bath model using Lagrangian descriptors. We first present the structure of the reactive islands for the two degrees of freedom system modelling isomerization in the absence of the bath using direct computation of cylindrical (tube) manifolds and verify the Lagrangian descriptor method for detecting the reactive islands. The hierarchy of the reactive islands as indicated in the recent work by Patra and Keshavamurthy is shown to be related to the temporal features in committor probabilities. Next, we investigate the influence of the solvent on the reactive islands that we previously revealed for the two degrees of freedom system and discuss the use of the Lagrangian descriptor in the high-dimensional phase space of the system-bath model.

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.

Neural network approach for the dynamics on the normally hyperbolic invariant manifold of periodically driven systems

Chemical reactions in multidimensional systems are often described by a rank-1 saddle, whose stable and unstable manifolds intersect in the normally hyperbolic invariant manifold (NHIM). Trajectories started on the NHIM in principle never leave this manifold when propagated forward or backward in time. However, the numerical investigation of the dynamics on the NHIM is difficult because of the instability of the motion. We apply a neural network to describe time-dependent NHIMs and use this network to stabilize the motion on the NHIM for a periodically driven model system with two degrees of freedom. The method allows us to analyze the dynamics on the NHIM via Poincaré surfaces of section (PSOS) and to determine the transition-state (TS) trajectory as a periodic orbit with the same periodicity as the driving saddle, viz. a fixed point of the PSOS surrounded by near-integrable tori. Based on transition state theory and a Floquet analysis of a periodic TS trajectory we compute the rate constant of the reaction with significantly reduced numerical effort compared to the propagation of a large trajectory ensemble.

Finding normally hyperbolic invariant manifolds in two and three degrees of freedom with Hénon-Heiles-type potential

We present a method based on a Lagrangian descriptor for revealing the high-dimensional phase space structures that are of interest in nonlinear Hamiltonian systems with index-1 saddle. These phase space structures include a normally hyperbolic invariant manifold and its stable and unstable manifolds, which act as codimension-1 barriers to phase space transport. In this article, finding the invariant manifolds in high-dimensional phase space will constitute identifying coordinates on these invariant manifolds. The method of Lagrangian descriptor is demonstrated by applying to classical two and three degrees of freedom Hamiltonian systems which have implications for myriad applications in chemistry, engineering, and physics.

Unveiling the chaotic structure in phase space of molecular systems using Lagrangian descriptors

We explore here the feasibility of using the recently introduced Lagrangian descriptors [A. M. Mancho et al., Commun. Nonlinear Sci. Numer. Simul. 18, 3530 (2013)1007-570410.1016/j.cnsns.2013.05.002] to unveil the usually rich dynamics taking place in the vibrations of molecular systems, especially if they are floppy. The principal novelty of our work is the inclusion of p norms in the definition of the descriptors in this kind of system, which greatly enhances their power to discern among the different structures existing in the phase space. As an illustration we use the LiCN molecule described by realistic potentials in two and three dimensions, which exhibits chaotic motion within a mixed phase space in the isomerization between the two wells corresponding to the linear isomer stable configurations, LiNC and LiCN. In particular, we pay special attention to the manifolds emerging from the unstable fixed point between the corresponding isomer wells, and also to the marginally stable structures around a parabolic point existing near the LiNC well.

Lagrangian descriptors of driven chemical reaction manifolds

The persistence of a transition state structure in systems driven by time-dependent environments allows the application of modern reaction rate theories to solution-phase and nonequilibrium chemical reactions. However, identifying this structure is problematic in driven systems and has been limited by theories built on series expansion about a saddle point. Recently, it has been shown that to obtain formally exact rates for reactions in thermal environments, a transition state trajectory must be constructed. Here, using optimized Lagrangian descriptors [G. T. Craven and R. Hernandez, Phys. Rev. Lett. 115, 148301 (2015)PRLTAO0031-900710.1103/PhysRevLett.115.148301], we obtain this so-called distinguished trajectory and the associated moving reaction manifolds on model energy surfaces subject to various driving and dissipative conditions. In particular, we demonstrate that this is exact for harmonic barriers in one dimension and this verification gives impetus to the application of Lagrangian descriptor-based methods in diverse classes of chemical reactions. The development of these objects is paramount in the theory of reaction dynamics as the transition state structure and its underlying network of manifolds directly dictate reactivity and selectivity.

Lagrangian descriptors in dissipative systems

The reaction dynamics of time-dependent systems can be resolved through a recrossing-free dividing surface associated with the transition state trajectory-that is, the unique trajectory which is bound to the barrier region for all time in response to a given time-dependent potential. A general procedure based on the minimization of Lagrangian descriptors has recently been developed by Craven and Hernandez [Phys. Rev. Lett., 2015, 115, 148301] to construct this particular trajectory without requiring perturbative expansions relative to the naive transition state point at the top of the barrier. The extension of the method to account for dissipation in the equations of motion requires additional considerations established in this paper because the calculation of the Lagrangian descriptor involves the integration of trajectories in forward and backward time. The two contributions are in general very different because the friction term can act as a source (in backward time) or sink (in forward time) of energy, leading to the possibility that information about the phase space structure may be lost due to the dominance of only one of the terms. To compensate for this effect, we introduce a weighting scheme within the Lagrangian descriptor and demonstrate that for thermal Langevin dynamics it preserves the essential phase space structures, while they are lost in the nonweighted case.

Performance of Lagrangian descriptors and their variants in incompressible flows

The method of Lagrangian Descriptors has been applied in many different contexts, specially in geophysical flows. In this paper, we analyze their performance in incompressible flows. We construct broad families of systems where this diagnostic fails in the detection of barriers to transport. Another aim of this manuscript is to illustrate the same deficiencies in the recent diagnostic proposed by Craven and Hernández.

Transition state geometry of driven chemical reactions on time-dependent double-well potentials

The minimum contour in the forward Lagrangian descriptor overlaps the invariant manifold (in green) dividing reactant and product regions.