Synchronised firing patterns in a random network of adaptive exponential integrate-and-fire

Synchronised firing patterns in a random network of adaptive exponential integrate-and-fire


We have studied neuronal synchronisation in a random network of adaptive exponential integrate-and-fire neurons. We study how spiking or bursting synchronous behaviour appears as a function of the coupling strength and the probability of connections, by constructing parameter spaces that identify these synchronous behaviours from measurements of the inter-spike interval and the calculation of the order parameter. Moreover, we verify the robustness of synchronisaton by applying an external perturbation to each neuron. The simulations show that bursting synchronisation is more robust than spike synchronisation.

synchronisation, integrate-and-fire, network

[cor]Corresponding author:

1 Introduction

The concept of synchronistion is based on the adjustment of rhythms of oscillating systems due to their interaction (1). Synchronisation phenomenon was recognised by Huygens in the 17th century, time when he performed experiments to understand this phenomenon (2). To date, several kinds of synchronisation among coupled systems were reported, such as complete (3), phase (4); (5), lag (6), and collective almost synchronisation (7).

Neuronal synchronous rhythms have been observed in a wide range of researches about cognitive functions (8); (9). Electroencephalography and magnetoencephalography studies have been suggested that neuronal synchronization in the gamma frequency plays a functional role for memories in humans (10); (11). Steinmetz et al. (12) investigated the synchronous behaviour of pairs of neurons in the secondary somatosensory cortex of monkey. They found that attention modulates oscillatory neuronal synchronisation in the somatosensory cortex. Moreover, in the literature it has been proposed that there is a relationship between conscious perception and synchronisation of neuronal activity (13).

We study spiking and bursting synchronisation between neuron in a neuronal network model. A spike refers to the action potential generated by a neuron that rapidly rises and falls (14), while bursting refers to a sequence of spikes that are followed by a quiescent time (15). It was demonstrated that spiking synchronisation is relevant to olfactory bulb (16) and is involved in motor cortical functions (17). The characteristics and mechanisms of bursting synchronisation were studied in cultured cortical neurons by means of planar electrode array (18). Jefferys Haas discovered synchronised bursting of CA1 hippocampal pyramidal cells (19).

There is a wide range of mathematical models used to describe neuronal activity, such as the cellular automaton (20), the Rulkov map (21), and differential equations (22); (23). One of the simplest mathematical models and that is widely used to depict neuronal behaviour is the integrate-and-fire (24), which is governed by a linear differential equation. A more realistic version of it is the adaptive exponential integrate-and-fire (aEIF) model which we consider in this work as the local neuronal activity of neurons in the network. The aEIF is a two-dimensional integrate-and-fire model introduced by Brette Gerstner (25). This model has an exponential spike mechanism with an adaptation current. Touboul Brette (26) studied the bifurcation diagram of the aEIF. They showed the existence of the Andronov-Hopf bifurcation and saddle-node bifurcations. The aEIF model can generate multiple firing patterns depending on the parameter and which fit experimental data from cortical neurons under current stimulation (27).

In this work, we focus on the synchronisation phenomenon in a randomly connected network. This kind of network, also called Erdös-Rényi network (28), has nodes where each pair is connected according to a probability. The random neuronal network was utilised to study oscillations in cortico-thalamic circuits (29) and dynamics of network with synaptic depression (30). We built a random neuronal network with unidirectional connections that represent chemical synapses.

We show that there are clearly separated ranges of parameters that lead to spiking or bursting synchronisation. In addition, we analyse the robustness to external perturbation of the synchronisation. We verify that bursting synchronisation is more robustness than spiking synchronisation. However, bursting synchronisation requires larger chemical synaptic strengths, and larger voltage potential relaxation reset to appear than those required for spiking synchronisation.

This paper is organised as follows: in Section II we present the adaptive exponential integrate-and-fire model. In Section III, we introduce the neuronal network with random features. In Section IV, we analyse the behaviour of spiking and bursting synchronisation. In the last Section, we draw our conclusions.

2 Adaptive exponential integrate-and-fire

As a local dynamics of the neuronal network, we consider the adaptive exponential integrate-and-fire (aEIF) model that consists of a system of two differential equations (25) given by


