Mean-field equations for stochastic firing-rate neural fields with delays: Derivation and noise-induced transitions

In this manuscript we analyze the collective behavior of mean-field limits of large-scale, spatially extended stochastic neuronal networks with delays. Rigorously, the asymptotic regime of such systems is characterized by a very intricate stochastic delayed integro-differential McKean-Vlasov equation that remain impenetrable, leaving the stochastic collective dynamics of such networks poorly understood. In order to study these macroscopic dynamics, we analyze networks of firing-rate neurons, i.e. with linear intrinsic dynamics and sigmoidal interactions. In that case, we prove that the solution of the mean-field equation is Gaussian, hence characterized by its two first moments, and that these two quantities satisfy a set of coupled delayed integro-differential equations. These equations are similar to usual neural field equations, and incorporate noise levels as a parameter, allowing analysis of noise-induced transitions. We identify through bifurcation analysis several qualitative transitions due to noise in the mean-field limit. In particular, stabilization of spatially homogeneous solutions, synchronized oscillations, bumps, chaotic dynamics, wave or bump splitting are exhibited and arise from static or dynamic Turing-Hopf bifurcations. These surprising phenomena allow further exploring the role of noise in the nervous system.Comment: Updated to the latest version published, and clarified the dependence in space of Brownian motion

Resonant spike propagation in coupled neurons with subthreshold activity

Màster en Biofísica, curs 2006-2007Chemical coupling between neurons is only active when the presynaptic neuron is firing, and thus it does not allow for the propagation of subthreshold activity. Electrical coupling via gap junctions, on the other hand, is also ubiquitous and, due to its diffusive nature, transmits both subthreshold and suprathreshold activity between neurons. We study theoretically the propagation of spikes between two neurons that exhibit subthreshold oscillations, and which are coupled via both chemical synapses and gap junctions. Due to the electrical coupling, the periodic subthreshold activity is synchronized in the two neurons, and affects propagation of spikes in such a way that for certain values of the delay in the synaptic coupling, propagation is not possible. This effect could provide a mechanism for the modulation of information transmission in neuronal networks

Synchronisation in networks of delay-coupled type-I excitable systems

We use a generic model for type-I excitability (known as the SNIPER or SNIC model) to describe the local dynamics of nodes within a network in the presence of non-zero coupling delays. Utilising the method of the Master Stability Function, we investigate the stability of the zero-lag synchronised dynamics of the network nodes and its dependence on the two coupling parameters, namely the coupling strength and delay time. Unlike in the FitzHugh-Nagumo model (a model for type-II excitability), there are parameter ranges where the stability of synchronisation depends on the coupling strength and delay time. One important implication of these results is that there exist complex networks for which the adding of inhibitory links in a small-world fashion may not only lead to a loss of stable synchronisation, but may also restabilise synchronisation or introduce multiple transitions between synchronisation and desynchronisation. To underline the scope of our results, we show using the Stuart-Landau model that such multiple transitions do not only occur in excitable systems, but also in oscillatory ones.Comment: 10 pages, 9 figure

Phase-responsiveness transmission in a network of quadratic integrate-and-fire neurons

In the study of the dynamics of neuronal networks, it is interesting to see how the interaction between neurons can elicit different behaviours in each individual one. Moreover, this can lead to the population exhibiting collective phenomena that is not intrinsic to a single cell, such as synchronization. In this project, we work with a large-scale network and a firing-rate model of quadratic integrate-and-fire (QIF) neurons. After studying the dynamics of the QIF model and computing its phase response curve (PRC), we propose an algorithm to describe the population through the PRCs. Our method is able to replicate the same dynamics we observe with the aforementioned models and it also serves us to gain more insight into the transmission of pulses and to explain how a network can maintain a state of synchronized firing

Dynamics of FitzHugh-Nagumo excitable systems with delayed coupling

Small lattices of $N$ nearest neighbor coupled excitable FitzHugh-Nagumo systems, with time-delayed coupling are studied, and compared with systems of FitzHugh-Nagumo oscillators with the same delayed coupling. Bifurcations of equilibria in N=2 case are studied analytically, and it is then numerically confirmed that the same bifurcations are relevant for the dynamics in the case $N>2$. Bifurcations found include inverse and direct Hopf and fold limit cycle bifurcations. Typical dynamics for different small time-lags and coupling intensities could be excitable with a single globally stable equilibrium, asymptotic oscillatory with symmetric limit cycle, bi-stable with stable equilibrium and a symmetric limit cycle, and again coherent oscillatory but non-symmetric and phase-shifted. For an intermediate range of time-lags inverse sub-critical Hopf and fold limit cycle bifurcations lead to the phenomenon of oscillator death. The phenomenon does not occur in the case of FitzHugh-Nagumo oscillators with the same type of coupling.Comment: accepted by Phys.Rev.