Source author record

Jonathan Touboul

Jonathan Touboul 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

36works
20topics
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

36 published item(s)

preprint2021arXiv

The Symmetry Basis of Pattern Formation in Reaction-Diffusion Networks

In networks of nonlinear oscillators, symmetries place hard constraints on the system that can be exploited to predict universal dynamical features and steady-states, providing a rare generic organizing principle for far-from-equilibrium systems. However, the robustness of this class of theories to symmetry-disrupting imperfections is untested. Here, we develop a model experimental reaction-diffusion network of chemical oscillators to test applications of this theory in the context of self-organizing systems relevant to biology and soft robotics. The network is a ring of 4 identical microreactors containing the oscillatory Belousov-Zhabotinsky reaction coupled to nearest neighbors via diffusion. Assuming perfect symmetry, theory predicts 4 categories of stable spatiotemporal phase-locked periodic states and 4 categories of invariant manifolds that guide and structure transitions between phase-locked states. In our experiments, we observed the predicted symmetry-derived synchronous clustered transients that occur when the dynamical trajectories coincide with invariant manifolds. However, we observe only 3 of the 4 phase-locked states that are predicted for the idealized homogeneous system. Quantitative agreement between experiment and numerical simulations is found by accounting for the small amount of experimentally determined heterogeneity. This work demonstrates that a surprising degree of the network's dynamics are constrained by symmetry in spite of the breakdown of the assumption of homogeneity and raises the question of why heterogeneity destabilizes some symmetry predicted states, but not others.

preprint2017arXiv

Large deviations for randomly connected neural networks: I. Spatially extended systems

In a series of two papers, we investigate the large deviations and asymptotic behavior of stochastic models of brain neural networks with random interaction coefficients. In this first paper, we take into account the spatial structure of the brain and consider (i) the presence of interaction delays that depend on the distance between cells and (ii) Gaussian random interaction amplitude whose mean and variance depend on the neurons positions and scale as the inverse of the network size. We show that the empirical measure satisfies a large-deviation principle with good rate function reaching its minimum at a unique spatially extended probability measure. This result implies averaged convergence of the empirical measure and propagation of chaos. The limit is characterized through complex non-Markovian implicit equation in which the network interaction term is replaced by a non-local Gaussian process whose statistics depend on the solution over the whole neural field.

preprint2017arXiv

Large deviations for randomly connected neural networks: II. State-dependent interactions

This work continues the analysis of large deviations for randomly connected neural networks models of the brain. The originality of the model relies on the fact that the directed impact of one particle onto another depends on the state of both particles, and (ii) have random Gaussian amplitude with mean and variance scaling as the inverse of the network size. Similarly to the spatially extended case, we show that under sufficient regularity assumptions, the empirical measure satisfies a large-deviation principle with good rate function achieving its minimum at a unique probability measure, implying in particular its convergence in both averaged and quenched cases, as well as a propagation of chaos property (in the averaged case only). The class of model we consider notably includes a stochastic version of Kuramoto model with random connections.

preprint2017arXiv

Wild oscillations in a nonlinear neuron model with resets: (II) Mixed-mode oscillations

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.

preprint2016arXiv

Propagation of chaos in neural fields

We consider the problem of the limit of bio-inspired spatially extended neuronal networks including an infinite number of neuronal types (space locations), with space-dependent propagation delays modeling neural fields. The propagation of chaos property is proved in this setting under mild assumptions on the neuronal dynamics, valid for most models used in neuroscience, in a mesoscopic limit, the neural-field limit, in which we can resolve the quite fine structure of the neuron's activity in space and where averaging effects occur. The mean-field equations obtained are of a new type: they take the form of well-posed infinite-dimensional delayed integro-differential equations with a nonlocal mean-field term and a singular spatio-temporal Brownian motion. We also show how these intricate equations can be used in practice to uncover mathematically the precise mesoscopic dynamics of the neural field in a particular model where the mean-field equations exactly reduce to deterministic nonlinear delayed integro-differential equations. These results have several theoretical implications in neuroscience we review in the discussion.

