Dynamics of globally coupled oscillators: Progress and perspectives

Synchronization transitions in phase oscillator populations with partial adaptive coupling

The adaptation underlying many realistic processes plays a pivotal role in shaping the collective dynamics of diverse systems. Here, we untangle the generic conditions for synchronization transitions in a system of coupled phase oscillators incorporating the adaptive scheme encoded by the feedback between the coupling and the order parameter via a power-law function with different weights. We mathematically argue that, in the subcritical and supercritical correlation scenarios, there exists no critical adaptive fraction for synchronization transitions converting from the first (second)-order to the second (first)-order. In contrast to the synchronization transitions previously deemed, the explosive and continuous phase transitions take place in the corresponding regions as long as the adaptive fraction is nonzero, respectively. Nevertheless, we uncover that, at the critical correlation, the routes toward synchronization depend crucially on the relative adaptive weights. In particular, we unveil that the emergence of a range of interrelated scaling behaviors of the order parameter near criticality, manifesting the subcritical and supercritical bifurcations, are responsible for various observed phase transitions. Our work, thus, provides profound insights for understanding the dynamical nature of phase transitions, and for better controlling and manipulating synchronization transitions in networked systems with adaptation.

Nature of the Volcano Transition in the Fully Disordered Kuramoto Model

Randomly coupled phase oscillators may synchronize into disordered patterns of collective motion. We analyze this transition in a large, fully connected Kuramoto model with symmetric but otherwise independent random interactions. Using the dynamical cavity method, we reduce the dynamics to a stochastic single-oscillator problem with self-consistent correlation and response functions that we study analytically and numerically. We clarify the nature of the volcano transition and elucidate its relation to the existence of an oscillator glass phase.

Synchronization dynamics of phase oscillators on power grid models

We investigate topological and spectral properties of models of European and US-American power grids and of paradigmatic network models as well as their implications for the synchronization dynamics of phase oscillators with heterogeneous natural frequencies. We employ the complex-valued order parameter-a widely used indicator for phase ordering-to assess the synchronization dynamics and observe the order parameter to exhibit either constant or periodic or non-periodic, possibly chaotic temporal evolutions for a given coupling strength but depending on initial conditions and the systems' disorder. Interestingly, both topological and spectral characteristics of the power grids point to a diminished capability of these networks to support a temporarily stable synchronization dynamics. We find non-trivial commonalities between the synchronization dynamics of oscillators on seemingly opposing topologies.

A robust balancing mechanism for spiking neural networks

Dynamical balance of excitation and inhibition is usually invoked to explain the irregular low firing activity observed in the cortex. We propose a robust nonlinear balancing mechanism for a random network of spiking neurons, which works also in the absence of strong external currents. Biologically, the mechanism exploits the plasticity of excitatory-excitatory synapses induced by short-term depression. Mathematically, the nonlinear response of the synaptic activity is the key ingredient responsible for the emergence of a stable balanced regime. Our claim is supported by a simple self-consistent analysis accompanied by extensive simulations performed for increasing network sizes. The observed regime is essentially fluctuation driven and characterized by highly irregular spiking dynamics of all neurons.

Integrability of a Globally Coupled Complex Riccati Array: Quadratic Integrate-and-Fire Neurons, Phase Oscillators, and All in Between

We present an exact dimensionality reduction for dynamics of an arbitrary array of globally coupled complex-valued Riccati equations. It generalizes the Watanabe-Strogatz theory [Integrability of a globally coupled oscillator array, Phys. Rev. Lett. 70, 2391 (1993).PRLTAO0031-900710.1103/PhysRevLett.70.2391] for sinusoidally coupled phase oscillators and seamlessly includes quadratic integrate-and-fire neurons as the real-valued special case. This simple formulation reshapes our understanding of a broad class of coupled systems-including a particular class of phase-amplitude oscillators-which newly fall under the category of integrable systems. Precise and rigorous analysis of complex Riccati arrays is now within reach, paving a way to a deeper understanding of emergent behavior of collective dynamics in coupled systems.

Deterministic correlations enhance synchronization in oscillator populations with heterogeneous coupling

Synchronization is a critical phenomenon that displays a pivotal role in a wealth of dynamical processes ranging from natural to artificial systems. Here, we untangle the synchronization optimization in a system of globally coupled phase oscillators incorporating heterogeneous interactions encoded by the deterministic-random coupling. We uncover that, within the given restriction, the added deterministic correlations can profoundly enhance the synchronizability in comparison with the uncorrelated scenario. The critical points manifesting the onset of synchronization and desynchronization transitions, as well as the level of phase coherence, are significantly shaped by the increment of deterministic correlations. In particular, we provide an analytical treatment to properly ground the mechanism underlying synchronization enhancement and substantiate that the analytical predictions are in fair agreement with the numerical simulations. This study is a step forward in highlighting the importance of heterogeneous coupling among dynamical agents, which provides insights for control strategies of synchronization in complex systems.

Chimeras in globally coupled oscillators: A review

The surprising phenomenon of chimera in an ensemble of identical oscillators is no more strange behavior of network dynamics and reality. By this time, this symmetry breaking self-organized collective dynamics has been established in many networks, a ring of non-locally coupled oscillators, globally coupled networks, a three-dimensional network, and multi-layer networks. A variety of coupling and dynamical models in addition to the phase oscillators has been used for a successful observation of chimera patterns. Experimental verification has also been done using metronomes, pendula, chemical, and opto-electronic systems. The phenomenon has also been shown to appear in small networks, and hence, it is not size-dependent. We present here a brief review of the origin of chimera patterns restricting our discussions to networks of globally coupled identical oscillators only. The history of chimeras in globally coupled oscillators is older than what has been reported in nonlocally coupled phase oscillators much later. We elaborate the story of the origin of chimeras in globally coupled oscillators in a chronological order, within our limitations, and with brief descriptions of the significant contributions, including our personal experiences. We first introduce chimeras in non-locally coupled and other network configurations, in general, and then discuss about globally coupled networks in more detail.

