The phase response curve (PRC) is a powerful tool to study the effect of a perturbation on the phase of an oscillator, assuming that all the dynamics can be explained by the phase variable. However, factors like the rate of convergence to the oscillator, strong forcing or high stimulation frequency may invalidate the above assumption and raise the question of how is the phase variation away from an attractor. The concept of isochrons turns out to be crucial to answer this question; from it, we have built up Phase Response Functions (PRF) and, in the present paper, we complete the extension of advancement functions to the transient states by defining the Amplitude Response Function (ARF) to control changes in the transversal variables. Based on the knowledge of both the PRF and the ARF, we study the case of a pulse-train stimulus, and compare the predictions given by the PRC-approach (a 1D map) to those given by the PRF-ARF-approach (a 2D map); we observe differences up to two orders of magnitude in favor of the 2D predictions, especially when the stimulation frequency is high or the strength of the stimulus is large. We also explore the role of hyperbolicity of the limit cycle as well as geometric aspects of the isochrons. Summing up, we aim at enlightening the contribution of transient effects in predicting the phase response and showing the limits of the phase reduction approach to prevent from falling into wrong predictions in synchronization problems.List of AbbreviationsPRC phase response curve, phase resetting curve.PRF phase response function.ARF amplitude response function.

In this paper we prove the breakdown of a heteroclinic connection in the analytic\ud versal unfoldings of the generic Hopf-zero singularity in an open set of the parameter space.\ud This heteroclinic orbit appears at any order if one performs the normal form around the\ud origin, therefore it is a phenomenon “beyond all orders”. In this paper we provide a formula\ud for the distance between the corresponding stable and unstable one-dimensional manifolds\ud which is given by an exponentially small function in the perturbation parameter. Our result\ud applies both for conservative and dissipative unfoldingsPeer ReviewedPostprint (published version

In this paper we study a beyond all order phenomenon which appears in the analytic unfoldings of the Hopf-zero singularity. It consists in the breakdown of a two-dimensional heteroclinic surface which exists in the truncated normal form of this singularity at any order. The results in this paper are twofold: on the one hand we give results for generic unfoldings which lead to sharp exponentially small upper bounds of the difference between these manifolds. On the other hand, we provide asymptotic formulas for this difference by means of the Melnikov function for some non-generic unfoldings.

In this paper we prove the breakdown of the two-dimensional stable and unstable manifolds associated to two saddle-focus points which appear in the unfoldings of the Hopf-zero singulariry. The method consists in obtaining an asymptotic formula for the difference between this manifolds which turns to be exponentially small respect to the unfolding parameter. The formula obtained is explicit but depends on the so-called Stokes constants, which arise in the study of original vector field and which corresponds to the so called inner equation in singular perturbation theory.

Phase response curves (PRCs) have been extensively used to control the phase of oscillators under perturbations. Their main advantage is the reduction of the whole model dynamics to a single variable (phase) dynamics. However, in some adverse situations (strong inputs, high-frequency stimuli, weak convergence,. . . ), the phase reduction does not provide enough information and, therefore, PRC lose predictive power. To overcome this shortcoming, in the last decade, new contributions have appeared that allow to reduce the system dynamics to the phase plus some transversal variable that controls the deviations from the asymptotic behaviour. We call this setting extended response functions. In particular, we single out the phase response function (PRF, a generalization of the PRC) and the amplitude response function (ARF) that account for the above-mentioned deviations from the oscillating attractor. It has been shown that in adverse situations, the PRC misestimate the phase dynamics whereas the PRF-ARF system provides accurate enough predictions. In this paper, we address the problem of studying the dynamics of the PRF-ARF systems under periodic pulsatile stimuli. This paradigm leads to a two-dimensional discrete dynamical system that we call 2D entrainment map. By using advanced methods to study invariant manifolds and the dynamics inside them, we construct an analytico-numerical method to track the invariant curves induced by the stimulus as two crucial parameters of the system increase (the strength of the input and its frequency). Our methodology also incorporates the computation of Arnold tongues associated to the 2D entrainment map. We apply the method developed to study inner dynamics of the invariant curves of a canonical type II oscillator model. We further compare the Arnold tongues of the 2D map with those obtained with the map induced only by the PRC, which give already noticeable differences. We also observe (via simulations) how high-frequency or strong enough stimuli break up the oscillatory dynamics and lead to phase-locking, which is well captured by the 2D entrainment map.