preprint2016arXiv

Real eigenvalues of non-symmetric random matrices: Transitions and Universality

In the past 20 years, the study of real eigenvalues of non-symmetric real random matrices has seen important progress. Notwithstanding, central questions still remain open, such as the characterization of their asymptotic statistics and the universality thereof. In this letter we show that for a wide class of matrices, the number $k_n$ of real eigenvalues of a matrix of size $n$ is asymptotically Gaussian with mean $\bar k_n=\mathcal{O}(\sqrt{n})$ and variance $\bar k_n(2-\sqrt{2})$. Moreover, we show that the limit distribution of real eigenvalues undergoes a transition between bimodal for $k_n=o(\sqrt{n})$ to unimodal for $k_n=\mathcal{O}(n)$, with a uniform distribution at the transition. We predict theoretically these behaviours in the Ginibre ensemble using a log-gas approach, and show numerically that they hold for a wide range of random matrices with independent entries beyond the universality class of the circular law.

preprint2015arXiv

Absorption properties of stochastic equations with Hölder diffusion coefficients

In this article, we address the absorption properties of a class of stochastic differ- ential equations around singular points where both the drift and diffusion functions vanish. According to the Hölder coefficient alpha of the diffusion function around the singular point, we identify different regimes. Stability of the absorbing state, large deviations for the absorption time, existence of stationary or quasi-stationary distributions are discussed. In particular, we show that quasi-stationary distributions only exist for alpha < 3/4, and for alpha in the interval (3/4, 1), no quasi-stationary distribution is found and numerical simulations tend to show that the process conditioned on not being absorbed initiates an almost sure exponential convergence towards the absorbing state (as is demonstrated to be true for alpha = 1). Applications of these results to stochastic bifurcations are discussed.

preprint2015arXiv

Complex oscillations in the delayed Fitzhugh-Nagumo equation

Motivated by the dynamics of neuronal responses, we analyze the dynamics of the Fitzhugh-Nagumo slow-fast system with delayed self-coupling. This system provides a canonical example of a canard explosion for sufficiently small delays. Beyond this regime, delays significantly enrich the dynamics, leading to mixed-mode oscillations, bursting and chaos. These behaviors emerge from a delay-induced subcritical Bogdanov-Takens instability arising at the fold points of the S-shaped critical manifold. Underlying the transition from canard-induced to delay-induced dynamics is an abrupt switch in the nature of the Hopf bifurcation.

preprint2015arXiv

Limits and dynamics of stochastic neuronal networks with random heterogeneous delays

Realistic networks display heterogeneous transmission delays. We analyze here the limits of large stochastic multi-populations networks with stochastic coupling and random interconnection delays. We show that depending on the nature of the delays distributions, a quenched or averaged propagation of chaos takes place in these networks, and that the network equations converge towards a delayed McKean-Vlasov equation with distributed delays. Our approach is mostly fitted to neuroscience applications. We instantiate in particular a classical neuronal model, the Wilson and Cowan system, and show that the obtained limit equations have Gaussian solutions whose mean and standard deviation satisfy a closed set of coupled delay differential equations in which the distribution of delays and the noise levels appear as parameters. This allows to uncover precisely the effects of noise, delays and coupling on the dynamics of such heterogeneous networks, in particular their role in the emergence of synchronized oscillations. We show in several examples that not only the averaged delay, but also the dispersion, govern the dynamics of such networks.

preprint2015arXiv

Noise-induced behaviors in neural mean field dynamics

