In a series of two papers, we investigate the mechanisms by which complex oscillations are generated in a class of nonlinear dynamical systems with resets modeling the voltage and adaptation of neurons. This first paper presents mathematical analysis showing that the system can support bursts of any period as a function of model parameters. In continuous dynamical systems with resets, period-incrementing structures are complex to analyze. In the present context, we use the fact that bursting patterns correspond to periodic orbits of the adaptation map that governs the sequence of values of the adaptation variable at the resets. Using a slow-fast approach, we show that this map converges towards a piecewise linear discontinuous map whose orbits are exactly characterized. That map shows a period-incrementing structure with instantaneous transitions. We show that the period-incrementing structure persists for the full system with non-constant adaptation, but the transitions are more complex. We investigate the presence of chaos at the transitions.

A. Map-based neuron models are an important tool in modelling neural dynamics and sometimes can be considered as an alternative to usually computationally costlier models based on continuous or hybrid dynamical systems. However, due to their discrete nature, rigorous mathematical analysis might be challenging. We study a discrete model of neuronal dynamics introduced by Chialvo [Chaos, Solitons & Fractals 5, 1995, 461-479]. In particular, we show that its reduced one-dimensional version can be treated as an independent simple model of neural activity where the input and the fixed value of the recovery variable are parameters. This one-dimensional model still displays very rich and varied dynamics. Using the fact that the map whose iterates define voltage dynamics is S-unimodal, we describe in detail both the periodic behaviour and the occurrence of different notions of chaos, indicating corresponding regions in parameter space. Our study is also complemented by a bifurcation analysis of the mentioned dynamical model.

This work continues the analysis of complex dynamics in a class of bidimensional nonlinear hybrid dynamical systems with resets modeling neuronal voltage dynamics with adaptation and spike emission. We show that these models can generically display a form of mixed-mode oscillations (MMOs), which are trajectories featuring an alternation of small oscillations with spikes or bursts (multiple consecutive spikes). The mechanism by which these are generated relies fundamentally on the hybrid structure of the flow: invariant manifolds of the continuous dynamics govern small oscillations, while discrete resets govern the emission of spikes or bursts, contrasting with classical MMO mechanisms in ordinary differential equations involving more than three dimensions and generally relying on a timescale separation. The decomposition of mechanisms reveals the geometrical origin of MMOs, allowing a relatively simple classification of points on the reset manifold associated to specific numbers of small oscillations. We show that the MMO pattern can be described through the study of orbits of a discrete adaptation map, which is singular as it features discrete discontinuities with unbounded left-and right-derivatives. We study orbits of the map via rotation theory for discontinuous circle maps and elucidate in detail complex behaviors arising in the case where MMOs display at most one small oscillation between each consecutive pair of spikes.

We conduct computer-assisted analysis of a two-dimensional model of a neuron introduced by Chialvo in 1995 [Chaos, Solitons Fractals 5, 461–479]. We apply the method of rigorous analysis of global dynamics based on a set-oriented topological approach, introduced by Arai et al. in 2009 [SIAM J. Appl. Dyn. Syst. 8, 757–789] and improved and expanded afterward. Additionally, we introduce a new algorithm to analyze the return times inside a chain recurrent set. Based on this analysis, together with the information on the size of the chain recurrent set, we develop a new method that allows one to determine subsets of parameters for which chaotic dynamics may appear. This approach can be applied to a variety of dynamical systems, and we discuss some of its practical aspects.