A COMPARATIVE STUDY OF THE HODGKIN–HUXLEY AND FITZHUGH–NAGUMO MODELS OF NEURON PULSE PROPAGATION

Abstract

The four-dimensional Hodgkin–Huxley equations are considered as the prototype for description of neural pulse propagation. Their mathematical complexity and sophistication prompted a simplified two-dimensional model, the FitzHugh–Nagumo equations, which display many of the former's dynamical features. Numerical and mathematical analysis are employed to demonstrate that the FitzHugh–Nagumo equations can provide quantitative predictions in close agreement with the Hodgkin–Huxley equations. The two most important parameters of a neural pulse are its speed c(T) and pulse height v max (T) and so numerical computations of these quantities predicted by the Hodgkin–Huxley equations are given over the entire temperature range T for stability of a neural pulse. Similarly, the FitzHugh–Nagumo equations are parameterized by two dimensionless quantities: a which determines the dynamics of the pulse front, and b whose departure from zero tailors the front to form the resultant pulse. Parallel computations are presented for the FitzHugh–Nagumo pulse whose relative simplicity permits analytic determination to close approximation of the dimensionless speed θ(a, b) and pulse height V max (a, b). It is shown that the two models are numerically identified by scaling according to c = 4904 θ cm/sec and v max = 115 V max mV where the numbers are a consequence of the experimental parameter values inherent to the Hodgkin–Huxley equations. With this connection, at a given temperature the Hodgkin–Huxley speed and pulse height determine unique values for the two FitzHugh–Nagumo parameters a and b. Approximate analytic solution for θ(a, b) allows construction of a three-dimensional [a, b, θ] state plot upon which a unique ridge defines, as a function of temperature, the speed and associated pulse height predicted by the Hodgkin–Huxley equations. The generality of the state plot suggests its application to other conductance models. Comparison of the Hodgkin–Huxley with the FitzHugh–Nagumo models highlight the quantitative limitations of the latter in the region of the minimum characterizing the back portion of the pulse. To overcome this limitation would require analytic extension of the FitzHugh–Nagumo dynamics to higher dimensionality.