Isostable reduction of periodic orbits

Reconstruction of phase-amplitude dynamics from signals of a network of oscillators

We present a novel method of reconstructing the phase-amplitude dynamics directly from signals of a network of oscillators to estimate the coupling between its nodes. For this purpose, we use the recent advances in the field of phase-amplitude reduction of oscillatory systems, which allow the representation of an uncoupled oscillatory system as a phase-amplitude oscillator in a unique form using transformations (parameterizations) related to the eigenfunctions of the Koopman operator. By combining the parameterization method and the Fourier-Laplace averaging method for finding the eigenfunctions of the Koopman operator, we developed a method of assessing the transformation functions from the signals of the interacting oscillatory systems. The resulting reconstructed dynamical system is a network of phase-amplitude oscillators with the interactions between them represented as coupling functions in phase and amplitude coordinates.

The efficiency of synchronization dynamics and the role of network syncreactivity

Synchronization of coupled oscillators is a fundamental process in both natural and artificial networks. While much work has investigated the asymptotic stability of the synchronous solution, the fundamental question of the transient behavior toward synchronization has received far less attention. In this work, we present the transverse reactivity as a metric to quantify the instantaneous rate of growth or decay of desynchronizing perturbations. We first use the transverse reactivity to design a coupling-efficient and energy-efficient synchronization strategy that involves varying the coupling strength dynamically according to the current state of the system. We find that our synchronization strategy is able to synchronize networks in both simulation and experiment over a significantly larger (often by orders of magnitude) range of coupling strengths than is possible when the coupling strength is constant. Then, we characterize the effects of network topology on the transient dynamics towards synchronization by introducing the concept of network syncreactivity: A network with a larger syncreactivity has a larger transverse reactivity at every point on the synchronization manifold, independent of the oscillator dynamics. We classify real-world examples of complex networks in terms of their syncreactivity.

Koopman analysis of the singularly perturbed van der Pol oscillator

The Koopman operator framework holds promise for spectral analysis of nonlinear dynamical systems based on linear operators. Eigenvalues and eigenfunctions of the Koopman operator, the so-called Koopman eigenvalues and Koopman eigenfunctions, respectively, mirror global properties of the system's flow. In this paper, we perform the Koopman analysis of the singularly perturbed van der Pol system. First, we show the spectral signature depending on singular perturbation: how two Koopman principal eigenvalues are ordered and what distinct shapes emerge in their associated Koopman eigenfunctions. Second, we discuss the singular limit of the Koopman operator, which is derived through the concatenation of Koopman operators for the fast and slow subsystems. From the spectral properties of the Koopman operator for the singularly perturbed system and the singular limit, we suggest that the Koopman eigenfunctions inherit geometric properties of the singularly perturbed system. These results are applicable to general planar singularly perturbed systems with stable limit cycles.

How to design optimal brain stimulation to modulate phase-amplitude coupling?

Objective.Phase-amplitude coupling (PAC), the coupling of the amplitude of a faster brain rhythm to the phase of a slower brain rhythm, plays a significant role in brain activity and has been implicated in various neurological disorders. For example, in Parkinson's disease, PAC between the beta (13-30 Hz) and gamma (30-100 Hz) rhythms in the motor cortex is exaggerated, while in Alzheimer's disease, PAC between the theta (4-8 Hz) and gamma rhythms is diminished. Modulating PAC (i.e. reducing or enhancing PAC) using brain stimulation could therefore open new therapeutic avenues. However, while it has been previously reported that phase-locked stimulation can increase PAC, it is unclear what the optimal stimulation strategy to modulate PAC might be. Here, we provide a theoretical framework to narrow down the experimental optimisation of stimulation aimed at modulating PAC, which would otherwise rely on trial and error.Approach.We make analytical predictions using a Stuart-Landau model, and confirm these predictions in a more realistic model of coupled neural populations.Main results.Our framework specifies the critical Fourier coefficients of the stimulation waveform which should be tuned to optimally modulate PAC. Depending on the characteristics of the amplitude response curve of the fast population, these components may include the slow frequency, the fast frequency, combinations of these, as well as their harmonics. We also show that the optimal balance of energy between these Fourier components depends on the relative strength of the endogenous slow and fast rhythms, and that the alignment of fast components with the fast rhythm should change throughout the slow cycle. Furthermore, we identify the conditions requiring to phase-lock stimulation to the fast and/or slow rhythms.Significance.Together, our theoretical framework lays the foundation for guiding the development of innovative and more effective brain stimulation aimed at modulating PAC for therapeutic benefit.