Multi-modal and multi-model interrogation of large-scale functional brain networks

Existing whole-brain models are generally tailored to the modelling of a particular data modality (e.g., fMRI or MEG/EEG). We propose that despite the differing aspects of neural activity each modality captures, they originate from shared network dynamics. Building on the universal principles of self-organising delay-coupled nonlinear systems, we aim to link distinct features of brain activity - captured across modalities - to the dynamics unfolding on a macroscopic structural connectome. To jointly predict connectivity, spatiotemporal and transient features of distinct signal modalities, we consider two large-scale models - the Stuart Landau and Wilson and Cowan models - which generate short-lived 40 Hz oscillations with varying levels of realism. To this end, we measure features of functional connectivity and metastable oscillatory modes (MOMs) in fMRI and MEG signals - and compare them against simulated data. We show that both models can represent MEG functional connectivity (FC), functional connectivity dynamics (FCD) and generate MOMs to a comparable degree. This is achieved by adjusting the global coupling and mean conduction time delay and, in the WC model, through the inclusion of balance between excitation and inhibition. For both models, the omission of delays dramatically decreased the performance. For fMRI, the SL model performed worse for FCD and MOMs, highlighting the importance of balanced dynamics for the emergence of spatiotemporal and transient patterns of ultra-slow dynamics. Notably, optimal working points varied across modalities and no model was able to achieve a correlation with empirical FC higher than 0.4 across modalities for the same set of parameters. Nonetheless, both displayed the emergence of FC patterns that extended beyond the constraints of the anatomical structure. Finally, we show that both models can generate MOMs with empirical-like properties such as size (number of brain regions engaging in a mode) and duration (continuous time interval during which a mode appears). Our results demonstrate the emergence of static and dynamic properties of neural activity at different timescales from networks of delay-coupled oscillators at 40 Hz. Given the higher dependence of simulated FC on the underlying structural connectivity, we suggest that mesoscale heterogeneities in neural circuitry may be critical for the emergence of parallel cross-modal functional networks and should be accounted for in future modelling endeavours.

A physics-based model of swarming jellyfish

We propose a model for the structure formation of jellyfish swimming based on active Brownian particles. We address the phenomena of counter-current swimming, avoidance of turbulent flow regions and foraging. We motivate corresponding mechanisms from observations of jellyfish swarming reported in the literature and incorporate them into the generic modelling framework. The model characteristics is tested in three paradigmatic flow environments.

Transient Phase Clusters in a Two-Population Network of Kuramoto Oscillators with Heterogeneous Adaptive Interaction

Adaptive interactions are an important property of many real-word network systems. A feature of such networks is the change in their connectivity depending on the current states of the interacting elements. In this work, we study the question of how the heterogeneous character of adaptive couplings influences the emergence of new scenarios in the collective behavior of networks. Within the framework of a two-population network of coupled phase oscillators, we analyze the role of various factors of heterogeneous interaction, such as the rules of coupling adaptation and the rate of their change in the formation of various types of coherent behavior of the network. We show that various schemes of heterogeneous adaptation lead to the formation of transient phase clusters of various types.

Proposing magnetoimpedance effect for neuromorphic computing

Oscillation of physical parameters in materials can result in a peak signal in the frequency spectrum of the voltage measured from the materials. This spectrum and its amplitude/frequency tunability, through the application of bias voltage or current, can be used to perform neuron-like cognitive tasks. Magnetic materials, after achieving broad distribution for data storage applications in classical Von Neumann computer architectures, are under intense investigation for their neuromorphic computing capabilities. A recent successful demonstration regards magnetisation oscillation in magnetic thin films by spin transfer or spin orbit torques accompanied by magnetoresistance (MR) effect that can give a voltage peak in the frequency spectrum of voltage with bias current dependence of both peak frequency and amplitude. Here we use classical magnetoimpedance (MI) effect in a magnetic wire to produce such a peak and manipulate its frequency and amplitude by means of the bias voltage. We applied a noise signal to a magnetic wire with high magnetic permeability and owing to the frequency dependence of the magnetic permeability we got frequency dependent impedance with a peak at the maximum permeability. Frequency dependence of the MI effect results in different changes in the voltage amplitude at each frequency when a bias voltage is applied and therefore a shift in the peak position and amplitude can be obtained. The presented method and material provide optimal features in structural simplicity, low-frequency operation (tens of MHz-order) and high robustness at different environmental conditions. Our universal approach can be applied to any system with frequency dependent bias responses.

Exact finite-dimensional description for networks of globally coupled spiking neurons

We consider large networks of globally coupled spiking neurons and derive an exact low-dimensional description of their collective dynamics in the thermodynamic limit. Individual neurons are described by the Ermentrout-Kopell canonical model that can be excitable or tonically spiking and interact with other neurons via pulses. Utilizing the equivalence of the quadratic integrate-and-fire and the theta-neuron formulations, we first derive the dynamical equations in terms of the Kuramoto-Daido order parameters (Fourier modes of the phase distribution) and relate them to two biophysically relevant macroscopic observables, the firing rate and the mean voltage. For neurons driven by Cauchy white noise or for Cauchy-Lorentz distributed input currents, we adapt the results by Cestnik and Pikovsky [Chaos 32, 113126 (2022)1054-150010.1063/5.0106171] and show that for arbitrary initial conditions the collective dynamics reduces to six dimensions. We also prove that in this case the dynamics asymptotically converges to a two-dimensional invariant manifold first discovered by Ott and Antonsen. For identical, noise-free neurons, the dynamics reduces to three dimensions, becoming equivalent to the Watanabe-Strogatz description. We illustrate the exact six-dimensional dynamics outside the invariant manifold by calculating nontrivial basins of different asymptotic regimes in a bistable situation.

Relaxation dynamics of phase oscillators with generic heterogeneous coupling