The collective behavior of cortical neurons is strongly affected by the presence of noise at the level of individual cells. In order to study these phenomena in large-scale assemblies of neurons, we consider networks of firing-rate neurons with linear intrinsic dynamics and nonlinear coupling, belonging to a few types of cell populations and receiving noisy currents. Asymptotic equations as the number of neurons tends to infinity (mean field equations) are rigorously derived based on a probabilistic approach. These equations are implicit on the probability distribution of the solutions which generally makes their direct analysis difficult. However, in our case, the solutions are Gaussian, and their moments satisfy a closed system of nonlinear ordinary differential equations (ODEs), which are much easier to study than the original stochastic network equations, and the statistics of the empirical process uniformly converge towards the solutions of these ODEs. Based on this description, we analytically and numerically study the influence of noise on the collective behaviors, and compare these asymptotic regimes to simulations of the network. We observe that the mean field equations provide an accurate description of the solutions of the network equations for network sizes as small as a few hundreds of neurons. In particular, we observe that the level of noise in the system qualitatively modifies its collective behavior, producing for instance synchronized oscillations of the whole network, desynchronization of oscillating regimes, and stabilization or destabilization of stationary solutions. These results shed a new light on the role of noise in shaping collective dynamics of neurons, and gives us clues for understanding similar phenomena observed in biological networks.

preprint2015arXiv

On a kinetic FitzHugh-Nagumo model of neuronal network

We investigate existence and uniqueness of solutions of a McKean-Vlasov evolution PDE representing the macroscopic behaviour of interacting Fitzhugh-Nagumo neurons. This equation is hypoelliptic, nonlocal and has unbounded coefficients. We prove existence of a solution to the evolution equation and non trivial stationary solutions. Moreover, we demonstrate uniqueness of the stationary solution in the weakly nonlinear regime. Eventually, using a semigroup factorisation method, we show exponential nonlinear stability in the small connectivity regime.

preprint2015arXiv

On the Dynamic Interplay between Positive and Negative Affects

Emotional disorders and psychological flourishing are the result of complex interactions between positive and negative affects that depend on external events and the subject's internal representations. Based on psychological data, we mathematically model the dynamical balance between positive and negative affects as a function of the response to external positive and negative events. This modeling allows the investigation of the relative impact of two leading forms of therapy on affect balance. The model uses a delay differential equation to analytically study the complete bifurcation diagram of the system. We compare the results of the model to psychological data on a single, recurrently depressed patient that was administered the two types of therapies considered (viz., coping-focused vs. affect-focused). The model leads to the prediction that stabilization at a normal state may rely on evaluating one's emotional state through an historical ongoing emotional state rather than in a narrow present window. The simple mathematical model proposed here offers a theoretically grounded quantitative framework for investigating the temporal process of change and parameters of resilience to relapse.

preprint2015arXiv

The real Ginibre ensemble with $k = O(n)$ real eigenvalues

We consider the ensemble of Real Ginibre matrices with a positive fraction $α>0$ of real eigenvalues. We demonstrate a large deviations principle for the joint eigenvalue density of such matrices and we introduce a two phase log-gas whose stationary distribution coincides with the spectral measure of the ensemble. Using these tools we provide an asymptotic expansion for the probability $p^n_{αn}$ that an $n\times n$ Ginibre matrix has $k=αn$ real eigenvalues and we characterize the spectral measures of these matrices.

preprint2014arXiv

Competition and boundary formation in heterogeneous media: Application to neuronal differentiation

We analyze an inhomogeneous system of coupled reaction-diffusion equations representing the dynamics of gene expression during differentiation of nerve cells. The outcome of this developmental phase is the formation of distinct functional areas separated by sharp and smooth boundaries. It proceeds through the competition between the expression of two genes whose expression is driven by monotonic gradients of chemicals, and the products of gene expression undergo local diffusion and drive gene expression in neighboring cells. The problem therefore falls in a more general setting of species in competition within a non-homogeneous medium. We show that in the limit of arbitrarily small diffusion, there exists a unique monotonic stationary solution, which splits the neural tissue into two winner-take-all parts at a precise boundary point: on both sides of the boundary, different neuronal types are present. In order to further characterize the location of this boundary, we use a blow-up of the system and define a traveling wave problem parametrized by the position within the monotonic gradient: the precise boundary location is given by the unique point in space at which the speed of the wave vanishes.

preprint2014arXiv

Index Distribution of the Ginibre Ensemble

