Hierarchical synchrony of phase oscillators in modular networks

Metastability of multi-population Kuramoto-Sakaguchi oscillators

An Ott-Antonsen reduced M-population of Kuramoto-Sakaguchi oscillators is investigated, focusing on the influence of the phase-lag parameter α on the collective dynamics. For oscillator populations coupled on a ring, we obtained a wide variety of spatiotemporal patterns, including coherent states, traveling waves, partially synchronized states, modulated states, and incoherent states. Back-and-forth transitions between these states are found, which suggest metastability. Linear stability analysis reveals the stable regions of coherent states with different winding numbers q. Within certain α ranges, the system settles into stable traveling wave solutions despite the coherent states also being linearly stable. For around α≈0.46π, the system displays the most frequent metastable transitions between coherent states and partially synchronized states, while for α closer to π/2, metastable transitions arise between partially synchronized states and modulated states. This model captures metastable dynamics akin to brain activity, offering insights into the synchronization of brain networks.

Modular subgraphs in large-scale connectomes underpin spontaneous co-fluctuation events in mouse and human brains

Previous studies have adopted an edge-centric framework to study fine-scale network dynamics in human fMRI. To date, however, no studies have applied this framework to data collected from model organisms. Here, we analyze structural and functional imaging data from lightly anesthetized mice through an edge-centric lens. We find evidence of "bursty" dynamics and events - brief periods of high-amplitude network connectivity. Further, we show that on a per-frame basis events best explain static FC and can be divided into a series of hierarchically-related clusters. The co-fluctuation patterns associated with each cluster centroid link distinct anatomical areas and largely adhere to the boundaries of algorithmically detected functional brain systems. We then investigate the anatomical connectivity undergirding high-amplitude co-fluctuation patterns. We find that events induce modular bipartitions of the anatomical network of inter-areal axonal projections. Finally, we replicate these same findings in a human imaging dataset. In summary, this report recapitulates in a model organism many of the same phenomena observed in previously edge-centric analyses of human imaging data. However, unlike human subjects, the murine nervous system is amenable to invasive experimental perturbations. Thus, this study sets the stage for future investigation into the causal origins of fine-scale brain dynamics and high-amplitude co-fluctuations. Moreover, the cross-species consistency of the reported findings enhances the likelihood of future translation.

Exponentially Long Transient Time to Synchronization of Coupled Chaotic Circle Maps in Dense Random Networks

We study the transition to synchronization in large, dense networks of chaotic circle maps, where an exact solution of the mean-field dynamics in the infinite network and all-to-all coupling limit is known. In dense networks of finite size and link probability of smaller than one, the incoherent state is meta-stable for coupling strengths that are larger than the mean-field critical coupling. We observe chaotic transients with exponentially distributed escape times and study the scaling behavior of the mean time to synchronization.

Multistability in coupled oscillator systems with higher-order interactions and community structure

We study synchronization dynamics in populations of coupled phase oscillators with higher-order interactions and community structure. We find that the combination of these two properties gives rise to a number of states unsupported by either higher-order interactions or community structure alone, including synchronized states with communities organized into clusters in-phase, anti-phase, and a novel skew-phase, as well as an incoherent-synchronized state. Moreover, the system displays strong multistability with many of these states stable at the same time. We demonstrate our findings by deriving the low dimensional dynamics of the system and examining the system's bifurcations using stability analysis and perturbation theory.

Chimeras on a ring of oscillator populations

Chimeras occur in networks of coupled oscillators and are characterized by coexisting groups of synchronous oscillators and asynchronous oscillators. We consider a network formed from N equal-sized populations at equally spaced points around a ring. We use the Ott/Antonsen ansatz to derive coupled ordinary differential equations governing the level of synchrony within each population and describe chimeras using a self-consistency argument. For N=2 and 3, our results are compared with previously known ones. We obtain new results for the cases of 4,5,…,12 populations and a numerically based conjecture resulting from the behavior of larger numbers of populations. We find macroscopic chaos when more than five populations are considered, but conjecture that this behavior vanishes as the number of populations is increased.

Modulations of local synchrony over time lead to resting-state functional connectivity in a parsimonious large-scale brain model

