Closed-form modeling of neuronal spike train statistics using multivariate Hawkes cumulants

We derive exact analytical expressions for the cumulants of any orders of neuronal membrane potentials driven by spike trains in a multivariate Hawkes process model with excitation and inhibition. Such expressions can be used for the prediction and sensitivity analysis of the statistical behavior of the model over time and to estimate the probability densities of neuronal membrane potentials using Gram-Charlier expansions. Our results are shown to provide a better alternative to Monte Carlo estimates via stochastic simulations and computer codes based on combinatorial recursions are included.

Fluid limit theorems for stochastic hybrid systems with application to neuron models

In this paper we establish limit theorems for a class of stochastic hybrid systems (continuous deterministic dynamics coupled with jump Markov processes) in the fluid limit (small jumps at high frequency), thus extending known results for jump Markov processes. We prove a functional law of large numbers with exponential convergence speed, derive a diffusion approximation, and establish a functional central limit theorem. We apply these results to neuron models with stochastic ion channels, as the number of channels goes to infinity, estimating the convergence to the deterministic model. In terms of neural coding, we apply our central limit theorems to numerically estimate the impact of channel noise both on frequency and spike timing coding.

Averaging for a Fully Coupled Piecewise-Deterministic Markov Process in Infinite Dimensions

In this paper we consider the generalized Hodgkin-Huxley model introduced in Austin (2008). This model describes the propagation of an action potential along the axon of a neuron at the scale of ion channels. Mathematically, this model is a fully coupled piecewise-deterministic Markov process (PDMP) in infinite dimensions. We introduce two time scales in this model in considering that some ion channels open and close at faster jump rates than others. We perform a slow-fast analysis of this model and prove that, asymptotically, this ‘two-time-scale’ model reduces to the so-called averaged model, which is still a PDMP in infinite dimensions, for which we provide effective evolution equations and jump rates.

Intrinsic variability of latency to first-spike

The assessment of the variability of neuronal spike timing is fundamental to gain understanding of latency coding. Based on recent mathematical results, we investigate theoretically the impact of channel noise on latency variability. For large numbers of ion channels, we derive the asymptotic distribution of latency, together with an explicit expression for its variance. Consequences in terms of information processing are studied with Fisher information in the Morris-Lecar model. A competition between sensitivity to input and precision is responsible for favoring two distinct regimes of latencies.