Phase and amplitude responses for delay equations using harmonic balance

Robust delay induced oscillations, common in nature, are often modeled by delay-differential equations (DDEs). Motivated by the success of phase-amplitude reductions for ordinary differential equations with limit cycle oscillations, there is now a growing interest in the development of analogous approaches for DDEs to understand their response to external forcing. When combined with Floquet theory, the fundamental quantities for this reduction are phase and amplitude response functions. Here, we develop a framework for their construction that utilizes the method of harmonic balance.

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.

Optimal phase-based control of strongly perturbed limit cycle oscillators using phase reduction techniques

Phase reduction is a well-established technique for analysis and control of weakly perturbed limit cycle oscillators. However, its accuracy is diminished in a strongly perturbed setting where information about the amplitude dynamics must also be considered. In this paper, we consider phase-based control of general limit cycle oscillators in both weakly and strongly perturbed regimes. For use at the strongly perturbed end of the continuum, we propose a strategy for optimal phase control of general limit cycle oscillators that uses an adaptive phase-amplitude reduced order model in conjunction with dynamic programming. This strategy can accommodate large magnitude inputs at the expense of requiring additional dimensions in the reduced order equations, thereby increasing the computational complexity. We apply this strategy to two biologically motivated prototype problems and provide direct comparisons to two related phase-based control algorithms. In situations where other commonly used strategies fail due to the application of large magnitude inputs, the adaptive phase-amplitude reduction provides a viable reduced order model while still yielding a computationally tractable control problem. These results highlight the need for discernment in reduced order model selection for limit cycle oscillators to balance the trade-off between accuracy and dimensionality.

Insights into oscillator network dynamics using a phase-isostable framework

Networks of coupled nonlinear oscillators can display a wide range of emergent behaviors under the variation of the strength of the coupling. Network equations for pairs of coupled oscillators where the dynamics of each node is described by the evolution of its phase and slowest decaying isostable coordinate have previously been shown to capture bifurcations and dynamics of the network, which cannot be explained through standard phase reduction. An alternative framework using isostable coordinates to obtain higher-order phase reductions has also demonstrated a similar descriptive ability for two oscillators. In this work, we consider the phase-isostable network equations for an arbitrary but finite number of identical coupled oscillators, obtaining conditions required for the stability of phase-locked states including synchrony. For the mean-field complex Ginzburg-Landau equation where the solutions of the full system are known, we compare the accuracy of the phase-isostable network equations and higher-order phase reductions in capturing bifurcations of phase-locked states. We find the former to be the more accurate and, therefore, employ this to investigate the dynamics of globally linearly coupled networks of Morris-Lecar neuron models (both two and many nodes). We observe qualitative correspondence between results from numerical simulations of the full system and the phase-isostable description demonstrating that in both small and large networks, the phase-isostable framework is able to capture dynamics that the first-order phase description cannot.

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 universal description of stochastic oscillators

