Phase autoencoder for limit-cycle oscillators.

We present a phase autoencoder that encodes the asymptotic phase of a limit-cycle oscillator, a fundamental quantity characterizing its synchronization dynamics. This autoencoder is trained in such a way that its latent variables directly represent the asymptotic phase of the oscillator. The trained autoencoder can perform two functions without relying on the mathematical model of the oscillator: first, it can evaluate the asymptotic phase and the phase sensitivity function of the oscillator; second, it can reconstruct the oscillator state on the limit cycle in the original space from the phase value as an input. Using several examples of limit-cycle oscillators, we demonstrate that the asymptotic phase and the phase sensitivity function can be estimated only from time-series data by the trained autoencoder. We also present a simple method for globally synchronizing two oscillators as an application of the trained autoencoder.

Higher-order interactions induce anomalous transitions to synchrony.

We analyze the simplest model of identical coupled phase oscillators subject to two-body and three-body interactions with permutation symmetry and phase lags. This model is derived from an ensemble of weakly coupled nonlinear oscillators by phase reduction, where the first and second harmonic interactions with phase lags naturally appear. Our study indicates that the higher-order interactions induce anomalous transitions to synchrony. Unlike the conventional Kuramoto model, higher-order interactions lead to anomalous phenomena such as multistability of full synchronization, incoherent, and two-cluster states, and transitions to synchrony through slow switching and clustering. Phase diagrams of the dynamical regimes are constructed theoretically and verified by direct numerical simulations. We also show that similar transition scenarios are observed even if a small heterogeneity in the oscillators' frequency is included.

Dynamic mode decomposition for Koopman spectral analysis of elementary cellular automata.

We apply dynamic mode decomposition (DMD) to elementary cellular automata (ECA). Three types of DMD methods are considered, and the reproducibility of the system dynamics and Koopman eigenvalues from observed time series is investigated. While standard DMD fails to reproduce the system dynamics and Koopman eigenvalues associated with a given periodic orbit in some cases, Hankel DMD with delay-embedded time series improves reproducibility. However, Hankel DMD can still fail to reproduce all the Koopman eigenvalues in specific cases. We propose an extended DMD method for ECA that uses nonlinearly transformed time series with discretized Walsh functions and show that it can completely reproduce the dynamics and Koopman eigenvalues. Linear-algebraic backgrounds for the reproducibility of the system dynamics and Koopman eigenvalues are also discussed.

Phase-amplitude reduction and optimal phase locking of collectively oscillating networks.

We present a phase-amplitude reduction framework for analyzing collective oscillations in networked dynamical systems. The framework, which builds on the phase reduction method, takes into account not only the collective dynamics on the limit cycle but also deviations from it by introducing amplitude variables and using them with the phase variable. The framework allows us to study how networks react to applied inputs or coupling, including their synchronization and phase locking, while capturing the deviations of the network states from the unperturbed dynamics. Numerical simulations are used to demonstrate the effectiveness of the framework for networks composed of FitzHugh-Nagumo elements. The resulting phase-amplitude equations can be used in deriving optimal periodic waveforms or introducing feedback control for achieving fast phase locking while stabilizing the collective oscillations.

A definition of the asymptotic phase for quantum nonlinear oscillators from the Koopman operator viewpoint.

We propose a definition of the asymptotic phase for quantum nonlinear oscillators from the viewpoint of the Koopman operator theory. The asymptotic phase is a fundamental quantity for the analysis of classical limit-cycle oscillators, but it has not been defined explicitly for quantum nonlinear oscillators. In this study, we define the asymptotic phase for quantum oscillatory systems by using the eigenoperator of the backward Liouville operator associated with the fundamental oscillation frequency. By using the quantum van der Pol oscillator with a Kerr effect as an example, we illustrate that the proposed asymptotic phase appropriately yields isochronous phase values in both semiclassical and strong quantum regimes.

