Skip to content

Advertisement

  • Poster presentation
  • Open Access

Numerical simulations in two-dimensional neural fields

BMC Neuroscience201516 (Suppl 1) :P22

https://doi.org/10.1186/1471-2202-16-S1-P22

  • Published:

Keywords

  • Membrane Potential
  • Positive Constant
  • Kernel Function
  • Firing Rate
  • Field Equation
We present and discuss a new numerical algorithm for solving the two-dimensional Neural Field Equation with space-dependent delays:
c t V ( x ¯ , t ) = I ( x ¯ , t ) V ( x ¯ , t ) + Ω K ( | x ¯ y ¯ | ) S ( V ( y , t τ ( x ¯ , y ¯ ) ) d y ¯ t [ 0 , T ] , x ¯ Ω 2

used in Neuroscience to describe the evolution of a population of neurons and the interactions between them. A similar equation (without delays) was first introduced by Wilson and Cowan [1], and then by Amari [2]. Here V( x ,t) represents the post-synaptic membrane potential at instant t and position x. The function I represents external sources of excitation and S describes the dependence between the firing rate of the neurons and their membrane potential (typically it is a function of sigmoidal type). The kernel function K(| x-y |) describes the connectivity between the neurons at positions × and y. The delay τ(x,y) takes into the consideration the time spent by an electrical signal to travel between these two positions.

The new numerical method presented in this work is based on an implicit second order method for discretisation in time and uses Gaussian quadratures for space integration. We use low-rank methods to reduce the computational effort, which enables to reduce very significantly the dimensions of the involved matrices, without affecting the final accuracy of the method.

This algorithm is targeted directly to the application in Neuroscience and Robotics. As an illustration of such applications, we consider a neural field described in [3], where the firing rate function has the form S(x)= 2/ (1+exp(- μ x)) (μ is a certain positive constant), and the connectivity function is given by K(x)= 1/(2 πξ12)1/2 exp(- ||x||2/ 2 ξ12) - A/(2 πξ22)1/2 exp(- ||x||2/ 2 ξ22), with A, ξ1, ξ2 given positive constants.

It is known [3] that the stability of the trivial solution of this equation depends on the value of μ and on the parameters of the connectivity function. In particular, for each set of values of these parameters, it is possible to compute a bifurcation value μbif, such that if μ < μbif the zero solution is stable, and otherwise it is unstable.

With the help of our algorithm, these bifurcation values can be efficiently computed. Some examples are displayed in Table 1 . The presented values are for the case τ(x,y)=0.

Table 1

A

ξ1

μbif

A

ξ1

ξ2

μbif

0

0.1

15.81

1

0.4

0.3

17.68

0

0.2

7.98

1

0.5

0.3

10.24

0

0.3

5.35

1

0.6

0.3

8.41

0

0.4

4.06

1

0.7

0.3

7.94

0

0.5

3.56

1

0.8

0.3

8.09

Declarations

Ackknowledgement

This research was supported by a Marie Curie Intra European Fellowship, PIEF-GA-2013-629496.

Authors’ Affiliations

(1)
Institute of Stochastics, Johannes Kepler University, 4040 Linz, Austria
(2)
CEMAT/Instituto Superior Técnico, University of Lisbon, 1049 Lisboa, Portugal

References

  1. Wilson HR, Cowan JD: Excitatory and inhibitory interactions in localized populations of model neurons. Biophys J. 1972, 12 (1): 1-24.PubMedPubMed CentralView ArticleGoogle Scholar
  2. Amari SL: Dynamics of pattern formation in lateral-inhibition type neural fields. Biol Cybernet. 1977, 27 (2): 77-87.View ArticleGoogle Scholar
  3. Faye G, Faugeras O: Some theoretical and numerical results for delayed neural field equations. Physica D. 2010, 239 (9): 561-578.View ArticleGoogle Scholar

Copyright

© Lima and Buckwar 2015

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Advertisement