Many systems in physics, chemistry, and biology exhibit oscillations with a pronounced random component. Such stochastic oscillations can emerge via different mechanisms, for example, linear dynamics of a stable focus with fluctuations, limit-cycle systems perturbed by noise, or excitable systems in which random inputs lead to a train of pulses. Despite their diverse origins, the phenomenology of random oscillations can be strikingly similar. Here, we introduce a nonlinear transformation of stochastic oscillators to a complex-valued function [Formula: see text](x) that greatly simplifies and unifies the mathematical description of the oscillator's spontaneous activity, its response to an external time-dependent perturbation, and the correlation statistics of different oscillators that are weakly coupled. The function [Formula: see text] (x) is the eigenfunction of the Kolmogorov backward operator with the least negative (but nonvanishing) eigenvalue λ1 = μ1 + iω1. The resulting power spectrum of the complex-valued function is exactly given by a Lorentz spectrum with peak frequency ω1 and half-width μ1; its susceptibility with respect to a weak external forcing is given by a simple one-pole filter, centered around ω1; and the cross-spectrum between two coupled oscillators can be easily expressed by a combination of the spontaneous power spectra of the uncoupled systems and their susceptibilities. Our approach makes qualitatively different stochastic oscillators comparable, provides simple characteristics for the coherence of the random oscillation, and gives a framework for the description of weakly coupled oscillators.

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.

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.

Control of coupled neural oscillations using near-periodic inputs

Deep brain stimulation (DBS) is a commonly used treatment for medication resistant Parkinson's disease and is an emerging treatment for other neurological disorders. More recently, phase-specific adaptive DBS (aDBS), whereby the application of stimulation is locked to a particular phase of tremor, has been proposed as a strategy to improve therapeutic efficacy and decrease side effects. In this work, in the context of these phase-specific aDBS strategies, we investigate the dynamical behavior of large populations of coupled neurons in response to near-periodic stimulation, namely, stimulation that is periodic except for a slowly changing amplitude and phase offset that can be used to coordinate the timing of applied input with a specified phase of model oscillations. Using an adaptive phase-amplitude reduction strategy, we illustrate that for a large population of oscillatory neurons, the temporal evolution of the associated phase distribution in response to near-periodic forcing can be captured using a reduced order model with four state variables. Subsequently, we devise and validate a closed-loop control strategy to disrupt synchronization caused by coupling. Additionally, we identify strategies for implementing the proposed control strategy in situations where underlying model equations are unavailable by estimating the necessary terms of the reduced order equations in real-time from observables.

Asymptotic Phase and Amplitude for Classical and Semiclassical Stochastic Oscillators via Koopman Operator Theory

The asymptotic phase is a fundamental quantity for the analysis of deterministic limit-cycle oscillators, and generalized definitions of the asymptotic phase for stochastic oscillators have also been proposed. In this article, we show that the asymptotic phase and also amplitude can be defined for classical and semiclassical stochastic oscillators in a natural and unified manner by using the eigenfunctions of the Koopman operator of the system. We show that the proposed definition gives appropriate values of the phase and amplitude for strongly stochastic limit-cycle oscillators, excitable systems undergoing noise-induced oscillations, and also for quantum limit-cycle oscillators in the semiclassical regime.

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.

Exploiting circadian memory to hasten recovery from circadian misalignment

Recent years have seen a sustained interest in the development of circadian reentrainment strategies to limit the deleterious effects of jet lag. Due to the dynamical complexity of many circadian models, phase-based model reduction techniques are often an imperative first step in the analysis. However, amplitude coordinates that capture lingering effects (i.e., memory) from past inputs are often neglected. In this work, we focus on these amplitude coordinates using an operational phase and an isostable coordinate framework in the context of the development of jet-lag amelioration strategies. By accounting for the influence of circadian memory, we identify a latent phase shift that can prime one's circadian cycle to reentrain more rapidly to an expected time-zone shift. A subsequent optimal control problem is proposed that balances the trade-off between control effort and the resulting latent phase shift. Data-driven model identification techniques for the inference of necessary reduced order, phase-amplitude-based models are considered in situations where the underlying model equations are unknown, and numerical results are illustrated in both a simple planar model and in a coupled population of circadian oscillators.

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 L₁ 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 L₁-norm optimization with additional constraints on the individual connection weights.

Leveraging deep learning to control neural oscillators

