A Global Bifurcation and the Appearance of a One-Dimensional Spiral Wave in Excitable Media

It is well known that excitable media can sustain "fast" stable traveling pulses of excitation. In models for which analysis is tractable, the fast pulse solution is known to appear through a saddle- node bifurcation accompanied by a "slow" unstable traveling pulse. Furthermore, the uniform rest state is also a stable solution. It is generally assumed that the boundary between the basins of attractions of the rest state and fast pulse (i.e., threshold) consists of the stable manifold of the slow pulse. We use numerical experiments to explore this issue. Our results indicate that, near the saddle-node bifurcation, the stable manifold of the slow pulse does indeed act as the threshold between the rest state and fast pulse. However, further away from the saddle-node bifurcation, a global bifurcation involving the heteroclinic connections between slow and fast pulses occurs. This bifurcation gives rise to an unstable periodic solution that has been referred to as a one-dimensional spiral wave. Beyond the bifurcation point, the stable manifold of the one-dimensional spiral wave appears to act as the threshold between the rest state and fast pulse.

[1]  Christopher Jones,et al.  Stability of the travelling wave solution of the FitzHugh-Nagumo system , 1984 .

[2]  R. FitzHugh Impulses and Physiological States in Theoretical Models of Nerve Membrane. , 1961, Biophysical journal.

[3]  J. Keener,et al.  Singular perturbation theory of traveling waves in excitable media (a review) , 1988 .

[4]  Nancy Kopell,et al.  Target Patterns and Horseshoes from a Perturbed Central-Force Problem: Some Temporally Periodic Solutions to Reaction-Diffusion Equations , 1981 .

[5]  C Antzelevitch,et al.  Characteristics of Reflection as a Mechanism of Reentrant Arrhythmias and Its Relationship to Parasystole , 1980, Circulation.

[6]  C. Luo,et al.  A dynamic model of the cardiac ventricular action potential. I. Simulations of ionic currents and concentration changes. , 1994, Circulation research.

[7]  A. Hodgkin,et al.  A quantitative description of membrane current and its application to conduction and excitation in nerve , 1952, The Journal of physiology.

[8]  James P. Keener,et al.  Mathematical physiology , 1998 .

[9]  G. W. Beeler,et al.  Reconstruction of the action potential of ventricular myocardial fibres , 1977, The Journal of physiology.

[10]  Björn Sandstede,et al.  Defects in Oscillatory Media: Toward a Classification , 2004, SIAM J. Appl. Dyn. Syst..

[11]  S. Hastings ON THE EXISTENCE OF HOMOCLINIC AND PERIODIC ORBITS FOR THE FITZHUGH-NAGUMO EQUATIONS , 1976 .

[12]  G. Carpenter A geometric approach to singular perturbation problems with applications to nerve impulse equations , 1977 .

[13]  A. Winfree Electrical instability in cardiac muscle: phase singularities and rotors. , 1989, Journal of theoretical biology.

[14]  Aoxiang Xu,et al.  Two forms of spiral-wave reentry in an ionic model of ischemic ventricular myocardium. , 1998, Chaos.

[15]  R C Barr,et al.  Reflection after delayed excitation in a computer model of a single fiber. , 1992, Circulation research.

[16]  G. Flores Stability analysis for the slow traveling pulse of the Fitzhugh-Nagumo system , 1991 .

[17]  C F Starmer,et al.  Vulnerability in an excitable medium: analytical and numerical studies of initiating unidirectional propagation. , 1993, Biophysical journal.

[18]  Excitation , 1964 .

[19]  Trine Krogh-Madsen,et al.  Resetting and termination of reentry in a loop-and-tail cardiac model. , 2008, Physical review. E, Statistical, nonlinear, and soft matter physics.

[20]  William H. Calvin,et al.  Impulses reflected from dorsal root ganglia and from focal nerve injuries , 1976, Brain Research.