Biophysical models of large-scale brain activity are a fundamental tool for understanding the mechanisms underlying the patterns observed with neuroimaging. These models combine a macroscopic description of the within- and between-ensemble dynamics of neurons within a single architecture. A challenge for these models is accounting for modulations of within-ensemble synchrony over time. Such modulations in local synchrony are fundamental for modeling behavioral tasks and resting-state activity. Another challenge comes from the difficulty in parametrizing large scale brain models which hinders researching principles related with between-ensembles differences. Here we derive a parsimonious large scale brain model that can describe fluctuations of local synchrony. Crucially, we do not reduce within-ensemble dynamics to macroscopic variables first, instead we consider within and between-ensemble interactions similarly while preserving their physiological differences. The dynamics of within-ensemble synchrony can be tuned with a parameter which manipulates local connectivity strength. We simulated resting-state static and time-resolved functional connectivity of alpha band envelopes in models with identical and dissimilar local connectivities. We show that functional connectivity emerges when there are high fluctuations of local and global synchrony simultaneously (i.e. metastable dynamics). We also show that for most ensembles, leaning towards local asynchrony or synchrony correlates with the functional connectivity with other ensembles, with the exception of some regions belonging to the default-mode network.

Optimizing synchrony with a minimal coupling strength of coupled phase oscillators on complex networks based on desynchronous clustering

Finding the globally optimal network is an unsolved problem in synchrony optimization. In this paper, an efficient edge-adding optimization method based on global information is proposed. The edge-adding scheme is obtained through the eigenvector corresponding to the maximum eigenvalue of the system's Jacobian matrix. With different frequency distributions, we find that the optimized networks have similar features, concluded as three conditions: (i) The deviations of nodes' frequencies from the mean value are linear with the nodes' degrees, (ii) the oscillators form a bipartite network divided according to the frequencies of the oscillators, and (iii) oscillators are only connected to those with sufficiently large frequency differences. An optimal network can be constructed directly based on these three conditions for a given distribution of natural frequencies. We show that the critical coupling strengths of these constructed networks approach the theoretical lower bound. The constructed networks are or at least close to the globally optimal ones for synchrony.

Thalamic bursts modulate cortical synchrony locally to switch between states of global functional connectivity in a cognitive task

Performing a cognitive task requires going through a sequence of functionally diverse stages. Although it is typically assumed that these stages are characterized by distinct states of cortical synchrony that are triggered by sub-cortical events, little reported evidence supports this hypothesis. To test this hypothesis, we first identified cognitive stages in single-trial MEG data of an associative recognition task, showing with a novel method that each stage begins with local modulations of synchrony followed by a state of directed functional connectivity. Second, we developed the first whole-brain model that can simulate cortical synchrony throughout a task. The model suggests that the observed synchrony is caused by thalamocortical bursts at the onset of each stage, targeted at cortical synapses and interacting with the structural anatomical connectivity. These findings confirm that cognitive stages are defined by distinct states of cortical synchrony and explains the network-level mechanisms necessary for reaching stage-dependent synchrony states.

Path-dependent dynamics induced by rewiring networks of inertial oscillators

In networks of coupled oscillators, it is of interest to understand how interaction topology affects synchronization. Many studies have gained key insights into this question by studying the classic Kuramoto oscillator model on static networks. However, new questions arise when the network structure is time varying or when the oscillator system is multistable, the latter of which can occur when an inertial term is added to the Kuramoto model. While the consequences of evolving topology and multistability on collective behavior have been examined separately, real-world systems such as gene regulatory networks and the brain may exhibit these properties simultaneously. It is thus relevant to ask how time-varying network connectivity impacts synchronization in systems that can exhibit multistability. To address this question, we study how the dynamics of coupled Kuramoto oscillators with inertia are affected when the topology of the underlying network changes in time. We show that hysteretic synchronization behavior in networks of coupled inertial oscillators can be driven by changes in connection topology alone. Moreover, we find that certain fixed-density rewiring schemes induce significant changes to the level of global synchrony that remain even after the network returns to its initial configuration, and we show that these changes are robust to a wide range of network perturbations. Our findings highlight that the specific progression of network topology over time, in addition to its initial or final static structure, can play a considerable role in modulating the collective behavior of systems evolving on complex networks.

Chimera states in coupled map lattices: Spatiotemporally intermittent behavior and an equivalent cellular automaton

