Dynamics of Oscillator Populations Globally Coupled with Distributed Phase Shifts.

We consider a population of globally coupled oscillators in which phase shifts in the coupling are random. We show that in the maximally disordered case, where the pairwise shifts are independent identically distributed random variables, the dynamics of a large population reduces to one without randomness in the shifts but with an effective coupling function, which is a convolution of the original coupling function with the distribution of the phase shifts. This result is valid for noisy oscillators and/or in the presence of a distribution of natural frequencies. We argue also, using the property of global asymptotic stability, that this reduction is valid in a partially disordered case, where random phase shifts are attributed to the forced units only. However, the reduction to an effective coupling in the partially disordered noise-free situation may fail if the coupling function is complex enough to ensure the multistability of locked states.

A unified quantification of synchrony in globally coupled populations with the Wiener order parameter.

We tackle the quantification of synchrony in globally coupled populations. Furthermore, we treat the problem of incomplete observations when the population mean field is unavailable, but only a small subset of units is observed. We introduce a new order parameter based on the integral of the squared autocorrelation function and demonstrate its efficiency for quantifying synchrony via monitoring general observables, regardless of whether the oscillations can be characterized in terms of the phases. Under condition of a significant irregularity in the dynamics of the coupled units, this order parameter provides a unified description of synchrony in populations of units of various complexities. The main examples include noise-induced oscillations, coupled strongly chaotic systems, and noisy periodic oscillations. Furthermore, we explore how this parameter works for the standard Kuramoto model of coupled regular-phase oscillators. The most significant advantage of our approach is its ability to infer and quantify synchrony from the observation of a small percentage of the units and even from a single unit, provided the observations are sufficiently long.

Dynamics of large oscillator populations with random interactions.

We explore large populations of phase oscillators interacting via random coupling functions. Two types of coupling terms, the Kuramoto-Daido coupling and the Winfree coupling, are considered. Under the assumption of statistical independence of the phases and the couplings, we derive reduced averaged equations with effective non-random coupling terms. As a particular example, we study interactions defined via the coupling functions that have the same shape but possess random coupling strengths and random phase shifts. While randomness in coupling strengths just renormalizes the interaction, a distribution of the phase shifts in coupling reshapes the coupling function.

Deterministic active particles in the overactive limit.

We consider two models of deterministic active particles in an external potential. In the limit where the speed of a particle is fixed, both models nearly coincide and can be formulated as a Hamiltonian system, but only if the potential is time-independent. If the particles are identical, their interaction via a potential force leads to conservative dynamics with a conserved phase volume. In contrast, the phase volume is shown to shrink for non-identical particles even if the confining potential is time-independent.

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.

Attractor-repeller collision and the heterodimensional dynamics.

We study the heterodimensional dynamics in a simple map on a three-dimensional torus. This map consists of a two-dimensional driving Anosov map and a one-dimensional driven Möbius map, and demonstrates the collision of a chaotic attractor with a chaotic repeller if parameters are varied. We explore this collision by following tangent bifurcations of the periodic orbits and establish a regime where periodic orbits with different numbers of unstable directions coexist in a chaotic set. For this situation, we construct a heterodimensional cycle connecting these periodic orbits. Furthermore, we discuss properties of the rotation number and of the nontrivial Lyapunov exponent at the collision and in the heterodimensional regime.

High-order phase reduction for coupled 2D oscillators.

Phase reduction is a general approach to describe coupled oscillatory units in terms of their phases, assuming that the amplitudes are enslaved. The coupling should be small for such reduction, but one also expects the reduction to be valid for finite coupling. This paper presents a general framework, allowing us to obtain coupling terms in higher orders of the coupling parameter for generic two-dimensional oscillators and arbitrary coupling terms. The theory is illustrated with an accurate prediction of Arnold's tongue for the van der Pol oscillator exploiting higher-order phase reduction.

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.

Confinement and Collective Escape of Active Particles.

Active matter broadly covers the dynamics of self-propelled particles. While the onset of collective behavior in homogenous active systems is relatively well understood, the effect of inhomogeneities such as obstacles and traps lacks overall clarity. Here, we study how interacting, self-propelled particles become trapped and released from a trap. We have found that captured particles aggregate into an orbiting condensate with a crystalline structure. As more particles are added, the trapped condensates escape as a whole. Our results shed light on the effects of confinement and quenched disorder in active matter.

Stochastic bursting in networks of excitable units with delayed coupling.

We investigate the phenomenon of stochastic bursting in a noisy excitable unit with multiple weak delay feedbacks, by virtue of a directed tree lattice model. We find statistical properties of the appearing sequence of spikes and expressions for the power spectral density. This simple model is extended to a network of three units with delayed coupling of a star type. We find the power spectral density of each unit and the cross-spectral density between any two units. The basic assumptions behind the analytical approach are the separation of timescales, allowing for a description of the spike train as a point process, and weakness of coupling, allowing for a representation of the action of overlapped spikes via the sum of the one-spike excitation probabilities.

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.

Mutual singularities of overlapping attractor and repeller.

