The importance of different timings of excitatory and inhibitory pathways in neural field models

Bumps and oscillons in networks of spiking neurons

We study localized patterns in an exact mean-field description of a spatially extended network of quadratic integrate-and-fire neurons. We investigate conditions for the existence and stability of localized solutions, so-called bumps, and give an analytic estimate for the parameter range, where these solutions exist in parameter space, when one or more microscopic network parameters are varied. We develop Galerkin methods for the model equations, which enable numerical bifurcation analysis of stationary and time-periodic spatially extended solutions. We study the emergence of patterns composed of multiple bumps, which are arranged in a snake-and-ladder bifurcation structure if a homogeneous or heterogeneous synaptic kernel is suitably chosen. Furthermore, we examine time-periodic, spatially localized solutions (oscillons) in the presence of external forcing, and in autonomous, recurrently coupled excitatory and inhibitory networks. In both cases, we observe period-doubling cascades leading to chaotic oscillations.

Kernel Reconstruction for Delayed Neural Field Equations

Understanding the neural field activity for realistic living systems is a challenging task in contemporary neuroscience. Neural fields have been studied and developed theoretically and numerically with considerable success over the past four decades. However, to make effective use of such models, we need to identify their constituents in practical systems. This includes the determination of model parameters and in particular the reconstruction of the underlying effective connectivity in biological tissues.In this work, we provide an integral equation approach to the reconstruction of the neural connectivity in the case where the neural activity is governed by a delay neural field equation. As preparation, we study the solution of the direct problem based on the Banach fixed-point theorem. Then we reformulate the inverse problem into a family of integral equations of the first kind. This equation will be vector valued when several neural activity trajectories are taken as input for the inverse problem. We employ spectral regularization techniques for its stable solution. A sensitivity analysis of the regularized kernel reconstruction with respect to the input signal u is carried out, investigating the Fréchet differentiability of the kernel with respect to the signal. Finally, we use numerical examples to show the feasibility of the approach for kernel reconstruction, including numerical sensitivity tests, which show that the integral equation approach is a very stable and promising approach for practical computational neuroscience.

Traveling waves and breathers in an excitatory-inhibitory neural field

We study existence and stability of traveling activity bump solutions in an excitatory-inhibitory (E-I) neural field with Heaviside firing rate functions by deriving existence conditions for traveling bumps and an Evans function to analyze their spectral stability. Subsequently, we show that these existence and stability results reduce, in the limit of wave speed c→0, to the equivalent conditions developed for the stationary bump case. Using the results for the stationary bump case, we show that drift bifurcations of stationary bumps serve as a mechanism for generating traveling bump solutions in the E-I neural field as parameters are varied. Furthermore, we explore the interrelations between stationary and traveling types of bumps and breathers (time-periodic oscillatory bumps) by bridging together analytical and simulation results for stationary and traveling bumps and their bifurcations in a region of parameter space. Interestingly, we find evidence for a codimension-2 drift-Hopf bifurcation occurring as two parameters, inhibitory time constant τ and I-to-I synaptic connection strength w[over ¯]_{ii}, are varied and show that the codimension-2 point serves as an organizing center for the dynamics of these four types of spatially localized solutions. Additionally, we describe a case involving subcritical bifurcations that lead to traveling waves and breathers as τ is varied.

Emergence of metastable state dynamics in interconnected cortical networks with propagation delays

The importance of the large number of thin-diameter and unmyelinated axons that connect different cortical areas is unknown. The pronounced propagation delays in these axons may prevent synchronization of cortical networks and therefore hinder efficient information integration and processing. Yet, such global information integration across cortical areas is vital for higher cognitive function. We hypothesized that delays in communication between cortical areas can disrupt synchronization and therefore enhance the set of activity trajectories and computations interconnected networks can perform. To evaluate this hypothesis, we studied the effect of long-range cortical projections with propagation delays in interconnected large-scale cortical networks that exhibited spontaneous rhythmic activity. Long-range connections with delays caused the emergence of metastable, spatio-temporally distinct activity states between which the networks spontaneously transitioned. Interestingly, the observed activity patterns correspond to macroscopic network dynamics such as globally synchronized activity, propagating wave fronts, and spiral waves that have been previously observed in neurophysiological recordings from humans and animal models. Transient perturbations with simulated transcranial alternating current stimulation (tACS) confirmed the multistability of the interconnected networks by switching the networks between these metastable states. Our model thus proposes that slower long-range connections enrich the landscape of activity states and represent a parsimonious mechanism for the emergence of multistability in cortical networks. These results further provide a mechanistic link between the known deficits in connectivity and cortical state dynamics in neuropsychiatric illnesses such as schizophrenia and autism, as well as suggest non-invasive brain stimulation as an effective treatment for these illnesses.

