Adjoint method provides phase response functions for delay-induced oscillations

Nonreciprocal synchronization in embryonic oscillator ensembles

Synchronization of coupled oscillators is a universal phenomenon encountered across different scales and contexts, e.g., chemical wave patterns, superconductors, and the unison applause we witness in concert halls. The existence of common underlying coupling rules defines universality classes, revealing a fundamental sameness between seemingly distinct systems. Identifying rules of synchronization in any particular setting is hence of paramount relevance. Here, we address the coupling rules within an embryonic oscillator ensemble linked to vertebrate embryo body axis segmentation. In vertebrates, the periodic segmentation of the body axis involves synchronized signaling oscillations in cells within the presomitic mesoderm (PSM), from which somites, the prevertebrae, form. At the molecular level, it is known that intact Notch-signaling and cell-to-cell contact are required for synchronization between PSM cells. However, an understanding of the coupling rules is still lacking. To identify these, we develop an experimental assay that enables direct quantification of synchronization dynamics within mixtures of oscillating cell ensembles, for which the initial input frequency and phase distribution are known. Our results reveal a "winner-takes-it-all" synchronization outcome, i.e., the emerging collective rhythm matches one of the input rhythms. Using a combination of theory and experimental validation, we develop a coupling model, the "Rectified Kuramoto" (ReKu) model, characterized by a phase-dependent, nonreciprocal interaction in the coupling of oscillatory cells. Such nonreciprocal synchronization rules reveal fundamental similarities between embryonic oscillators and a class of collective behaviors seen in neurons and fireflies, where higher-level computations are performed and linked to nonreciprocal synchronization.

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.

Reconstruction of phase dynamics from macroscopic observations based on linear and nonlinear response theories

We propose a method to reconstruct the phase dynamics in rhythmical interacting systems from macroscopic responses to weak inputs by developing linear and nonlinear response theories, which predict the responses in a given system. By solving an inverse problem, the method infers an unknown system: the natural frequency distribution, the coupling function, and the time delay which is inevitable in real systems. In contrast to previous methods, our method requires neither strong invasiveness nor microscopic observations. We demonstrate that the method reconstructs two phase systems from observed responses accurately. The qualitative methodological advantages demonstrated by our quantitative numerical examinations suggest its broad applicability in various fields, including brain systems, which are often observed through macroscopic signals such as electroencephalograms and functional magnetic response imaging.

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.

Controlling fluidic oscillator flow dynamics by elastic structure vibration

In this study, we introduce a design of a feedback-type fluidic oscillator with elastic structures surrounding its feedback channel. By employing phase reduction theory, we extract the phase sensitivity function of the complex fluid-structure coupled system, which represents the system's oscillatory characteristics. We show that the frequency of the oscillating flow inside the fluidic oscillator can be modulated by inducing synchronization with the weak periodic forcing from the elastic structure vibration. This design approach adds controllability to the fluidic oscillator, where conventionally, the intrinsic oscillatory characteristics of such device were highly determined by its geometry. The synchronization-induced control also changes the physical characteristics of the oscillatory fluid flow, which can be beneficial for practical applications, such as promoting better fluid mixing without changing the overall geometry of the device. Furthermore, by analyzing the phase sensitivity function, we demonstrate how the use of phase reduction theory gives good estimation of the synchronization condition with minimal number of experiments, allowing for a more efficient control design process. Finally, we show how an optimal control signal can be designed to reach the fastest time to synchronization.

Arnold tongue entrainment reveals dynamical principles of the embryonic segmentation clock

Living systems exhibit an unmatched complexity, due to countless, entangled interactions across scales. Here, we aim to understand a complex system, that is, segmentation timing in mouse embryos, without a reference to these detailed interactions. To this end, we develop a coarse-grained approach, in which theory guides the experimental identification of the segmentation clock entrainment responses. We demonstrate period- and phase-locking of the segmentation clock across a wide range of entrainment parameters, including higher-order coupling. These quantifications allow to derive the phase response curve (PRC) and Arnold tongues of the segmentation clock, revealing its essential dynamical properties. Our results indicate that the somite segmentation clock has characteristics reminiscent of a highly non-linear oscillator close to an infinite period bifurcation and suggests the presence of long-term feedbacks. Combined, this coarse-grained theoretical-experimental approach reveals how we can derive simple, essential features of a highly complex dynamical system, providing precise experimental control over the pace and rhythm of the somite segmentation clock.

