DSNov 4, 2017Code
Computational Method for Phase Space Transport with Applications to Lobe Dynamics and Rate of EscapeShibabrat Naik, Francois Lekien, Shane D. Ross
Lobe dynamics and escape from a potential well are general frameworks introduced to study phase space transport in chaotic dynamical systems. While the former approach studies how regions of phase space are transported by reducing the flow to a two-dimensional map, the latter approach studies the phase space structures that lead to critical events by crossing periodic orbit around saddles. Both of these frameworks require computation with curves represented by millions of points-computing intersection points between these curves and area bounded by the segments of these curves-for quantifying the transport and escape rate. We present a theory for computing these intersection points and the area bounded between the segments of these curves based on a classification of the intersection points using equivalence class. We also present an alternate theory for curves with nontransverse intersections and a method to increase the density of points on the curves for locating the intersection points accurately.The numerical implementation of the theory presented herein is available as an open source software called Lober. We used this package to demonstrate the application of the theory to lobe dynamics that arises in fluid mechanics, and rate of escape from a potential well that arises in ship dynamics.
CHEM-PHJul 21, 2021
Predicting trajectory behaviour via machine-learned invariant manifoldsVladimír Krajňák, Shibabrat Naik, Stephen Wiggins
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.
DSJul 18, 2021
Support vector machines for learning reactive islandsShibabrat Naik, Vladimír Krajňák, Stephen Wiggins
We develop a machine learning framework that can be applied to data sets derived from the trajectories of Hamilton's equations. The goal is to learn the phase space structures that play the governing role for phase space transport relevant to particular applications. Our focus is on learning reactive islands in two degrees-of-freedom Hamiltonian systems. Reactive islands are constructed from the stable and unstable manifolds of unstable periodic orbits and play the role of quantifying transition dynamics. We show that support vector machines (SVM) is an appropriate machine learning framework for this purpose as it provides an approach for finding the boundaries between qualitatively distinct dynamical behaviors, which is in the spirit of the phase space transport framework. We show how our method allows us to find reactive islands directly in the sense that we do not have to first compute unstable periodic orbits and their stable and unstable manifolds. We apply our approach to the Hénon-Heiles Hamiltonian system, which is a benchmark system in the dynamical systems community. We discuss different sampling and learning approaches and their advantages and disadvantages.