Pseudo-Lyapunov exponents and predictability of Hodgkin-Huxley neuronal network dynamics

Yi Sun, Doug Zhou, Aaditya Rangan, David Cai

Research output: Contribution to journalArticle

Abstract

We present a numerical analysis of the dynamics of all-to-all coupled Hodgkin-Huxley (HH) neuronal networks with Poisson spike inputs. It is important to point out that, since the dynamical vector of the system contains discontinuous variables, we propose a so-called pseudo-Lyapunov exponent adapted from the classical definition using only continuous dynamical variables, and apply it in our numerical investigation. The numerical results of the largest Lyapunov exponent using this new definition are consistent with the dynamical regimes of the network. Three typical dynamical regimes - asynchronous, chaotic and synchronous, are found as the synaptic coupling strength increases from weak to strong. We use the pseudo-Lyapunov exponent and the power spectrum analysis of voltage traces to characterize the types of the network behavior. In the nonchaotic (asynchronous or synchronous) dynamical regimes, i.e., the weak or strong coupling limits, the pseudo-Lyapunov exponent is negative and there is a good numerical convergence of the solution in the trajectory-wise sense by using our numerical methods. Consequently, in these regimes the evolution of neuronal networks is reliable. For the chaotic dynamical regime with an intermediate strong coupling, the pseudo-Lyapunov exponent is positive, and there is no numerical convergence of the solution and only statistical quantifications of the numerical results are reliable. Finally, we present numerical evidence that the value of pseudo-Lyapunov exponent coincides with that of the standard Lyapunov exponent for systems we have been able to examine.

Original languageEnglish (US)
Pages (from-to)247-266
Number of pages20
JournalJournal of Computational Neuroscience
Volume28
Issue number2
DOIs
StatePublished - Apr 2010

Fingerprint

Spectrum Analysis

Keywords

  • Chaos
  • Hodgkin-Huxley neuron
  • Lyapunov exponents
  • Neuronal network
  • Numerical analysis

ASJC Scopus subject areas

  • Cellular and Molecular Neuroscience
  • Cognitive Neuroscience
  • Sensory Systems

Cite this

Pseudo-Lyapunov exponents and predictability of Hodgkin-Huxley neuronal network dynamics. / Sun, Yi; Zhou, Doug; Rangan, Aaditya; Cai, David.

In: Journal of Computational Neuroscience, Vol. 28, No. 2, 04.2010, p. 247-266.

Research output: Contribution to journalArticle

@article{5d3a639ad1874b478f072efb213a993e,
title = "Pseudo-Lyapunov exponents and predictability of Hodgkin-Huxley neuronal network dynamics",
abstract = "We present a numerical analysis of the dynamics of all-to-all coupled Hodgkin-Huxley (HH) neuronal networks with Poisson spike inputs. It is important to point out that, since the dynamical vector of the system contains discontinuous variables, we propose a so-called pseudo-Lyapunov exponent adapted from the classical definition using only continuous dynamical variables, and apply it in our numerical investigation. The numerical results of the largest Lyapunov exponent using this new definition are consistent with the dynamical regimes of the network. Three typical dynamical regimes - asynchronous, chaotic and synchronous, are found as the synaptic coupling strength increases from weak to strong. We use the pseudo-Lyapunov exponent and the power spectrum analysis of voltage traces to characterize the types of the network behavior. In the nonchaotic (asynchronous or synchronous) dynamical regimes, i.e., the weak or strong coupling limits, the pseudo-Lyapunov exponent is negative and there is a good numerical convergence of the solution in the trajectory-wise sense by using our numerical methods. Consequently, in these regimes the evolution of neuronal networks is reliable. For the chaotic dynamical regime with an intermediate strong coupling, the pseudo-Lyapunov exponent is positive, and there is no numerical convergence of the solution and only statistical quantifications of the numerical results are reliable. Finally, we present numerical evidence that the value of pseudo-Lyapunov exponent coincides with that of the standard Lyapunov exponent for systems we have been able to examine.",
keywords = "Chaos, Hodgkin-Huxley neuron, Lyapunov exponents, Neuronal network, Numerical analysis",
author = "Yi Sun and Doug Zhou and Aaditya Rangan and David Cai",
year = "2010",
month = "4",
doi = "10.1007/s10827-009-0202-2",
language = "English (US)",
volume = "28",
pages = "247--266",
journal = "Journal of Computational Neuroscience",
issn = "0929-5313",
publisher = "Springer Netherlands",
number = "2",

}

