- Research
- Open access
- Published:
Distributed Nonlocal Feedback Delays May Destabilize Fronts in Neural Fields, Distributed Transmission Delays Do Not
The Journal of Mathematical Neuroscience volume 3, Article number: 9 (2013)
Abstract
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.
1 Introduction
The spatio-temporal dynamics of extended neuronal networks has attracted much attention in recent years [3, 9, 48, 49]. They are powerful models to reproduce encephalographic data [40], to explain phenomena observed in medicine [39] such as general anaesthesia [26, 29, 34, 44] and describe experimental spatio-temporal propagation of electric activity in neural tissues [23, 38, 41].
The spatially-extended neural network under study implies axonal connections with finite transmission speeds, which essentially leads to transmission delays. This delay depends strongly on the axonal branching architecture [2] and the degree of myelination of axonal branches [42]. For instance, unmyelinated axons exhibit a small transmission speed in the range of 0.1–1.0 m/s and occur mainly in short-range intracortical connections. In long-range axonal fibers such as cortico-cortical connections, the axons are myelinated leading to a faster transmission speed in the range of 1 m/s–100 m/s. Consequently, the resulting transmission delay between two spatial locations depend on the axonal paths and varies between 0.5 ms and 100 ms. Since these delay times are in the same range as time constants of synaptic responses of tens of milliseconds, effects of finite transmission delays on the spatio-temporal evolution of activity occur [6, 30, 46]. Although one may estimate such mean transmission delay times along axonal fibers, physiological studies point out that the transmission speed depends on the specific path the action potential takes, and hence varies in a single axonal branching structure from one neuron to another [42]. In addition, the detailed branching structure of neural tissue changes on a time scale of few months [45] or even few days [21], and hence changes the transmission speed. Consequently, it is not reasonable to study the effect of a single transmission speed, but rather a distribution of speeds. The present work considers such a distribution of transmission speeds and extends previous studies assuming a single speed only [13, 30, 37, 38, 46].
The model under study considers two types of axonal pathways. In intra-cortical connections, the lengths of axonal paths may vary due to the absence of fiber bundles of fixed length yielding a transmission delay proportional to the distance between two spatial locations. In contrast, cortico-cortical feedback connections may exhibit fiber bundles with fixed length yielding a constant feedback delay. By virtue of the distribution of transmission speeds, these two pathways are modeled by a distribution of transmission speeds and distributed feedback delays.
Most recent studies of extended neural networks considered either transmission speeds in intraarea connections or delay in feedback connections, although experimental findings indicate the presence of both connections [8], a distribution of transmission speeds [20] and feedback delays [8]. Only few previous studies consider both a single transmission delay in intraarea connections and the feedback delay [24, 28, 36, 51]. The present work extends these studies by a detailed spectral stability analysis.
It is well known that transmission delays and feedback delays may destabilize stationary activity yielding oscillatory phenomena [4, 7, 30]. To better understand the spatio-temporal dynamics of neural populations, it is essential to comprehend the role of delays and their effects on the activity propagation. The present work undertakes the analysis of stationary front for the general case of large classes of delay distributions and spatial interactions and, therefore, aims to reveal answers for this problem.
2 Methods
This section introduces the model equation and gives the major previous results, which represent the basis of the novel results given in Sect. 3. To not expand the section too much, the discussion of the already published findings is kept short and concise.
2.1 The Model Equation
The model under study describes the activity in neural populations on a mesoscopic spatial and temporal scale with typical spatial range of 500 μm and temporal time constants of 5–10 ms. The corresponding spatial domain is coarse-grained and exhibits spatial patches. This structure reflects the macrocolumnar structure observed in primary sensory areas [19]. Such a neural field model allows the successful reproduction of electroencephalographic activity on the head [35] and the successful description of spiral waves in neural tissue [23]. More precisely, the neural field model describes a rate-coding neural population involving synaptic interactions whose spatio-temporal evolution obeys the following nonlinear scalar integro-differential equation:
where represents the mean membrane potential of a spatial patch at position x and time t [4] and is the initial activity. This model neglects units since they can be considered by appropriate scaling of time and space [27]. The prefactors and are nonnegative constants and reflect the synaptic weights of intraarea connections and feedback connections, respectively, and . The function represents the probability density distribution of axonal transmission speeds. Similarly, the function represents the probability density function for feedback delays. We would like to point out that Eq. (1) describes the evolution of a minimal scalar model which, however, considers a large number of features. In spite of its reduced form, i.e., scalar and activity evolves in a one-dimensional space, and it promises to give insights into the effect of various delay types.
To learn more about the wave speed, at some point in the work, we will investigate the dependence of the wave speed of the traveling wave front on the transmission speeds and delays. To this end, the authors assume the axonal transmission speed distribution and the feedback delay distribution to a sum of two terms
i.e., the presence of two axonal transmission speeds and two feedback delays. This choice reflects the presence of short intracortical connections showing a small transmission speed [31] and small constant feedback delay and long-range cortico-cortical connections, which exhibit large axonal transmission speeds and large feedback delays.
In addition, denotes the transfer function of the model and is chosen to the Heaviside step function: for all , , and for all . This assumption is valid for identical neurons in the population [27]. Although this is a strong approximation, it gives first insights into the possible dynamics of homogeneous neural populations. The parameter θ is constant and represents the mean firing threshold of the neurons.
The spatial kernel functions K and W are real-valued and reflect two different nonlocal axonal connectivities between neurons and synapses: K represents intracortical axonal interaction in the neural population, which exhibits finite transmission speeds, whereas W denotes the axonal connectivity along fibers that leave the neural population and reenter it with a constant delay. These functions may be seen as weighted sums of probability density functions of axonal connections of subnetworks where the weights represent the synaptic strengths in each subnetwork [27]. We require that at least in a small open interval , where is a positive constant. This assumption is reasonable since the resulting traveling wave front propagates with a wave speed , that is equal to or smaller than the transmission speed due to physical reasoning. It is assumed that
for two positive constants and .
To investigate various superpositions of subnetworks, the paper considers the following three general classes of synaptic interactions:
-
(A)
This class consists of all nonnegative kernel functions, reflecting global excitation in the neuronal population, i.e., on ℝ.
-
(B)
This class consists of all Mexican hat kernel functions, reflecting lateral inhibition and local excitation in the neuronal population, where each kernel function satisfies the conditions on and on , for a positive constant . This neural interaction is a successful model for interactions in visual receptive fields, e.g., to explain orientation tuning [5, 47].
-
(C)
This class consists of all upside down Mexican hat kernel functions, reflecting lateral excitation and local inhibition, where each kernel function satisfies the conditions on and on , for a positive constant . This interaction is motivated by the physiological finding [32] that networks of inhibitory neurons act locally only and excitatory neurons rarely exhibit local connections, but rather long-range interactions.
The synaptic feedback coupling W satisfies the same assumptions as K and it belongs to one of the three classes. Certainly, K and W are not necessarily in the same class. Moreover, K and W are not necessarily symmetric functions.
2.2 The Standing Wave Front
In the presence of a single firing threshold θ as in Eq. (1), a traveling wave front connects two constant states, i.e., one above and one below the threshold θ. Let and choose in Eq. (1), then
with two exponentially stable stationary states and . They are stable irrespective of the spatial interactions and the delay distributions.
Suppose that and . Now let us consider a stationary standing wave front with profile . Without loss of generality, suppose that the front crosses the threshold θ at the point , that is, on , , and on . Moreover, suppose that , and . Plugging such a solution back into the equation, noting that the solution is independent of time, we get
Amari had studied solutions similar to (5) for general stationary states in his celebrated work [3]. Section 3.1 will elaborate in some detail the spectral stability of the stationary standing wave front given by (5).
It has been shown in many previous studies [4, 27] that certain spatial interactions may destabilize spatially constant stationary states subject to the nonlinear gain function . In the present model, the functional derivative of the nonlinear gain function is given by . In other words, spatial interactions at a certain spatial location x with activity contribute to the stability of a stationary state only if the location is close to the threshold θ.
2.3 The Traveling Wave Front and Its Stability
After the study of the existence of a standing wave front, this subsection focuses on the existence of a travelling wave front. Let , , be constants such that . As shown in the previous section, there are two constant solutions and to Eq. (1).
2.3.1 The Front Shape
Suppose that is a traveling wave front of Eq. (1), where represents the wave speed and represents a moving coordinate. Due to translation invariance, suppose that the traveling wave front satisfies the conditions on , , and on . Suppose that the front satisfies the boundary conditions , and . For physical reasons, suppose that the wave speed satisfies the conditions , where is the smallest occurring transmission speed. Then the traveling wave front and the wave speed μ satisfy the equation
After a series of change of variables (such as and , etc.), this above equation becomes
where is the sign function, which is defined by for all , and for all . By using the integrating factor idea and integration by parts, we find the representation of the front
The derivative of is given by
These expressions are useful in a later part of the work.
2.3.2 The Wave Speed of the Traveling Wave Front
To compute the wave speed, setting and in (8) and making some simple change of variables yields
This equation may be rewritten as
where and are called speed index functions, defined on , by
Define the sub speed index functions and on by
There holds . We will see that there exists a unique positive solution (that is, the wave speed ) to the equation .
Interestingly, the speed index functions and allow us to express the slope of the front at the threshold θ by utilizing (9) and (10)
These expressions will be useful in later discussions of the spectral stability of the traveling wave front.
Moreover, for single transmission speed and single feedback delay (i.e., and , where and are positive constants) and identical spatial kernel functions , we have
Thus, there is no distinction between intracortical and feedback interactions. This is obvious from Eq. (1) where the two integrals may be written as a single one.
2.3.3 The Spectral Stability of the Traveling Wave Front
Real neural structures exhibit a certain level of background noise, which may disturb the propagation of activity. Hence, it is important to study the stability of the stationary traveling wave front with respect to small perturbations. This kind of analysis has already been performed before for similar equations in [4, 14, 16, 43].
Previous studies on integral and/or partial differential equations involving infinite delays have found criteria for the existence, uniqueness, and stability of traveling wave solutions [1, 18]. Following the successful extrapolation approach for infinite delays in integral and/or partial differential equations based on finite delays [12, 18] and recalling the close relationship of partial differential equations and integro-differential equations of the type discussed here [25], we assume in the following that mathematical analysis based on finite delays are applicable. An additional discussion of this approximation may appear necessary, but is neglected since it would exceed the major aim of the present work. Motivated by a previous study of Coombes and Owen [10], we study the spectral stability of the traveling wave front of (1) by constructing Evans functions which are well known from the literature of integral and/or partial differential equations [43].
To study the spectral stability of the front, we must rewrite the equation under consideration in moving coordinate and linearize the new equation with respect to the traveling wave front to obtain a linear equation. Then we seek solutions of the form to separate time and the moving coordinate z and to obtain an eigenvalue problem with the complex eigenvalue λ and the eigenfunction ψ; see also the Appendix. Solving the eigenvalue problem yields the definition of the Evans function. It turns out that a complex number is an eigenvalue of the eigenvalue problem if and only if is a zero of the Evans function. Note that is always an eigenvalue, reflecting the translation invariance of the travelling wave front. Previous studies have shown that the equation determines the isolated spectrum . As before, the stationary front is unstable if there exists some eigenvalue with positive real part or if the neutral eigenvalue is not simple.
The spectrum of the linear differential operator consists of two parts: the essential spectrum and the isolated spectrum (point spectrum comprising the eigenvalues). Please see the Appendix for the definitions of the linear differential operators and . Are there any other spectrum in other than the essential spectrum and the isolated spectrum to the linear differential operator ? By assumption, the kernel functions K and W converge to zero exponentially fast as . Therefore, the operator is a compact operator in . The residual spectrum does not exist in our model.
The essential spectrum is easy to calculate by following the original ideas of John Evans [14]. The complex number belongs to the essential spectrum if and only if λ is a complex number such that the solutions of the differential equation
is bounded on ℝ. The solution of this equation is given by . This solution is bounded on ℝ if and only if , for some real number r. Since the subsequent sections show that linear deviations from stationary traveling wave front obey differential equations of the type (17), we find that the essential spectrum contains those values of λ with only, i.e., the essential spectrum does not threaten the stability of the traveling wave front.
It remains to investigate whether the isolated spectrum (the eigenvalues) threatens the stability, i.e., whether there are eigenvalues with positive real part , threatening the stability of the stationary front.
Define the Evans function for the traveling wave front of Eq. (1) by
where and are also called Evans functions, defined by
in the open domain . Here, the Evans function is also called the stability index function and represents a sum of two single stability index functions and , reflecting the contribution of transmission delay and constant delay, respectively.
Properties of the Evans functions: A complex number is an eigenvalue of the eigenvalue problem (details to be given in the Appendix) if and only if is a zero of the Evans function. Moreover,
2.4 The Numerical Simulations
To complement the analytical study, subsequent sections show numerical integrations of the model equation (1). To this end, if not stated otherwise, the spatial kernel functions are chosen to
with the real parameters , , . Equation (1) is integrated numerically in time by an Euler-forward method with time step , the spatial integral has been computed on a circular spatial grid of length and 600 intervals by the Monte Carlo-type VEGAS-algorithm [22] with 5000 random draws. If not stated differently, the initial values of the neural activity have been chosen to the analytical stationary front perturbed by random values taken from a uniform distribution in the initial interval .
3 Results
This section presents the new findings on the effect of distributed delays in both standing wave front and traveling wave front. They extend previous results mentioned in the previous section.
3.1 The Standing Wave Front
Derivation of an eigenvalue problem. Recall that the standing wave front is a time independent solution of the nonlinear scalar integro-differential equation
To study the spectral stability of the stationary standing wave front, we linearize this equation about the front and we find that obeys the differential equation (neglecting higher order terms)
where , because . Making a change of variable and using (7), we find the equation
where . We understand that the contribution of the spatial interactions to the activity is maximum if and are large. For kernel functions with a maximum value at the origin, the largest contribution is expected at , i.e., close to the threshold θ. Hence, one expects a large change of the activity close to the threshold θ. This is different from other spatial interactions, such as gamma-distributed kernel functions [27], where the strongest contribution occurs away from the threshold.
A standing wave front is translation invariant. We may let . Suppose that is a solution of the above equation, where λ is a complex constant and is a bounded continuous function defined on ℝ. After canceling out , we obtain the eigenvalue problem
The spectral stability of the stationary front (5) is determined by the eigenvalues of this eigenvalue problem. Letting in this equation and canceling out because it is not equal to zero, we get
At first, we observe that the spectral stability of the standing wave front just depends on the spatial self-interactions and , i.e., lateral interactions do not determine the stability of the stationary standing wave front. Moreover, finite transmission speeds do not affect the stability, but feedback delays may do. This result originates from the different nature of the delays: the transmission delay does not affect the dynamics at the threshold θ since it vanishes at zero distance from the threshold θ whereas the feedback delay is independent of the distance and affects the dynamics at all spatial locations.
Spectral analysis. Obviously, the neutral eigenvalue is a trivial solution of this equation. To study the spectral stability of the standing wave front with respect to feedback delays, let us define the auxiliary function on the interval :
For all kernel functions K and W, for all probability density functions ξ and η, f is an increasing of λ. There exists a unique solution to the equation . Moreover, is a simple solution (that is, is a simple eigenvalue). Therefore, the standing wave front is spectrally stable. Let us consider a very special case: absent intraarea connections (i.e., ) and focus on nonlocal feedback connections with a single feedback delay (i.e. , where is a positive constant). Then . That is . Given the positive constant , it is very easy to see that there exists a unique solution to .
Moreover, Fig. 1 shows the space-time activity of oscillatory unstable standing front gained by a numerical simulation of Eq. (1) for a single delay (a) and a distribution of two delays (b). We observe oscillatory activity close to the threshold value θ consistent with the reasoning in Sect. 2.2.
3.2 Traveling Wave Front
This subsection investigates analytically the uniqueness of the stationary travelling wave front and the spectral stability of the front. Numerical simulations validate the analytical findings.
3.2.1 The Uniqueness
Taking a close look at the implicit equation (10), the question arises whether its solution, i.e., the wave speed , is unique, or whether there are several possible traveling wave fronts with different wave speeds. Previous studies on propagating front in neural fields involving a single axonal finite transmission speed [11] have established the uniqueness of the traveling wave front. The present work extends these studies by considering distributed transmission speeds and distributed feedback delays.
If the speed index functions , , and defined in (11), (12), and (10) are monotonic in μ, then the wave speed of the traveling wave front is unique. By using rigorous mathematical analysis (to keep the paper from too long, details are not to be given here), we find that
-
For all kernel functions in classes (A) and (B), for all speed and delay distributions , for all , there hold the following estimates:
(23) -
For all kernel functions in class (C), for all speed and delay distributions , for all , there hold the following estimates:
(24)
for two positive constants and , where and .
In addition,
and hence and . Since , condition (10) stipulates . Consequently, it is necessary that and in at least a small neighborhood of the wave speed .
Overall, there exists a unique wave speed and there is a unique traveling wave front to Eq. (1).
3.2.2 Dependence on Transmission Speed and Feedback Delays
Now we investigate the change of the wave speed while changing the speed and delay distributions assuming the distributions (2). For simplicity, we choose , where κ is a fixed positive constant. Then (23), (24), and (25) and the definitions (11), (12) yield for general synaptic interactions K and W that
and for excitatory delayed feedback interaction, i.e., ,
Figure 2 shows the wave speed given implicitly by (10) with respect to different transmission speed distributions and feedback delay distributions corresponding to (2) for the kernel functions in the three classes (A), (B), and (C).
3.2.3 Stability Analysis
Analysis of the eigenvalues of the eigenvalue problem. The introduction of the speed index functions and the stability index functions given in Sect. 2.3.3 reveals novel relationship between the uniqueness of the wave speed and the stability of the traveling wave front. Comparing the stability index functions given in Eqs. (19), (20) and the speed index functions in Eqs. (11), (12), we can obtain the following results.
-
(I)
For intracortical interactions only (i.e., , , and ), the speed index function from Eq. (11) and the stability index function from (18) are related by
(26)
Let , then
For the eigenvalue , we find that , and hence (27) holds. This eigenvalue reflects the translation invariance of the stationary traveling wave front. To study the spectral stability of the front, the question arises whether there is another complex eigenvalue with positive real part for which (27) holds. Since the kernel function K is real, it follows that the eigenvalue (if it exists) is real and if . Then Eq. (27) breaks down to the question whether there is another solution or, equivalently, whether is monotonic in . This question has already been answered in the context of the uniqueness of the wave speed in Sect. 2.3.2. Since is the unique solution of Eq. (10), Eq. (27) has a unique solution so is the only eigenvalue.
-
(II)
Now let us consider feedback connections only (i.e., , , and ). The stability index function (18) reads
Then (10), (13), (14), and (15) yield
Since instabilities imply , , the left-hand side of (28) is smaller than its right-hand side, and hence (28) holds true only for , i.e., the front is spectrally stable.
-
(III)
For both intracortical and feedback connections (i.e., , , and ), the stability index function reads
Similar to cases (I) and (II), we obtain the relationship
By virtue of the monotonicity of the speed index function , is the only solution of (29), i.e., is the only eigenvalue and the front is spectrally stable.
Now let us consider distributed feedback delays, the stability index function (18) and Eqs. (10), (15) yield
Since for all λ with , the same reasoning as in (I) and (II) applies and the right side of Eq. (30) is smaller than its left side. Hence, the only solution of Eq. (30) is . Moreover, the derivative of the Evans function at is positive, that is . In another words, the neutral eigenvalue is a simple eigenvalue.
On the other hand, note that
and
We write and . Note that both the real part and the imaginary part of the Evans function are real harmonic functions of x and y. Hence, they satisfy the mean value formula; see Evans [14–17]. As a result, also satisfies the mean value formula. Thus, cannot attain a local maximum inside any open domain of ℂ. By using a strong maximum principle of in Ω, we find that , for all with . The spectral stability follows immediately. To illustrate this finding, Fig. 3 shows numerical simulations of Eq. (1) and we observe stationary propagating front for kernel functions in the three classes (A), (B), and (C).
The simulations consider both distributed axonal transmission speeds and distributed feedback delays. We observe that starting from a noise-perturbed stationary front, the activity approaches the smooth stationary propagating front after some time, i.e., the front is also exponentially stable in the presence of feedback delay. This numerical result for nonvanishing feedback connections complements the analytical result for vanishing feedback delay above.
4 Discussion
4.1 Stability of the Standing Wave Front
The results in Sect. 3.1 show that transmission delays do not destabilize standing wave front. In contrast, nonlocal feedback delays may destabilize standing wave front for certain delays resulting to oscillatory activity. The corresponding analytical conditions for a single delay hold for feedback connections in class (C), i.e., for local inhibition–lateral excitation interactions. Since such synaptic couplings as well as feedback delays are omnipresent in the brain, oscillatory standing wave front may occur frequently in real neural structures. For distributed delays, such instability may be possible, but no analytical method is known up to today proving oscillatory instability. Nevertheless, numerical simulations of standing wave front subject to distributed feedback delays reveal oscillatory instability as well.
4.2 Wave Speed of the Traveling Wave Front
Section 3.2 shows that the speed index function is monotonically increasing, i.e., on the interval where the wave speed exist. Consequently, the wave speed is unique. Moreover, the same subsection shows analytically that the increase of the transmission speed, i.e., the decrease of the transmission delay, increases the wave speed. Similarly, the increase of the feedback delays decreases the wave speed. Summarizing these results, increasing delays slows down the wave speed of the traveling wave front.
4.3 Stability of the Travelling Wave Front
The results in Sect. 3.2 indicate that traveling wave front is spectrally stable in the absence of delayed feedback due to the uniqueness of the wave speed. To our best knowledge, this relationship between the uniqueness and the stability has not been found before, although claimed in previous studies [50, 51].
Focusing on the effect of delayed feedback, the stability analysis in Sect. 3.2.3 reveals the stability of traveling wave front for distributed delays. To support the analytical results, numerical simulations were performed showing stable traveling wave front, cf. Fig. 3. Additional extensive numerical studies (not shown) on the effect of distributed feedback delays have found stable traveling wave front only.
Summarizing the latter results, distributed transmission delays do not destabilize traveling wave front, but feedback delays may induce oscillatory instability.
4.4 Summary and Outlook
We find that the standing wave front and the stationary traveling wave front involving distributed transmission speeds exhibit a unique stable traveling wave front to Eq. (1) given any pair of synaptic couplings (K, W), probability density functions (ξ, η), synaptic rate constants (α, β), any threshold θ and if assumptions (3) hold. In addition, the standing wave front and the stationary traveling wave front are spectrally stable in the presence of small external perturbations due to the found relationship between uniqueness and stability. In contrast, the additional presence of nonlocal feedback delays may render the stationary front instability. We find that the delay-induced loss of stability is oscillatory. This is valid for the present model involving a single synaptic time scale. It is worth mentioning, however, that neural fields involving multiple synaptic time scales are sensitive to transmission delays [33].
Future work may further investigate the relation of uniqueness and stability in other wave phenomena, such as traveling pulse solutions or global waves. Moreover, it is interesting to further study the impact of nonlocal delayed feedback on these wave phenomena. These studies may permit deeper insight into the role of delayed feedback loops which are omnipresent in neural systems.
Appendix
Let represent the traveling wave front of the nonlinear scalar integro-differential equation (1).
Definition 1 (Eigenvalue problems, eigenvalues, and eigenfunctions)
Let and . Define a linear differential operator by
Define a family of linear differential operators by using the eigenvalue problem
The eigenvalue problem is given explicitly by
If there exists a complex number and there exists a nontrivial bounded continuously differentiable function defined on ℝ, such that , then is called an eigenvalue and is called an eigenfunction of the eigenvalue problem.
Definition 2 (The spectral stability)
-
(1)
If and if is a simple eigenvalue of the eigenvalue problem , then we say that the traveling wave front is spectrally stable, where represents the spectrum of the operator , a positive constant.
-
(2)
If there exists an eigenvalue with positive real part to the eigenvalue problem , then we say that the traveling wave front is exponentially unstable.
-
(3)
If , that is, there exists no nonzero eigenvalue to the eigenvalue problem in the right half plane , but is not a simple eigenvalue of the eigenvalue problem , then we say that the traveling wave front is algebraically unstable.
Therefore, we have finished the proof of the spectral stability.
References
Adimy M, Bouzahir H, Ezzinbi K: Local existence and stability for some partial functional differential equations with infinite delay. Nonlinear Anal, Theory Methods Appl 2002, 48: 323–348. 10.1016/S0362-546X(00)00184-X
Agmon-Sir H, Segev I: Signal delay and input synchronisation in passive dendritic structures. J Neurophysiol 1993, 70: 2066–2085.
Amari S: Dynamics of pattern formation in lateral-inhibition type neural fields. Biol Cybern 1977, 27: 77–87. 10.1007/BF00337259
Atay FM, Hutt A: Neural fields with distributed transmission speeds and long-range feedback delays. SIAM J Appl Dyn Syst 2006, 5: 670–698. 10.1137/050629367
Ben-Yishai R, Bar-Or RL, Sompolinsky H: Theory of orientation tuning in visual cortex. Proc Natl Acad Sci USA 1995, 92: 3844–3848. 10.1073/pnas.92.9.3844
Bojak I, Liley DTJ: Axonal velocity distributions in neural field equations. PLoS Comput Biol 2010., 6: Article ID e1000653 Article ID e1000653
Bressloff PC, Folias SE: Front bifurcations in an excitatory neural network. SIAM J Appl Math 2004, 65: 131–151. 10.1137/S0036139903434481
Bullier J, Hupe JM, James AC, Girard P: The role of feedback connections in shaping the responses of visual cortical neurons. Prog Brain Res 2001, 134: 193–204.
Coombes S: Waves, bumps, and patterns in neural field theories. Biol Cybern 2005, 93: 91–108. 10.1007/s00422-005-0574-y
Coombes S, Owen MR: Evans functions for integral neural field equations with Heaviside firing rate function. SIAM J Appl Dyn Syst 2004, 3: 574–600. 10.1137/040605953
Coombes S, Lord GJ, Owen MR: Waves and bumps in neuronal networks with axo-dendritic synaptic interactions. Physica D 2003, 178: 219–241. 10.1016/S0167-2789(03)00002-2
Engel K, Nagel R Graduate Texts in Mathematics 194. In One-Parameter Semigroups for Linear Evolution Equations. Springer, Berlin; 2001.
Ermentrout GB, McLeod JB: Existence and uniqueness of travelling waves for a neural network. Proc R Soc Edinb, Sect A, Math 1993, 123(3):461–478. 10.1017/S030821050002583X
Evans JW: Nerve axon equations, I Linear approximations. Indiana Univ Math J 1971, 21: 877–885.
Evans JW: Nerve axon equations, II Stability at rest. Indiana Univ Math J 1972, 22: 75–90. 10.1512/iumj.1973.22.22009
Evans JW: Nerve axon equations, III Stability of the nerve impulse. Indiana Univ Math J 1972, 22: 577–593. 10.1512/iumj.1973.22.22048
Evans JW: Nerve axon equations, IV The stable and the unstable impulse. Indiana Univ Math J 1975, 24: 1169–1190. 10.1512/iumj.1975.24.24096
Ezzinbi K: Existence and stability for some partial functional differential equations with infinite delay. Electron J Differ Equ 2003., 2003: Article ID 116 Article ID 116
Gerstner W, Kistler W: Spiking Neuron Models. Cambridge University Press, Cambridge; 2002.
Girard P, Hupe JM, Bullier J: Feedforward and feedback connections between areas V1 and V2 of the monkey have similar rapid conduction velocities. J Neurophysiol 2001, 85: 1328–1331.
Gogolla N, Galimberti I, Caroni P: Structural plasticity of axon terminals in the adult. Curr Opin Neurobiol 2007, 17(5):516–524. 10.1016/j.conb.2007.09.002
Gough B: GNU Scientific Library Reference Manual. 2nd edition. Network Theory Ltd, Godalming; 2003.
Huang X, Troy WC, Yang Q, Ma H, Laing CR, Schiff SJ, Wu JY: Spiral waves in disinhibited mammalian neocortex. J Neurosci 2004, 24: 9897–9902. 10.1523/JNEUROSCI.2705-04.2004
Hutt A: Effects of nonlocal feedback on travelling fronts in neural fields subject to transmission delay. Phys Rev E 2004., 70: Article ID 052902 Article ID 052902
Hutt A: Generalization of the reaction-diffusion, Swift-Hohenberg, and Kuramoto-Sivashinsky equations and effects of finite propagation speeds. Phys Rev E 2007., 75: Article ID 026214 Article ID 026214
Hutt A (Ed): Sleep and Anesthesia: Neural Correlates in Theory and Experiment. Springer, New York; 2011.
Hutt A, Atay FM: Analysis of nonlocal neural fields for both general and gamma-distributed connectivities. Physica D 2005, 203: 30–54. 10.1016/j.physd.2005.03.002
Hutt A, Atay FM: Effects of distributed transmission speeds on propagating activity in neural populations. Phys Rev E 2006., 73: Article ID 021906 Article ID 021906
Hutt A, Longtin A: Effects of the anesthetic agent propofol on neural populations. Cogn Neurodyn 2009, 4: 37–59.
Hutt A, Rougier N: Activity spread and breathers induced by finite transmission speeds in two-dimensional neural fields. Phys Rev E 2010., 82: Article ID R055701 Article ID R055701
Hutt A, Bestehorn M, Wennekers T: Pattern formation in intracortical neuronal fields. Netw Comput Neural Syst 2003, 14: 351–368. 10.1088/0954-898X/14/2/310
Koch C: Biophysics of Computation: Information Processing in Single Neurons. Oxford University Press, London; 1999.
Laing C, Coombes S: The importance of different timings of excitatory and inhibitory pathways in neural field models. Network 2006, 17: 151–172. 10.1080/09548980500533461
Liley DTJ, Bojak I: Understanding the transition to seizure by modeling the epileptiform activity of general anaesthetic agents. J Clin Neurophysiol 2005, 22: 300–313.
Liley DTJ, Cadusch PJ, Dafilis MP: A spatially continuous mean field theory of electrocortical activity. Netw Comput Neural Syst 2002, 13: 67–113.
Magpantay FMG, Zou X: Wave fronts in neuronal fields with nonlocal post-synaptic axonal connections and delayed nonlocal feedback connections. Math Biosci Eng 2010, 7: 421–442.
Pinto DJ, Ermentrout GB: Spatially structured activity in synaptically coupled neuronal networks. I. Travelling fronts and pulses. SIAM J Appl Math 2001, 62: 206–225. 10.1137/S0036139900346453
Prechtl JC, Cohen LB, Pesaran B, Mitra PP, Kleinfeld D: Visual stimuli induce waves of electrical activity in turtle cortex. Proc Natl Acad Sci USA 1997, 94: 7621–7626. 10.1073/pnas.94.14.7621
Reggia JA, Montgomery D: A computational model of visual hallucinations in migraine. Comput Biol Med 1996, 26: 133–139. 10.1016/0010-4825(95)00051-8
Rennie CJ, Robinson PA, Wright JJ: Unified neurophysical model of EEG spectra and evoked potentials. Biol Cybern 2002, 86: 457–471. 10.1007/s00422-002-0310-9
Richardson KA, Schiff SJ, Gluckman BJ: Control of travelling waves in the mammalian cortex. Phys Rev Lett 2005., 94: Article ID 028103 Article ID 028103
Ringo JL, Doty RW, Dementer S, Simard PY: Time is of the essence: a conjecture that hemispheric specialization arises from interhemispheric conduction delay. Cereb Cortex 1994, 4: 331–343. 10.1093/cercor/4.4.331
Sandstede B: Evans functions and nonlinear stability of travelling waves in neuronal network models. Int J Bifurc Chaos 2007, 17: 2693–2704. 10.1142/S0218127407018695
Steyn-Ross ML, Steyn-Ross DA, Sleigh JW, Wilcocks LC: Toward a theory of the general-anesthetic-induced phase transition of the cerebral cortex: I. A thermodynamic analogy. Phys Rev E 2001., 64: Article ID 011917 Article ID 011917
Swadlow HA: Impulse conduction in the mammalian brain: physiological properties of individual axons monitored for several months. Science 1982, 218: 911–913. 10.1126/science.7134984
Veltz R, Faugeras O: Stability of the stationary solutions of neural field equations with propagation delays. J Math Neurosci 2011, 1: 1–28. 10.1186/2190-8567-1-1
Wennekers T: Orientation tuning properties of simple cells in area V1 derived from an approximate analysis of nonlinear neural field models. Neural Comput 2001, 13: 1721–1747. 10.1162/08997660152469323
Wilson HR, Cowan JD: A mathematical theory of the functional dynamics of cortical and thalamic nervous tissue. Kybernetic 1973, 13: 55–80. 10.1007/BF00288786
Zhang L: How do synaptic coupling and spatial temporal delay influence travelling waves in nonlinear nonlocal neuronal networks? SIAM J Appl Dyn Syst 2007, 6: 597–644. 10.1137/06066789X
Zhang L, Wu P-S, Stoner MA: Influence of sodium currents on speeds of travelling wave fronts in synaptically coupled neuronal networks. Physica D 2010, 239: 9–32. 10.1016/j.physd.2009.09.022
Zhang L, Wu P-S, Stoner MA: Influence of neurobiological mechanisms on speeds of travelling wave fronts in mathematical neuroscience. Discrete Contin Dyn Syst, Ser B 2011, 16: 1003–1037.
Acknowledgements
AH acknowledges funding from the European Research Council for support under the European Union’s Seventh Framework Programme (FP7/2007-2013)/ERC grant agreement no. 257253. LZ acknowledges funding from Faculty Research Grant and Faculty Grant for International Connections of Lehigh University.
Author information
Authors and Affiliations
Corresponding authors
Additional information
Competing Interests
The authors declare that no competing interests exist.
Authors’ Contributions
AH carried out the numerical simulations, AH and LZ designed the study and performed the analytical computations. Both authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License (https://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
About this article
Cite this article
Hutt, A., Zhang, L. Distributed Nonlocal Feedback Delays May Destabilize Fronts in Neural Fields, Distributed Transmission Delays Do Not. J. Math. Neurosc. 3, 9 (2013). https://doi.org/10.1186/2190-8567-3-9
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/2190-8567-3-9