Poster presentation | Open | Published:
Smale horseshoe structure in the firing rate model
BMC Neurosciencevolume 14, Article number: P152 (2013)
The firing rate model in the form of nonlinear integrodifferential equations can characterize spatiotemporal patterns of a continuum neural field. These patterns are associated with a wide range of neurobiological phenomena, such as persistent activity and propagating waves in neural networks.
To understand the substrates of neural circuitry that supports the localized stationary patterns, we study the existence of multi-bump pulse solutions of an integral equation that is the equilibrium equation of the firing rate model. If the integral coupling function, which describes the spatial connection among the network of neurons, is even and its positive half solves a second order linear ordinary differential equation (ODE), then the multi-bump pulse solutions of the integral equation are homoclinic solutions of a reversible fourth order ODE. It was known previously that the corresponding ODEs are conservative for a class of oscillatory and decaying coupling functions . We show that for the Amari-type firing rate model, which assumes the Heaviside firing rate function and "Mexican hat" coupling functions that decay exponentially, the corresponding ODEs are also conservative. Then, by analyzing the configurations of stable and unstable manifolds (see the figure below) in the corresponding reversible, conservative ODE system, we establish the existence of the Smale horseshoe for an open set of model parameters such as the decay rate of the coupling and the threshold of the Heaviside firing rate function. Consequently, even for the Amari-type firing rate model there are countably many symmetric and asymmetric multi-bump stationary pulse solutions as well as spatially chaotic stationary solutions. Furthermore, the robustness of the Smale horseshoe implies that similar solutions also exist for nonsaturating piecewise-linear firing rate functions with small gains and for smooth firing rate functions that are "close" to the Heaviside or piecewise-linear case.
Laing CR, Troy WC, Gutkin B, Ermentrout GB: Multiple bumps in a neuronal model of working memory. SIAM J Appl Math. 2002, 63: 62-97. 10.1137/S0036139901389495.