Source author record

Eva Löcherbach

Eva Löcherbach 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

16works
3topics
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

16 published item(s)

preprint2022arXiv

Fluctuations for mean field limits of interacting systems of spiking neurons

We consider a system of $N$ neurons, each spiking randomly with rate depending on its membrane potential. When a neuron spikes, its potential is reset to $0$ and all other neurons receive an additional amount $h/N$ of potential, where $ h > 0$ is some fixed parameter. In between successive spikes, each neuron's potential undergoes some leakage at constant rate $ α. $ While the propagation of chaos of the system, as $N \to \infty$, to a limit nonlinear jumping stochastic differential equation has already been established in a series of papers, see De Masi et al. (2015) and Fournier and Löcherbach (2016), the present paper is devoted to the associated central limit theorem. More precisely we study the measure valued process of fluctuations at scale $ N^{-1/2}$ of the empirical measures of the membrane potentials, centered around the associated limit. We show that this fluctuation process, interpreted as càdlàg process taking values in a suitable weighted Sobolev space, converges in law to a limit process characterized by a system of stochastic differential equations driven by Gaussian white noise. We complete this picture by studying the fluctuations, at scale $ N^{-1/2}, $ of a fixed number of membrane potential processes around their associated limit quantities, giving rise to a mesoscopic approximation of the membrane potentials that take into account the correlations within the finite system.

preprint2022arXiv

Kalikow decomposition for counting processes with stochastic intensity and application to simulation algorithms

We propose a new Kalikow decomposition for continuous time multivariate counting processes, on potentially infinite networks. We prove the existence of such a decomposition in various cases. This decomposition allows us to derive simulation algorithms that hold either for stationary processes with potentially infinite network but bounded intensities, or for processes with unbounded intensities in a finite network and with empty past before 0. The Kalikow decomposition is not unique and we discuss the choice of the decomposition in terms of algorithmic efficiency in certain cases. We apply these methods on several examples: linear Hawkes process, age dependent Hawkes process, exponential Hawkes process, Galves-Löcherbach process.

preprint2022arXiv

On a finite-size neuronal population equation

Population equations for infinitely large networks of spiking neurons have a long tradition in theoretical neuroscience. In this work, we analyze a recent generalization of these equations to populations of finite size, which takes the form of a nonlinear stochastic integral equation. We prove that, in the case of leaky integrate-and-fire (LIF) neurons with escape noise and for a slightly simplified version of the model, the equation is well-posed and stable in the sense of Brémaud-Massoulié. The proof combines methods from Markov processes taking values in the space of positive measures and nonlinear Hawkes processes. For applications, we also provide efficient simulation algorithms.

preprint2021arXiv

White-noise driven conditional McKean-Vlasov limits for systems of particles with simultaneous and random jumps

We study the convergence of $N-$particle systems described by SDEs driven by Brownian motion and Poisson random measure, where the coefficients depend on the empirical measure of the system. Every particle jumps with a jump rate depending on its position and on the empirical measure of the system. Jumps are simultaneous, that is, at each jump time, all particles of the system are affected by this jump and receive a random jump height that is centred and scaled in $N^{-1/2}.$ This particular scaling implies that the limit of the empirical measures of the system is random, describing the conditional distribution of one particle in the limit system. We call such limits {\it conditional McKean-Vlasov limits}. The conditioning in the limit measure reflects the dependencies between coexisting particles in the limit system such that we are dealing with a {\it conditional propagation of chaos property}. As a consequence of the scaling in $N^{-1/2}$ and of the fact that the limit of the empirical measures is not deterministic the limit system turns out to be solution of a non-linear SDE, where not independent martingale measures and white noises appear having an intensity that depends on the conditional law of the process.

preprint2020arXiv

Conditional propagation of chaos for mean field systems of interacting neurons

We study the stochastic system of interacting neurons introduced in De Masi et al. (2015) and in Fournier and Löcherbach (2016) in a diffusive scaling. The system consists of $N$ neurons, each spiking randomly with rate depending on its membrane potential. At its spiking time, the potential of the spiking neuron is reset to $0$ and all other neurons receive an additional amount of potential which is a centred random variable of order $ 1 / \sqrt{N}.$ In between successive spikes, each neuron's potential follows a deterministic flow. We prove the convergence of the system, as $N \to \infty$, to a limit nonlinear jumping stochastic differential equation driven by Poisson random measure and an additional Brownian motion $W$ which is created by the central limit theorem. This Brownian motion is underlying each particle's motion and induces a common noise factor for all neurons in the limit system. Conditionally on $W,$ the different neurons are independent in the limit system. This is the {\it conditional propagation of chaos} property. We prove the well-posedness of the limit equation by adapting the ideas of Graham (1992) to our frame. To prove the convergence in distribution of the finite system to the limit system, we introduce a new martingale problem that is well suited for our framework. The uniqueness of the limit is deduced from the exchangeability of the underlying system.

