Remote synchronization reveals network symmetries and functional modules

Scalable synchronization cluster in networked chaotic oscillators

Cluster synchronization in synthetic networks of coupled chaotic oscillators is investigated. It is found that despite the asymmetric nature of the network structure, a subset of the oscillators can be synchronized as a cluster while the other oscillators remain desynchronized. Interestingly, with the increase in the coupling strength, the cluster is expanding gradually by recruiting the desynchronized oscillators one by one. This new synchronization phenomenon, which is named "scalable synchronization cluster," is explored theoretically by the method of eigenvector-based analysis, and it is revealed that the scalability of the cluster is attributed to the unique feature of the eigenvectors of the network coupling matrix. The transient dynamics of the cluster in response to random perturbations are also studied, and it is shown that in restoring to the synchronization state, oscillators inside the cluster are stabilized in sequence, illustrating again the hierarchy of the oscillators. The findings shed new light on the collective behaviors of networked chaotic oscillators and are helpful for the design of real-world networks where scalable synchronization clusters are concerned.

The transition to synchronization of networked systems

We study the synchronization properties of a generic networked dynamical system, and show that, under a suitable approximation, the transition to synchronization can be predicted with the only help of eigenvalues and eigenvectors of the graph Laplacian matrix. The transition comes out to be made of a well defined sequence of events, each of which corresponds to a specific clustered state. The network's nodes involved in each of the clusters can be identified, and the value of the coupling strength at which the events are taking place can be approximately ascertained. Finally, we present large-scale simulations which show the accuracy of the approximation made, and of our predictions in describing the synchronization transition of both synthetic and real-world large size networks, and we even report that the observed sequence of clusters is preserved in heterogeneous networks made of slightly non-identical systems.

Multi-type synchronization for coupled van der Pol oscillator systems with multiple coupling modes

In this paper, we investigate synchronous solutions of coupled van der Pol oscillator systems with multiple coupling modes using the theory of rotating periodic solutions. Multiple coupling modes refer to two or three types of coupling modes in van der Pol oscillator networks, namely, position, velocity, and acceleration. Rotating periodic solutions can represent various types of synchronous solutions corresponding to different phase differences of coupled oscillators. When matrices representing the topology of different coupling modes have symmetry, the overall symmetry of the oscillator system depends on the intersection of the symmetries of the different topologies, determining the type of synchronous solutions for the coupled oscillator network. When matrices representing the topology of different coupling modes lack symmetry, if the adjacency matrices representing different coupling modes can be simplified into structurally identical quotient graphs (where weights can be proportional) through the same external equitable partition, the symmetry of the quotient graph determines the synchronization type of the original system. All these results are consistent with multi-layer networks where connections between different layers are one-to-one.

Does the brain behave like a (complex) network? I. Dynamics

Graph theory is now becoming a standard tool in system-level neuroscience. However, endowing observed brain anatomy and dynamics with a complex network structure does not entail that the brain actually works as a network. Asking whether the brain behaves as a network means asking whether network properties count. From the viewpoint of neurophysiology and, possibly, of brain physics, the most substantial issues a network structure may be instrumental in addressing relate to the influence of network properties on brain dynamics and to whether these properties ultimately explain some aspects of brain function. Here, we address the dynamical implications of complex network, examining which aspects and scales of brain activity may be understood to genuinely behave as a network. To do so, we first define the meaning of networkness, and analyse some of its implications. We then examine ways in which brain anatomy and dynamics can be endowed with a network structure and discuss possible ways in which network structure may be shown to represent a genuine organisational principle of brain activity, rather than just a convenient description of its anatomy and dynamics.

Breathing cluster in complex neuron-astrocyte networks

Brain activities are featured by spatially distributed neural clusters of coherent firings and a spontaneous slow switching of the clusters between the coherent and incoherent states. Evidences from recent in vivo experiments suggest that astrocytes, a type of glial cell regarded previously as providing only structural and metabolic supports to neurons, participate actively in brain functions by regulating the neural firing activities, yet the underlying mechanism remains unknown. Here, introducing astrocyte as a reservoir of the glutamate released from the neuron synapses, we propose the model of the complex neuron-astrocyte network, and investigate the roles of astrocytes in regulating the cluster synchronization behaviors of networked chaotic neurons. It is found that a specific set of neurons on the network are synchronized and form a cluster, while the remaining neurons are kept as desynchronized. Moreover, during the course of network evolution, the cluster is switching between the synchrony and asynchrony states in an intermittent fashion, henceforth the phenomenon of "breathing cluster." By the method of symmetry-based analysis, we conduct a theoretical investigation on the synchronizability of the cluster. It is revealed that the contents of the cluster are determined by the network symmetry, while the breathing of the cluster is attributed to the interplay between the neural network and the astrocyte. The phenomenon of breathing cluster is demonstrated in different network models, including networks with different sizes, nodal dynamics, and coupling functions. The findings shed light on the cellular mechanism of astrocytes in regulating neural activities and give insights into the state-switching of the neocortex.

Endowing networks with desired symmetries and modular behavior

Symmetries in a network regulate its organization into functional clustered states. Given a generic ensemble of nodes and a desirable cluster (or group of clusters), we exploit the direct connection between the elements of the eigenvector centrality and the graph symmetries to generate a network equipped with the desired cluster(s), with such a synthetical structure being furthermore perfectly reflected in the modular organization of the network's functioning. Our results solve a relevant problem of designing a desired set of clusters and are of generic application in all cases where a desired parallel functioning needs to be blueprinted.

Noisy image segmentation based on synchronous dynamics of coupled photonic spiking neurons