We construct an equivalent cellular automaton (CA) for a system of globally coupled sine circle maps with two populations and distinct values for intergroup and intragroup coupling. The phase diagram of the system shows that the coupled map lattice can exhibit chimera states with synchronized and spatiotemporally intermittent subgroups after evolution from random initial conditions in some parameter regimes, as well as to other kinds of solutions in other parameter regimes. The CA constructed by us reflects the global nature and the two population structure of the coupled map lattice and is able to reproduce the phase diagram accurately. The CA depends only on the total number of laminar and burst sites and shows a transition from co-existing deterministic and probabilistic behavior in the chimera region to fully probabilistic behavior at the phase boundaries. This identifies the characteristic signature of the transition of a cellular automaton to a chimera state. We also construct an evolution equation for the average number of laminar/burst sites from the CA, analyze its behavior and solutions, and correlate these with the behavior seen for the coupled map lattice. Our CA and methods of analysis can have relevance in wider contexts.

Anti-phase collective synchronization with intrinsic in-phase coupling of two groups of electrochemical oscillators

The synchronization of two groups of electrochemical oscillators is investigated during the electrodissolution of nickel in sulfuric acid. The oscillations are coupled through combined capacitance and resistance, so that in a single pair of oscillators (nearly) in-phase synchronization is obtained. The internal coupling within each group is relatively strong, but there is a phase difference between the fast and slow oscillators. The external coupling between the two groups is weak. The experiments show that the two groups can exhibit (nearly) anti-phase collective synchronization. Such synchronization occurs only when the external coupling is weak, and the interactions are delayed by the capacitance. When the external coupling is restricted to those between the fast and the slow elements, the anti-phase synchronization is more prominent. The results are interpreted with phase models. The theory predicts that, for anti-phase collective synchronization, there must be a minimum internal phase difference for a given shift in the phase coupling function. This condition is less stringent with external fast-to-slow coupling. The results provide a framework for applications of collective phase synchronization in modular networks where weak coupling between the groups can induce synchronization without rearrangements of the phase dynamics within the groups. This article is part of the theme issue 'Coupling functions: dynamical interaction mechanisms in the physical, biological and social sciences'.

Low-dimensional dynamics of the Kuramoto model with rational frequency distributions

The Kuramoto model is a paradigmatic tool for studying the dynamics of collective behavior in large ensembles of coupled dynamical systems. Over the past decade a great deal of progress has been made in analytical descriptions of the macroscopic dynamics of the Kuramoto model, facilitated by the discovery of Ott and Antonsen's dimensionality reduction method. However, the vast majority of these works relies on a critical assumption where the oscillators' natural frequencies are drawn from a Cauchy, or Lorentzian, distribution, which allows for a convenient closure of the evolution equations from the dimensionality reduction. In this paper we investigate the low-dimensional dynamics that emerge from a broader family of natural frequency distributions, in particular, a family of rational distribution functions. We show that, as the polynomials that characterize the frequency distribution increase in order, the low-dimensional evolution equations become more complicated, but nonetheless the system dynamics remain simple, displaying a transition from incoherence to partial synchronization at a critical coupling strength. Using the low-dimensional equations we analytically calculate the critical coupling strength corresponding to the onset of synchronization and investigate the scaling properties of the order parameter near the onset of synchronization. These results agree with calculations from Kuramoto's original self-consistency framework, but we emphasize that the low-dimensional equations approach used here allows for a true stability analysis categorizing the bifurcations.

Structural disconnection is responsible for increased functional connectivity in multiple sclerosis

Increased synchrony within neuroanatomical networks is often observed in neurophysiologic studies of human brain disease. Most often, this phenomenon is ascribed to a compensatory process in the face of injury, though evidence supporting such accounts is limited. Given the known dependence of resting-state functional connectivity (rsFC) on underlying structural connectivity (SC), we examine an alternative hypothesis: that topographical changes in SC, specifically particular patterns of disconnection, contribute to increased network rsFC. We obtain measures of rsFC using fMRI and SC using probabilistic tractography in 50 healthy and 28 multiple sclerosis subjects. Using a computational model of neuronal dynamics, we simulate BOLD using healthy subject SC to couple regions. We find that altering the model by introducing structural disconnection patterns observed in those multiple sclerosis subjects with high network rsFC generates simulations with high rsFC as well, suggesting that disconnection itself plays a role in producing high network functional connectivity. We then examine SC data in individuals. In multiple sclerosis subjects with high network rsFC, we find a preferential disconnection between the relevant network and wider system. We examine the significance of such network isolation by introducing random disconnection into the model. As observed empirically, simulated network rsFC increases with removal of connections bridging a community with the remainder of the brain. We thus show that structural disconnection known to occur in multiple sclerosis contributes to network rsFC changes in multiple sclerosis and further that community isolation is responsible for elevated network functional connectivity.

