Numerical investigation of a neural field model including dendritic processing

Daniele Avitabile, Stephen Coombes, Pedro M. Lima

Research output: Contribution to JournalArticleAcademicpeer-review

Abstract

We consider a simple neural field model in which the state variable is dendritic voltage, and in which somas form a continuous one-dimensional layer. This neural field model with dendritic processing is formulated as an integro-differential equation. We introduce a computational method for approximating solutions to this nonlocal model, and use it to perform numerical simulations for neuro-biologically realistic choices of anatomical connectivity and nonlinear firing rate function. For the time discretisation we adopt an Implicit-Explicit (IMEX) scheme; the space discretisation is based on a finitedi fference scheme to approximate the diffusion term and uses the trapezoidal rule to approximate integrals describing the nonlocal interactions in the model. We prove that the scheme is of first-order in time and second order in space, and can be efficiently implemented if the factorisation of a small, banded matrix is precomputed. By way of validation we compare the outputs of a numerical realisation to theoretical predictions for the onset of a Turing pattern, and to the speed and shape of a travelling front for a specific choice of Heaviside firing rate. We find that theory and numerical simulations are in excellent agreement.

Original languageEnglish
Pages (from-to)271-290
Number of pages20
JournalJournal of Computational Dynamics
Volume7
Issue number2
DOIs
Publication statusPublished - Dec 2020

Keywords

  • Cable equation
  • Finite difference scheme
  • Integral equations
  • Neural fields
  • Travelling waves

Fingerprint Dive into the research topics of 'Numerical investigation of a neural field model including dendritic processing'. Together they form a unique fingerprint.

  • Cite this