The collective dynamics in neural networks is essential for information processing and has attracted much interest on the application in artificial intelligence. Synchronization is one of the most dominant phenomenon in the collective dynamics of neural network. Here, we propose to use the spiking dynamics and collective synchronization of coupled photonic spiking neurons for noisy image segmentation. Based on the synchronization mechanism and synchronization control, the noised pattern segmentation is demonstrated numerically. This work provides insight into the possible application based on the collective dynamics of large-scale photonic networks and opens a way for ultra-high speed image processing.

The co-presence of polystyrene nanoplastics and ofloxacin demonstrates combined effects on the structure, assembly, and metabolic activities of marine microbial community

Nanoplastic is increasing in environments and can address toxic effects on various organisms. Particle size, concentration, and surface functionalization most influence nanoplastic toxicity. Besides, nanoplastic can adsorb other contaminants (e.g., antibiotics) to aggravate its adverse effects. The combined effects of nanoplastics and antibiotics on planktonic/benthic microbial communities, however, are still largely unknown. In this study, the combined effects of polystyrene nanoplastic and ofloxacin on the structure, assembly, and metabolic activities of marine microbial communities were investigated based on amplicon sequencing data. The results mainly demonstrate that: (1) nanoplastic and ofloxacin have greater impacts on prokaryotic communities than eukaryotic ones; (2) niche breadths of planktonic prokaryotes and benthic eukaryotes were shrank with both high nanoplastic and ofloxacin concentrations; (3) increased ofloxacin mainly reduces nodes/edges of co-occurrence networks, while nanoplastic centralizes network modularity; (4) increased nanoplastic under high ofloxacin concentration induces more differential prokaryotic pathways in planktonic communities, while benthic communities are less influenced. The present work indicates that co-presence of nanoplastics and ofloxacin has synergistic combined effects on community structure shifts, niche breadth shrinking, network simplifying, and differential prokaryotic pathways inducing in marine microbial communities, suggesting nanoplastics and its combined impacts with other pollutions should be paid with more concerns.

General optimization framework for accurate and efficient reconstruction of symmetric complex networks from dynamical data

The challenging problem of network reconstruction from dynamical data can in general be formulated as an optimization task of solving multiple linear equations. Existing approaches are of the two types: Point-by-point (PBP) and global methods. The local PBP method is computationally efficient, but the accuracies of its solutions are somehow low, while a global method has the opposite traits: High accuracy and high computational cost. Taking advantage of the network symmetry, we develop a novel framework integrating the advantages of both the PBP and global methods while avoiding their shortcomings: i.e., high reconstruction accuracy is guaranteed, but the computational cost is orders of magnitude lower than that of the global methods in the literature. The mathematical principle underlying our framework is block coordinate descent (BCD) for solving optimization problems, where the various blocks are determined by the network symmetry. The reconstruction framework is validated by numerical examples with a variety of network structures (i.e., sparse and dense networks) and dynamical processes. Our success is a demonstration that the general principle of exploiting symmetry can be extended to tackling the challenging inverse problem or reverse engineering of complex networks. Since solving a large number of linear equations is key to a plethora of problems in science and engineering, our BCD-based network reconstruction framework will find broader applications.

Detecting partial synchrony in a complex oscillatory network using pseudovortices

Partial synchronization is an important dynamical process of coupled oscillators on various natural and artificial networks, which can remain undetected due to the system complexity. With an analogy between pairwise asynchrony of oscillators and topological defects, i.e., vortices, in the two-dimensional XY model, we propose a robust and data-driven method to identify the partial synchronization on complex networks. The proposed method is based on an integer matrix whose element is pseudovorticity that discretely quantifies asynchronous phase dynamics in every two oscillators, which results in graphical and entropic representations of partial synchrony. As a first trial, we apply our method to 200 FitzHugh-Nagumo neurons on a complex small-world network. Partially synchronized chimera states are revealed by discriminating synchronized states even with phase lags. Such phase lags also appear in partial synchronization in chimera states. Our topological, graphical, and entropic method is implemented solely with measurable phase dynamics data, which will lead to a straightforward application to general oscillatory networks including neural networks in the brain.

Transitional cluster dynamics in a model for delay-coupled chemical oscillators

Cluster synchronization is a fundamental phenomenon in systems of coupled oscillators. Here, we investigate clustering patterns that emerge in a unidirectional ring of four delay-coupled electrochemical oscillators. A voltage parameter in the experimental setup controls the onset of oscillations via a Hopf bifurcation. For a smaller voltage, the oscillators exhibit simple, so-called primary, clustering patterns, where all phase differences between each set of coupled oscillators are identical. However, upon increasing the voltage, secondary states, where phase differences differ, are detected, in addition to the primary states. Previous work on this system saw the development of a mathematical model that explained how the existence, stability, and common frequency of the experimentally observed cluster states could be accurately controlled by the delay time of the coupling. In this study, we revisit the mathematical model of the electrochemical oscillators in order to address open questions by means of bifurcation analysis. Our analysis reveals how the stable cluster states, corresponding to experimental observations, lose their stability via an assortment of bifurcation types. The analysis further reveals complex interconnectedness between branches of different cluster types. We find that each secondary state provides a continuous transition between certain primary states. These connections are explained by studying the phase space and parameter symmetries of the respective states. Furthermore, we show that it is only for a larger value of the voltage parameter that the branches of secondary states develop intervals of stability. For a smaller voltage, all the branches of secondary states are completely unstable and are, therefore, hidden to experimentalists.

Collective behavior of identical Stuart-Landau oscillators in a star network with coupling asymmetry effects

