Low-dimensional dynamics of phase oscillators driven by Cauchy noise

Collective synchronization through noise cancellation

After decades of study, there are only two known mechanisms to induce global synchronization in a population of oscillators: Deterministic coupling and common forcing. The inclusion of independent noise in these models typically serves to drive disorder, increasing the stability of the incoherent state. Here we show that the reverse is also possible. We propose and analyze a simple general model of purely noise coupled oscillators. In the first explicit choice of noise coupling, we find the linear response around incoherence is identical to that of the paradigmatic Kuramoto model but exhibits binary phase locking instead of full coherence. We characterize the phase diagram, stationary states, and approximate low-dimensional dynamics for the model, revealing the curious behavior of this mechanism of synchronization. In the second minimal case we connect the final synchronized state to the initial conditions of the system.

Exact low-dimensional description for fast neural oscillations with low firing rates

Recently, low-dimensional models of neuronal activity have been exactly derived for large networks of deterministic, quadratic integrate-and-fire (QIF) neurons. Such firing rate models (FRM) describe the emergence of fast collective oscillations (>30 Hz) via the frequency locking of a subset of neurons to the global oscillation frequency. However, the suitability of such models to describe realistic neuronal states is seriously challenged by the fact that during episodes of fast collective oscillations, neuronal discharges are often very irregular and have low firing rates compared to the global oscillation frequency. Here we extend the theory to derive exact FRM for QIF neurons to include noise and show that networks of stochastic neurons displaying irregular discharges at low firing rates during episodes of fast oscillations are governed by exactly the same evolution equations as deterministic networks. Our results reconcile two traditionally confronted views on neuronal synchronization and upgrade the applicability of exact FRM to describe a broad range of biologically realistic neuronal states.

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.

Exact finite-dimensional reduction for a population of noisy oscillators and its link to Ott-Antonsen and Watanabe-Strogatz theories

Populations of globally coupled phase oscillators are described in the thermodynamic limit by kinetic equations for the distribution densities or, equivalently, by infinite hierarchies of equations for the order parameters. Ott and Antonsen [Chaos 18, 037113 (2008)] have found an invariant finite-dimensional subspace on which the dynamics is described by one complex variable per population. For oscillators with Cauchy distributed frequencies or for those driven by Cauchy white noise, this subspace is weakly stable and, thus, describes the asymptotic dynamics. Here, we report on an exact finite-dimensional reduction of the dynamics outside of the Ott-Antonsen subspace. We show that the evolution from generic initial states can be reduced to that of three complex variables, plus a constant function. For identical noise-free oscillators, this reduction corresponds to the Watanabe-Strogatz system of equations [Watanabe and Strogatz, Phys. Rev. Lett. 70, 2391 (1993)]. We discuss how the reduced system can be used to explore the transient dynamics of perturbed ensembles.

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.

Hierarchy of Exact Low-Dimensional Reductions for Populations of Coupled Oscillators

We consider an ensemble of phase oscillators in the thermodynamic limit, where it is described by a kinetic equation for the phase distribution density. We propose an Ansatz for the circular moments of the distribution (Kuramoto-Daido order parameters) that allows for an exact truncation at an arbitrary number of modes. In the simplest case of one mode, the Ansatz coincides with that of Ott and Antonsen [Chaos 18, 037113 (2008)CHAOEH1054-150010.1063/1.2930766]. Dynamics on the extended manifolds facilitate higher-dimensional behavior such as chaos, which we demonstrate with a simulation of a Josephson junction array. The findings are generalized for oscillators with a Cauchy-Lorentzian distribution of natural frequencies.

Lévy on-off intermittency

We present an alternative form of intermittency, Lévy on-off intermittency, which arises from multiplicative α-stable white noise close to an instability threshold. We study this problem in the linear and nonlinear regimes, both theoretically and numerically, for the case of a pitchfork bifurcation with fluctuating growth rate. We compute the stationary distribution analytically and numerically from the associated fractional Fokker-Planck equation in the Stratonovich interpretation. We characterize the system in the parameter space (α,β) of the noise, with stability parameter α∈(0,2) and skewness parameter β∈[-1,1]. Five regimes are identified in this parameter space, in addition to the well-studied Gaussian case α=2. Three regimes are located at 1<α<2, where the noise has finite mean but infinite variance. They are differentiated by β and all display a critical transition at the deterministic instability threshold, with on-off intermittency close to onset. Critical exponents are computed from the stationary distribution. Each regime is characterized by a specific form of the density and specific critical exponents, which differ starkly from the Gaussian case. A finite or infinite number of integer-order moments may converge, depending on parameters. Two more regimes are found at 0<α≤1. There, the mean of the noise diverges, and no critical transition occurs. In one case, the origin is always unstable, independently of the distance μ from the deterministic threshold. In the other case, the origin is conversely always stable, independently of μ. We thus demonstrate that an instability subject to nonequilibrium, power-law-distributed fluctuations can display substantially different properties than for Gaussian thermal fluctuations, in terms of statistics and critical behavior.

Low-dimensional description for ensembles of identical phase oscillators subject to Cauchy noise

We study ensembles of globally coupled or forced identical phase oscillators subject to independent white Cauchy noise. We demonstrate that if the oscillators are forced in several harmonics, stationary synchronous regimes can be exactly described with a finite number of complex order parameters. The corresponding distribution of phases is a product of wrapped Cauchy distributions. For sinusoidal forcing, the Ott-Antonsen low-dimensional reduction is recovered.