Modulation of the firing times of neural oscillators has long been an important control objective, with applications including Parkinson's disease, Tourette's syndrome, epilepsy, and learning. One common goal for such modulation is desynchronization, wherein two or more oscillators are stimulated to transition from firing in phase with each other to firing out of phase. The optimization of such stimuli has been well studied, but this typically relies on either a reduction of the dimensionality of the system or complete knowledge of the parameters and state of the system. This limits the applicability of results to real problems in neural control. Here, we present a trained artificial neural network capable of accurately estimating the effects of square-wave stimuli on neurons using minimal output information from the neuron. We then apply the results of this network to solve several related control problems in desynchronization, including desynchronizing pairs of neurons and achieving clustered subpopulations of neurons in the presence of coupling and noise.

Degenerate isostable reduction for fixed-point and limit-cycle attractors with defective linearizations

Isostable coordinates provide a convenient framework for understanding the transient behavior of dynamical systems with stable attractors. These isostable coordinates are often used to characterize the slowest decaying eigenfunctions of the Koopman operator; by neglecting the rapidly decaying Koopman eigenfunctions a reduced order model can be obtained. Existing work has focused primarily on nondegenerate isostable coordinates, that is, isostable coordinates that are associated with eigenvalues that have identical algebraic and geometric multiplicities. Current isostable reduction methods cannot be applied to characterize the decay associated with a defective eigenvalue. In this work, a degenerate isostable framework is proposed for use when eigenvalues are defective. These degenerate isostable coordinates are investigated in the context of various reduced order modeling frameworks that retain many of the important properties of standard (nondegenerate) isostable reduced modeling strategies. Reduced order modeling examples that require the use of degenerate isostable coordinates are presented with relevance to both circadian physiology and nonlinear fluid flows.

Shape versus timing: linear responses of a limit cycle with hard boundaries under instantaneous and static perturbation

When dynamical systems that produce rhythmic behaviors operate within hard limits, they may exhibit limit cycles with sliding components, that is, closed isolated periodic orbits that make and break contact with a constraint surface. Examples include heel-ground interaction in locomotion, firing rate rectification in neural networks, and stick-slip oscillators. In many rhythmic systems, robustness against external perturbations involves response of both the shape and the timing of the limit cycle trajectory. The existing methods of infinitesimal phase response curve (iPRC) and variational analysis are well established for quantifying changes in timing and shape, respectively, for smooth systems. These tools have recently been extended to nonsmooth dynamics with transversal crossing boundaries. In this work, we further extend the iPRC method to nonsmooth systems with sliding components, which enables us to make predictions about the synchronization properties of weakly coupled stick-slip oscillators. We observe a new feature of the isochrons in a planar limit cycle with hard sliding boundaries: a nonsmooth kink in the asymptotic phase function, originating from the point at which the limit cycle smoothly departs the constraint surface, and propagating away from the hard boundary into the interior of the domain. Moreover, the classical variational analysis neglects timing information and is restricted to instantaneous perturbations. By defining the "infinitesimal shape response curve" (iSRC), we incorporate timing sensitivity of an oscillator to describe the shape response of this oscillator to parametric perturbations. In order to extract timing information, we also develop a "local timing response curve" (lTRC) that measures the timing sensitivity of a limit cycle within any given region. We demonstrate in a specific example that taking into account local timing sensitivity in a nonsmooth system greatly improves the accuracy of the iSRC over global timing analysis given by the iPRC.

Global phase-amplitude description of oscillatory dynamics via the parameterization method

