Source author record

Mathieu Desroches

Mathieu Desroches appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

15works
11topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

15 published item(s)

preprint2022arXiv

Cross-scale excitability in networks of synaptically-coupled quadratic integrate-and-fire neurons

From the action potentials of neurons and cardiac cells to the amplification of calcium signals in oocytes, excitability is a hallmark of many biological signalling processes. In recent years, excitability in single cells has been related to multiple-timescale dynamics through canards, special solutions which determine the effective thresholds of the all-or-none responses. However, the emergence of excitability in large populations remains an open problem. Here, we show that the mechanism of excitability in large networks and mean-field descriptions of coupled quadratic integrate-and-fire (QIF) cells mirrors that of the individual components. We initially exploit the Ott-Antonsen ansatz to derive low-dimensional dynamics for the coupled network and use it to describe the structure of canards via slow periodic forcing. We demonstrate that the thresholds for onset and offset of population firing can be found in the same way as those of the single cell. We combine theoretical analysis and numerical computations to develop a novel and comprehensive framework for excitability in large populations, applicable not only to models amenable to Ott-Antonsen reduction, but also to networks without a closed-form mean-field limit, in particular sparse networks.

preprint2021arXiv

Spike-adding and reset-induced canard cycles in adaptive integrate and fire models

We study a class of planar integrate and fire (IF) models called adaptive integrate and fire (AIF) models, which possesses an adaptation variable on top of membrane potential, and whose subthreshold dynamics is piecewise linear (PWL). These AIF models therefore have two reset conditions, which enable bursting dynamics to emerge for suitable parameter values. Such models can be thought of as hybrid dynamical systems. We consider a particular slow dynamics within AIF models and prove the existence of bursting cycles with $N$ resets, for any integer $N$. Furthermore, we study the transition between $N$- and $(N+1)$-reset cycles upon vanishingly small parameter variations and prove (for $N=2$) that such transitions are organised by canard cycles. Finally, using numerical continuation we compute branches of bursting cycles, including canard-explosive branches, in these AIF models, by suitably recasting the periodic problem as a two-point boundary-value problem.

preprint2020arXiv

Towards a new classification of bursting patterns: review & extensions

The mathematical classification of complex bursting oscillations in multiscale excitable systems, seen for example in physics and neuroscience, has been the subject of active enquiry since the early 1980s. This classification problem is fundamental as it also establishes analytical and numerical foundations for studying complex temporal behaviours in multiple timescale models. This manuscript begins by reviewing the seminal works of Rinzel and Izhikevich in classifying bursting patterns of excitable cell models. Moreover, we recall an alternative, yet complementary, mathematical classification approach by Golubitsky, which together with the Rinzel-Izhikevich proposals provide the state-of-the-art foundations to the classification problem. Unexpectedly, while keeping within the Rinzel-Izhikevich framework, we find novel cases of bursting mechanisms not considered before. Moving beyond the state-of-the-art, we identify novel bursting mechanisms that fall outside the current classifications. This leads us towards a new classification, which requires the analysis of both the fast and the slow subsystems of an underlying slow-fast model. This new classification allows the dynamical dissection of a larger class of bursters. To substantiate this, we add a new class of bursters with at least two slow variables, which we denote folded-node bursters, to convey the idea that the bursts are initiated or annihilated via a folded-node singularity. In fact, there are two main families of folded-node bursters, depending upon the phase of the bursting cycle during which folded-node dynamics occurs. If it occurs during the silent phase, we obtain the classical folded-node bursting (or simply folded-node bursting). If it occurs during the active phase, we have cyclic folded-node bursting. We classify both families and give examples of minimal systems displaying these novel types of bursting behaviour.

preprint2020arXiv

Why we should use Topological Data Analysis in Ageing: towards defining the "Topological shape of ageing"

