Jasmine Otto

  [CV]  

Ph.D. Computational Media, UC Santa Cruz

[back]Fitz-Hugh Nagumo in one spatial dimension (with p5.js and vega-lite)

Click (tap) to depolarize the left end. (Advanced: Modify parameter dictionary 'par' from the JS console.)

Voltage may be greater (or less than) resting potential. Light borders mark inactivated segments.


$$ \begin{align*} \dot{v}=&v\left(v-a\right)\left(1-v\right)-w\\ \dot{w}=&\epsilon\left(v-bw\right) \end{align*} $$

The action potential (a depolarization of the membrane) travels down an unmyelinated axon by diffusion. (Since myelinated axon experiences fast saltatory conduction, we might pretend each segment of our model is a node of Ranvier, since the ion channels of each segment recover independently.)

The action potential dies out unless its shape is self-reinforcing (stereotyped), so the axon robustly transmits a signal. (The frequency of spikes, not their amplitude, encodes qualitative information.)

The generalized FitzHugh-Nagumo equations consider only 'excitation' and 'inactivation' at a point along the axon. Although the literal ion populations, ligand-gated channels, and osmotic gradients are all neglected, this simple nonlinear dynamics produces the desired qualia: a stereotyped response (regular spiking at fixed frequency, seen as an excursion in the phase-plane) under fixed excitatory input.


10min animated explanations of action potential and resting potential exist on Osmosis.


The action potential is generated when the voltage across the membrane rises from resting potential (about -70mV due to active K+ transport out of the axon) past a threshold (which triggers Na+ channels to open, leading to fast equilibriation).

As more ion channels slowly open to let K+ exit, others become inactivated to prevent more Na+ from entering. Thus the downstroke, where the neuron cannot respond to further stimuli (the absolute refractory period), ends in an overshoot (the relative refractory period).

Continuous stimulation causes the neuron to spike with a certain frequency. In this simple model, that frequency is fixed.


The change in excitation is a cubic function of excitation. It is negative for both low and high excitation, but positive for excitation between the threshold and peak potential. It turns out that excitation always trends towards zero (absent external stimulus), but if it is pushed past threshold potential, it takes a long detour towards peak potential.

The change in inactivation is a linear function of both exictation and inactivation. It is positive for high excitation, but quickly becomes negative as inactivation rises. So inactivation tends to peak slightly after excitation does.

Notice that no specific count is kept of the two most important species of charge carrier, K+ and Na+. In this model, diffusion between segments is treated as a linear function of the difference in excitement between them, which does not seem to be a disastrous oversimplification.

Theme derived from Minimal (orderedlist) for Jekyll with flat icons (alexisdoreau).