In this paper, we use the parameterization method to provide a complete description of the dynamics of an n-dimensional oscillator beyond the classical phase reduction. The parameterization method allows us, via efficient algorithms, to obtain a parameterization of the attracting invariant manifold of the limit cycle in terms of the phase-amplitude variables. The method has several advantages. It provides analytically a Fourier-Taylor expansion of the parameterization up to any order, as well as a simplification of the dynamics that allows for a numerical globalization of the manifolds. Thus, one can obtain the local and global isochrons and isostables, including the slow attracting manifold, up to high accuracy, which offer a geometrical portrait of the oscillatory dynamics. Furthermore, it provides straightforwardly the infinitesimal phase and amplitude response functions, that is, the extended infinitesimal phase and amplitude response curves, which monitor the phase and amplitude shifts beyond the asymptotic state. Thus, the methodology presented yields an accurate description of the phase dynamics for perturbations not restricted to the limit cycle but to its attracting invariant manifold. Finally, we explore some strategies to reduce the dimension of the dynamics, including the reduction of the dynamics to the slow stable submanifold. We illustrate our methods by applying them to different three-dimensional single neuron and neural population models in neuroscience.

Stabilization of Weakly Unstable Fixed Points as a Common Dynamical Mechanism of High-Frequency Electrical Stimulation

While high-frequency electrical stimulation often used to treat various biological diseases, it is generally difficult to understand its dynamical mechanisms of action. In this work, high-frequency electrical stimulation is considered in the context of neurological and cardiological systems. Despite inherent differences between these systems, results from both theory and computational modeling suggest identical dynamical mechanisms responsible for desirable qualitative changes in behavior in response to high-frequency stimuli. Specifically, desynchronization observed in a population of periodically firing neurons and reversible conduction block that occurs in cardiomyocytes both result from bifurcations engendered by stimulation that modifies the stability of unstable fixed points. Using a reduced order phase-amplitude modeling framework, this phenomenon is described in detail from a theoretical perspective. Results are consistent with and provide additional insight for previously published experimental observations. Also, it is found that sinusoidal input is energy-optimal for modifying the stability of weakly unstable fixed points using periodic stimulation.

Phase-amplitude reduction far beyond the weakly perturbed paradigm

While phase reduction is a well-established technique for the analysis of perturbed limit cycle oscillators, practical application requires perturbations to be sufficiently weak thereby limiting its utility in many situations. Here, a general strategy is developed for constructing a set of phase-amplitude reduced equations that is valid to arbitrary orders of accuracy in the amplitude coordinates. This reduction framework can be used to investigate the behavior of oscillatory dynamical systems far beyond the weakly perturbed paradigm. Additionally, a patchwork phase-amplitude reduction method is suggested that is useful when exceedingly large magnitude perturbations are considered. This patchwork method incorporates the high-accuracy phase-amplitude reductions of multiple nearby periodic orbits that result from modifications to nominal parameters. The proposed method of high-accuracy phase-amplitude reduction can be readily implemented numerically and examples are provided where reductions are computed up to fourteenth order accuracy.

A data-driven phase and isostable reduced modeling framework for oscillatory dynamical systems

Phase-amplitude reduction is of growing interest as a strategy for the reduction and analysis of oscillatory dynamical systems. Augmentation of the widely studied phase reduction with amplitude coordinates can be used to characterize transient behavior in directions transverse to a limit cycle to give a richer description of the dynamical behavior. Various definitions for amplitude coordinates have been suggested, but none are particularly well suited for implementation in experimental systems where output recordings are readily available but the underlying equations are typically unknown. In this work, a reduction framework is developed for inferring a phase-amplitude reduced model using only the observed model output from an arbitrarily high-dimensional system. This framework employs a proper orthogonal reduction strategy to identify important features of the transient decay of solutions to the limit cycle. These features are explicitly related to previously developed phase and isostable coordinates and used to define so-called data-driven phase and isostable coordinates that are valid in the entire basin of attraction of a limit cycle. The utility of this reduction strategy is illustrated in examples related to neural physiology and is used to implement an optimal control strategy that would otherwise be computationally intractable. The proposed data-driven phase and isostable coordinate system and associated reduced modeling framework represent a useful tool for the study of nonlinear dynamical systems in situations where the underlying dynamical equations are unknown and in particularly high-dimensional or complicated numerical systems for which standard phase-amplitude reduction techniques are not computationally feasible.