Phase response approaches to neural activity models with distributed delay

In weakly coupled neural oscillator networks describing brain dynamics, the coupling delay is often distributed. We present a theoretical framework to calculate the phase response curve of distributed-delay induced limit cycles with infinite-dimensional phase space. Extending previous works, in which non-delayed or discrete-delay systems were investigated, we develop analytical results for phase response curves of oscillatory systems with distributed delay using Gaussian and log-normal delay distributions. We determine the scalar product and normalization condition for the linearized adjoint of the system required for the calculation of the phase response curve. As a paradigmatic example, we apply our technique to the Wilson-Cowan oscillator model of excitatory and inhibitory neuronal populations under the two delay distributions. We calculate and compare the phase response curves for the Gaussian and log-normal delay distributions. The phase response curves obtained from our adjoint calculations match those compiled by the direct perturbation method, thereby proving that the theory of weakly coupled oscillators can be applied successfully for distributed-delay-induced limit cycles. We further use the obtained phase response curves to derive phase interaction functions and determine the possible phase locked states of multiple inter-coupled populations to illuminate different synchronization scenarios. In numerical simulations, we show that the coupling delay distribution can impact the stability of the synchronization between inter-coupled gamma-oscillatory networks.

Unstable delayed feedback control to change sign of coupling strength for weakly coupled limit cycle oscillators

Weakly coupled limit cycle oscillators can be reduced into a system of weakly coupled phase models. These phase models are helpful to analyze the synchronization phenomena. For example, a phase model of two oscillators has a one-dimensional differential equation for the evolution of the phase difference. The existence of fixed points determines frequency-locking solutions. By treating each oscillator as a black-box possessing a single input and a single output, one can investigate various control algorithms to change the synchronization of the oscillators. In particular, we are interested in a delayed feedback control algorithm. Application of this algorithm to the oscillators after a subsequent phase reduction should give the same phase model as in the control-free case, but with a rescaled coupling strength. The conventional delayed feedback control is limited to the change of magnitude but does not allow the change of sign of the coupling strength. In this work, we present a modification of the delayed feedback algorithm supplemented by an additional unstable degree of freedom, which is able to change the sign of the coupling strength. Various numerical calculations performed with Landau-Stuart and FitzHugh-Nagumo oscillators show successful switching between an in-phase and anti-phase synchronization using the provided control algorithm. Additionally, we show that the control force becomes non-invasive if our objective is stabilization of an unstable phase difference for two coupled oscillators.

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.

Analytical approximations for the speed of pacemaker-generated waves

In oscillatory media, waves can be generated by pacemaker regions which oscillate faster than their surroundings. In many chemical and biological systems, such waves can synchronize the whole medium and as such they are a means of transmitting information at a fixed speed over large distances. In this paper, we apply analytical tools to investigate the factors that determine the speed of these waves. More precisely, we apply singular perturbation and phase reduction methods to two types of negative-feedback oscillators, one built on underlying bistability and one including a time delay in the negative feedback. In both systems, we investigate the influence of time-scale separation on the resulting wave speed, as well as the effect of size and frequency of the pacemaker region. We compare our analytical estimates to numerical simulations which we described previously [J. Rombouts and L. Gelens, Phys. Rev. Research 2, 043038 (2020)2643-156410.1103/PhysRevResearch.2.043038].

From local resynchronization to global pattern recovery in the zebrafish segmentation clock

