Kuramoto Model for Excitation-Inhibition-Based Oscillations

Solving the kinetic Ising model with nonreciprocity

Nonreciprocal interactions are a generic feature of nonequilibrium systems. We define a nonreciprocal generalization of the kinetic Ising model in one spatial dimension. We solve the model exactly using two different approaches for infinite, semi-infinite, and finite systems with either periodic or open boundary conditions. The exact solution allows us to explore a range of novel phenomena tied to nonreciprocity like nonreciprocity induced frustration and wave phenomena with interesting parity-dependence for finite systems of size N. We study dynamical questions like the approach to equilibrium with various boundary conditions. We find different regimes, separated by Nth-order exceptional points, which can be classified as overdamped, underdamped, or critically damped phases. Despite these different regimes, long-time order is only present at zero temperature. Additionally, we explore the low-energy behavior of the system in various limits, including the aging and spatiotemporal Porod regimes, demonstrating that nonreciprocity induces unique scaling behavior at zero temperature. Lastly, we present general results for systems where spins interact with no more than two spins, outlining the conditions under which long-time order may exist.

Deviation of a system of nonreciprocally coupled harmonic oscillators from a conservative system

Discrete systems of coupled linear mechanical oscillators with nonreciprocal interaction are a model for a variety of physical systems. In general, the presence of nonreciprocal interactions renders their dynamics nonconservative, but under certain conditions it remains conservative. In this paper we show which thermodynamic properties induced by nonreciprocity can be observed in conservative systems and which are specific to nonconservative systems. To this end, we formulate a criterion for identifying conservative systems and construct a measure to quantify the deviation from conservativity.

Analyzing top-down visual attention in the context of gamma oscillations: a layer- dependent network-of- networks approach

Top-down visual attention is a fundamental cognitive process that allows individuals to selectively attend to salient visual stimuli in the environment. Recent empirical findings have revealed that gamma oscillations participate in the modulation of visual attention. However, computational studies face challenges when analyzing the attentional process in the context of gamma oscillation due to the unstable nature of gamma oscillations and the complexity induced by the layered fashion in the visual cortex. In this study, we propose a layer-dependent network-of-networks approach to analyze such attention with gamma oscillations. The model is validated by reproducing empirical findings on orientation preference and the enhancement of neuronal response due to top-down attention. We perform parameter plane analysis to classify neuronal responses into several patterns and find that the neuronal response to sensory and attention signals was modulated by the heterogeneity of the neuronal population. Furthermore, we revealed a counter-intuitive scenario that the excitatory populations in layer 2/3 and layer 5 exhibit opposite responses to the attentional input. By modification of the original model, we confirmed layer 6 plays an indispensable role in such cases. Our findings uncover the layer-dependent dynamics in the cortical processing of visual attention and open up new possibilities for further research on layer-dependent properties in the cerebral cortex.

How synaptic function controls critical transitions in spiking neuron networks: insight from a Kuramoto model reduction

The dynamics of synaptic interactions within spiking neuron networks play a fundamental role in shaping emergent collective behavior. This paper studies a finite-size network of quadratic integrate-and-fire neurons interconnected via a general synaptic function that accounts for synaptic dynamics and time delays. Through asymptotic analysis, we transform this integrate-and-fire network into the Kuramoto-Sakaguchi model, whose parameters are explicitly expressed via synaptic function characteristics. This reduction yields analytical conditions on synaptic activation rates and time delays determining whether the synaptic coupling is attractive or repulsive. Our analysis reveals alternating stability regions for synchronous and partially synchronous firing, dependent on slow synaptic activation and time delay. We also demonstrate that the reduced microscopic model predicts the emergence of synchronization, weakly stable cyclops states, and non-stationary regimes remarkably well in the original integrate-and-fire network and its theta neuron counterpart. Our reduction approach promises to open the door to rigorous analysis of rhythmogenesis in networks with synaptic adaptation and plasticity.

Nonreciprocal Dicke Model

We investigate the physics of an open two-component Dicke model, where the light field mediates nonreciprocal interactions between two spin species. We show that the model, which we dub nonreciprocal Dicke model, exhibits a discrete parity-time (PT) symmetry and we characterize the emergence of a nonstationary phase, so far explained in terms of dissipation-induced instability, as spontaneous breaking of PT symmetry. We further show that such PT symmetry breaking embodies an instance of a nonreciprocal phase transition, a concept recently introduced by Fruchart et al. [Nature (London) 592, 363 (2021)NATUAS0028-083610.1038/s41586-021-03375-9]. Remarkably, the phase transition in our model does not necessitate the presence of any underlying broken symmetry or exceptional points in the spectrum, both believed to be essential requirements for nonreciprocal phase transitions. Our results establish driven-dissipative light-matter systems as a new avenue for exploring nonreciprocal phase transitions and contribute to the theory of nonreciprocal collective phenomena.

