BZPEERReading, trust and collaboration for researchDiscover
Workspaces
Network
Opportunities
Account
Researcher profile
Stephen Wiggins contributes to research discovery and scholarly infrastructure.
Trust snapshot
Actions
Identity and collaboration
Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.
Log in to claimDirect collaboration
Claim this author entity first to unlock direct invitations.
Research graph
Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph 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.
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 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.
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.