Living systems are subject to the arrow of time; from birth, they undergo complex transformations (self-organization) in a constant battle for survival, but inevitably ageing and disease trap them to death. Can ageing be understood and eventually reversed? What tools can be employed to further our understanding of ageing? The present article is an invitation for biologists and clinicians to consider key conceptual ideas and computational tools (known to mathematicians and physicists), which potentially may help dissect some of the underlying processes of ageing and disease. Specifically, we first discuss how to classify and analyze complex systems, as well as highlight critical theoretical difficulties that make complex systems hard to study. Subsequently, we introduce Topological Data Analysis - a novel Big Data tool - which may help in the study of complex systems since it extracts knowledge from data in a holistic approach via topological considerations. These conceptual ideas and tools are discussed in a rather informal way to pave future discussions and collaborations between mathematicians and biologists studying ageing.

preprint2016arXiv

A modular architecture for transparent computation in Recurrent Neural Networks

Computation is classically studied in terms of automata, formal languages and algorithms; yet, the relation between neural dynamics and symbolic representations and operations is still unclear in traditional eliminative connectionism. Therefore, we suggest a unique perspective on this central issue, to which we would like to refer as to transparent connectionism, by proposing accounts of how symbolic computation can be implemented in neural substrates. In this study we first introduce a new model of dynamics on a symbolic space, the versatile shift, showing that it supports the real-time simulation of a range of automata. We then show that the Goedelization of versatile shifts defines nonlinear dynamical automata, dynamical systems evolving on a vectorial space. Finally, we present a mapping between nonlinear dynamical automata and recurrent artificial neural networks. The mapping defines an architecture characterized by its granular modularity, where data, symbolic operations and their control are not only distinguishable in activation space, but also spatially localizable in the network itself, while maintaining a distributed encoding of symbolic representations. The resulting networks simulate automata in real-time and are programmed directly, in absence of network training. To discuss the unique characteristics of the architecture and their consequences, we present two examples: i) the design of a Central Pattern Generator from a finite-state locomotive controller, and ii) the creation of a network simulating a system of interactive automata that supports the parsing of garden-path sentences as investigated in psycholinguistics experiments.

preprint2015arXiv

Canards and curvature: nonsmooth approximation by pinching

In multiple time-scale (singularly perturbed) dynamical systems, canards are counterintuitive solutions that evolve along both attracting and repelling invariant manifolds. In two dimensions, canards result in periodic oscillations whose amplitude and period grow in a highly nonlinear way: they are slowly varying with respect to a control parameter, except for an exponentially small range of values where they grow extremely rapidly. This sudden growth, called a canard explosion, has been encountered in many applications ranging from chemistry to neuronal dynamics, aerospace engineering and ecology. Canards were initially studied using nonstandard analysis, and later the same results were proved by standard techniques such as matched asymptotics, invariant manifold theory and parameter blow-up. More recently, canard-like behaviour has been linked to surfaces of discontinuity in piecewise-smooth dynamical systems. This paper provides a new perspective on the canard phenomenon by showing that the nonstandard analysis of canard explosions can be recast into the framework of piecewise-smooth dynamical systems. An exponential coordinate scaling is applied to a singularly perturbed system of ordinary differential equations. The scaling acts as a lens that resolves dynamics across all time-scales. The changes of local curvature that are responsible for canard explosions are then analysed. Regions where different time-scales dominate are separated by hypersurfaces, and these are pinched together to obtain a piecewise-smooth system, in which curvature changes manifest as discontinuity-induced bifurcations. The method is used to classify canards in arbitrary dimensions, and to derive the parameter values over which canards form either small cycles (canards without head) or large cycles (canards with head).

preprint2015arXiv

Nonsmooth analogues of slow-fast dynamics -- pinching at a folded node

The folded node is a singularity associated with loss of normal hyperbolicity in systems where mixtures of slow and fast timescales arise due to singular perturbations. Canards are special solutions that reveal a counteractive feature of the local dynamics, namely trajectories that flow from attractive regions of space into repulsive regions. An alternative way to model switches between timescales is using piecewise-smooth differential equations. There is presently no adequate theory or method for relating slow-fast and piecewise-smooth models. Here we derive the analogous piecewise-smooth system for the folded node by pinching phase space to sharpen the switch between timescales. The corresponding piecewise-smooth system contains a so-called two-fold singularity, and exhibits the same topology and number of canards as the slow-fast system. Thus pinching provides a piecewise-smooth approximation to a slow-fast system.

