Abstract
Neural Field Equations (NFE) are intended to model the synaptic interactions between neurons in a continuous neural network, called a neural field. This kind of integro-differential equations proved to be a useful tool to describe the spatiotemporal neuronal activity from a macroscopic point of view, allowing the study of a wide variety of neurobiological phenomena, such as the sensory stimuli processing. The present article aims to study the effects of additive noise in one- and two-dimensional neural fields, while taking into account finite axonal velocity and an external stimulus. A Galerkin-type method is presented, which applies Fast Fourier Transforms to optimise the computational effort required to solve these equations. The explicit Euler-Maruyama scheme is implemented to obtain the stochastic numerical solution. An open-source numerical solver written in Julia was developed to simulate the neural fields in study.






Similar content being viewed by others
Code availability
A numerical solver written in Julia, whose purpose is to solve the types of NFEs discussed in this paper, is published in Julia's library and available for installation, more details at https://github.com/tiagoseq/NeuralFieldEq.jl. Technical details of code usage are described in Sequeira (2021). Julia proved to be a great tool to simulate the addressed NF, and allied to the \(\mathcal {RFFT}\) routines it was possible to work in strongly delayed scenarios, when the computational effort increases significantly. The solver is efficient enough to not run out of memory, even in the stochastic case, where we had to compute a large number of trajectories for each noise level. In all simulations carried out, the computing time never exceeded 40 minutes. Note that the computations were performed in a Laptop with a 1.30GHz Intel(R) Core(TM) i7 CPU processor and 16GB memory RAM.
References
Amari, S. L. (1977). Dynamics of pattern formation in lateral.inhibition type neural fields. Biological Cybernetics, 27, 77–87
Bell, J., & Cosner, C. (1983). Threshold conditions for a diffusive model of a myelinated axon. Journal of Mathematical Biology 18(1), 39–52. https://doi.org/10.1007/BF00275909.
Bell, J. (1984). Behaviour of some models of myelinated axons. Mathematical Medicine and Biology: A Journal of the IMA 1(2), 149-167. https://doi.org/10.1093/imammb/1.2.149
Bressloff, P., Cowan, J., Golubitsky, M., Thomas, P., & Wiener, M. (2002). What geometric visual hallucinations tell us about the visual cortex. Neural computation, 14, 473–91. https://doi.org/10.1162/089976602317250861.
Chi, H., Bell, J., & Hassard, B. (1986). Numerical solution of a nonlinear advance-delay-differential equation from nerve conduction theory. Journal of mathematical biology, 24, 583–601. https://doi.org/10.1007/BF00275686.
Clearfield, M., Dineva, E., Smith, L., Diedrich, F., & Thelen, E. (2009). Cue salience and infant perseverative reaching: Tests of the dynamic field theory. Developmental science, 12, 26–40b. https://doi.org/10.1111/j.1467-7687.2008.00769.x.
Coombes, S. (2010). Large-scale neural dynamics: Simple and complex. NeuroImage, 52, 731–9. https://doi.org/10.1016/j.neuroimage.2010.01.045.
Desmedt, J. E., & Cheron, G. (1980). Central somatosensory conduction in man: Neural generators and interpeak latencies of the far-field components recorded from neck and right or left scalp and earlobes. Electroencephalography and Clinical Neurophysiology, 50(5), 382–403. https://doi.org/10.1016/0013-4694(80)90006-1.
Erlhagen, W., & Bicho, E. (2006). The dynamic neural field approach to cognitive robotics. Journal of neural engineering, 3, 36–54. https://doi.org/10.1088/1741-2560/3/3/R02.
Faye, G., & Faugeras, O. (2010). Some theoretical and numerical results for delayed neural field equations. Physica D: Nonlinear Phenomena, 239, 561–578. https://doi.org/10.1016/j.physd.2010.01.010.
Faye, G., & Faugeras, O. (2010). New results for delayed neural field equations. Cinquiéme conférence pléniére française de Neurosciences omputationnelles. Neurocomp, 10
Folias, S. E., & Bressloff, P. C. (2005). Breathers in two-dimensional neural media. Physical Review Letters, 95, 208107.
Funahashi, S., Bruce, C., & Goldman-Rakic, P. S. (1989). Mnemonic coding of visual space in the monkeys dorsolateral prefrontal cortex. Journal of Neurophysiology, 61, 1–9. https://doi.org/10.1152/jn.1989.61.2.331.
Gerstner, W., & Kistler, W. (2002). Spiking Neuron Models. Single Neurons, Populations, Plasticity. Cambridge Univeristy Press, Cambridge.
Hahnloser, R., Xie, X., & Seung, H. (2001). A theory of neural integration in the head-direction system. Advances in Neural Information Processing Systems, 14.
Higham, D. J. (2001). An algorithmic introduction to numerical simulation of stochastic differential equations. SIAM Review, 43(3), 525-546. https://doi.org/10.1137/S0036144500378302
Hutt, A., & Rougier, N. (2010). Activity spread and breathers induced by finite transmission speeds in two-dimensional neural fields. Physical Review E, 82(5), 055701. https://doi.org/10.1103/PhysRevE.82.055701
Hutt, A., & Rougier, N. (2013). Numerical simulation scheme of one-and two dimensional neural fields involving space-dependent delays. In Neural Fields (pp. 175-185). Springer, Berlin, Heidelberg. https://doi.org/10.1007/978-3-642-54593-16
Kilpatrick, Z., & Ermentrout, B. (2012). Wandering bumps in stochastic neural fields. SIAM Journal on Applied Dynamical Systems, 12, https://doi.org/10.1137/120877106.
Kuehn, C. G. & Riedler, M. (2014). Large deviations for nonlocal stochastic neural fields. Journal of Mathematical Neuroscience 4, 1 (2014). https://doi.org/10.1186/2190-8567-4-1
Laing, C. R., Troy, W. C., Gutkin, B., & Ermentrout, G. B. (2002). Multiple bumps in a neuronal model of working memory. SIAM Journal on Applied Mathematics, 63, https://doi.org/10.1137/S0036139901389495
Liley, D., Cadusch, P. J., & Dafilis, M. P. (2009). A spatially continuous mean field theory of electrocortical activity. Network: Computation in Neural Systems, 13, https://doi.org/10.1088/0954-898X/14/2/601
Lima, P. M. (2019). Numerical investigation of stochastic neural field equations. In: Sing, V.K. (ed.). Advances in Mathematical Methods and High Performance Computing, 3 51-67. New York.
Lund, J., Angelucci, A., & Bressloff, P. (2003). Anatomical substrates for functional columns in macaque monkey primary visual cortex. Cerebral cortex (New York, N.Y. : 1991) 13, 15-24. https://doi.org/10.1093/cercor/13.1.15
Lima, P., & Buckwar, E. (2015). Numerical solution of the neural field equation in the two-dimensional case. SIAM Journal on Scientific Computing, 37, 962–979. https://doi.org/10.1137/15M1022562.
Nichols, E., & Hutt, A. (2015). Neural field simulator: two-dimensional spatio-temporal dynamics involving finite transmission speed. Frontiers in Neuroinformatics, 9, 25. https://doi.org/10.3389/fninf.2015.00025.
Parsons, B., Novich, S., & Eagleman, D. (2013). Motor-sensory recalibration modulates perceived simultaneity of cross-modal events at different distances. Frontiers in psychology, 4, 46. https://doi.org/10.3389/fpsyg.2013.00046.
Shepherd, G. M. (1994). Neurobiology. New York: Oxford University Press.
Sequeira, T. (2021). Neuralfieldeq.jl: A flexible solver to compute neural field equations in several scenarios - in press. Journal of Open Source Software.
Thelen, E., Schöner, G., Scheier, C., & Smith, L. B. (2001). The dynamics of embodiment: A field theory of infant perseverative reaching. Behavioral and Brain Sciences, 24(1), 1–34 https://doi.org/10.1017/s0140525x01003910
Wilson, H., & Cowan, J. (1973). A mathematical theory of the functional dynamics of cortical and thalamic nervous tissue. Kybernetik, 13, 55-80. https://link.springer.com/article/10.1007/BF00288786
Acknowledgements
The second author acknowledges the financial support of the portuguese FCT (Fundação para a Ciência e Tecnologia), through projects UIDB/04621/2020, UIDP/04621/2020 and PTDC/MAT-APL/31393/2017. Both authors are grateful to the anonymous reviewers, whose comments and suggestions helped to improve the paper.
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflicts of interest
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
Additional information
Action Editor: Alain Destexhe.
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
About this article
Cite this article
Sequeira, T.F., Lima, P.M. Numerical simulations of one- and two-dimensional stochastic neural field equations with delay. J Comput Neurosci 50, 299–311 (2022). https://doi.org/10.1007/s10827-022-00816-w
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10827-022-00816-w