The phase sensitivity curve or phase response curve (PRC) quantifies the oscillator's reaction to stimulation at a specific phase and is a primary characteristic of a self-sustained oscillatory unit. Knowledge of this curve yields a phase dynamics description of the oscillator for arbitrary weak forcing. Similar, though much less studied characteristic, is the amplitude response that can be defined either using an ad hoc approach to amplitude estimation or via the isostable variables. Here, we discuss the problem of the phase and amplitude response inference from observations using test stimulation. Although PRC determination for noise-free neuronal-like oscillators perturbed by narrow pulses is a well-known task, the general case remains a challenging problem. Even more challenging is the inference of the amplitude response. This characteristic is crucial, e.g., for controlling the amplitude of the collective mode in a network of interacting units; this task is relevant for neuroscience. Here, we compare the performance of different techniques suitable for inferring the phase and amplitude response, particularly with application to macroscopic oscillators. We suggest improvements to these techniques, e.g., demonstrating how to obtain the PRC in case of stimuli of arbitrary shape. Our main result is a novel technique denoted by IPID-1, based on the direct reconstruction of the Winfree equation and the analogous first-order equation for isostable dynamics. The technique works for signals with or without well-pronounced marker events and pulses of arbitrary shape; in particular, we consider charge-balanced pulses typical in neuroscience applications. Moreover, this technique is superior for noisy and high-dimensional systems. Additionally, we describe an error measure that can be computed solely from data and complements any inference technique.

Synchronization of two or more self-sustained oscillators is a well-known and studied phenomenon, appearing both in natural and designed systems. In some cases, the synchronized state is undesired, and the aim is to destroy synchrony by external intervention. In this paper, we focus on desynchronizing two self-sustained oscillators by short pulses delivered to the system in a phase-specific manner. We analyze a non-trivial case when we cannot access both oscillators but stimulate only one. The following restriction is that we can monitor only one unit, be it a stimulated or non-stimulated one. First, we use a system of two coupled Rayleigh oscillators to demonstrate how a loss of synchrony can be induced by stimulating a unit once per period at a specific phase and detected by observing consecutive inter-pulse durations. Next, we exploit the phase approximation to develop a rigorous theory formulating the problem in terms of a map. We derive exact expressions for the phase -isostable coordinates of this coupled system and show a relation between the phase and isostable response curves to the phase response curve of the uncoupled oscillator. Finally, we demonstrate how to obtain phase response information from the system using time series and discuss the differences between observing the stimulated and unstimulated oscillator.

Synchronization of two or more self-sustained oscillators is a well-known and studied phenomenon, appearing both in natural and designed systems. In some cases, the synchronized state is undesired, and the aim is to destroy synchrony by external intervention. In this paper, we focus on desynchronizing two self-sustained oscillators by short pulses delivered to the system in a phase-specific manner. We analyze a non-trivial case when we cannot access both oscillators but stimulate only one. The following restriction is that we can monitor only one unit, be it a stimulated or non-stimulated one. First, we use a system of two coupled Rayleigh oscillators to demonstrate how a loss of synchrony can be induced by stimulating a unit once per period at a specific phase and detected by observing consecutive inter-pulse durations. Next, we exploit the phase approximation to develop a rigorous theory formulating the problem in terms of a map. We derive exact expressions for the phase–isostable coordinates of this coupled system and show a relation between the phase and isostable response curves to the phase response curve of the uncoupled oscillator. Finally, we demonstrate how to obtain phase response information from the system using time series and discuss the differences between observing the stimulated and unstimulated oscillator.

The phase sensitivity curve or phase response curve (PRC) quantifies the oscillator's reaction to stimulation at a specific phase and is a primary characteristic of a self-sustained oscillatory unit. Knowledge of this curve yields a phase dynamics description of the oscillator for arbitrary weak forcing. Similar, though much less studied characteristic, is the amplitude response that can be defined either using an ad hoc approach to amplitude estimation or via the isostable variables. Here, we discuss the problem of the phase and amplitude response inference from observations using test stimulation. Although PRC determination for noise-free neuronal-like oscillators perturbed by narrow pulses is a well-known task, the general case remains a challenging problem. Even more challenging is the inference of the amplitude response. This characteristic is crucial, e.g., for controlling the amplitude of the collective mode in a network of interacting units; this task is relevant for neuroscience. Here, we compare the performance of different techniques suitable for inferring the phase and amplitude response, particularly with application to macroscopic oscillators. We suggest improvements to these techniques, e.g., demonstrating how to obtain the PRC in case of stimuli of arbitrary shape. Our main result is a novel technique denoted by IPID-1, based on the direct reconstruction of the Winfree equation and the analogous first-order equation for isostable dynamics. The technique works for signals with or without well-pronounced marker events and pulses of arbitrary shape; in particular, we consider charge-balanced pulses typical in neuroscience applications. Moreover, this technique is superior for noisy and high-dimensional systems. Additionally, we describe an error measure that can be computed solely from data and complements any inference technique.

Collective synchronization in a large population of self-sustained units appears both in natural and engineered systems. Sometimes this effect is in demand, while in some cases, it is undesirable, which calls for control techniques. In this paper, we focus on pulsatile control, with the goal to either increase or decrease the level of synchrony. We quantify this level by the entropy of the phase distribution. Motivated by possible applications in neuroscience, we consider pulses of a realistic shape. Exploiting the noisy Kuramoto–Winfree model, we search for the optimal pulse profile and the optimal stimulation phase. For this purpose, we derive an expression for the change of the phase distribution entropy due to the stimulus. We relate this change to the properties of individual units characterized by generally different natural frequencies and phase response curves and the population’s state. We verify the general result by analyzing a two-frequency population model and demonstrating a good agreement of the theory and numerical simulations.