Math. Model. Nat. Phenom.
Volume 17, 2022
Systems with Hysteresis and Switching
|Number of page(s)||26|
|Published online||09 June 2022|
A multiple timescale network model of intracellular calcium concentrations in coupled neurons: Insights from ROM simulations
Ecuaciones Diferenciales y Análisis Numérico, Universidad de Sevilla, Calle Tarfia s/n, Seville 41012, Spain
2 IMUS, Universidad de Sevilla, Calle Tarfia s/n, Seville 41012, Spain
3 Laboratoire de Mathématiques et Modélisation d’Évry (LAMME), Univ Evry, CNRS, Université Paris-Saclay, IBGBI, 23 Bld de France, Evry 91037, France
* Corresponding author: firstname.lastname@example.org
Accepted: 11 April 2022
In Fernández-García and Vidal [Physica D 401 (2020) 132129], the authors have analyzed the synchronization features between two identical 3D slow-fast oscillators, symmetrically coupled, and built as an extension of the FitzHugh—Nagumo dynamics generating Mixed-Mode Oscillations. The third variable in each oscillator aims at representing the time-varying intracellular calcium concentration in neurons. The global model is therefore six-dimensional with two fast variables and four slow variables with strong symmetry properties. In the present article, we consider an extension of this model in two different directions. First, we consider heterogeneity among cells and analyze the coupling of two oscillators with different values for one parameter which tunes the intrinsic frequency of the output. We therefore identify new patterns of antiphasic synchronization, with non trivial signatures and that exhibit a Devil’s Staircase phenomenon in signature transitions when varying the coupling gain parameter value. Second, we introduce a network of N cells divided into two clusters: the coupling between neurons in each cluster is excitatory, while the coupling between the two clusters is inhibitory. Such system aims at modelling the interactions between neurons tending to synchronization in each of two subpopulations that inhibit each other, like ipsi- and contra-lateral motoneurons assemblies. To perform the numerical simulations in this case when N is large, as an initial step towards the network analysis, we consider Reduced Order Models in order to save computational costs. We present the numerical reduction results in a network of 100 cells. For the sake of validation of the numerical reduction method, we both compare the outputs and CPU times obtained with the original and the reduced models in different cases of network coupling structures.
Mathematics Subject Classification: 34C15 / 34C25 / 34C27 / 34C28 / 92B20 / 34K28 / 37M05 / 70K70
Key words: Slow-fast dynamics / coupled oscillators / mixed-mode oscillations / synchronization / network neuron model / reduced order models
© The authors. Published by EDP Sciences, 2022
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.