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.