Complex systems, and in particular random neural networks, are often described by randomly interacting dynamical systems with no specific symmetry. In that context, characterizing the number of relevant directions necessitates fine estimates on the Ginibre ensemble. In this Letter, we compute analytically the probability distribution of the number of eigenvalues $N_R$ with modulus greater than $R$ (the index) of a large $N\times N$ random matrix in the real or complex Ginibre ensemble. We show that the fraction $N_R/N=p$ has a distribution scaling as $\exp(-βN^2 ψ_R(p))$ with $β=1$ (respectively $β=1/2$) for the complex (resp. real) Ginibre ensemble. For any $p\in[0,1]$, the equilibrium spectral densities as well as the rate function $ψ_R(p)$ are explicitly derived. This function displays a third order phase transition at the critical (minimum) value $p^*_R=1-R^2$, associated to a phase transition of the Coulomb gas. We deduce that, in the central regime, the fluctuations of the index $N_R$ around its typical value $p^*_R N$ scale as $N^{1/3}$.

preprint2014arXiv

Limits and dynamics of randomly connected neuronal networks

Networks of the brain are composed of a very large number of neurons connected through a random graph and interacting after random delays that both depend on the anatomical distance between cells. In order to comprehend the role of these random architectures on the dynamics of such networks, we analyze the mesoscopic and macroscopic limits of networks with random correlated connectivity weights and delays. We address both averaged and quenched limits, and show propagation of chaos and convergence to a complex integral McKean-Vlasov equations with distributed delays. We then instantiate a completely solvable model illustrating the role of such random architectures in the emerging macroscopic activity. We particularly focus on the role of connectivity levels in the emergence of periodic solutions.

preprint2014arXiv

Pulsatile localized dynamics in delayed neural-field equations in arbitrary dimension

Neural field equations are integro-differential systems describing the macroscopic activity of spatially extended pieces of cortex. In such cortical assemblies, the propagation of information and the transmission machinery induce communication delays, due to the transport of information (propagation delays) and to the synaptic machinery (constant delays). We investigate the role of these delays on the formation of structured spatiotemporal patterns for these systems in arbitrary dimensions. We focus on localized activity, either induced by the presence of a localized stimulus (pulses) or by transitions between two levels of activity (fronts). Linear stability analysis allows to reveal the existence of Hopf bifurcation curves induced by the delays, along different modes that may be symmetric or asymmetric. We show that instabilities strongly depend on the dimension, and in particular may exhibit transversal instabilities along invariant directions. These instabilities yield pulsatile localized activity, and depending on the symmetry of the destabilized modes, either produce spatiotemporal breathing or sloshing patterns.

preprint2014arXiv

Spatially extended networks with singular multi-scale connectivity patterns

The cortex is a very large network characterized by a complex connectivity including at least two scales: a microscopic scale at which the interconnections are non-specific and very dense, while macroscopic connectivity patterns connecting different regions of the brain at larger scale are extremely sparse. This motivates to analyze the behavior of networks with multiscale coupling, in which a neuron is connected to its $v(N)$ nearest-neighbors where $v(N)=o(N)$, and in which the probability of macroscopic connection between two neurons vanishes. These are called singular multi-scale connectivity patterns. We introduce a class of such networks and derive their continuum limit. We show convergence in law and propagation of chaos in the thermodynamic limit. The limit equation obtained is an intricate non-local McKean-Vlasov equation with delays which is universal with respect to the type of micro-circuits and macro-circuits involved.

preprint2014arXiv

The heterogeneous gas with singular interaction: Generalized circular law and heterogeneous renormalized energy

