1
|
Sequeira TF, Lima PM. Numerical simulations of one- and two-dimensional stochastic neural field equations with delay. J Comput Neurosci 2022; 50:299-311. [PMID: 35618864 DOI: 10.1007/s10827-022-00816-w] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2021] [Revised: 02/16/2022] [Accepted: 03/09/2022] [Indexed: 11/29/2022]
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.
Collapse
Affiliation(s)
- Tiago F Sequeira
- Center for Computational and Stochastic Mathematics, Instituto Superior Técnico, Universidade de Lisboa, Av. Rovisco Pais, Lisboa, 1049-001, Portugal
| | - Pedro M Lima
- Center for Computational and Stochastic Mathematics, Instituto Superior Técnico, Universidade de Lisboa, Av. Rovisco Pais, Lisboa, 1049-001, Portugal.
| |
Collapse
|
2
|
Kulikova MV, Lima PM, Kulikov GY. Sequential method for fast neural population activity reconstruction in the cortex from incomplete noisy measurements. Comput Biol Med 2021; 141:105103. [PMID: 34959112 DOI: 10.1016/j.compbiomed.2021.105103] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2021] [Revised: 11/29/2021] [Accepted: 11/29/2021] [Indexed: 11/03/2022]
Abstract
During recent years there has been a growing interest in stochastic dynamic neural fields employed for modeling and predictions in biomedical and technical systems. In this paper, given some incomplete noisy data available from sensors, we propose and explore a state estimation method for fast restorations of membrane potential in the cortex based on such measurements and the Amari equation used for simulations of neural population activity in a stochastic setting. Our novel technique relies upon a Galerkin-type spectral approximation utilized within the conventional state-space approach. Translating a stochastic system into its state-space form creates a straightforward and fruitful way to the data-driven parameter estimation, filtering, prediction and smoothing. The present study is particularly focused on establishing a nonlinear stochastic Galerkin-spectral-approximation-induced system of large size, which is further estimated by the traditional extended Kalman filter (EKF). The efficiency of calculations is the main purpose of our research. That is why the fast filtering solution devised is based on processing the incoming data incrementally, that is, by processing measurements one at a time, rather than handling them as a unified high-dimensional vector. Such sequential filters suit well for dealing with large data sets as well as with real-time on-line computations. Also, their derivation and substantiation is of great interest in the context of neural network training because of large stochastic systems arisen there. In comparison to the batch filtering, our novel algorithm reduces the computational cost of membrane potential reconstructions in terms of the amount of grid nodes N accepted in the underlying spacial discretization, significantly. Apart from its computation efficiency, this sequential method is more robust to round-off errors committed within a computer-based finite precision arithmetics than the classical EKF because of the (N × N)-matrix inversion elimination from such membrane potential calculations. The superior performance of our technique is examined and confirmed in comparison to the batch one on two known scenarios in the dynamic neural field modeling.
Collapse
Affiliation(s)
- M V Kulikova
- CEMAT, Instituto Superior Técnico, Universidade de Lisboa, Av. Rovisco Pais, 1049-001, Lisboa, Portugal.
| | - P M Lima
- CEMAT, Instituto Superior Técnico, Universidade de Lisboa, Av. Rovisco Pais, 1049-001, Lisboa, Portugal.
| | - G Yu Kulikov
- CEMAT, Instituto Superior Técnico, Universidade de Lisboa, Av. Rovisco Pais, 1049-001, Lisboa, Portugal.
| |
Collapse
|
3
|
Lima PM, Erlhagen W, Kulikov GY, Kulikova MV. Mathematical Modeling of Working Memory in the Presence of Random Disturbance using Neural Field Equations. EPJ WEB OF CONFERENCES 2021. [DOI: 10.1051/epjconf/202124801021] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this paper, we describe a neural field model which explains how a population of cortical neurons may encode in its firing pattern simultaneously the nature and time of sequential stimulus events. Moreover, we investigate how noise-induced perturbations may affect the coding process. This is obtained by means of a two-dimensional neural field equation, where one dimension represents the nature of the event (for example, the color of a light signal) and the other represents the moment when the signal has occurred. The additive noise is represented by a Q-Wiener process. Some numerical experiments reported are carried out using a computational algorithm for two-dimensional stochastic neural field equations.
Collapse
|
4
|
Romeo A, Supèr H. Bump competition and lattice solutions in two-dimensional neural fields. Neural Netw 2017; 94:141-158. [PMID: 28779599 DOI: 10.1016/j.neunet.2017.07.003] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/28/2016] [Revised: 05/19/2017] [Accepted: 07/02/2017] [Indexed: 10/19/2022]
Abstract
Some forms of competition among activity bumps in a two-dimensional neural field are studied. First, threshold dynamics is included and rivalry evolutions are considered. The relations between parameters and dominance durations can match experimental observations about ageing. Next, the threshold dynamics is omitted from the model and we focus on the properties of the steady-state. From noisy inputs, hexagonal grids are formed by a symmetry-breaking process. Particular issues about solution existence and stability conditions are considered. We speculate that they affect the possibility of producing basis grids which may be combined to form feature maps.
Collapse
Affiliation(s)
- August Romeo
- Department of Cognition, Development and Educational Psychology, Faculty of Psychology, University of Barcelona, Spain
| | - Hans Supèr
- Department of Cognition, Development and Educational Psychology, Faculty of Psychology, University of Barcelona, Spain; Institut de Neurociències, University of Barcelona, Spain; Catalan Institution for Research and Advanced Studies (ICREA), Spain.
| |
Collapse
|
5
|
Visser S, Nicks R, Faugeras O, Coombes S. Standing and travelling waves in a spherical brain model: The Nunez model revisited. PHYSICA D. NONLINEAR PHENOMENA 2017; 349:27-45. [PMID: 28626276 PMCID: PMC5421190 DOI: 10.1016/j.physd.2017.02.017] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/09/2016] [Revised: 02/27/2017] [Accepted: 02/28/2017] [Indexed: 06/07/2023]
Abstract
The Nunez model for the generation of electroencephalogram (EEG) signals is naturally described as a neural field model on a sphere with space-dependent delays. For simplicity, dynamical realisations of this model either as a damped wave equation or an integro-differential equation, have typically been studied in idealised one dimensional or planar settings. Here we revisit the original Nunez model to specifically address the role of spherical topology on spatio-temporal pattern generation. We do this using a mixture of Turing instability analysis, symmetric bifurcation theory, centre manifold reduction and direct simulations with a bespoke numerical scheme. In particular we examine standing and travelling wave solutions using normal form computation of primary and secondary bifurcations from a steady state. Interestingly, we observe spatio-temporal patterns which have counterparts seen in the EEG patterns of both epileptic and schizophrenic brain conditions.
Collapse
Affiliation(s)
- S. Visser
- School of Mathematical Sciences, University of Nottingham, NG7 2RD, UK
- Wellcome Trust Centre for Biomedical Modelling and Analysis, RILD Building, University of Exeter, EX2 5DW, UK
| | - R. Nicks
- School of Mathematical Sciences, University of Nottingham, NG7 2RD, UK
| | - O. Faugeras
- INRIA Sophia Antipolis Mediterannee, 2004 Route Des Lucioles, Sophia Antipolis, 06410, France
| | - S. Coombes
- School of Mathematical Sciences, University of Nottingham, NG7 2RD, UK
| |
Collapse
|
6
|
Nichols EJ, Hutt A. Neural field simulator: two-dimensional spatio-temporal dynamics involving finite transmission speed. Front Neuroinform 2015; 9:25. [PMID: 26539105 PMCID: PMC4611063 DOI: 10.3389/fninf.2015.00025] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/20/2015] [Accepted: 10/02/2015] [Indexed: 12/14/2022] Open
Abstract
Neural Field models (NFM) play an important role in the understanding of neural population dynamics on a mesoscopic spatial and temporal scale. Their numerical simulation is an essential element in the analysis of their spatio-temporal dynamics. The simulation tool described in this work considers scalar spatially homogeneous neural fields taking into account a finite axonal transmission speed and synaptic temporal derivatives of first and second order. A text-based interface offers complete control of field parameters and several approaches are used to accelerate simulations. A graphical output utilizes video hardware acceleration to display running output with reduced computational hindrance compared to simulators that are exclusively software-based. Diverse applications of the tool demonstrate breather oscillations, static and dynamic Turing patterns and activity spreading with finite propagation speed. The simulator is open source to allow tailoring of code and this is presented with an extension use case.
Collapse
Affiliation(s)
- Eric J. Nichols
- Team Neurosys, Loria, Centre National de la Recherche Scientifique, INRIA, UMR no. 7503, Université de LorraineNancy, France
| | | |
Collapse
|
7
|
Detorakis GI, Rougier NP. Structure of receptive fields in a computational model of area 3b of primary sensory cortex. Front Comput Neurosci 2014; 8:76. [PMID: 25120461 PMCID: PMC4112916 DOI: 10.3389/fncom.2014.00076] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/06/2014] [Accepted: 06/29/2014] [Indexed: 11/24/2022] Open
Abstract
In a previous work, we introduced a computational model of area 3b which is built upon the neural field theory and receives input from a simplified model of the index distal finger pad populated by a random set of touch receptors (Merkell cells). This model has been shown to be able to self-organize following the random stimulation of the finger pad model and to cope, to some extent, with cortical or skin lesions. The main hypothesis of the model is that learning of skin representations occurs at the thalamo-cortical level while cortico-cortical connections serve a stereotyped competition mechanism that shapes the receptive fields. To further assess this hypothesis and the validity of the model, we reproduced in this article the exact experimental protocol of DiCarlo et al. that has been used to examine the structure of receptive fields in area 3b of the primary somatosensory cortex. Using the same analysis toolset, the model yields consistent results, having most of the receptive fields to contain a single region of excitation and one to several regions of inhibition. We further proceeded our study using a dynamic competition that deeply influences the formation of the receptive fields. We hypothesized this dynamic competition to correspond to some form of somatosensory attention that may help to precisely shape the receptive fields. To test this hypothesis, we designed a protocol where an arbitrary region of interest is delineated on the index distal finger pad and we either (1) instructed explicitly the model to attend to this region (simulating an attentional signal) (2) preferentially trained the model on this region or (3) combined the two aforementioned protocols simultaneously. Results tend to confirm that dynamic competition leads to shrunken receptive fields and its joint interaction with intensive training promotes a massive receptive fields migration and shrinkage.
Collapse
Affiliation(s)
| | - Nicolas P Rougier
- INRIA Bordeaux Sud-Ouest Bordeaux, France ; Institut des Maladies Neurodégénératives, Université de Bordeaux, Centre National de la Recherche Scientifique, UMR 5293 Bordeaux, France ; LaBRI, Université de Bordeaux, Institut Polytechnique de Bordeaux, Centre National de la Recherche Scientifique, UMR 5800 Talence, France
| |
Collapse
|
8
|
Nichols E, Green K, Hutt A, van Veen L. Two-dimensional patterns in neural fields subject to finite transmission speed. BMC Neurosci 2014. [PMCID: PMC4125084 DOI: 10.1186/1471-2202-15-s1-p16] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022] Open
|
9
|
Rougier NP, Fix J. DANA: distributed numerical and adaptive modelling framework. NETWORK (BRISTOL, ENGLAND) 2012; 23:237-253. [PMID: 22994650 DOI: 10.3109/0954898x.2012.721573] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/01/2023]
Abstract
DANA is a python framework ( http://dana.loria.fr ) whose computational paradigm is grounded on the notion of a unit that is essentially a set of time dependent values varying under the influence of other units via adaptive weighted connections. The evolution of a unit's value are defined by a set of differential equations expressed in standard mathematical notation which greatly ease their definition. The units are organized into groups that form a model. Each unit can be connected to any other unit (including itself) using a weighted connection. The DANA framework offers a set of core objects needed to design and run such models. The modeler only has to define the equations of a unit as well as the equations governing the training of the connections. The simulation is completely transparent to the modeler and is handled by DANA. This allows DANA to be used for a wide range of numerical and distributed models as long as they fit the proposed framework (e.g. cellular automata, reaction-diffusion system, decentralized neural networks, recurrent neural networks, kernel-based image processing, etc.).
Collapse
Affiliation(s)
- Nicolas P Rougier
- INRIA Bordeaux - Sud Ouest, 351, Cours de la Libération, 33405 Talence Cedex, France.
| | | |
Collapse
|
10
|
Detorakis GI, Rougier NP. A neural field model of the somatosensory cortex: formation, maintenance and reorganization of ordered topographic maps. PLoS One 2012; 7:e40257. [PMID: 22808127 PMCID: PMC3395710 DOI: 10.1371/journal.pone.0040257] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/28/2012] [Accepted: 06/04/2012] [Indexed: 11/18/2022] Open
Abstract
We investigate the formation and maintenance of ordered topographic maps in the primary somatosensory cortex as well as the reorganization of representations after sensory deprivation or cortical lesion. We consider both the critical period (postnatal) where representations are shaped and the post-critical period where representations are maintained and possibly reorganized. We hypothesize that feed-forward thalamocortical connections are an adequate site of plasticity while cortico-cortical connections are believed to drive a competitive mechanism that is critical for learning. We model a small skin patch located on the distal phalangeal surface of a digit as a set of 256 Merkel ending complexes (MEC) that feed a computational model of the primary somatosensory cortex (area 3b). This model is a two-dimensional neural field where spatially localized solutions (a.k.a. bumps) drive cortical plasticity through a Hebbian-like learning rule. Simulations explain the initial formation of ordered representations following repetitive and random stimulations of the skin patch. Skin lesions as well as cortical lesions are also studied and results confirm the possibility to reorganize representations using the same learning rule and depending on the type of the lesion. For severe lesions, the model suggests that cortico-cortical connections may play an important role in complete recovery.
Collapse
Affiliation(s)
- Georgios Is. Detorakis
- INRIA CNRS: UMR 7503 Université Henri Poincaré - Nancy I Université Nancy II Institut National Polytechnique de Lorraine, Nancy, France
| | - Nicolas P. Rougier
- INRIA CNRS: UMR 7503 Université Henri Poincaré - Nancy I Université Nancy II Institut National Polytechnique de Lorraine, Nancy, France
- * E-mail:
| |
Collapse
|