Fast optimal entrainment of limit-cycle oscillators by strong periodic inputs via phase-amplitude reduction and Floquet theory.

Optimal entrainment of limit-cycle oscillators by strong periodic inputs is studied on the basis of the phase-amplitude reduction and Floquet theory. Two methods for deriving the input waveforms that keep the system state close to the original limit cycle are proposed, which enable the use of strong inputs for entrainment. The first amplitude-feedback method uses feedback control to suppress deviations of the system state from the limit cycle, while the second amplitude-penalty method seeks an input waveform that does not excite large deviations from the limit cycle in the feedforward framework. Optimal entrainment of the van der Pol and Willamowski-Rössler oscillators with real or complex Floquet exponents is analyzed as examples. It is demonstrated that the proposed methods can achieve considerably faster entrainment and provide wider entrainment ranges than the conventional method that relies only on phase reduction.

Sparse optimization of mutual synchronization in collectively oscillating networks.

We consider a pair of collectively oscillating networks of dynamical elements and optimize their internetwork coupling for efficient mutual synchronization based on the phase reduction theory developed by Nakao et al. [Chaos 28, 045103 (2018)]. The dynamical equations describing a pair of weakly coupled networks are reduced to a pair of coupled phase equations, and the linear stability of the synchronized state between the networks is represented as a function of the internetwork coupling matrix. We seek the optimal coupling by minimizing the Frobenius and L1 norms of the internetwork coupling matrix for the prescribed linear stability of the synchronized state. Depending on the norm, either a dense or sparse internetwork coupling yielding efficient mutual synchronization of the networks is obtained. In particular, a sparse yet resilient internetwork coupling is obtained by L1-norm optimization with additional constraints on the individual connection weights.

Koopman spectral analysis of elementary cellular automata.

We perform a Koopman spectral analysis of elementary cellular automata (ECA). By lifting the system dynamics using a one-hot representation of the system state, we derive a matrix representation of the Koopman operator as the transpose of the adjacency matrix of the state-transition network. The Koopman eigenvalues are either zero or on the unit circle in the complex plane, and the associated Koopman eigenfunctions can be explicitly constructed. From the Koopman eigenvalues, we can judge the reversibility, determine the number of connected components in the state-transition network, evaluate the period of asymptotic orbits, and derive the conserved quantities for each system. We numerically calculate the Koopman eigenvalues of all rules of ECA on a one-dimensional lattice of 13 cells with periodic boundary conditions. It is shown that the spectral properties of the Koopman operator reflect Wolfram's classification of ECA.

Spectral analysis of the Koopman operator for partial differential equations.

We provide an overview of the Koopman-operator analysis for a class of partial differential equations describing relaxation of the field variable to a stable stationary state. We introduce Koopman eigenfunctionals of the system and use the notion of conjugacy to develop spectral expansion of the Koopman operator. For linear systems such as the diffusion equation, the Koopman eigenfunctionals can be expressed as linear functionals of the field variable. The notion of inertial manifolds is shown to correspond to joint zero level sets of Koopman eigenfunctionals, and the notion of isostables is defined as the level sets of the slowest decaying Koopman eigenfunctional. Linear diffusion equation, nonlinear Burgers equation, and nonlinear phase-diffusion equation are analyzed as examples.

On the concept of dynamical reduction: the case of coupled oscillators.

An overview is given on two representative methods of dynamical reduction known as centre-manifold reduction and phase reduction. These theories are presented in a somewhat more unified fashion than the theories in the past. The target systems of reduction are coupled limit-cycle oscillators. Particular emphasis is placed on the remarkable structural similarity existing between these theories. While the two basic principles, i.e. (i) reduction of dynamical degrees of freedom and (ii) transformation of reduced evolution equation to a canonical form, are shared commonly by reduction methods in general, it is shown how these principles are incorporated into the above two reduction theories in a coherent manner. Regarding the phase reduction, a new formulation of perturbative expansion is presented for discrete populations of oscillators. The style of description is intended to be so informal that one may digest, without being bothered with technicalities, what has been done after all under the word reduction. This article is part of the theme issue 'Coupling functions: dynamical interaction mechanisms in the physical, biological and social sciences'.