We introduce and analyze $d$ dimensional Coulomb gases with random charge distribution and general external confining potential. We show that these gases satisfy a large deviations principle. The analysis of the minima of the rate function (which is the leading term of the energy) reveals that at equilibrium, the particle distribution is a generalized circular law (i.e. with spherical support but non-necessarily uniform distribution). In the classical electrostatic external potential, there are infinitely many minimizers of the rate function. The most likely macroscopic configuration is a disordered distribution in which particles are uniformly distributed (for $d=2$, the circular law), and charges are independent of the positions of the particles. General charge-dependent confining potentials unfold this degenerate situation: in contrast, the particle density is not uniform, and particles spontaneously organize according to their charge. In that picture the classical electrostatic potential appears as a transition at which order is lost. Sub-leading terms of the energy are derived: we show that these are related to an operator, generalizing the Coulomb renormalized energy, which incorporates the heterogeneous nature of the charges. This heterogeneous renormalized energy informs us about the microscopic arrangements of the particles, which are non-standard, strongly depending on the charges, and include progressive and irregular lattices.

preprint2013arXiv

Large deviations, dynamics and phase transitions in large stochastic heterogeneous neural networks

We analyze the macroscopic behavior of multi-populations randomly connected neural networks with interaction delays. Similar to cases occurring in spin glasses, we show that the sequences of empirical measures satisfy a large deviation principle, and converge towards a self-consistent non-Markovian process. The proof differs in that we are working in infinite-dimensional spaces (interaction delays), non-centered interactions and multiple cell types. The limit equation is qualitatively analyzed, and we identify a number of phase transitions in such systems upon changes in delays, connectivity patterns and dispersion, particularly focusing on the emergence of non-equilibrium states involving synchronized oscillations.

preprint2013arXiv

Macroscopic equations governing noisy spiking neuronal populations

At functional scales, cortical behavior results from the complex interplay of a large number of excitable cells operating in noisy environments. Such systems resist to mathematical analysis, and computational neurosciences have largely relied on heuristic partial (and partially justified) macroscopic models, which successfully reproduced a number of relevant phenomena. The relationship between these macroscopic models and the spiking noisy dynamics of the underlying cells has since then been a great endeavor. Based on recent mean-field reductions for such spiking neurons, we present here {a principled reduction of large biologically plausible neuronal networks to firing-rate models, providing a rigorous} relationship between the macroscopic activity of populations of spiking neurons and popular macroscopic models, under a few assumptions (mainly linearity of the synapses). {The reduced model we derive consists of simple, low-dimensional ordinary differential equations with} parameters and {nonlinearities derived from} the underlying properties of the cells, and in particular the noise level. {These simple reduced models are shown to reproduce accurately the dynamics of large networks in numerical simulations}. Appropriate parameters and functions are made available {online} for different models of neurons: McKean, Fitzhugh-Nagumo and Hodgkin-Huxley models.

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

Synchronization in random balanced networks

Characterizing the in uence of network properties on the global emerging behavior of interacting elements constitutes a central question in many areas, from physical to social sciences. In this article we study a primary model of disordered neuronal networks with excitatory-inhibitory structure and balance constraints. We show how the interplay between structure and disorder in the connectivity leads to a universal transition from trivial to synchronized stationary or periodic states. This transition cannot be explained only through the analysis of the spectral density of the connectivity matrix. We provide a low dimensional approximation that shows the role of both the structure and disorder in the dynamics.

preprint2013arXiv

Topological and Dynamical Complexity of Random Neural Networks

Random neural networks are dynamical descriptions of randomly interconnected neural units. These show a phase transition to chaos as a disorder parameter is increased. The microscopic mechanisms underlying this phase transition are unknown, and similarly to spin-glasses, shall be fundamentally related to the behavior of the system. In this Letter we investigate the explosion of complexity arising near that phase transition. We show that the mean number of equilibria undergoes a sharp transition from one equilibrium to a very large number scaling exponentially with the dimension on the system. Near criticality, we compute the exponential rate of divergence, called topological complexity. Strikingly, we show that it behaves exactly as the maximal Lyapunov exponent, a classical measure of dynamical complexity. This relationship unravels a microscopic mechanism leading to chaos which we further demonstrate on a simpler class of disordered systems, suggesting a deep and underexplored link between topological and dynamical complexity.

preprint2012arXiv

Controllability of the heat and wave equations and their finite difference approximations by the shape of the domain