In this study, we investigated the impact of the asymmetry of a coupling scheme on oscillator dynamics in a star network. We obtained stability conditions for the collective behavior of the systems, ranging from an equilibrium point over complete synchronization (CS) and quenched hub incoherence to remote synchronization states using both numerical and analytical methods. The coupling asymmetry factor α significantly influences and determines the stable parameter region of each state. For α ≠ 1, the equilibrium point can emerge when the Hopf bifurcation parameter a is positive, which is impossible for diffusive coupling. However, CS can occur even if a is negative under α < 1. Unlike diffusive coupling, we observe more behavior when α ≠ 1, including additional in-phase remote synchronization. These results are supported by theoretical analysis and validated through numerical simulations and independent of network size. The findings may offer practical methods for controlling, restoring, or obstructing specific collective behavior.

Percolation critical exponents in cluster kinetics of pulse-coupled oscillators

Transient dynamics leading to the synchrony of a type of pulse-coupled oscillators, so-called scrambler oscillators, has previously been studied as an aggregation process of synchronous clusters, and a rate equation for the cluster size distribution has been proposed. However, the evolution of the cluster size distribution for general cluster sizes has not been fully understood yet. In this paper, we study the evolution of the cluster size distribution from the perspective of a percolation model by regarding the number of aggregations as the number of attached bonds. Specifically, we derive the scaling form of the cluster size distribution with specific values of the critical exponents using the property that the characteristic cluster size diverges as the percolation threshold is approached from below. Through simulation, it is confirmed that the scaling form well explains the evolution of the cluster size distribution. Based on the distribution behavior, we find that a giant cluster of all oscillators is formed discontinuously at the threshold and also that further aggregation does not occur like in a one-dimensional bond percolation model. Finally, we discuss the origin of the discontinuous formation of the giant cluster from the perspective of global suppression in explosive percolation models. For this, we approximate the aggregation process as a cluster-cluster aggregation with a given collision kernel. We believe that the theoretical approach presented in this paper can be used to understand the transient dynamics of a broad range of synchronizations.

Phase frustration induced remote synchronization

Remote synchronization (RS) may take an important role in brain functioning and its study has attracted much attention in recent years. So far, most studies of RS are focused on the Stuart-Landau oscillators with mean-field coupling. However, realistic cases may have more complicated couplings and behaviors, such as the brain networks. To make the study of RS a substantial progress toward realistic situations, we here present a model of RS with phase frustration and show that RS can be induced for those systems where no RS exists when there is no phase frustration. By numerical simulations on both the Stuart-Landau and Kuramoto oscillators, we find that the optimal range of RS depends on the match of phase frustrations between the hub and leaf nodes and a fixed relationship of this match is figured out. While for the non-optimal range of RS, we find that RS exists only in a linear band between the phase frustrations of the hub and leaf nodes. A brief theoretical analysis is provided to explain these results.

A multiplex analysis of phonological and orthographic networks

The study of natural language using a network approach has made it possible to characterize novel properties ranging from the level of individual words to phrases or sentences. A natural way to quantitatively evaluate similarities and differences between spoken and written language is by means of a multiplex network defined in terms of a similarity distance between words. Here, we use a multiplex representation of words based on orthographic or phonological similarity to evaluate their structure. We report that from the analysis of topological properties of networks, there are different levels of local and global similarity when comparing written vs. spoken structure across 12 natural languages from 4 language families. In particular, it is found that differences between the phonetic and written layers is markedly higher for French and English, while for the other languages analyzed, this separation is relatively smaller. We conclude that the multiplex approach allows us to explore additional properties of the interaction between spoken and written language.

Structural position vectors and symmetries in complex networks

Symmetries, due to their fundamental importance to dynamical processes on networks, have attracted a great deal of current research. Finding all symmetric nodes in large complex networks typically relies on automorphism groups from algebraic-group theory, which are solvable in quasipolynomial time. We articulate a conceptually appealing and computationally extremely efficient approach to finding and characterizing all symmetric nodes by introducing a structural position vector (SPV) for each node in networks. We establish the mathematical result that symmetric nodes must have the same SPV value and demonstrate, using six representative complex networks from the real world, that all symmetric nodes in these networks can be found in linear time. Furthermore, the SPVs not only characterize the similarity of nodes but also quantify the nodal influences in propagation dynamics. A caveat is that the proved mathematical result relating the SPV values to nodal symmetries is not sufficient; i.e., nodes having the same SPV values may not be symmetric, which arises in regular networks or networks with a dominant regular component. We point out with an analysis that this caveat is, in fact, shared by the known existing approaches to finding symmetric nodes in the literature. We further argue, with the aid of a mathematical analysis, that our SPV method is generally effective for finding the symmetric nodes in real-world networks that typically do not have a dominant regular component. Our SPV-based framework, therefore, provides a physically intuitive and computationally efficient way to uncover, understand, and exploit symmetric structures in complex networks arising from real-world applications.

Modelling the perception of music in brain network dynamics

We analyze the influence of music in a network of FitzHugh-Nagumo oscillators with empirical structural connectivity measured in healthy human subjects. We report an increase of coherence between the global dynamics in our network and the input signal induced by a specific music song. We show that the level of coherence depends crucially on the frequency band. We compare our results with experimental data, which also describe global neural synchronization between different brain regions in the gamma-band range in a time-dependent manner correlated with musical large-scale form, showing increased synchronization just before transitions between different parts in a musical piece (musical high-level events). The results also suggest a separation in musical form-related brain synchronization between high brain frequencies, associated with neocortical activity, and low frequencies in the range of dance movements, associated with interactivity between cortical and subcortical regions.

Supermodal Decomposition of the Linear Swing Equation for Multilayer Networks