We apply the concepts of relative dimensions and mutual singularities to characterize the fractal properties of overlapping attractor and repeller in chaotic dynamical systems. We consider one analytically solvable example (a generalized baker's map); two other examples, the Anosov-Möbius and the Chirikov-Möbius maps, which possess fractal attractor and repeller on a two-dimensional torus, are explored numerically. We demonstrate that although for these maps the stable and unstable directions are not orthogonal to each other, the relative Rényi and Kullback-Leibler dimensions as well as the mutual singularity spectra for the attractor and repeller can be well approximated under orthogonality assumption of two fractals.

Some elements for a history of the dynamical systems theory.

Writing a history of a scientific theory is always difficult because it requires to focus on some key contributors and to "reconstruct" some supposed influences. In the 1970s, a new way of performing science under the name "chaos" emerged, combining the mathematics from the nonlinear dynamical systems theory and numerical simulations. To provide a direct testimony of how contributors can be influenced by other scientists or works, we here collected some writings about the early times of a few contributors to chaos theory. The purpose is to exhibit the diversity in the paths and to bring some elements-which were never published-illustrating the atmosphere of this period. Some peculiarities of chaos theory are also discussed.

Solitary phase waves in a chain of autonomous oscillators.

In the present paper, we study phase waves of self-sustained oscillators with a nearest-neighbor dispersive coupling on an infinite lattice. To analyze the underlying dynamics, we approximate the lattice with a quasi-continuum (QC). The resulting partial differential model is then further reduced to the Gardner equation, which predicts many properties of the underlying solitary structures. Using an iterative procedure on the original lattice equations, we determine the shapes of solitary waves, kinks, and the flat-like solitons that we refer to as flatons. Direct numerical experiments reveal that the interaction of solitons and flatons on the lattice is notably clean. All in all, we find that both the QC and the Gardner equation predict remarkably well the discrete patterns and their dynamics.

Kantorovich-Rubinstein-Wasserstein distance between overlapping attractor and repeller.

We consider several examples of dynamical systems demonstrating overlapping attractor and repeller. These systems are constructed via introducing controllable dissipation to prototypic models with chaotic dynamics (Anosov cat map, Chirikov standard map, and incompressible three-dimensional flow of the ABC-type on a three-torus) and ergodic non-chaotic behavior (skew-shift map). We employ the Kantorovich-Rubinstein-Wasserstein distance to characterize the difference between the attractor and the repeller, in dependence on the dissipation level.

Nonlinear phase coupling functions: a numerical study.

Phase reduction is a general tool widely used to describe forced and interacting self-sustained oscillators. Here, we explore the phase coupling functions beyond the usual first-order approximation in the strength of the force. Taking the periodically forced Stuart-Landau oscillator as the paradigmatic model, we determine and numerically analyse the coupling functions up to the fourth order in the force strength. We show that the found nonlinear phase coupling functions can be used for predicting synchronization regions of the forced oscillator. This article is part of the theme issue 'Coupling functions: dynamical interaction mechanisms in the physical, biological and social sciences'.

Numerical phase reduction beyond the first order approximation.

We develop a numerical approach to reconstruct the phase dynamics of driven or coupled self-sustained oscillators. Employing a simple algorithm for computation of the phase of a perturbed system, we construct numerically the equation for the evolution of the phase. Our simulations demonstrate that the description of the dynamics solely by phase variables can be valid for rather strong coupling strengths and large deviations from the limit cycle. Coupling functions depend crucially on the coupling and are generally non-decomposable in phase response and forcing terms. We also discuss the limitations of the approach.

Stochastic bursting in unidirectionally delay-coupled noisy excitable systems.

We show that "stochastic bursting" is observed in a ring of unidirectional delay-coupled noisy excitable systems, thanks to the combinational action of time-delayed coupling and noise. Under the approximation of timescale separation, i.e., when the time delays in each connection are much larger than the characteristic duration of the spikes, the observed rather coherent spike pattern can be described by an idealized coupled point process with a leader-follower relationship. We derive analytically the statistics of the spikes in each unit, the pairwise correlations between any two units, and the spectrum of the total output from the network. Theory is in good agreement with the simulations with a network of theta-neurons.

Blinking chimeras in globally coupled rotators.

In globally coupled ensembles of identical oscillators so-called chimera states can be observed. The chimera state is a symmetry-broken regime, where a subset of oscillators forms a cluster, a synchronized population, while the rest of the system remains a collection of nonsynchronized, scattered units. We describe here a blinking chimera regime in an ensemble of seven globally coupled rotators (Kuramoto oscillators with inertia). It is characterized by a death-birth process, where a long-term stable cluster of four oscillators suddenly dissolves and is very quickly reborn with a new reshuffled configuration. We identify three different kinds of rare blinking events and give a quantitative characterization by applying stability analysis to the long-lived chaotic state and to the short-lived regular regimes that arise when the cluster dissolves.

Repulsively coupled Kuramoto-Sakaguchi phase oscillators ensemble subject to common noise.

We consider the Kuramoto-Sakaguchi model of identical coupled phase oscillators with a common noisy forcing. While common noise always tends to synchronize the oscillators, a strong repulsive coupling prevents the fully synchronous state and leads to a nontrivial distribution of oscillator phases. In previous numerical simulations, the formation of stable multicluster states has been observed in this regime. However, we argue here that because identical phase oscillators in the Kuramoto-Sakaguchi model form a partially integrable system according to the Watanabe-Strogatz theory, the formation of clusters is impossible. Integrating with various time steps reveals that clustering is a numerical artifact, explained by the existence of higher order Fourier terms in the errors of the employed numerical integration schemes. By monitoring the induced change in certain integrals of motion, we quantify these errors. We support these observations by showing, on the basis of the analysis of the corresponding Fokker-Planck equation, that two-cluster states are non-attractive. On the other hand, in ensembles of general limit cycle oscillators, such as Van der Pol oscillators, due to an anharmonic phase response function as well as additional amplitude dynamics, multiclusters can occur naturally.