Exploiting the Hamiltonian structure of a neural field model. (English) Zbl 1203.37034

This work is concerned with the following model \[ \frac{\partial u(x,t)}{\partial t} \, = \, - u(x,t) + \int_{-\infty}^{+\infty} w(x-y) f[u(y,t)] dy, \] where \(u(x,t)\) is the average voltage, or activity level, of a neuronal population of spatial position \(x\) and time \(t\). The coupling function \(w(x)\) is the distance-dependent strength of connectivity between neuronal elements and is given by \[ w(x) \, = \, e^{-b|x|} \left( b \sin |x| + \cos x \right), \] with \(b\) a parameter which governs the rate at which oscillations in \(w\) decay with distance. The firing rate function \(f(u)\) models neurons firing once threshold is reached and tends to a maximal limit as the stimulus is increased. It is given by: \[ f(u) \, = \, 2 \, \exp \left[-r/(u-\theta)^2 \right] \Theta(u-\theta), \] where the parameter \(\theta\) is the firing threshold, \(r\) is the steepness parameter and \(\Theta\) is the Heaviside step function. We recall that the Heaviside step function is defined by \[ \Theta(x) \, = \, \int_{-\infty}^{x} \delta(t) \, dt, \] where \(\delta(t)\) is the Dirac delta function. \(\Theta(x)\) is zero for \(x<0\) and it is \(1\) for \(x>0\). This model has been studied in [C. R. Laing et al., SIAM J. Appl. Math., 63, No.1, 62–97 (2002; Zbl 1017.45006)] and [C. R. Laing, W.C. Troy, SIAM J. Appl. Dyn. Syst., 2, No. 3, 487–516 (2003; Zbl 1088.34011)].
The authors consider time-independent solutions of the equation which are stationary and spatially-localised, that is, \[ \lim_{|x| \to \infty} (u,u',u'',u''') \, = \, (0,0,0,0). \] These solutions satisfy the integral equation \[ u(x) \, = \, \int_{-\infty}^{\infty} w(x-y) f[u(y)] dy, \] which can be transformed to the following differential equation by use of Fourier transforms: \[ u'''' - 2(b^2-1) u'' + (b^2+1)^2 u \, = \, 4b(b^2+1) f(u). \tag{eq9} \] This equation can be written as a system of four first-order ODEs with the following properties:
1) the linearisation of the system at the origin has the four eigenvalues \(\pm b\pm i\), i.e., the origin is a bifocus with a two-dimensional unstable manifold and a two-dimensional stable manifold.
2) the system can be written as a Hamiltonian system.
3) this Hamiltonian system is invariant under a space-reversing symmetry.
Using all these properties the authors derive a real scalar function whose zeros give rise to the single-bump homoclinic orbits to the fixed point at the origin of (eq9).
The authors numerically integrate the equations so as to examine the unexpected disappearance of stable homoclinic orbits in certain regions of the parameter space. The solution curve of stable homoclinic orbits breaks when the firing rate function is sufficiently steep.


37C29 Homoclinic and heteroclinic orbits for dynamical systems
92C20 Neural biology


HomCont; AUTO
Full Text: DOI Link


[1] Champneys, A.R., Homoclinic orbits in reversible systems and their applications in mechanics, fluids and optics, Physica D, 112, 158-186, (1998) · Zbl 1194.37154
[2] Hunt, G.W.; Peletier, M.A.; Champneys, A.R.; Woods, P.D.; Ahmer Wadee, M.; Budd, C.J.; Lord, G.J., Cellular buckling in long structures, Nonlinear dyn., 21, 3-29, (2000) · Zbl 0974.74024
[3] Tlidi, M.; Mandel, P.; Lefever, R., Localized structures and localized patterns in optical bistability, Phys. rev. lett., 73, 640-643, (1994)
[4] Umbanhowar, P.B.; Melo, F.; Swinney, H.L., Localized excitations in a vertically vibrated granular layer, Nature, 382, 793, (1996)
[5] Dawes, J.H.P., Localized convection cells in the presence of a vertical magnetic field, J. fluid mech., 570, 385-406, (2007) · Zbl 1105.76026
[6] Guo, Y.; Chow, C., Existence and stability of standing pulses in neural networks: I. existence, SIAM J. appl. dyn. syst., 4, 217-248, (2005) · Zbl 1109.34002
[7] Laing, C.R.; Troy, W.C.; Gutkin, B.; Ermentrout, G.B., Multiple bumps in a neuronal model of working memory, SIAM J. appl. math., 63, 1, 62-97, (2002) · Zbl 1017.45006
[8] Amari, S., Dynamics of pattern formation in lateral-inhibition type neural fields, Biol. cybern., 27, 2, 77-87, (1977) · Zbl 0367.92005
[9] Ermentrout, G.B., Neural networks as spatio – temporal pattern – forming systems, Rep. prog. phys., 61, 353-430, (1998)
[10] Pinto, D.J.; Ermentrout, G.B., Spatially structured activity in synaptically coupled neuronal networks: II. lateral inhibition and standing pulses, SIAM J. appl. math., 62, 1, 226-243, (2001) · Zbl 1070.92506
[11] Burke, J.; Knobloch, E., Localized states in the generalised swift – hohenberg equation, Phys. rev. E, 73, 056211, (2006)
[12] Burke, J.; Knobloch, E., Homoclinic snaking: structure and stability, Chaos, 17, 037102, (2007) · Zbl 1163.37317
[13] Burke, J.; Knobloch, E., Snakes and ladders: localized states in the swift – hohenberg equation, Phys. lett. A, 360, 681-688, (2007) · Zbl 1236.35144
[14] Knobloch, E., Spatially localized structures in dissipative systems: open problems, Nonlinearity, 21, T45-T60, (2008) · Zbl 1139.35002
[15] Peletier, L.A.; Troy, W.C., Spatial patterns: higher order models in physics and mechanics, (2001), Birkhauser Boston · Zbl 0872.34032
[16] Lloyd, D.J.B; Sandstede, B.; Avitabile, D.; Champneys, A.R., Localized hexagon patterns of the planar swift – hohenberg equation, SIAM J. appl. dyn. syst., 7, 3, 1049-1100, (2008) · Zbl 1168.35311
[17] Laing, C.R.; Troy, W.C., PDE methods for nonlocal models, SIAM J. appl. dyn. syst., 2, 3, 487-516, (2003) · Zbl 1088.34011
[18] Gutkin, B.S.; Ermentrout, G.B.; O’Sullivan, J., Layer 3 patchy recurrent excitatory connections may determine the spatial organization of sustained activity in the primate prefrontal cortex, Neurocomputing, 32-33, 391-400, (2000)
[19] Levitt, J.B.; Lewis, D.A.; Yoshioka, T.; Lund, J.S., Topography of pyramidal neuron intrinsic connections in macaque monkey prefrontal cortex (areas 9 and 46), J. comput. neuro., 338, 360-376, (1993)
[20] Lewis, D.A.; Anderson, S.A., The functional architecture of the prefrontal cortex and schizophrenia, Psych. med., 25, 887-894, (1995)
[21] A.J. Elvin, Pattern formation in a neural field model. Ph.D. Thesis, Massey University, 2008. http://hdl.handle.net/10179/717
[22] Elvin, A.J.; Laing, C.R.; Roberts, M.G., Transient Turing patterns in a neural field model, Phys. rev. E, 79, 011911, (2009)
[23] Coombes, S.; Lord, G.J.; Owen, M.R., Waves and bumps in neuronal networks with axo-dendritic synaptic interactions, Physica D, 178, 3-4, 219-241, (2003) · Zbl 1013.92006
[24] Coombes, S., Waves, bumps, and patterns in neural field theories, Biol. cybern., 93, 2, 91-108, (2005) · Zbl 1116.92012
[25] Coombes, S.; Venkov, N.A.; Shiau, L.; Bojak, I.; Liley, D.T.J.; Laing, C.R., Modeling electrocortical activity through improved local approximations of integral neural field equations, Phys. rev. E, 76, 5, 51901, (2007)
[26] Jirsa, V.K.; Haken, H., Field theory of electromagnetic brain activity, Phys. rev. lett., 77, 5, 960-963, (1996)
[27] E.J. Doedel, R.C. Paffenroth, A.R. Champneys, T.F. Fairgrieve, Y.A. Kuznetsov, B.E. Oldeman, B. Sandstede, X. Wang, AUTO 2000: Continuation and bifurcation software for ordinary differential equations (with HomCont), 2000
[28] Drazin, P.G., Nonlinear systems, (2002), Cambridge University Press · Zbl 0753.34001
[29] Glendinning, P.; Sparrow, C., T-points: A codimension two heteroclinic bifurcation, J. stat. phys., 43, 3-4, (1986) · Zbl 0635.58031
[30] Coombes, S.; Owen, M.R., Bumps, breathers, and waves in a neural network with spike frequency adaptation, Phys. rev. lett., 94, (2005)
[31] Kilpatrick, Z.P.; Bressloff, P.C., Effects of synaptic depression and adaptation on spatiotemporal dynamics of an excitatory neuronal network, Physica D, 239, 9, 547-560, (2010) · Zbl 1186.92010
[32] Coombes, S.; Owen, M.R., Evans functions for integral neural field equations with heaviside firing rate function, SIAM J. appl. dyn. syst., 3, 4, 574-600, (2004) · Zbl 1067.92019
[33] Laing, C.; Coombes, S., The importance of different timings of excitatory and inhibitory pathways in neural field models, Network: comp. neural sys., 17, 2, 151-172, (2006)
[34] Owen, M.R.; Laing, C.R.; Coombes, S., Bumps, rings, and spots in a two-dimensional neural field: splitting and rotational instabilities, New J. phys., 9, 378, (2007)
[35] Traub, R.D.; Jefferys, J.G.; Miles, R., Analysis of the propagation of disinhibition-induced after-discharges along the guinea-pig hippocampal slice in vitro, J. physiol., 472, 267-287, (1993)
[36] Wu, J.-Y.; Huang, X.; Zhang, C., Propagating waves of activity in the neocortex: what they are, what they do, The neuroscientist, 14, 487-502, (2008)
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.