Two-Dimensional Sodium plus PotassiumNeuron Model

Requires a Wolfram Notebook System

Interact on desktop, mobile and cloud with the free Wolfram CDF Player or other Wolfram Language products.

Requires a Wolfram Notebook System

Edit on desktop, mobile and cloud with any Wolfram Language product.

This Demonstration simulates the two-dimensional persistent sodium plus potassium () neuron model described in [1].






where is the capacitance density (for the squid axon, , so is fixed and set to one); is the voltage across a cell membrane in mV; is the total current density flowing across a patch of the cell membrane in ; , , are the equilibrium potentials (reverse potentials) in mV; , , are conductance densities in ; and are the activation variables (probabilities of activation gates to be open) and when with the subscript , they are voltage-sensitive steady-state activation functions; and is the time constant. Note that stands for leak (ohmic leak current density is given by and it is mostly carried by ions); and are sodium and potassium ions, respectively. Finally, the Boltzmann functions and each have two parameters: and , where satisfies and similarly for sodium. The , are slope factors; that is, they adjust the slope of the Boltzmann functions at (smaller values of result in a steeper Boltzmann function).

We plot four graphs. The main graph is the phase plane of the potassium conductance versus membrane voltage with the vector field (black), nullclines (orange and green), and voltage trajectory (red) superimposed. We also plot the membrane voltage versus time, sodium () and potassium () conductances versus time, and membrane currents versus time. Blue corresponds to sodium and green to potassium; nullcline is green as well.

There are four bookmarks, one for each bifurcation of the resting state. The parameter values for each bifurcation are listed below.


Contributed by: Joshua Ryan and Jan Fiala (August 2015)
(UVa Wise)
Open content licensed under CC BY-NC-SA




[1] E. M. Izhikevich, Dynamical Systems in Neuroscience: The Geometry of Excitability and Bursting, Cambridge, MA: MIT Press, 2007.

Feedback (field required)
Email (field required) Name
Occupation Organization
Note: Your message & contact information may be shared with the author of any specific Demonstration for which you give feedback.