Analysis of inverse stochastic resonance and long-term firing in Hodgkin-Huxley neurons
Henry Tuckwell and Jürgen Jost
Contact the author: Please use for correspondence this email.
Submission date: 25. Feb. 2012 (revised version: April 2012)
published in: Physica / A, 391 (2012) 22, p. 5311-5325
DOI number (of the published article): 10.1016/j.physa.2012.06.019
with the following different title: Analysis of inverse stochastic resonance and the long-term firing of Hodgkin-Huxley neurons with Gaussian white noise
Keywords and phrases: Hodgkin-Huxley, Stochastic analysis
Download full preprint: PDF (1601 kB)
In order to explain the occurrence of a minimum in firing rate which occurs for certain mean input levels μ as noise level σ increases (inverse stochastic resonance, ISR) in Hodgkin-Huxley (HH) systems, we analyse the underlying transitions from a stable equilibrium point to limit cycle and vice-versa. For a value of μ at which ISR is pronounced, properties of the corresponding stable equilibrium point are found. A linearized approximation around this point has oscillatory solutions from whose maxima spikes tend to occur. A one dimensional diffusion is also constructed for small noise. Properties of the basin of attraction of the limit cycle (spike) are investigated heuristically. Long term trials of duration 500000 ms are carried out for values of σ from 0 to 2.0. The graph of mean spike count versus σ is divided into 4 regions R1,...,R4, where R3 contains the minimum associated with ISR. In R1 transitions to the basin of attraction of the rest point are not observed until a small critical value of σ = σc1 is reached, at the beginning of R2. The sudden decline in firing rate when σ is just greater than σc1 implies that there is only a small range of noise levels 0 < σ < σc1 where repetitive spiking is safe from annihilation by noise. The firing rate remains small throughout R3. At a larger critical value σ = σc2 which signals the beginning of R4, the probability of transitions from the basin of attraction of the equilibrium point to that of the limit cycle apparently becomes greater than zero and the spike rate thereafter increases with increasing σ. The quantitative scheme underlying the ISR curve is outlined in terms of the properties of exit time random variables. In the final subsection, several statistical properties of the main random variables associated with long term spiking activity are given, including distributions of exit times from the two relevant basins of attraction and the interspike interval.