Phase synchronization of coupled bursting neurons and the generalized Kuramoto model

Kuramoto variables as eigenvalues of unitary matrices

We generalize the Kuramoto model by interpreting the N variables on the unit circle as eigenvalues of a N-dimensional unitary matrix U in three versions: general unitary, symmetric unitary, and special orthogonal. The time evolution is generated by N^{2} coupled differential equations for the matrix elements of U, and synchronization happens when U evolves into a multiple of the identity. The Ott-Antonsen ansatz is related to the Poisson kernels that are so useful in quantum transport, and we prove it in the case of identical natural frequencies. When the coupling constant is a matrix, we find some surprising new dynamical behaviors.

Generalized frustration in the multidimensional Kuramoto model

The Kuramoto model describes how coupled oscillators synchronize their phases as the intensity of the coupling increases past a threshold. The model was recently extended by reinterpreting the oscillators as particles moving on the surface of unit spheres in a D-dimensional space. Each particle is then represented by a D-dimensional unit vector; for D=2 the particles move on the unit circle and the vectors can be described by a single phase, recovering the original Kuramoto model. This multidimensional description can be further extended by promoting the coupling constant between the particles to a matrix K that acts on the unit vectors. As the coupling matrix changes the direction of the vectors, it can be interpreted as a generalized frustration that tends to hinder synchronization. In a recent paper we studied in detail the role of the coupling matrix for D=2. Here we extend this analysis to arbitrary dimensions. We show that, for identical particles, when the natural frequencies are set to zero, the system converges either to a stationary synchronized state, given by one of the real eigenvectors of K, or to an effective two-dimensional rotation, defined by one of the complex eigenvectors of K. The stability of these states depends on the set eigenvalues and eigenvectors of the coupling matrix, which controls the asymptotic behavior of the system, and therefore, can be used to manipulate these states. When the natural frequencies are not zero, synchronization depends on whether D is even or odd. In even dimensions the transition to synchronization is continuous and rotating states are replaced by active states, where the module of the order parameter oscillates while it rotates. If D is odd the phase transition is discontinuous and active states can be suppressed for some distributions of natural frequencies.

The brainstem connectome database

Connectivity data of the nervous system and subdivisions, such as the brainstem, cerebral cortex and subcortical nuclei, are necessary to understand connectional structures, predict effects of connectional disorders and simulate network dynamics. For that purpose, a database was built and analyzed which comprises all known directed and weighted connections within the rat brainstem. A longterm metastudy of original research publications describing tract tracing results form the foundation of the brainstem connectome (BC) database which can be analyzed directly in the framework neuroVIISAS. The BC database can be accessed directly by connectivity tables, a web-based tool and the framework. Analysis of global and local network properties, a motif analysis, and a community analysis of the brainstem connectome provides insight into its network organization. For example, we found that BC is a scale-free network with a small-world connectivity. The Louvain modularity and weighted stochastic block matching resulted in partially matching of functions and connectivity. BC modeling was performed to demonstrate signal propagation through the somatosensory pathway which is affected in Multiple sclerosis.

Measurement(s)	brainstem
Technology Type(s)	tract tracing metastudy
Factor Type(s)	brain region
Sample Characteristic - Organism	Rattus rattus
Sample Characteristic - Environment	Experimental setup
Sample Characteristic - Location	Germany

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.

A Model for Evolutionary Structural Plasticity and Synchronization of a Network of Neurons

A model of time-dependent structural plasticity for the synchronization of neuron networks is presented. It is known that synchronized oscillations reproduce structured communities, and this synchronization is transient since it can be enhanced or suppressed, and the proposed model reproduces this characteristic. The evolutionary behavior of the couplings is comparable to those of a network of biological neurons. In the structural network, the physical connections of axons and dendrites between neurons are modeled, and the evolution in the connections depends on the neurons' potential. Moreover, it is shown that the coupling force's function behaves as an adaptive controller that leads the neurons in the network to synchronization. The change in the node's degree shows that the network exhibits time-dependent structural plasticity, achieved through the evolutionary or adaptive change of the coupling force between the nodes. The coupling force function is based on the computed magnitude of the membrane potential deviations with its neighbors and a threshold that determines the neuron's connections. These rule the functional network structure along the time.

A model integrating multiple processes of synchronization and coherence for information instantiation within a cortical area