Recent advances in coupled oscillator theory

We review the theory of weakly coupled oscillators for smooth systems. We then examine situations where application of the standard theory falls short and illustrate how it can be extended. Specific examples are given to non-smooth systems with applications to the Izhikevich neuron. We then introduce the idea of isostable reduction to explore behaviours that the weak coupling paradigm cannot explain. In an additional example, we show how bifurcations that change the stability of phase-locked solutions in a pair of identical coupled neurons can be understood using the notion of isostable reduction. This article is part of the theme issue 'Coupling functions: dynamical interaction mechanisms in the physical, biological and social sciences'.

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'.

Dynamic causal modelling of phase-amplitude interactions

Models of coupled phase oscillators are used to describe a wide variety of phenomena in neuroimaging. These models typically rest on the premise that oscillator dynamics do not evolve beyond their respective limit cycles, and hence that interactions can be described purely in terms of phase differences. Whilst mathematically convenient, the restrictive nature of phase-only models can limit their explanatory power. We therefore propose a generalisation of dynamic causal modelling that incorporates both phase and amplitude. This allows for the separate quantifications of phase and amplitude contributions to the connectivity between neural regions. We show, using model-generated data and simulations of coupled pendula, that phase-amplitude models can describe strongly coupled systems more effectively than their phase-only counterparts. We relate our findings to four metrics commonly used in neuroimaging: the Kuramoto order parameter, cross-correlation, phase-lag index, and spectral entropy. We find that, with the exception of spectral entropy, the phase-amplitude model is able to capture all metrics more effectively than the phase-only model. We then demonstrate, using local field potential recordings in rodents and functional magnetic resonance imaging in macaque monkeys, that amplitudes in oscillator models play an important role in describing neural dynamics in anaesthetised brain states.

Synchronization of Electrically Coupled Resonate-and-Fire Neurons

Electrical coupling between neurons is broadly present across brain areas and is typically assumed to synchronize network activity. However, intrinsic properties of the coupled cells can complicate this simple picture. Many cell types with electrical coupling show a diversity of post-spike subthreshold fluctuations, often linked to subthreshold resonance, which are transmitted through electrical synapses in addition to action potentials. Using the theory of weakly coupled oscillators, we explore the effect of both subthreshold and spike-mediated coupling on synchrony in small networks of electrically coupled resonate-and-fire neurons, a hybrid neuron model with damped subthreshold oscillations and a range of post-spike voltage dynamics. We calculate the phase response curve using an extension of the adjoint method that accounts for the discontinuous post-spike reset rule. We find that both spikes and subthreshold fluctuations can jointly promote synchronization. The subthreshold contribution is strongest when the voltage exhibits a significant post-spike elevation in voltage, or plateau potential. Additionally, we show that the geometry of trajectories approaching the spiking threshold causes a "reset-induced shear" effect that can oppose synchrony in the presence of network asymmetry, despite having no effect on the phase-locking of symmetrically coupled pairs.

Optimal phase control of biological oscillators using augmented phase reduction

We develop a novel optimal control algorithm to change the phase of an oscillator using a minimum energy input, which also minimizes the oscillator's transversal distance to the uncontrolled periodic orbit. Our algorithm uses a two-dimensional reduction technique based on both isochrons and isostables. We develop a novel method to eliminate cardiac alternans by connecting our control algorithm with the underlying physiological problem. We also describe how the devised algorithm can be used for spike timing control which can potentially help with motor symptoms of essential and parkinsonian tremor, and aid in treating jet lag. To demonstrate the advantages of this algorithm, we compare it with a previously proposed optimal control algorithm based on standard phase reduction for the Hopf bifurcation normal form, and models for cardiac pacemaker cells, thalamic neurons, and circadian gene regulation cycle in the suprachiasmatic nucleus. We show that our control algorithm is effective even when a large phase change is required or when the nontrivial Floquet multiplier is close to unity; in such cases, the previously proposed control algorithm fails.