The coupled phase oscillator model serves as a paradigm that has been successfully used to shed light on the collective dynamics occurring in large ensembles of interacting units. It was widely known that the system experiences a continuous (second-order) phase transition to synchronization by gradually increasing the homogeneous coupling among the oscillators. As the interest in exploring synchronized dynamics continues to grow, the heterogeneous patterns between phase oscillators have received ample attention during the past years. Here, we consider a variant of the Kuramoto model with quenched disorder in their natural frequencies and coupling. Correlating these two types of heterogeneity via a generic weighted function, we systematically investigate the impacts of the heterogeneous strategies, the correlation function, and the natural frequency distribution on the emergent dynamics. Importantly, we develop an analytical treatment for capturing the essential dynamical properties of the equilibrium states. In particular, we uncover that the critical threshold corresponding to the onset of synchronization is unaffected by the location of the inhomogeneity, which, however, does depend crucially on the value of the correlation function at its center. Furthermore, we reveal that the relaxation dynamics of the incoherent state featuring the responses to external perturbations is significantly shaped by all the considered effects, thereby leading to various decaying mechanisms of the order parameters in the subcritical region. Moreover, we untangle that synchronization is facilitated by the out-coupling strategy in the supercritical region. Our study is a step forward in highlighting the potential importance of the inhomogeneous patterns involved in the complex systems, and could thus provide theoretical insights for profoundly understanding the generic statistical mechanical properties of the steady states toward synchronization.

Introduction to Focus Issue: Dynamics of oscillator populations

Even after about 50 years of intensive research, the dynamics of oscillator populations remain one of the most popular topics in nonlinear science. This Focus Issue brings together studies on such diverse aspects of the problem as low-dimensional description, effects of noise and disorder on synchronization transition, control of synchrony, the emergence of chimera states and chaotic regimes, stability of power grids, etc.

Coevolutionary Dynamics with Global Fields

We investigate the effects of external and autonomous global interaction fields on an adaptive network of social agents with an opinion formation dynamics based on a simple imitation rule. We study the competition between global fields and adaptive rewiring on the space of parameters of the system. The model represents an adaptive society subject to global mass media such as a directed opinion influence or feedback of endogenous cultural trends. We show that, in both situations, global mass media contribute to consensus and to prevent the fragmentation of the social network induced by the coevolutionary dynamics. We present a discussion of these results in the context of dynamical systems and opinion formation dynamics.

Symmetry breaking yields chimeras in two small populations of Kuramoto-type oscillators

Despite their simplicity, networks of coupled phase oscillators can give rise to intriguing collective dynamical phenomena. However, the symmetries of globally and identically coupled identical units do not allow solutions where distinct oscillators are frequency-unlocked-a necessary condition for the emergence of chimeras. Thus, forced symmetry breaking is necessary to observe chimera-type solutions. Here, we consider the bifurcations that arise when full permutational symmetry is broken for the network to consist of coupled populations. We consider the smallest possible network composed of four phase oscillators and elucidate the phase space structure, (partial) integrability for some parameter values, and how the bifurcations away from full symmetry lead to frequency-unlocked weak chimera solutions. Since such solutions wind around a torus they must arise in a global bifurcation scenario. Moreover, periodic weak chimeras undergo a period-doubling cascade leading to chaos. The resulting chaotic dynamics with distinct frequencies do not rely on amplitude variation and arise in the smallest networks that support chaos.

A global synchronization theorem for oscillators on a random graph

Consider n identical Kuramoto oscillators on a random graph. Specifically, consider Erdős-Rényi random graphs in which any two oscillators are bidirectionally coupled with unit strength, independently and at random, with probability 0 ≤ p ≤ 1. We say that a network is globally synchronizing if the oscillators converge to the all-in-phase synchronous state for almost all initial conditions. Is there a critical threshold for p above which global synchrony is extremely likely but below which it is extremely rare? It is suspected that a critical threshold exists and is close to the so-called connectivity threshold, namely, p ∼ log ⁡ ( n ) / n for n ≫ 1. Ling, Xu, and Bandeira made the first progress toward proving a result in this direction: they showed that if p ≫ log ⁡ ( n ) / ⁿ, then Erdős-Rényi networks of Kuramoto oscillators are globally synchronizing with high probability as n → ∞. Here, we improve that result by showing that p ≫ ^log ⁡ ( n ) / n suffices. Our estimates are explicit: for example, we can say that there is more than a 99.9996 % chance that a random network with n = ¹⁰ and p > 0.011 17 is globally synchronizing.

Intensity-Varied Closed-Loop Noise Stimulation for Oscillation Suppression in the Parkinsonian State

This work explores the effectiveness of the intensity-varied closed-loop noise stimulation on the oscillation suppression in the Parkinsonian state. Deep brain stimulation (DBS) is the standard therapy for Parkinson's disease (PD), but its effects need to be improved. The noise stimulation has compelling results in alleviating the PD state. However, in the open-loop control scheme, the noise stimulation parameters cannot be self-adjusted to adapt to the amplitude of the synchronized neuronal activities in real time. Thus, based on the delayed-feedback control algorithm, an intensity-varied closed-loop noise stimulation strategy is proposed. Based on a computational model of the basal ganglia (BG) that can present the intrinsic properties of the BG neurons and their interactions with the thalamic neurons, the proposed stimulation strategy is tested. Simulation results show that the noise stimulation suppresses the pathological beta (12-35 Hz) oscillations without any new rhythms in other bands compared with traditional high-frequency DBS. The intensity-varied closed-loop noise stimulation has a more profound role in removing the pathological beta oscillations and improving the thalamic reliability than open-loop noise stimulation, especially for different PD states. And the closed-loop noise stimulation enlarges the parameter space of the delayed-feedback control algorithm due to the randomness of noise signals. We also provide a theoretical analysis of the effective parameter domain of the delayed-feedback control algorithm by simplifying the BG model to an oscillator model. This exploration may guide a new approach to treating PD by optimizing the noise-induced improvement of the BG dysfunction.

