Greater accuracy and broadened applicability of phase reduction using isostable coordinates.

*(English)*Zbl 1392.92007Summary: The applicability of phase models is generally limited by the constraint that the dynamics of a perturbed oscillator must stay near its underlying periodic orbit. Consequently, external perturbations must be sufficiently weak so that these assumptions remain valid. Using the notion of isostables of periodic orbits to provide a simplified coordinate system from which to understand the dynamics transverse to a periodic orbit, we devise a strategy to correct for changing phase dynamics for locations away from the limit cycle. Consequently, these corrected phase dynamics allow for perturbations of larger magnitude without invalidating the underlying assumptions of the reduction. The proposed reduction strategy yields a closed set of equations and can be applied to periodic orbits embedded in arbitrarily high dimensional spaces. We illustrate the utility of this strategy in two models with biological relevance. In the first application, we find that an optimal control strategy for modifying the period of oscillation can be improved with the corrected phase reduction. In the second, the corrected phase reduced dynamics are used to understand adaptation and memory effects resulting from past perturbations.

##### MSC:

92B05 | General biology and biomathematics |

34C15 | Nonlinear oscillations and coupled oscillators for ordinary differential equations |

##### Software:

XPPAUT
PDF
BibTeX
XML
Cite

\textit{D. Wilson} and \textit{B. Ermentrout}, J. Math. Biol. 76, No. 1--2, 37--66 (2018; Zbl 1392.92007)

Full Text:
DOI

##### References:

[1] | Ascher UM, Petzold LR (1998) Computer methods for ordinary differential equations and differential-algebraic equations, vol 61. SIAM, Philadelphia · Zbl 0908.65055 |

[2] | Azodi-Avval, R; Gharabaghi, A, Phase-dependent modulation as a novel approach for therapeutic brain stimulation, Front Comput Neurosci, 9, 26, (2015) |

[3] | Bagheri, N; Stelling, J; Doyle, FJ, Circadian phase resetting via single and multiple control targets, PLoS Comput Biol, 4, e1000104, (2008) |

[4] | Brown, E; Moehlis, J; Holmes, P, On the phase reduction and response dynamics of neural oscillator populations, Neural Comput, 16, 673-715, (2004) · Zbl 1054.92006 |

[5] | Budišić, M; Mohr, R; Mezić, I, Applied koopmanism. chaos: an interdisciplinary, J Nonlinear Sci, 22, 047510, (2012) · Zbl 1319.37013 |

[6] | Castejón, O; Guillamon, A; Huguet, G, Phase-amplitude response functions for transient-state stimuli, J Math Neurosci, 3, 13, (2013) · Zbl 1291.92026 |

[7] | Cherry, EM; Fenton, FH, Suppression of alternans and conduction blocks despite steep APD restitution: electrotonic, memory, and conduction velocity restitution effects, Am J Physiol Heart Circ Physiol, 55, h2332-h2341, (2004) |

[8] | Constantin P, Foias C, Nicolaenko B, Temam R (1989) Integral manifolds and inertial manifolds for dissipative partial differential equations, vol 70. Springer, New York · Zbl 0683.58002 |

[9] | Crowley, SJ; Acebo, C; Carskadon, MA, Sleep, Circadian rhythms, and delayed phase in adolescence, Sleep Med, 8, 602-612, (2007) |

[10] | Cui, J; Canavier, CC; Butera, RJ, Functional phase response curves: a method for understanding synchronization of adapting neurons, J Neurophysiol, 102, 387-398, (2009) |

[11] | Dean, DA; Forger, DB; Klerman, EB, Taking the lag out of jet lag through model-based schedule design, PLoS Comput Biol, 5, e1000418, (2009) |

[12] | Detrixhe, M; Doubeck, M; Moehlis, J; Gibou, F, A fast Eulerian approach for computation of global isochrons in high dimensions, SIAM J Appl Dyn Syst, 15, 1501-1527, (2016) · Zbl 1372.65331 |

[13] | Ermentrout, B; Pascal, M; Gutkin, B, The effects of spike frequency adaptation and negative feedback on the synchronization of neural oscillators, Neural Comput, 13, 1285-1310, (2001) · Zbl 0963.68647 |

[14] | Ermentrout GB (2002) Simulating, analyzing and animating dynamical systems: a guide to XPPAUT for researchers and students. SIAM, Philadelphia · Zbl 1003.68738 |

