FitzHugh–Nagumo model

The FitzHugh–Nagumo model (FHN) describes a prototype of an excitable system (e.g., a neuron).

FitzHugh-Nagumo model in phase space, with . The green line is the cubic nullcline and the red line is the linear nullcline. The black lines are integral curves.
Graph of v with parameters I=0.5, a=0.7, b=0.8, and τ=12.5

It is an example of a relaxation oscillator because, if the external stimulus exceeds a certain threshold value, the system will exhibit a characteristic excursion in phase space, before the variables and relax back to their rest values.

This behaviour is a sketch for neural spike generations, with a short, nonlinear elevation of membrane voltage , diminished over time by a slower, linear recovery variable representing sodium channel reactivation and potassium channel deactivation, after stimulation by an external input current.[1]

The equations for this dynamical system read

The FitzHugh–Nagumo model is a simplified 2D version of the Hodgkin–Huxley model which models in a detailed manner activation and deactivation dynamics of a spiking neuron.

In turn, the Van der Pol oscillator is a special case of the FitzHugh–Nagumo model, with .

History

edit

It was named after Richard FitzHugh (1922–2007)[2] who suggested the system in 1961[3] and Jinichi Nagumo et al. who created the equivalent circuit the following year.[4]

In the original papers of FitzHugh, this model was called Bonhoeffer–Van der Pol oscillator (named after Karl-Friedrich Bonhoeffer and Balthasar van der Pol) because it contains the Van der Pol oscillator as a special case for  . The equivalent circuit was suggested by Jin-ichi Nagumo, Suguru Arimoto, and Shuji Yoshizawa.[5]

Qualitative analysis

edit

Qualitatively, the dynamics of this system is determined by the relation between the three branches of the cubic nullcline and the linear nullcline.

The cubic nullcline is defined by  .

The linear nullcline is defined by  .

In general, the two nullclines intersect at one or three points, each of which is an equilibrium point. At large values of  , far from origin, the flow is a clockwise circular flow, consequently the sum of the index for the entire vector field is +1. This means that when there is one equilibrium point, it must be a clockwise spiral point or a node. When there are three equilibrium points, they must be two clockwise spiral points and one saddle point.

  • If the linear nullcline pierces the cubic nullcline from downwards then it is a clockwise spiral point or a node.
  • If the linear nullcline pierces the cubic nullcline from upwards in the middle branch, then it is a saddle point.

The type and stability of the index +1 can be numerically computed by computing the trace and determinant of its Jacobian:  The point is stable iff the trace is negative. That is,  .

The point is a spiral point iff  . That is,  .

The limit cycle is born when a stable spiral point becomes unstable by Hopf bifurcation.[1]

Only when the linear nullcline pierces the cubic nullcline at three points, the system has a separatrix, being the two branches of the stable manifold of the saddle point in the middle.

  • If the separatrix is a curve, then trajectories to the left of the separatrix converge to the left sink, and similarly for the right.
  • If the separatrix is a cycle around the left intersection, then trajectories inside the separatrix converge to the left spiral point. Trajectories outside the separatrix converge to the right sink. The separatrix itself is the limit cycle of the lower branch of the stable manifold for the saddle point in the middle. Similarly for the case where the separatrix is a cycle around the right intersection.
  • Between the two cases, the system undergoes a homoclinic bifurcation.

Gallery figures: FitzHugh-Nagumo model, with  , and varying  . (They are animated. Open them to see the animation.)

See also

edit

References

edit
  1. ^ a b Sherwood, William Erik (2013), "FitzHugh–Nagumo Model", in Jaeger, Dieter; Jung, Ranu (eds.), Encyclopedia of Computational Neuroscience, New York, NY: Springer, pp. 1–11, doi:10.1007/978-1-4614-7320-6_147-1, ISBN 978-1-4614-7320-6, retrieved 2023-04-15
  2. ^ "Richard FitzHugh at the National Institute of Health - CHM Revolution". www.computerhistory.org. Archived from the original on 25 Mar 2023. Retrieved 2023-06-20.
  3. ^ FitzHugh, Richard (July 1961). "Impulses and Physiological States in Theoretical Models of Nerve Membrane". Biophysical Journal. 1 (6): 445–466. Bibcode:1961BpJ.....1..445F. doi:10.1016/S0006-3495(61)86902-6. PMC 1366333. PMID 19431309.
  4. ^ Nagumo, J.; Arimoto, S.; Yoshizawa, S. (October 1962). "An Active Pulse Transmission Line Simulating Nerve Axon". Proceedings of the IRE. 50 (10): 2061–2070. doi:10.1109/jrproc.1962.288235. ISSN 0096-8390. S2CID 51648050.
  5. ^ "SIAM: Obituaries: Jin-Ichi Nagumo". www.siam.org.

Further reading

edit
  • 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