Phase reduction and synchronization of a network of coupled dynamical elements exhibiting collective oscillations.

A general phase reduction method for a network of coupled dynamical elements exhibiting collective oscillations, which is applicable to arbitrary networks of heterogeneous dynamical elements, is developed. A set of coupled adjoint equations for phase sensitivity functions, which characterize the phase response of the collective oscillation to small perturbations applied to individual elements, is derived. Using the phase sensitivity functions, collective oscillation of the network under weak perturbation can be described approximately by a one-dimensional phase equation. As an example, mutual synchronization between a pair of collectively oscillating networks of excitable and oscillatory FitzHugh-Nagumo elements with random coupling is studied.

Phase-amplitude reduction of transient dynamics far from attractors for limit-cycling systems.

Phase reduction framework for limit-cycling systems based on isochrons has been used as a powerful tool for analyzing the rhythmic phenomena. Recently, the notion of isostables, which complements the isochrons by characterizing amplitudes of the system state, i.e., deviations from the limit-cycle attractor, has been introduced to describe the transient dynamics around the limit cycle [Wilson and Moehlis, Phys. Rev. E 94, 052213 (2016)]. In this study, we introduce a framework for a reduced phase-amplitude description of transient dynamics of stable limit-cycling systems. In contrast to the preceding study, the isostables are treated in a fully consistent way with the Koopman operator analysis, which enables us to avoid discontinuities of the isostables and to apply the framework to system states far from the limit cycle. We also propose a new, convenient bi-orthogonalization method to obtain the response functions of the amplitudes, which can be interpreted as an extension of the adjoint covariant Lyapunov vector to transient dynamics in limit-cycling systems. We illustrate the utility of the proposed reduction framework by estimating the optimal injection timing of external input that efficiently suppresses deviations of the system state from the limit cycle in a model of a biochemical oscillator.

Phase reduction method for strongly perturbed limit cycle oscillators.

The phase reduction method for limit cycle oscillators subjected to weak perturbations has significantly contributed to theoretical investigations of rhythmic phenomena. We here propose a generalized phase reduction method that is also applicable to strongly perturbed limit cycle oscillators. The fundamental assumption of our method is that the perturbations can be decomposed into a slowly varying component as compared to the amplitude relaxation time and remaining weak fluctuations. Under this assumption, we introduce a generalized phase parameterized by the slowly varying component and derive a closed equation for the generalized phase describing the oscillator dynamics. The proposed method enables us to explore a broader class of rhythmic phenomena, in which the shape and frequency of the oscillation may vary largely because of the perturbations. We illustrate our method by analyzing the synchronization dynamics of limit cycle oscillators driven by strong periodic signals. It is shown that the proposed method accurately predicts the synchronization properties of the oscillators, while the conventional method does not.

Collective phase description of oscillatory convection.

We formulate a theory for the collective phase description of oscillatory convection in Hele-Shaw cells. It enables us to describe the dynamics of the oscillatory convection by a single degree of freedom which we call the collective phase. The theory can be considered as a phase reduction method for limit-cycle solutions in infinite-dimensional dynamical systems, namely, stable time-periodic solutions to partial differential equations, representing the oscillatory convection. We derive the phase sensitivity function, which quantifies the phase response of the oscillatory convection to weak perturbations applied at each spatial point, and analyze the phase synchronization between two weakly coupled Hele-Shaw cells exhibiting oscillatory convection on the basis of the derived phase equations.

Adjoint method provides phase response functions for delay-induced oscillations.