What is the form of dynamic, e.g., sensory, information in the mammalian cortex? Information in the cortex is modeled as a coherence map of a mixed chimera state of synchronous, phasic, and disordered minicolumns. The theoretical model is built on neurophysiological evidence. Complex spatiotemporal information is instantiated through a system of interacting biological processes that generate a synchronized cortical area, a coherent aperture. Minicolumn elements are grouped in macrocolumns in an array analogous to a phased-array radar, modeled as an aperture, a "hole through which radiant energy flows." Coherence maps in a cortical area transform inputs from multiple sources into outputs to multiple targets, while reducing complexity and entropy. Coherent apertures can assume extremely large numbers of different information states as coherence maps, which can be communicated among apertures with corresponding very large bandwidths. The coherent aperture model incorporates considerable reported research, integrating five conceptually and mathematically independent processes: 1) a damped Kuramoto network model, 2) a pumped area field potential, 3) the gating of nearly coincident spikes, 4) the coherence of activity across cortical lamina, and 5) complex information formed through functions in macrocolumns. Biological processes and their interactions are described in equations and a functional circuit such that the mathematical pieces can be assembled the same way the neurophysiological ones are. The model can be conceptually convolved over the specifics of local cortical areas within and across species. A coherent aperture becomes a node in a graph of cortical areas with a corresponding distribution of information.

Noninvasive Continuous Monitoring of Vital Signs With Wearables: Fit for Medical Use?

BACKGROUND

Wearables (= wearable computer) enable continuous and noninvasive monitoring of a range of vital signs. Mobile and cost-effective devices, combined with powerful data analysis tools, open new dimensions in assessing body functions ("digital biomarkers").

METHODS

To answer the question whether wearables are ready for use in the medical context, a PubMed literature search and analysis for their clinical-scientific use using publications from the years 2008 to 2018 was performed.

RESULTS

A total of 79 out of 314 search hits were publications on clinical trials with wearables, of which 16 were randomized controlled trials. Motion sensors were most frequently used to measure defined movements, movement disorders, or general physical activity. Approximately 20% of the studies used sensors to detect cardiovascular parameters. As for the sensor location, the wrist was chosen in most studies (22.8%).

CONCLUSION

Wearables can be used in a precisely defined medical context, when taking into account complex influencing factors.

Phase and amplitude dynamics of coupled oscillator systems on complex networks

We investigated locking behaviors of coupled limit-cycle oscillators with phase and amplitude dynamics. We focused on how the dynamics are affected by inhomogeneous coupling strength and by angular and radial shifts in coupling functions. We performed mean-field analyses of oscillator systems with inhomogeneous coupling strength, testing Gaussian, power-law, and brain-like degree distributions. Even for oscillators with identical intrinsic frequencies and intrinsic amplitudes, we found that the coupling strength distribution and the coupling function generated a wide repertoire of phase and amplitude dynamics. These included fully and partially locked states in which high-degree or low-degree nodes would phase-lead the network. The mean-field analytical findings were confirmed via numerical simulations. The results suggest that, in oscillator systems in which individual nodes can independently vary their amplitude over time, qualitatively different dynamics can be produced via shifts in the coupling strength distribution and the coupling form. Of particular relevance to information flows in oscillator networks, changes in the non-specific drive to individual nodes can make high-degree nodes phase-lag or phase-lead the rest of the network.

The role of timescale separation in oscillatory ensembles with competitive coupling

We study a heterogeneous population consisting of two groups of oscillatory elements, one with attractive and one with repulsive coupling. Moreover, we set different internal timescales for the oscillators of the two groups and concentrate on the role of this timescale separation in the collective behavior. Our results demonstrate that it may significantly modify synchronization properties of the system, and the implications are fundamentally different depending on the ratio between the group timescales. For the slower attractive group, synchronization properties are similar to the case of equal timescales. However, when the attractive group is faster, these properties significantly change and bistability appears. The other collective regimes such as frozen states and solitary states are also shown to be crucially influenced by timescale separation.

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.

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.

Switchable slow cellular conductances determine robustness and tunability of network states

Neuronal information processing is regulated by fast and localized fluctuations of brain states. Brain states reliably switch between distinct spatiotemporal signatures at a network scale even though they are composed of heterogeneous and variable rhythms at a cellular scale. We investigated the mechanisms of this network control in a conductance-based population model that reliably switches between active and oscillatory mean-fields. Robust control of the mean-field properties relies critically on a switchable negative intrinsic conductance at the cellular level. This conductance endows circuits with a shared cellular positive feedback that can switch population rhythms on and off at a cellular resolution. The switch is largely independent from other intrinsic neuronal properties, network size and synaptic connectivity. It is therefore compatible with the temporal variability and spatial heterogeneity induced by slower regulatory functions such as neuromodulation, synaptic plasticity and homeostasis. Strikingly, the required cellular mechanism is available in all cell types that possess T-type calcium channels but unavailable in computational models that neglect the slow kinetics of their activation.

