Chaotic phase synchronization in scale-free networks of bursting neurons

Synchronization in scale-free neural networks under electromagnetic radiation

The functional networks of the human brain exhibit the structural characteristics of a scale-free topology, and these neural networks are exposed to the electromagnetic environment. In this paper, we consider the effects of magnetic induction on synchronous activity in biological neural networks, and the magnetic effect is evaluated by the four-stable discrete memristor. Based on Rulkov neurons, a scale-free neural network model is established. Using the initial value and the strength of magnetic induction as control variables, numerical simulations are carried out. The research reveals that the scale-free neural network exhibits multiple coexisting behaviors, including resting state, period-1 bursting synchronization, asynchrony, and chimera states, which are dependent on the different initial values of the multi-stable discrete memristor. In addition, we observe that the strength of magnetic induction can either enhance or weaken the synchronization in the scale-free neural network when the parameters of Rulkov neurons in the network vary. This investigation is of significant importance in understanding the adaptability of organisms to their environment.

Analyzing bursting synchronization in structural connectivity matrix of a human brain under external pulsed currents

Cognitive tasks in the human brain are performed by various cortical areas located in the cerebral cortex. The cerebral cortex is separated into different areas in the right and left hemispheres. We consider one human cerebral cortex according to a network composed of coupled subnetworks with small-world properties. We study the burst synchronization and desynchronization in a human neuronal network under external periodic and random pulsed currents. With and without external perturbations, the emergence of bursting synchronization is observed. Synchronization can contribute to the processing of information, however, there are evidences that it can be related to some neurological disorders. Our results show that synchronous behavior can be suppressed by means of external pulsed currents.

Role of assortativity in predicting burst synchronization using echo state network

In this study, we use a reservoir computing based echo state network (ESN) to predict the collective burst synchronization of neurons. Specifically, we investigate the ability of ESN in predicting the burst synchronization of an ensemble of Rulkov neurons placed on a scale-free network. We have shown that a limited number of nodal dynamics used as input in the machine can capture the real trend of burst synchronization in this network. Further, we investigate the proper selection of nodal inputs of degree-degree (positive and negative) correlated networks. We show that for a disassortative network, selection of different input nodes based on degree has no significant role in the machine's prediction. However, in the case of assortative network, training the machine with the information (i.e., time series) of low degree nodes gives better results in predicting the burst synchronization. The results are found to be consistent with the investigation carried out with a continuous time Hindmarsh-Rose neuron model. Furthermore, the role of hyperparameters like spectral radius and leaking parameter of ESN on the prediction process has been examined. Finally, we explain the underlying mechanism responsible for observing these differences in the prediction in a degree correlated network.

In-vitro engineered human cerebral tissues mimic pathological circuit disturbances in 3D

In-vitro modeling of brain network disorders such as epilepsy remains a major challenge. A critical step is to develop an experimental approach that enables recapitulation of in-vivo-like three-dimensional functional complexity while allowing local modulation of the neuronal networks. Here, by promoting matrix-supported active cell reaggregation, we engineered multiregional cerebral tissues with intact 3D neuronal networks and functional interconnectivity characteristic of brain networks. Furthermore, using a multi-chambered tissue-culture chip, we show that our separated but interconnected cerebral tissues can mimic neuropathological signatures such as the propagation of epileptiform discharges.

A method is developed to engineer cerebral tissues with intact 3D neuronal networks, mimicking neuropathological signatures such as the propagation of epileptiform discharges, using a multi-chambered tissue culture chip.

Computational neuronal correlation with enhanced synchronized activity in the basal ganglia and the slowing of thalamic theta and alpha rhythms in Parkinson's disease

The aim of this work is computationally to correlate the synchronized neuronal activity of basal ganglia and slowing in theta and alpha rhythms in electroencephalogram (EEG) signal in thalamic region in case of dopamine depletion and decrease of synaptic connections. The used network topology is a scale-free network with constant node degree. The dopamine-modulated type Izikhevich neuron model is used for modeling the striatal region, consisting of fast-spiking interneurons, D1 and D2 type dopamine expressing medium spiny neurons. On the other hand, the ordinary Izikhevich neuron model is used in the modeling of extrastriatal basal ganglia (BG) regions where globus pallidus (GP) subregion neurons have also dopamine-dependent parameters. The thalamic region of the network is mass modeled including inhibitory input from basal ganglia. Depending on the decrease of synaptic connections and dopamine level, the synchronization among basal ganglia neuron populations is investigated. The effect of synaptic delay on synchronization is also considered. It is observed that the decrease of dopamine neurotransmitter and decrease in the number of synaptic connections cause an increased synchronous activity in BG. Also, slowing in theta and alpha bands in thalamus EEG signals is observed. This shows the causal relation between synchronization and power shifting to lower frequency components in the case of neurodegenerative diseases such as Parkinson's disease (PD).

Cluster synchronization of delayed coupled neural networks: Delay-dependent distributed impulsive control

This paper investigates the issue of cluster synchronization (CS) for the coupled neural networks (CNNs) with time-varying delays via the delay-dependent distributed impulsive control. A new Halanay-like inequality, where delayed impulses are taken into consideration, is proposed. Based on the Lyapunov theory and the new differential inequality, sufficient conditions of CS for delayed CNNs with fixed and switching coupling topology are obtained, respectively. Moreover, delay-dependent distributed impulsive controllers with fixed or switching topology are designed thereby. Finally, we present a numerical example of CNNs with fixed or switching coupling to verify the effectiveness of our results, respectively.

Influence of Delayed Conductance on Neuronal Synchronization

In the brain, the excitation-inhibition balance prevents abnormal synchronous behavior. However, known synaptic conductance intensity can be insufficient to account for the undesired synchronization. Due to this fact, we consider time delay in excitatory and inhibitory conductances and study its effect on the neuronal synchronization. In this work, we build a neuronal network composed of adaptive integrate-and-fire neurons coupled by means of delayed conductances. We observe that the time delay in the excitatory and inhibitory conductivities can alter both the state of the collective behavior (synchronous or desynchronous) and its type (spike or burst). For the weak coupling regime, we find that synchronization appears associated with neurons behaving with extremes highest and lowest mean firing frequency, in contrast to when desynchronization is present when neurons do not exhibit extreme values for the firing frequency. Synchronization can also be characterized by neurons presenting either the highest or the lowest levels in the mean synaptic current. For the strong coupling, synchronous burst activities can occur for delays in the inhibitory conductivity. For approximately equal-length delays in the excitatory and inhibitory conductances, desynchronous spikes activities are identified for both weak and strong coupling regimes. Therefore, our results show that not only the conductance intensity, but also short delays in the inhibitory conductance are relevant to avoid abnormal neuronal synchronization.

Synchronous patterns and intermittency in a network induced by the rewiring of connections and coupling

The connection architecture plays an important role in the synchronization of networks, where the presence of local and nonlocal connection structures are found in many systems, such as the neural ones. Here, we consider a network composed of chaotic bursting oscillators coupled through a Watts-Strogatz-small-world topology. The influence of coupling strength and rewiring of connections is studied when the network topology is varied from regular to small-world to random. In this scenario, we show two distinct nonstationary transitions to phase synchronization: one induced by the increase in coupling strength and another resulting from the change from local connections to nonlocal ones. Besides this, there are regions in the parameter space where the network depicts a coexistence of different bursting frequencies where nonstationary zig-zag fronts are observed. Regarding the analyses, we consider two distinct methodological approaches: one based on the phase association to the bursting activity where the Kuramoto order parameter is used and another based on recurrence quantification analysis where just a time series of the network mean field is required.

Cluster burst synchronization in a scale-free network of inhibitory bursting neurons

We consider a scale-free network of inhibitory Hindmarsh-Rose (HR) bursting neurons, and make a computational study on coupling-induced cluster burst synchronization by varying the average coupling strength J 0 . For sufficiently small J 0 , non-cluster desynchronized states exist. However, when passing a critical point J c ∗ ( ≃ 0.16 ) , the whole population is segregated into 3 clusters via a constructive role of synaptic inhibition to stimulate dynamical clustering between individual burstings, and thus 3-cluster desynchronized states appear. As J 0 is further increased and passes a lower threshold J l ∗ ( ≃ 0.78 ) , a transition to 3-cluster burst synchronization occurs due to another constructive role of synaptic inhibition to favor population synchronization. In this case, HR neurons in each cluster make burstings every 3rd cycle of the instantaneous burst rate R w ( t ) of the whole population, and exhibit burst synchronization. However, as J 0 passes an intermediate threshold J m ∗ ( ≃ 5.2 ) , HR neurons fire burstings intermittently at a 4th cycle of R w ( t ) via burst skipping rather than at its 3rd cycle, and hence they begin to make intermittent hoppings between the 3 clusters. Due to such intermittent intercluster hoppings via burst skippings, the 3 clusters become broken up (i.e., the 3 clusters are integrated into a single one). However, in spite of such break-up (i.e., disappearance) of the 3-cluster states, (non-cluster) burst synchronization persists in the whole population, which is well visualized in the raster plot of burst onset times where bursting stripes (composed of burst onset times and indicating burst synchronization) appear successively. With further increase in J 0 , intercluster hoppings are intensified, and bursting stripes also become dispersed more and more due to a destructive role of synaptic inhibition to spoil the burst synchronization. Eventually, when passing a higher threshold J h ∗ ( ≃ 17.8 ) a transition to desynchronization occurs via complete overlap between the bursting stripes. Finally, we also investigate the effects of stochastic noise on both 3-cluster burst synchronization and intercluster hoppings.

Role of time scales and topology on the dynamics of complex networks

The interplay between time scales and structural properties of complex networks of nonlinear oscillators can generate many interesting phenomena, like amplitude death, cluster synchronization, frequency synchronization, etc. We study the emergence of such phenomena and their transitions by considering a complex network of dynamical systems in which a fraction of systems evolves on a slower time scale on the network. We report the transition to amplitude death for the whole network and the scaling near the transitions as the connectivity pattern changes. We also discuss the suppression and recovery of oscillations and the crossover behavior as the number of slow systems increases. By considering a scale free network of systems with multiple time scales, we study the role of heterogeneity in link structure on dynamical properties and the consequent critical behaviors. In this case with hubs made slow, our main results are the escape time statistics for loss of complete synchrony as the slowness spreads on the network and the self-organization of the whole network to a new frequency synchronized state. Our results have potential applications in biological, physical, and engineering networks consisting of heterogeneous oscillators.

Synchronization of Rulkov neuron networks coupled by excitatory and inhibitory chemical synapses

This paper takes into account a neuron network model in which the excitatory and the inhibitory Rulkov neurons interact each other through excitatory and inhibitory chemical coupling, respectively. Firstly, for two or more identical or non-identical Rulkov neurons, the existence conditions of the synchronization manifold of the fixed points are investigated, which have received less attention over the past decades. Secondly, the master stability equation of the arbitrarily connected neuron network under the existence conditions of the synchronization manifold is discussed. Thirdly, taking three identical Rulkov neurons as an example, some new results are presented: (1) topological structures that can make the synchronization manifold exist are given, (2) the stability of synchronization when different parameters change is discussed, and (3) the roles of the control parameters, the ratio, as well as the size of the coupling strength and sigmoid function are analyzed. Finally, for the chemical coupling between two non-identical neurons, the transversal system is given and the effect of two coupling strengths on synchronization is analyzed.

Neuron dynamics variability and anomalous phase synchronization of neural networks

Anomalous phase synchronization describes a synchronization phenomenon occurring even for the weakly coupled network and characterized by a non-monotonous dependence of the synchronization strength on the coupling strength. Its existence may support a theoretical framework to some neurological diseases, such as Parkinson's and some episodes of seizure behavior generated by epilepsy. Despite the success of controlling or suppressing the anomalous phase synchronization in neural networks applying external perturbations or inducing ambient changes, the origin of the anomalous phase synchronization as well as the mechanisms behind the suppression is not completely known. Here, we consider networks composed of N = 2000 coupled neurons in a small-world topology for two well known neuron models, namely, the Hodgkin-Huxley-like and the Hindmarsh-Rose models, both displaying the anomalous phase synchronization regime. We show that the anomalous phase synchronization may be related to the individual behavior of the coupled neurons; particularly, we identify a strong correlation between the behavior of the inter-bursting-intervals of the neurons, what we call neuron variability, to the ability of the network to depict anomalous phase synchronization. We corroborate the ideas showing that external perturbations or ambient parameter changes that eliminate anomalous phase synchronization and at the same time promote small changes in the individual dynamics of the neurons, such that an increasing individual variability of neurons implies a decrease of anomalous phase synchronization. Finally, we demonstrate that this effect can be quantified using a well known recurrence quantifier, the "determinism." Moreover, the results obtained by the determinism are based on only the mean field potential of the network, turning these measures more suitable to be used in experimental situations.

Burst synchronization in a scale-free neuronal network with inhibitory spike-timing-dependent plasticity

We are concerned about burst synchronization (BS), related to neural information processes in health and disease, in the Barabási-Albert scale-free network (SFN) composed of inhibitory bursting Hindmarsh-Rose neurons. This inhibitory neuronal population has adaptive dynamic synaptic strengths governed by the inhibitory spike-timing-dependent plasticity (iSTDP). In previous works without considering iSTDP, BS was found to appear in a range of noise intensities for fixed synaptic inhibition strengths. In contrast, in our present work, we take into consideration iSTDP and investigate its effect on BS by varying the noise intensity. Our new main result is to find occurrence of a Matthew effect in inhibitory synaptic plasticity: good BS gets better via LTD, while bad BS get worse via LTP. This kind of Matthew effect in inhibitory synaptic plasticity is in contrast to that in excitatory synaptic plasticity where good (bad) synchronization gets better (worse) via LTP (LTD). We note that, due to inhibition, the roles of LTD and LTP in inhibitory synaptic plasticity are reversed in comparison with those in excitatory synaptic plasticity. Moreover, emergences of LTD and LTP of synaptic inhibition strengths are intensively investigated via a microscopic method based on the distributions of time delays between the pre- and the post-synaptic burst onset times. Finally, in the presence of iSTDP we investigate the effects of network architecture on BS by varying the symmetric attachment degree l ∗ and the asymmetry parameter Δ l in the SFN.

Recurrence quantification analysis for the identification of burst phase synchronisation

In this work, we apply the spatial recurrence quantification analysis (RQA) to identify chaotic burst phase synchronisation in networks. We consider one neural network with small-world topology and another one composed of small-world subnetworks. The neuron dynamics is described by the Rulkov map, which is a two-dimensional map that has been used to model chaotic bursting neurons. We show that with the use of spatial RQA, it is possible to identify groups of synchronised neurons and determine their size. For the single network, we obtain an analytical expression for the spatial recurrence rate using a Gaussian approximation. In clustered networks, the spatial RQA allows the identification of phase synchronisation among neurons within and between the subnetworks. Our results imply that RQA can serve as a useful tool for studying phase synchronisation even in networks of networks.

Mean-field dynamics of a population of stochastic map neurons

We analyze the emergent regimes and the stimulus-response relationship of a population of noisy map neurons by means of a mean-field model, derived within the framework of cumulant approach complemented by the Gaussian closure hypothesis. It is demonstrated that the mean-field model can qualitatively account for stability and bifurcations of the exact system, capturing all the generic forms of collective behavior, including macroscopic excitability, subthreshold oscillations, periodic or chaotic spiking, and chaotic bursting dynamics. Apart from qualitative analogies, we find a substantial quantitative agreement between the exact and the approximate system, as reflected in matching of the parameter domains admitting the different dynamical regimes, as well as the characteristic properties of the associated time series. The effective model is further shown to reproduce with sufficient accuracy the phase response curves of the exact system and the assembly's response to external stimulation of finite amplitude and duration.

Effect of spike-timing-dependent plasticity on stochastic burst synchronization in a scale-free neuronal network

We consider an excitatory population of subthreshold Izhikevich neurons which cannot fire spontaneously without noise. As the coupling strength passes a threshold, individual neurons exhibit noise-induced burstings. This neuronal population has adaptive dynamic synaptic strengths governed by the spike-timing-dependent plasticity (STDP). However, STDP was not considered in previous works on stochastic burst synchronization (SBS) between noise-induced burstings of sub-threshold neurons. Here, we study the effect of additive STDP on SBS by varying the noise intensity D in the Barabási-Albert scale-free network (SFN). One of our main findings is a Matthew effect in synaptic plasticity which occurs due to a positive feedback process. Good burst synchronization (with higher bursting measure) gets better via long-term potentiation (LTP) of synaptic strengths, while bad burst synchronization (with lower bursting measure) gets worse via long-term depression (LTD). Consequently, a step-like rapid transition to SBS occurs by changing D, in contrast to a relatively smooth transition in the absence of STDP. We also investigate the effects of network architecture on SBS by varying the symmetric attachment degree [Formula: see text] and the asymmetry parameter [Formula: see text] in the SFN, and Matthew effects are also found to occur by varying [Formula: see text] and [Formula: see text]. Furthermore, emergences of LTP and LTD of synaptic strengths are investigated in details via our own microscopic methods based on both the distributions of time delays between the burst onset times of the pre- and the post-synaptic neurons and the pair-correlations between the pre- and the post-synaptic instantaneous individual burst rates (IIBRs). Finally, a multiplicative STDP case (depending on states) with soft bounds is also investigated in comparison with the additive STDP case (independent of states) with hard bounds. Due to the soft bounds, a Matthew effect with some quantitative differences is also found to occur for the case of multiplicative STDP.

Delay-induced locking in bursting neuronal networks

In this paper, the collective behaviors for ring structured bursting neuronal networks with electrical couplings and distance-dependent delays are studied. Each neuron is modeled by the Hindmarsh-Rose neuron. Through changing time delays between connected neurons, different spatiotemporal patterns are obtained. These patterns can be explained by calculating the ratios between the bursting period and the delay which exhibit clear locking relations. The holding and the failure of the lockings are investigated via bifurcation analysis. In particular, the bursting death phenomenon is observed for large coupling strengths and small time delays which is in fact the result of the partial amplitude death in the fast subsystem. These results indicate that the collective behaviors of bursting neurons critically depend on the bifurcation structure of individual ones and thus the variety of bifurcation types for bursting neurons may create diverse behaviors in similar neuronal ensembles.

A chaotic bursting-spiking transition in a pancreatic beta-cells system: Observation of an interior glucose-induced crisis

Nonlinear systems are commonly able to display abrupt qualitative changes (or transitions) in the dynamics. A particular type of these transitions occurs when the size of a chaotic attractor suddenly changes. In this article, we present such a transition through the observation of a chaotic interior crisis in the Deng bursting-spiking model for the glucose-induced electrical activity of pancreatic β-cells. To this chaos-chaos transition corresponds precisely the change between the bursting and spiking dynamics, which are central and key dynamical regimes that the Deng model is able to perform. We provide a description of the crisis mechanism at the bursting-spiking transition point in terms of time series variations and based on certain amplitudes of invariant intervals associated with return maps. Using symbolic dynamics, we are able to accurately compute the points of a curve representing the transition between the bursting and spiking regimes in a biophysical meaningfully parameter space. The analysis of the chaotic interior crisis is complemented by means of topological invariants with the computation of the topological entropy and the maximum Lyapunov exponent. Considering very recent developments in the literature, we construct analytical solutions triggering the bursting-spiking transition in the Deng model. This study provides an illustration of how an integrated approach, involving numerical evidences and theoretical reasoning within the theory of dynamical systems, can directly enhance our understanding of biophysically motivated models.

Dynamical responses to external stimuli for both cases of excitatory and inhibitory synchronization in a complex neuronal network

For studying how dynamical responses to external stimuli depend on the synaptic-coupling type, we consider two types of excitatory and inhibitory synchronization (i.e., synchronization via synaptic excitation and inhibition) in complex small-world networks of excitatory regular spiking (RS) pyramidal neurons and inhibitory fast spiking (FS) interneurons. For both cases of excitatory and inhibitory synchronization, effects of synaptic couplings on dynamical responses to external time-periodic stimuli S(t) (applied to a fraction of neurons) are investigated by varying the driving amplitude A of S(t). Stimulated neurons are phase-locked to external stimuli for both cases of excitatory and inhibitory couplings. On the other hand, the stimulation effect on non-stimulated neurons depends on the type of synaptic coupling. The external stimulus S(t) makes a constructive effect on excitatory non-stimulated RS neurons (i.e., it causes external phase lockings in the non-stimulated sub-population), while S(t) makes a destructive effect on inhibitory non-stimulated FS interneurons (i.e., it breaks up original inhibitory synchronization in the non-stimulated sub-population). As results of these different effects of S(t), the type and degree of dynamical response (e.g., synchronization enhancement or suppression), characterized by the dynamical response factor [Formula: see text] (given by the ratio of synchronization degree in the presence and absence of stimulus), are found to vary in a distinctly different way, depending on the synaptic-coupling type. Furthermore, we also measure the matching degree between the dynamics of the two sub-populations of stimulated and non-stimulated neurons in terms of a "cross-correlation" measure [Formula: see text]. With increasing A, based on [Formula: see text], we discuss the cross-correlations between the two sub-populations, affecting the dynamical responses to S(t).

Effects of partial time delays on phase synchronization in Watts-Strogatz small-world neuronal networks

In this paper, we study effects of partial time delays on phase synchronization in Watts-Strogatz small-world neuronal networks. Our focus is on the impact of two parameters, namely the time delay τ and the probability of partial time delay p_delay, whereby the latter determines the probability with which a connection between two neurons is delayed. Our research reveals that partial time delays significantly affect phase synchronization in this system. In particular, partial time delays can either enhance or decrease phase synchronization and induce synchronization transitions with changes in the mean firing rate of neurons, as well as induce switching between synchronized neurons with period-1 firing to synchronized neurons with period-2 firing. Moreover, in comparison to a neuronal network where all connections are delayed, we show that small partial time delay probabilities have especially different influences on phase synchronization of neuronal networks.

Emergence of ultrafast sparsely synchronized rhythms and their responses to external stimuli in an inhomogeneous small-world complex neuronal network

We consider an inhomogeneous small-world network (SWN) composed of inhibitory short-range (SR) and long-range (LR) interneurons, and investigate the effect of network architecture on emergence of synchronized brain rhythms by varying the fraction of LR interneurons p_long. The betweenness centralities of the LR and SR interneurons (characterizing the potentiality in controlling communication between other interneurons) are distinctly different. Hence, in view of the betweenness, SWNs we consider are inhomogeneous, unlike the "canonical" Watts-Strogatz SWN with nearly the same betweenness centralities. For small p_long, the load of communication traffic is much concentrated on a few LR interneurons. However, as p_long is increased, the number of LR connections (coming from LR interneurons) increases, and then the load of communication traffic is less concentrated on LR interneurons, which leads to better efficiency of global communication between interneurons. Sparsely synchronized rhythms are thus found to emerge when passing a small critical value p_long^(c)(≃0.16). The population frequency of the sparsely synchronized rhythm is ultrafast (higher than 100 Hz), while the mean firing rate of individual interneurons is much lower (∼30 Hz) due to stochastic and intermittent neural discharges. These dynamical behaviors in the inhomogeneous SWN are also compared with those in the homogeneous Watts-Strogatz SWN, in connection with their network topologies. Particularly, we note that the main difference between the two types of SWNs lies in the distribution of betweenness centralities. Unlike the case of the Watts-Strogatz SWN, dynamical responses to external stimuli vary depending on the type of stimulated interneurons in the inhomogeneous SWN. We consider two cases of external time-periodic stimuli applied to sub-populations of the LR and SR interneurons, respectively. Dynamical responses (such as synchronization suppression and enhancement) to these two cases of stimuli are studied and discussed in relation to the betweenness centralities of stimulated interneurons, representing the effectiveness for transfer of stimulation effect in the whole network.

Locking induced by distance-dependent delay in neuronal networks

In the present paper, the locking phenomenon induced by distance-dependent delay in ring structured neuronal networks is investigated, wherein each neuron is modeled by a FitzHugh-Nagumo neuron. Through increasing the element time delay, the different spatiotemporal patterns are observed. By calculating the interspike interval and its value that is divided by the delay of the nearest neurons, it is found that these patterns are actually the lockings between the period of spiking and the distance-dependent delay of the connected neurons. The lockings could also be revealed by the mean time lag of the neurons and in different connection topologies. Furthermore, the influences of the network size and the coupling strength are investigated, wherein the former seems to play a negligible role on these locking patterns; in contrast, too small coupling strengths will blur the boundaries of different patterns and too large ones may destroy the high ratio locking patterns. Finally, one may predict the locking order which determines the emergence order of the patterns in the networks.

Modeling of synchronization behavior of bursting neurons at nonlinearly coupled dynamical networks

Synchronization behaviors of bursting neurons coupled through electrical and dynamic chemical synapses are investigated. The Izhikevich model is used with random and small world network of bursting neurons. Various currents which consist of diffusive electrical and time-delayed dynamic chemical synapses are used in the simulations to investigate the influences of synaptic currents and couplings on synchronization behavior of bursting neurons. The effects of parameters, such as time delay, inhibitory synaptic strengths, and decay time on synchronization behavior are investigated. It is observed that in random networks with no delay, bursting synchrony is established with the electrical synapse alone, single spiking synchrony is observed with hybrid coupling. In small world network with no delay, periodic bursting behavior with multiple spikes is observed when only chemical and only electrical synapse exist. Single-spike and multiple-spike bursting are established with hybrid couplings. A decrease in the synchronization measure is observed with zero time delay, as the decay time is increased in random network. For synaptic delays which are above active phase period, synchronization measure increases with an increase in synaptic strength and time delay in small world network. However, in random network, it increases with only an increase in synaptic strength.

Effects of distance-dependent delay on small-world neuronal networks

We study firing behaviors and the transitions among them in small-world noisy neuronal networks with electrical synapses and information transmission delay. Each neuron is modeled by a two-dimensional Rulkov map neuron. The distance between neurons, which is a main source of the time delay, is taken into consideration. Through spatiotemporal patterns and interspike intervals as well as the interburst intervals, the collective behaviors are revealed. It is found that the networks switch from resting state into intermittent firing state under Gaussian noise excitation. Initially, noise-induced firing behaviors are disturbed by small time delays. Periodic firing behaviors with irregular zigzag patterns emerge with an increase of the delay and become progressively regular after a critical value is exceeded. More interestingly, in accordance with regular patterns, the spiking frequency doubles compared with the former stage for the spiking neuronal network. A growth of frequency persists for a larger delay and a transition to antiphase synchronization is observed. Furthermore, it is proved that these transitions are generic also for the bursting neuronal network and the FitzHugh-Nagumo neuronal network. We show these transitions due to the increase of time delay are robust to the noise strength, coupling strength, network size, and rewiring probability.

Development and function of human cerebral cortex neural networks from pluripotent stem cells in vitro

A key aspect of nervous system development, including that of the cerebral cortex, is the formation of higher-order neural networks. Developing neural networks undergo several phases with distinct activity patterns in vivo, which are thought to prune and fine-tune network connectivity. We report here that human pluripotent stem cell (hPSC)-derived cerebral cortex neurons form large-scale networks that reflect those found in the developing cerebral cortex in vivo. Synchronised oscillatory networks develop in a highly stereotyped pattern over several weeks in culture. An initial phase of increasing frequency of oscillations is followed by a phase of decreasing frequency, before giving rise to non-synchronous, ordered activity patterns. hPSC-derived cortical neural networks are excitatory, driven by activation of AMPA- and NMDA-type glutamate receptors, and can undergo NMDA-receptor-mediated plasticity. Investigating single neuron connectivity within PSC-derived cultures, using rabies-based trans-synaptic tracing, we found two broad classes of neuronal connectivity: most neurons have small numbers (<10) of presynaptic inputs, whereas a small set of hub-like neurons have large numbers of synaptic connections (>40). These data demonstrate that the formation of hPSC-derived cortical networks mimics in vivo cortical network development and function, demonstrating the utility of in vitro systems for mechanistic studies of human forebrain neural network biology.

Summary: Human PSC-derived cerebral cortex neurons form large-scale functional networks that change over time and mimic those found in the developing cerebral cortex in vivo.

Fast sparsely synchronized brain rhythms in a scale-free neural network

We consider a directed version of the Barabási-Albert scale-free network model with symmetric preferential attachment with the same in- and out-degrees and study the emergence of sparsely synchronized rhythms for a fixed attachment degree in an inhibitory population of fast-spiking Izhikevich interneurons. Fast sparsely synchronized rhythms with stochastic and intermittent neuronal discharges are found to appear for large values of J (synaptic inhibition strength) and D (noise intensity). For an intensive study we fix J at a sufficiently large value and investigate the population states by increasing D. For small D, full synchronization with the same population-rhythm frequency fp and mean firing rate (MFR) fi of individual neurons occurs, while for large D partial synchronization with fp>〈fi〉 (〈fi〉: ensemble-averaged MFR) appears due to intermittent discharge of individual neurons; in particular, the case of fp>4〈fi〉 is referred to as sparse synchronization. For the case of partial and sparse synchronization, MFRs of individual neurons vary depending on their degrees. As D passes a critical value D* (which is determined by employing an order parameter), a transition to unsynchronization occurs due to the destructive role of noise to spoil the pacing between sparse spikes. For D<D*, population synchronization emerges in the whole population because the spatial correlation length between the neuronal pairs covers the whole system. Furthermore, the degree of population synchronization is also measured in terms of two types of realistic statistical-mechanical measures. Only for the partial and sparse synchronization do contributions of individual neuronal dynamics to population synchronization change depending on their degrees, unlike in the case of full synchronization. Consequently, dynamics of individual neurons reveal the inhomogeneous network structure for the case of partial and sparse synchronization, which is in contrast to the case of statistically homogeneous random graphs and small-world networks. Finally, we investigate the effect of network architecture on sparse synchronization for fixed values of J and D in the following three cases: (1) variation in the degree of symmetric attachment, (2) asymmetric preferential attachment of new nodes with different in- and out-degrees, and (3) preferential attachment between pre-existing nodes (without addition of new nodes). In these three cases, both relation between network topology (e.g., average path length and betweenness centralization) and sparse synchronization and contributions of individual dynamics to the sparse synchronization are discussed.

Frequency-domain order parameters for the burst and spike synchronization transitions of bursting neurons

We are interested in characterization of synchronization transitions of bursting neurons in the frequency domain. Instantaneous population firing rate (IPFR) [Formula: see text], which is directly obtained from the raster plot of neural spikes, is often used as a realistic collective quantity describing population activities in both the computational and the experimental neuroscience. For the case of spiking neurons, a realistic time-domain order parameter, based on [Formula: see text], was introduced in our recent work to characterize the spike synchronization transition. Unlike the case of spiking neurons, the IPFR [Formula: see text] of bursting neurons exhibits population behaviors with both the slow bursting and the fast spiking timescales. For our aim, we decompose the IPFR [Formula: see text] into the instantaneous population bursting rate [Formula: see text] (describing the bursting behavior) and the instantaneous population spike rate [Formula: see text] (describing the spiking behavior) via frequency filtering, and extend the realistic order parameter to the case of bursting neurons. Thus, we develop the frequency-domain bursting and spiking order parameters which are just the bursting and spiking "coherence factors" [Formula: see text] and [Formula: see text] of the bursting and spiking peaks in the power spectral densities of [Formula: see text] and [Formula: see text] (i.e., "signal to noise" ratio of the spectral peak height and its relative width). Through calculation of [Formula: see text] and [Formula: see text], we obtain the bursting and spiking thresholds beyond which the burst and spike synchronizations break up, respectively. Consequently, it is shown in explicit examples that the frequency-domain bursting and spiking order parameters may be usefully used for characterization of the bursting and the spiking transitions, respectively.

Noise-induced burst and spike synchronizations in an inhibitory small-world network of subthreshold bursting neurons

We are interested in noise-induced firings of subthreshold neurons which may be used for encoding environmental stimuli. Noise-induced population synchronization was previously studied only for the case of global coupling, unlike the case of subthreshold spiking neurons. Hence, we investigate the effect of complex network architecture on noise-induced synchronization in an inhibitory population of subthreshold bursting Hindmarsh-Rose neurons. For modeling complex synaptic connectivity, we consider the Watts-Strogatz small-world network which interpolates between regular lattice and random network via rewiring, and investigate the effect of small-world connectivity on emergence of noise-induced population synchronization. Thus, noise-induced burst synchronization (synchrony on the slow bursting time scale) and spike synchronization (synchrony on the fast spike time scale) are found to appear in a synchronized region of the [Formula: see text]-[Formula: see text] plane ([Formula: see text]: synaptic inhibition strength and [Formula: see text]: noise intensity). As the rewiring probability [Formula: see text] is decreased from 1 (random network) to 0 (regular lattice), the region of spike synchronization shrinks rapidly in the [Formula: see text]-[Formula: see text] plane, while the region of the burst synchronization decreases slowly. We separate the slow bursting and the fast spiking time scales via frequency filtering, and characterize the noise-induced burst and spike synchronizations by employing realistic order parameters and statistical-mechanical measures introduced in our recent work. Thus, the bursting and spiking thresholds for the burst and spike synchronization transitions are determined in terms of the bursting and spiking order parameters, respectively. Furthermore, we also measure the degrees of burst and spike synchronizations in terms of the statistical-mechanical bursting and spiking measures, respectively.

Phase synchronization of coupled bursting neurons and the generalized Kuramoto model

Bursting neurons fire rapid sequences of action potential spikes followed by a quiescent period. The basic dynamical mechanism of bursting is the slow currents that modulate a fast spiking activity caused by rapid ionic currents. Minimal models of bursting neurons must include both effects. We considered one of these models and its relation with a generalized Kuramoto model, thanks to the definition of a geometrical phase for bursting and a corresponding frequency. We considered neuronal networks with different connection topologies and investigated the transition from a non-synchronized to a partially phase-synchronized state as the coupling strength is varied. The numerically determined critical coupling strength value for this transition to occur is compared with theoretical results valid for the generalized Kuramoto model.

ADAPTATION-DEPENDENT SYNCHRONIZATION TRANSITIONS AND BURST GENERATIONS IN ELECTRICALLY COUPLED NEURAL NETWORKS

