Bifurcation analysis of a periodically forced relaxation oscillator: differential model versus phase-resetting map.

We compare the dynamics of the periodically forced FitzHugh-Nagumo oscillator in its relaxation regime to that of a one-dimensional discrete map of the circle derived from the phase-resetting response of this oscillator (the "phase-resetting map"). The forcing is a periodic train of Gaussian-shaped pulses, with the width of the pulses much shorter than the intrinsic period of the oscillator. Using numerical continuation techniques, we compute bifurcation diagrams for the periodic solutions of the full differential equations, with the stimulation period being the bifurcation parameter. The period-1 solutions, which belong either to isolated loops or to an everywhere-unstable branch in the bifurcation diagram at sufficiently small stimulation amplitudes, merge together to form a single branch at larger stimulation amplitudes. As a consequence of the fast-slow nature of the oscillator, this merging occurs at virtually the same stimulation amplitude for all the period-1 loops. Again using continuation, we show that this stimulation amplitude corresponds, in the circle map, to a change of topological degree from one to zero. We explain the origin of this coincidence, and also discuss the translational symmetry properties of the bifurcation diagram.

Continuation and bifurcation analysis of a periodically forced excitable system.

The response of an excitable cell to periodic electrical stimulation is modeled using the FitzHugh-Nagumo (FHN) system submitted to a gaussian-shaped pacing, the width of which is small compared with the action potential duration. The influence of the amplitude and the period of the stimulation is studied using numerical continuation and bifurcation techniques (AUTO97 software). Results are discussed in the light of prior experimental and theoretical findings. In particular, agreement with the documented behavior of periodically stimulated cardiac cells and squid axons is discussed. As previously reported, we find many different "M:N" periodic solutions, period-doubling sequences leading to seemingly chaotic regimes, and bistability phenomena. In addition, the use of continuation techniques has allowed us to track unstable solutions of the system and thus to determine how the different stable rhythms are connected with each other in a bifurcation diagram. Depending on the stimulus amplitude, the aspect of the bifurcation diagram with the stimulus period as main varying parameter can vary from very simple to very complex. In its most developed structure, this bifurcation diagram consists of a main "tree" of period-2(P) branches, where the 1:1, 1:0, 2:2, 2:1,... rhythms are located, and of several closed loops made up of period-{N x 2(P)} branches (N>2), isolated from each other and from the main tree. It is mainly on such loops that N:1 rhythms (N>2) on one hand, and N:N-1 or Wenckebach rhythms (N>2) on the other hand, are located. Stable M:N and M:N-1 rhythms (M>or=N) can be found on the same branch of solutions. They are separated by a region of unstable solutions at small stimulus amplitudes, but this region shrinks gradually as the stimulus amplitude is raised, until it finally disappears. We believe that this property is related to the excitability characteristics of the FHN system. It would be interesting to know if it has any correspondence in the behavior of real excitable cells.

Standing waves in the FitzHugh-Nagumo model of cardiac electrical activity.

When excitable media are submitted to appropriate time dependent boundary conditions, a standing wavelike pattern can be observed in the system, as shown in recent experiments. In the present analysis, the physical mechanism explaining the occurrence of such space-time patterns is shown to be a competition between Ohmic diffusion and an action potential propagation across the system, coupled with the existence of refractory states for excitable media.

Rolling and slipping motion of Euler's disk.

We present an experimental study of the motion of a circular disk spun onto a table. With the help of a high speed video system, the temporal evolution of (i) the inclination angle alpha, (ii) the angular velocity omega, and (iii) the precession rate Omega are studied. The influence of the mass of the disk as well as the friction between the disk and the supporting surface are considered. Both inclination angle and angular velocity are observed to decrease according to a power law. We also show that the precession rate diverges as the motion stops. Measurements are performed very near the collapse as well as on long range times. Times to collapse have been also measured. Results are compared with previous theoretical and experimental works. The major source of energy dissipation is found to be the slipping of the disk on the plane.