Library-based numerical reduction of the Hodgkin - Huxley neuron for network simulation

Yi Sun, Doug Zhou, Aaditya Rangan, David Cai

Research output: Contribution to journalArticle

Abstract

We present an efficient library-based numerical method for simulating the Hodgkin - Huxley (HH) neuronal networks. The key components in our numerical method involve (i) a pre-computed high resolution data library which contains typical neuronal trajectories (i.e., the time-courses of membrane potential and gating variables) during the interval of an action potential (spike), thus allowing us to avoid resolving the spikes in detail and to use large numerical time steps for evolving the HH neuron equations; (ii) an algorithm of spike-spike corrections within the groups of strongly coupled neurons to account for spike-spike interactions in a single large time step. By using the library method, we can evolve the HH networks using time steps one order of magnitude larger than the typical time steps used for resolving the trajectories without the library, while achieving comparable resolution in statistical quantifications of the network activity, such as average firing rate, interspike interval distribution, power spectra of voltage traces. Moreover, our large time steps using the library method can break the stability requirement of standard methods (such as Runge-Kutta (RK) methods) for the original dynamics. We compare our library-based method with RK methods, and find that our method can capture very well phase-locked, synchronous, and chaotic dynamics of HH neuronal networks. It is important to point out that, in essence, our library-based HH neuron solver can be viewed as a numerical reduction of the HH neuron to an integrate-and-fire (I&F) neuronal representation that does not sacrifice the gating dynamics (as normally done in the analytical reduction to an I&F neuron).

Original languageEnglish (US)
Pages (from-to)369-390
Number of pages22
JournalJournal of Computational Neuroscience
Volume27
Issue number3
DOIs
StatePublished - 2009

Fingerprint

Libraries
Neurons
Membrane Potentials
Action Potentials

Keywords

  • Chaos
  • Hodgkin - Huxley neuron
  • Integrate-and-fire neurons
  • Library method
  • Neuronal network
  • Numerical algorithm

ASJC Scopus subject areas

  • Cellular and Molecular Neuroscience
  • Cognitive Neuroscience
  • Sensory Systems

Cite this

Library-based numerical reduction of the Hodgkin - Huxley neuron for network simulation. / Sun, Yi; Zhou, Doug; Rangan, Aaditya; Cai, David.

In: Journal of Computational Neuroscience, Vol. 27, No. 3, 2009, p. 369-390.

Research output: Contribution to journalArticle

@article{4a29c881a2064261aa77ea7c86e075af,
title = "Library-based numerical reduction of the Hodgkin - Huxley neuron for network simulation",
abstract = "We present an efficient library-based numerical method for simulating the Hodgkin - Huxley (HH) neuronal networks. The key components in our numerical method involve (i) a pre-computed high resolution data library which contains typical neuronal trajectories (i.e., the time-courses of membrane potential and gating variables) during the interval of an action potential (spike), thus allowing us to avoid resolving the spikes in detail and to use large numerical time steps for evolving the HH neuron equations; (ii) an algorithm of spike-spike corrections within the groups of strongly coupled neurons to account for spike-spike interactions in a single large time step. By using the library method, we can evolve the HH networks using time steps one order of magnitude larger than the typical time steps used for resolving the trajectories without the library, while achieving comparable resolution in statistical quantifications of the network activity, such as average firing rate, interspike interval distribution, power spectra of voltage traces. Moreover, our large time steps using the library method can break the stability requirement of standard methods (such as Runge-Kutta (RK) methods) for the original dynamics. We compare our library-based method with RK methods, and find that our method can capture very well phase-locked, synchronous, and chaotic dynamics of HH neuronal networks. It is important to point out that, in essence, our library-based HH neuron solver can be viewed as a numerical reduction of the HH neuron to an integrate-and-fire (I&F) neuronal representation that does not sacrifice the gating dynamics (as normally done in the analytical reduction to an I&F neuron).",
keywords = "Chaos, Hodgkin - Huxley neuron, Integrate-and-fire neurons, Library method, Neuronal network, Numerical algorithm",
author = "Yi Sun and Doug Zhou and Aaditya Rangan and David Cai",
year = "2009",
doi = "10.1007/s10827-009-0151-9",
language = "English (US)",
volume = "27",
pages = "369--390",
journal = "Journal of Computational Neuroscience",
issn = "0929-5313",
publisher = "Springer Netherlands",
number = "3",

}

