An official website of the United States government
Phase reduction is a well-established method to study weakly driven and weakly perturbed oscillators. Traditional phase-reduction approaches characterize the perturbed system dynamics solely in terms of the timing of the oscillations. In the case of large perturbations, the introduction of amplitude (isostable) coordinates improves the accuracy of the phase description by providing a sense of distance from the underlying limit cycle. Importantly, phase-amplitude coordinates allow for the study of both the timing and shape of system oscillations. A parallel tool is the infinitesimal shape response curve (iSRC), a variational method that characterizes the shape change of a limit-cycle oscillator under sustained perturbation. Despite the importance of oscillation amplitude in a wide range of physical systems, systematic studies on the shape change of oscillations remain scarce. Both phase-amplitude coordinates and the iSRC represent methods to analyze oscillation shape change, yet a relationship between the two has not been previously explored. In this work, we establish the iSRC and phase-amplitude coordinates as complementary tools to study oscillation amplitude. We extend existing iSRC theory and specify conditions under which a general class of systems can be analyzed by the joint iSRC phase-amplitude approach. We show that the iSRC takes on a dramatically simple form in phase-amplitude coordinates, and directly relate the phase and isostable response curves to the iSRC. We apply our theory to weakly perturbed single oscillators, and to study the synchronization and entrainment of coupled oscillators.
more »
« less
Quantitative comparison of the mean–return-time phase and the stochastic asymptotic phase for noisy oscillators
Abstract Seminal work by A. Winfree and J. Guckenheimer showed that a deterministic phase variable can be defined either in terms of Poincaré sections or in terms of the asymptotic (long-time) behaviour of trajectories approaching a stable limit cycle. However, this equivalence between the deterministic notions of phase is broken in the presence of noise. Different notions of phase reduction for a stochastic oscillator can be defined either in terms of mean–return-time sections or as the argument of the slowest decaying complex eigenfunction of the Kolmogorov backwards operator. Although both notions of phase enjoy a solid theoretical foundation, their relationship remains unexplored. Here, we quantitatively compare both notions of stochastic phase. We derive an expression relating both notions of phase and use it to discuss differences (and similarities) between both definitions of stochastic phase for (i) a spiral sink motivated by stochastic models for electroencephalograms, (ii) noisy limit-cycle systems-neuroscience models, and (iii) a stochastic heteroclinic oscillator inspired by a simple motor-control system.
more »
« less
- Award ID(s):
- 2052109
- PAR ID:
- 10351835
- Date Published:
- Journal Name:
- Biological Cybernetics
- Volume:
- 116
- Issue:
- 2
- ISSN:
- 1432-0770
- Page Range / eLocation ID:
- 219 to 234
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)Abstract The characteristics of El-Niño-Southern Oscillation (ENSO) phase-locking in observations and CMIP5 and CMIP6 models are examined in this study. Two metrics based on the peaking month histogram for all El Niño and La Niña events are adopted to delineate the basic features of ENSO phase-locking in terms of the preferred calendar month and strength of this preference. It turns out that most models are poor at simulating the ENSO phase-locking, either showing little peak strengths or peaking at the wrong seasons. By deriving ENSO’s linear dynamics based on the conceptual recharge oscillator (RO) framework through the seasonal linear inverse model (sLIM) approach, various simulated phase-locking behaviors of CMIP models are systematically investigated in comparison with observations. In observations, phase-locking is mainly attributed to the seasonal modulation of ENSO’s SST growth rate. In contrast, in a significant portion of CMIP models, phase-locking is co-determined by the seasonal modulations of both SST growth and phase-transition rates. Further study of the joint effects of SST growth and phase-transition rates suggests that for simulating realistic winter peak ENSO phase-locking with the right dynamics, climate models need to have four key factors in the right combination: (1) correct phase of SST growth rate modulation peaking at the fall; (2) large enough amplitude for the annual cycle in growth rate; (3) amplitude of semi-annual cycle in growth rate needs to be small; and (4) amplitude of seasonal modulation in SST phase-transition rate needs to be small.more » « less
-
Abstract Stochastic oscillations can be characterized by a corresponding point process; this is a common practice in computational neuroscience, where oscillations of the membrane voltage under the influence of noise are often analyzed in terms of the interspike interval statistics, specifically the distribution and correlation of intervals between subsequent threshold-crossing times. More generally, crossing times and the corresponding interval sequences can be introduced for different kinds of stochastic oscillators that have been used to model variability of rhythmic activity in biological systems. In this paper we show that if we use the so-called mean-return-time (MRT) phase isochrons (introduced by Schwabedal and Pikovsky) to count the cycles of a stochastic oscillator with Markovian dynamics, the interphase interval sequence does not show any linear correlations, i.e., the corresponding sequence of passage times forms approximately a renewal point process. We first outline the general mathematical argument for this finding and illustrate it numerically for three models of increasing complexity: (i) the isotropic Guckenheimer–Schwabedal–Pikovsky oscillator that displays positive interspike interval (ISI) correlations if rotations are counted by passing the spoke of a wheel; (ii) the adaptive leaky integrate-and-fire model with white Gaussian noise that shows negative interspike interval correlations when spikes are counted in the usual way by the passage of a voltage threshold; (iii) a Hodgkin–Huxley model with channel noise (in the diffusion approximation represented by Gaussian noise) that exhibits weak but statistically significant interspike interval correlations, again for spikes counted when passing a voltage threshold. For all these models, linear correlations between intervals vanish when we count rotations by the passage of an MRT isochron. We finally discuss that the removal of interval correlations does not change the long-term variability and its effect on information transmission, especially in the neural context.more » « less
-
El Niño–Southern Oscillation (ENSO) events tend to peak at the end of the calendar year, a phenomenon called ENSO phase locking. This phase locking is a fundamental ENSO property that is determined by its basic dynamics. The conceptual ENSO recharge oscillator (RO) model is adopted to examine the ENSO phase-locking behavior in terms of its peak time, strength of phase locking, and asymmetry between El Niño and La Niña events. The RO model reproduces the main phase-locking characteristics found in observations, and the results show that the phase locking of ENSO is mainly dominated by the seasonal modulation of ENSO growth/decay rate. In addition, the linear/nonlinear mechanism of ENSO phase preference/phase locking is investigated using RO model. The difference between the nonlinear phase-locking mechanism and linear phase-preference mechanism is largely smoothed out in the presence of noise forcing. Further, the impact on ENSO phase locking from annual cycle modulation of the growth/decay rate, stochastic forcing, nonlinearity, and linear frequency are examined in the RO model. The preferred month of ENSO peak time depends critically on the phase and strength of the seasonal modulation of the ENSO growth/decay rate. Furthermore, the strength of phase locking is mainly controlled by the linear growth/decay rate, the amplitude of seasonal modulation of growth/decay rate, the amplitude of noise, the SST-dependent factor of multiplicative noise, and the linear frequency. The asymmetry of the sharpness of ENSO phase locking is induced by the asymmetric effect of state-dependent noise forcing in El Niño and La Niña events.more » « less
-
null (Ed.)There is a growing recognition that ecological systems can spend extended periods of time far away from an asymptotic state, and that ecological understanding will therefore require a deeper appreciation for how long ecological transients arise. Recent work has defined classes of deterministic mechanisms that can lead to long transients. Given the ubiquity of stochasticity in ecological systems, a similar systematic treatment of transients that includes the influence of stochasticity is important. Stochasticity can of course promote the appearance of transient dynamics by preventing systems from settling permanently near their asymptotic state, but stochasticity also interacts with deterministic features to create qualitatively new dynamics. As such, stochasticity may shorten, extend or fundamentally change a system’s transient dynamics. Here, we describe a general framework that is developing for understanding the range of possible outcomes when random processes impact the dynamics of ecological systems over realistic time scales. We emphasize that we can understand the ways in which stochasticity can either extend or reduce the lifetime of transients by studying the interactions between the stochastic and deterministic processes present, and we summarize both the current state of knowledge and avenues for future advances.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1007/s00422-022-00929-6
