A leaky integrate-and-fire model with adaptation for the generation of a spike train.

A model is proposed to describe the spike-frequency adaptation observed in many neuronal systems. We assume that adaptation is mainly due to a calcium-activated potassium current, and we consider two coupled stochastic differential equations for which an analytical approach combined with simulation techniques and numerical methods allow to obtain both qualitative and quantitative results about asymptotic mean firing rate, mean calcium concentration and the firing probability density. A related algorithm, based on the Hazard Rate Method, is also devised and described.

A simple algorithm to generate firing times for leaky integrate-and-fire neuronal model.

A method to generate first passage times for a class of stochastic processes is proposed. It does not require construction of the trajectories as usually needed in simulation studies, but is based on an integral equation whose unknown quantity is the probability density function of the studied first passage times and on the application of the hazard rate method. The proposed procedure is particularly efficient in the case of the Ornstein-Uhlenbeck process, which is important for modeling spiking neuronal activity.

A non-autonomous stochastic predator-prey model.

The aim of this paper is to consider a non-autonomous predator-prey-like system, with a Gompertz growth law for the prey. By introducing random variations in both prey birth and predator death rates, a stochastic model for the predator-prey-like system in a random environment is proposed and investigated. The corresponding Fokker-Planck equation is solved to obtain the joint probability density for the prey and predator populations and the marginal probability densities. The asymptotic behavior of the predator-prey stochastic model is also analyzed.

Gauss-diffusion processes for modeling the dynamics of a couple of interacting neurons.

With the aim to describe the interaction between a couple of neurons a stochastic model is proposed and formalized. In such a model, maintaining statements of the Leaky Integrate-and-Fire framework, we include a random component in the synaptic current, whose role is to modify the equilibrium point of the membrane potential of one of the two neurons and when a spike of the other one occurs it is turned on. The initial and after spike reset positions do not allow to identify the inter-spike intervals with the corresponding first passage times. However, we are able to apply some well-known results for the first passage time problem for the Ornstein-Uhlenbeck process in order to obtain (i) an approximation of the probability density function of the inter-spike intervals in one-way-type interaction and (ii) an approximation of the tail of the probability density function of the inter-spike intervals in the mutual interaction. Such an approximation is admissible for small instantaneous firing rates of both neurons.

On the evaluation of firing densities for periodically driven neuron models.

The leaky integrate-and-fire model for neuronal spiking events driven by a periodic stimulus is studied by using the Fokker-Planck formulation. To this purpose, an essential use is made of the asymptotic behavior of the first-passage-time probability density function of a time homogeneous diffusion process through an asymptotically periodic threshold. Numerical comparisons with some recently published results derived by a different approach are performed. Use of a new asymptotic approximation is then made in order to design a numerical algorithm of predictor-corrector type to solve the integral equation in the unknown first-passage-time probability density function. Such algorithm, characterized by a reduced (linear) computation time, is seen to provide a high computation accuracy. Finally, it is shown that such an approach yields excellent approximations to the firing probability density function for a wide range of parameters, including the case of high stimulus frequencies.