- Poster presentation
- Open Access
Extensive chaotic dynamics in neural networks in the balanced state
© Kreissl et al; licensee BioMed Central Ltd. 2009
- Published: 13 July 2009
- Lyapunov Exponent
- Chaotic Dynamic
- Entropy Production
- Balance State
- Neuron Model
Neurons embedded in operational cortical networks fire action potentials in highly irregular sequences . The balanced state of cortical networks provides an attractive explanation of this temporally irregular behaviour . In this state, neurons are driven by large input fluctuations, resulting from a dynamically achieved balance of excitation and inhibition. It has been shown that the balanced state can robustly emerge from the collective dynamics of spiking neuron networks . The detailed dynamics of the balanced state, however, is not well understood. Initially, Vreeswijk and Sompolinsky  found a kind of "hyper"-chaotic dynamics in networks of binary neurons, characterized by an infinite positive Lyapunov exponent, which is hard to reconcile with classical notions of nonlinear dynamics. More recently, Zillmer et al.  and Jahnke et al.  reported stable dynamics in networks of leaky integrate and fire neurons. Notably, all of these previous works neglect a major cellular mechanism of dynamic instability, the instability of the membrane voltage near the spike threshold.
To clarify which kind of balanced state dynamics is expected in networks of neuron models describing such dynamic spike generation, we analyzed networks of pulse-coupled theta neurons in sparse, random networks in the balanced state. The theta neuron model is a standard model of type I neuronal excitability . Based on the analytic solution of the single neuron dynamics, we derived a map that was used for event based simulations and obtained the exact single spike Jacobian matrix for the entire network. This Jacobian was used to calculate the whole set of Lyapunov exponents, following the standard procedure . From the Lyapunov exponents we acquired the entropy production via Pesin's formula and the attractor dimension via the Kaplan-Yorke conjecture.
In general, the studied networks in the balanced state exhibited chaotic dynamics. The largest Lyapunov exponent was positive and finite, giving evidence to conventional chaos. The entropy production and dimension of the chaotic attractor both scaled linearly with the number of neurons in the networks, hence chaos in the balanced state appears to be extensive spatio-temporal. Furthermore, we found a high attractor dimension, corresponding to many effective degrees of freedom that may potentially encode information in the network. Yet, this is accompanied by a rapid loss of information of 0.5 bits per spike. We conclude that extensive chaos should be expected in general for networks of type I neurons in the balanced state. If information loss in cortical networks is as high as in our model, information could hardly be encoded in detailed spiking patterns beyond the immediate stimulus response.
- Softky WR, Koch C: The highly irregular firing of cortical cells is inconsistent with temporal integration of random EPSPs. J Neurosci. 1993, 13: 334-350.PubMedGoogle Scholar
- Van Vreeswijk C, Sompolinsky H: Chaotic balanced state in a model of cortical circuits. Neural Comput. 1998, 10: 1321-1371. 10.1162/089976698300017214.PubMedView ArticleGoogle Scholar
- Zillmer R, Livi R, Politi A, Torcini A: Desynchronization in diluted neural networks. Phys Rev E. 2006, 74: 036203-10.1103/PhysRevE.74.036203.View ArticleGoogle Scholar
- Jahnke S, Memmesheimer RM, Timme M: Stable irregular dynamics in complex neural networks. Phys Rev Lett. 2008, 100: 048102-10.1103/PhysRevLett.100.048102.PubMedView ArticleGoogle Scholar
- Gutkin BS, Ermentrout GB: Dynamics of membrane excitability determine interspike interval variability: a link between spike generation mechanisms and cortical spike train statistics. Neural Comput. 1998, 10: 1047-1065. 10.1162/089976698300017331.PubMedView ArticleGoogle Scholar
- Benettin G, Galgani L, Giorgilli A, Strelcyn JM: Lyapunov characteristic exponents for smooth dynamical systems and for hamiltonian systems; a method for computing all of them, Part I & Part II. Meccanica. 1980, 15: 9-30. 10.1007/BF02128236.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd.