Regularization of synchronized chaotic bursts

Unfolding the distribution of periodicity regions and diversity of chaotic attractors in the Chialvo neuron map

We performed an exhaustive numerical analysis of the two-dimensional Chialvo map by obtaining the parameter planes based on the computation of periodicities and Lyapunov exponents. Our results allowed us to determine the different regions of dynamical behavior, identify regularities in the distribution of periodicities in regions indicating regular behavior, find some pseudofractal structures, identify regions such as the "eyes of chaos" similar to those obtained in parameter planes of continuous systems, and, finally, characterize the statistical properties of chaotic attractors leading to possible hyperchaotic behavior.

Emergence of complex oscillatory dynamics in the neuronal networks with long activity time of inhibitory synapses

The brain displays complex dynamics, including collective oscillations, and extensive research has been conducted to understand their generation. However, our understanding of how biological constraints influence these oscillations is incomplete. This study investigates the essential properties of neuronal networks needed to generate oscillations resembling those in the brain. A simple discrete-time model of interconnected excitable elements is developed, capable of closely resembling the complex oscillations observed in biological neural networks. In the model, synaptic connections remain active for a duration exceeding individual neuron activity. We show that the inhibitory synapses must exhibit longer activity than excitatory synapses to produce a diverse range of the dynamical states, including biologically plausible oscillations. Upon meeting this condition, the transition between different dynamical states can be controlled by external stochastic input to the neurons. The study provides a comprehensive explanation for the emergence of distinct dynamical states in neural networks based on specific parameters.

Analysis of dynamics of a map-based neuron model via Lorenz maps

Modeling nerve cells can facilitate formulating hypotheses about their real behavior and improve understanding of their functioning. In this paper, we study a discrete neuron model introduced by Courbage et al. [Chaos 17, 043109 (2007)], where the originally piecewise linear function defining voltage dynamics is replaced by a cubic polynomial, with an additional parameter responsible for varying the slope. Showing that on a large subset of the multidimensional parameter space, the return map of the voltage dynamics is an expanding Lorenz map, we analyze both chaotic and periodic behavior of the system and describe the complexity of spiking patterns fired by a neuron. This is achieved by using and extending some results from the theory of Lorenz-like and expanding Lorenz mappings.

Complete synchronization of three-layer Rulkov neuron network coupled by electrical and chemical synapses

This paper analyzes the complete synchronization of a three-layer Rulkov neuron network model connected by electrical synapses in the same layers and chemical synapses between adjacent layers. The outer coupling matrix of the network is not Laplacian as in linear coupling networks. We develop the master stability function method, in which the invariant manifold of the master stability equations (MSEs) does not correspond to the zero eigenvalues of the connection matrix. After giving the existence conditions of the synchronization manifold about the nonlinear chemical coupling, we investigate the dynamics of the synchronization manifold, which will be identical to that of a synchronous network by fixing the same parameters and initial values. The waveforms show that the transient chaotic windows and the transient approximate periodic windows with increased or decreased periods occur alternatively before asymptotic behaviors. Furthermore, the Lyapunov exponents of the MSEs indicate that the network with a periodic synchronization manifold can achieve complete synchronization, while the network with a chaotic synchronization manifold can not. Finally, we simulate the effects of small perturbations on the asymptotic regimes and the evolution routes for the synchronous periodic and the non-synchronous chaotic network.

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.

Adaptive exponential integrate-and-fire model with fractal extension

The description of neuronal activity has been of great importance in neuroscience. In this field, mathematical models are useful to describe the electrophysical behavior of neurons. One successful model used for this purpose is the Adaptive Exponential Integrate-and-Fire (Adex), which is composed of two ordinary differential equations. Usually, this model is considered in the standard formulation, i.e., with integer order derivatives. In this work, we propose and study the fractal extension of Adex model, which in simple terms corresponds to replacing the integer derivative by non-integer. As non-integer operators, we choose the fractal derivatives. We explore the effects of equal and different orders of fractal derivatives in the firing patterns and mean frequency of the neuron described by the Adex model. Previous results suggest that fractal derivatives can provide a more realistic representation due to the fact that the standard operators are generalized. Our findings show that the fractal order influences the inter-spike intervals and changes the mean firing frequency. In addition, the firing patterns depend not only on the neuronal parameters but also on the order of respective fractal operators. As our main conclusion, the fractal order below the unit value increases the influence of the adaptation mechanism in the spike firing patterns.

Expecting the Unexpected: Entropy and Multifractal Systems in Finance

Entropy serves as a measure of chaos in systems by representing the average rate of information loss about a phase point's position on the attractor. When dealing with a multifractal system, a single exponent cannot fully describe its dynamics, necessitating a continuous spectrum of exponents, known as the singularity spectrum. From an investor's point of view, a rise in entropy is a signal of abnormal and possibly negative returns. This means he has to expect the unexpected and prepare for it. To explore this, we analyse the New York Stock Exchange (NYSE) U.S. Index as well as its constituents. Through this examination, we assess their multifractal characteristics and identify market conditions (bearish/bullish markets) using entropy, an effective method for recognizing fluctuating fractal markets. Our findings challenge conventional beliefs by demonstrating that price declines lead to increased entropy, contrary to some studies in the literature that suggest that reduced entropy in market crises implies more determinism. Instead, we propose that bear markets are likely to exhibit higher entropy, indicating a greater chance of unexpected extreme events. Moreover, our study reveals a power-law behaviour and indicates the absence of variance.

Both electrical and metabolic coupling shape the collective multimodal activity and functional connectivity patterns in beta cell collectives: A computational model perspective

Pancreatic beta cells are coupled excitable oscillators that synchronize their activity via different communication pathways. Their oscillatory activity manifests itself on multiple timescales and consists of bursting electrical activity, subsequent oscillations in the intracellular Ca^{2+}, as well as oscillations in metabolism and exocytosis. The coordination of the intricate activity on the multicellular level plays a key role in the regulation of physiological pulsatile insulin secretion and is incompletely understood. In this paper, we investigate theoretically the principles that give rise to the synchronized activity of beta cell populations by building up a phenomenological multicellular model that incorporates the basic features of beta cell dynamics. Specifically, the model is composed of coupled slow and fast oscillatory units that reflect metabolic processes and electrical activity, respectively. Using a realistic description of the intercellular interactions, we study how the combination of electrical and metabolic coupling generates collective rhythmicity and shapes functional beta cell networks. It turns out that while electrical coupling solely can synchronize the responses, the addition of metabolic interactions further enhances coordination, the spatial range of interactions increases the number of connections in the functional beta cell networks, and ensures a better consistency with experimental findings. Moreover, our computational results provide additional insights into the relationship between beta cell heterogeneity, their activity profiles, and functional connectivity, supplementing thereby recent experimental results on endocrine networks.

Reliability and robustness of oscillations in some slow-fast chaotic systems

A variety of nonlinear models of biological systems generate complex chaotic behaviors that contrast with biological homeostasis, the observation that many biological systems prove remarkably robust in the face of changing external or internal conditions. Motivated by the subtle dynamics of cell activity in a crustacean central pattern generator (CPG), this paper proposes a refinement of the notion of chaos that reconciles homeostasis and chaos in systems with multiple timescales. We show that systems displaying relaxation cycles while going through chaotic attractors generate chaotic dynamics that are regular at macroscopic timescales and are, thus, consistent with physiological function. We further show that this relative regularity may break down through global bifurcations of chaotic attractors such as crises, beyond which the system may also generate erratic activity at slow timescales. We analyze these phenomena in detail in the chaotic Rulkov map, a classical neuron model known to exhibit a variety of chaotic spike patterns. This leads us to propose that the passage of slow relaxation cycles through a chaotic attractor crisis is a robust, general mechanism for the transition between such dynamics. We validate this numerically in three other models: a simple model of the crustacean CPG neural network, a discrete cubic map, and a continuous flow.

Body as First Teacher: The Role of Rhythmic Visceral Dynamics in Early Cognitive Development

Embodied cognition-the idea that mental states and processes should be understood in relation to one's bodily constitution and interactions with the world-remains a controversial topic within cognitive science. Recently, however, increasing interest in predictive processing theories among proponents and critics of embodiment alike has raised hopes of a reconciliation. This article sets out to appraise the unificatory potential of predictive processing, focusing in particular on embodied formulations of active inference. Our analysis suggests that most active-inference accounts invoke weak, potentially trivial conceptions of embodiment; those making stronger claims do so independently of the theoretical commitments of the active-inference framework. We argue that a more compelling version of embodied active inference can be motivated by adopting a diachronic perspective on the way rhythmic physiological activity shapes neural development in utero. According to this visceral afferent training hypothesis, early-emerging physiological processes are essential not only for supporting the biophysical development of neural structures but also for configuring the cognitive architecture those structures entail. Focusing in particular on the cardiovascular system, we propose three candidate mechanisms through which visceral afferent training might operate: (a) activity-dependent neuronal development, (b) periodic signal modeling, and (c) oscillatory network coordination.

Topological-numerical analysis of a two-dimensional discrete neuron model

We conduct computer-assisted analysis of a two-dimensional model of a neuron introduced by Chialvo in 1995 [Chaos, Solitons Fractals 5, 461-479]. We apply the method of rigorous analysis of global dynamics based on a set-oriented topological approach, introduced by Arai et al. in 2009 [SIAM J. Appl. Dyn. Syst. 8, 757-789] and improved and expanded afterward. Additionally, we introduce a new algorithm to analyze the return times inside a chain recurrent set. Based on this analysis, together with the information on the size of the chain recurrent set, we develop a new method that allows one to determine subsets of parameters for which chaotic dynamics may appear. This approach can be applied to a variety of dynamical systems, and we discuss some of its practical aspects.

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.

Variations of the spontaneous electrical activities of the neuronal networks imposed by the exposure of electromagnetic radiations using computational map-based modeling

The interaction between neurons in a neuronal network develops spontaneous electrical activities. But the effects of electromagnetic radiation on these activities have not yet been well explored. In this study, a ring of three coupled 1-dimensional Rulkov neurons and the generated electromagnetic field (EMF) are considered to investigate how the spontaneous activities might change regarding the EMF exposure. By employing the bifurcation analysis and time series, a comprehensive view of neuronal behavioral changes due to electromagnetic inductions is provided. The main findings of this study are as follows: 1) When a neuronal network is showing a spontaneous chaotic firing manner (without any external stimuli), a generated magnetic field inhibits this type of behavior. In fact, EMF completely eliminated the chaotic intrinsic behaviors of the neuronal loop. 2) When the network is exhibiting regular period-3 spiking patterns, the generated magnetic field changes its firing pattern to chaotic spiking, which is similar to epileptic seizures. 3) With weak synaptic connections, electromagnetic radiation inhibits and suppresses neuronal activities. 4) If the external magnetic flux has a high amplitude, it can change the shape of the induction current according to its shape 5) when there are weak synaptic connections in the network, a high-frequency external magnetic flux engenders high-frequency fluctuations in the membrane voltages. On the whole, electromagnetic radiation changes the pattern of the spontaneous activities of neuronal networks in the brain according to synaptic strengths and initial states of the neurons.

Exploiting deterministic features in apparently stochastic data

Many processes in nature are the result of many coupled individual subsystems (like population dynamics or neurosystems). Not always such systems exhibit simple stable behaviors that in the past science has mostly focused on. Often, these systems are characterized by bursts of seemingly stochastic activity, interrupted by quieter periods. The hypothesis is that the presence of a strong deterministic ingredient is often obscured by the stochastic features. We test this by modeling classically stochastic considered real-world data from both, the stochastic as well as the deterministic approaches to find that the deterministic approach's results level with those from the stochastic side. Moreover, the deterministic approach is shown to reveal the full dynamical systems landscape, which can be exploited for steering the dynamics into a desired regime.

Rulkov neural network coupled with discrete memristors

The features of memristive-coupled neural networks have been studied extensively in the continuous field. However, the particularities of the discrete domain are rarely mentioned. This paper constructs a discrete memristor with sine-type conductance and applies the discrete memristor to coupling the Rulkov neuron maps for the first time. The properties of the proposed memristive-coupled bi-neuron Rulkov map and multi-neuron Rulkov neural network are probed. In order to better characterize the discrete system, many numerical simulation methods are used. Such as the normalized mean synchronization error, bifurcation diagrams, phase portraits, spatiotemporal patterns and so on. Numerical studies have shown that in discrete memristor-coupled neural networks, both parameters and coupling factors affect the dynamics of the system, resulting in complex and interesting behavioural changes.

Financial markets' deterministic aspects modeled by a low-dimensional equation

We ask whether empirical finance market data (Financial Stress Index, swap and equity, emerging and developed, corporate and government, short and long maturity), with their recently observed alternations between calm periods and financial turmoil, could be described by a low-dimensional deterministic model, or whether this requests a stochastic approach. We find that a deterministic model performs at least as well as one of the best stochastic models, but may offer additional insight into the essential mechanisms that drive financial markets.

Intermittent evolution routes to the periodic or the chaotic orbits in Rulkov map

This paper concerns the intermittent evolution routes to the asymptotic regimes in the Rulkov map. That is, the windows with transient approximate periodic and transient chaotic behaviors occur alternatively before the system reaches the periodic or the chaotic orbits. Meanwhile, the evolution routes to chaotic orbits can be classified into different types according to the windows before reaching asymptotic chaotic states. In addition, the initial values can be regarded as a key factor affecting the asymptotic behaviors and the evolution routes. The effects of the initial values are given by parameter planes, bifurcation diagrams, and waveforms. In order to investigate whether the intermittent evolution routes can be learned by machine learning, some experiments are given to understanding the differences between the trajectories of the Rulkov map generated by the numerical simulations and predicted by the neural networks. These results show that there is about 60% accuracy rate of successfully predicting both the evolution routes and the asymptotic period-3 orbits using a three-layer feedforward neural network, while the bifurcation diagrams can be reconstructed using reservoir computing except a few parameter conditions.

Suppression of chaotic bursting synchronization in clustered scale-free networks by an external feedback signal

Oscillatory activities in the brain, detected by electroencephalograms, have identified synchronization patterns. These synchronized activities in neurons are related to cognitive processes. Additionally, experimental research studies on neuronal rhythms have shown synchronous oscillations in brain disorders. Mathematical modeling of networks has been used to mimic these neuronal synchronizations. Actually, networks with scale-free properties were identified in some regions of the cortex. In this work, to investigate these brain synchronizations, we focus on neuronal synchronization in a network with coupled scale-free networks. The networks are connected according to a topological organization in the structural cortical regions of the human brain. The neuronal dynamic is given by the Rulkov model, which is a two-dimensional iterated map. The Rulkov neuron can generate quiescence, tonic spiking, and bursting. Depending on the parameters, we identify synchronous behavior among the neurons in the clustered networks. In this work, we aim to suppress the neuronal burst synchronization by the application of an external perturbation as a function of the mean-field of membrane potential. We found that the method we used to suppress synchronization presents better results when compared to the time-delayed feedback method when applied to the same model of the neuronal network.

Synchronization in populations of electrochemical bursting oscillators with chaotic slow dynamics

We investigate the synchronization of coupled electrochemical bursting oscillators using the electrodissolution of iron in sulfuric acid. The dynamics of a single oscillator consisted of slow chaotic oscillations interrupted by a burst of fast spiking, generating a multiple time-scale dynamical system. A wavelet analysis first decomposed the time series data from each oscillator into a fast and a slow component, and the corresponding phases were also obtained. The phase synchronization of the fast and slow dynamics was analyzed as a function of electrical coupling imposed by an external coupling resistance. For two oscillators, a progressive transition was observed: With increasing coupling strength, first, the fast bursting intervals overlapped, which was followed by synchronization of the fast spiking, and finally, the slow chaotic oscillations synchronized. With a population of globally coupled 25 oscillators, the coupling eliminated the fast dynamics, and only the synchronization of the slow dynamics can be observed. The results demonstrated the complexities of synchronization with bursting oscillations that could be useful in other systems with multiple time-scale dynamics, in particular, in neuronal networks.

A tale of two rhythms: Locked clocks and chaos in biology

The fundamental mechanisms that control and regulate biological organisms exhibit a surprising level of complexity. Oscillators are perhaps the simplest motifs that produce time-varying dynamics and are ubiquitous in biological systems. It is also known that such biological oscillators interact with each other-for instance, circadian oscillators affect the cell cycle, and somitogenesis clock proteins in adjacent cells affect each other in developing embryos. Therefore, it is vital to understand the effects that can emerge from non-linear interaction between oscillations. Here, we show how oscillations typically arise in biology and take the reader on a tour through the great variety in dynamics that can emerge even from a single pair of coupled oscillators. We explain how chaotic dynamics can emerge and outline the methods of detecting this in experimental time traces. Finally, we discuss the potential role of such complex dynamical features in biological systems.

Pattern selection in coupled neurons under high-low frequency electric field

Biological neurons exposed to an external electric field can induce polarization and charge fluctuation. Indeed, the exchange of calcium, sodium and potassium ions across the cell membrane can induce an electric field as a result of a time-varying electromagnetic field set up. This field could further modulate the cell electrical activity by inducing multiple firing modes. Based on the physical law of electric field, an improved model which includes an additive electrical field variable is constructed for a network of neurons to study wave propagation and mode transition by exploring the longtime dynamics of slightly perturbed plane waves in the network. Wave pattern and mode transition dependence on the different parameters of external electric field are discussed. It is found that the plane wave propagating in the coupled system breaks down to localized structures under the activation of modulational instability. The network under high-low external electric field supports bursting synchronization. This could be a fruitful avenue to discern the occurrence of paroxysmal epilepsy.

FRACTAL DIMENSION ESTIMATION WITH PERSISTENT HOMOLOGY: A COMPARATIVE STUDY

We propose that the recently defined persistent homology dimensions are a practical tool for fractal dimension estimation of point samples. We implement an algorithm to estimate the persistent homology dimension, and compare its performance to classical methods to compute the correlation and box-counting dimensions in examples of self-similar fractals, chaotic attractors, and an empirical dataset. The performance of the 0-dimensional persistent homology dimension is comparable to that of the correlation dimension, and better than box-counting.

Proper spatial heterogeneities expand the regime of scale-free behavior in a lattice of excitable elements

Signatures of criticality, such as power law scaling of observables, have been empirically found in a plethora of real-life settings, including biological systems. The presence of critical states is believed to have many functional advantages and is associated with optimal operational abilities. Typically, critical dynamics arises in the proximity of phase transition points between absorbing disordered states (subcriticality) and ordered active regimes (supercriticality) and requires a high degree of fine tuning to emerge, which is unlikely to occur in real biological systems. In the present study we propose a rather simple, and biologically relevant mechanism that profoundly expands the critical-like region. In particular, by means of numerical simulation we show that incorporating spatial heterogeneities into the square lattice of map-based excitable oscillators broadens the parameter space in which the distribution of excitation wave sizes follows closely a power law. Most importantly, this behavior is only observed if the spatial profile exhibits intermediate-sized patches with similar excitability levels, whereas for large and small spatial clusters only marginal widening of the critical state is detected. Furthermore, it turned out that the presence of spatial disorder in general amplifies the size of excitation waves, whereby the relatively highest contributions are observed in the proximity of the critical point. We argue that the reported mechanism is of particular importance for excitable systems with local interactions between individual elements.

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.

On the accuracy and computational cost of spiking neuron implementation

Since more than a decade ago, three statements about spiking neuron (SN) implementations have been widely accepted: 1) Hodgkin and Huxley (HH) model is computationally prohibitive, 2) Izhikevich (IZH) artificial neuron is as efficient as Leaky Integrate-and-Fire (LIF) model, and 3) IZH model is more efficient than HH model (Izhikevich, 2004). As suggested by Hodgkin and Huxley (1952), their model operates in two modes: by using the α's and β's rate functions directly (HH model) and by storing them into tables (HHT model) for computational cost reduction. Recently, it has been stated that: 1) HHT model (HH using tables) is not prohibitive, 2) IZH model is not efficient, and 3) both HHT and IZH models are comparable in computational cost (Skocik & Long, 2014). That controversy shows that there is no consensus concerning SN simulation capacities. Hence, in this work, we introduce a refined approach, based on the multiobjective optimization theory, describing the SN simulation capacities and ultimately choosing optimal simulation parameters. We have used normalized metrics to define the capacity levels of accuracy, computational cost, and efficiency. Normalized metrics allowed comparisons between SNs at the same level or scale. We conducted tests for balanced, lower, and upper boundary conditions under a regular spiking mode with constant and random current stimuli. We found optimal simulation parameters leading to a balance between computational cost and accuracy. Importantly, and, in general, we found that 1) HH model (without using tables) is the most accurate, computationally inexpensive, and efficient, 2) IZH model is the most expensive and inefficient, 3) both LIF and HHT models are the most inaccurate, 4) HHT model is more expensive and inaccurate than HH model due to α's and β's table discretization, and 5) HHT model is not comparable in computational cost to IZH model. These results refute the theory formulated over a decade ago (Izhikevich, 2004) and go more in-depth in the statements formulated by Skocik and Long (2014). Our statements imply that the number of dimensions or FLOPS in the SNs are theoretical but not practical indicators of the true computational cost. The metric we propose for the computational cost is more precise than FLOPS and was found to be invariant to computer architecture. Moreover, we found that the firing frequency used in previous works is a necessary but an insufficient metric to evaluate the simulation accuracy. We also show that our results are consistent with the theory of numerical methods and the theory of SN discontinuity. Discontinuous SNs, such LIF and IZH models, introduce a considerable error every time a spike is generated. In addition, compared to the constant input current, the random input current increases the computational cost and inaccuracy. Besides, we found that the search for optimal simulation parameters is problem-specific. That is important because most of the previous works have intended to find a general and unique optimal simulation. Here, we show that this solution could not exist because it is a multiobjective optimization problem that depends on several factors. This work sets up a renewed thesis concerning the SN simulation that is useful to several related research areas, including the emergent Deep Spiking Neural Networks.

Random temporal connections promote network synchronization

We report a phenomenon of collective dynamics on discrete-time complex networks: a random temporal interaction matrix even of zero or/and small average is able to significantly enhance synchronization with probability one. According to current knowledge, there is no verifiably sufficient criterion for the phenomenon. We use the standard method of synchronization analytics and the theory of stochastic processes to establish a criterion, by which we rigorously and accurately depict how synchronization occurring with probability one is affected by the statistical characteristics of the random temporal connections such as the strength and topology of the connections as well as their probability distributions. We also illustrate the enhancement phenomenon using physical and biological complex dynamical networks.

The emergence of synchrony in networks of mutually inferring neurons

This paper considers the emergence of a generalised synchrony in ensembles of coupled self-organising systems, such as neurons. We start from the premise that any self-organising system complies with the free energy principle, in virtue of placing an upper bound on its entropy. Crucially, the free energy principle allows one to interpret biological systems as inferring the state of their environment or external milieu. An emergent property of this inference is synchronisation among an ensemble of systems that infer each other. Here, we investigate the implications of neuronal dynamics by simulating neuronal networks, where each neuron minimises its free energy. We cast the ensuing ensemble dynamics in terms of inference and show that cardinal behaviours of neuronal networks - both in vivo and in vitro - can be explained by this framework. In particular, we test the hypotheses that (i) generalised synchrony is an emergent property of free energy minimisation; thereby explaining synchronisation in the resting brain: (ii) desynchronisation is induced by exogenous input; thereby explaining event-related desynchronisation and (iii) structure learning emerges in response to causal structure in exogenous input; thereby explaining functional segregation in real neuronal systems.

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.

Phase synchronization and intermittent behavior in healthy and Alzheimer-affected human-brain-based neural network

We study the dynamical proprieties of phase synchronization and intermittent behavior of neural systems using a network of networks structure based on an experimentally obtained human connectome for healthy and Alzheimer-affected brains. We consider a network composed of 78 neural subareas (subnetworks) coupled with a mean-field potential scheme. Each subnetwork is characterized by a small-world topology, composed of 250 bursting neurons simulated through a Rulkov model. Using the Kuramoto order parameter we demonstrate that healthy and Alzheimer-affected brains display distinct phase synchronization and intermittence properties as a function of internal and external coupling strengths. In general, for the healthy case, each subnetwork develops a substantial level of internal synchronization before a global stable phase-synchronization state has been established. For the unhealthy case, despite the similar internal subnetwork synchronization levels, we identify higher levels of global phase synchronization occurring even for relatively small internal and external coupling. Using recurrence quantification analysis, namely the determinism of the mean-field potential, we identify regions where the healthy and unhealthy networks depict nonstationary behavior, but the results denounce the presence of a larger region or intermittent dynamics for the case of Alzheimer-affected networks. A possible theoretical explanation based on two locally stable but globally unstable states is discussed.

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.

How synaptic plasticity influences spike synchronization and its transitions in complex neuronal network

There is evidence that synaptic plasticity is a vital feature of realistic neuronal systems. This study, describing synaptic plasticity by a modified Oja learning rule, focuses on the effect of synapse learning rate on spike synchronization and its relative transitions in a Newman-Watts small-world neuronal network. The individual dynamics of each neuron is modeled by a simple Rulkov map that produces spiking behavior. Numerical results have indicated that large coupling can lead to a spatiotemporally synchronous pattern of spiking neurons; in addition, this kind of spike synchronization can emerge intermittently by turning information transmission delay between coupled neurons. Interestingly, with the advent of synaptic plasticity, spike synchronization is gradually destroyed by increasing synapse learning rate; moreover, the phenomenon of intermittent synchronization transitions becomes less and less obvious and it even disappears for relative larger learning rate. Further simulations confirm that spike synchronization as well as synchronization transitions is largely independent of network size. Meanwhile, we detect that large shortcuts probability can facilitate spike synchronization, but it is disadvantageous for delay-induced synchronization transitions.

Chimera states in two-dimensional networks of locally coupled oscillators

Chimera state is defined as a mixed type of collective state in which synchronized and desynchronized subpopulations of a network of coupled oscillators coexist and the appearance of such anomalous behavior has strong connection to diverse neuronal developments. Most of the previous studies on chimera states are not extensively done in two-dimensional ensembles of coupled oscillators by taking neuronal systems with nonlinear coupling function into account while such ensembles of oscillators are more realistic from a neurobiological point of view. In this paper, we report the emergence and existence of chimera states by considering locally coupled two-dimensional networks of identical oscillators where each node is interacting through nonlinear coupling function. This is in contrast with the existence of chimera states in two-dimensional nonlocally coupled oscillators with rectangular kernel in the coupling function. We find that the presence of nonlinearity in the coupling function plays a key role to produce chimera states in two-dimensional locally coupled oscillators. We analytically verify explicitly in the case of a network of coupled Stuart-Landau oscillators in two dimensions that the obtained results using Ott-Antonsen approach and our analytical finding very well matches with the numerical results. Next, we consider another type of important nonlinear coupling function which exists in neuronal systems, namely chemical synaptic function, through which the nearest-neighbor (locally coupled) neurons interact with each other. It is shown that such synaptic interacting function promotes the emergence of chimera states in two-dimensional lattices of locally coupled neuronal oscillators. In numerical simulations, we consider two paradigmatic neuronal oscillators, namely Hindmarsh-Rose neuron model and Rulkov map for each node which exhibit bursting dynamics. By associating various spatiotemporal behaviors and snapshots at particular times, we study the chimera states in detail over a large range of coupling parameter. The existence of chimera states is confirmed by instantaneous angular frequency, order parameter and strength of incoherence.

Partial coupling delay induced multiple spatiotemporal orders in a modular neuronal network

The influence of partial coupling delay on the spatiotemporal spiking dynamics is explored in a modular neuronal network. The modular neuronal network is composed of two subnetworks which present the small-world property and scale-free property, respectively. Numerical results show that spatiotemporal order that the modular network is most coherent in time and nearly synchronized in space can emerge intermittently when the coupling delays among neurons are appropriately tuned. The appropriately tuned delays are further detected to be integer multiples of the intrinsic spiking period of the modular neuronal network, which implies that the phenomenon of multiple spatiotemporal orders could be the result of a locking between the length of coupling delay and the intrinsic spiking period of the modular neuronal network. Moreover, the multiple spatiotemporal orders are verified to be robust against variations of the fraction of delayed connection as well as the key parameters of network architecture such as the rewiring probability, the average degree of small-world subnetwork, the initial nodes of scale-free subnetwork and the total size of the modular network.

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

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 synchronisation by applying an external perturbation to each neuron. The simulations show that bursting synchronisation is more robust than spike synchronisation.

Phase diagrams and dynamics of a computationally efficient map-based neuron model

We introduce a new map-based neuron model derived from the dynamical perceptron family that has the best compromise between computational efficiency, analytical tractability, reduced parameter space and many dynamical behaviors. We calculate bifurcation and phase diagrams analytically and computationally that underpins a rich repertoire of autonomous and excitable dynamical behaviors. We report the existence of a new regime of cardiac spikes corresponding to nonchaotic aperiodic behavior. We compare the features of our model to standard neuron models currently available in the literature.

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.

Suppression of phase synchronisation in network based on cat's brain

We have studied the effects of perturbations on the cat's cerebral cortex. According to the literature, this cortex structure can be described by a clustered network. This way, we construct a clustered network with the same number of areas as in the cat matrix, where each area is described as a sub-network with a small-world property. We focus on the suppression of neuronal phase synchronisation considering different kinds of perturbations. Among the various controlling interventions, we choose three methods: delayed feedback control, external time-periodic driving, and activation of selected neurons. We simulate these interventions to provide a procedure to suppress undesired and pathological abnormal rhythms that can be associated with many forms of synchronisation. In our simulations, we have verified that the efficiency of synchronisation suppression by delayed feedback control is higher than external time-periodic driving and activation of selected neurons of the cat's cerebral cortex with the same coupling strengths.

Inverse period-doubling bifurcations determine complex structure of bursting in a one-dimensional non-autonomous map

We propose a simple one-dimensional non-autonomous map, in which some novel bursting patterns (e.g., "fold/double inverse flip" bursting, "fold/multiple inverse flip" bursting, and "fold/a cascade of inverse flip" bursting) can be observed. Typically, these bursting patterns exhibit complex structures containing a chain of inverse period-doubling bifurcations. The active states related to these bursting can be period-2(n) (n = 1, 2, 3,…) attractors or chaotic attractors, which may evolve to quiescence by a chain of inverse period-doubling bifurcations when the slow excitation decreases through period-doubling bifurcation points of the map. This accounts for the complex inverse period-doubling bifurcation structures observed in bursting patterns. Our findings enrich the possible routes to bursting as well as the underlying mechanisms of bursting.

Chimera states in bursting neurons

We study the existence of chimera states in pulse-coupled networks of bursting Hindmarsh-Rose neurons with nonlocal, global, and local (nearest neighbor) couplings. Through a linear stability analysis, we discuss the behavior of the stability function in the incoherent (i.e., disorder), coherent, chimera, and multichimera states. Surprisingly, we find that chimera and multichimera states occur even using local nearest neighbor interaction in a network of identical bursting neurons alone. This is in contrast with the existence of chimera states in populations of nonlocally or globally coupled oscillators. A chemical synaptic coupling function is used which plays a key role in the emergence of chimera states in bursting neurons. The existence of chimera, multichimera, coherent, and disordered states is confirmed by means of the recently introduced statistical measures and mean phase velocity.

Basin stability for burst synchronization in small-world networks of chaotic slow-fast oscillators

The impact of connectivity and individual dynamics on the basin stability of the burst synchronization regime in small-world networks consisting of chaotic slow-fast oscillators is studied. It is shown that there are rewiring probabilities corresponding to the largest basin stabilities, which uncovers a reason for finding small-world topologies in real neuronal networks. The impact of coupling density and strength as well as the nodal parameters of relaxation or excitability are studied. Dynamic mechanisms are uncovered that most strongly influence basin stability of the burst synchronization regime.

Local and global synchronization transitions induced by time delays in small-world neuronal networks with chemical synapses

Effects of time delay on the local and global synchronization in small-world neuronal networks with chemical synapses are investigated in this paper. Numerical results show that, for both excitatory and inhibitory coupling types, the information transmission delay can always induce synchronization transitions of spiking neurons in small-world networks. In particular, regions of in-phase and out-of-phase synchronization of connected neurons emerge intermittently as the synaptic delay increases. For excitatory coupling, all transitions to spiking synchronization occur approximately at integer multiples of the firing period of individual neurons; while for inhibitory coupling, these transitions appear at the odd multiples of the half of the firing period of neurons. More importantly, the local synchronization transition is more profound than the global synchronization transition, depending on the type of coupling synapse. For excitatory synapses, the local in-phase synchronization observed for some values of the delay also occur at a global scale; while for inhibitory ones, this synchronization, observed at the local scale, disappears at a global scale. Furthermore, the small-world structure can also affect the phase synchronization of neuronal networks. It is demonstrated that increasing the rewiring probability can always improve the global synchronization of neuronal activity, but has little effect on the local synchronization of neighboring neurons.

Two distinct desynchronization processes caused by lesions in globally coupled neurons

To accomplish a task, the brain works like a synchronized neuronal network where all the involved neurons work together. When a lesion spreads in the brain, depending on its evolution, it can reach a significant portion of relevant area. As a consequence, a phase transition might occur: the neurons desynchronize and cannot perform a certain task anymore. Lesions are responsible for either disrupting the neuronal connections or, in some cases, for killing the neuron. In this work, we will use a simplified model of neuronal network to show that these two types of lesions cause different types of desynchronization.Received: 20 November 2014,  Accepted: 10 March 2015; Edited by: C. A. Condat, G. J. Sibona; DOI: http://dx.doi.org/10.4279/PIP.070002Cite as: F A S Ferrari, R L Viana, Papers in Physics 7, 070002 (2015)  This paper, by Fabiano A. S. Ferrari, Ricardo L. Viana , is licensed under the Creative Commons Attribution License 3.0.

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.

Control for a synchronization-desynchronization switch

How to freely enhance or suppress synchronization of networked dynamical systems is of great importance in many disciplines. A unified precise control method for a synchronization-desynchronization switch, called the pull-push control method, is suggested. Namely, synchronization can be achieved when the original systems are desynchronous by pulling (or protecting) one node or a certain subset of nodes, whereas desynchronization can be accomplished when the systems are already synchronous by pushing (or kicking) one node or a certain subset of nodes. With this method, the controlled nodes should be chosen by the generalized eigenvector centrality of the critical synchronization mode of the Laplacian matrix. Compared with existing control methods for synchronization, it displays high efficiency, flexibility, and precision as well.

Modular networks with delayed coupling: synchronization and frequency control

We study the collective dynamics of modular networks consisting of map-based neurons which generate irregular spike sequences. Three types of intramodule topology are considered: a random Erdös-Rényi network, a small-world Watts-Strogatz network, and a scale-free Barabási-Albert network. The interaction between the neurons of different modules is organized by relatively sparse connections with time delay. For all the types of the network topology considered, we found that with increasing delay two regimes of module synchronization alternate with each other: inphase and antiphase. At the same time, the average rate of collective oscillations decreases within each of the time-delay intervals corresponding to a particular synchronization regime. A dual role of the time delay is thus established: controlling a synchronization mode and degree and controlling an average network frequency. Furthermore, we investigate the influence on the modular synchronization by other parameters: the strength of intermodule coupling and the individual firing rate.

Synchronization and stochastic resonance of the small-world neural network based on the CPG

According to biological knowledge, the central nervous system controls the central pattern generator (CPG) to drive the locomotion. The brain is a complex system consisting of different functions and different interconnections. The topological properties of the brain display features of small-world network. The synchronization and stochastic resonance have important roles in neural information transmission and processing. In order to study the synchronization and stochastic resonance of the brain based on the CPG, we establish the model which shows the relationship between the small-world neural network (SWNN) and the CPG. We analyze the synchronization of the SWNN when the amplitude and frequency of the CPG are changed and the effects on the CPG when the SWNN's parameters are changed. And we also study the stochastic resonance on the SWNN. The main findings include: (1) When the CPG is added into the SWNN, there exists parameters space of the CPG and the SWNN, which can make the synchronization of the SWNN optimum. (2) There exists an optimal noise level at which the resonance factor Q gets its peak value. And the correlation between the pacemaker frequency and the dynamical response of the network is resonantly dependent on the noise intensity. The results could have important implications for biological processes which are about interaction between the neural network and the CPG.

Noise induced complexity: patterns and collective phenomena in a small-world neuronal network

The effects of noise on patterns and collective phenomena are studied in a small-world neuronal network with the dynamics of each neuron being described by a two-dimensional Rulkov map neuron. It is shown that for intermediate noise levels, noise-induced ordered patterns emerge spatially, which supports the spatiotemporal coherence resonance. However, the inherent long range couplings of small-world networks can effectively disrupt the internal spatial scale of the media at small fraction of long-range couplings. The temporal order, characterized by the autocorrelation of a firing rate function, can be greatly enhanced by the introduction of small-world connectivity. There exists an optimal fraction of randomly rewired links, where the temporal order and synchronization can be optimized.