TY - JOUR

T1 - Pseudo-Lyapunov exponents and predictability of Hodgkin-Huxley neuronal network dynamics

AU - Sun, Yi

AU - Zhou, Doug

AU - Rangan, Aaditya

AU - Cai, David

PY - 2010/4

Y1 - 2010/4

N2 - We present a numerical analysis of the dynamics of all-to-all coupled Hodgkin-Huxley (HH) neuronal networks with Poisson spike inputs. It is important to point out that, since the dynamical vector of the system contains discontinuous variables, we propose a so-called pseudo-Lyapunov exponent adapted from the classical definition using only continuous dynamical variables, and apply it in our numerical investigation. The numerical results of the largest Lyapunov exponent using this new definition are consistent with the dynamical regimes of the network. Three typical dynamical regimes - asynchronous, chaotic and synchronous, are found as the synaptic coupling strength increases from weak to strong. We use the pseudo-Lyapunov exponent and the power spectrum analysis of voltage traces to characterize the types of the network behavior. In the nonchaotic (asynchronous or synchronous) dynamical regimes, i.e., the weak or strong coupling limits, the pseudo-Lyapunov exponent is negative and there is a good numerical convergence of the solution in the trajectory-wise sense by using our numerical methods. Consequently, in these regimes the evolution of neuronal networks is reliable. For the chaotic dynamical regime with an intermediate strong coupling, the pseudo-Lyapunov exponent is positive, and there is no numerical convergence of the solution and only statistical quantifications of the numerical results are reliable. Finally, we present numerical evidence that the value of pseudo-Lyapunov exponent coincides with that of the standard Lyapunov exponent for systems we have been able to examine.

AB - We present a numerical analysis of the dynamics of all-to-all coupled Hodgkin-Huxley (HH) neuronal networks with Poisson spike inputs. It is important to point out that, since the dynamical vector of the system contains discontinuous variables, we propose a so-called pseudo-Lyapunov exponent adapted from the classical definition using only continuous dynamical variables, and apply it in our numerical investigation. The numerical results of the largest Lyapunov exponent using this new definition are consistent with the dynamical regimes of the network. Three typical dynamical regimes - asynchronous, chaotic and synchronous, are found as the synaptic coupling strength increases from weak to strong. We use the pseudo-Lyapunov exponent and the power spectrum analysis of voltage traces to characterize the types of the network behavior. In the nonchaotic (asynchronous or synchronous) dynamical regimes, i.e., the weak or strong coupling limits, the pseudo-Lyapunov exponent is negative and there is a good numerical convergence of the solution in the trajectory-wise sense by using our numerical methods. Consequently, in these regimes the evolution of neuronal networks is reliable. For the chaotic dynamical regime with an intermediate strong coupling, the pseudo-Lyapunov exponent is positive, and there is no numerical convergence of the solution and only statistical quantifications of the numerical results are reliable. Finally, we present numerical evidence that the value of pseudo-Lyapunov exponent coincides with that of the standard Lyapunov exponent for systems we have been able to examine.

KW - Chaos

KW - Hodgkin-Huxley neuron

KW - Lyapunov exponents

KW - Neuronal network

KW - Numerical analysis

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

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

U2 - 10.1007/s10827-009-0202-2

DO - 10.1007/s10827-009-0202-2

M3 - Article

VL - 28

SP - 247

EP - 266

JO - Journal of Computational Neuroscience

JF - Journal of Computational Neuroscience

SN - 0929-5313

IS - 2

ER -