Limit-cycle oscillations induced by time delay are widely observed in various systems, but a systematic phase-reduction theory for them has yet to be developed. Here we present a practical theoretical framework to calculate the phase response function Z(θ), a fundamental quantity for the theory, of delay-induced limit cycles with infinite-dimensional phase space. We show that Z(θ) can be obtained as a zero eigenfunction of the adjoint equation associated with an appropriate bilinear form for the delay differential equations. We confirm the validity of the proposed framework for two biological oscillators and demonstrate that the derived phase equation predicts intriguing multimodal locking behavior.

Noise-tuned bursting in a Hedgehog burster.

Noise can shape the firing behaviors of neurons. Here, we show that noise acting on the fast variable of the Hedgehog burster can tune the spike counts of bursts via the self-induced stochastic resonance (SISR) phenomenon. Using the distance matching condition, the critical transition positions on the slow manifolds can be predicted and the stochastic periodic orbits for various noise strengths are obtained. The critical transition positions on the slow manifold with non-monotonic potential differences exhibit a staircase-like dependence on the noise strength, which is also revealed by the stepwise change in the period of the stochastic periodic orbit. The noise-tuned bursting is more coherent within each step while displaying mixed-mode oscillations near the boundaries between the steps. When noise is large enough, noise-induced trapping of the slow variable can be observed, where the number of coexisting traps increases with the noise strength. It is argued that the robustness of SISR underlies the generality of the results discovered in this paper.

Turing instability in quantum activator-inhibitor systems.

Turing instability is a fundamental mechanism of nonequilibrium self-organization. However, despite the universality of its essential mechanism, Turing instability has thus far been investigated mostly in classical systems. In this study, we show that Turing instability can occur in a quantum dissipative system and analyze its quantum features such as entanglement and the effect of measurement. We propose a degenerate parametric oscillator with nonlinear damping in quantum optics as a quantum activator-inhibitor unit and demonstrate that a system of two such units can undergo Turing instability when diffusively coupled with each other. The Turing instability induces nonuniformity and entanglement between the two units and gives rise to a pair of nonuniform states that are mixed due to quantum noise. Further performing continuous measurement on the coupled system reveals the nonuniformity caused by the Turing instability. Our results extend the universality of the Turing mechanism to the quantum realm and may provide a novel perspective on the possibility of quantum nonequilibrium self-organization and its application in quantum technologies.

Estimating asymptotic phase and amplitude functions of limit-cycle oscillators from time series data.

We propose a method for estimating the asymptotic phase and amplitude functions of limit-cycle oscillators using observed time series data without prior knowledge of their dynamical equations. The estimation is performed by polynomial regression and can be solved as a convex optimization problem. The validity of the proposed method is numerically illustrated by using two-dimensional limit-cycle oscillators as examples. As an application, we demonstrate data-driven fast entrainment with amplitude suppression using the optimal periodic input derived from the estimated phase and amplitude functions.

Dynamics of limit-cycle oscillators subject to general noise.

The phase description is a powerful tool for analyzing noisy limit-cycle oscillators. The method, however, has found only limited applications so far, because the present theory is applicable only to Gaussian noise while noise in the real world often has non-Gaussian statistics. Here, we provide the phase reduction method for limit-cycle oscillators subject to general, colored and non-Gaussian, noise including a heavy-tailed one. We derive quantifiers like mean frequency, diffusion constant, and the Lyapunov exponent to confirm consistency of the results. Applying our results, we additionally study a resonance between the phase and noise.

Effective long-time phase dynamics of limit-cycle oscillators driven by weak colored noise.

An effective white-noise Langevin equation is derived that describes long-time phase dynamics of a limit-cycle oscillator driven by weak stationary colored noise. Effective drift and diffusion coefficients are given in terms of the phase sensitivity of the oscillator and the correlation function of the noise, and are explicitly calculated for oscillators with sinusoidal phase sensitivity functions driven by two typical colored Gaussian processes. The results are verified by numerical simulations using several types of stochastic or chaotic noise. The drift and diffusion coefficients of oscillators driven by chaotic noise exhibit anomalous dependence on the oscillator frequency, reflecting the peculiar power spectrum of the chaotic noise.