Integrity of rhythmic spatial gene expression patterns in the vertebrate segmentation clock requires local synchronization between neighboring cells by Delta-Notch signaling and its inhibition causes defective segment boundaries. Whether deformation of the oscillating tissue complements local synchronization during patterning and segment formation is not understood. We combine theory and experiment to investigate this question in the zebrafish segmentation clock. We remove a Notch inhibitor, allowing resynchronization, and analyze embryonic segment recovery. We observe unexpected intermingling of normal and defective segments, and capture this with a new model combining coupled oscillators and tissue mechanics. Intermingled segments are explained in the theory by advection of persistent phase vortices of oscillators. Experimentally observed changes in recovery patterns are predicted in the theory by temporal changes in tissue length and cell advection pattern. Thus, segmental pattern recovery occurs at two length and time scales: rapid local synchronization between neighboring cells, and the slower transport of the resulting patterns across the tissue through morphogenesis.

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

Optimization of linear and nonlinear interaction schemes for stable synchronization of weakly coupled limit-cycle oscillators

Optimization of mutual synchronization between a pair of limit-cycle oscillators with weak symmetric coupling is considered in the framework of the phase-reduction theory. By generalizing our previous study [S. Shirasaka, N. Watanabe, Y. Kawamura, and H. Nakao, Optimizing stability of mutual synchronization between a pair of limit-cycle oscillators with weak cross coupling, Phys. Rev. E 96, 012223 (2017)2470-004510.1103/PhysRevE.96.012223] on the optimization of cross-diffusion coupling matrices between the oscillators, we consider optimization of mutual coupling signals to maximize the linear stability of the synchronized state, which are functionals of the past time sequences of the oscillator states. For the case of linear coupling, optimization of the delay time and linear filtering of coupling signals are considered. For the case of nonlinear coupling, general drive-response coupling is considered and the optimal response and driving functions are derived. The theoretical results are illustrated by numerical simulations.

Macroscopic phase resetting-curves determine oscillatory coherence and signal transfer in inter-coupled neural circuits

Macroscopic oscillations of different brain regions show multiple phase relationships that are persistent across time and have been implicated in routing information. While multiple cellular mechanisms influence the network oscillatory dynamics and structure the macroscopic firing motifs, one of the key questions is to identify the biophysical neuronal and synaptic properties that permit such motifs to arise. A second important issue is how the different neural activity coherence states determine the communication between the neural circuits. Here we analyse the emergence of phase-locking within bidirectionally delayed-coupled spiking circuits in which global gamma band oscillations arise from synaptic coupling among largely excitable neurons. We consider both the interneuronal (ING) and the pyramidal-interneuronal (PING) population gamma rhythms and the inter coupling targeting the pyramidal or the inhibitory neurons. Using a mean-field approach together with an exact reduction method, we reduce each spiking network to a low dimensional nonlinear system and derive the macroscopic phase resetting-curves (mPRCs) that determine how the phase of the global oscillation responds to incoming perturbations. This is made possible by the use of the quadratic integrate-and-fire model together with a Lorentzian distribution of the bias current. Depending on the type of gamma (PING vs. ING), we show that incoming excitatory inputs can either speed up the macroscopic oscillation (phase advance; type I PRC) or induce both a phase advance and a delay (type II PRC). From there we determine the structure of macroscopic coherence states (phase-locking) of two weakly synaptically-coupled networks. To do so we derive a phase equation for the coupled system which links the synaptic mechanisms to the coherence states of the system. We show that a synaptic transmission delay is a necessary condition for symmetry breaking, i.e. a non-symmetric phase lag between the macroscopic oscillations. This potentially provides an explanation to the experimentally observed variety of gamma phase-locking modes. Our analysis further shows that symmetry-broken coherence states can lead to a preferred direction of signal transfer between the oscillatory networks where this directionality also depends on the timing of the signal. Hence we suggest a causal theory for oscillatory modulation of functional connectivity between cortical circuits.