preprint2015arXiv

Turing Computation with Recurrent Artificial Neural Networks

We improve the results by Siegelmann & Sontag (1995) by providing a novel and parsimonious constructive mapping between Turing Machines and Recurrent Artificial Neural Networks, based on recent developments of Nonlinear Dynamical Automata. The architecture of the resulting R-ANNs is simple and elegant, stemming from its transparent relation with the underlying NDAs. These characteristics yield promise for developments in machine learning methods and symbolic computation with continuous time dynamical systems. A framework is provided to directly program the R-ANNs from Turing Machine descriptions, in absence of network training. At the same time, the network can potentially be trained to perform algorithmic tasks, with exciting possibilities in the integration of approaches akin to Google DeepMind's Neural Turing Machines.

preprint2013arXiv

A method for the reconstruction of unknown non-monotonic growth functions in the chemostat

We propose an adaptive control law that allows one to identify unstable steady states of the open-loop system in the single-species chemostat model without the knowledge of the growth function. We then show how one can use this control law to trace out (reconstruct) the whole graph of the growth function. The process of tracing out the graph can be performed either continuously or step-wise. We present and compare both approaches. Even in the case of two species in competition, which is not directly accessible with our approach due to lack of controllability, feedback control improves identifiability of the non-dominant growth rate.

preprint2013arXiv

Extremum seeking via continuation techniques for optimizing biogas production in the chemostat

We consider the chemostat model with the substrate concentration as the single measurement. We propose a control strategy that drives the system at a steady state maximizing the gas production without the knowledge of the specific growth rate. Our approach separates the extremum seeking problem from the feedback control problem such that each of the two subproblems can be solved with relatively simple algorithms. We are then free to choose any numerical optimization algorithm. We give a demonstration for two choices: one is based on slow-fast dynamics and numerical continuation, the other is a combination of golden-section and Newton iteration. The method copes with non-monotonic growth functions.

preprint2013arXiv

Noise-induced canard and mixed-mode oscillations in large stochastic networks with multiple timescales

We investigate the dynamics of large stochastic networks with different timescales and nonlinear mean-field interactions. After deriving the limit equations for a general class of network models, we apply our results to the celebrated Wilson-Cowan system with two populations with or without slow adaptation, paradigmatic example of nonlinear mean-field network. This system has the property that the dynamics of the mean of the solution exactly satisfies an ODE. This reduction allows to show that in the mean-field limit and in multiple populations with multiple timescales, noise induces canard explosions and Mixed-Mode Oscillations on the mean of the solution. This sheds new light on the qualitative effects of noise and sensitivity to precise noise values in large stochastic networks. We further investigate finite-sized networks and show that systematic differences with the mean-field limits arise in bistable regimes (where random switches between different attractors occur) or in mixed-mode oscillations, were the finite-size effects induce early jumps due to the sensitivity of the attractor.

preprint2013arXiv

Short-term synaptic plasticity in the deterministic Tsodyks-Markram model leads to unpredictable network dynamics

Short-Term Synaptic Plasticity (STSP) strongly affects the neural dynamics of cortical networks. The Tsodyks and Markram (TM) model for STSP accurately accounts for a wide range of physiological responses at different types of cortical synapses. Here, we report for the first time a route to chaotic behavior via a Shilnikov homoclinic bifurcation that dynamically organises some of the responses in the TM model. In particular, the presence of such a homoclinic bifurcation strongly affects the shape of the trajectories in the phase space and induces highly irregular transient dynamics; indeed, in the vicinity of the Shilnikov homoclinic bifurcation, the number of population spikes and their precise timing are unpredictable and highly sensitive to the initial conditions. Such an irregular deterministic dynamics has its counterpart in stochastic/network versions of the TM model: the existence of the Shilnikov homoclinic bifurcation generates complex and irregular spiking patterns and --acting as a sort of springboard-- facilitates transitions between the down-state and unstable periodic orbits. The interplay between the (deterministic) homoclinic bifurcation and stochastic effects may give rise to some of the complex dynamics observed in neural systems.