We study the swing equation in the case of a multilayer network in which generators and motors are modeled differently; namely, the model for each generator is given by second order dynamics and the model for each motor is given by first order dynamics. We also remove the commonly used assumption of equal damping coefficients in the second order dynamics. Under these general conditions, we are able to obtain a decomposition of the linear swing equation into independent modes describing the propagation of small perturbations. In the process, we identify symmetries affecting the structure and dynamics of the multilayer network and derive an essential model based on a 'quotient network.' We then compare the dynamics of the full network and that of the quotient network and obtain a modal decomposition of the error dynamics. We also provide a method to quantify the steady-state error and the maximum overshoot error. Two case studies are presented to illustrate application of our method.

Chimera Patterns of Synchrony in a Frustrated Array of Hebb Synapses

The union of the Kuramoto–Sakaguchi model and the Hebb dynamics reproduces the Lisman switch through a bistability in synchronized states. Here, we show that, within certain ranges of the frustration parameter, the chimera pattern can emerge, causing a different, time-evolving, distribution in the Hebbian synaptic strengths. We study the stability range of the chimera as a function of the frustration (phase-lag) parameter. Depending on the range of the frustration, two different types of chimeras can appear spontaneously, i.e., from randomized initial conditions. In the first type, the oscillators in the coherent region rotate, on average, slower than those in the incoherent region; while in the second type, the average rotational frequencies of the two regions are reversed, i.e., the coherent region runs, on average, faster than the incoherent region. We also show that non-stationary behavior at finite N can be controlled by adjusting the natural frequency of a single pacemaker oscillator. By slowly cycling the frequency of the pacemaker, we observe hysteresis in the system. Finally, we discuss how we can have a model for learning and memory.

Interlayer Connectivity Affects the Coherence Resonance and Population Activity Patterns in Two-Layered Networks of Excitatory and Inhibitory Neurons

The firing patterns of neuronal populations often exhibit emergent collective oscillations, which can display substantial regularity even though the dynamics of individual elements is very stochastic. One of the many phenomena that is often studied in this context is coherence resonance, where additional noise leads to improved regularity of spiking activity in neurons. In this work, we investigate how the coherence resonance phenomenon manifests itself in populations of excitatory and inhibitory neurons. In our simulations, we use the coupled FitzHugh-Nagumo oscillators in the excitable regime and in the presence of neuronal noise. Formally, our model is based on the concept of a two-layered network, where one layer contains inhibitory neurons, the other excitatory neurons, and the interlayer connections represent heterotypic interactions. The neuronal activity is simulated in realistic coupling schemes in which neurons within each layer are connected with undirected connections, whereas neurons of different types are connected with directed interlayer connections. In this setting, we investigate how different neurophysiological determinants affect the coherence resonance. Specifically, we focus on the proportion of inhibitory neurons, the proportion of excitatory interlayer axons, and the architecture of interlayer connections between inhibitory and excitatory neurons. Our results reveal that the regularity of simulated neural activity can be increased by a stronger damping of the excitatory layer. This can be accomplished with a higher proportion of inhibitory neurons, a higher fraction of inhibitory interlayer axons, a stronger coupling between inhibitory axons, or by a heterogeneous configuration of interlayer connections. Our approach of modeling multilayered neuronal networks in combination with stochastic dynamics offers a novel perspective on how the neural architecture can affect neural information processing and provide possible applications in designing networks of artificial neural circuits to optimize their function via noise-induced phenomena.

Metamorphoses and explosively remote synchronization in dynamical networks

We uncover a phenomenon in coupled nonlinear networks with a symmetry: as a bifurcation parameter changes through a critical value, synchronization among a subset of nodes can deteriorate abruptly, and, simultaneously, perfect synchronization emerges suddenly among a different subset of nodes that are not directly connected. This is a synchronization metamorphosis leading to an explosive transition to remote synchronization. The finding demonstrates that an explosive onset of synchrony and remote synchronization, two phenomena that have been studied separately, can arise in the same system due to symmetry, providing another proof that the interplay between nonlinear dynamics and symmetry can lead to a surprising phenomenon in physical systems.

Matryoshka and disjoint cluster synchronization of networks

The main motivation for this paper is to characterize network synchronizability for the case of cluster synchronization (CS), in an analogous fashion to Barahona and Pecora [Phys. Rev. Lett. 89, 054101 (2002)] for the case of complete synchronization. We find this problem to be substantially more complex than the original one. We distinguish between the two cases of networks with intertwined clusters and no intertwined clusters and between the two cases that the master stability function is negative either in a bounded range or in an unbounded range of its argument. Our proposed definition of cluster synchronizability is based on the synchronizability of each individual cluster within a network. We then attempt to generalize this definition to the entire network. For CS, the synchronous solution for each cluster may be stable, independent of the stability of the other clusters, which results in possibly different ranges in which each cluster synchronizes (isolated CS). For each pair of clusters, we distinguish between three different cases: Matryoshka cluster synchronization (when the range of the stability of the synchronous solution for one cluster is included in that of the other cluster), partially disjoint cluster synchronization (when the ranges of stability of the synchronous solutions partially overlap), and complete disjoint cluster synchronization (when the ranges of stability of the synchronous solutions do not overlap).

Looking beyond community structure leads to the discovery of dynamical communities in weighted networks

A fundamental question is whether groups of nodes of a complex network can possibly display long-term cluster-synchronized behavior. While this question has been addressed for the restricted classes of unweighted and labeled graphs, it remains an open problem for the more general class of weighted networks. The emergence of coordinated motion of nodes in natural and technological networks is directly related to the network structure through the concept of an equitable partition, which determines which nodes can show long-term synchronized behavior and which nodes cannot. We provide a method to detect the presence of nearly equitable partitions in weighted networks, based on minimal information about the network structure. With this approach we are able to discover the presence of dynamical communities in both synthetic and real technological, biological, and social networks, to a statistically significant level. We show that our approach based on dynamical communities is better at predicting the emergence of synchronized behavior than existing methods to detect community structure.

