Figure 1: A canard (with hat) of the van der Pol oscillator
Canards were discovered and first analyzed by French mathematicians (Benoît et al. [1981]) who studied 2D
relaxation oscillators, in particular the
van der Pol oscillator. There the
classical canard phenomenon explains the very fast transition upon variation of a parameter from a small amplitude
limit cycle via
canard cycles to a large amplitude
relaxation cycle. This very fast transition called
canard explosion happens within an exponentially small range of a control parameter. Thus this phenomenon is very hard to detect and it seems to be more like a
canard in a newspaper. Furthermore, the shape of a canard cycle in the
phase space resembles that of a duck (Figure 1). So the notion
canard was born and the chase on these creatures began with either nonstandard (Benoît et al [1981]) or standard (Eckhaus [1983]) methods.
[edit] Methods
The main mathematical methods to analyze canards
are nonstandard analysis (see, e.g., Benoît et al. [1981], Diener [1984]), matched asymptotic expansions (see, e.g., Eckhaus [1983], Mishchenko et al. [1994]) and the blow-up technique
(see, e.g., Dumortier and Roussarie [1996], Krupa and Szmolyan [2001],
Szmolyan and Wechselberger [2001]), which extends
geometric singular perturbation theory known as
Fenichel theory to non-hyperbolic
points. Moreover, with complex analysis, the canard phenomenon can be
understood with either a geometrical point of view (see, e.g. Callot
[1993] or Benoît et al. [1998]) or the Gevrey theory which is
convenient to study the divergent asymptotic expansions of canards
(see, e.g. Canalis-Durand et al. [2000] or Benoît [2001]).
[edit] Definition
Canards are a phenomenon occurring in singularly perturbed systems (also known as slow-fast systems), i.e. they occur in systems of the form
with fast variable slow variable sufficiently smooth functions and small parameter
This system evolves on a slow time scale The limiting problem on this slow time scale is called the reduced problem and describes the evolution of the slow variable The phase space of the reduced problem is the critical manifold defined by
On the other hand, the limiting problem on the fast time scale
is called the layer problem and describes the evolution of the fast variable
for fixed
By Fenichel theory (Fenichel [1979]), normally hyperbolic subsets
of perturb to nearby slow invariant manifolds
of the singularly perturbed system with the
flow given approximately by the flow of the reduced system. The most
common case where normal hyperbolicity breaks down is given by a
folded critical manifold where
denotes the attracting part of the critical manifold
the repelling part of
and the fold of along which normal
hyperbolicity is lost via a saddle-node bifurcation of the layer
problem. Of course, any bifurcation of the layer problem leads to a
loss of normal hyperbolicity of Canards are a special
class of solutions of singularly perturbed systems where normal
hyperbolicity is lost.
Figure 2: Examples of canard solutions (green) near critical manifold S; normal hyperbolicity is lost at L, e.g., via saddle-node bifurcation of layer problem (Fig. A), or via transcritical bifurcation of layer problem (Fig. B)
Definition: A canard is a solution of a singularly perturbed system which follows an attracting slow manifold passes close to a bifurcation point of the critical manifold, and then follows a repelling slow manifold for a considerable amount of time.
In geometric terms a canard solution corresponds to the intersection of an attracting and repelling slow manifold near a non-hyperbolic point This geometric object is called a maximal canard.
In the case of a Hopf bifurcation of the layer
problem, the phenomenon is quite different: there is a delayed loss of stability
but there exist no canards (see Wallet [1986], Neishtadt [1987]). This case will not
be explained here.
[edit] Codimension of the Canard Phenomenon
Canards in singularly perturbed systems with just one slow variable
() and one fast variable
() are non generic, since maximal canards in such systems occur only for discrete values of a control parameter (Krupa and Szmolyan [2001]), i.e. a one parameter family of singularly perturbed systems is needed to unfold this canard problem. In general, if canards are considered in singular perturbation problems with and where normal hyperbolicity is lost via a saddle-node bifurcation of the layer problem, then the singularly perturbed system can be reduced to a one slow and one fast variable system by a center manifold reduction.
[edit] The Classical Canard Phenomenon
[edit] Canard Cycles and their Explosion
The classical canard phenomenon (discovered by Benoît et al. [1981])
occurs in singularly perturbed systems with and Its prototypical example is the van der Pol oscillator given by
where denotes a control parameter (external forcing).
The critical manifold is cubic shaped consisting of two attracting outer branches and a middle repelling branch connected via two fold-points where normal hyperbolicity is lost.
Figure 3: Canard explosion of the
van der Pol oscillator () within an exponentially small neighbourhood of a=0.998740451245 where the transition from relaxation
oscillations to small amplitude limit cycles happens via canard cycles. The small amplitude limit cycles then terminate at a=1 via a Hopf bifurcation.
Here, the nature of the classical canard phenomenon is the transition from a small amplitude oscillatory state to a (large amplitude) relaxation oscillatory state within an exponentially small range of the control parameter
(see Figure 3). This transition, also called canard explosion, occurs through a sequence of canard cycles which can be asymptotically stable, but they are hard to observe in an experiment or simulation because of sensitivity to the control parameter and also because of sensitivity to noise. This is well known in chemical literature where a canard explosion is classified as a hard transition, because, for practical purposes, the transition from a small amplitude oscillation to a relaxation oscillation occurs immediately (see, e.g., Brøns and Bar-Eli [1991], Peng et al. [1991]).
[edit] Generic Canards
Canards in singularly perturbed systems with two or more slow variables () and one fast variable () are robust, since maximal canards generically persist
under small parameter changes (Benoît [1983], Mishchenko et
al. [1994], Szmolyan and Wechselberger [2001], Benoît [2001], Wechselberger [2005]). In general, if canards are considered in singular perturbation problems with and where normal hyperbolicity is lost via a saddle-node bifurcation of the layer problem, then the singularly perturbed system can be reduced to an slow and one fast variable system by a center manifold reduction. Therefore, a generic example with minimal dimension is given by a 3D singularly perturbed system with a 2D folded critical manifold.
Figure 4: A: The reduced flow near a saddle singularity, B-C: The corresponding reduced flow near a folded saddle singularity (B) and the same reduced flow shown on the 2D folded critical manifold (C)
[edit] Classification
There exist different types of canards in 3D systems with 2D folded critical manifolds (Benoît [1983], Szmolyan and Wechselberger [2001]). The classification of these canards is based on the analysis of the corresponding 2D reduced system. There exist (discrete) folded singularities on the fold called canard points, where the reduced flow crosses from the attracting to the repelling branch of the critical manifold. Generically, these folded singularities are either folded saddles (see Figure 4) or folded nodes (see Figure 5), resembling similar phase portraits as ordinary saddles or nodes in 2D phase space, besides that the reduced flow allow trajectories to cross the fold at these canard points. Note that in the folded saddle case (Figure 4B-4C) two solutions of the reduced flow cross via the folded saddle from the attracting branch to the repelling branch of the critical manifold or vice versa. More strikingly, a whole family of solutions (Figure 5B, shadowed sector)
crosses via the folded node singularity from the attracting to the repelling branch. This is possible since existence and uniqueness results of ordinary differential equations are violated along the fold
Figure 5: The reduced flow near a node singularity, B-C: The corresponding reduced flow near a folded node singularity (B) and the same reduced flow shown on the 2D folded critical manifold (C)
The folded singularity corresponding to a saddle-node bifurcation of a folded saddle and a folded node is called a folded saddle-node type I. Another possible bifurcation scenario is a transcritical bifurcation of a folded singularity and an ordinary singularity. In this case, the ordinary singularity crosses via a transcritical bifurcation from the attracting branch to the repelling branch or vice versa. The corresponding folded singularity at this transcritical bifurcation is called a folded saddle-node type II.
[edit] The Generalized Canard Phenomenon
[edit] Mixed Mode Oscillations (MMOs)
MMOs correspond to switching between small amplitude oscillations and
relaxation oscillations. These patterns were first discovered in the
famous Belousov-Zhabotinsky reaction and, since then, have been
frequently observed in experiments and models of chemical and
biological rhythms. One way to explain these patterns is based on
canards of folded node type. The reason is that canards of folded node
type can be responsible for small amplitude oscillations
(Wechselberger [2005]). A good intuition for MMOs is that a system
moves dynamically from a small amplitude oscillatory state to a
relaxation oscillatory state and the feature of the large relaxation
oscillation is to bring the system back to the basin of attraction of
the small amplitude oscillatory state. Other proposed mechanisms for
MMOs are break-up/loss of stability of a
Shilnikov homoclinic orbit (Koper [1995]), break-up of an invariant
torus (Larter and Steinmetz [1991]) or slow passage through a
delayed Hopf bifurcation (Larter et al. [1988]).
Figure 6: Explanation of MMOs via generalized canard phenomenon: the
upper center figure shows a time trace of a MMO pattern (); the
lower right figure summarizes Assumptions 1-3 in the phase space such that a singularly perturbed system possesses MMOs; the
lower left figure shows the 4 small amplitude oscillations in the phase space near the corresponding folded node singularity.
Figure 6 shows a MMO pattern consisting of 1 large amplitude oscillation and 4 small amplitude oscillations. In general, the symbol is assigned to a MMO pattern with large and small oscillations. The observed MMO pattern in Figure 6 can be explained as follows. Given a singularly perturbed system:
-
Assumption 1: The critical manifold of the singularly perturbed system is (locally) a folded surface.
-
Assumption 2: The corresponding reduced problem possesses a folded node singularity.
-
Assumption 3: There exist a singular periodic orbit (Figure 6, lower right) which consists of a segment on (blue) within the singular funnel (shadowed region) with the folded node singularity (black circle) as an endpoint, fast fibers (red) of the layer problem and a global return mechanism (green). The global return mechanism has to be specified.
Theorem (Brøns et al. [2006]): Given the above singularly
perturbed system under Assumptions 1-3. Then, for sufficiently small
there exist MMOs of type
It is possible to calculate the number
of small oscillations.
The number of small oscillations is given by
with where
and are the corresponding eigenvalues
of the folded node singularity whenever (except resonance cases).
[edit] References
-
E. Benoît [1983]; Systèmes lents-rapides dans et leurs
canards; Astérisque, 109-110(2), 159-191.
-
E. Benoît [2001]; Perturbation singulière en dimension trois : Canards
en un point pseudo-singulier nœud; Bulletin de la Société
Mathématique de France, 129-1, 91-113.
-
E. Benoît, J. F. Callot, F. Diener, M. Diener [1981]; Chasse au canard; Collectanea Mathematica, 31-32 (1-3), 37-119.
-
E. Benoît, A. Fruchard, R. Schaefke, G. Wallet [1998];
Solutions surstables des équations différentielles lentes-rapides
à point tournant; Annales de la Faculté des Sciences de Toulouse,
VII, 4, 627-658.
-
K. Bold, C. Edwards, J. Guckenheimer, S. Guharay, K. Hoffman, J. Hubbard, R. Oliva, W. Weckesser [2003];
The forced van der Pol oscillator II: Canards in the reduced flow; SIAM J. Appl. Dyn. Syst. 2, 570-608.
-
M. Brøns, K. Bar-Eli [1991]; Canard explosion and excitation in a model of the Belousov-Zhabotinskii reaction; J. Phys. Chem 95, 8706-8713.
-
M. Brøns, M. Krupa, M. Wechselberger [2006];
Mixed Mode Oscillations Due to the Generalized Canard Phenomenon;
Fields Institute Communications 49, 39-63.
-
J.L. Callot [1993]; Champs lents-rapides complexes à une dimension
lente; Annales scientifiques de l'Ecole Normale Supérieure, 4,26,149-173.
-
M. Canalis-Durand, J.P. Ramis, R. Schaefke, Y. Sibuya [2000]; Gevrey
solutions of singularly perturbed differential equations; J. Reine
Angew. Math., 518, 95-129.
- M. Diener [1984]; The canard unchained or how fast/slow dynamical problems bifurcate; The Mathematical Intelligencer, 6, 38-49.
-
J. Drover, J. Rubin, J. Su and B. Ermentrout [2004]; Analysis of a canard mechanism by which excitatory synaptic coupling can synchronize neurons at low firing frequencies; SIAM J. Appl. Math. 65, 69-92.
-
F. Dumortier, R. Roussarie [1996]; Canard cycles and center manifolds; Memoires of the AMS, 557.
-
W. Eckhaus [1983]; Relaxation oscillations including a standard chase on French ducks; in Asymptotic Analysis II, Springer Lecture Notes Math. 985, 449-494.
-
N. Fenichel [1979]; Geometric singular perturbation theory; JDE 31, 53-98.
-
J. Guckenheimer, R. Haiduc [2005]; Canards at folded nodes; Mosc. Math. J. 5, 91--103.
-
M. Koper [1995];
Bifurcations of mixed-mode oscillations in a three-variable autonomous van der Pol-Duffing model with a cross-shaped phase diagram; Physica D 80, 72-94.
-
M. Krupa, P. Szmolyan [2001]; Extending geometric singular perturbation theory to nonhyperbolic points - fold and canard points in two dimensions; SIAM J. of Math. Anal. 33, 286-314.
-
M. Krupa, P. Szmolyan [2001]; Relaxation oscillations and canard explosion; J. Diff. Eq. 174, 312--368.
-
R. Larter, C. Steinmetz [1991]; Chaos via mixed-mode oscillations; Phil. Trans. R. Soc. Lond. A 337, 291-298.
-
R. Larter, C. Steinmetz, B. Aguda [1988]; Fast-slow variable analysis of the transition to mixed-mode oscillations and chaos in the peroxidase reaction; J. Phys. Chem. 89, 6506-6514.
-
E.F. Mishchenko, Yu.S. Kolesov, A.Yu. Kolesov, N.Kh. Rhozov [1994]; Asymptotic methods in singularly perturbed systems, Monographs in Contemporary Mathematics; Consultants Bureau, New York, A Division of Plenum Publishing Cooperation 233 Spring Street, New York, N.Y. 10013.
-
J. Moehlis [2002]; Canards in a Surface Oxidation Reaction; J. Nonlinear Sci. 12, 319-345.
-
A.I. Neishtadt [1987]; Persistence of stability loss for dynamical bifurcations I, Diff. Eq. 23, 1385-1391.
-
B. Peng, V. Gaspar, K. Showalter [1991]; False bifurcations in chemical systems: Canards; Phil. Trans. R. Soc. Lond. A 337, 275-289.
-
H. G. Rotstein, N. Kopell, A. M. Zhabotinsky, I. R. Epstein [2003];
Canard phenomenon and localization of oscillations in the Belousov-Zhabotinsky reaction with global feedback; J. Chem. Phys. 119, 8824-8832.
-
P. Szmolyan, M. Wechselberger [2001]; Canards in JDE 177, 419-453.
-
G. Wallet [1986]; Entrée-sortie dans un tourbillon; Annales de l'Institut Fourier, 36, 4, 157-184.
-
M. Wechselberger [2005];
Existence and Bifurcation of Canards in in the case of a Folded Node;
SIAM J. Applied Dynamical Systems 4, 101-139.
Internal references
- Yuri A. Kuznetsov (2006) Andronov-Hopf bifurcation. Scholarpedia, 1(10):1858.
- Edward Ott (2006) Basin of attraction. Scholarpedia, 1(8):1701.
- Anatol M. Zhabotinsky (2007) Belousov-Zhabotinsky reaction. Scholarpedia, 2(9):1435.
- John Guckenheimer (2007) Bifurcation. Scholarpedia, 2(6):1517.
- Jack Carr (2006) Center manifold. Scholarpedia, 1(12):1826.
- James Meiss (2007) Dynamical systems. Scholarpedia, 2(2):1629.
- Eugene M. Izhikevich and Richard FitzHugh (2006) FitzHugh-Nagumo model. Scholarpedia, 1(9):1349.
- Jeff Moehlis, Kresimir Josic, Eric T. Shea-Brown (2006) Periodic orbit. Scholarpedia, 1(7):1358.
- Yuri A. Kuznetsov (2006) Saddle-node bifurcation. Scholarpedia, 1(10):1859.
- Leonid Pavlovich Shilnikov and Andrey Shilnikov (2007) Shilnikov bifurcation. Scholarpedia, 2(8):1891.
- Philip Holmes and Eric T. Shea-Brown (2006) Stability. Scholarpedia, 1(10):1838.
- Takashi Kanamaru (2007) Van der Pol oscillator. Scholarpedia, 2(1):2202.
[edit] External Links
[edit] See Also
Normal Hyperbolicity, Singularly perturbed systems, Relaxation Oscillator, FitzHugh-Nagumo Model, Bifurcations, Dynamical Systems