Large scale brain oscillations emerge from synaptic interactions within neuronal circuits. Over the past years, such macroscopic rhythms have been suggested to play a crucial role in routing the flow of information across cortical regions, resulting in a functional connectome. The underlying mechanism is cortical oscillations that bind together following a well-known motif called phase-locking. While there is significant experimental support for multiple phase-locking modes in the brain, it is still unclear what is the underlying mechanism that permits macroscopic rhythms to phase lock. In the present paper we take up with this issue, and to show that, one can study the emergent macroscopic phase-locking within the mathematical framework of weakly coupled oscillators. We find that under synaptic delays, fully symmetrically coupled networks can display symmetry-broken states of activity, where one network starts to lead in phase the second (also sometimes known as stuttering states). When we analyse how incoming transient signals affect the coupled system, we find that in the symmetry-broken state, the effect depends strongly on which network is targeted (the leader or the follower) as well as the timing of the input. Hence we show how the dynamics of the emergent phase-locked activity imposes a functional directionality on how signals are processed. We thus offer clarification on the synaptic and circuit properties responsible for the emergence of multiple phase-locking patterns and provide support for its functional implication in information transfer.

Response of an oscillatory differential delay equation to a periodic stimulus

Periodic hematological diseases such as cyclical neutropenia or cyclical thrombocytopenia, with their characteristic oscillations of circulating neutrophils or platelets, may pose grave problems for patients. Likewise, periodically administered chemotherapy has the unintended side effect of establishing periodic fluctuations in circulating white cells, red cell precursors and/or platelets. These fluctuations, either spontaneous or induced, often have serious consequences for the patient (e.g. neutropenia, anemia, or thrombocytopenia respectively) which exogenously administered cytokines can partially correct. The question of when and how to administer these drugs is a difficult one for clinicians and not easily answered. In this paper we use a simple model consisting of a delay differential equation with a piecewise linear nonlinearity, that has a periodic solution, to model the effect of a periodic disease or periodic chemotherapy. We then examine the response of this toy model to both single and periodic perturbations, meant to mimic the drug administration, as a function of the drug dose and the duration and frequency of its administration to best determine how to avoid side effects.

Optimizing mutual synchronization of rhythmic spatiotemporal patterns in reaction-diffusion systems

Optimization of the stability of synchronized states between a pair of symmetrically coupled reaction-diffusion systems exhibiting rhythmic spatiotemporal patterns is studied in the framework of the phase reduction theory. The optimal linear filter that maximizes the linear stability of the in-phase synchronized state is derived for the case in which the two systems are nonlocally coupled. The optimal nonlinear interaction function that theoretically gives the largest linear stability of the in-phase synchronized state is also derived. The theory is illustrated by using typical rhythmic patterns in FitzHugh-Nagumo systems as examples.

Macroscopic phase-resetting curves for spiking neural networks

The study of brain rhythms is an open-ended, and challenging, subject of interest in neuroscience. One of the best tools for the understanding of oscillations at the single neuron level is the phase-resetting curve (PRC). Synchronization in networks of neurons, effects of noise on the rhythms, effects of transient stimuli on the ongoing rhythmic activity, and many other features can be understood by the PRC. However, most macroscopic brain rhythms are generated by large populations of neurons, and so far it has been unclear how the PRC formulation can be extended to these more common rhythms. In this paper, we describe a framework to determine a macroscopic PRC (mPRC) for a network of spiking excitatory and inhibitory neurons that generate a macroscopic rhythm. We take advantage of a thermodynamic approach combined with a reduction method to simplify the network description to a small number of ordinary differential equations. From this simplified but exact reduction, we can compute the mPRC via the standard adjoint method. Our theoretical findings are illustrated with and supported by numerical simulations of the full spiking network. Notably our mPRC framework allows us to predict the difference between effects of transient inputs to the excitatory versus the inhibitory neurons in the network.

Optimizing stability of mutual synchronization between a pair of limit-cycle oscillators with weak cross coupling