where is the membrane potential when a current is injected, is the membrane capacitance, is the leak conductance, is the resting potential, is the slope factor, is the threshold potential, is an adaptation variable, is the time constant, and is the level of subthreshold adaptation. If reaches the threshold , a reset condition is applied: and . In our simulations, we consider pF, nS, mV, mV, mV, pA, ms, nS, and mV (27).

The firing pattern depends on the reset parameters and . Table 1 exhibits some values that generate five different firing patterns (Fig. 1). In Fig. 1 we represent each firing pattern with a different colour in the parameter space : adaptation in red, tonic spiking in blue, initial bursting in green, regular bursting in yellow, and irregular in black. In Figs. 1a, 1b, and 1c we observe adaptation, tonic spiking, and initial burst pattern, respectively, due to a step current stimulation. Adaptation pattern has increasing inter-spike interval during a sustained stimulus, tonic spiking pattern is the simplest regular discharge of the action potential, and the initial bursting pattern starts with a group of spikes presenting a frequency larger than the steady state frequency. The membrane potential evolution with regular bursting is showed in Fig. 1d, while Fig. 1e displays irregular pattern.

Firing patterns Fig. b (pA) (mV) Layout
adaptation 1(a) 60.0 -68.0 red
tonic spiking 1(b) 5.0 -65.0 blue
initial burst 1(c) 35.0 -48.8 green
regular bursting 1(d) 40.0 -45.0 yellow
irregular 1(e) 41.2 -47.4 black
Table 1: Reset parameters.
Figure 1: (Colour online) Parameter space for the firing patterns as a function of the reset parameters and . (a) Adaptation in red, (b) tonic spiking in blue, (c) initial bursting in green, (d) regular bursting in yellow, and (e) irregular in black.

As we have interest in spiking and bursting synchronisation, we separate the parameter space into a region with spike and another with bursting patterns (Fig. 2). To identify these two regions of interest, we use the coefficient of variation (CV) of the neuronal inter-spike interval (ISI), that is given by


where is the standard deviation of the ISI normalised by the mean (31). Spiking patterns produce . Parameter regions that represent the neurons firing with spiking pattern are denoted by gray colour in Fig. 2. Whereas, the black region represents the bursting patterns, which results in .

Figure 2: Parameter space for the firing patterns as a function of the reset parameters and . Spike pattern in region I () and bursting pattern in region II () are separated by white circles.

3 Spiking or bursting synchronisation

In this work, we constructed a network where the neurons are randomly connected (28). Our network is given by


where is the membrane potential of the neuron , is the synaptic conductance, is the synaptic reversal potential, is the synaptic time constant, is the synaptic weight, is the adjacency matrix, is the external perturbation, and is randomly distributed in the interval .

The schematic representation of the neuronal network that we have considered is illustrated in Fig 3. Each neuron is randomly linked to other neurons with a probability by means of directed connections. When is equal to 1, the neuronal network becames an all-to-all network. A network with this topology was used by Borges et al. (32) to study the effects of the spike timing-dependent plasticity on the synchronisation in a Hodgkin-Huxley neuronal network.

Figure 3: Schematic representation of the neuronal network where the neurons are connected according to a probability .

A useful diagnostic tool to determine synchronous behaviour is the complex phase order parameter defined as (33)


where and are the amplitude and angle of a centroid phase vector, respectively, and the phase is given by


where corresponds to the time when a spike () of a neuron happens (). We have considered the beginning of the spike when mV. The value of the order parameter magnitude goes to 1 in a totally synchronised state. To study the neuronal synchronisation of the network, we have calculated the time-average order-parameter, that is given by


where is the time window for calculating .

Figs. 4a, 4b, and 4c show the raster plots for nS, nS, and nS, respectively, considering mV, , and pA, where the dots correspond to the spiking activities generated by neurons. For nS (Fig. 4a) the network displays a desynchonised state, and as a result, the order parameter values are very small (black line in Fig. 4d). Increasing the synaptic conductance for nS, the neuronal network exhibits spike synchronisation (Fig. 4b) and the order parameter values are near unity (red line in Fig. 4d). When the network presents bursting synchronisation (Fig. 4c), the order parameter values vary between and (blue line in Fig. 4d). to the time when the neuron are firing.