On the origin of chaotic attractors with two zero Lyapunov exponents in a system of five biharmonically coupled phase oscillators

We study chaotic dynamics in a system of four differential equations describing the interaction of five identical phase oscillators coupled via biharmonic function. We show that this system exhibits strange spiral attractors (Shilnikov attractors) with two zero (indistinguishable from zero in numerics) Lyapunov exponents in a wide region of the parameter space. We explain this phenomenon by means of bifurcation analysis of a three-dimensional Poincaré map for the system under consideration. We show that chaotic dynamics develop here near a codimension three bifurcation, when a periodic orbit (fixed point of the Poincaré map) has the triplet of multipliers ( 1 , 1 , 1 ). As it is known, the flow normal form for such bifurcation is the well-known three-dimensional Arneodó-Coullet-Spiegel-Tresser (ACST) system, which exhibits spiral attractors. According to this, we conclude that the additional zero Lyapunov exponent for orbits in the observed attractors appears due to the fact that the corresponding three-dimensional Poincaré map is very close to the time-shift map of the ACST-system.

Synchronization in the Kuramoto model in presence of stochastic resetting

What happens when the paradigmatic Kuramoto model involving interacting oscillators of distributed natural frequencies and showing spontaneous collective synchronization in the stationary state is subject to random and repeated interruptions of its dynamics with a reset to the initial condition? While resetting to a synchronized state, it may happen between two successive resets that the system desynchronizes, which depends on the duration of the random time interval between the two resets. Here, we unveil how such a protocol of stochastic resetting dramatically modifies the phase diagram of the bare model, allowing, in particular, for the emergence of a synchronized phase even in parameter regimes for which the bare model does not support such a phase. Our results are based on an exact analysis invoking the celebrated Ott-Antonsen ansatz for the case of the Lorentzian distribution of natural frequencies and numerical results for Gaussian frequency distribution. Our work provides a simple protocol to induce global synchrony in the system through stochastic resetting.

Partial locking in phase-oscillator populations with heterogenous coupling

We consider a variant of the mean-field model of coupled phase oscillators with uniform distribution of natural frequencies. By establishing correlations between the quenched disorder of intrinsic frequencies and coupling strength with both in- and out-coupling heterogeneities, we reveal a generic criterion for the onset of partial locking that takes place in a domain with the coexistence of phase-locked oscillators and drifters. The critical points manifesting the instability of the stationary states are obtained analytically. In particular, the bifurcation mechanism of the equilibrium states is uncovered by the use of frequency-dependent version of the Ott-Antonsen reduction consistently with the analysis based on the self-consistent approach. We demonstrate that both the manner of coupling heterogeneity and correlation exponent have influence on the emergent patterns of partial locking. Our research could find applicability in better understanding the phase transitions and related collective phenomena involving synchronization control in networked systems.

Mean-field analysis of Stuart-Landau oscillator networks with symmetric coupling and dynamical noise

This paper presents analyses of networks composed of homogeneous Stuart-Landau oscillators with symmetric linear coupling and dynamical Gaussian noise. With a simple mean-field approximation, the original system is transformed into a surrogate system that describes uncorrelated oscillation/fluctuation modes of the original system. The steady-state probability distribution for these modes is described using an exponential family, and the dynamics of the system are mainly determined by the eigenvalue spectrum of the coupling matrix and the noise level. The variances of the modes can be expressed as functions of the eigenvalues and noise level, yielding the relation between the covariance matrix and the coupling matrix of the oscillators. With decreasing noise, the leading mode changes from fluctuation to oscillation, generating apparent synchrony of the coupled oscillators, and the condition for such a transition is derived. Finally, the approximate analyses are examined via numerical simulation of the oscillator networks with weak coupling to verify the utility of the approximation in outlining the basic properties of the considered coupled oscillator networks. These results are potentially useful for the modeling and analysis of indirectly measured data of neurodynamics, e.g., via functional magnetic resonance imaging and electroencephalography, as a counterpart of the frequently used Ising model.

A self-consistent analytical theory for rotator networks under stochastic forcing: Effects of intrinsic noise and common input

Despite the incredible complexity of our brains' neural networks, theoretical descriptions of neural dynamics have led to profound insights into possible network states and dynamics. It remains challenging to develop theories that apply to spiking networks and thus allow one to characterize the dynamic properties of biologically more realistic networks. Here, we build on recent work by van Meegen and Lindner who have shown that "rotator networks," while considerably simpler than real spiking networks and, therefore, more amenable to mathematical analysis, still allow one to capture dynamical properties of networks of spiking neurons. This framework can be easily extended to the case where individual units receive uncorrelated stochastic input, which can be interpreted as intrinsic noise. However, the assumptions of the theory do not apply anymore when the input received by the single rotators is strongly correlated among units. As we show, in this case, the network fluctuations become significantly non-Gaussian, which calls for reworking of the theory. Using a cumulant expansion, we develop a self-consistent analytical theory that accounts for the observed non-Gaussian statistics. Our theory provides a starting point for further studies of more general network setups and information transmission properties of these networks.

Efficient moment-based approach to the simulation of infinitely many heterogeneous phase oscillators

The dynamics of ensembles of phase oscillators are usually described considering their infinite-size limit. In practice, however, this limit is fully accessible only if the Ott-Antonsen theory can be applied, and the heterogeneity is distributed following a rational function. In this work, we demonstrate the usefulness of a moment-based scheme to reproduce the dynamics of infinitely many oscillators. Our analysis is particularized for Gaussian heterogeneities, leading to a Fourier-Hermite decomposition of the oscillator density. The Fourier-Hermite moments obey a set of hierarchical ordinary differential equations. As a preliminary experiment, the effects of truncating the moment system and implementing different closures are tested in the analytically solvable Kuramoto model. The moment-based approach proves to be much more efficient than the direct simulation of a large oscillator ensemble. The convenience of the moment-based approach is exploited in two illustrative examples: (i) the Kuramoto model with bimodal frequency distribution, and (ii) the "enlarged Kuramoto model" (endowed with nonpairwise interactions). In both systems, we obtain new results inaccessible through direct numerical integration of populations.