[21]  S. Yoshizawa,et al.  An Active Pulse Transmission Line Simulating Nerve Axon , 1962, Proceedings of the IRE.

[22]  Bard Ermentrout,et al.  Reflected Waves in an Inhomogeneous Excitable Medium , 1996, SIAM J. Appl. Math..

[23]  Natalia A Trayanova,et al.  Virtual Electrode‐Induced Positive and Negative Graded Responses: , 2003, Journal of cardiovascular electrophysiology.

[24]  James P. Keener,et al.  Stability conditions for the traveling pulse: Modifying the restitution hypothesis. , 2002, Chaos.

[25]  W Rall,et al.  Changes of action potential shape and velocity for changing core conductor geometry. , 1974, Biophysical journal.

[26]  C. Morris,et al.  Voltage oscillations in the barnacle giant muscle fiber. , 1981, Biophysical journal.

[27]  G. Flores The stable manifold of the standing wave of the Nagumo equation , 1989 .

[28]  J Rinzel,et al.  Traveling wave solutions of a nerve conduction equation. , 1973, Biophysical journal.

[29]  H. McKean Nagumo's equation , 1970 .

[30]  J Jalife,et al.  Excitation, Conduction, and Reflection Impulses in Isolated Bovine and Canine Cardiac Purkinje Fibers , 1981, Circulation research.

[31]  C. Luo,et al.  A dynamic model of the cardiac ventricular action potential. II. Afterdepolarizations, triggered activity, and potentiation. , 1994, Circulation research.

[32]  J W Moore,et al.  Propagation of action potentials in inhomogeneous axon regions. , 1975, Federation proceedings.

[33]  E. A. Ermakova,et al.  Rotating spiral waves in a modified Fitz-Hugh-Nagumo model , 1984 .

[34]  R. G. Casten,et al.  Perturbation analysis of an approximation to the Hodgkin-Huxley theory , 1975 .

[35]  J. McLeod,et al.  The approach of solutions of nonlinear diffusion equations to travelling front solutions , 1977 .

[36]  J. Bell,et al.  Study of propagation along nonuniform excitable fibers. , 1994, Mathematical biosciences.

[37]  J. Rinzel Excitation dynamics: insights from simplified membrane models. , 1985, Federation proceedings.

[38]  S. Baccus,et al.  Synaptic facilitation by reflected action potentials: enhancement of transmission when nerve impulses reverse direction at axon branch points. , 1998, Proceedings of the National Academy of Sciences of the United States of America.

[39]  Calculation of the threshold surface for nerve equations , 1990 .

[40]  G. Bard Ermentrout,et al.  One-dimensional λ-ω target patterns: Empirical stability tests , 1980 .

[41]  John J. Tyson,et al.  When Time Breaks Down: The Three‐Dimensional Dynamics of Electrochemical Waves and Cardiac Arrhythmias , 1988 .

[42]  J. Keener,et al.  The effect of spatial scale of resistive inhomogeneity on defibrillation of cardiac tissue. , 2003, Journal of theoretical biology.

[43]  J. Keener,et al.  The biphasic mystery: why a biphasic shock is more effective than a monophasic shock for defibrillation. , 1999, Journal of theoretical biology.

[44]  James P Keener,et al.  The topology of defibrillation. , 2004, Journal of theoretical biology.

[45]  Wanda Krassowska,et al.  Initiation of propagation in a one-dimensional excitable medium , 1997 .

[46]  Sitabhra Sinha,et al.  Termination of reentry in an inhomogeneous ring of model cardiac cells. , 2002, Physical review. E, Statistical, nonlinear, and soft matter physics.

[47]  Glass,et al.  Resetting and Annihilation of Reentrant Abnormally Rapid Heartbeat. , 1995, Physical review letters.

[48]  J. Jalife,et al.  Reflected reentry in nonhomogeneous ventricular muscle as a mechanism of cardiac arrhythmias. , 1984, Circulation.