Brain information processing involves electrophysiological signals at multiple temporal and spatial timescales, from the single neuron level to whole brain areas. A fast and local control of these signals by neurochemicals called neuromodulators is essential in complex tasks such as movement initiation and attentional focus. The neuromodulators act at the cellular scale to control signals that propagate at potentially much larger scales. The present paper highlights the critical role of a cellular switch of excitability for the fast and localized control of cellular and network states. By turning ON and OFF the cellular switch, neuromodulators can robustly switch large populations between distinct network states. We stress the importance of controlling the switch at a cellular level and independently of the connectivity to allow for tunable spatiotemporal signatures of the network states.

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.

Sampled-Data Consensus of Linear Multi-agent Systems With Packet Losses

In this paper, the consensus problem is studied for a class of multi-agent systems with sampled data and packet losses, where random and deterministic packet losses are considered, respectively. For random packet losses, a Bernoulli-distributed white sequence is used to describe packet dropouts among agents in a stochastic way. For deterministic packet losses, a switched system with stable and unstable subsystems is employed to model packet dropouts in a deterministic way. The purpose of this paper is to derive consensus criteria, such that linear multi-agent systems with sampled-data and packet losses can reach consensus. By means of the Lyapunov function approach and the decomposition method, the design problem of a distributed controller is solved in terms of convex optimization. The interplay among the allowable bound of the sampling interval, the probability of random packet losses, and the rate of deterministic packet losses are explicitly derived to characterize consensus conditions. The obtained criteria are closely related to the maximum eigenvalue of the Laplacian matrix versus the second minimum eigenvalue of the Laplacian matrix, which reveals the intrinsic effect of communication topologies on consensus performance. Finally, simulations are given to show the effectiveness of the proposed results.In this paper, the consensus problem is studied for a class of multi-agent systems with sampled data and packet losses, where random and deterministic packet losses are considered, respectively. For random packet losses, a Bernoulli-distributed white sequence is used to describe packet dropouts among agents in a stochastic way. For deterministic packet losses, a switched system with stable and unstable subsystems is employed to model packet dropouts in a deterministic way. The purpose of this paper is to derive consensus criteria, such that linear multi-agent systems with sampled-data and packet losses can reach consensus. By means of the Lyapunov function approach and the decomposition method, the design problem of a distributed controller is solved in terms of convex optimization. The interplay among the allowable bound of the sampling interval, the probability of random packet losses, and the rate of deterministic packet losses are explicitly derived to characterize consensus conditions. The obtained criteria are closely related to the maximum eigenvalue of the Laplacian matrix versus the second minimum eigenvalue of the Laplacian matrix, which reveals the intrinsic effect of communication topologies on consensus performance. Finally, simulations are given to show the effectiveness of the proposed results.

Chaotic Resonance in Typical Routes to Chaos in the Izhikevich Neuron Model

Chaotic resonance (CR), in which a system responds to a weak signal through the effects of chaotic activities, is a known function of chaos in neural systems. The current belief suggests that chaotic states are induced by different routes to chaos in spiking neural systems. However, few studies have compared the efficiency of signal responses in CR across the different chaotic states in spiking neural systems. We focused herein on the Izhikevich neuron model, comparing the characteristics of CR in the chaotic states arising through the period-doubling or tangent bifurcation routes. We found that the signal response in CR had a unimodal maximum with respect to the stability of chaotic orbits in the tested chaotic states. Furthermore, the efficiency of signal responses at the edge of chaos became especially high as a result of synchronization between the input signal and the periodic component in chaotic spiking activity.

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.

On synchronization in power-grids modelled as networks of second-order Kuramoto oscillators

This work concerns analytical results on the role of coupling strength in the phenomenon of onset of complete frequency locking in power-grids modelled as a network of second-order Kuramoto oscillators. Those results allow estimation of the coupling strength for the onset of complete frequency locking and to assess the features of network and oscillators that favor synchronization. The analytical results are evaluated using an order parameter defined as the normalized sum of absolute values of phase deviations of the oscillators over time. The investigation of the frequency synchronization within the subsets of the parameter space involved in the synchronization problem is also carried out. It is shown that the analytical results are in good agreement with those observed in the numerical simulations. In order to illustrate the methodology, a case study is presented, involving the Brazilian high-voltage transmission system under a load peak condition to study the effect of load on the syncronizability of the grid. The results show that both the load and the centralized generation might have concurred to the 2014 blackout.