In this article we study a controllability problem for a parabolic and a hyperbolic partial differential equations in which the control is the shape of the domain where the equation holds. The quantity to be controlled is the trace of the solution into an open subdomain and at a given time, when the right hand side source term is known. The mapping that associates this trace to the shape of the domain is nonlinear. We show (i) an approximate controllability property for the linearized parabolic problem and (ii) an exact local controllability property for the linearized and the nonlinear equations in the hyperbolic case. We then address the same questions in the context of a finite difference spatial semi-discretization in both the parabolic and hyperbolic problems. In this discretized case again we prove a local controllability result for the parabolic problem, and an exact controllability for the hyperbolic case, applying a local surjectivity theorem together with a unique continuation property of the underlying adjoint discrete system.

preprint2012arXiv

Heterogeneous connections induce oscillations in large scale networks

Realistic large-scale networks display an heterogeneous distribution of connectivity weights, that might also randomly vary in time. We show that depending on the level of heterogeneity in the connectivity coefficients, different qualitative macroscopic and microscopic regimes emerge. We evidence in particular generic transitions from stationary to perfectly periodic phase-locked regimes as the disorder parameter is increased, both in a simple model treated analytically and in a biologically relevant network made of excitable cells.

preprint2011arXiv

A Markovian event-based framework for stochastic spiking neural networks

In spiking neural networks, the information is conveyed by the spike times, that depend on the intrinsic dynamics of each neuron, the input they receive and on the connections between neurons. In this article we study the Markovian nature of the sequence of spike times in stochastic neural networks, and in particular the ability to deduce from a spike train the next spike time, and therefore produce a description of the network activity only based on the spike times regardless of the membrane potential process. To study this question in a rigorous manner, we introduce and study an event-based description of networks of noisy integrate-and-fire neurons, i.e. that is based on the computation of the spike times. We show that the firing times of the neurons in the networks constitute a Markov chain, whose transition probability is related to the probability distribution of the interspike interval of the neurons in the network. In the cases where the Markovian model can be developed, the transition probability is explicitly derived in such classical cases of neural networks as the linear integrate-and-fire neuron models with excitatory and inhibitory interactions, for different types of synapses, possibly featuring noisy synaptic integration, transmission delays and absolute and relative refractory period. This covers most of the cases that have been investigated in the event-based description of spiking deterministic neural networks.

preprint2011arXiv

Finite-size and correlation-induced effects in Mean-field Dynamics

The brain's activity is characterized by the interaction of a very large number of neurons that are strongly affected by noise. However, signals often arise at macroscopic scales integrating the effect of many neurons into a reliable pattern of activity. In order to study such large neuronal assemblies, one is often led to derive mean-field limits summarizing the effect of the interaction of a large number of neurons into an effective signal. Classical mean-field approaches consider the evolution of a deterministic variable, the mean activity, thus neglecting the stochastic nature of neural behavior. In this article, we build upon two recent approaches that include correlations and higher order moments in mean-field equations, and study how these stochastic effects influence the solutions of the mean-field equations, both in the limit of an infinite number of neurons and for large yet finite networks. We introduce a new model, the infinite model, which arises from both equations by a rescaling of the variables and, which is invertible for finite-size networks, and hence, provides equivalent equations to those previously derived models. The study of this model allows us to understand qualitative behavior of such large-scale networks. We show that, though the solutions of the deterministic mean-field equation constitute uncorrelated solutions of the new mean-field equations, the stability properties of limit cycles are modified by the presence of correlations, and additional non-trivial behaviors including periodic orbits appear when there were none in the mean field. The origin of all these behaviors is then explored in finite-size networks where interesting mesoscopic scale effects appear. This study leads us to show that the infinite-size system appears as a singular limit of the network equations, and for any finite network, the system will differ from the infinite system.

preprint2011arXiv

Mean Field description of and propagation of chaos in recurrent multipopulation networks of Hodgkin-Huxley and Fitzhugh-Nagumo neurons