[15] | Ermentrout, GB; Kopell, N, Multiple pulse interactions and averaging in systems of coupled neural oscillators, J Math Biol, 29, 195-217, (1991) · Zbl 0718.92004 |

[16] | Ermentrout GB, Terman DH (2010) Mathematical foundations of neuroscience, vol 35. Springer, New York · Zbl 1320.92002 |

[17] | Foias, C; Sell, GR; Temam, R, Inertial manifolds for nonlinear evolutionary equations, J Differ Equ, 73, 309-353, (1988) · Zbl 0643.58004 |

[18] | Gonze, D; Bernard, S; Waltermann, C; Kramer, A; Herzel, H, Spontaneous synchronization of coupled Circadian oscillators, Biophys J, 89, 120-129, (2005) |

[19] | Govaerts, W; Sautois, B, Computation of the phase response curve: a direct numerical approach, Neural Comput, 18, 817-847, (2006) · Zbl 1087.92001 |

[20] | Grimshaw R (1993) Nonlinear ordinary differential equations, vol 2. CRC Press, Baca Raton · Zbl 0743.34002 |

[21] | Guckenheimer, J, Isochrons and phaseless sets, J Math Biol, 1, 259-273, (1975) · Zbl 0345.92001 |

[22] | Guckenheimer J, Holmes P (1983) Nonlinear oscillations, dynamical systems, and bifurcations of vector fields, vol 42. Springer, New York · Zbl 0515.34001 |

[23] | Guillamon, A; Huguet, G, A computational and geometric approach to phase resetting curves and surfaces, SIAM J Appl Dyn Syst, 8, 1005-1042, (2009) · Zbl 1216.34030 |

[24] | Hoppensteadt FC, Izhikevich EM (1997) Weakly connected neural networks. Springer, New York · Zbl 0887.92003 |

[25] | Huguet, G; Llave, R, Computation of limit cycles and their isochrons: fast algorithms and their convergence, SIAM J Appl Dyn Syst, 12, 1763-1802, (2013) · Zbl 1291.37108 |

[26] | Ichinose, N; Aihara, K; Judd, K, Extending the concept of isochrons from oscillatory to excitable systems for modeling an excitable neuron, Int J Bifurc Chaos, 8, 2375-2385, (1998) · Zbl 0935.92012 |

[27] | Izhikevich EM (2007) Dynamical systems in neuroscience: the geometry of excitability and bursting. MIT Press, London |

[28] | Jolly, MS; Kevrekidis, IG; Titi, ES, Approximate inertial manifolds for the Kuramoto-Sivashinsky equation: analysis and computations, Phys D Nonlinear Phenom, 44, 38-60, (1990) · Zbl 0704.58030 |

[29] | Kirk D (1998) Optimal control theory. Dover Publications, New York |

[30] | Klausmeier, CA, Floquet theory: a useful tool for understanding nonequilibrium dynamics, Theor Ecol, 1, 153-161, (2008) |

[31] | Kuramoto Y (1984) Chemical oscillations, waves, and turbulence. Springer, Berlin · Zbl 0558.76051 |

[32] | Kuramoto, Y, Phase-and center-manifold reductions for large populations of coupled oscillators with application to non-locally coupled systems, Int J Bifurc Chaos, 7, 789-805, (1997) · Zbl 0901.92015 |

[33] | Kurebayashi, W; Shirasaka, S; Nakao, H, Phase reduction method for strongly perturbed limit cycle oscillators, Phys Rev Lett, 111, 214101, (2013) |

[34] | Maran, SK; Canavier, CC, Using phase resetting to predict 1:1 and 2:2 locking in two neuron networks in which firing order is not always preserved, J Comput Neurosci, 24, 37-55, (2008) |

[35] | Mauroy, A; Mezić, I; Moehlis, J, Isostables, isochrons, and koopman spectrum for the action-angle representation of stable fixed point dynamics, Phys D Nonlinear Phenom, 261, 19-30, (2013) · Zbl 1284.37047 |

[36] | Nabi, A; Mirzadeh, M; Gibou, F; Moehlis, J, Minimum energy desynchronizing control for coupled neurons, J Comput Neurosci, 34, 259-271, (2013) · Zbl 1276.92017 |