preprint2012arXiv

Mixed-mode oscillations in a multiple time scale phantom bursting system

In this work we study mixed mode oscillations in a model of secretion of GnRH (Gonadotropin Releasing Hormone). The model is a phantom burster consisting of two feedforward coupled FitzHugh-Nagumo systems, with three time scales. The forcing system (Regulator) evolves on the slowest scale and acts by moving the slow nullcline of the forced system (Secretor). There are three modes of dynamics: pulsatility (transient relaxation oscillation), surge (quasi steady state) and small oscillations related to the passage of the slow nullcline through a fold point of the fast nullcline. We derive a variety of reductions, taking advantage of the mentioned features of the system. We obtain two results; one on the local dynamics near the fold in the parameter regime corresponding to the presence of small oscillations and the other on the global dynamics, more specifically on the existence of an attracting limit cycle. Our local result is a rigorous characterization of small canards and sectors of rotation in the case of folded node with an additional time scale, a feature allowing for a clear geometric argument. The global result gives the existence of an attracting unique limit cycle, which, in some parameter regimes, remains attracting and unique even during passages through a canard explosion.

preprint2011arXiv

A showcase of torus canards in neuronal bursters

Rapid action potential generation --- spiking --- and alternating intervals of spiking and quiescence --- bursting --- are two dynamic patterns observed in neuronal activity. In computational models of neuronal systems, the transition from spiking to bursting often exhibits complex bifurcation structure. One type of transition involves the torus canard, which was originally observed in a simple biophysical model of a Purkinje cell. In this article, we expand on that original result by showing that torus canards arise in a broad array of well-known computational neuronal models with three different classes of bursting dynamics: sub-Hopf/fold cycle bursting, circle/fold cycle bursting, and fold/fold cycle bursting. The essential features that these models share are multiple time scales leading naturally to decomposition into slow and fast systems, a saddle-node of periodic orbits in the fast system, and a torus bifurcation in the full system. We show that the transition from spiking to bursting in each model system is given by an explosion of torus canards. Based on these examples, as well as on emerging theory, we propose that torus canards are a common dynamic phenomenon separating the regimes of spiking and bursting activity.

preprint2011arXiv

Codimension-two homoclinic bifurcations underlying spike adding in the Hindmarsh-Rose burster

The well-studied Hindmarsh-Rose model of neural action potential is revisited from the point of view of global bifurcation analysis. This slow-fast system of three paremeterised differential equations is arguably the simplest reduction of Hodgkin-Huxley models capable of exhibiting all qualitatively important distinct kinds of spiking and bursting behaviour. First, keeping the singular perturbation parameter fixed, a comprehensive two-parameter bifurcation diagram is computed by brute force. Of particular concern is the parameter regime where lobe-shaped regions of irregular bursting undergo a transition to stripe-shaped regions of periodic bursting. The boundary of each stripe represents a fold bifurcation that causes a smooth spike-adding transition where the number of spikes in each burst is increased by one. Next, numerical continuation studies reveal that the global structure is organised by various curves of homoclinic bifurcations. In particular the lobe to stripe transition is organised by a sequence of codimension-two orbit- and inclination-flip points that occur along {\em each} homoclinic branch. Each branch undergoes a sharp turning point and hence approximately has a double-cover of the same curve in parameter space. The sharp turn is explained in terms of the interaction between a two-dimensional unstable manifold and a one-dimensional slow manifold in the singular limit. Finally, a new local analysis is undertaken using approximate Poincaré maps to show that the turning point on each homoclinic branch in turn induces an inclination flip that gives birth to the fold curve that organises the spike-adding transition. Implications of this mechanism for explaining spike-adding behaviour in other excitable systems are discussed.