Abstract

In this paper, we analyze the invasion and extinction of activity in heterogeneous neural fields. We first consider the effects of spatial heterogeneities on the propagation of an invasive activity front. In contrast to previous studies of front propagation in neural media, we assume that the front propagates into an unstable rather than a metastable zero-activity state. For sufficiently localized initial conditions, the asymptotic velocity of the resulting pulled front is given by the linear spreading velocity, which is determined by linearizing about the unstable state within the leading edge of the front. One of the characteristic features of these so-called pulled fronts is their sensitivity to perturbations inside the leading edge. This means that standard perturbation methods for studying the effects of spatial heterogeneities or external noise fluctuations break down. We show how to extend a partial differential equation method for analyzing pulled fronts in slowly modulated environments to the case of neural fields with slowly modulated synaptic weights. The basic idea is to rescale space and time so that the front becomes a sharp interface whose location can be determined by solving a corresponding local Hamilton-Jacobi equation. We use steepest descents to derive the Hamilton-Jacobi equation from the original nonlocal neural field equation. In the case of weak synaptic heterogenities, we then use perturbation theory to solve the corresponding Hamilton equations and thus determine the time-dependent wave speed. In the second part of the paper, we investigate how time-dependent heterogenities in the form of extrinsic multiplicative noise can induce rare noise-driven transitions to the zero-activity state, which now acts as an absorbing state signaling the extinction of all activity. In this case, the most probable path to extinction can be obtained by solving the classical equations of motion that dominate a path integral representation of the stochastic neural field in the weak noise limit. These equations take the form of nonlocal Hamilton equations in an infinite-dimensional phase space.

Highlights

  • Reaction-diffusion equations based on the Fisher-Kolmogorov-Petrovskii-Piskunov (F-KPP) model and its generalizations have been used extensively to describe the spatial spread of invading species including plants, insects, diseases, and genes in terms of propagating fronts [1,2,3,4,5,6,7]

  • We have explored some of the consequences of having an unstable zero-activity state in a scalar neural field model

  • We showed how time-dependent corrections to the wave speed could be estimated using a Hamilton-Jacobi theory of sharp interface dynamics combined with perturbation theory

Read more

Summary

Introduction

Reaction-diffusion equations based on the Fisher-Kolmogorov-Petrovskii-Piskunov (F-KPP) model and its generalizations have been used extensively to describe the spatial spread of invading species including plants, insects, diseases, and genes in terms of propagating fronts [1,2,3,4,5,6,7]. One fundamental result in the theory of deterministic fronts is the difference between fronts propagating into a linearly unstable (zero) state and those propagating into a metastable state (a state that is linearly stable but nonlinearly unstable) In the latter case, the front has a unique velocity that is obtained by solving the associated partial differential equation (PDE) in traveling wave coordinates. As in the case of reaction diffusion systems, neural fields can exhibit a rich repertoire of wave phenomena, including solitary traveling fronts, pulses, and spiral waves [27,28,29] They have been used to model wave propagation in cortical slices [30, 31] and in vivo [32]. We begin by considering a neural field model that supports propagating pulled fronts, and determine the asymptotic wave speed by linearizing about the unstable state within the leading edge of the front (Section 2). The most probable path to extinction can be obtained by solving the classical equations of motion that dominate the path integral representation in the weak noise limit; these equations take the form of nonlocal Hamilton equations in an infinite-dimensional phase space (Section 6)

Neural fields with propagating pulled fronts
Hamilton-Jacobi dynamics and slow spatial heterogeneities
Calculation of wave speed
Numerical results
Path integral formulation of a stochastic neural field
Hamiltonian-Jacobi dynamics and population extinction in the weak-noise limit
Discussion
Fisher RA
Murray JD
40. Bressloff PC
49. Gourley SA
56. Gardiner CW
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call