A typical feature of neurons is their ability to encode neural information dynamically through spike frequency adaptation (SFA). Previous studies of SFA on neuronal synchronization were mainly concentrated on the correlated firing between neuron pairs, while the synchronization of neuron populations in the presence of SFA is still unclear. In this study, the influence of SFA on the population synchronization of neurons was numerically explored in electrically coupled networks, with regular, small-world, and random connectivity, respectively. The simulation results indicate that cross-correlation indices decrease significantly when the neurons have adaptation compared with those of nonadapting neurons, similar to previous experimental observations. However, the synchronous activity of population neurons exhibits a rather complex adaptation-dependent manner. Specifically, synchronization strength of neuron populations changes nonmonotonically, depending on the degree of adaptation. In addition, single neurons in the networks can switch from regular spiking to bursting with the increase of adaptation degree. Furthermore, the connection probability among neurons exhibits significant influence on the population synchronous activity, but has little effect on the burst generation of single neurons. Accordingly, the results may suggest that synchronous activity and burst firing of population neurons are both adaptation-dependent.

On the capabilities and computational costs of neuron models

We review the Hodgkin-Huxley, Izhikevich, and leaky integrate-and-fire neuron models in regular spiking modes solved with the forward Euler, fourth-order Runge-Kutta, and exponential Euler methods and determine the necessary time steps and corresponding computational costs required to make the solutions accurate. We conclude that the leaky integrate-and-fire needs the least number of computations, and that the Hodgkin-Huxley and Izhikevich models are comparable in computational cost.

Bursting frequency versus phase synchronization in time-delayed neuron networks

We investigate the dependence of the average bursting frequency on time delay for neuron networks with randomly distributed time-delayed chemical synapses. The result is compared with the corresponding curve for the phase synchronization and it turns out that, in some intervals, these have a very similar shape and appear as almost mirror images of each other. We have analyzed both the map-based chaotic Rulkov model and the continuous Hindmarsh-Rose model, yielding the same conclusions. In order to gain further insight, we also analyzed time-delayed Kuramoto models displaying an overall behavior similar to that observed on the neuron network models. For the Kuramoto models, we were able to derive analytical formulas providing an implicit functional relationship between the mean frequency and the phase synchronization. These formulas suggest a strong dependence between those two measures, which could explain the similarities in shape between the curves.

Control of bursting synchronization in networks of Hodgkin-Huxley-type neurons with chemical synapses

Thermally sensitive neurons present bursting activity for certain temperature ranges, characterized by fast repetitive spiking of action potential followed by a short quiescent period. Synchronization of bursting activity is possible in networks of coupled neurons, and it is sometimes an undesirable feature. Control procedures can suppress totally or partially this collective behavior, with potential applications in deep-brain stimulation techniques. We investigate the control of bursting synchronization in small-world networks of Hodgkin-Huxley-type thermally sensitive neurons with chemical synapses through two different strategies. One is the application of an external time-periodic electrical signal and another consists of a time-delayed feedback signal. We consider the effectiveness of both strategies in terms of protocols of applications suitable to be applied by pacemakers.

Suppression of bursting synchronization in clustered scale-free (rich-club) neuronal networks

Functional brain networks are composed of cortical areas that are anatomically and functionally connected. One of the cortical networks for which more information is available in the literature is the cat cerebral cortex. Statistical analyses of the latter suggest that its structure can be described as a clustered network, in which each cluster is a scale-free network possessing highly connected hubs. Those hubs are, on their hand, connected together in a strong fashion ("rich-club" network). We have built a clustered scale-free network inspired in the cat cortex structure so as to study their dynamical properties. In this article, we focus on the synchronization of bursting activity of the cortical areas and how it can be suppressed by means of neuron deactivation through suitably applied light pulses. We show that it is possible to effectively suppress bursting synchronization by acting on a single, yet suitably chosen neuron, as long as it is highly connected, thanks to the "rich-club" structure of the network.

Collective almost synchronisation in complex networks

This work introduces the phenomenon of Collective Almost Synchronisation (CAS), which describes a universal way of how patterns can appear in complex networks for small coupling strengths. The CAS phenomenon appears due to the existence of an approximately constant local mean field and is characterised by having nodes with trajectories evolving around periodic stable orbits. Common notion based on statistical knowledge would lead one to interpret the appearance of a local constant mean field as a consequence of the fact that the behaviour of each node is not correlated to the behaviours of the others. Contrary to this common notion, we show that various well known weaker forms of synchronisation (almost, time-lag, phase synchronisation, and generalised synchronisation) appear as a result of the onset of an almost constant local mean field. If the memory is formed in a brain by minimising the coupling strength among neurons and maximising the number of possible patterns, then the CAS phenomenon is a plausible explanation for it.

Onset of chaotic phase synchronization in complex networks of coupled heterogeneous oscillators

Existing studies on network synchronization focused on complex networks possessing either identical or nonidentical but simple nodal dynamics. We consider networks of both complex topologies and heterogeneous but chaotic oscillators, and investigate the onset of global phase synchronization. Based on a heuristic analysis and by developing an efficient numerical procedure to detect the onset of phase synchronization, we uncover a general scaling law, revealing that chaotic phase synchronization can be facilitated by making the network more densely connected. Our methodology can find applications in probing the fundamental network dynamics in realistic situations, where both complex topology and complicated, heterogeneous nodal dynamics are expected.

Phase synchronization of bursting neurons in clustered small-world networks

We investigate the collective dynamics of bursting neurons on clustered networks. The clustered network model is composed of subnetworks, each of them presenting the so-called small-world property. This model can also be regarded as a network of networks. In each subnetwork a neuron is connected to other ones with regular as well as random connections, the latter with a given intracluster probability. Moreover, in a given subnetwork each neuron has an intercluster probability to be connected to the other subnetworks. The local neuron dynamics has two time scales (fast and slow) and is modeled by a two-dimensional map. In such small-world network the neuron parameters are chosen to be slightly different such that, if the coupling strength is large enough, there may be synchronization of the bursting (slow) activity. We give bounds for the critical coupling strength to obtain global burst synchronization in terms of the network structure, that is, the probabilities of intracluster and intercluster connections. We find that, as the heterogeneity in the network is reduced, the network global synchronizability is improved. We show that the transitions to global synchrony may be abrupt or smooth depending on the intercluster probability.