Mean-field equations for neural populations with q-Gaussian heterogeneities

Describing the collective dynamics of large neural populations using low-dimensional models for averaged variables has long been an attractive task in theoretical neuroscience. Recently developed reduction methods make it possible to derive such models directly from the microscopic dynamics of individual neurons. To simplify the reduction, the Cauchy distribution is usually assumed for heterogeneous network parameters. Here we extend the reduction method for a wider class of heterogeneities defined by the q-Gaussian distribution. The shape of this distribution depends on the Tsallis index q and gradually changes from the Cauchy distribution to the normal Gaussian distribution as this index changes. We derive the mean-field equations for an inhibitory network of quadratic integrate-and-fire neurons with a q-Gaussian-distributed excitability parameter. It is shown that the dynamic modes of the network significantly depend on the form of the distribution determined by the Tsallis index. The results obtained from the mean-field equations are confirmed by numerical simulation of the microscopic model.

Sparsity-driven synchronization in oscillator networks

The emergence of synchronized behavior is a direct consequence of networking dynamical systems. Naturally, strict instances of this phenomenon, such as the states of complete synchronization, are favored or even ensured in networks with a high density of connections. Conversely, in sparse networks, the system state-space is often shared by a variety of coexistent solutions. Consequently, the convergence to complete synchronized states is far from being certain. In this scenario, we report the surprising phenomenon in which completely synchronized states are made the sole attractor of sparse networks by removing network links, the sparsity-driven synchronization. This phenomenon is observed numerically for nonlocally coupled Kuramoto networks and verified analytically for locally coupled ones. In addition, we unravel the bifurcation scenario underlying the network transition to completely synchronized behavior. Furthermore, we present a simple procedure, based on the bifurcations in the thermodynamic limit, that determines the minimum number of links to be removed in order to ensure complete synchronization. Finally, we propose an application of the reported phenomenon as a control scheme to drive complete synchronization in high connectivity networks.

Coherent oscillations in balanced neural networks driven by endogenous fluctuations

We present a detailed analysis of the dynamical regimes observed in a balanced network of identical quadratic integrate-and-fire neurons with sparse connectivity for homogeneous and heterogeneous in-degree distributions. Depending on the parameter values, either an asynchronous regime or periodic oscillations spontaneously emerge. Numerical simulations are compared with a mean-field model based on a self-consistent Fokker-Planck equation (FPE). The FPE reproduces quite well the asynchronous dynamics in the homogeneous case by either assuming a Poissonian or renewal distribution for the incoming spike trains. An exact self-consistent solution for the mean firing rate obtained in the limit of infinite in-degree allows identifying balanced regimes that can be either mean- or fluctuation-driven. A low-dimensional reduction of the FPE in terms of circular cumulants is also considered. Two cumulants suffice to reproduce the transition scenario observed in the network. The emergence of periodic collective oscillations is well captured both in the homogeneous and heterogeneous setups by the mean-field models upon tuning either the connectivity or the input DC current. In the heterogeneous situation, we analyze also the role of structural heterogeneity.

Kuramoto model for populations of quadratic integrate-and-fire neurons with chemical and electrical coupling

We derive the Kuramoto model (KM) corresponding to a population of weakly coupled, nearly identical quadratic integrate-and-fire (QIF) neurons with both electrical and chemical coupling. The ratio of chemical to electrical coupling determines the phase lag of the characteristic sine coupling function of the KM and critically determines the synchronization properties of the network. We apply our results to uncover the presence of chimera states in two coupled populations of identical QIF neurons. We find that the presence of both electrical and chemical coupling is a necessary condition for chimera states to exist. Finally, we numerically demonstrate that chimera states gradually disappear as coupling strengths cease to be weak.

Collective behavior of swarmalators on a ring

We study the collective behavior of swarmalators, generalizations of phase oscillators that both sync and swarm, confined to move on a one-dimensional (1D) ring. This simple model captures the essence of movement in two or three dimensions, but has the benefit of being solvable: most of the collective states and their bifurcations can be specified exactly. The model also captures the behavior of real-world swarmalators which swarm in quasi-1D rings such as bordertaxic vinegar eels and sperm.

Collective dynamics of phase oscillator populations with three-body interactions

Many-body interactions between dynamical agents have caught particular attention in recent works that found wide applications in physics, neuroscience, and sociology. In this paper we investigate such higher order (nonadditive) interactions on collective dynamics in a system of globally coupled heterogeneous phase oscillators. We show that the three-body interactions encoded microscopically in nonlinear couplings give rise to added dynamic phenomena occurring beyond the pairwise interactions. The system in general displays an abrupt desynchronization transition characterized by irreversible explosive synchronization via an infinite hysteresis loop. More importantly, we give a mathematical argument that such an abrupt dynamic pattern is a universally expected effect. Furthermore, the origin of this abrupt transition is uncovered by performing a rigorous stability analysis of the equilibrium states, as well as by providing a detailed description of the spectrum structure of linearization around the steady states. Our work reveals a self-organized phenomenon that is responsible for the rapid switching to synchronization in diverse complex systems exhibiting critical transitions with nonpairwise interactions.

Contrarians Synchronize beyond the Limit of Pairwise Interactions

We give evidence that a population of pure contrarian globally coupled D-dimensional Kuramoto oscillators reaches a collective synchronous state when the interplay between the units goes beyond the limit of pairwise interactions. Namely, we will show that the presence of higher-order interactions may induce the appearance of a coherent state even when the oscillators are coupled negatively to the mean field. An exact solution for the description of the microscopic dynamics for forward and backward transitions is provided, which entails imperfect symmetry breaking of the population into a frequency-locked state featuring two clusters of different instantaneous phases. Our results contribute to a better understanding of the powerful potential of group interactions entailing multidimensional choices and novel dynamical states in many circumstances, such as in social systems.