We consider optimization of the linear stability of synchronized states between a pair of weakly coupled limit-cycle oscillators with cross coupling, where different components of state variables of the oscillators are allowed to interact. On the basis of the phase reduction theory, we derive the coupling matrix between different components of the oscillator states that maximizes the linear stability of the synchronized state under given constraints on the overall coupling intensity and the stationary phase difference. The improvement in the linear stability is illustrated by using several types of limit-cycle oscillators as examples.

Critical dynamics of Hopf bifurcations in the corticothalamic system: Transitions from normal arousal states to epileptic seizures

A physiologically based corticothalamic model of large-scale brain activity is used to analyze critical dynamics of transitions from normal arousal states to epileptic seizures, which correspond to Hopf bifurcations. This relates an abstract normal form quantitatively to underlying physiology that includes neural dynamics, axonal propagation, and time delays. Thus, a bridge is constructed that enables normal forms to be used to interpret quantitative data. The normal form of the Hopf bifurcations with delays is derived using Hale's theory, the center manifold theorem, and normal form analysis, and it is found to be explicitly expressed in terms of transfer functions and the sensitivity matrix of a reduced open-loop system. It can be applied to understand the effect of each physiological parameter on the critical dynamics and determine whether the Hopf bifurcation is supercritical or subcritical in instabilities that lead to absence and tonic-clonic seizures. Furthermore, the effects of thalamic and cortical nonlinearities on the bifurcation type are investigated, with implications for the roles of underlying physiology. The theoretical predictions about the bifurcation type and the onset dynamics are confirmed by numerical simulations and provide physiologically based criteria for determining bifurcation types from first principles. The results are consistent with experimental data from previous studies, imply that new regimes of seizure transitions may exist in clinical settings, and provide a simplified basis for control-systems interventions. Using the normal form, and the full equations from which it is derived, more complex dynamics, such as quasiperiodic cycles and saddle cycles, are discovered near the critical points of the subcritical Hopf bifurcations.

Deriving theoretical phase locking values of a coupled cortico-thalamic neural mass model using center manifold reduction

Cognitive functions such as sensory processing and memory processes lead to phase synchronization in the electroencephalogram or local field potential between different brain regions. There are a lot of computational researches deriving phase locking values (PLVs), which are an index of phase synchronization intensity, from neural models. However, these researches derive PLVs numerically. To the best of our knowledge, there have been no reports on the derivation of a theoretical PLV. In this study, we propose an analytical method for deriving theoretical PLVs from a cortico-thalamic neural mass model described by a delay differential equation. First, the model for generating neural signals is transformed into a normal form of the Hopf bifurcation using center manifold reduction. Second, the normal form is transformed into a phase model that is suitable for analyzing synchronization phenomena. Third, the Fokker-Planck equation of the phase model is derived and the phase difference distribution is obtained. Finally, the PLVs are calculated from the stationary distribution of the phase difference. The validity of the proposed method is confirmed via numerical simulations. Furthermore, we apply the proposed method to a working memory process, and discuss the neurophysiological basis behind the phase synchronization phenomenon. The results demonstrate the importance of decreasing the intensity of independent noise during the working memory process. The proposed method will be of great use in various experimental studies and simulations relevant to phase synchronization, because it enables the effect of neurophysiological changes on PLVs to be analyzed from a mathematical perspective.

Phase reduction theory for hybrid nonlinear oscillators

Hybrid dynamical systems characterized by discrete switching of smooth dynamics have been used to model various rhythmic phenomena. However, the phase reduction theory, a fundamental framework for analyzing the synchronization of limit-cycle oscillations in rhythmic systems, has mostly been restricted to smooth dynamical systems. Here we develop a general phase reduction theory for weakly perturbed limit cycles in hybrid dynamical systems that facilitates analysis, control, and optimization of nonlinear oscillators whose smooth models are unavailable or intractable. On the basis of the generalized theory, we analyze injection locking of hybrid limit-cycle oscillators by periodic forcing and reveal their characteristic synchronization properties, such as ultrafast and robust entrainment to the periodic forcing and logarithmic scaling at the synchronization transition. We also illustrate the theory by analyzing the synchronization dynamics of a simple physical model of biped locomotion.

