A Universal Description of Stochastic Oscillators

Abstract

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 new complex-valued function Q * 1 (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 Q * 1 (x) is the eigenfunction of the Kolmogorov backward operator with the least negative (but non-vanishing) 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.