[37] | Netoff T, Schwemmer MA, Lewis TJ (2012) Experimentally estimating phase response curves of neurons: theoretical and practical issues. In: Phase response curves in neuroscience. Springer, pp 95-129 |

[38] | Oprisan, SA; Prinz, AA; Canavier, CC, Phase resetting and phase locking in hybrid circuits of one model and one biological neuron, Biophys J, 87, 2283-2298, (2004) |

[39] | Osinga, HM; Moehlis, J, Continuation-based computation of global isochrons, SIAM J Appl Dyn Syst, 9, 1201-1228, (2010) · Zbl 1232.37014 |

[40] | Park, Y; Ermentrout, B, Weakly coupled oscillators in a slowly varying world, J Comput Neurosci, 40, 269-281, (2016) · Zbl 1382.93017 |

[41] | Pyragas, K; Novičenko, V, Phase reduction of a limit cycle oscillator perturbed by a strong amplitude-modulated high-frequency force, Phys Rev E, 92, 012910, (2015) |

[42] | Rabinovitch, A; Thieberger, R; Friedman, M, Forced bonhoeffer-van der Pol oscillator in its excited mode, Phys Rev E, 50, 1572, (1994) |

[43] | Roberts, AJ, Appropriate initial conditions for asymptotic descriptions of the long term evolution of dynamical systems. the journal of the Australian mathematical society. series B, Appl Math, 31, 48-75, (1989) · Zbl 0679.34044 |

[44] | Rosen, MR; Cohen, IS, Cardiac memory... new insights into molecular mechanisms, J Physiol, 570, 209-218, (2006) |

[45] | Rubin, J; Terman, D, High frequency stimulation of the subthalamic nucleus eliminates pathological thalamic rhythmicity in a computational model, J Comput Neurosci, 16, 211-235, (2004) |

[46] | Serkh, K; Forger, DB, Optimal schedules of light exposure for rapidly correcting Circadian misalignment, PLoS Comput Biol, 10, e1003523, (2014) |

[47] | Shirasaka, S; Kurebayashi, W; Nakao, H, Phase-amplitude reduction of transient dynamics far from attractors for limit-cycling systems. chaos: an interdisciplinary, J Nonlinear Sci, 27, 023119, (2017) · Zbl 1387.37045 |

[48] | Strogatz SH (1994) Nonlinear dynamics and chaos: with applications to physics, biology, chemistry, and engineering. Addison-Wesley, Boston · Zbl 1343.37001 |

[49] | Tass, PA; Qin, L; Hauptmann, C; Dovero, S; Bezard, E; Boraud, T; Meissner, WG, Coordinated reset has sustained aftereffects in Parkinsonian monkeys, Ann Neurol, 72, 816-820, (2012) |

[50] | Tolkacheva, EG; Schaeffer, DG; Gauthier, DJ; Krassowska, W, Condition for alternans and stability of the 1:1 response pattern in a model of paced cardiac dynamics, Phys Rev E, 67, 031904, (2003) |

[51] | Vreeswijk, C; Hansel, D, Patterns of synchrony in neural networks with spike adaptation, Neural Comput, 13, 959-992, (2001) · Zbl 1004.92011 |

[52] | Wedgwood, KCA; Lin, KK; Thul, R; Coombes, S, Phase-amplitude descriptions of neural oscillator models, J Math Neurosci, 3, 2, (2013) · Zbl 1291.92052 |

[53] | Wilson, D; Moehlis, J, An energy-optimal approach for entrainment of uncertain Circadian oscillators, Biophys J, 107, 1744-1755, (2014) |

[54] | Wilson, D; Moehlis, J, Clustered desynchronization from high-frequency deep brain stimulation, PLoS Comput Biol, 11, e1004673, (2015) |

[55] | Wilson, D; Moehlis, J, Isostable reduction of periodic orbits, Phys Rev E, 94, 052213, (2016) |

[56] | Winfree A (2001) The geometry of biological time, 2nd edn. Springer, New York · Zbl 1014.92001 |

[57] | Xu, Y; Padiath, QS; Shapiro, RE; Jones, CR; Wu, SC; Saigoh, N; Saigoh, K; Ptáček, LJ; Fu, YH, Functional consequences of a CKI\(δ \) mutation causing familial advanced sleep phase syndrome, Nature, 434, 640-644, (2005) |

This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.