Mathematical Frameworks for Oscillatory Network Dynamics in Neuroscience

The tools of weakly coupled phase oscillator theory have had a profound impact on the neuroscience community, providing insight into a variety of network behaviours ranging from central pattern generation to synchronisation, as well as predicting novel network states such as chimeras. However, there are many instances where this theory is expected to break down, say in the presence of strong coupling, or must be carefully interpreted, as in the presence of stochastic forcing. There are also surprises in the dynamical complexity of the attractors that can robustly appear-for example, heteroclinic network attractors. In this review we present a set of mathematical tools that are suitable for addressing the dynamics of oscillatory neural networks, broadening from a standard phase oscillator perspective to provide a practical framework for further successful applications of mathematics to understanding network dynamics in neuroscience.

Response of an oscillatory differential delay equation to a single stimulus

Here we analytically examine the response of a limit cycle solution to a simple differential delay equation to a single pulse perturbation of the piecewise linear nonlinearity. We construct the unperturbed limit cycle analytically, and are able to completely characterize the perturbed response to a pulse of positive amplitude and duration with onset at different points in the limit cycle. We determine the perturbed minima and maxima and period of the limit cycle and show how the pulse modifies these from the unperturbed case.

Delayed feedback control of synchronization in weakly coupled oscillator networks

We study control of synchronization in weakly coupled oscillator networks by using a phase-reduction approach. Starting from a general class of limit-cycle oscillators we derive a phase model, which shows that delayed feedback control changes effective coupling strengths and effective frequencies. We derive the analytical condition for critical control gain, where the phase dynamics of the oscillator becomes extremely sensitive to any perturbations. As a result the network can attain phase synchronization even if the natural interoscillatory couplings are small. In addition, we demonstrate that delayed feedback control can disrupt the coherent phase dynamic in synchronized networks. The validity of our results is illustrated on networks of diffusively coupled Stuart-Landau and FitzHugh-Nagumo models.

Phase reduction of a limit cycle oscillator perturbed by a strong amplitude-modulated high-frequency force

The phase reduction method for a limit cycle oscillator subjected to a strong amplitude-modulated high-frequency force is developed. An equation for the phase dynamics is derived by introducing a new, effective phase response curve. We show that if the effective phase response curve is everywhere positive (negative), then an entrainment of the oscillator to an envelope frequency is possible only when this frequency is higher (lower) than the natural frequency of the oscillator. Also, by using the Pontryagin maximum principle, we have derived an optimal waveform of the perturbation that ensures an entrainment of the oscillator with minimal power. The theoretical results are demonstrated with the Stuart-Landau oscillator and model neurons.

Delayed feedback control and phase reduction of unstable quasi-periodic orbits

The delayed feedback control (DFC) is applied to stabilize unstable quasi-periodic orbits (QPOs) in discrete-time systems. The feedback input is given by the difference between the current state and a time-delayed state in the DFC. However, there is an inevitable time-delay mismatch in QPOs. To evaluate the influence of the time-delay mismatch on the DFC, we propose a phase reduction method for QPOs and construct a phase response curve (PRC) from unstable QPOs directly. Using the PRC, we estimate the rotation number of QPO stabilized by the DFC. We show that the orbit of the DFC is consistent with the unstable QPO perturbed by a small state difference resulting from the time-delay mismatch, implying that the DFC can certainly stabilize the unstable QPO.

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.

Time-delayed feedback control design beyond the odd-number limitation

We present an algorithm for a time-delayed feedback control design to stabilize periodic orbits with an odd number of positive Floquet exponents in autonomous systems. Due to the so-called odd-number theorem such orbits have been considered as uncontrollable by time-delayed feedback methods. However, this theorem has been refuted by a counterexample and recently a corrected version of the theorem has been proved. In our algorithm, the control matrix is designed using a relationship between Floquet multipliers of the systems controlled by time-delayed and proportional feedback. The efficacy of the algorithm is demonstrated with the Lorenz and Chua systems.