Figure 4: (Colour online) Raster plot for (a) nS, (b) nS, and (c) nS, considering mV, , and pA. In (d) the order parameter is computed for nS (black line), nS (red line), and nS (blue line).

In Fig. 5a we show as a function of for , pA (black line), pA (red line), and pA (blue line). The three results exhibit strong synchronous behaviour () for many values of when nS . However, for nS, it is possible to see synchronous behaviour only for pA in the range . In addition, we calculate the coefficient of variation (CV) to determine the range in where the neurons of the network have spiking or bursting behaviour (Fig. 5b). We consider that for CV (black dashed line) the neurons exhibit spiking behaviour, while for CV the neurons present bursting behaviour. We observe that in the range for pA there is spiking sychronisation, and bursting synchronisation for nS.

Figure 5: (Colour online) (a) Time-average order parameter and (b) CV for mV, , pA (black line), pA (red line), and pA (blue line).

4 Parameter space of synchronisation

The synchronous behaviour depends on the synaptic conductance and the probability of connections. Fig. 6 exhibits the time-averaged order parameter in colour scale as a function of and . We verify a large parameter region where spiking and bursting synchronisation is strong, characterised by . The regions I and II correspond to spiking and bursting patterns, respectively, and these regions are separated by a white line with circles. We obtain the regions by means of the coefficient of variation (CV). There is a transition between region I and region II, where neurons initially synchronous in the spike, loose spiking synchronicity to give place to a neuronal network with a regime of bursting synchronisation.

Figure 6: (Colour online) for mV and pA, where the colour bar represents the time-average order parameter. The regions I (spike patterns) and II (bursting patterns) are separated by the white line with circles.

We investigate the dependence of spiking and bursting synchronisation on the control parameters and . To do that, we use the time average order parameter and the coefficient of variation. Figure 7 shows that the spike patterns region (region I) decreases when increases. This way, the region I for pA and mV of parameters leading to no synchronous behaviour (Fig. 7a), becomes a region of parameters that promote synchronised bursting (Fig. 7b and 7c). However, a large region of desynchronised bursting appears for nS about mV and pA in the region II (Fig. 7b). For nS, we see, in Fig. 7c, three regions of desynchronous behaviour, one in the region I for pA, other in region II for pA, and another one is located around the border (white line with circles) between regions I and II for pA.

Figure 7: (Colour online) Parameter space for , (a) nS, (b) nS, and (c) nS, where the colour bar represents the time-average order parameter. The regions I (spike patterns) and II (bursting patterns) are separated by white circles.

It has been found that external perturbations on neuronal networks not only can induce synchronous behaviour (34); (35), but also can suppress synchronisation (36). Aiming to study the robustness to perturbations of the synchronous behaviour, we consider an external perturbation (3). It is applied on each neuron with an average time interval of about ms and with a constant intensity during ms.

Figure 8 shows the plots for , where the regions I and II correspond to spiking and bursting patterns, respectively, separated by white line with circles, and the colour bar indicates the time-average order parameter values. In this Figure, we consider mV, pA, (a) pA, (b) pA, and (c) pA. For pA (Fig. 8a) the perturbation does not suppress spike synchronisation, whereas for pA the synchronisation is completely suppressed in region I (Fig. 8b). In Fig. 8c, we see that increasing further the constant intensity for pA, the external perturbation suppresses also bursting synchronisation in region II. Therefore,the synchronous behavior in region II is more robustness to perturbations than in the region I, due to the fact that the region II is in a range with high and values, namely strong coupling and high connectivity.

Figure 8: (Colour online) for mV, pA, (a) pA, (b) pA, and (c) pA.

In order to understand the perturbation effect on the spike and bursting patterns, we consider the same values of and as Fig. 7a. Figure 9 exhibits the space parameter , where is equal to pA. The external perturbation suppresses synchronisation in the region I, whereas we observe synchronisation in region II. The synchronous behaviour in region II can be suppressed if the constant intensity is increased. Therefore, bursting synchronisation is more robustness to perturbations than spike synchronisation.

