FitzHugh-Nagumo Model

From Scholarpedia

From Scholarpedia, the free peer-reviewed encyclopedia p.20223

Curator: Eugene M. Izhikevich, The Neurosciences Institute, San Diego, California
Curator: Richard FitzHugh, Maryland, U.S.A.

Figure 1: Circuit diagram of the tunnel-diode nerve model of Nagumo et al. (1962).

The FitzHugh-Nagumo model

\begin{matrix} \dot{V} & = & V-V^3/3 - W + I \\ \dot{W} & = & 0.08(V+0.7 - 0.8W) \end{matrix}

is a two-dimensional simplification of the Hodgkin-Huxley model of spike generation in squid giant axons. Here,

  • V is the membrane potential,
  • W is a recovery variable,
  • I is the magnitude of stimulus current.

This system was suggested by FitzHugh (1961), who called it "Bonhoeffer-van der Pol model", and the equivalent circuit by Nagumo et al. (1962; see Fig.9).

Contents

[edit]

Principal Assumptions

The motivation for the FitzHugh-Nagumo model was to isolate conceptually the essentially mathematical properties of excitation and propagation from the electrochemical properties of sodium and potassium ion flow. The model consists of

The model is sometimes written in the abstract form

\begin{matrix} \dot{V} & = & f(V) - W + I \\ \dot{W} & = & a(bV -cW) \end{matrix}

where f(V) is a polynomial of third degree, and a, b, and c are constant parameters.

[edit]

Explained Phenomena

Enlarge
Figure 2: Phase portrait and physiological state diagram of FitzHugh-Nagumo model (modified from FitzHugh 1961).

While the Hodgkin-Huxley Model is more realistic and biophysically sound, only projections of its four-dimensional phase trajectories can be observed. The simplicity of the FitzHugh-Nagumo model permits the entire solution to be viewed at once. This allows a geometrical explanation of important biological phenomena related to neuronal excitability and spike-generating mechanism.

The phase portrait of the FitzHugh-Nagumo model in the figure depicts

  • V-nullcline, which is the N-shaped curve obtained from the condition \dot{V}=0, where the sign of \dot{V} passes through zero,
  • W-nullcline, which is a straight line obtained from the condition \dot{W}=0, where the sign of \dot{W} passes through zero, and
  • some typical trajectories starting with various initial conditions.

The intersection of nullclines is an equilibrium (because \dot{V}=\dot{W}=0), which may be unstable if it is on the middle branch of the V-nullcline, i.e., when I is strong enough. In this case, the model exhibits periodic (tonic spiking) activity.

Most of the labels in the figure are explained in the text. The "no man's land" region of the phase space is highly unstable, containing trajectories starting very close to the quasi-threshold.

Figure 3: Absence of all-or-none spikes in the FitzHugh-Nagumo model.
[edit]

Absence of All-or-None Spikes

The FitzHugh-Nagumo model explained the absence of all-or-none spikes in the HH model in response to stimuli, i.e., pulses of the injected current I. Weak stimuli (small pulses of I) result in small-amplitude trajectories that correspond to subthreshold responses; stronger stimuli result in intermediate-amplitude trajectories that correspond to partial-amplitude spikes; and strong stimuli result in large-amplitude trajectories that correspond to suprathreshold response -- firing a spike.

[edit]

Absence of Threshold

Similarly to the HH model, FitzHugh-Nagumo model does not have a well-defined firing threshold. This feature is the consequence of the absence of all-or-none responses, and it is related, from the mathematical point of view, to the absence of a saddle equilibrium (FitzHugh 1955). The apparent illusion of threshold dynamics and all-or-none responses in both models is due to the existence of the ‘’quasi-threshold’’, which is a canard trajectory that follows the unstable (middle) branch of the N-shaped V-nullcline. Nearby trajectories diverge sharply away from the canard trajectory to the left or right, producing an apparently "all-or-none" response and threshold-like behavior. (A point moving along a canard trajectory is like a tightrope walker walking slowly along a rope; if he loses his balance, he quickly falls away from the rope to one side or the other.)

Figure 4: Excitation block in the FitzHugh-Nagumo model.
[edit]

Excitation Block

The FitzHugh-Nagumo model explains the excitation block phenomenon, i.e., the cessation of repetitive spiking as the amplitude of the stimulus current increases. When I is weak or zero, the equilibrium (intersection of nullclines) is on the left (stable) branch of V-nullcline, and the model is resting. Increasing I shifts the nullcline upward and the equilibrium slides onto the middle (unstable) branch of the nullcline. The model exhibits periodic spiking activity in this case. Increasing the stimulus further shifts the equilibrium to the right (stable) branch of the N-shaped nullcline, and the oscillations are blocked (by excitation!). The precise mathematical mechanism involves appearance and disappearance of a limit cycle attractor, and it is reviewed in detail by Izhikevich (2007).