Distributed nonlocal feedback delays may destabilize fronts in neural fields, distributed transmission delays do not

The spread of activity in neural populations is a well-known phenomenon. To understand the propagation speed and the stability of stationary fronts in neural populations, the present work considers a neural field model that involves intracortical and cortico-cortical synaptic interactions. This includes distributions of axonal transmission speeds and nonlocal feedback delays as well as general classes of synaptic interactions. The work proves the spectral stability of standing and traveling fronts subject to general transmission speeds for large classes of spatial interactions and derives conditions for the front instabilities subjected to nonlocal feedback delays. Moreover, it turns out that the uniqueness of the stationary traveling fronts guarantees its exponential stability for vanishing feedback delay. Numerical simulations complement the analytical findings.

Transmission efficiency in ring, brain inspired neuronal networks. Information and energetic aspects

Organisms often evolve as compromises, and many of these compromises can be expressed in terms of energy efficiency. Thus, many authors analyze energetic costs processes during information transmission in the brain. In this paper we study information transmission rate per energy used in a class of ring, brain inspired neural networks, which we assume to involve components like excitatory and inhibitory neurons or long-range connections. Choosing model of neuron we followed a probabilistic approach proposed by Levy and Baxter (2002), which contains all essential qualitative mechanisms participating in the transmission process and provides results consistent with physiologically observed values. Our research shows that all network components, in broad range of conditions, significantly improve the information-energetic efficiency. It turned out that inhibitory neurons can improve the information-energetic transmission efficiency by 50%, while long-range connections can improve the efficiency even by 70%. We also found that the most effective is the network with the smallest size: we observed that two times increase of the size can cause even three times decrease of the information-energetic efficiency. This article is part of a Special Issue entitled Neural Coding 2012.

Reduced dynamics for delayed systems with harmonic or stochastic forcing

The analysis of nonlinear delay-differential equations (DDEs) subjected to external forcing is difficult due to the infinite dimensionality of the space in which they evolve. To simplify the analysis of such systems, the present work develops a non-homogeneous center manifold (CM) reduction scheme, which allows the derivation of a time-dependent order parameter equation in finite dimension. This differential equation captures the major dynamical features of the delayed system. The forcing is assumed to be small compared to the amplitude of the autonomous system, in order to cause only small variations of the fixed points and of the autonomous CM. The time-dependent CM is shown to satisfy a non-homogeneous partial differential equation. We first briefly review CM theory for DDEs. Then we show, for the general scalar case, how an ansatz that separates the CM into one for the autonomous problem plus an additional time-dependent order-two correction leads to satisfying results. The paper then details the application to a transcritical bifurcation subjected to single or multiple periodic forcings. The validity limits of the reduction scheme are also highlighted. Finally, we characterize the specific case of additive stochastic driving of the transcritical bifurcation, where additive white noise shifts the mode of the probability density function of the state variable to larger amplitudes.

Neural population densities shape network correlations

The way sensory microcircuits manage cellular response correlations is a crucial question in understanding how such systems integrate external stimuli and encode information. Most sensory systems exhibit heterogeneities in terms of population sizes and features, which all impact their dynamics. This work addresses how correlations between the dynamics of neural ensembles depend on the relative size or density of excitatory and inhibitory populations. To do so, we study an apparently symmetric system of coupled stochastic differential equations that model the evolution of the populations' activities. Excitatory and inhibitory populations are connected by reciprocal recurrent connections, and both receive different stimuli exhibiting a certain level of correlation with each other. A stability analysis is performed, which reveals an intrinsic asymmetry in the distribution of the fixed points with respect to the threshold of the nonlinearities. Based on this, we show how the cross correlation between the population responses depends on the density of the inhibitory population, and that a specific ratio between both population sizes leads to a state of zero correlation. We show that this so-called asynchronous state subsists, despite the presence of stimulus correlation, and most importantly, that it occurs only in asymmetrical systems where one population outnumbers the other. Using linear approximations, we derive analytical expressions for the root of the cross-correlation function and study how the asynchronous state is impacted by the model's parameters. This work suggests a possible explanation for why inhibitory cells outnumber excitatory cells in the visual system.

ON THE ORIGIN AND PROPERTIES OF TWO-POPULATION NEURAL FIELD MODELS - A TUTORIAL INTRODUCTION

