No Arabic abstract
We show that the stochastic Morris-Lecar neuron, in a neighborhood of its stable point, can be approximated by a two-dimensional Ornstein-Uhlenbeck (OU) modulation of a constant circular motion. The associated radial OU process is an example of a leaky integrate-and-fire (LIF) model prior to firing. A new model constructed from a radial OU process together with a simple firing mechanism based on detailed Morris-Lecar firing statistics reproduces the Morris-Lecar Interspike Interval (ISI) distribution, and has the computational advantages of a LIF. The result justifies the large amount of attention paid to the LIF models.
In this paper, we provide a complete mathematical construction for a stochastic leaky-integrate-and-fire model (LIF) mimicking the interspike interval (ISI) statistics of a stochastic FitzHugh-Nagumo neuron model (FHN) in the excitable regime, where the unique fixed point is stable. Under specific types of noises, we prove that there exists a global random attractor for the stochastic FHN system. The linearization method is then applied to estimate the firing time and to derive the associated radial equation representing a LIF equation. This result confirms the previous prediction in [Ditlevsen, S. and Greenwood, P. (2013). The Morris-Lecar neuron model embeds a leaky integrate-and-fire model. Journal of Mathematical Biology, 67(2):239-259] for the Morris-Lecar neuron model in the bistability regime consisting of a stable fixed point and a stable limit cycle.
We study the dynamics of identical leaky integrate-and-fire neurons with symmetric non-local coupling. Upon varying control parameters (coupling strength, coupling range, refractory period) we investigate the systems behaviour and highlight the formation of chimera states. We show that the introduction of a refractory period enlarges the parameter region where chimera states appear and affects the chimera multiplicity.
Neural network dynamics emerge from the interaction of spiking cells. One way to formulate the problem is through a theoretical framework inspired by ideas coming from statistical physics, the so-called mean-field theory. In this document, we investigate different issues related to the mean-field description of an excitatory network made up of leaky integrate-and-fire neurons. The description is written in the form a nonlinear partial differential equation which is known to blow up in finite time when the network is strongly connected. We prove that in a moderate coupling regime the equation is globally well-posed in the space of measures, and that there exist stationary solutions. In the case of weak connectivity we also demonstrate the uniqueness of the steady state and its global exponential stability. The method to show those mathematical results relies on a contraction argument of Doeblins type in the linear case, which corresponds to a population of non-interacting units.
The classical biophysical Morris-Lecar model of neuronal excitability predicts that upon stimulation of the neuron with a sufficiently large constant depolarizing current there exists a finite interval of the current values where periodic spike generation occurs. Above the upper boundary of this interval, there is four-stage damping of the spike amplitude: 1) minor primary damping, which reflects a typical transient to stationary dynamic state, 2) plateau of nearly undamped periodic oscillations, 3) strong damping, and 4) reaching a constant asymptotic value of the neuron potential. We have shown that in the vicinity of the asymptote the Morris-Lecar equations can be reduced to the standard equation for exponentially damped harmonic oscillations. Importantly, all coefficients of this equation can be explicitly expressed through parameters of the original Morris-Lecar model, enabling direct comparison of the numerical and analytical solutions for the neuron potential dynamics at later stages of the spike amplitude damping.
In the mean field integrate-and-fire model, the dynamics of a typical neuron within a large network is modeled as a diffusion-jump stochastic process whose jump takes place once the voltage reaches a threshold. In this work, the main goal is to establish the convergence relationship between the regularized process and the original one where in the regularized process, the jump mechanism is replaced by a Poisson dynamic, and jump intensity within the classically forbidden domain goes to infinity as the regularization parameter vanishes. On the macroscopic level, the Fokker-Planck equation for the process with random discharges (i.e. Poisson jumps) are defined on the whole space, while the equation for the limit process is on the half space. However, with the iteration scheme, the difficulty due to the domain differences has been greatly mitigated and the convergence for the stochastic process and the firing rates can be established. Moreover, we find a polynomial-order convergence for the distribution by a re-normalization argument in probability theory. Finally, by numerical experiments, we quantitatively explore the rate and the asymptotic behavior of the convergence for both linear and nonlinear models.