Computation of phase response curves via a direct method adapted to infinitesimal perturbations

A new numerical algorithm for computation of phase response curves of stable limit cycle oscillators is proposed. The idea of the algorithm originates from a direct method that is based on computation of the oscillator response to short finite pulses delivered at different phases of oscillations. Here we adapt the direct method to the case of infinitesimal perturbations and compare our algorithm with the standard algorithm based on the backward integration of the adjoint equations. In contrast to the standard algorithm, our algorithm does not require any backward integration and it is easier to program since a necessity of numerical interpolation for the Jacobian matrix is avoided. In addition, we demonstrate by examples that our algorithm is faster than the standard algorithm and this advantage is especially notable for weakly stable limit cycle oscillators.

[1]  Yoshiki Kuramoto,et al.  Chemical Oscillations, Waves, and Turbulence , 1984, Springer Series in Synergetics.

[2]  John Guckenheimer,et al.  Dissecting the Phase Response of a Model Bursting Neuron , 2009, SIAM J. Appl. Dyn. Syst..

[3]  A. Winfree.  The geometry of biological time , 1991 .

[4]  W. Govaerts,et al.  Computation of the Phase Response Curve: A Direct Numerical Approach , 2006, Neural Computation.

[5]  Bard Ermentrout,et al.  Simulating, analyzing, and animating dynamical systems - a guide to XPPAUT for researchers and students , 2002, Software, environments, tools.

[6]  C. Johnson,et al.  Forty years of PRCs--what have we learned? , 1999, Chronobiology international.

[7]  E. Izhikevich,et al.  Weakly connected neural networks , 1997 .

[8]  A. Hodgkin,et al.  A quantitative description of membrane current and its application to conduction and excitation in nerve , 1990 .

[9]  G Bard Ermentrout,et al.  Efficient estimation of phase-resetting curves in real neurons and its significance for neural-network modeling. , 2005, Physical review letters.

[10]  A. Reyes,et al.  Two modes of interspike interval shortening by brief transient depolarizations in cat neocortical neurons. , 1993, Journal of neurophysiology.

[11]  A. Winfree.  Patterns of phase compromise in biological cycles , 1974 .

[12]  G. Ermentrout,et al.  Multiple pulse interactions and averaging in systems of coupled neural oscillators , 1991 .

[13]  I. G. Malkin,et al.  Some problems in the theory of nonlinear oscillations , 1959 .

[14]  Eric Shea-Brown,et al.  On the Phase Reduction and Response Dynamics of Neural Oscillator Populations , 2004, Neural Computation.

[15]  P. McClintock.  Phase resetting in medicine and biology , 2008 .

[16]  Jürgen Kurths,et al.  Synchronization - A Universal Concept in Nonlinear Sciences , 2001, Cambridge Nonlinear Science Series.

[17]  H. Robinson,et al.  Phase resetting curves and oscillatory stability in interneurons of rat somatosensory cortex. , 2007, Biophysical journal.

[18]  N. Ikeda,et al.  Model of bidirectional interaction between myocardial pacemakers based on the phase response curve , 1982, Biological Cybernetics.

[19]  Jürgen Kurths,et al.  Synchronization: Phase locking and frequency entrainment , 2001 .

[20]  Gemma Huguet,et al.  A Computational and Geometric Approach to Phase Resetting Curves and Surfaces , 2009, SIAM J. Appl. Dyn. Syst..

[21]  Jeff Moehlis,et al.  Linear control of neuronal spike timing using phase response curves , 2009, 2009 Annual International Conference of the IEEE Engineering in Medicine and Biology Society.

[22]  Dimitrios I. Fotiadis,et al.  Phase response characteristics of sinoatrial node cells , 2007, Comput. Biol. Medicine.

[23]  C. Morris,et al.  Voltage oscillations in the barnacle giant muscle fiber. , 1981, Biophysical journal.

[24]  J. W. Hastings,et al.  A PERSISTENT DIURNAL RHYTHM OF LUMINESCENCE IN GONYAULAX POLYEDRA , 1958 .

[25]  Bard Ermentrout,et al.  Type I Membranes, Phase Resetting Curves, and Synchrony , 1996, Neural Computation.

[26]  J. Hindmarsh,et al.  A model of neuronal bursting using three coupled first order differential equations , 1984, Proceedings of the Royal Society of London. Series B. Biological Sciences.

[27]  Alexander P. Kuznetsov,et al.  Coupled van der Pol–Duffing oscillators: Phase dynamics and structure of synchronization tongues , 2009 .

[28]  P. McClintock.  Synchronization:a universal concept in nonlinear science , 2003 .

[29]  Jeff Moehlis,et al.  Linear Control of Neuronal Spike Timing Using Phase Response Curves , 2010 .

[30]  A. Winfree.  When time breaks down , 1987 .

[31]  Alvin Shrier,et al.  Chaos in neurobiology , 1983, IEEE Transactions on Systems, Man, and Cybernetics.

[32]  J. Guckenheimer,et al.  Isochrons and phaseless sets , 1975, Journal of mathematical biology.

[33]  Eugene M. Izhikevich,et al.  Dynamical Systems in Neuroscience: The Geometry of Excitability and Bursting , 2006 .

[34]  G. Ermentrout.  Stable periodic solutions to discrete and continuum arrays of weakly coupled nonlinear oscillators , 1992 .