Neural field models have a long tradition in mathematical neuroscience, and in the present tutorial paper we outline the neurobiological and biophysical origin of such models, in particular two-population field models describing excitatory and inhibitory neurons interacting via nonlocal spatial connections. Results from investigations of such models on the existence and stability of stationary localized activity pulses ('bumps') and generation of stationary spatial and spatiotemporal oscillations through Turing-type instabilities are described.

Neural adaptation facilitates oscillatory responses to static inputs in a recurrent network of ON and OFF cells

We investigate the role of adaptation in a neural field model, composed of ON and OFF cells, with delayed all-to-all recurrent connections. As external spatially profiled inputs drive the network, ON cells receive inputs directly, while OFF cells receive an inverted image of the original signals. Via global and delayed inhibitory connections, these signals can cause the system to enter states of sustained oscillatory activity. We perform a bifurcation analysis of our model to elucidate how neural adaptation influences the ability of the network to exhibit oscillatory activity. We show that slow adaptation encourages input-induced rhythmic states by decreasing the Andronov-Hopf bifurcation threshold. We further determine how the feedback and adaptation together shape the resonant properties of the ON and OFF cell network and how this affects the response to time-periodic input. By introducing an additional frequency in the system, adaptation alters the resonance frequency by shifting the peaks where the response is maximal. We support these results with numerical experiments of the neural field model. Although developed in the context of the circuitry of the electric sense, these results are applicable to any network of spontaneously firing cells with global inhibitory feedback to themselves, in which a fraction of these cells receive external input directly, while the remaining ones receive an inverted version of this input via feedforward di-synaptic inhibition. Thus the results are relevant beyond the many sensory systems where ON and OFF cells are usually identified, and provide the backbone for understanding dynamical network effects of lateral connections and various forms of ON/OFF responses.

Bifurcations of lurching waves in a thalamic neuronal network

We consider a two-layer, one-dimensional lattice of neurons; one layer consists of excitatory thalamocortical neurons, while the other is comprised of inhibitory reticular thalamic neurons. Such networks are known to support "lurching" waves, for which propagation does not appear smooth, but rather progresses in a saltatory fashion; these waves can be characterized by different spatial widths (different numbers of neurons active at the same time). We show that these lurching waves are fixed points of appropriately defined Poincaré maps, and follow these fixed points as parameters are varied. In this way, we are able to explain observed transitions in behavior, and, in particular, to show how branches with different spatial widths are linked with each other. Our computer-assisted analysis is quite general and could be applied to other spatially extended systems which exhibit this non-trivial form of wave propagation.

Responses of recurrent nets of asymmetric ON and OFF cells

A neural field model of ON and OFF cells with all-to-all inhibitory feedback is investigated. External spatiotemporal stimuli drive the ON and OFF cells with, respectively, direct and inverted polarity. The dynamic differences between networks built of ON and OFF cells ("ON/OFF") and those having only ON cells ("ON/ON") are described for the general case where ON and OFF cells can have different spontaneous firing rates; this asymmetric case is generic. Neural responses to nonhomogeneous static and time-periodic inputs are analyzed in regimes close to and away from self-oscillation. Static stimuli can cause oscillatory behavior for certain asymmetry levels. Time-periodic stimuli expose dynamical differences between ON/OFF and ON/ON nets. Outside the stimulated region, we show that ON/OFF nets exhibit frequency doubling, while ON/ON nets cannot. On the other hand, ON/ON networks show antiphase responses between stimulated and unstimulated regions, an effect that does not rely on specific receptive field circuitry. An analysis of the resonance properties of both net types reveals that ON/OFF nets exhibit larger response amplitude. Numerical simulations of the neural field models agree with theoretical predictions for localized static and time-periodic forcing. This is also the case for simulations of a network of noisy integrate-and-fire neurons. We finally discuss the application of the model to the electrosensory system and to frequency-doubling effects in retina.

Oscillatory response in a sensory network of ON and OFF cells with instantaneous and delayed recurrent connections

A neural field model with multiple cell-to-cell feedback connections is investigated. Our model incorporates populations of ON and OFF cells, receiving sensory inputs with direct and inverted polarity, respectively. Oscillatory responses to spatially localized stimuli are found to occur via Andronov-Hopf bifurcations of stationary activity. We explore the impact of multiple delayed feedback components as well as additional excitatory and/or inhibitory non-delayed recurrent signals on the instability threshold. Paradoxically, instantaneous excitatory recurrent terms are found to enhance network responsiveness by reducing the oscillatory response threshold, allowing smaller inputs to trigger oscillatory activity. Instantaneous inhibitory components do the opposite. The frequency of these response oscillations is further shaped by the polarity of the non-delayed terms.

Large-scale neural dynamics: simple and complex