Pattern Discovery in Multilayer Networks

MOTIVATION

In bioinformatics, complex cellular modeling and behavior simulation to identify significant molecular interactions is considered a relevant problem. Traditional methods model such complex systems using single and binary network. However, this model is inadequate to represent biological networks as different sets of interactions can simultaneously take place for different interaction constraints (such as transcription regulation and protein interaction). Furthermore, biological systems may exhibit varying interaction topologies even for the same interaction type under different developmental stages or stress conditions. Therefore, models which consider biological systems as solitary interactions are inaccurate as they fail to capture the complex behavior of cellular interactions within organisms. Identification and counting of recurrent motifs within a network is one of the fundamental problems in biological network analysis. Existing methods for motif counting on single network topologies are inadequate to capture patterns of molecular interactions that have significant changes in biological expression when identified across different organisms that are similar, or even time-varying networks within the same organism. That is, they fail to identify recurrent interactions as they consider a single snapshot of a network among a set of multiple networks. Therefore, we need methods geared towards studying multiple network topologies and the pattern conservation among them. Contributions: In this paper, we consider the problem of counting the number of instances of a user supplied motif topology in a given multilayer network. We model interactions among a set of entities (e.g., genes)describing various conditions or temporal variation as multilayer networks. Thus a separate network as each layer shows the connectivity of the nodes under a unique network state. Existing motif counting and identification methods are limited to single network topologies, and thus cannot be directly applied on multilayer networks. We apply our model and algorithm to study frequent patterns in cellular networks that are common in varying cellular states under different stress conditions, where the cellular network topology under each stress condition describes a unique network layer.

RESULTS

We develop a methodology and corresponding algorithm based on the proposed model for motif counting in multilayer networks. We performed experiments on both real and synthetic datasets. We modeled the synthetic datasets under a wide spectrum of parameters, such as network size, density, motif frequency. Results on synthetic datasets demonstrate that our algorithm finds motif embeddings with very high accuracy compared to existing state-of-the-art methods such as G-tries, ESU (FANMODE)and mfinder. Furthermore, we observe that our method runs from several times to several orders of magnitude faster than existing methods. For experiments on real dataset, we consider Escherichia coli (E. coli)transcription regulatory network under different experimental conditions. We observe that the genes selected by our method conserves functional characteristics under various stress conditions with very low false discovery rates. Moreover, the method is scalable to real networks in terms of both network size and number of layers.

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.

Symmetry-breaking mechanism for the formation of cluster chimera patterns

The emergence of order in collective dynamics is a fascinating phenomenon that characterizes many natural systems consisting of coupled entities. Synchronization is such an example where individuals, usually represented by either linear or nonlinear oscillators, can spontaneously act coherently with each other when the interactions' configuration fulfills certain conditions. However, synchronization is not always perfect, and the coexistence of coherent and incoherent oscillators, broadly known in the literature as chimera states, is also possible. Although several attempts have been made to explain how chimera states are created, their emergence, stability, and robustness remain a long-debated question. We propose an approach that aims to establish a robust mechanism through which cluster synchronization and chimera patterns originate. We first introduce a stability-breaking method where clusters of synchronized oscillators can emerge. At variance with the standard approach where synchronization arises as a collective behavior of coupled oscillators, in our model, the system initially sets on a homogeneous fixed-point regime, and, only due to a global instability principle, collective oscillations emerge. Following a combination of the network modularity and the model's parameters, one or more clusters of oscillators become incoherent within yielding a particular class of patterns that we here name cluster chimera states.

Two mechanisms of remote synchronization in a chain of Stuart-Landau oscillators

Remote synchronization implies that oscillators interacting not directly but via an additional unit (hub) adjust their frequencies and exhibit frequency locking while the hub remains asynchronous. In this paper, we analyze the mechanisms of remote synchrony in a small network of three coupled Stuart-Landau oscillators using recent results on higher-order phase reduction. We analytically demonstrate the role of two factors promoting remote synchrony. These factors are the nonisochronicity of oscillators and the coupling terms appearing in the second-order phase approximation. We show a good correspondence between our theory and numerical results for small and moderate coupling strengths.

Cluster synchronization of networks via a canonical transformation for simultaneous block diagonalization of matrices

We study cluster synchronization of networks and propose a canonical transformation for simultaneous block diagonalization of matrices that we use to analyze the stability of the cluster synchronous solution. Our approach has several advantages as it allows us to: (1) decouple the stability problem into subproblems of minimal dimensionality while preserving physically meaningful information, (2) study stability of both orbital and equitable partitions of the network nodes, and (3) obtain a parameterization of the problem in a small number of parameters. For the last point, we show how the canonical transformation decouples the problem into blocks that preserve key physical properties of the original system. We also apply our proposed algorithm to analyze several real networks of interest, and we find that it runs faster than alternative algorithms from the literature.

Attracting Poisson chimeras in two-population networks