Figure 9: (Colour online) for nS, , and pA, where the colour bar represents the time-average order parameter. The regions I (spike patterns) and II (bursting patterns) are separated by white line with circles.

5 Conclusion

In this paper, we studied the spiking and bursting synchronous behaviour in a random neuronal network where the local dynamics of the neurons is given by the adaptive exponential integrate-and-fire (aEIF) model. The aEIF model can exhibit different firing patterns, such as adaptation, tonic spiking, initial burst, regular bursting, and irregular bursting.

In our network, the neurons are randomly connected according to a probability. The larger the probability of connection, and the strength of the synaptic connection, the more likely is to find bursting synchronisation.

It is possible to suppress synchronous behaviour by means of an external perturbation. However, synchronous behaviour with higher values of and , which typically promotes bursting synchronisation, are more robust to perturbations, then spike synchronous behaviour appearing for smaller values of these parameters. We concluded that bursting synchronisation provides a good environment to transmit information when neurons are strongly perturbed (large ).


This study was possible by partial financial support from the following Brazilian government agencies: CNPq, CAPES, and FAPESP (2011/19296-1 and 2015/07311-7). We also wish thank Newton Fund and COFAP.


  1. journal: Arxiv


  1. Pikovsky, A., Rosenblum, M., Kurths, J. (2001). Synchronization: A universal concept in nonlinear sciences. Cambridge University Press.
  2. Bennet, M., Schatz, M. F., Rockwood, H., Wiesenfeld, K. (2002). Huygens’s clocks. Proceedings: Mathematical, Physical and Engineering Sciences, 458, 563-579.
  3. Li, Y. Li, C. (2016). Complete synchronization of delayed chaotic neural networks by intermittent control with two switches in a control period. Neurocomputing, 173, 1341-1347.
  4. Pereira, T., Baptista, M. S., Kurths, J. (2007). Onset of phase synchronisation in neurons with chemical synapse. International Journal of Bifurcation and Chaos, 17, 3545-3549.
  5. Batista, C. A. S., Lopes, S. R., Viana, R. L., Batista, A. M. (2010). Delayed feedback control of bursting synchronization in a scale-free neuronal network. Neural Networks, 23, 114-124.
  6. Huang, J., Li, C., Huang, T., He, X. (2014). Finite-time lag synchronization of delayed neural networks. Neurocomputing, 139, 145-149.
  7. Baptista, M. S., Ren, H.-P., Swarts, J. C. M., Carareto, R., Nijmeijer, H., Grebogi, C. (2012). Collective Almost Synchronisation in Complex Networks. Plos One, 7, e48118.
  8. Wang, X.-J. (2010). Neurophysiological and Computational Principles of Cortical Rhythms in Cognition. Physiological Reviews, 90, 1195-1268.
  9. Hutcheon, B., Yarom, Y. (2000). Resonance oscillation and the intrinsic frequency preferences of neurons. Trends In NeuroScience, 23, 216-222.
  10. Axmacher, N., Mormann, F., Fernandez, G., Elger, C. E., Fell, J. (2006). Memory formation by neuronal synchronization. Brain Research Reviews, 52, 170-182.
  11. Fell, J. Axmacher, N. (2011). The role of phase synchronization in memory processes. Nature Reviews Neuroscience, 12, 105-118.
  12. Steinmetz, P. N., Roy, A., Fitzgerald, P. J., Hsiao, S. S., Johnson, K. O., Niebur, E. (2000). Attention modulates synchronized neuronal firing in primate somatosensory cortex. Nature, 404, 187-190.
  13. Hipp, J. F., Engel, A. K., Siegel, M. (2011). Oscillatory synchronization in large-scale cortical networks predicts perception. Neuron, 69, 387-396.
  14. de Lange, E. Hasler, M. (2008). Predicting single spikes and spike patterns with the Hindmarsh-Rose model. Biological Cybernetics, 99, 349-360.
  15. Wu, Y., Lu, W., Lin, W., Leng, G., Feng, J. (2012). Bifurcations of emergent bursting in a neuronal network. Plos One, 7, e38402.
  16. Davison, A. P., Feng, J., Brown, D. (2001). Spike synchronization in a biophysically-detailed model of the olfactory bulb. Neurocomputing, 38, 515-521.
  17. Riehle, A., Grun, S., Diesmann, M., Aertsen, A. (1997). Spike synchronization and rate modulation differentially involved in motor cortical function. Science, 278, 1950-1953.
  18. Maeda, E., Robinson P. C., Kawana A. (1995). The mechanisms of generation and propagation of synchronized bursting in developing networks of cortical neurons. The Journal of Neuroscience, 15, 6834-6845.
  19. Jefferys, J. G. R. Haas, H. L. (1982). Synchronized bursting of CA1 hippocampal pyramidal cells in the absence of synaptic transmission. Nature, 300, 448-450.
  20. Viana, R. L., Borges, F. S., Iarosz, K. C., Batista, A. M., Lopes, S. R., Caldas, I. L. (2014). Dynamic range in a neuron network with electrical and chemical synapses. Communications in Nonlinear Science and Numerical Simulation, 19, 164-172.
  21. Rulkov, N. F. (2001). Regularization of synchronized chaotic bursts. Physical Review Letters, 86, 183-186.
  22. Hodgkin, A. L., Huxley, A. F. (1952). A quantitative description of membrane current and its application to conduction and excitation in nerve. The Journal of Physiology, 117, 500-544.
  23. Hindmarsh, L. J., Rose, R. M. (1984). A model of neuronal bursting using three coupled first order differential equations. Proceeding of the Royal Society of London Series B, 221, 87-102.
  24. Lapicque, L. (1907). Recherches quantitatives sur l’excitation electrique des nerfs traitée comme une polarization. Journal de Physiologie et de Pathologie Générale, 9, 620-635.
  25. Brette, R. Gerstner, W. (2005). Adaptive exponential integrate-and-fire model as an effective description of neuronal activity. Journal of Neurophysiology, 94, 3637-3642.
  26. Touboul, J. Brette, R. (2008). Dynamics and bifurcations of the adaptive exponential integrate-and-fire model. Biological Cybernetics, 99, 319-334.
  27. Naud, R., Marcille, N., Clopath, C., Gerstner, G. (2008). Firing patterns in the adaptive exponential integrate-and-fire model. Biological Cybernetics, 99, 335-347.
  28. Erdös, P. Rényi, A. (1959). On random graphs, I. Publicationes Mathematicae, 6, 290-297.
  29. Gelenbe, E., Cramer, C. (1998). Oscillatory corticothalamic response to somatosensory input. Biosystems, 48, 67-75.
  30. Senn, W., Wyler, K., Streit, J., Larkum, M., Lüscher H.-R., Mey, H., Müller, L., Stainhauser, D., Vogt, K., Wannier, Th. (1996). Dynamics of a random neural network with synaptic depression. Neural Networks, 9, 575-588.
  31. Gabbiani, F. Koch, C. (1998). Principles of spike train analysis, in Methods in Neuronal Modeling: From Ions to Networks, 2nd Edn, eds C. Koch and I. Segev (Cambridge, MA: MIT Press), 313-360.
  32. Borges, R. R., Borges, F. S., Lameu, E. L., Batista, A. M., Iarosz, K. C., Caldas, I. L., Viana, R. L., Sanjuán (2016). Effects of the spike timing-dependent plasticity on the synchronisation in a random Hodgkin-Huxley neuronal network. Communications in Nonlinear Science Numerical Simulation, 34, 12-22.
  33. Kuramoto, Y. (2003). Chemical oscillations, waves and turbulence. Dover, New York.
  34. Baptista, M. S., Zhou, C., Kurths, J. (2006). Information transmission in phase synchronous chaotic arrays. Chinese Physics Letters, 23, 560-563.
  35. Zhang, W., Li, C., Huang, T., Xiao, M. (2015). Synchronization of neural networks with stochastic perturbation via aperiodically intermittent control. Neural Networks, 71, 105-111.
  36. Lameu, E. L., Borges, F. S., Borges, R. R., Iarosz, K. C., Caldas, I. L., Batista, A. M., Viana, R. L., Kurths, J. (2016). Suppression of phase synchronisation in network based on cat’s brains. Chaos, 26, 043107.
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minumum 40 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description