Phase reduction and phase-based optimal control for biological systems: a tutorial

A powerful technique for the analysis of nonlinear oscillators is the rigorous reduction to phase models, with a single variable describing the phase of the oscillation with respect to some reference state. An analog to phase reduction has recently been proposed for systems with a stable fixed point, and phase reduction for periodic orbits has recently been extended to take into account transverse directions and higher-order terms. This tutorial gives a unified treatment of such phase reduction techniques and illustrates their use through mathematical and biological examples. It also covers the use of phase reduction for designing control algorithms which optimally change properties of the system, such as the phase of the oscillation. The control techniques are illustrated for example neural and cardiac systems.

Isostable reduction of oscillators with piecewise smooth dynamics and complex Floquet multipliers

Phase-amplitude reduction is a widely applied technique in the study of limit cycle oscillators with the ability to represent a complicated and high-dimensional dynamical system in a more analytically tractable set of coordinates. Recent work has focused on the use of isostable coordinates, which characterize the transient decay of solutions toward a periodic orbit, and can ultimately be used to increase the accuracy of these reduced models. The breadth of systems to which this phase-amplitude reduction strategy can be applied, however, is still rather limited. In this work, the theory of phase-amplitude reduction using isostable coordinates is further developed to accommodate a broader set of dynamical systems. In the first part, limit cycles of piecewise smooth dynamical systems are considered and strategies are developed to compute the associated reduced equations. In the second part, the notion of isostable coordinates for complex-valued Floquet multipliers is introduced, resulting in one phaselike coordinate and one amplitudelike coordinate for each pair of complex conjugate Floquet multipliers. Examples are given with relevance to piecewise smooth representations of excitable cardiomyocytes and the relationship between the reduced coordinate system and the emergence of cardiac alternans is discussed. Also, phase-amplitude reduction is implemented for a chaotic, externally forced pendulum with complex Floquet multipliers and a resulting control strategy for the stabilization of its periodic solution is investigated.

Global computation of phase-amplitude reduction for limit-cycle dynamics

Recent years have witnessed increasing interest in phase-amplitude reduction of limit-cycle dynamics. Adding an amplitude coordinate to the phase coordinate allows us to take into account the dynamics transversal to the limit cycle and thereby overcome the main limitations of classic phase reduction (strong convergence to the limit cycle and weak inputs). While previous studies, mostly focus on local quantities such as infinitesimal responses, a major and limiting challenge of phase-amplitude reduction is to compute amplitude coordinates globally, in the basin of attraction of the limit cycle. In this paper, we propose a method to compute the full set of phase-amplitude coordinates in the large. Our method is based on the so-called Koopman (composition) operator and aims at computing the eigenfunctions of the operator through Laplace averages (in combination with the harmonic balance method). This yields a forward integration method that is not limited to two-dimensional systems. We illustrate the method by computing the so-called isostables of limit cycles in two-, three-, and four-dimensional state spaces, as well as their responses to strong external inputs.

Greater accuracy and broadened applicability of phase reduction using isostable coordinates

The applicability of phase models is generally limited by the constraint that the dynamics of a perturbed oscillator must stay near its underlying periodic orbit. Consequently, external perturbations must be sufficiently weak so that these assumptions remain valid. Using the notion of isostables of periodic orbits to provide a simplified coordinate system from which to understand the dynamics transverse to a periodic orbit, we devise a strategy to correct for changing phase dynamics for locations away from the limit cycle. Consequently, these corrected phase dynamics allow for perturbations of larger magnitude without invalidating the underlying assumptions of the reduction. The proposed reduction strategy yields a closed set of equations and can be applied to periodic orbits embedded in arbitrarily high dimensional spaces. We illustrate the utility of this strategy in two models with biological relevance. In the first application, we find that an optimal control strategy for modifying the period of oscillation can be improved with the corrected phase reduction. In the second, the corrected phase reduced dynamics are used to understand adaptation and memory effects resulting from past perturbations.

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.