Chimera states, i.e., dynamical states composed of coexisting synchronous and asynchronous oscillations, have been reported to exist in diverse topologies of oscillators in simulations and experiments. Two-population networks with distinct intra- and inter-population coupling have served as simple model systems for chimera states since they possess an invariant synchronized manifold in contrast to networks on a spatial structure. Here, we study dynamical and spectral properties of finite-sized chimeras on two-population networks. First, we elucidate how the Kuramoto order parameter of the finite-sized globally coupled two-population network of phase oscillators is connected to that of the continuum limit. These findings suggest that it is suitable to classify the chimera states according to their order parameter dynamics, and therefore, we define Poisson and non-Poisson chimera states. We then perform a Lyapunov analysis of these two types of chimera states, which yields insight into the full stability properties of the chimera trajectories as well as of collective modes. In particular, our analysis also confirms that Poisson chimeras are neutrally stable. We then introduce two types of "perturbation" that act as small heterogeneities and render Poisson chimeras attracting: A topological variation via the simplest nonlocal intra-population coupling that keeps the network symmetries and the allowance of amplitude variations in the globally coupled two-population network; i.e., we replace the phase oscillators by Stuart-Landau oscillators. The Lyapunov spectral properties of chimera states in the two modified networks are investigated, exploiting an approach based on network symmetry-induced cluster pattern dynamics of the finite-size network.

Rapid exponential stabilization of Lotka-McKendrick's equation via event-triggered impulsive control

This paper investigates the problem of rapid exponential stabilization for linear Lotka-McKendrick's equation. Based on a new event-triggered impulsive control (ETIC) method, an impulsive control is designed to solve the rapid exponential stabilization of the dynamic population Lotka-McKendrick's equation. The effectiveness of our control is verified through a numerical example.

Finite-size scaling of geometric renormalization flows in complex networks

Some characteristics of complex networks need to be derived from global knowledge of the network topologies, which challenges the practice for studying many large-scale real-world networks. Recently, the geometric renormalization technique has provided a good approximation framework to significantly reduce the size and complexity of a network while retaining its "slow" degrees of freedom. However, due to the finite-size effect of real networks, excessive renormalization iterations will eventually cause these important "slow" degrees of freedom to be filtered out. In this paper, we systematically investigate the finite-size scaling of structural and dynamical observables in geometric renormalization flows of both synthetic and real evolutionary networks. Our results show that these observables can be well characterized by a certain scaling function. Specifically, we show that the critical exponent implied by the scaling function is independent of these observables but depends only on the structural properties of the network. To a certain extent, the results of this paper are of great significance for predicting the observable quantities of large-scale real systems and further suggest that the potential scale invariance of many real-world networks is often masked by finite-size effects.

Dynamics of coupled Kuramoto oscillators with distributed delays

This paper studies the effects of two different types of distributed-delay coupling in the system of two mutually coupled Kuramoto oscillators: one where the delay distribution is considered inside the coupling function and the other where the distribution enters outside the coupling function. In both cases, the existence and stability of phase-locked solutions is analyzed for uniform and gamma distribution kernels. The results show that while having the distribution inside the coupling function only changes parameter regions where phase-locked solutions exist, when the distribution is taken outside the coupling function, it affects both the existence, as well as stability properties of in- and anti-phase states. For both distribution types, various branches of phase-locked solutions are computed, and regions of their stability are identified for uniform, weak, and strong gamma distributions.

Asymmetry-induced order in multilayer networks

Symmetries naturally occur in real-world networks and can significantly influence the observed dynamics. For instance, many synchronization patterns result from the underlying network symmetries, and high symmetries are known to increase the stability of synchronization. Yet here we find that general macroscopic features of network solutions such as regularity can be induced by breaking their symmetry of interactions. We demonstrate this effect in an ecological multilayer network where the topological asymmetries occur naturally. These asymmetries rescue the system from chaotic oscillations by establishing stable periodic orbits and equilibria. We call this phenomenon asymmetry-induced order and uncover its mechanism by analyzing both analytically and numerically the absence of dynamics on the system's synchronization manifold. Moreover, the bifurcation scenario describing the route from chaos to order is also disclosed. We demonstrate that this result also holds for generic node dynamics by analyzing coupled paradigmatic Rössler and Lorenz systems.

Synchronization scenarios in three-layer networks with a hub

We study various relay synchronization scenarios in a three-layer network, where the middle (relay) layer is a single node, i.e., a hub. The two remote layers consist of non-locally coupled rings of FitzHugh-Nagumo oscillators modeling neuronal dynamics. All nodes of the remote layers are connected to the hub. The role of the hub and its importance for the existence of chimera states are investigated in dependence on the inter-layer coupling strength and inter-layer time delay. Tongue-like regions in the parameter plane exhibiting double chimeras, i.e., chimera states in the remote layers whose coherent cores are synchronized with each other, and salt-and-pepper states are found. At very low intra-layer coupling strength, when chimera states do not exist in single layers, these may be induced by the hub. Also, the influence of the dilution of links between the remote layers and the hub upon the dynamics is investigated. The greatest effect of dilution is observed when links to the coherent domain of the chimeras are removed.

Generalized splay states in phase oscillator networks

Networks of coupled phase oscillators play an important role in the analysis of emergent collective phenomena. In this article, we introduce generalized m-splay states constituting a special subclass of phase-locked states with vanishing mth order parameter. Such states typically manifest incoherent dynamics, and they often create high-dimensional families of solutions (splay manifolds). For a general class of phase oscillator networks, we provide explicit linear stability conditions for splay states and exemplify our results with the well-known Kuramoto-Sakaguchi model. Importantly, our stability conditions are expressed in terms of just a few observables such as the order parameter or the trace of the Jacobian. As a result, these conditions are simple and applicable to networks of arbitrary size. We generalize our findings to phase oscillators with inertia and adaptively coupled phase oscillator models.

Network isolators inhibit failure spreading in complex networks

In our daily lives, we rely on the proper functioning of supply networks, from power grids to water transmission systems. A single failure in these critical infrastructures can lead to a complete collapse through a cascading failure mechanism. Counteracting strategies are thus heavily sought after. In this article, we introduce a general framework to analyse the spreading of failures in complex networks and demostrate that not only decreasing but also increasing the connectivity of the network can be an effective method to contain damages. We rigorously prove the existence of certain subgraphs, called network isolators, that can completely inhibit any failure spreading, and we show how to create such isolators in synthetic and real-world networks. The addition of selected links can thus prevent large scale outages as demonstrated for power transmission grids.

