An iterative method for the canard explosion in general planar systems

The canard explosion is the change of amplitude and period of a limit cycle born in a Hopf bifurcation in a very narrow parameter interval. The phenomenon is well understood in singular perturbation problems where a small parameter controls the slow/fast dynamics. However, canard explosions are also observed in systems where no such parameter is present. Here we show how the iterative method of Roussel and Fraser, devised to construct regular slow manifolds, can be used to determine a canard point in a general planar system of nonlinear ODEs. We demonstrate the method on the van der Pol equation, showing that the asymptotics of the method is correct, and on a templator model for a self-replicating system.Comment: Paper presented at the 9th AIMS Conference on Dynamical Systems, Differential Equations and Applications, Orlando, Florida, USA July 1 - 5, 201

An iterative method for the approximation of fibers in slow-fast systems

In this paper we extend a method for iteratively improving slow manifolds so that it also can be used to approximate the fiber directions. The extended method is applied to general finite dimensional real analytic systems where we obtain exponential estimates of the tangent spaces to the fibers. The method is demonstrated on the Michaelis-Menten-Henri model and the Lindemann mechanism. The latter example also serves to demonstrate the method on a slow-fast system in non-standard slow-fast form. Finally, we extend the method further so that it also approximates the curvature of the fibers.Comment: To appear in SIAD

Canards in stiction: on solutions of a friction oscillator by regularization

We study the solutions of a friction oscillator subject to stiction. This discontinuous model is non-Filippov, and the concept of Filippov solution cannot be used. Furthermore some Carath\'eodory solutions are unphysical. Therefore we introduce the concept of stiction solutions: these are the Carath\'eodory solutions that are physically relevant, i.e. the ones that follow the stiction law. However, we find that some of the stiction solutions are forward non-unique in subregions of the slip onset. We call these solutions singular, in contrast to the regular stiction solutions that are forward unique. In order to further the understanding of the non-unique dynamics, we introduce a regularization of the model. This gives a singularly perturbed problem that captures the main features of the original discontinuous problem. We identify a repelling slow manifold that separates the forward slipping to forward sticking solutions, leading to a high sensitivity to the initial conditions. On this slow manifold we find canard trajectories, that have the physical interpretation of delaying the slip onset. We show with numerics that the regularized problem has a family of periodic orbits interacting with the canards. We observe that this family has a saddle stability and that it connects, in the rigid body limit, the two regular, slip-stick branches of the discontinuous problem, that were otherwise disconnected.Comment: Submitted to: SIADS. 28 pages, 12 figure

Mixed-mode oscillations in a multiple time scale phantom bursting system

In this work we study mixed mode oscillations in a model of secretion of GnRH (Gonadotropin Releasing Hormone). The model is a phantom burster consisting of two feedforward coupled FitzHugh-Nagumo systems, with three time scales. The forcing system (Regulator) evolves on the slowest scale and acts by moving the slow nullcline of the forced system (Secretor). There are three modes of dynamics: pulsatility (transient relaxation oscillation), surge (quasi steady state) and small oscillations related to the passage of the slow nullcline through a fold point of the fast nullcline. We derive a variety of reductions, taking advantage of the mentioned features of the system. We obtain two results; one on the local dynamics near the fold in the parameter regime corresponding to the presence of small oscillations and the other on the global dynamics, more specifically on the existence of an attracting limit cycle. Our local result is a rigorous characterization of small canards and sectors of rotation in the case of folded node with an additional time scale, a feature allowing for a clear geometric argument. The global result gives the existence of an attracting unique limit cycle, which, in some parameter regimes, remains attracting and unique even during passages through a canard explosion.Comment: 38 pages, 16 figure

Chaos at the border of criticality

The present paper points out to a novel scenario for formation of chaotic attractors in a class of models of excitable cell membranes near an Andronov-Hopf bifurcation (AHB). The mechanism underlying chaotic dynamics admits a simple and visual description in terms of the families of one-dimensional first-return maps, which are constructed using the combination of asymptotic and numerical techniques. The bifurcation structure of the continuous system (specifically, the proximity to a degenerate AHB) endows the Poincare map with distinct qualitative features such as unimodality and the presence of the boundary layer, where the map is strongly expanding. This structure of the map in turn explains the bifurcation scenarios in the continuous system including chaotic mixed-mode oscillations near the border between the regions of sub- and supercritical AHB. The proposed mechanism yields the statistical properties of the mixed-mode oscillations in this regime. The statistics predicted by the analysis of the Poincare map and those observed in the numerical experiments of the continuous system show a very good agreement.Comment: Chaos: An Interdisciplinary Journal of Nonlinear Science (tentatively, Sept 2008