BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
Stephen Wiggins appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
This is the preface to a special issue of the Journal of Physical Organic Chemistry dedicated to an outstanding physical organic chemist, mentor and friend, Barry Carpenter on the occasion of his official retirement.
In this paper, we use support vector machines (SVM) to develop a machine learning framework to discover phase space structures that distinguish between distinct reaction pathways. The SVM model is trained using data from trajectories of Hamilton's equations and works well even with relatively few trajectories. Moreover, this framework is specifically designed to require minimal a priori knowledge of the dynamics in a system. This makes our approach computationally better suited than existing methods for high-dimensional systems and systems where integrating trajectories is expensive. We benchmark our approach on Chesnavich's CH$_4^+$ Hamiltonian.
In this paper we study an asymmetric valley-ridge inflection point (VRI) potential, whose energy surface (PES) features two sequential index-1 saddles (the upper and the lower), with one saddle having higher energy than the other and two potential wells separated by the lower index-1 saddle. We show how the depth and the flatness of our potential changes as we modify the parameter that controls the asymmetry as well as how the branching ratio (ratio of the trajectories that enter each well) is changing as we modify the same parameter and its correlation with the area of the lobes as they have formulated by the stable and unstable manifolds that have been extracted from the gradient of the LD scalar fields.
We used for the first time the method of periodic orbit dividing surfaces in a non-integrable Hamiltonian system with three degrees of freedom. We have studied the structure of these four dimensional objects in the five dimensional phase space. This method enabled us to detect the reactive and non-reactive trajectories in a three dimensional Caldera potential energy surface. We distinguished four distinct types of trajectory behavior. Two of the types of trajectories could only occur in a three dimensional Caldera potential energy surface and not in a two dimensional surface, and we have shown that this is a result of homoclinic intersections. These homoclinic intersections were analyzed with the method of Lagrangian descriptors. Finally, we were able to detect and describe the phenomenon of dynamical matching in the three dimensional Caldera potential energy surface, which is an important mechanism for understanding the reaction dynamics of organic molecules.
In this article, we review the analytical and numerical approaches for computing the phase space structures in two degrees-of-freedom Hamiltonian systems that arise in chemical reactions. In particular, these phase space structures are the unstable periodic orbit associated with an index-1 saddle, the periodic orbit dividing surface, and the stable and unstable invariant manifolds of the unstable periodic orbit. We review the approaches in the context of a two degrees-of-freedom Hamiltonian with a quartic potential coupled with a quadratic potential. We derive the analytical form of the phase space structures for the integrable case and visualize their geometry on the three dimensional energy surface. We then investigate the bifurcation of the dividing surface due to the changes in the parameters of the potential energy. We also review the numerical method of \emph{turning point} and present its new modification called the \emph{turning point based on configuration difference} for computing the unstable periodic orbit in two degrees-of-freedom systems. These methods are implemented in the open-source python package, UPOsHam~\cite{Lyu2020}.
In this paper we explore the phase space structures governing isomerization dynamics on a potential energy surface with four wells and an index-2 saddle. For this model, we analyze the influence that coupling both degrees of freedom of the system and breaking the symmetry of the problem have on the geometrical template of phase space structures that characterizes reaction. To achieve this goal we apply the method of Lagrangian descriptors, a technique with the capability of unveiling the key invariant manifolds that determine transport processes in nonlinear dynamical systems. This approach reveals with extraordinary detail the intricate geometry of the isomerization routes interconnecting the different potential wells, and provides us with valuable information to distinguish between initial conditions that undergo sequential and concerted isomerization.
We study the phase space structures that control the transport in a classical Hamiltonian model for a chemical reaction. This model has been proposed to study the yield of products in an ultracold exothermic reaction. In the considered model, two elements determine the evolution of the system: a Van der Waals force and short-range force associated with the many-body interactions. In the previous work has been used small random periodic changes in the direction of the momentum to simulate the short-range many-body interactions. In the present work, random Gaussian bumps have been added to the Van der Waals potential energy simulate the short-range effects between the particles in the system. We compare both variants of the model and explain their differences similarities and differences from a phase space perspective. In order to visualize the structures that direct the dynamics in the phase space, we construct a natural Lagrangian descriptor for Hamiltonian systems based on the Maupertuis action $S_0$.
In this paper, we analyse the phase space structure of the roaming dynamics in a two degree of freedom potential energy surface consisting of two identical planar Morse potentials separated by a distance. This potential energy surface was previously studied in [1], and it has two potential wells surrounded by an unbounded flat region containing no critical points. We study the phase space mechanism for the transference between the wells using the method of Lagrangian descriptors.
In light of the rapid recent retreat of Arctic sea ice, the extreme weather events triggering the variability in Arctic ice cover has drawn increasing attention. A non-Gaussian $α$-stable Lévy process is thought to be an appropriate model to describe such extreme event. The maximal likely trajectory, based on the nonlocal Fokker-Planck equation, is applied to a nonautonomous Arctic sea ice system under $α$-stable Lévy noise. Two types of tipping times, the early-warning tipping time and the disaster-happening tipping time, are used to predict the critical time for the maximal likely transition from a perennially ice-covered state to a seasonally ice-free one, and from a seasonally ice-free state to a perennially ice-free one, respectively. We find that the increased intensity of extreme events results in shorter warning time for sea ice melting, and that an enhanced greenhouse effect will intensify this influence, making the arrival of warning time significantly earlier. Meanwhile, for the enhanced greenhouse effect, we discover that increased intensity and frequency of extreme events will advance the disaster-happening tipping time, in which an ice-free state is maintained throughout the year in the Arctic Ocean. Finally, we identify values of Lévy index $α$ and noise intensity $ε$ in $αε$-space that can trigger a transition between the Arctic sea ice state. These results provide an effective theoretical framework for studying Arctic sea ice variations under the influence of extreme events.
In this article we present the influence of a Hamiltonian saddle-node bifurcation on the high-dimensional phase space structures that mediate reaction dynamics. To achieve this goal, we identify the phase space invariant manifolds using Lagrangian descriptors, which is a trajectory-based diagnostic suitable for the construction of a complete ``phase space tomography'' by means of analyzing dynamics on low-dimensional slices. First, we build a Hamiltonian system with one degree-of-freedom (DoF) that models reaction, and study the effect of adding a parameter to the potential energy function that controls the depth of the well. Then, we extend this framework to a saddle-node bifurcation for a two DoF Hamiltonian, constructed by coupling a harmonic oscillator, i.e. a bath mode, to the other reactive DoF in the system. For this problem, we describe the phase space structures associated with the rank-1 saddle equilibrium point in the bottleneck region, which is a Normally Hyperbolic Invariant Manifold (NHIM) and its stable and unstable manifolds. Finally, we address the qualitative changes in the reaction dynamics of the Hamiltonian system due to changes in the well depth of the potential energy surface that gives rise to the saddle-node bifurcation.
The paper deals with the problem of existence of a convergent "strong" normal form in the neighbourhood of an equilibrium, for a finite dimensional system of differential equations with analytic and time-dependent non-linear term. The problem can be solved either under some non-resonance hypotheses on the spectrum of the linear part or if the non-linear term is assumed to be (slowly) decaying in time. This paper "completes" a pioneering work of Pustil'nikov in which, despite under weaker non-resonance hypotheses, the nonlinearity is required to be asymptotically autonomous. The result is obtained as a consequence of the existence of a strong normal form for a suitable class of real-analytic Hamiltonians with non-autonomous perturbations.
The paper deals with the problem of the existence of a normal form for a nearly-integrable real-analytic Hamiltonian with aperiodically time-dependent perturbation decaying (slowly) in time. In particular, in the case of an isochronous integrable part, the system can be cast in an exact normal form, regardless of the properties of the frequency vector. The general case is treated by a suitable adaptation of the finite order normalization techniques usually used for Nekhoroshev arguments. The key point is that the so called "geometric part" is not necessary in this case. As a consequence, no hypotheses on the integrable part are required, apart from analyticity. The work, based on two different perturbative approaches developed by A.Giorgilli et al., is a generalisation of the techniques used by the same authors to treat more specific aperiodically time-dependent problems.
The classical theorem of Moser, on the existence of a normal form in the neighbourhood of a hyperbolic equilibrium, is extended to a class of real-analytic Hamiltonians with aperiodically time-dependent perturbations. A stronger result is obtained in the case in which the perturbing function exhibits a time decay.
We examine the phase space structures that govern reaction dynamics in the absence of critical points on the potential energy surface. We show that in the vicinity of hyperbolic invariant tori it is possible to define phase space dividing surfaces that are analogous to the dividing surfaces governing transition from reactants to products near a critical point of the potential energy surface. We investigate the problem of capture of an atom by a diatomic molecule and show that a normally hyperbolic invariant manifold exists at large atom-diatom distances, away from any critical points on the potential. This normally hyperbolic invariant manifold is the anchor for the construction of a dividing surface in phase space, which defines the outer or loose transition state governing capture dynamics. We present an algorithm for sampling an approximate capture dividing surface, and apply our methods to the recombination of the ozone molecule. We treat both 2 and 3 degree of freedom models with zero total angular momentum. We have located the normally hyperbolic invariant manifold from which the orbiting (outer) transition state is constructed. This forms the basis for our analysis of trajectories for ozone in general, but with particular emphasis on the roaming trajectories.
The aim of this paper is to prove the Kolmogorov theorem of persistence of Diophantine flows for nearly-integrable Poisson systems associated to a real analytic Hamiltonian with aperiodic time dependence, provided that the perturbation is asymptotically vanishing. The paper is an extension of an analogous result by the same authors for canonical Hamiltonian systems; the flexibility of the Lie series method developed by A. Giorgilli et al., is profitably used in the present generalisation.
We explore both classical and quantum dynamics of a model potential exhibiting a caldera: that is, a shallow potential well with two pairs of symmetry related index one saddles associated with entrance/exit channels. Classical trajectory simulations at several different energies confirm the existence of the `dynamical matching' phenomenon originally proposed by Carpenter, where the momentum direction associated with an incoming trajectory initiated at a high energy saddle point determines to a considerable extent the outcome of the reaction (passage through the diametrically opposing exit channel). By studying a `stretched' version of the caldera model, we have uncovered a generalized dynamical matching: bundles of trajectories can reflect off a hard potential wall so as to end up exiting predominantly through the transition state opposite the reflection point. We also investigate the effects of dissipation on the classical dynamics. In addition to classical trajectory studies, we examine the dynamics of quantum wave packets on the caldera potential (stretched and unstretched). These computations reveal a quantum mechanical analogue of the `dynamical matching' phenomenon, where the initial expectation value of the momentum direction for the wave packet determines the exit channel through which most of the probability density passes to product.
The aim of this paper is to prove a Kolmogorov-type result for a nearly-integrable Hamiltonian, quadratic in the actions, with an aperiodic time dependence. The existence of a torus with a prefixed Diophantine frequency is shown in the forced system, provided that the perturbation is real-analytic and (exponentially) decaying with time. The advantage consists of the possibility to choose an arbitrarily small decaying coefficient, consistently with the perturbation size.
A reduced two dimensional model is used to study Ketene isomerization reaction. In light of recent results by Ulusoy \textit{et al.} [J.\ Phys.\ Chem.\ A {\bf 117}, 7553 (2013)], the present work focuses on the generalization of the roaming mechanism to the Ketene isomerization reaction by applying our phase space approach previously used to elucidate the roaming phenomenon in ion-molecule reactions. Roaming is again found be associated with the trapping of trajectories in a phase space region between two dividing surfaces; trajectories are classified as reactive or nonreactive, and are further naturally classified as direct or non-direct (roaming). The latter long-lived trajectories are trapped in the region of non-linear mechanical resonances, which in turn define alternative reaction pathways in phase space. It is demonstrated that resonances associated with periodic orbits provide a dynamical explanation of the quantum mechanical resonances found in the isomerization rate constant calculations by Gezelter and Miller [J.\ Chem.\ Phys.\ {\bf 103}, 7868-7876 (1995)]. Evidence of the trapping of trajectories by `sticky' resonant periodic orbits is provided by plotting Poincaré surfaces of section, and a gap time analysis is carried out in order to investigate the statistical assumption inherent in transition state theory for Ketene isomerization.
We develop a framework to study the role of variability in transport across a streamline of a reference flow. Two complementary schemes are presented: a graphical approach for individual cases, and an analytical approach for general properties. The spatially nonlinear interaction of dynamic variability and the reference flow results in flux variability. The characteristic time-scale of the dynamic variability and the length-scale of the flux variability in a unit of flight-time govern the spatio-temporal interaction that leads to transport. The non-dimensional ratio of the two characteristic scales is shown to be a a critical parameter. The pseudo-lobe sequence along the reference streamline describes spatial coherency and temporal evolution of transport. For finite-time transport from an initial time up to the present, the characteristic length-scale of the flux variability regulates the width of the pseudo-lobes. The phase speed of pseudo-lobe propagation averages the reference flow and the flux variability. In contrast, for definite transport over a fixed time interval and spatial segment, the characteristic time-scale of the dynamic variability regulates the width of the pseudo-lobes. Generation of the pseudo-lobe sequence appears to be synchronous with the dynamic variability, although it propagates with the reference flow. In either case, the critical characteristic ratio is found to be one, corresponding to a resonance of the flux variability with the reference flow. Using a kinematic model, we demonstrate the framework for two types of transport in a blocked flow of the mid-latitude atmosphere: across the meandering jet axis and between the jet and recirculating cell.
In this paper we develop new techniques for revealing geometrical structures in phase space that are valid for aperiodically time dependent dynamical systems, which we refer to as Lagrangian descriptors. These quantities are based on the integration, for a finite time, along trajectories of an intrinsic bounded, positive geometrical and/or physical property of the trajectory itself. We discuss a general methodology for constructing Lagrangian descriptors, and we discuss a "heuristic argument" that explains why this method is successful for revealing geometrical structures in the phase space of a dynamical system. We support this argument by explicit calculations on a benchmark problem having a hyperbolic fixed point with stable and unstable manifolds that are known analytically. Several other benchmark examples are considered that allow us the assess the performance of Lagrangian descriptors in revealing invariant tori and regions of shear. Throughout the paper "side-by-side" comparisons of the performance of Lagrangian descriptors with both finite time Lyapunov exponents (FTLEs) and finite time averages of certain components of the vector field ("time averages") are carried out and discussed. In all cases Lagrangian descriptors are shown to be both more accurate and computationally efficient than these methods. We also perform computations for an explicitly three dimensional, aperiodically time-dependent vector field and an aperiodically time dependent vector field defined as a data set. Comparisons with FTLEs and time averages for these examples are also carried out, with similar conclusions as for the benchmark examples.
We provide a dynamical interpretation of the recently identified `roaming' mechanism for molecular dissociation reactions in terms of geometrical structures in phase space. These are NHIMs (Normally Hyperbolic Invariant Manifolds) and their stable/unstable manifolds that define transition states for ion-molecule association or dissociation reactions. The associated dividing surfaces rigorously define a roaming region of phase space, in which both reactive and nonreactive trajectories can be trapped for arbitrarily long times.
We study reaction dynamics on a model potential energy surface exhibiting post-transition state bifurcation in the vicinity of a valley ridge inflection point. We compute fractional yields of products reached after the VRI region is traversed, both with and without dissipation. It is found that apparently minor variations in the potential lead to significant changes in the reaction dynamics. Moreover, when dissipative effects are incorporated, the product ratio depends in a complicated and highly non-monotonic fashion on the dissipation parameter. Dynamics in the vicinity of the VRI point itself play essentially no role in determining the product ratio, except in the highly dissipative regime.
The aim of this paper is to extend the results of Giorgilli and Zehnder for aperiodic time dependent systems to a case of general nearly-integrable convex analytic Hamiltonians. The existence of a normal form and then a stability result are shown in the case of a slow aperiodic time dependence that, under some smallness conditions, is independent on the size of the perturbation.
We analyze the dynamics of a model of a nanobeam under compression. The model is a two mode truncation of the Euler-Bernoulli beam equation subject to compressive stress. We consider parameter regimes where the first mode is unstable and the second mode can be either stable or unstable, and the remaining modes (neglected) are always stable. Material parameters used correspond to silicon. The two mode model Hamiltonian is the sum of a (diagonal) kinetic energy term and a potential energy term. The form of the potential energy function suggests an analogy with isomerisation reactions in chemistry. We therefore study the dynamics of the buckled beam using the conceptual framework established for the theory of isomerisation reactions. When the second mode is stable the potential energy surface has an index one saddle and when the second mode is unstable the potential energy surface has an index two saddle and two index one saddles. Symmetry of the system allows us to construct a phase space dividing surface between the two "isomers" (buckled states). The energy range is sufficiently wide that we can treat the effects of the index one and index two saddles in a unified fashion. We have computed reactive fluxes, mean gap times and reactant phase space volumes for three stress values at several different energies. In all cases the phase space volume swept out by isomerizing trajectories is considerably less than the reactant density of states, proving that the dynamics is highly nonergodic. The associated gap time distributions consist of one or more `pulses' of trajectories. Computation of the reactive flux correlation function shows no sign of a plateau region; rather, the flux exhibits oscillatory decay, indicating that, for the 2-mode model in the physical regime considered, a rate constant for isomerization does not exist.
In this paper we continue our studies of the phase space geometry and dynamics associated with index k saddles (k > 1) of the potential energy surface. Using normal form theory, we give an explicit formula for a "dividing surface" in phase space, i.e. a co-dimension one surface (within the energy shell) through which all trajectories that "cross" the region of the index k saddle must pass. With a generic non-resonance assumption, the normal form provides k (approximate) integrals that describe the saddle dynamics in a neighborhood of the index k saddle. These integrals provide a symbolic description of all trajectories that pass through a neighborhood of the saddle. We give a parametrization of the dividing surface which is used as the basis for a numerical method to sample the dividing surface. Our techniques are applied to isomerization dynamics on a potential energy surface having 4 minima; two symmetry related pairs of minima are connected by low energy index one saddles, with the pairs themselves connected via higher energy index one saddles and an index two saddle at the origin. We compute and sample the dividing surface and show that our approach enables us to distinguish between concerted crossing ("hilltop crossing") isomerizing trajectories and those trajectories that are not concerted crossing (potentially sequentially isomerizing trajectories). We then consider the effect of additional "bath modes" on the dynamics, which is a four degree-of-freedom system. For this system we show that the normal form and dividing surface can be realized and sampled and that, using the approximate integrals of motion and our symbolic description of trajectories, we are able to choose initial conditions corresponding to concerted crossing isomerizing trajectories and (potentially) sequentially isomerizing trajectories.
In this paper we develop a method for the estimation of {\bf T}ransport {\bf I}nduced by the {\bf M}ean-{\bf E}ddy interaction (TIME) in two-dimensional unsteady flows. The method is built on the dynamical systems approach and can be viewed as a hybrid combination of Lagrangian and Eulerian methods. The (Eulerian) boundaries across which we consider (Lagrangian) transport are kinematically defined by appropriately chosen streamlines of the mean flow. By evaluating the impact of the mean-eddy interaction on transport, the TIME method can be used as a diagnostic tool for transport processes that occur during a specified time interval along a specified boundary segment. We introduce two types of TIME functions: one that quantifies the accumulation of flow properties and another that measures the displacement of the transport geometry. The spatial geometry of transport is described by the so-called pseudo-lobes, and temporal evolution of transport by their dynamics. In the case where the TIME functions are evaluated along a separatrix, the pseudo-lobes have a relationship to the lobes of Lagrangian transport theory. In fact, one of the TIME functions is identical to the Melnikov function that is used to measure the distance, at leading order in a small parameter, between the two invariant manifolds that define the Lagrangian lobes. We contrast the similarities and differences between the TIME and Lagrangian lobe dynamics in detail. An application of the TIME method is carried out for inter-gyre transport in the wind-driven oceanic circulation model and a comparison with the Lagrangian transport theory is made.
We present a framework for the analysis of transport processes resulting from the mean-eddy interaction in a flow. The framework is based on the {\bf T}ransport {\bf I}nduced by the {\bf M}ean-{\bf E}ddy {\bf I}nteraction (TIME) method presented in a companion paper \cite{ide_wiggins_pd06a}. The TIME method estimates the (Lagrangian) transport across stationary (Eulerian) boundaries defined by chosen streamlines of the mean flow. Our framework proceeds after first carrying out a sequence of preparatory steps that link the flow dynamics to the transport processes. This includes the construction of the so-called "instantaneous flux" as the Hovmöller diagram. Transport processes are studied by linking the signals of the instantaneous flux field to the dynamical variability of the flow. This linkage also reveals how the variability of the flow contributes to the transport. The spatio-temporal analysis of the flux diagram can be used to assess the efficiency of the variability in transport processes. We apply the method to the double-gyre ocean circulation model in the situation where the Rossby-wave mode dominates the dynamic variability. The spatio-temporal analysis shows that the inter-gyre transport is controlled by the circulating eddy vortices in the fast eastward jet region, whereas the basin-scale Rossby waves have very little impact.
Transition State Theory forms the basis of computing reaction rates in chemical and other systems. Recently it has been shown how transition state theory can rigorously be realized in phase space using an explicit algorithm. The quantization has been demonstrated to lead to an efficient procedure to compute cumulative reaction probabilities and the associated Gamov-Siegert resonances. In this letter these results are used to express the cumulative reaction probability as an absolutely convergent sum over periodic orbits contained in the transition state.
We consider the existence of invariant manifolds in phase space governing reaction dynamics in situations where there are no saddle points on the potential energy surface in the relevant regions of configuration space. We point out that such situations occur in a number of important classes of chemical reactions, and we illustrate this concretely by considering a model for transition state switching in an ion-molecule association reaction due to Chesnavich (J. Chem. Phys. {\bf 84}, 2615 (1986)). For this model we show that, in the region of configuration space relevant to the reaction, there are no saddle points on the potential energy surface, but that in phase space there is a normally hyperbolic invariant manifold (NHIM) bounding a dividing surface having the property that the reactive flux through this dividing surface is a minimum. We then describe two methods for finding NHIMs and their associated phase space structures in systems with more than two degrees-of-freedom. These methods do not rely on the existence of saddle points, or any other particular feature, of the potential energy surface.
We investigate the phase space structure and dynamics of a Hamiltonian isokinetic thermostat, for which ergodic thermostat trajectories at fixed (zero) energy generate a canonical distribution in configuration space. Model potentials studied consist of a single bistable mode plus transverse harmonic modes. Interpreting the bistable mode as a reaction (isomerization) coordinate, we establish connections with the theory of unimolecular reaction rates, in particular the formulation of isomerization rates in terms of gap times. The distribution of gap times (or associated lifetimes) for a microcanonical ensemble initiated on the dividing surface is of great dynamical significance; an exponential lifetime distribution is usually taken to be an indicator of `statistical' behavior. Moreover, comparison of the magnitude of the phase space volume swept out by reactive trajectories as they pass through the reactant region with the total phase space volume (classical density of states) for the reactant region provides a necessary condition for ergodic dynamics. We compute gap times, associated lifetime distributions, mean gap times, reactive fluxes, reactive volumes and total reactant phase space volumes for model systems with 3 degrees of freedom, at three different temperatures. At all three temperatures, the necessary condition for ergodicity is approximately satisfied. At high temperatures a non-exponential lifetime distribution is found, while at low temperatures the lifetime is more nearly exponential. The degree of exponentiality of the lifetime distribution is quantified by computing the information entropy deficit with respect to pure exponential decay. The efficacy of the Hamiltonian isokinetic thermostat is examined by computing coordinate distributions averaged over single long trajectories initiated on the dividing surface.
We review recent results on quantum reactive scattering from a phase space perspective. The approach uses classical and quantum versions of normal form theory and the perspective of dynamical systems theory. Over the past ten years the classical normal form theory has provided a method for realizing the phase space structures that are responsible for determining reactions in high dimensional Hamiltonian systems. This has led to the understanding that a new (to reaction dynamics) type of phase space structure, a {\em normally hyperbolic invariant manifold} (or, NHIM) is the "anchor" on which the phase space structures governing reaction dynamics are built. The quantum normal form theory provides a method for quantizing these phase space structures through the use of the Weyl quantization procedure. We show that this approach provides a solution of the time-independent Schrödinger equation leading to a (local) S-matrix in a neighborhood of the saddle point governing the reaction. It follows easily that the quantization of the directional flux through the dividing surface with the properties noted above is a flux operator that can be expressed in a "closed form". Moreover, from the local S-matrix we easily obtain an expression for the cumulative reactio probability (CRP). Significantly, the expression for the CRP can be evaluated without the need to compute classical trajectories. The quantization of the NHIM is shown to lead to the activated complex, and the lifetimes of quantum states initialized on the NHIM correspond to the Gamov-Siegert resonances. We apply these results to the collinear nitrogen exchange reaction and a three degree-of-freedom system corresponding to an Eckart barrier coupled to two Morse oscillators.
The flux-flux correlation function formalism is a standard and widely used approach for the computation of reaction rates. In this paper we introduce a method to compute the classical and quantum flux-flux correlation functions for anharmonic barriers essentially analytically through the use of the classical and quantum normal forms. In the quantum case we show that the quantum normal form reduces the computation of the flux-flux correlation function to that of an effective one dimensional anharmonic barrier. The example of the computation of the quantum flux-flux correlation function for a fourth order anharmonic barrier is worked out in detail, and we present an analytical expression for the quantum mechanical microcanonical flux-flux correlation function. We then give a discussion of the short-time and harmonic limits.
We consider issues associated with the Lagrangian characterisation of flow structures arising in aperiodically time-dependent vector fields that are only known on a finite time interval. A major motivation for the consideration of this problem arises from the desire to study transport and mixing problems in geophysical flows where the flow is obtained from a numerical solution, on a finite space-time grid, of an appropriate partial differential equation model for the velocity field. Of particular interest is the characterisation, location, and evolution of "transport barriers" in the flow, i.e. material curves and surfaces. We argue that a general theory of Lagrangian transport has to account for the effects of transient flow phenomena which are not captured by the infinite-time notions of hyperbolicity even for flows defined for all time. Notions of finite-time hyperbolic trajectories, their finite time stable and unstable manifolds, as well as finite-time Lyapunov exponent (FTLE) fields and associated Lagrangian coherent structures have been the main tools for characterizing transport barriers in the time-aperiodic situation. In this paper we consider a variety of examples, some with explicit solutions, that illustrate, in a concrete manner, the issues and phenomena that arise in the setting of finite-time dynamical systems. Of particular significance for geophysical applications is the notion of "flow transition" which occurs when finite-time hyperbolicity is lost, or gained. The phenomena discovered and analysed in our examples point the way to a variety of directions for rigorous mathematical research in this rapidly developing, and important, new area of dynamical systems theory.
Hamiltonian dynamical systems possessing equilibria of ${saddle} \times {centre} \times...\times {centre}$ stability type display \emph{reaction-type dynamics} for energies close to the energy of such equilibria; entrance and exit from certain regions of the phase space is only possible via narrow \emph{bottlenecks} created by the influence of the equilibrium points. In this paper we provide a thorough pedagogical description of the phase space structures that are responsible for controlling transport in these problems. Of central importance is the existence of a \emph{Normally Hyperbolic Invariant Manifold (NHIM)}, whose \emph{stable and unstable manifolds} have sufficient dimensionality to act as separatrices, partitioning energy surfaces into regions of qualitatively distinct behavior. This NHIM forms the natural (dynamical) equator of a (spherical) \emph{dividing surface} which locally divides an energy surface into two components (`reactants' and `products'), one on either side of the bottleneck. This dividing surface has all the desired properties sought for in \emph{transition state theory} where reaction rates are computed from the flux through a dividing surface. In fact, the dividing surface that we construct is crossed exactly once by reactive trajectories, and not crossed by nonreactive trajectories, and related to these properties, minimizes the flux upon variation of the dividing surface. We discuss three presentations of the energy surface and the phase space structures contained in it for 2-degree-of-freedom (DoF) systems in the threedimensional space $\R^3$, and two schematic models which capture many of the essential features of the dynamics for $n$-DoF systems. In addition, we elucidate the structure of the NHIM.
The general approach to classical unimolecular reaction rates due to Thiele is revisited in light of recent advances in the phase space formulation of transition state theory for multidimensional systems. We analyze in detail the gap time distribution and associated reactant lifetime distribution for the isomerization reaction HCN $\rightleftharpoons$ CNH. Both algebraic (power law) and exponential decay regimes have been identified. Statistical estimates of the isomerization rate are compared with the numerically determined decay rate. Examination of the decay properties of subsensembles of trajectories that exit the HCN well through either of 2 available symmetry related product channels shows that the complete trajectory ensemble effectively attains the full symmetry of the system phase space on a short timescale $t \lesssim 0.5$ ps, after which the product branching ratio is 1:1, the "statistical" value. At intermediate times, this statistical product ratio is accompanied by nonexponential (nonstatistical) decay. We point out close parallels between the dynamical behavior inferred from the gap time distribution for HCN and nonstatistical behavior recently identified in reactions of some organic molecules.
We study the phase space geometry associated with index 2 saddles of a potential energy surface and its influence on reaction dynamics for $n$ degree-of-freedom (DoF) Hamiltonian systems. For index 1 saddles of potential energy surfaces (the case of classical transition state theory), the existence of a normally hyperbolic invariant manifold (NHIM) of saddle stability type has been shown, where the NHIM serves as the "anchor" for the construction of dividing surfaces having the no-recrossing property and minimal flux. For the index 1 saddle case the stable and unstable manifolds of the NHIM are co-dimension one in the energy surface, and act as conduits for reacting trajectories in phase space. The situation for index 2 saddles is quite different. We show that NHIMs with their stable and unstable manifolds still exist, but that these manifolds by themselves lack sufficient dimension to act as barriers in the energy surface. Rather, there are different types of invariant manifolds, containing the NHIM and its stable and unstable manifolds, that act as co-dimension one barriers in the energy surface. These barriers divide the energy surface in the vicinity of the index 2 saddle into regions of qualitatively different trajectories exhibiting a wider variety of dynamical behavior than for the case of index 1 saddles. In particular, we can identify a class of trajectories, which we refer to as "roaming trajectories", which are not associated with reaction along the classical minimum energy path (MEP). We illustrate the significance of our analysis of the index 2 saddle for reaction dynamics with two examples.
The quantum normal form approach to quantum transition state theory is used to compute the cumulative reaction probability for collinear exchange reactions. It is shown that for heavy atom systems like the nitrogen exchange reaction the quantum normal form approach gives excellent results and has major computational benefits over full reactive scattering approaches. For light atom systems like the hydrogen exchange reaction however the quantum normal approach is shown to give only poor results. This failure is attributed to the importance of tunnelling trajectories in light atom reactions that are not captured by the quantum normal form as indicated by the only very slow convergence of the quantum normal form for such systems.