Stability of synchronization in simplicial complexes

Various systems in physics, biology, social sciences and engineering have been successfully modeled as networks of coupled dynamical systems, where the links describe pairwise interactions. This is, however, too strong a limitation, as recent studies have revealed that higher-order many-body interactions are present in social groups, ecosystems and in the human brain, and they actually affect the emergent dynamics of all these systems. Here, we introduce a general framework to study coupled dynamical systems accounting for the precise microscopic structure of their interactions at any possible order. We show that complete synchronization exists as an invariant solution, and give the necessary condition for it to be observed as a stable state. Moreover, in some relevant instances, such a necessary condition takes the form of a Master Stability Function. This generalizes the existing results valid for pairwise interactions to the case of complex systems with the most general possible architecture.

Optimal cost tuning of frustration: Achieving desired states in the Kuramoto-Sakaguchi model

There are numerous examples of studied real-world systems that can be described as dynamical systems characterized by individual phases and coupled in a networklike structure. Within the framework of oscillatory models, much attention has been devoted to the Kuramoto model, which considers a collection of oscillators interacting through a sinus function of the phase differences. In this paper, we draw on an extension of the Kuramoto model, called the Kuramoto-Sakaguchi model, which adds a phase lag parameter to each node. We construct a general formalism that allows us to compute the set of lag parameters that may lead to any phase configuration within a linear approximation. In particular, we devote special attention to the cases of full synchronization and symmetric configurations. We show that the set of natural frequencies, phase lag parameters, and phases at the steady state is coupled by an equation and a continuous spectra of solutions is feasible. In order to quantify the system's strain to achieve that particular configuration, we define a cost function and compute the optimal set of parameters that minimizes it. Despite considering a linear approximation of the model, we show that the obtained tuned parameters for the case of full synchronization enhance frequency synchronization in the nonlinear model as well.

Comparing phonological and orthographic networks: A multiplex analysis

The complexity of natural language can be explored by means of multiplex analyses at different scales, from single words to groups of words or sentence levels. Here, we plan to investigate a multiplex word-level network, which comprises an orthographic and a phonological network defined in terms of distance similarity. We systematically compare basic structural network properties to determine similarities and differences between them, as well as their combination in a multiplex configuration. As a natural extension of our work, we plan to evaluate the preservation of the structural network properties and information-based quantities from the following perspectives: (i) presence of similarities across 12 natural languages from 4 linguistic families (Romance, Germanic, Slavic and Uralic), (ii) increase of the size of the number of words (corpus) from 104 to 50 × 103, and (iii) robustness of the networks. Our preliminary findings reinforce the idea of common organizational properties among natural languages. Once concluded, will contribute to the characterization of similarities and differences in the orthographic and phonological perspectives of language networks at a word-level.

Delays induced cluster synchronization in chaotic networks

We study networks of coupled oscillators and analyze the role of coupling delays in determining the emergence of cluster synchronization. Given a network topology and a particular arrangement of the coupling delays over the network connections, different patterns of cluster synchronization may emerge. We focus on a simple ring network of six bidirectionally coupled identical oscillators, for which with two different values of the delays, a total of eight cluster synchronization patterns may emerge, depending on the assignment of the delays to the ring connections. We analyze stability of each of the patterns and find that for large enough coupling strength and specific values of the delays, they can all be stabilized. We construct an experimental ring of six bidirectionally coupled Colpitts oscillators, with delayed connections obtained by coupling the oscillators via RF cables of appropriate length. We find that experimental observations of cluster synchronization are in essential agreement with theoretical predictions. We also verify our theory in a fully connected network of fifty nodes for which connections are randomly assigned to be either undelayed or delayed with a given probability.

Synchronization within synchronization: transients and intermittency in ecological networks

Transients are fundamental to ecological systems with significant implications to management, conservation and biological control. We uncover a type of transient synchronization behavior in spatial ecological networks whose local dynamics are of the chaotic, predator–prey type. In the parameter regime where there is phase synchronization among all the patches, complete synchronization (i.e. synchronization in both phase and amplitude) can arise in certain pairs of patches as determined by the network symmetry—henceforth the phenomenon of  ‘synchronization within synchronization.’ Distinct patterns of complete synchronization coexist but, due to intrinsic instability or noise, each pattern is a transient and there is random, intermittent switching among the patterns in the course of time evolution. The probability distribution of the transient time is found to follow an algebraic scaling law with a divergent average transient lifetime. Based on symmetry considerations, we develop a stability analysis to understand these phenomena. The general principle of symmetry can also be exploited to explain previously discovered, counterintuitive synchronization behaviors in ecological networks.

Symmetries and cluster synchronization in multilayer networks

Real-world systems in epidemiology, social sciences, power transportation, economics and engineering are often described as multilayer networks. Here we first define and compute the symmetries of multilayer networks, and then study the emergence of cluster synchronization in these networks. We distinguish between independent layer symmetries, which occur in one layer and are independent of the other layers, and dependent layer symmetries, which involve nodes in different layers. We study stability of the cluster synchronous solution by decoupling the problem into a number of independent blocks and assessing stability of each block through a Master Stability Function. We see that blocks associated with dependent layer symmetries have a different structure to the other blocks, which affects the stability of clusters associated with these symmetries. Finally, we validate the theory in a fully analog experiment in which seven electronic oscillators of three kinds are connected with two kinds of coupling.

Effect of topology upon relay synchronization in triplex neuronal networks