Delayed Feedback-Based Suppression of Pathological Oscillations in a Neural Mass Model

Suppression of excessively synchronous beta frequency (12-35 Hz) oscillatory activity in the basal ganglia is believed to correlate with the alleviation of hypokinetic motor symptoms of the Parkinson's disease. Delayed feedback is an effective strategy to interrupt the synchronization and has been used in the design of closed-loop neuromodulation methods computationally. Although tremendous efforts in this are being made by optimizing delayed feedback algorithm and stimulation waveforms, there are still remaining problems in the selection of effective parameters in the delayed feedback control schemes. In most delayed feedback neuromodulation strategies, the stimulation signal is obtained from the local field potential (LFP) of the excitatory subthalamic nucleus (STN) neurons and then is administered back to STN itself only. The inhibitory external globus pallidus (GPe) nucleus in the excitatory-inhibitory STN-GPe reciprocal network has not been involved in the design of the delayed feedback control strategies. Thus, considering the role of GPe, this paper proposes three schemes involving GPe in the design of the delayed feedback strategies and compared their effectiveness to the traditional paradigm using STN only. Based on a neural mass model of STN-GPe network having capability of simulating the LFP directly, the proposed stimulation strategies are tested and compared. Our simulation results show that the four types of delayed feedback control schemes are all effective, even if with a simple linear delayed feedback algorithm. But the three new control strategies we propose here further improve the control performance by enlarging the oscillatory suppression space and reducing the energy expenditure, suggesting that they may be more effective in applications. This paper may guide a new approach to optimize the closed-loop deep brain stimulation treatment to alleviate the Parkinsonian state by retargeting the measurement and stimulation nucleus.

The Kuramoto model on a sphere: Explaining its low-dimensional dynamics with group theory and hyperbolic geometry

We study a system of N identical interacting particles moving on the unit sphere in d-dimensional space. The particles are self-propelled and coupled all to all, and their motion is heavily overdamped. For d=2, the system reduces to the classic Kuramoto model of coupled oscillators; for d=3, it has been proposed to describe the orientation dynamics of swarms of drones or other entities moving about in three-dimensional space. Here, we use group theory to explain the recent discovery that the model shows low-dimensional dynamics for all N≥3 and to clarify why it admits the analog of the Ott-Antonsen ansatz in the continuum limit N→∞. The underlying reason is that the system is intimately connected to the natural hyperbolic geometry on the unit ball B^d. In this geometry, the isometries form a Lie group consisting of higher-dimensional generalizations of the Möbius transformations used in complex analysis. Once these connections are realized, the reduced dynamics and the generalized Ott-Antonsen ansatz follow immediately. This framework also reveals the seamless connection between the finite and infinite- N cases. Finally, we show that special forms of coupling yield gradient dynamics with respect to the hyperbolic metric and use that fact to obtain global stability results about convergence to the synchronized state.

Sufficiently dense Kuramoto networks are globally synchronizing

Consider any network of n identical Kuramoto oscillators in which each oscillator is coupled bidirectionally with unit strength to at least μ(n-1) other oscillators. There is a critical value of the connectivity, μ_c, such that whenever μ>μ_c, the system is guaranteed to converge to the all-in-phase synchronous state for almost all initial conditions, but when μ<μ_c, there are networks with other stable states. The precise value of the critical connectivity remains unknown, but it has been conjectured to be μ_c=0.75. In 2020, Lu and Steinerberger proved that μ_c≤0.7889, and Yoneda, Tatsukawa, and Teramae proved in 2021 that μ_c>0.6838. This paper proves that μ_c≤0.75 and explain why this is the best upper bound that one can obtain by a purely linear stability analysis.

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.

Transient circulant clusters in two-population network of Kuramoto oscillators with different rules of coupling adaptation

We considered a network consisting of two populations of phase oscillators, the interaction of which is determined by different rules for the coupling adaptation. The introduction of various adaptation rules leads to the suppression of splay states and the emergence of each population complex non-stationary behavior called transient circulant clusters. In such states, each population contains a pair of anti-phase clusters whose size and composition slowly change over time as a result of successive transitions of oscillators between clusters. We show that an increase in the mismatch of the adaptation rules makes it possible to stop the process of rearrangement of clusters in one or both populations of the network. Transitions to such modes are always preceded by the appearance of solitary states in one of the populations.

Discordant synchronization patterns on directed networks of identical phase oscillators with attractive and repulsive couplings

We study the collective dynamics of identical phase oscillators on globally coupled networks whose interactions are asymmetric and mediated by positive and negative couplings. We split the set of oscillators into two interconnected subpopulations. In this setup, oscillators belonging to the same group interact via symmetric couplings while the interaction between subpopulations occurs in an asymmetric fashion. By employing the dimensional reduction scheme of the Ott-Antonsen (OA) theory, we verify the existence of traveling wave and π-states, in addition to the classical fully synchronized and incoherent states. Bistability between all collective states is reported. Analytical results are generally in excellent agreement with simulations; for some parameters and initial conditions, however, we numerically detect chimera-like states which are not captured by the OA theory.

Modeling synchronization in globally coupled oscillatory systems using model order reduction

We construct reduced order models for two classes of globally coupled multi-component oscillatory systems, selected as prototype models that exhibit synchronization. These are the Kuramoto model, considered both in its original formulation and with a suitable change of coordinates, and a model for the circadian clock. The systems of interest possess strong reduction properties, as their dynamics can be efficiently described with a low-dimensional set of coordinates. Specifically, the solution and selected quantities of interest are well approximated at the reduced level, and the reduced models recover the expected transition to synchronized states as the coupling strengths vary. Assuming that the interactions depend only on the averages of the system variables, the surrogate models ensure a significant computational speedup for large systems.