preprint2020arXiv

Exponential ergodicity for diffusions with jumps driven by a Hawkes process

In this paper, we introduce a new class of processes which are diffusions with jumps driven by a multivariate nonlinear Hawkes process. Our goal is to study their long-time behavior. In the case of exponential memory kernels for the underlying Hawkes process we establish conditions for the positive Harris recurrence of the couple (X, Y), where X denotes the diffusion process and Y the piecewise deterministic Markov process (PDMP) defining the stochastic intensity of the driving Hawkes. As a direct consequence of the Harris recurrence, we obtain the ergodic theorem for X. Furthermore, we provide sufficient conditions under which the process is exponentially $β$--mixing.

preprint2020arXiv

Stability for Hawkes processes with inhibition

We consider a multivariate non-linear Hawkes process in a multi-class setup where particles are organised within two populations of possibly different sizes, such that one of the populations acts excitatory on the system while the other population acts inhibitory on the system. The goal of this note is to present a class of Hawkes Processes with stable dynamics without assumptions on the spectral radius of the associated weight function matrix. This illustrates how inhibition in a Hawkes system significantly affects the stability properties of the system.

preprint2019arXiv

A system of interacting neurons with short term synaptic facilitation

In this paper we present a simple microscopic stochastic model describing short term plasticity within a large homogeneous network of interacting neurons. Each neuron is represented by its membrane potential and by the residual calcium concentration within the cell at a given time. Neurons spike at a rate depending on their membrane potential. When spiking, the residual calcium concentration of the spiking neuron increases by one unit. Moreover, an additional amount of potential is given to all other neurons in the system. This amount depends linearly on the current residual calcium concentration within the cell of the spiking neuron. In between successive spikes, the potentials and the residual calcium concentrations of each neuron decrease at a constant rate. We show that in this framework, short time memory can be described as the tendency of the system to keep track of an initial stimulus by staying within a certain region of the space of configurations during a short but macroscopic amount of time before finally being kicked out of this region and relaxing to equilibrium. The main technical tool is a rigorous justification of the passage to a large population limit system and a thorough study of the limit equation.

preprint2016arXiv

Absolute continuity of the invariant measure in Piecewise Deterministic Markov Processes having degenerate jumps

We consider piecewise deterministic Markov processes with degenerate transition kernels of the "house-of-cards"-type. We use a splitting scheme based on jump times to prove the absolute continuity, as well as some regularity, of the invariant measure of the process. Finally, we obtain finer results on the regularity of the one-dimensional marginals of the invariant measure, using integration by parts with respect to the jump times.

preprint2016arXiv

Multi-class oscillating systems of interacting neurons

We consider multi-class systems of interacting nonlinear Hawkes processes modeling several large families of neurons and study their mean field limits. As the total number of neurons goes to infinity we prove that the evolution within each class can be described by a nonlinear limit differential equation driven by a Poisson random measure, and state associated central limit theorems. We study situations in which the limit system exhibits oscillatory behavior, and relate the results to certain piecewise deterministic Markov processes and their diffusion approximations.

preprint2016arXiv

Non-parametric estimation of the spiking rate in systems of interacting neurons

We consider a model of interacting neurons where the membrane potentials of the neurons are described by a multidimensional piecewise deterministic Markov process (PDMP) with values in ${\mathbb R}^N, $ where $ N$ is the number of neurons in the network. A deterministic drift attracts each neuron's membrane potential to an equilibrium potential $m.$ When a neuron jumps, its membrane potential is reset to $0,$ while the other neurons receive an additional amount of potential $\frac{1}{N}.$ We are interested in the estimation of the jump (or spiking) rate of a single neuron based on an observation of the membrane potentials of the $N$ neurons up to time $t.$ We study a Nadaraya-Watson type kernel estimator for the jump rate and establish its rate of convergence in $L^2 .$ This rate of convergence is shown to be optimal for a given Hölder class of jump rate functions. We also obtain a central limit theorem for the error of estimation. The main probabilistic tools are the uniform ergodicity of the process and a fine study of the invariant measure of a single neuron.