We review the use of neural field models for modelling the brain at the large scales necessary for interpreting EEG, fMRI, MEG and optical imaging data. Albeit a framework that is limited to coarse-grained or mean-field activity, neural field models provide a framework for unifying data from different imaging modalities. Starting with a description of neural mass models, we build to spatially extend cortical models of layered two-dimensional sheets with long range axonal connections mediating synaptic interactions. Reformulations of the fundamental non-local mathematical model in terms of more familiar local differential (brain wave) equations are described. Techniques for the analysis of such models, including how to determine the onset of spatio-temporal pattern forming instabilities, are reviewed. Extensions of the basic formalism to treat refractoriness, adaptive feedback and inhomogeneous connectivity are described along with open challenges for the development of multi-scale models that can integrate macroscopic models at large spatial scales with models at the microscopic scale.

Dynamics of driven recurrent networks of ON and OFF cells

A globally coupled network of ON and OFF cells is studied using neural field theory. ON cells increase their activity when the amplitude of an external stimulus increases, while OFF cells do the opposite given the same stimulus. Theory predicts that, without input, multiple transitions to oscillations can occur depending on feedback delay and the difference between ON and OFF resting states. Static spatial stimuli can induce or suppress global oscillations via a Andronov-Hopf bifurcation. This is the case for either polarity of such stimuli. In contrast, only excitatory inputs can induce or suppress oscillations in an equivalent network built of ON cells only even though oscillations are more prevalent in such systems. Nonmonotonic responses to local stimuli occur where responses lateral to the stimulus switch from excitatory to inhibitory as the input amplitude increases. With local time-periodic forcing, the unforced cells oscillate at twice the driving frequency via full-wave rectification mediated by the feedback. Our results agree with simulations of the neural field model, and further, qualitative agreement is found with the behavior of a network of spiking stochastic integrate-and-fire model neurons.

Delays in activity-based neural networks

In this paper, we study the effect of two distinct discrete delays on the dynamics of a Wilson-Cowan neural network. This activity-based model describes the dynamics of synaptically interacting excitatory and inhibitory neuronal populations. We discuss the interpretation of the delays in the language of neurobiology and show how they can contribute to the generation of network rhythms. First, we focus on the use of linear stability theory to show how to destabilize a fixed point, leading to the onset of oscillatory behaviour. Next, we show for the choice of a Heaviside nonlinearity for the firing rate that such emergent oscillations can be either synchronous or anti-synchronous, depending on whether inhibition or excitation dominates the network architecture. To probe the behaviour of smooth (sigmoidal) nonlinear firing rates, we use a mixture of numerical bifurcation analysis and direct simulations, and uncover parameter windows that support chaotic behaviour. Finally, we comment on the role of delays in the generation of bursting oscillations, and discuss natural extensions of the work in this paper.

Modeling electrocortical activity through improved local approximations of integral neural field equations

Neural field models of firing rate activity typically take the form of integral equations with space-dependent axonal delays. Under natural assumptions on the synaptic connectivity we show how one can derive an equivalent partial differential equation (PDE) model that properly treats the axonal delay terms of the integral formulation. Our analysis avoids the so-called long-wavelength approximation that has previously been used to formulate PDE models for neural activity in two spatial dimensions. Direct numerical simulations of this PDE model show instabilities of the homogeneous steady state that are in full agreement with a Turing instability analysis of the original integral model. We discuss the benefits of such a local model and its usefulness in modeling electrocortical activity. In particular, we are able to treat "patchy" connections, whereby a homogeneous and isotropic system is modulated in a spatially periodic fashion. In this case the emergence of a "lattice-directed" traveling wave predicted by a linear instability analysis is confirmed by the numerical simulation of an appropriate set of coupled PDEs.

Generalization of the reaction-diffusion, Swift-Hohenberg, and Kuramoto-Sivashinsky equations and effects of finite propagation speeds

The work proposes and studies a model for one-dimensional spatially extended systems, which involve nonlocal interactions and finite propagation speed. It shows that the general reaction-diffusion equation, the Swift-Hohenberg equation, and the general Kuramoto-Sivashinsky equation represent special cases of the proposed model for limited spatial interaction ranges and for infinite propagation speeds. Moreover, the Swift-Hohenberg equation is derived from a general energy functional. After a detailed validity study on the generalization conditions, the three equations are extended to involve finite propagation speeds. Moreover, linear stability studies of the extended equations reveal critical propagation speeds and unusual types of instabilities in all three equations. In addition, an extended diffusion equation is derived and studied in some detail with respect to finite propagation speeds. The extended model allows for the explanation of recent experimental results on non-Fourier heat conduction in nonhomogeneous material.