Stability and bifurcation of collective dynamics in phase oscillator populations with general coupling

The Kuramoto model serves as an illustrative paradigm for studying the synchronization transitions and collective behaviors in large ensembles of coupled dynamical units. In this paper, we present a general framework for analytically capturing the stability and bifurcation of the collective dynamics in oscillator populations by extending the global coupling to depend on an arbitrary function of the Kuramoto order parameter. In this generalized Kuramoto model with rotation and reflection symmetry, we show that all steady states characterizing the long-term macroscopic dynamics can be expressed in a universal profile given by the frequency-dependent version of the Ott-Antonsen reduction, and the introduced empirical stability criterion for each steady state degenerates to a remarkably simple expression described by the self-consistent equation [Iatsenko et al., Phys. Rev. Lett. 110, 064101 (2013)PRLTAO0031-900710.1103/PhysRevLett.110.064101]. Here, we provide a detailed description of the spectrum structure in the complex plane by performing a rigorous stability analysis of various steady states in the reduced system. More importantly, we uncover that the empirical stability criterion for each steady state involved in the system is completely equivalent to its linear stability condition that is determined by the nontrivial eigenvalues (discrete spectrum) of the linearization. Our study provides a new and widely applicable approach for exploring the stability properties of collective synchronization, which we believe improves the understanding of the underlying mechanisms of phase transitions and bifurcations in coupled dynamical networks.

Explosive synchronization in interlayer phase-shifted Kuramoto oscillators on multiplex networks

Different methods have been proposed in the past few years to incite explosive synchronization (ES) in Kuramoto phase oscillators. In this work, we show that the introduction of a phase shift α in interlayer coupling terms of a two-layer multiplex network of Kuramoto oscillators can also instigate ES in the layers. As α→π/2, ES emerges along with hysteresis. The width of hysteresis depends on the phase shift α, interlayer coupling strength, and natural frequency mismatch between mirror nodes. A mean-field analysis is performed to justify the numerical results. Similar to earlier works, the suppression of synchronization is accountable for the occurrence of ES. The robustness of ES against changes in network topology and natural frequency distribution is tested. Finally, taking a suggestion from the synchronized state of the multiplex networks, we extend the results to classical single networks where some specific links are assigned phase-shifted interactions.

Critical exponents in coupled phase-oscillator models on small-world networks

A coupled phase-oscillator model consists of phase oscillators, each of which has the natural frequency obeying a probability distribution and couples with other oscillators through a given periodic coupling function. This type of model is widely studied since it describes the synchronization transition, which emerges between the nonsynchronized state and partially synchronized states. The synchronization transition is characterized by several critical exponents, and we focus on the critical exponent defined by coupling strength dependence of the order parameter for revealing universality classes. In a typical interaction represented by the perfect graph, an infinite number of universality classes is yielded by dependency on the natural frequency distribution and the coupling function. Since the synchronization transition is also observed in a model on a small-world network, whose number of links is proportional to the number of oscillators, a natural question is whether the infinite number of universality classes remains in small-world networks irrespective of the order of links. Our numerical results suggest that the number of universality classes is reduced to one and the critical exponent is shared in the considered models having coupling functions up to second harmonics with unimodal and symmetric natural frequency distributions.

Quasi phase reduction of all-to-all strongly coupled λ-ω oscillators near incoherent states

The dynamics of an ensemble of N weakly coupled limit-cycle oscillators can be captured by their N phases using standard phase reduction techniques. However, it is a phenomenological fact that all-to-all strongly coupled limit-cycle oscillators may behave as "quasiphase oscillators," evidencing the need of novel reduction strategies. We introduce, here, quasi phase reduction (QPR), a scheme suited for identical oscillators with polar symmetry (λ-ω systems). By applying QPR, we achieve a reduction to N+2 degrees of freedom: N phase oscillators interacting through one independent complex variable. This "quasi phase model" is asymptotically valid in the neighborhood of incoherent states, irrespective of the coupling strength. The effectiveness of QPR is illustrated in a particular case, an ensemble of Stuart-Landau oscillators, obtaining exact stability boundaries of uniform and nonuniform incoherent states for a variety of couplings. An extension of QPR beyond the neighborhood of incoherence is also explored. Finally, a general QPR model with N+2M degrees of freedom is obtained for coupling through the first M harmonics.

Polaritons and excitons: Hamiltonian design for enhanced coherence

The primary questions motivating this report are: Are there ways to increase coherence and delocalization of excitation among many molecules at moderate electronic coupling strength? Coherent delocalization of excitation in disordered molecular systems is studied using numerical calculations. The results are relevant to molecular excitons, polaritons, and make connections to classical phase oscillator synchronization. In particular, it is hypothesized that it is not only the magnitude of electronic coupling relative to the standard deviation of energetic disorder that decides the limits of coherence, but that the structure of the Hamiltonian-connections between sites (or molecules) made by electronic coupling-is a significant design parameter. Inspired by synchronization phenomena in analogous systems of phase oscillators, some properties of graphs that define the structure of different Hamiltonian matrices are explored. The report focuses on eigenvalues and ensemble density matrices of various structured, random matrices. Some reasons for the special delocalization properties and robustness of polaritons in the single-excitation subspace (the star graph) are discussed. The key result of this report is that, for some classes of Hamiltonian matrix structure, coherent delocalization is not easily defeated by energy disorder, even when the electronic coupling is small compared to disorder.

Low-dimensional dynamics of phase oscillators driven by Cauchy noise

Phase oscillator systems with global sine coupling are known to exhibit low-dimensional dynamics. In this paper, such characteristics are extended to phase oscillator systems driven by Cauchy noise. The low-dimensional dynamics solution agreed well with the numerical simulations of noise-driven phase oscillators in the present study. The low-dimensional dynamics of identical oscillators with Cauchy noise coincided with those of heterogeneous oscillators with Cauchy-distributed natural frequencies. This allows for the study of noise-driven identical oscillator systems through heterogeneous oscillators without noise and vice versa.

