Mixed-mode oscillations for slow-fast perturbed systems

This article concerns the dynamics of mixed-mode oscillations (MMOs) emerging from the calcium-based inner hair cells (IHCs) model in the auditory cortex. The paper captures the MMOs generation mechanism based on the geometric singular perturbation theory (GSPT) after exploiting the average analysis for reducing the full model. Our analysis also finds that the critical manifold and folded surface are central to the mechanism of the existence of MMOs at the folded saddle for the perturbed system. The system parameters, such like the maximal calcium channels conductance, controls the firing patterns, and many new oscillations occur for the IHCs model. Tentatively, we conduct dynamic analysis combined with dynamic method based on GSPT by giving slow-fast analysis for the singular perturbed models and bifurcation analysis. In particular, we explore the two-slow-two-fast and three-slow-one-fast IHCs perturbed systems with layer and reduced problems so that differential-algebraic equations are obtained. This paper reveals the underlying dynamic properties of perturbed systems under singular perturbation theory.

Tipping Cascades in a Multi-patch System with Noise and Spatial Coupling

Forecasting tipping points in spatially extended systems is a key area of interest to ecologists. A slowly declining spatially distributed population is an important example of an ecological system that could exhibit a cascade of tipping points. Here, we develop a spatial two-patch model with environmental stochasticity that is slowly forced through population collapse, in the presence of changing environmental conditions. We begin with a basic spatial model, then introduce a fast–slow version of the model using geometric singular perturbation theory, followed by the inclusion of stochasticity. Using the spectral density of the fluctuating subpopulation in each patch, we derive analytic expressions for candidate indicators of population extinction and evaluate their performance through a simulation study. We find that coupling and spatial heterogeneity decrease the magnitude of the proposed indicators in coupled populations relative to isolated populations. Moreover, the degree of coupling dictates the trends in summary statistics. We conclude that this theory may be applied to other contexts, including the control of invasive species.

Canard solutions in neural mass models: consequences on critical regimes

Mathematical models at multiple temporal and spatial scales can unveil the fundamental mechanisms of critical transitions in brain activities. Neural mass models (NMMs) consider the average temporal dynamics of interconnected neuronal subpopulations without explicitly representing the underlying cellular activity. The mesoscopic level offered by the neural mass formulation has been used to model electroencephalographic (EEG) recordings and to investigate various cerebral mechanisms, such as the generation of physiological and pathological brain activities. In this work, we consider a NMM widely accepted in the context of epilepsy, which includes four interacting neuronal subpopulations with different synaptic kinetics. Due to the resulting three-time-scale structure, the model yields complex oscillations of relaxation and bursting types. By applying the principles of geometric singular perturbation theory, we unveil the existence of the canard solutions and detail how they organize the complex oscillations and excitability properties of the model. In particular, we show that boundaries between pathological epileptic discharges and physiological background activity are determined by the canard solutions. Finally we report the existence of canard-mediated small-amplitude frequency-specific oscillations in simulated local field potentials for decreased inhibition conditions. Interestingly, such oscillations are actually observed in intracerebral EEG signals recorded in epileptic patients during pre-ictal periods, close to seizure onsets.

Controlling Canard Cycles

Canard cycles are periodic orbits that appear as special solutions of fast-slow systems (or singularly perturbed ordinary differential equations). It is well known that canard cycles are difficult to detect, hard to reproduce numerically, and that they are sensible to exponentially small changes in parameters. In this paper, we combine techniques from geometric singular perturbation theory, the blow-up method, and control theory, to design controllers that stabilize canard cycles of planar fast-slow systems with a folded critical manifold. As an application, we propose a controller that produces stable mixed-mode oscillations in the van der Pol oscillator.

Geometric analysis of mixed-mode oscillations in a model of electrical activity in human beta-cells

Human pancreatic beta-cells may exhibit complex mixed-mode oscillatory electrical activity, which underlies insulin secretion. A recent biophysical model of human beta-cell electrophysiology can simulate such bursting behavior, but a mathematical understanding of the model’s dynamics is still lacking. Here we exploit time-scale separation to simplify the original model to a simpler three-dimensional model that retains the behavior of the original model and allows us to apply geometric singular perturbation theory to investigate the origin of mixed-mode oscillations. Changing a parameter modeling the maximal conductance of a potassium current, we find that the reduced model possesses a singular Hopf bifurcation that results in small-amplitude oscillations, which go through a period-doubling sequence and chaos until the birth of a large-scale return mechanism and bursting dynamics. The theory of folded node singularities provide insight into the bursting dynamics further away from the singular Hopf bifurcation and the eventual transition to simple spiking activity. Numerical simulations confirm that the insight obtained from the analysis of the reduced model can be lifted back to the original model.

Geometric analysis of a two-body problem with quick loss of mass

We consider a two-body problem with quick loss of mass which was formulated by Verhulst (Verhulst in J Inst Math Appl 18: 87–98, 1976). The corresponding dynamical system is singularly perturbed due to the presence of a small parameter in the governing equations which corresponds to the reciprocal of the initial rate of loss of mass, resulting in a boundary layer in the asymptotics. Here, we showcase a geometric approach which allows us to derive asymptotic expansions for the solutions of that problem via a combination of geometric singular perturbation theory (Fenichel in J Differ Equ 31: 53–98, 1979) and the desingularization technique known as “blow-up” (Dumortier, in: Bifurcations and Periodic Orbits of Vector Fields, Springer, Dordrecht, 1993). In particular, we justify the unexpected dependence of those expansions on fractional powers of the singular perturbation parameter; moreover, we show that the occurrence of logarithmic (“switchback”) terms therein is due to a resonance phenomenon that arises in one of the coordinate charts after blow-up.

Geometric Analysis of Mixed-mode Oscillations in a Model of Electrical Activity in Human Beta-cells

Human pancreatic beta-cells may exhibit complex mixed-mode oscillatory electrical activity, which underlies insulin secretion. A recent biophysical model of human beta-cell electrophysiology can simulate such bursting behavior, but a mathematical understanding of the model's dynamics is still lacking. Here we exploit time-scale separation to simplify the original model to a simpler three-dimensional model that retains the behavior of the original model and allows us to apply geometric singular perturbation theory to investigate the origin of mixed-mode oscillations. Changing a parameter modeling the maximal conductance of a potassium current, we nd that the reduced model possesses a singular Hopf bifurcation that results in small-amplitude oscillations, which go through a period-doubling sequence and chaos until the birth of a large-scale return mechanism and bursting dynamics. The theory of folded node singularities provide insight into the bursting dynamics further away from the singular Hopf bifurcation and the eventual transition to simple spiking activity. Numerical simulations confirm that the insight obtained from the analysis of the reduced model can be lifted back to the original model.