TY - JOUR

T1 - Library-based numerical reduction of the Hodgkin - Huxley neuron for network simulation

AU - Sun, Yi

AU - Zhou, Doug

AU - Rangan, Aaditya

AU - Cai, David

PY - 2009

Y1 - 2009

N2 - We present an efficient library-based numerical method for simulating the Hodgkin - Huxley (HH) neuronal networks. The key components in our numerical method involve (i) a pre-computed high resolution data library which contains typical neuronal trajectories (i.e., the time-courses of membrane potential and gating variables) during the interval of an action potential (spike), thus allowing us to avoid resolving the spikes in detail and to use large numerical time steps for evolving the HH neuron equations; (ii) an algorithm of spike-spike corrections within the groups of strongly coupled neurons to account for spike-spike interactions in a single large time step. By using the library method, we can evolve the HH networks using time steps one order of magnitude larger than the typical time steps used for resolving the trajectories without the library, while achieving comparable resolution in statistical quantifications of the network activity, such as average firing rate, interspike interval distribution, power spectra of voltage traces. Moreover, our large time steps using the library method can break the stability requirement of standard methods (such as Runge-Kutta (RK) methods) for the original dynamics. We compare our library-based method with RK methods, and find that our method can capture very well phase-locked, synchronous, and chaotic dynamics of HH neuronal networks. It is important to point out that, in essence, our library-based HH neuron solver can be viewed as a numerical reduction of the HH neuron to an integrate-and-fire (I&F) neuronal representation that does not sacrifice the gating dynamics (as normally done in the analytical reduction to an I&F neuron).

AB - We present an efficient library-based numerical method for simulating the Hodgkin - Huxley (HH) neuronal networks. The key components in our numerical method involve (i) a pre-computed high resolution data library which contains typical neuronal trajectories (i.e., the time-courses of membrane potential and gating variables) during the interval of an action potential (spike), thus allowing us to avoid resolving the spikes in detail and to use large numerical time steps for evolving the HH neuron equations; (ii) an algorithm of spike-spike corrections within the groups of strongly coupled neurons to account for spike-spike interactions in a single large time step. By using the library method, we can evolve the HH networks using time steps one order of magnitude larger than the typical time steps used for resolving the trajectories without the library, while achieving comparable resolution in statistical quantifications of the network activity, such as average firing rate, interspike interval distribution, power spectra of voltage traces. Moreover, our large time steps using the library method can break the stability requirement of standard methods (such as Runge-Kutta (RK) methods) for the original dynamics. We compare our library-based method with RK methods, and find that our method can capture very well phase-locked, synchronous, and chaotic dynamics of HH neuronal networks. It is important to point out that, in essence, our library-based HH neuron solver can be viewed as a numerical reduction of the HH neuron to an integrate-and-fire (I&F) neuronal representation that does not sacrifice the gating dynamics (as normally done in the analytical reduction to an I&F neuron).

KW - Chaos

KW - Hodgkin - Huxley neuron

KW - Integrate-and-fire neurons

KW - Library method

KW - Neuronal network

KW - Numerical algorithm

UR - http://www.scopus.com/inward/record.url?scp=70350568843&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=70350568843&partnerID=8YFLogxK

U2 - 10.1007/s10827-009-0151-9

DO - 10.1007/s10827-009-0151-9

M3 - Article

VL - 27

SP - 369

EP - 390

JO - Journal of Computational Neuroscience

JF - Journal of Computational Neuroscience

SN - 0929-5313

IS - 3

ER -