Diffusion dynamics and synchronizability of hierarchical products of networks

The hierarchical product of networks represents a natural tool for building large networks out of two smaller subnetworks: a primary subnetwork and a secondary subnetwork. Here we study the dynamics of diffusion and synchronization processes on hierarchical products. We apply techniques previously used for approximating the eigenvalues of the adjacency matrix to the Laplacian matrix, allowing us to quantify the effects that the primary and secondary subnetworks have on diffusion and synchronization in terms of a coupling parameter that weighs the secondary subnetwork relative to the primary subnetwork. Diffusion processes are separated into two regimes: for small coupling the diffusion rate is determined by the structure of the secondary network, scaling with the coupling parameter, while for large coupling it is determined by the primary network and saturates. Synchronization, on the other hand, is separated into three regimes, for both small and large coupling hierarchical products have poor synchronization properties, but is optimized at an intermediate value. Moreover, the critical coupling value that optimizes synchronization is shaped by the relative connectivities of the primary and secondary subnetworks, compensating for significant differences between the two subnetworks.

Uncovering low dimensional macroscopic chaotic dynamics of large finite size complex systems

In the last decade, it has been shown that a large class of phase oscillator models admit low dimensional descriptions for the macroscopic system dynamics in the limit of an infinite number N of oscillators. The question of whether the macroscopic dynamics of other similar systems also have a low dimensional description in the infinite N limit has, however, remained elusive. In this paper, we show how techniques originally designed to analyze noisy experimental chaotic time series can be used to identify effective low dimensional macroscopic descriptions from simulations with a finite number of elements. We illustrate and verify the effectiveness of our approach by applying it to the dynamics of an ensemble of globally coupled Landau-Stuart oscillators for which we demonstrate low dimensional macroscopic chaotic behavior with an effective 4-dimensional description. By using this description, we show that one can calculate dynamical invariants such as Lyapunov exponents and attractor dimensions. One could also use the reconstruction to generate short-term predictions of the macroscopic dynamics.

Development of structural correlations and synchronization from adaptive rewiring in networks of Kuramoto oscillators

Synchronization of non-identical oscillators coupled through complex networks is an important example of collective behavior, and it is interesting to ask how the structural organization of network interactions influences this process. Several studies have explored and uncovered optimal topologies for synchronization by making purposeful alterations to a network. On the other hand, the connectivity patterns of many natural systems are often not static, but are rather modulated over time according to their dynamics. However, this co-evolution and the extent to which the dynamics of the individual units can shape the organization of the network itself are less well understood. Here, we study initially randomly connected but locally adaptive networks of Kuramoto oscillators. In particular, the system employs a co-evolutionary rewiring strategy that depends only on the instantaneous, pairwise phase differences of neighboring oscillators, and that conserves the total number of edges, allowing the effects of local reorganization to be isolated. We find that a simple rule-which preserves connections between more out-of-phase oscillators while rewiring connections between more in-phase oscillators-can cause initially disordered networks to organize into more structured topologies that support enhanced synchronization dynamics. We examine how this process unfolds over time, finding a dependence on the intrinsic frequencies of the oscillators, the global coupling, and the network density, in terms of how the adaptive mechanism reorganizes the network and influences the dynamics. Importantly, for large enough coupling and after sufficient adaptation, the resulting networks exhibit interesting characteristics, including degree-frequency and frequency-neighbor frequency correlations. These properties have previously been associated with optimal synchronization or explosive transitions in which the networks were constructed using global information. On the contrary, by considering a time-dependent interplay between structure and dynamics, this work offers a mechanism through which emergent phenomena and organization can arise in complex systems utilizing local rules.

Collective phase reduction of globally coupled noisy dynamical elements

We formulate a theory for the collective phase reduction of globally coupled noisy dynamical elements exhibiting macroscopic rhythms. We first transform the Langevin-type equation that represents a group of globally coupled noisy dynamical elements into the corresponding nonlinear Fokker-Planck equation and then develop the phase reduction method for limit-cycle solutions to the nonlinear Fokker-Planck equation. The theory enables us to describe the collective dynamics of a group of globally coupled noisy dynamical elements by a single degree of freedom called the collective phase. As long as the group collectively exhibits macroscopic rhythms, the theory is applicable even when the coupling and noise are strong; it is also independent of the assumption that each element of the group is a self-sustained oscillator. We also provide a simple and accurate numerical algorithm for the collective phase description method and numerically illustrate the theory using a group of globally coupled noisy FitzHugh-Nagumo elements.

