- Poster presentation
- Open Access
- Published:

# Multistability in large scale models of brain activity

*BMC Neuroscience*
**volume 14**, Article number: P84 (2013)

Rich spatiotemporal dynamical patterns, observed in the brain at rest, reveal several large-scale functional networks, presumably involved in different brain functions. In parallel, structural networks obtained by Diffusion Spectrum Imaging ("connectome") identify several interconnected sub-networks that overlap with the functional networks [3]. Neural mass model simulations aim at developing realistic models of the brain activity. In particular, multiple fixed-point attractors can be identified and spontaneous alternation between several brain "states" can be obtained [2]. The number and the organization of the "fixed-point" attractors strongly vary depending of the parameters used in the model.

In order to provide a general view of connectome-based dynamical systems, we developed a simplified neural mass model inspired by the continuous Hopfield network [4], to which we add a stochastic component, and a dynamic threshold (eq.2) that prevents reaching the two trivial attractors of system (eq.1).

where *W* is the N = 66 nodes human connectome [3], *x* is the node potential, *φ* is the node output, *θ* is the adaptive threshold, *η*(*t*) is a white noise and σ the diffusion parameter. Three main parameters (excitation strength ||*W*||, inhibition strength *Ω* and node excitability γ) are systematically varied in order to identify the fixed points attractors of the dynamics (Figure 1).

The adaptive threshold plays the role of a global inhibition and allows stabilizing many attractors, where the Ω / ||W|| ratio controls the sparsity of the activity. Then multistable dynamical systems are obtained in a large region of the parameter space, allowing to identify many different attractor patterns. The number of attractors is found to increase with the value of gamma (node excitability), at the expense of their intrinsic stability. Itinerant dynamics is obtained when a significant noise level is introduced in the system. Then, several "attractive" patterns can be reached on a single trajectory, where the duration of the visit reflects the stability of the pattern.

The simple dynamical system we have implemented allows exploring a large region of the parameters space, but is also capable of reproducing some aspects of the brain's spontaneous activity (switching between different attractors). The set of attractors we find in some regions of the parameters space share similarities with the different functional modes observed in the resting-state activity [1].

## References

- 1.
Damoiseaux JS, Rombouts S, Barkhof F, Scheltens P, Stam C, Smith P, Beckmann C: Consistent resting-state networks across healthy subjects. PNAS. 2006

- 2.
Deco G, Jirsa VK: Ongoing cortical activity at rest: criticality, multistability and ghost attractors. J Neuroscience. 2012

- 3.
Hagmann P, Cammoun L, Gigandet X, Meuli R, Honey C, Wedeen V, Sporns O: Mapping the structural core of human cerebral cortex. PLoS Biol. 2008

- 4.
Hopfield JJ: Neural networks and physical systems with emergent collective computational abilities. Proc Natl Acad Sci USA. 1982

## Author information

### Affiliations

### Corresponding author

## Rights and permissions

**Open Access**
This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License (
https://creativecommons.org/licenses/by/2.0
), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## About this article

### Cite this article

Golos, M., Jirsa, V.K. & Daucé, E. Multistability in large scale models of brain activity.
*BMC Neurosci* **14, **P84 (2013). https://0-doi-org.brum.beds.ac.uk/10.1186/1471-2202-14-S1-P84

Published:

DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/1471-2202-14-S1-P84

### Keywords

- Adaptive Threshold
- Spontaneous Alternation
- Hopfield Network
- Node Potential
- Neural Mass Model