Figure 5: Anodal break excitation (post-inhibitory rebound spike) in the FitzHugh-Nagumo model.
[edit]

Anodal Break Excitation

The FitzHugh-Nagumo model explained the phenomenon of post-inhibitory (rebound) spikes, called anodal break excitation at that time. As the stimulus I becomes negative (hyperpolarization), the resting state shifts to the left. As the system is released from hyperpolarization (anodal break), the trajectory starts from a point far below the resting state (outside the quasi-threshold, see the first figure), makes a large-amplitude excursion, i.e., fires a transient spike, and then returns to the resting state.

Figure 6: Spike accommodation to slowly increasing stimulus in the FitzHugh-Nagumo model.
[edit]

Spike Accommodation

The FitzHugh-Nagumo model explained the dynamical mechanism of spike accommodation in HH-type models. When stimulation strength I increases slowly, the neuron remains quiescent. The resting equilibrium of the FitzHugh-Nagumo model shifts slowly to the right, and the state of the system follows it smoothly without firing spikes. In contrast, when the stimulation is increased abruptly, even by a smaller amount, the trajectory could not go directly to the new resting state, but fires a transient spike; see figure. Geometrically, this phenomenon is similar to the post-inhibitory (rebound) response.

Figure 7: Traveling pulse the FitzHugh-Nagumo reaction-diffusion model. All spatial points are projected onto their V and W coordinates, so that the traveling pulse looks like a circle on the phase plane (notice that because of the diffusion term, the points do not exhibit relaxation oscillations).
[edit]

Traveling Waves

The FitzHugh-Nagumo equations became a favorite model for reaction-diffusion systems

\begin{matrix} \dot{V} & = & f(V) - W + I + V_{xx} \\ \dot{W} & = & a(bV -cW) \qquad \qquad \end{matrix}

which simulate propagation of waves in excitable media, such as heart tissue or nerve fiber. Here, the diffusion term V_{xx} is the second derivative with respect to spatial variable x. Its success is mostly due to the fact that the model is analytically tractable, and hence it allows derivation of many important properties of traveling pulses without resort to computer simulations.

Enlarge
Figure 8: Richard FitzHugh with analog computer at NIH, ca. 1960.
[edit]

History

Shortly after the publication of Hodgkin and Huxley's equations for the squid giant axon, Richard FitzHugh was working at the Biophysics Laboratory of the National Institutes of Health (NIH) in Bethesda, Maryland. He undertook an analysis of the mathematical properties of their equations. He used the new techniques of nonlinear mechanics which had been developed by Russian mathematicians led by A. A. Andronov. This was before digital computers became easily accessible. John Moore and FitzHugh started by planning how to programm an analog computer which could be used to solve the Hodgkin-Huxley equations. The equipment needed included operational amplifiers, function generators, multipliers, and an ink pen plotter. The laboratory purchased the computer, which occupied four floor-to-ceiling relay racks, full of vacuum tubes. These were continually failing, and FitzHugh had to find and replace several tubes a week, requiring some detective work. The heat from all these tubes sometimes overloaded the air conditioning, so that on hot summer days he had to take off his shirt and wear shorts to be comfortable.

With this computer he plotted solutions of the HH equations. The operation of the analog computer required the skill of an electronic engineer as well as those of a mathematician.

In this analog computer, the variables in the HH equations, V, m, h, n, are represented by voltages. Each variable was transformed into a voltage with a separate scale factor. These voltage signals were passed from one unit to another.

One of the basic units of the computer was the operational amplifier (Op Amp). Each one occupied a metal box about six inches long with six vacuum tubes on top. (Today one can buy a tiny solid-state chip with several Op Amps in it.) To these Op Amps one could connect highly accurate resistors and capacitors to perform the mathematical operations of addition, subtraction and also integration of the signals representing first derivatives with respect to time. Another type of unit was the function generator. There were six of these, for the alpha and beta functions in the HH equations. The voltage value of each function was set into the unit at the evenly spaced points, at intervals of ten volts, between which the function was approximated by straight line segments. This was not accurate enough for the computation. To provide smoothing of the function curve, a high frequency zigzag signal of ten volts peak-to-peak amplitude from a signal generator was added to the input. This was of too high a frequency to appear in the slowly changing voltage passed to the plotter. The effect was to average the function locally, to produce a smooth curve, which however no longer exactly fitted accurately at the set points. Finally, after readjusting the output to accurately fit the function curve at the original set points, an accurate fit to all the function was obtained.

All the units were connected by the maze of wires shown in the photo of FitzHugh operating the computer. The wires protrude from an insulated board, underneath which they made contact with terminals connected to the various units. To change the connections on the board, for a different problem, it was disengaged from its position, exposing the terminals behind. The computer power supply delivered voltages from -100 volts to +100 volts, and if one of the terminals was accidentally touched, one might receive a nasty shock.

The first thing to do in the morning was to turn on the computer and let it warm up until the voltages from the power supply stabilized. Then computation could begin.

In order to distinguish between the physical basis of the HH equations, in terms of the flow through the axon membrane of sodium and potassium ions, on the one hand, and the phenomena of excitation above a threshold value of stimulus, and propagation along the axon, on the other, it seemed that it would be useful to simplify their equations in order to isolate these properties from each other. At the suggestion of his lab chief, Dr. Kenneth S. (Kacy) Cole, FitzHugh modified the Van der Pol equations for the nonlinear relaxation oscillator. The result had a stable resting state, from which it could be excited by a sufficiently large electrical stimulus to produce an impulse. A large enough constant current stimulus produced a train of impulses (FitzHugh 1961,1969).

Enlarge
Figure 9: The original Nagumo circuit is held as a trust in the laboratory of Dr. K. Aihara in University of Tokyo. Photo provided by H. Suetani and K. Aihara with permissions.

These equations were similar to those describing the electronic circuit called a monostable multivibrator. At about the same time, an electronic circuit was built by the Japanese engineer Jin-Ichi Nagumo, using tunnel (Esaki) diodes; see Fig.9. These diodes have a current-voltage curve similar to the cubic shape used in FitzHugh's equations. These equations have since become known as the FitzHugh-Nagumo equations, though they were originally called "Bonhoeffer-van der Pol model" by FitzHugh. Reprogramming the analog computer for the FitzHugh-Nagumo equations was much simpler. Only two multipliers and no function generators were needed.

Aside from mathematicians in the computer Division and in John Rinzel's Mathematical Biology Section, there were few people at NIH who were interested in the application of mathematics to biological problems. Most did experiments only. FitzHugh fulfilled his obligations in such work by going each summer to the Marine Biological Laboratory at Woods Hole Massachusetts, where he assisted Cole and Moore in their experiments, using Cole's voltage clamp technique on the squid which were caught at sea there. Otherwise, at NIH there was not much interest in mathematical models.

Within 30-40 years of this first forage by Fitzhugh and co-workers into mathematical neuroscience whole communities of mathematicians, physicists, and engineers studying nonlinear dynamics of biological systems have started and matured. Over that time the Fitzhugh-Nagumo model has remained the prototypical example of an excitable system and new discoveries with this system are still being made.

[edit]

References

  • FitzHugh R. (1955) Mathematical models of threshold phenomena in the nerve membrane. Bull. Math. Biophysics, 17:257--278
  • FitzHugh R. (1961) Impulses and physiological states in theoretical models of nerve membrane. Biophysical J. 1:445-466
  • FitzHugh R. (1969) Mathematical models of excitation and propagation in nerve. Chapter 1 (pp. 1-85 in H.P. Schwan, ed. Biological Engineering, McGraw-Hill Book Co., N.Y.)
  • Nagumo J., Arimoto S., and Yoshizawa S. (1962) An active pulse transmission line simulating nerve axon. Proc IRE. 50:2061–2070.
[edit]

Further reading

  • Rocsoreanu C., Georgescu A., and Giurgiteanu N. (2000) The FitzHugh-Nagumo Model: Bifurcation and Dynamics. Kluwer Academic Publishers, Boston, 2000.
[edit]

See also

Hodgkin-Huxley model, Relaxation Oscillator, Canards, Periodic Orbit, Reaction-Diffusion Systems, Traveling Wave, Van der Pol Oscillator

[edit]

External links


Eugene M. Izhikevich, Richard FitzHugh (2006) FitzHugh-Nagumo Model. Scholarpedia, p.3193
Created: 1 March 2006, reviewed: 18 August 2006, accepted: 23 September 2006
Editor: Frances K. Skinner, Toronto Western Research Institute, University Health Network, and University of Toronto
Retrieved from "http://www.scholarpedia.org/article/FitzHugh-Nagumo_Model"