Mathematical Frameworks for Oscillatory Network Dynamics in Neuroscience

The tools of weakly coupled phase oscillator theory have had a profound impact on the neuroscience community, providing insight into a variety of network behaviours ranging from central pattern generation to synchronisation, as well as predicting novel network states such as chimeras. However, there are many instances where this theory is expected to break down, say in the presence of strong coupling, or must be carefully interpreted, as in the presence of stochastic forcing. There are also surprises in the dynamical complexity of the attractors that can robustly appear-for example, heteroclinic network attractors. In this review we present a set of mathematical tools that are suitable for addressing the dynamics of oscillatory neural networks, broadening from a standard phase oscillator perspective to provide a practical framework for further successful applications of mathematics to understanding network dynamics in neuroscience.

Spectral properties of the hierarchical product of graphs

The hierarchical product of two graphs represents a natural way to build a larger graph out of two smaller graphs with less regular and therefore more heterogeneous structure than the Cartesian product. Here we study the eigenvalue spectrum of the adjacency matrix of the hierarchical product of two graphs. Introducing a coupling parameter describing the relative contribution of each of the two smaller graphs, we perform an asymptotic analysis for the full spectrum of eigenvalues of the adjacency matrix of the hierarchical product. Specifically, we derive the exact limit points for each eigenvalue in the limits of small and large coupling, as well as the leading-order relaxation to these values in terms of the eigenvalues and eigenvectors of the two smaller graphs. Given its central roll in the structural and dynamical properties of networks, we study in detail the Perron-Frobenius, or largest, eigenvalue. Finally, as an example application we use our theory to predict the epidemic threshold of the susceptible-infected-susceptible model on a hierarchical product of two graphs.

Ott-Antonsen attractiveness for parameter-dependent oscillatory systems

The Ott-Antonsen (OA) ansatz [Ott and Antonsen, Chaos 18, 037113 (2008); Chaos 19, 023117 (2009)] has been widely used to describe large systems of coupled phase oscillators. If the coupling is sinusoidal and if the phase dynamics does not depend on the specific oscillator, then the macroscopic behavior of the systems can be fully described by a low-dimensional dynamics. Does the corresponding manifold remain attractive when introducing an intrinsic dependence between an oscillator's phase and its dynamics by additional, oscillator specific parameters? To answer this, we extended the OA ansatz and proved that parameter-dependent oscillatory systems converge to the OA manifold given certain conditions. Our proof confirms recent numerical findings that already hinted at this convergence. Furthermore, we offer a thorough mathematical underpinning for networks of so-called theta neurons, where the OA ansatz has just been applied. In a final step, we extend our proof by allowing for time-dependent and multi-dimensional parameters as well as for network topologies other than global coupling. This renders the OA ansatz an excellent starting point for the analysis of a broad class of realistic settings.

Optimal synchronization of directed complex networks

We study optimal synchronization of networks of coupled phase oscillators. We extend previous theory for optimizing the synchronization properties of undirected networks to the important case of directed networks. We derive a generalized synchrony alignment function that encodes the interplay between the network structure and the oscillators' natural frequencies and serves as an objective measure for the network's degree of synchronization. Using the generalized synchrony alignment function, we show that a network's synchronization properties can be systematically optimized. This framework also allows us to study the properties of synchrony-optimized networks, and in particular, investigate the role of directed network properties such as nodal in- and out-degrees. For instance, we find that in optimally rewired networks, the heterogeneity of the in-degree distribution roughly matches the heterogeneity of the natural frequency distribution, but no such relationship emerges for out-degrees. We also observe that a network's synchronization properties are promoted by a strong correlation between the nodal in-degrees and the natural frequencies of oscillators, whereas the relationship between the nodal out-degrees and the natural frequencies has comparatively little effect. This result is supported by our theory, which indicates that synchronization is promoted by a strong alignment of the natural frequencies with the left singular vectors corresponding to the largest singular values of the Laplacian matrix.

Dynamics of two populations of phase oscillators with different frequency distributions

