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

Integration of continuous-time dynamics in a spiking neural network simulator

Contemporary modeling approaches to the dynamics of neural networks consider two main classes of models: biologically grounded spiking neurons and functionally inspired rate-based units. The unified simulation framework presented here supports the combination of the two for multi-scale modeling approaches, the quantitative validation of mean-field approaches by spiking network simulations, and an increase in reliability by usage of the same simulation code and the same network model specifications for both model classes. While most efficient spiking simulations rely on the communication of discrete events, rate models require time-continuous interactions between neurons. Exploiting the conceptual similarity to the inclusion of gap junctions in spiking network simulations, we arrive at a reference implementation of instantaneous and delayed interactions between rate-based models in a spiking network simulator. The separation of rate dynamics from the general connection and communication infrastructure ensures flexibility of the framework. We further demonstrate the broad applicability of the framework by considering various examples from the literature ranging from random networks to neural field models. The study provides the prerequisite for interactions between rate-based and spiking models in a joint simulation

Limits and dynamics of stochastic neuronal networks with random heterogeneous delays

Realistic networks display heterogeneous transmission delays. We analyze here the limits of large stochastic multi-populations networks with stochastic coupling and random interconnection delays. We show that depending on the nature of the delays distributions, a quenched or averaged propagation of chaos takes place in these networks, and that the network equations converge towards a delayed McKean-Vlasov equation with distributed delays. Our approach is mostly fitted to neuroscience applications. We instantiate in particular a classical neuronal model, the Wilson and Cowan system, and show that the obtained limit equations have Gaussian solutions whose mean and standard deviation satisfy a closed set of coupled delay differential equations in which the distribution of delays and the noise levels appear as parameters. This allows to uncover precisely the effects of noise, delays and coupling on the dynamics of such heterogeneous networks, in particular their role in the emergence of synchronized oscillations. We show in several examples that not only the averaged delay, but also the dispersion, govern the dynamics of such networks.Comment: Corrected misprint (useless stopping time) in proof of Lemma 1 and clarified a regularity hypothesis (remark 1

Invariant Measures for Dissipative Dynamical Systems: Abstract Results and Applications

In this work we study certain invariant measures that can be associated to the time averaged observation of a broad class of dissipative semigroups via the notion of a generalized Banach limit. Consider an arbitrary complete separable metric space $X$ which is acted on by any continuous semigroup $\{S(t)\}_{t \geq 0}$. Suppose that $\S(t)\}_{t \geq 0}$ possesses a global attractor $\mathcal{A}$. We show that, for any generalized Banach limit $\underset{T \rightarrow \infty}{\rm{LIM}}$ and any distribution of initial conditions $\mathfrak{m}_0$, that there exists an invariant probability measure $\mathfrak{m}$, whose support is contained in $\mathcal{A}$, such that $$\int_{X} \phi(x) d\mathfrak{m} (x) = \underset{T\to \infty}{\rm{LIM}} \frac{1}{T}\int_0^T \int_X \phi(S(t) x) d \mathfrak{m}_0(x) d t,$$ for all observables $\phi$ living in a suitable function space of continuous mappings on $X$. This work is based on a functional analytic framework simplifying and generalizing previous works in this direction. In particular our results rely on the novel use of a general but elementary topological observation, valid in any metric space, which concerns the growth of continuous functions in the neighborhood of compact sets. In the case when $\{S(t)\}_{t \geq 0}$ does not possess a compact absorbing set, this lemma allows us to sidestep the use of weak compactness arguments which require the imposition of cumbersome weak continuity conditions and limits the phase space $X$ to the case of a reflexive Banach space. Two examples of concrete dynamical systems where the semigroup is known to be non-compact are examined in detail.Comment: To appear in Communications in Mathematical Physic

Limits and dynamics of randomly connected neuronal networks

Networks of the brain are composed of a very large number of neurons connected through a random graph and interacting after random delays that both depend on the anatomical distance between cells. In order to comprehend the role of these random architectures on the dynamics of such networks, we analyze the mesoscopic and macroscopic limits of networks with random correlated connectivity weights and delays. We address both averaged and quenched limits, and show propagation of chaos and convergence to a complex integral McKean-Vlasov equations with distributed delays. We then instantiate a completely solvable model illustrating the role of such random architectures in the emerging macroscopic activity. We particularly focus on the role of connectivity levels in the emergence of periodic solutions

Stability of stochastic impulsive differential equations: integrating the cyber and the physical of stochastic systems

According to Newton's second law of motion, we humans describe a dynamical system with a differential equation, which is naturally discretized into a difference equation whenever a computer is used. The differential equation is the physical model in human brains and the difference equation the cyber model in computers for the dynamical system. The physical model refers to the dynamical system itself (particularly, a human-designed system) in the physical world and the cyber model symbolises it in the cyber counterpart. This paper formulates a hybrid model with impulsive differential equations for the dynamical system, which integrates its physical model in real world/human brains and its cyber counterpart in computers. The presented results establish a theoretic foundation for the scientific study of control and communication in the animal/human and the machine (Norbert Wiener) in the era of rise of the machines as well as a systems science for cyber-physical systems (CPS)

Linear feedback stabilization of a dispersively monitored qubit

The state of a continuously monitored qubit evolves stochastically, exhibiting competition between coherent Hamiltonian dynamics and diffusive partial collapse dynamics that follow the measurement record. We couple these distinct types of dynamics together by linearly feeding the collected record for dispersive energy measurements directly back into a coherent Rabi drive amplitude. Such feedback turns the competition cooperative, and effectively stabilizes the qubit state near a target state. We derive the conditions for obtaining such dispersive state stabilization and verify the stabilization conditions numerically. We include common experimental nonidealities, such as energy decay, environmental dephasing, detector efficiency, and feedback delay, and show that the feedback delay has the most significant negative effect on the feedback protocol. Setting the measurement collapse timescale to be long compared to the feedback delay yields the best stabilization.Comment: 16 pages, 7 figure