We derive the mean-field equations arising as the limit of a network of interacting spiking neurons, as the number of neurons goes to infinity. The neurons belong to a fixed number of populations and are represented either by the Hodgkin-Huxley model or by one of its simplified version, the Fitzhugh-Nagumo model. The synapses between neurons are either electrical or chemical. The network is assumed to be fully connected. The maximum conductances vary randomly. Under the condition that all neurons initial conditions are drawn independently from the same law that depends only on the population they belong to, we prove that a propagation of chaos phenomenon takes places, namely that in the mean-field limit, any finite number of neurons become independent and, within each population, have the same probability distribution. This probability distribution is solution of a set of implicit equations, either nonlinear stochastic differential equations resembling the McKean-Vlasov equations, or non-local partial differential equations resembling the McKean-Vlasov-Fokker- Planck equations. We prove the well-posedness of these equations, i.e. the existence and uniqueness of a solution. We also show the results of some preliminary numerical experiments that indicate that the mean-field equations are a good representation of the mean activity of a finite size network, even for modest sizes. These experiment also indicate that the McKean-Vlasov-Fokker- Planck equations may be a good way to understand the mean-field dynamics through, e.g., a bifurcation analysis.

preprint2010arXiv

Can power-law scaling and neuronal avalanches arise from stochastic dynamics?

The presence of self-organized criticality in biology is often evidenced by a power-law scaling of event size distributions, which can be measured by linear regression on logarithmic axes. We show here that such a procedure does not necessarily mean that the system exhibits self-organized criticality. We first provide an analysis of multisite local field potential (LFP) recordings of brain activity and show that event size distributions defined as negative LFP peaks can be close to power-law distributions. However, this result is not robust to change in detection threshold, or when tested using more rigorous statistical analyses such as the Kolmogorov-Smirnov test. Similar power-law scaling is observed for surrogate signals, suggesting that power-law scaling may be a generic property of thresholded stochastic processes. We next investigate this problem analytically, and show that, indeed, stochastic processes can produce spurious power-law scaling without the presence of underlying self-organized criticality. However, this power-law is only apparent in logarithmic representations, and does not survive more rigorous analysis such as the Kolmogorov-Smirnov test. The same analysis was also performed on an artificial network known to display self-organized criticality. In this case, both the graphical representations and the rigorous statistical analysis reveal with no ambiguity that the avalanche size is distributed as a power-law. We conclude that logarithmic representations can lead to spurious power-law scaling induced by the stochastic nature of the phenomenon. This apparent power-law scaling does not constitute a proof of self-organized criticality, which should be demonstrated by more stringent statistical tests.

preprint2010arXiv

Multi-Resolution Schauder Approach to Multidimensional Gauss-Markov Processes

The study of multidimensional stochastic processes involves complex computations in intricate functional spaces. In particular, the diffusion processes, which include the practically important Gauss-Markov processes, are ordinarily defined through the theory of stochastic integration. Here, inspired by the Lévy-Cieselski construction of the Wiener process, we propose an alternative representation of multidimensional Gauss-Markov processes as expansions on well-chosen Schauder bases, with independent random coefficients of normal law with zero mean and unitary variance. We thereby offer a natural multi-resolution description of Gauss-Markov processes as limits of the finite-dimensional partial sums of the expansion, that are strongly almost-surely convergent. Moreover, such finite-dimensional random processes constitute an optimal approximation of the process, in the sense of minimizing the associated Dirichlet energy under interpolating constraints. This approach allows simpler treatment in many applied and theoretical fields and we provide a short overview of applications we are currently developing.

preprint2010arXiv

On the simulation of nonlinear bidimensional spiking neuron models