A large variety of rhythms are observed in nature. Rhythms such as electroencephalogram signals in the brain can often be regarded as interacting. In this study, we investigate the dynamical properties of rhythmic systems in two populations of phase oscillators with different frequency distributions. We assume that the average frequency ratio between two populations closely approximates some small integer. Most importantly, we adopt a specific coupling function derived from phase reduction theory. Under some additional assumptions, the system of two populations of coupled phase oscillators reduces to a low-dimensional system in the continuum limit. Consequently, we find chimera states in which clustering and incoherent states coexist. Finally, we confirm consistent behaviors of the derived low-dimensional model and the original model.

Erosion of synchronization: Coupling heterogeneity and network structure

We study the dynamics of network-coupled phase oscillators in the presence of coupling frustration. It was recently demonstrated that in heterogeneous network topologies, the presence of coupling frustration causes perfect phase synchronization to become unattainable even in the limit of infinite coupling strength. Here, we consider the important case of heterogeneous coupling functions and extend previous results by deriving analytical predictions for the total erosion of synchronization. Our analytical results are given in terms of basic quantities related to the network structure and coupling frustration. In addition to fully heterogeneous coupling, where each individual interaction is allowed to be distinct, we also consider partially heterogeneous coupling and homogeneous coupling in which the coupling functions are either unique to each oscillator or identical for all network interactions, respectively. We demonstrate the validity of our theory with numerical simulations of multiple network models, and highlight the interesting effects that various coupling choices and network models have on the total erosion of synchronization. Finally, we consider some special network structures with well-known spectral properties, which allows us to derive further analytical results.

Hamiltonian mean field model: Effect of network structure on synchronization dynamics

The Hamiltonian mean field model of coupled inertial Hamiltonian rotors is a prototype for conservative dynamics in systems with long-range interactions. We consider the case where the interactions between the rotors are governed by a network described by a weighted adjacency matrix. By studying the linear stability of the incoherent state, we find that the transition to synchrony begins when the coupling constant K is inversely proportional to the largest eigenvalue of the adjacency matrix. We derive a closed system of equations for a set of local order parameters to study the effect of network heterogeneity on the synchronization of the rotors. When K is just beyond the transition to synchronization, we find that the degree of synchronization is highly dependent on the network's heterogeneity, but that for large K the degree of synchronization is robust to changes in the degree distribution. Our results are illustrated with numerical simulations on Erdös-Renyi networks and networks with power-law degree distributions.

Dynamics of globally coupled oscillators: Progress and perspectives

In this paper, we discuss recent progress in research of ensembles of mean field coupled oscillators. Without an ambition to present a comprehensive review, we outline most interesting from our viewpoint results and surprises, as well as interrelations between different approaches.

Control of coupled oscillator networks with application to microgrid technologies

The control of complex systems and network-coupled dynamical systems is a topic of vital theoretical importance in mathematics and physics with a wide range of applications in engineering and various other sciences. Motivated by recent research into smart grid technologies, we study the control of synchronization and consider the important case of networks of coupled phase oscillators with nonlinear interactions-a paradigmatic example that has guided our understanding of self-organization for decades. We develop a method for control based on identifying and stabilizing problematic oscillators, resulting in a stable spectrum of eigenvalues, and in turn a linearly stable synchronized state. The amount of control, that is, number of oscillators, required to stabilize the network is primarily dictated by the coupling strength, dynamical heterogeneity, and mean degree of the network, and depends little on the structural heterogeneity of the network itself.

Frequency assortativity can induce chaos in oscillator networks

We investigate the effect of preferentially connecting oscillators with similar frequency to each other in networks of coupled phase oscillators (i.e., frequency assortativity). Using the network Kuramoto model as an example, we find that frequency assortativity can induce chaos in the macroscopic dynamics. By applying a mean-field approximation in combination with the dimension reduction method of Ott and Antonsen, we show that the dynamics can be described by a low dimensional system of equations. We use the reduced system to characterize the macroscopic chaos using Lyapunov exponents, bifurcation diagrams, and time-delay embeddings. Finally, we show that the emergence of chaos stems from the formation of multiple groups of synchronized oscillators, i.e., meta-oscillators.

Chimera states in time-varying complex networks