A general fractional-order dynamical network: synchronization behavior and state tuning

A general fractional-order dynamical network model for synchronization behavior is proposed. Different from previous integer-order dynamical networks, the model is made up of coupled units described by fractional differential equations, where the connections between individual units are nondiffusive and nonlinear. We show that the synchronous behavior of such a network cannot only occur, but also be dramatically different from the behavior of its constituent units. In particular, we find that simple behavior can emerge as synchronized dynamics although the isolated units evolve chaotically. Conversely, individually simple units can display chaotic attractors when the network synchronizes. We also present an easily checked criterion for synchronization depending only on the eigenvalues distribution of a decomposition matrix and the fractional orders. The analytic results are complemented with numerical simulations for two networks whose nodes are governed by fractional-order Lorenz dynamics and fractional-order Rössler dynamics, respectively.

Vibrational resonance in excitable neuronal systems

In this paper, we investigate the effect of a high-frequency driving on the dynamical response of excitable neuronal systems to a subthreshold low-frequency signal by numerical simulation. We demonstrate the occurrence of vibrational resonance in spatially extended neuronal networks. Different network topologies from single small-world networks to modular networks of small-world subnetworks are considered. It is shown that an optimal amplitude of high-frequency driving enhances the response of neuron populations to a low-frequency signal. This effect of vibrational resonance of neuronal systems depends extensively on the network structure and parameters, such as the coupling strength between neurons, network size, and rewiring probability of single small-world networks, as well as the number of links between different subnetworks and the number of subnetworks in the modular networks. All these parameters play a key role in determining the ability of the network to enhance the outreach of the localized subthreshold low-frequency signal. Considering that two-frequency signals are ubiquity in brain dynamics, we expect the presented results could have important implications for the weak signal detection and information propagation across neuronal systems.

Chaotic phase synchronization in a modular neuronal network of small-world subnetworks

We investigate the onset of chaotic phase synchronization of bursting oscillators in a modular neuronal network of small-world subnetworks. A transition to mutual phase synchronization takes place on the bursting time scale of coupled oscillators, while on the spiking time scale, they behave asynchronously. It is shown that this bursting synchronization transition can be induced not only by the variations of inter- and intra-coupling strengths but also by changing the probability of random links between different subnetworks. We also analyze the effect of external chaotic phase synchronization of bursting behavior in this clustered network by an external time-periodic signal applied to a single neuron. Simulation results demonstrate a frequency locking tongue in the driving parameter plane, where bursting synchronization is maintained, even with the external driving. The width of this synchronization region increases with the signal amplitude and the number of driven neurons but decreases rapidly with the network size. Considering that the synchronization of bursting neurons is thought to play a key role in some pathological conditions, the presented results could have important implications for the role of externally applied driving signal in controlling bursting activity in neuronal ensembles.

Stochastic resonance on a modular neuronal network of small-world subnetworks with a subthreshold pacemaker

We study the phenomenon of stochastic resonance on a modular neuronal network consisting of several small-world subnetworks with a subthreshold periodic pacemaker. Numerical results show that the correlation between the pacemaker frequency and the dynamical response of the network is resonantly dependent on the intensity of additive spatiotemporal noise. This effect of pacemaker-driven stochastic resonance of the system depends extensively on the local and the global network structure, such as the intra- and inter-coupling strengths, rewiring probability of individual small-world subnetwork, the number of links between different subnetworks, and the number of subnetworks. All these parameters play a key role in determining the ability of the network to enhance the noise-induced outreach of the localized subthreshold pacemaker, and only they bounded to a rather sharp interval of values warrant the emergence of the pronounced stochastic resonance phenomenon. Considering the rather important role of pacemakers in real-life, the presented results could have important implications for many biological processes that rely on an effective pacemaker for their proper functioning.

Chaotic phase synchronization in small-world networks of bursting neurons

We investigate the chaotic phase synchronization in a system of coupled bursting neurons in small-world networks. A transition to mutual phase synchronization takes place on the bursting time scale of coupled oscillators, while on the spiking time scale, they behave asynchronously. It is shown that phase synchronization is largely facilitated by a large fraction of shortcuts, but saturates when it exceeds a critical value. We also study the external chaotic phase synchronization of bursting oscillators in the small-world network by a periodic driving signal applied to a single neuron. It is demonstrated that there exists an optimal small-world topology, resulting in the largest peak value of frequency locking interval in the parameter plane, where bursting synchronization is maintained, even with the external driving. The width of this interval increases with the driving amplitude, but decrease rapidly with the network size. We infer that the externally applied driving parameters outside the frequency locking region can effectively suppress pathologically synchronized rhythms of bursting neurons in the brain.

Delay-induced intermittent transition of synchronization in neuronal networks with hybrid synapses

We study the dependence of synchronization transitions in scale-free networks of bursting neurons with hybrid synapses on the information transmission delay and the probability of inhibitory synapses. It is shown that, irrespective of the probability of inhibitory synapses, the delay always plays a subtle role during synchronization transition of the scale-free neuronal networks. In particular, regions of irregular and regular propagating excitatory fronts appear intermittently as the delay increases. These delay-induced synchronization transitions are manifested as well-expressed minima in the measure for spatiotemporal synchrony. In addition, it is found that, for smaller and larger probability of inhibitory synapses, intermittent synchronization transition is relatively profound, while for the moderate probability of inhibitory synapses, synchronization transition seems less profound. More interestingly, it is found that as the probability of inhibitory synapses is large, regions of synchronization are upscattering.