Volcano transition in populations of phase oscillators with random nonreciprocal interactions

Populations of heterogeneous phase oscillators with frustrated random interactions exhibit a quasiglassy state in which the distribution of local fields is volcanoshaped. In a recent work [Phys. Rev. Lett. 120, 264102 (2018)10.1103/PhysRevLett.120.264102], the volcano transition was replicated in a solvable model using a low-rank, random coupling matrix M. We extend here that model including tunable nonreciprocal interactions, i.e., M^{T}≠M. More specifically, we formulate two different solvable models. In both of them the volcano transition persists if matrix elements M_{jk} and M_{kj} are enough correlated. Our numerical simulations fully confirm the analytical results. To put our work in a wider context, we also investigate numerically the volcano transition in the analogous model with a full-rank random coupling matrix.

Collective motion of active particles exhibiting non-reciprocal orientational interactions

We present a Brownian dynamics study of a 2D bath of active particles interacting among each other through usual steric interactions and, additionally, via non-reciprocal avoidant orientational interactions. We motivate them by the fact that the two flagella of the alga Chlamydomonas interact sterically with nearby surfaces such that a torque acts on the alga. As expected, in most cases such interactions disrupt the motility-induced particle clustering in active baths. Surprisingly, however, we find that the active particles can self-organize into collectively moving flocks if the range of non-reciprocal interactions is close to that of steric interactions. We observe that the flocking motion can manifest itself through a variety of structural forms, spanning from single dense bands to multiple moderately-dense stripes, which are highly dynamic. The flocking order parameter is found to be only weakly dependent on the underlying flock structure. Together with the variance of the local-density distribution, one can clearly group the flocking motion into the two separate band and dynamic-stripes states.

Interacting Oscillators with Fluctuating Coupling: Mode Mixing without Cross-Correlations

Coupled oscillators are one of the basic models in nonlinear dynamics. Here, we study numerically and analytically the spectra of two harmonic oscillators with stochastically fluctuating coupling and driving forces reproducing a thermostat. We show that, even at small coupling, vanishing on average, the oscillation spectra exhibit mixing, even though no cross-correlations exists between the oscillators. Our results reveal a new mechanism of mode mixing for stochastically uncorrelated systems that is crucial for analysis of spectra in various systems, from simple liquids to living systems.

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.

Sandpile cascades on oscillator networks: The BTW model meets Kuramoto

Cascading failures abound in complex systems and the Bak-Tang-Weisenfeld (BTW) sandpile model provides a theoretical underpinning for their analysis. Yet, it does not account for the possibility of nodes having oscillatory dynamics, such as in power grids and brain networks. Here, we consider a network of Kuramoto oscillators upon which the BTW model is unfolding, enabling us to study how the feedback between the oscillatory and cascading dynamics can lead to new emergent behaviors. We assume that the more out-of-sync a node is with its neighbors, the more vulnerable it is and lower its load-carrying capacity accordingly. Also, when a node topples and sheds load, its oscillatory phase is reset at random. This leads to novel cyclic behavior at an emergent, long timescale. The system spends the bulk of its time in a synchronized state where load builds up with minimal cascades. Yet, eventually, the system reaches a tipping point where a large cascade triggers a "cascade of larger cascades," which can be classified as a dragon king event. The system then undergoes a short transient back to the synchronous, buildup phase. The coupling between capacity and synchronization gives rise to endogenous cascade seeds in addition to the standard exogenous ones, and we show their respective roles. We establish the phenomena from numerical studies and develop the accompanying mean-field theory to locate the tipping point, calculate the load in the system, determine the frequency of the long-time oscillations, and find the distribution of cascade sizes during the buildup phase.

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.

The Sakaguchi-Kuramoto model in presence of asymmetric interactions that break phase-shift symmetry

The celebrated Kuramoto model provides an analytically tractable framework to study spontaneous collective synchronization and comprises globally coupled limit-cycle oscillators interacting symmetrically with one another. The Sakaguchi-Kuramoto model is a generalization of the basic model that considers the presence of a phase lag parameter in the interaction, thereby making it asymmetric between oscillator pairs. Here, we consider a further generalization by adding an interaction that breaks the phase-shift symmetry of the model. The highlight of our study is the unveiling of a very rich bifurcation diagram comprising of both oscillatory and non-oscillatory synchronized states as well as an incoherent state: There are regions of two-state as well as an interesting and hitherto unexplored three-state coexistence arising from asymmetric interactions in our model.

Non-reciprocal phase transitions