preprint2015arXiv

Hawkes processes with variable length memory and an infinite number of components

In this paper, we build a model for biological neural nets where the activity of the network is described by Hawkes processes having a variable length memory. The particularity of this paper is to deal with an infinite number of components. We propose a graphical construction of the process and we build, by means of a perfect simulation algorithm, a stationary version of the process. To carry out this algorithm, we make use of a Kalikow-type decomposition technique. Two models are described in this paper. In the first model, we associate to each edge of the interaction graph a saturation threshold that controls the influence of a neuron on another. In the second model, we impose a structure on the interaction graph leading to a cascade of spike trains. Such structures, where neurons are divided into layers can be found in retina.

preprint2015arXiv

On a toy model of interacting neurons

We continue the study of a stochastic system of interacting neurons introduced in De Masi-Galves-Löcherbach-Presutti (2014). The system consists of N neurons, each spiking randomly with rate depending on its membrane potential. At its spiking time, the neuron potential is reset to 0 and all other neurons receive an additional amount 1/N of potential. Moreover, electrical synapses induce a deterministic drift of the system towards its center of mass. We prove propagation of chaos of the system, as N tends to infinity, to a limit nonlinear jumping stochastic differential equation. We consequently improve on the results of De Masi-Galves-Löcherbach-Presutti (2014), since (i) we remove the compact support condition on the initial datum, (ii) we get a rate of convergence in $1/\sqrt N$. Finally, we study the limit equation: we describe the shape of its time-marginals, we prove the existence of a unique non-trivial invariant distribution, we show that the trivial invariant distribution is not attractive, and in a special case, we establish the convergence to equilibrium.

preprint2014arXiv

Hydrodynamic limit for interacting neurons

This paper studies the hydrodynamic limit of a stochastic process describing the time evolution of a system with N neurons with mean-field interactions produced both by chemical and by electrical synapses. This system can be informally described as follows. Each neuron spikes randomly following a point process with rate depending on its membrane potential. At its spiking time, the membrane potential of the spiking neuron is reset to the value 0 and, simultaneously, the membrane potentials of the other neurons are increased by an amount of potential 1/N . This mimics the effect of chemical synapses. Additionally, the effect of electrical synapses is represented by a deterministic drift of all the membrane potentials towards the average value of the system. We show that, as the system size N diverges, the distribution of membrane potentials becomes deterministic and is described by a limit density which obeys a non linear PDE which is a conservation law of hyperbolic type.

preprint2013arXiv

Infinite systems of interacting chains with memory of variable length - a stochastic model for biological neural nets

We consider a new class of non Markovian processes with a countable number of interacting components. At each time unit, each component can take two values, indicating if it has a spike or not at this precise moment. The system evolves as follows. For each component, the probability of having a spike at the next time unit depends on the entire time evolution of the system after the last spike time of the component. This class of systems extends in a non trivial way both the interacting particle systems, which are Markovian, and the stochastic chains with memory of variable length which have finite state space. These features make it suitable to describe the time evolution of biological neural systems. We construct a stationary version of the process by using a probabilistic tool which is a Kalikow-type decomposition either in random environment or in space-time. This construction implies uniqueness of the stationary process. Finally we consider the case where the interactions between components are given by a critical directed Erdös-Rényi-type random graph with a large but finite number of components. In this framework we obtain an explicit upper-bound for the correlation between successive inter-spike intervals which is compatible with previous empirical findings.

preprint2012arXiv

Transition densities for stochastic Hodgkin-Huxley models

We consider a stochastic Hodgkin-Huxley model driven by a periodic signal as model for the membrane potential of a pyramidal neuron. The associated five dimensional diffusion process is a time inhomogeneous highly degenerate diffusion for which the weak Hoermander condition holds only locally. Using a technique which is based on estimates of the Fourier transform, inspired by Fournier 2008, Bally 2007 and De Marco 2011, we show that the process admits locally a strictly positive continuous transition density. Moreover, we show that the presence of noise enables the stochastic system to imitate any possible deterministic spiking behavior, i.e. mixtures of regularly spiking and non-spiking time periods are possible features of the stochastic model. This is a fundamental difference between stochastic and deterministic Hodgkin-Huxley models.