"We describe a simple conductance-based model neuron that includes intra and
extracellular ion concentration dynamics and show that this model exhibits periodic
bursting. The bursting arises as the fast-spiking behavior of the neuron is modulated
by the slow oscillatory behavior in the ion concentration variables and vice versa. By
separating these time scales and studying the bifurcation structure of the neuron, we catalog
several qualitatively different bursting profiles that are strikingly similar to those seen in
experimental preparations. Our work suggests that ion concentration dynamics may play an
important role in modulating neuronal excitability in real biological systems."
These are the source codes used in
E. Barreto and J.R. Cressman, "Ion Concentration Dynamics as a
Mechanism for Neuronal Bursting",
Journal of Biological Physics 37, 361-373 (2011).
Link to the paper:
http://www.springerlink.com/content/v52215p195159211/
Author-generated version available at: http://arxiv.org/abs/1012.3124
The .ode files run with XPP.
The .c file is in C and uses headers and the rk4 routine from
Numerical Recipies.
The parameters that recreate the subplots of Figure 4 are shown using
XPP in a collection of jpgs included in this archive. Note that the
XPP code is written so as to measure time in milliseconds.
The authors recommend that modelers use the more efficient C code,
which outputs time in seconds, and which includes in the comments the
parameter values used for all five burst types.
Parameter values for all figures are also listed in the paper itself.