Out of equilibrium, a lack of reciprocity is the rule rather than the exception. Non-reciprocity occurs, for instance, in active matter^1-6, non-equilibrium systems^7-9, networks of neurons^10,11, social groups with conformist and contrarian members¹², directional interface growth phenomena^13-15 and metamaterials^16-20. Although wave propagation in non-reciprocal media has recently been closely studied^1,16-20, less is known about the consequences of non-reciprocity on the collective behaviour of many-body systems. Here we show that non-reciprocity leads to time-dependent phases in which spontaneously broken continuous symmetries are dynamically restored. We illustrate this mechanism with simple robotic demonstrations. The resulting phase transitions are controlled by spectral singularities called exceptional points²¹. We describe the emergence of these phases using insights from bifurcation theory^22,23 and non-Hermitian quantum mechanics^24,25. Our approach captures non-reciprocal generalizations of three archetypal classes of self-organization out of equilibrium: synchronization, flocking and pattern formation. Collective phenomena in these systems range from active time-(quasi)crystals to exceptional-point-enforced pattern formation and hysteresis. Our work lays the foundation for a general theory of critical phenomena in systems whose dynamics is not governed by an optimization principle.

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

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

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.

The Winfree model with non-infinitesimal phase-response curve: Ott-Antonsen theory

A novel generalization of the Winfree model of globally coupled phase oscillators, representing phase reduction under finite coupling, is studied analytically. We consider interactions through a non-infinitesimal (or finite) phase-response curve (PRC), in contrast to the infinitesimal PRC of the original model. For a family of non-infinitesimal PRCs, the global dynamics is captured by one complex-valued ordinary differential equation resorting to the Ott-Antonsen ansatz. The phase diagrams are thereupon obtained for four illustrative cases of non-infinitesimal PRC. Bistability between collective synchronization and full desynchronization is observed in all cases.

Spatial correlations in a finite-range Kuramoto model

We study spatial correlations between oscillator phases in the steady state of a Kuramoto model, in which phase oscillators that are randomly distributed in space interact with constant strength but within a limited range. Such a model could be relevant, for example, in the synchronization of gene expression oscillations in cells, where only oscillations of neighboring cells are coupled through cell-cell contacts. We analytically infer spatial phase-phase correlation functions from the known steady-state distribution of oscillators for the case of homogenous frequencies and show that these can contain information about the range and strength of interactions, provided the noise in the system can be estimated. We suggest a method for the latter, and also explore when correlations appear to be ergodic in this model, which would enable an experimental measurement of correlation functions to utilize temporal averages. Simulations show that our techniques also give qualitative results for the model with heterogenous frequencies. We illustrate our results by comparison with experimental data on genetic oscillations in the segmentation clock of zebrafish embryos.

The role of timescale separation in oscillatory ensembles with competitive coupling

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

Low-dimensional dynamics for higher-order harmonic, globally coupled phase-oscillator ensembles

The Kuramoto model, despite its popularity as a mean-field theory for many synchronization phenomenon of oscillatory systems, is limited to a first-order harmonic coupling of phases. For higher-order coupling, there only exists a low-dimensional theory in the thermodynamic limit. In this paper, we extend the formulation used by Watanabe and Strogatz to obtain a low-dimensional description of a system of arbitrary size of identical oscillators coupled all-to-all via their higher-order modes. To demonstrate an application of the formulation, we use a second harmonic globally coupled model, with a mean-field equal to the square of the Kuramoto mean-field. This model is known to exhibit asymmetrical clustering in previous numerical studies. We try to explain the phenomenon of asymmetrical clustering using the analytical theory developed here, as well as discuss certain phenomena not observed at the level of first-order harmonic coupling.

Oscillations and collective excitability in a model of stochastic neurons under excitatory and inhibitory coupling

We study a model with excitable neurons modeled as stochastic units with three states, representing quiescence, firing, and refractoriness. The transition rates between quiescence and firing depend exponentially on the number of firing neighbors, whereas all other rates are kept constant. This model class was shown to exhibit collective oscillations (synchronization) if neurons are spiking autonomously, but not if neurons are in the excitable regime. In both cases, neurons were restricted to interact through excitatory coupling. Here we show that a plethora of collective phenomena appear if inhibitory coupling is added. Besides the usual transition between an absorbing and an active phase, the model with excitatory and inhibitory neurons can also undergo reentrant transitions to an oscillatory phase. In the mean-field description, oscillations can emerge through supercritical or subcritical Hopf bifurcations, as well as through infinite period bifurcations. The model has bistability between active and oscillating behavior, as well as collective excitability, a regime where the system can display a peak of global activity when subject to a sufficiently strong perturbation. We employ a variant of the Shinomoto-Kuramoto order parameter to characterize the phase transitions and their system-size dependence.