Chimera states have been recently found in a variety of different coupling schemes and geometries. In most cases, the underlying coupling structure is considered to be static, while many realistic systems display significant temporal changes in the pattern of connectivity. In this work we investigate a time-varying network made of two coupled populations of Kuramoto oscillators, where the links between the two groups are considered to vary over time. As a main result we find that the network may support stable, breathing, and alternating chimera states. We also find that, when the rate of connectivity changes is fast, compared to the oscillator dynamics, the network may be described by a low-dimensional system of equations. Unlike in the static heterogeneous case, the onset of alternating chimera states is due to the presence of fluctuations, which may be induced either by the finite size of the network or by large switching times.

Erosion of synchronization in networks of coupled oscillators

We report erosion of synchronization in networks of coupled phase oscillators, a phenomenon where perfect phase synchronization is unattainable in steady state, even in the limit of infinite coupling. An analysis reveals that the total erosion is separable into the product of terms characterizing coupling frustration and structural heterogeneity, both of which amplify erosion. The latter, however, can differ significantly from degree heterogeneity. Finally, we show that erosion is marked by the reorganization of oscillators according to their node degrees rather than their natural frequencies.

Optimal synchronization of complex networks

We study optimal synchronization in networks of heterogeneous phase oscillators. Our main result is the derivation of a synchrony alignment function that encodes the interplay between network structure and oscillators' frequencies and that can be readily optimized. We highlight its utility in two general problems: constrained frequency allocation and network design. In general, we find that synchronization is promoted by strong alignments between frequencies and the dominant Laplacian eigenvectors, as well as a matching between the heterogeneity of frequencies and network structure.

Synchronization of bursting Hodgkin-Huxley-type neurons in clustered networks

We considered a clustered network of bursting neurons described by the Huber-Braun model. In the upper level of the network we used the connectivity matrix of the cat cerebral cortex network, and in the lower level each cortex area (or cluster) is modelled as a small-world network. There are two different coupling strengths related to inter- and intracluster dynamics. Each bursting cycle is composed of a quiescent period followed by a rapid chaotic sequence of spikes, and we defined a geometric phase which enables us to investigate the onset of synchronized bursting, as the state in which the neuron start bursting at the same time, whereas their spikes may remain uncorrelated. The bursting synchronization of a clustered network has been investigated using an order parameter and the average field of the network in order to identify regimes in which each cluster may display synchronized behavior, whereas the overall network does not. We introduce quantifiers to evaluate the relative contribution of each cluster in the partial synchronized behavior of the whole network. Our main finding is that we typically observe in the clustered network not a complete phase synchronized regime but instead a complex pattern of partial phase synchronization in which different cortical areas may be internally synchronized at distinct phase values, hence they are not externally synchronized, unless the coupling strengths are too large.

Frustrated hierarchical synchronization and emergent complexity in the human connectome network

The spontaneous emergence of coherent behavior through synchronization plays a key role in neural function, and its anomalies often lie at the basis of pathologies. Here we employ a parsimonious (mesoscopic) approach to study analytically and computationally the synchronization (Kuramoto) dynamics on the actual human-brain connectome network. We elucidate the existence of a so-far-uncovered intermediate phase, placed between the standard synchronous and asynchronous phases, i.e. between order and disorder. This novel phase stems from the hierarchical modular organization of the connectome. Where one would expect a hierarchical synchronization process, we show that the interplay between structural bottlenecks and quenched intrinsic frequency heterogeneities at many different scales, gives rise to frustrated synchronization, metastability, and chimera-like states, resulting in a very rich and complex phenomenology. We uncover the origin of the dynamic freezing behind these features by using spectral graph theory and discuss how the emerging complex synchronization patterns relate to the need for the brain to access –in a robust though flexible way– a large variety of functional attractors and dynamical repertoires without ad hoc fine-tuning to a critical point.

Glassy states and super-relaxation in populations of coupled phase oscillators

Large networks of coupled oscillators appear in many branches of science, so that the kinds of phenomena they exhibit are not only of intrinsic interest but also of very wide importance. In 1975, Kuramoto proposed an analytically tractable model to describe these systems, which has since been successfully applied in many contexts and remains a subject of intensive research. Some related problems, however, remain unclarified for decades, such as the existence and properties of the oscillator glass state. Here we present a detailed analysis of a very general form of the Kuramoto model. In particular, we find the conditions when it can exhibit glassy behaviour, which represents a kind of synchronous disorder in the present case. Furthermore, we discover a new and intriguing phenomenon that we refer to as super-relaxation where the oscillators feel no interaction at all while relaxing to incoherence. Our findings offer the possibility of creating glassy states and observing super-relaxation in real systems.

Synchronization of oscillators in a Kuramoto-type model with generic coupling