Coupled Möbius maps as a tool to model Kuramoto phase synchronization

We propose Möbius maps as a tool to model synchronization phenomena in coupled phase oscillators. Not only does the map provide fast computation of phase synchronization, it also reflects the underlying group structure of the sinusoidally coupled continuous phase dynamics. We study map versions of various known continuous-time collective dynamics, such as the synchronization transition in the Kuramoto-Sakaguchi model of nonidentical oscillators, chimeras in two coupled populations of identical phase oscillators, and Kuramoto-Battogtokh chimeras on a ring, and demonstrate similarities and differences between the iterated map models and their known continuous-time counterparts.

Controlling collective synchrony in oscillatory ensembles by precisely timed pulses

We present an efficient technique for control of synchrony in a globally coupled ensemble by pulsatile action. We assume that we can observe the collective oscillation and can stimulate all elements of the ensemble simultaneously. We pay special attention to the minimization of intervention into the system. The key idea is to stimulate only at the most sensitive phase. To find this phase, we implement an adaptive feedback control. Estimating the instantaneous phase of the collective mode on the fly, we achieve efficient suppression using a few pulses per oscillatory cycle. We discuss the possible relevance of the results for neuroscience, namely, for the development of advanced algorithms for deep brain stimulation, a medical technique used to treat Parkinson's disease.

Model reduction for the collective dynamics of globally coupled oscillators: From finite networks to the thermodynamic limit

Model reduction techniques have been widely used to study the collective behavior of globally coupled oscillators. However, most approaches assume that there are infinitely many oscillators. Here, we propose a new ansatz, based on the collective coordinate approach, that reproduces the collective dynamics of the Kuramoto model for finite networks to high accuracy, yields the same bifurcation structure in the thermodynamic limit of infinitely many oscillators as previous approaches, and additionally captures the dynamics of the order parameter in the thermodynamic limit, including critical slowing down that results from a cascade of saddle-node bifurcations.

Dense networks that do not synchronize and sparse ones that do

Consider any network of n identical Kuramoto oscillators in which each oscillator is coupled bidirectionally with unit strength to at least μ(n-1) other oscillators. Then, there is a critical value of μ above which the system is guaranteed to converge to the in-phase synchronous state for almost all initial conditions. The precise value of μ remains unknown. In 2018, Ling, Xu, and Bandeira proved that if each oscillator is coupled to at least 79.29% of all the others, global synchrony is ensured. In 2019, Lu and Steinerberger improved this bound to 78.89%. Here, we find clues that the critical connectivity may be exactly 75%. Our methods yield a slight improvement on the best known lower bound on the critical connectivity from 68.18% to 68.28%. We also consider the opposite end of the connectivity spectrum, where the networks are sparse rather than dense. In this regime, we ask how few edges one needs to add to a ring of n oscillators to turn it into a globally synchronizing network. We prove a partial result: all the twisted states in a ring of size n=2^m can be destabilized by adding just O(nlog₂⁡n) edges. To finish the proof, one needs to rule out all other candidate attractors. We have done this for n≤8 but the problem remains open for larger n. Thus, even for systems as simple as Kuramoto oscillators, much remains to be learned about dense networks that do not globally synchronize and sparse ones that do.

Understanding the dynamics of biological and neural oscillator networks through exact mean-field reductions: a review

Many biological and neural systems can be seen as networks of interacting periodic processes. Importantly, their functionality, i.e., whether these networks can perform their function or not, depends on the emerging collective dynamics of the network. Synchrony of oscillations is one of the most prominent examples of such collective behavior and has been associated both with function and dysfunction. Understanding how network structure and interactions, as well as the microscopic properties of individual units, shape the emerging collective dynamics is critical to find factors that lead to malfunction. However, many biological systems such as the brain consist of a large number of dynamical units. Hence, their analysis has either relied on simplified heuristic models on a coarse scale, or the analysis comes at a huge computational cost. Here we review recently introduced approaches, known as the Ott-Antonsen and Watanabe-Strogatz reductions, allowing one to simplify the analysis by bridging small and large scales. Thus, reduced model equations are obtained that exactly describe the collective dynamics for each subpopulation in the oscillator network via few collective variables only. The resulting equations are next-generation models: Rather than being heuristic, they exactly link microscopic and macroscopic descriptions and therefore accurately capture microscopic properties of the underlying system. At the same time, they are sufficiently simple to analyze without great computational effort. In the last decade, these reduction methods have become instrumental in understanding how network structure and interactions shape the collective dynamics and the emergence of synchrony. We review this progress based on concrete examples and outline possible limitations. Finally, we discuss how linking the reduced models with experimental data can guide the way towards the development of new treatment approaches, for example, for neurological disease.

Synaptic Diversity Suppresses Complex Collective Behavior in Networks of Theta Neurons

Comprehending how the brain functions requires an understanding of the dynamics of neuronal assemblies. Previous work used a mean-field reduction method to determine the collective dynamics of a large heterogeneous network of uniformly and globally coupled theta neurons, which are a canonical formulation of Type I neurons. However, in modeling neuronal networks, it is unreasonable to assume that the coupling strength between every pair of neurons is identical. The goal in the present work is to analytically examine the collective macroscopic behavior of a network of theta neurons that is more realistic in that it includes heterogeneity in the coupling strength as well as in neuronal excitability. We consider the occurrence of dynamical structures that give rise to complicated dynamics via bifurcations of macroscopic collective quantities, concentrating on two biophysically relevant cases: (1) predominantly excitable neurons with mostly excitatory connections, and (2) predominantly spiking neurons with inhibitory connections. We find that increasing the synaptic diversity moves these dynamical structures to distant extremes of parameter space, leaving simple collective equilibrium states in the physiologically relevant region. We also study the node vs. focus nature of stable macroscopic equilibrium solutions and discuss our results in the context of recent literature.