Relay synchronization in complex networks is characterized by the synchronization of remote parts of the network due to their interaction via a relay. In multilayer networks, distant layers that are not connected directly can synchronize due to signal propagation via relay layers. In this work, we investigate relay synchronization of partial synchronization patterns like chimera states in three-layer networks of interacting FitzHugh-Nagumo oscillators. We demonstrate that the phenomenon of relay synchronization is robust to topological random inhomogeneities of small-world type in the layer networks. We show that including randomness in the connectivity structure either of the remote network layers or of the relay layer increases the range of interlayer coupling strength where relay synchronization can be observed.

Effect of repulsive links on frustration in attractively coupled networks

We investigate the impact of attractive-repulsive interaction in networks of limit cycle oscillators. Mainly we focus on the design principle for generating an antiphase state between adjacent nodes in a complex network. We establish that a partial negative control throughout the branches of a spanning tree inside the positively coupled limit cycle oscillators works efficiently well in comparison with randomly chosen negative links to establish zero frustration (antiphase synchronization) in bipartite graphs. Based on the emergence of zero frustration, we develop a universal 0-π rule to understand the antiphase synchronization in a bipartite graph. Further, this rule is used to construct a nonbipartite graph for a given nonzero frustrated value. We finally show the generality of 0-π rule by implementing it in arbitrary undirected nonbipartite graphs of attractive-repulsively coupled limit cycle oscillators and successfully calculate the nonzero frustration value, which matches with numerical data. The validation of the rule is checked through the bifurcation analysis of small networks. Our work may unveil the underlying mechanism of several synchronization phenomena that exist in a network of oscillators having a mixed type of coupling.

Network-induced multistability through lossy coupling and exotic solitary states

The stability of synchronised networked systems is a multi-faceted challenge for many natural and technological fields, from cardiac and neuronal tissue pacemakers to power grids. For these, the ongoing transition to distributed renewable energy sources leads to a proliferation of dynamical actors. The desynchronisation of a few or even one of those would likely result in a substantial blackout. Thus the dynamical stability of the synchronous state has become a leading topic in power grid research. Here we uncover that, when taking into account physical losses in the network, the back-reaction of the network induces new exotic solitary states in the individual actors and the stability characteristics of the synchronous state are dramatically altered. These effects will have to be explicitly taken into account in the design of future power grids. We expect the results presented here to transfer to other systems of coupled heterogeneous Newtonian oscillators.

Functionability in complex networks: Leading nodes for the transition from structural to functional networks through remote asynchronization

Complex networks are essentially heterogeneous not only in the basic properties of the constituent nodes, such as their degree, but also in the effects that these have on the global dynamical properties of the network. Networks of coupled identical phase oscillators are good examples for analyzing these effects, since an overall synchronized state can be considered a reference state. A small variation of intrinsic node parameters may cause the system to move away from synchronization, and a new phase-locked stationary state can be achieved. We propose a measure of phase dispersion that quantifies the functional response of the system to a given local perturbation. As a particular implementation, we propose a variation of the standard Kuramoto model in which the nodes of a complex network interact with their neighboring nodes, by including a node-dependent frustration parameter. The final stationary phase-locked state now depends on the particular frustration parameter at each node and also on the network topology. We exploit this scenario by introducing individual frustration parameters and measuring what their effect on the whole network is, measured in terms of the phase dispersion, which depends only on the topology of the network and on the choice of the particular node that is perturbed. This enables us to define a characteristic of the node, its functionability, that can be computed analytically in terms of the network topology. Finally, we provide a thorough comparison with other centrality measures.

Effective degrees of freedom of the Pearson's correlation coefficient under autocorrelation

The dependence between pairs of time series is commonly quantified by Pearson's correlation. However, if the time series are themselves dependent (i.e. exhibit temporal autocorrelation), the effective degrees of freedom (EDF) are reduced, the standard error of the sample correlation coefficient is biased, and Fisher's transformation fails to stabilise the variance. Since fMRI time series are notoriously autocorrelated, the issue of biased standard errors - before or after Fisher's transformation - becomes vital in individual-level analysis of resting-state functional connectivity (rsFC) and must be addressed anytime a standardised Z-score is computed. We find that the severity of autocorrelation is highly dependent on spatial characteristics of brain regions, such as the size of regions of interest and the spatial location of those regions. We further show that the available EDF estimators make restrictive assumptions that are not supported by the data, resulting in biased rsFC inferences that lead to distorted topological descriptions of the connectome on the individual level. We propose a practical "xDF" method that accounts not only for distinct autocorrelation in each time series, but instantaneous and lagged cross-correlation. We find the xDF correction varies substantially over node pairs, indicating the limitations of global EDF corrections used previously. In addition to extensive synthetic and real data validations, we investigate the impact of this correction on rsFC measures in data from the Young Adult Human Connectome Project, showing that accounting for autocorrelation dramatically changes fundamental graph theoretical measures relative to no correction.

Common-injection-induced isolated desynchronization in delay-coupled VCSELs networks with variable-polarization optical feedback

We have investigated the cluster isolated desynchronization, a symmetry-breaking state, in the delay-coupled vertical-cavity surface-emitting lasers (VCSELs) networks subject to variable-polarization optical feedback (VPOF). It is shown that, in the VPOF-VCSELs networks, the elusive isolated desynchronization phenomenon could emerge out of the cluster synchronization by the common-signal injection approach from an additional auxiliary VCSEL. The influences of parameters in VPOF-VCSELs networks on the existence and stability of isolated desynchronization are systematically investigated. Moreover, the generality of the proposed scheme is validated in the VPOF-VCSELs network with real-world network topology (Nepal power grid network). Our results offer a new insight to manage the synchronization patterns of a VCSELs network.