We study synchronization properties of coupled oscillators on networks that allow description in terms of global mean field coupling. These models generalize the standard Kuramoto-Sakaguchi model, allowing for different contributions of oscillators to the mean field and to different forces from the mean field on oscillators. We present the explicit solutions of self-consistency equations for the amplitude and frequency of the mean field in a parametric form, valid for noise-free and noise-driven oscillators. As an example, we consider spatially spreaded oscillators for which the coupling properties are determined by finite velocity of signal propagation.

Solitary state at the edge of synchrony in ensembles with attractive and repulsive interactions

We discuss the desynchronization transition in networks of globally coupled identical oscillators with attractive and repulsive interactions. We show that, if attractive and repulsive groups act in antiphase or close to that, a solitary state emerges with a single repulsive oscillator split up from the others fully synchronized. With further increase of the repulsing strength, the synchronized cluster becomes fuzzy and the dynamics is given by a variety of stationary states with zero common forcing. Intriguingly, solitary states represent the natural link between coherence and incoherence. The phenomenon is described analytically for phase oscillators with sine coupling and demonstrated numerically for more general amplitude models.

Phase synchronization between collective rhythms of fully locked oscillator groups

A system of coupled oscillators can exhibit a rich variety of dynamical behaviors. When we investigate the dynamical properties of the system, we first analyze individual oscillators and the microscopic interactions between them. However, the structure of a coupled oscillator system is often hierarchical, so that the collective behaviors of the system cannot be fully clarified by simply analyzing each element of the system. For example, we found that two weakly interacting groups of coupled oscillators can exhibit anti-phase collective synchronization between the groups even though all microscopic interactions are in-phase coupling. This counter-intuitive phenomenon can occur even when the number of oscillators belonging to each group is only two, that is, when the total number of oscillators is only four. In this paper, we clarify the mechanism underlying this counter-intuitive phenomenon for two weakly interacting groups of two oscillators with global sinusoidal coupling.

Stationary and traveling wave states of the Kuramoto model with an arbitrary distribution of frequencies and coupling strengths

We consider the Kuramoto model of an ensemble of interacting oscillators allowing for an arbitrary distribution of frequencies and coupling strengths. We define a family of traveling wave states as stationary in a rotating frame, and derive general equations for their parameters. We suggest empirical stability conditions which, for the case of incoherence, become exact. In addition to making new theoretical predictions, we show that many earlier results follow naturally from our general framework. The results are applicable in scientific contexts ranging from physics to biology.

Mean-field approximation of two coupled populations of excitable units

The analysis on stability and bifurcations in the macroscopic dynamics exhibited by the system of two coupled large populations composed of N stochastic excitable units each is performed by studying an approximate system, obtained by replacing each population with the corresponding mean-field model. In the exact system, one has the units within an ensemble communicating via the time-delayed linear couplings, whereas the interensemble terms involve the nonlinear time-delayed interaction mediated by the appropriate global variables. The aim is to demonstrate that the bifurcations affecting the stability of the stationary state of the original system, governed by a set of 4N stochastic delay-differential equations for the microscopic dynamics, can accurately be reproduced by a flow containing just four deterministic delay-differential equations which describe the evolution of the mean-field based variables. In particular, the considered issues include determining the parameter domains where the stationary state is stable, the scenarios for the onset, and the time-delay induced suppression of the collective mode, as well as the parameter domains admitting bistability between the equilibrium and the oscillatory state. We show how analytically tractable bifurcations occurring in the approximate model can be used to identify the characteristic mechanisms by which the stationary state is destabilized under different system configurations, like those with symmetrical or asymmetrical interpopulation couplings.

Disorder-induced dynamics in a pair of coupled heterogeneous phase oscillator networks

We consider a pair of coupled heterogeneous phase oscillator networks and investigate their dynamics in the continuum limit as the intrinsic frequencies of the oscillators are made more and more disparate. The Ott/Antonsen Ansatz is used to reduce the system to three ordinary differential equations. We find that most of the interesting dynamics, such as chaotic behaviour, can be understood by analysing a gluing bifurcation of periodic orbits; these orbits can be thought of as "breathing chimeras" in the limit of identical oscillators. We also add Gaussian white noise to the oscillators' dynamics and derive a pair of coupled Fokker-Planck equations describing the dynamics in this case. Comparison with simulations of finite networks of oscillators is used to confirm many of the results.