Computing One-Dimensional Global Manifolds of Poincaré Maps by Continuation

Abstract

We present an algorithm to compute one-dimensional stable and unstable manifolds of saddle periodic orbits in a Poincaré section. The computation is set up as a boundary value problem by restricting the beginning and end points of orbit segments to the section. Starting from the periodic orbit itself, we use collocation routines from AUTO to continue the solutions of the boundary value problem such that one end point of the orbit segment varies along a part of the manifold that was already computed. In this way, the other end point of the orbit segment traces out a new piece of the manifold. As opposed to standard methods that use shooting to compute the Poincaré map as the k-th return map, our approach defines the Poincaré map as the solution to a boundary value problem. This enables us to compute global manifolds through points where the flow is tangent to the section-a situation that is typically encountered unless one is dealing with a periodically forced system. Another major advantage of our approach is that it deals effectively with the problem of extreme sensitivity of the Poincaré map on its argument, which is a typical feature in the important class of slow-fast systems. We illustrate and test our algorithm by computing stable and unstable manifold in three examples: the forced Van der Pol oscillator, a model of a semiconductor laser with optical injection, and a slow-fast chemical oscillator. All examples are accompanied by animations of how the manifolds grow during the computation.

Solving Winfree's puzzle: The isochrons in the FitzHugh-Nagumo model

“…This approach is quite general [1,20], and its significant advantage is that each isochron in [25] is found directly as a parametrised curve. As well as producing a good mesh resolution for the isochron, their approach is not affected by areas of extreme phase sensitivity that occur in systems of multiple time-scales; see also [8].…”

Mixed-Mode Oscillations with Multiple Time Scales

“…Despite its simplicity, this sweeping method fails to produce satisfactory results in some cases. In particular, strong convergence or divergence of trajectories toward one another makes the choice of the initial mesh problematic and can produce very nonuniform "coverage" of the desired manifold; see [61,62]. In multiple-time-scale systems, the fast exponential instability of Fenichel manifolds that are not attracting makes initial value solvers incapable of tracking these manifolds by forward integration.…”

High-Order Bisection Method for Computing Invariant Manifolds of Two-Dimensional Maps

“…Similar methods for computing the unstable manifold of planar maps are developed in [England et al, 2005;Krauskopf & Osinga, 1998a, 1998bKrauskopf et al, 2004]. These are based on similar conditions to Eqs.…”