Bidimensional spiking models currently gather a lot of attention for their simplicity and their ability to reproduce various spiking patterns of cortical neurons, and are particularly used for large network simulations. These models describe the dynamics of the membrane potential by a nonlinear differential equation that blows up in finite time, coupled to a second equation for adaptation. Spikes are emitted when the membrane potential blows up or reaches a cutoff value. The precise simulation of the spike times and of the adaptation variable is critical for it governs the spike pattern produced, and is hard to compute accurately because of the exploding nature of the system at the spike times. We thoroughly study the precision of fixed time-step integration schemes for this type of models and demonstrate that these methods produce systematic errors that are unbounded, as the cutoff value is increased, in the evaluation of the two crucial quantities: the spike time and the value of the adaptation variable at this time. Precise evaluation of these quantities therefore involve very small time steps and long simulation times. In order to achieve a fixed absolute precision in a reasonable computational time, we propose here a new algorithm to simulate these systems based on a variable integration step method that either integrates the original ordinary differential equation or the equation of the orbits in the phase plane, and compare this algorithm with fixed time-step Euler scheme and other more accurate simulation algorithms.

preprint2010arXiv

Stochastic firing rate models

We review a recent approach to the mean-field limits in neural networks that takes into account the stochastic nature of input current and the uncertainty in synaptic coupling. This approach was proved to be a rigorous limit of the network equations in a general setting, and we express here the results in a more customary and simpler framework. We propose a heuristic argument to derive these equations providing a more intuitive understanding of their origin. These equations are characterized by a strong coupling between the different moments of the solutions. We analyse the equations, present an algorithm to simulate the solutions of these mean-field equations, and investigate numerically the equations. In particular, we build a bridge between these equations and Sompolinsky and collaborators approach (1988, 1990), and show how the coupling between the mean and the covariance function deviates from customary approaches.

preprint2009arXiv

A constructive mean field analysis of multi population neural networks with random synaptic weights and stochastic inputs

We deal with the problem of bridging the gap between two scales in neuronal modeling. At the first (microscopic) scale, neurons are considered individually and their behavior described by stochastic differential equations that govern the time variations of their membrane potentials. They are coupled by synaptic connections acting on their resulting activity, a nonlinear function of their membrane potential. At the second (mesoscopic) scale, interacting populations of neurons are described individually by similar equations. The equations describing the dynamical and the stationary mean field behaviors are considered as functional equations on a set of stochastic processes. Using this new point of view allows us to prove that these equations are well-posed on any finite time interval and to provide a constructive method for effectively computing their unique solution. This method is proved to converge to the unique solution and we characterize its complexity and convergence rate. We also provide partial results for the stationary problem on infinite time intervals. These results shed some new light on such neural mass models as the one of Jansen and Rit \cite{jansen-rit:95}: their dynamics appears as a coarse approximation of the much richer dynamics that emerges from our analysis. Our numerical experiments confirm that the framework we propose and the numerical methods we derive from it provide a new and powerful tool for the exploration of neural behaviors at different scales.

preprint2008arXiv

Sensitivity to the cutoff value in the quadratic adaptive integrate-and-fire model

The quadratic adaptive integrate-and-fire model (Izhikecih 2003, 2007) is recognized as very interesting for its computational efficiency and its ability to reproduce many behaviors observed in cortical neurons. For this reason it is currently widely used, in particular for large scale simulations of neural networks. This model emulates the dynamics of the membrane potential of a neuron together with an adaptation variable. The subthreshold dynamics is governed by a two-parameter differential equation, and a spike is emitted when the membrane potential variable reaches a given cutoff value. Subsequently the membrane potential is reset, and the adaptation variable is added a fixed value called the spike-triggered adaptation parameter. We show in this note that when the system does not converge to an equilibrium point, both variables of the subthreshold dynamical system blow up in finite time whatever the parameters of the dynamics. The cutoff is therefore essential for the model to be well defined and simulated. The divergence of the adaptation variable makes the system very sensitive to the cutoff: changing this parameter dramatically changes the spike patterns produced. Furthermore from a computational viewpoint, the fact that the adaptation variable blows up and the very sharp slope it has when the spike is emitted implies that the time step of the numerical simulation needs to be very small (or adaptive) in order to catch an accurate value of the adaptation at the time of the spike. It is not the case for the similar quartic (Touboul 2008) and exponential (Brette and Gerstner 2005) models whose adaptation variable does not blow up in finite time, and which are therefore very robust to changes in the cutoff value.