Journal of Mathematical Neuroscience最新文献

Continuum neural field equations model the large-scale spatio-temporal dynamics of interacting neurons on a cortical surface. They have been extensively studied, both analytically and numerically, on bounded as well as unbounded domains. Neural field models do not require the specification of boundary conditions. Relatively little attention has been paid to the imposition of neural activity on the boundary, or to its role in inducing patterned states. Here we redress this imbalance by studying neural field models of Amari type (posed on one- and two-dimensional bounded domains) with Dirichlet boundary conditions. The Amari model has a Heaviside nonlinearity that allows for a description of localised solutions of the neural field with an interface dynamics. We show how to generalise this reduced but exact description by deriving a normal velocity rule for an interface that encapsulates boundary effects. The linear stability analysis of localised states in the interface dynamics is used to understand how spatially extended patterns may develop in the absence and presence of boundary conditions. Theoretical results for pattern formation are shown to be in excellent agreement with simulations of the full neural field model. Furthermore, a numerical scheme for the interface dynamics is introduced and used to probe the way in which a Dirichlet boundary condition can limit the growth of labyrinthine structures.

We examine a family of random firing-rate neural networks in which we enforce the neurobiological constraint of Dale's Law-each neuron makes either excitatory or inhibitory connections onto its post-synaptic targets. We find that this constrained system may be described as a perturbation from a system with nontrivial symmetries. We analyze the symmetric system using the tools of equivariant bifurcation theory and demonstrate that the symmetry-implied structures remain evident in the perturbed system. In comparison, spectral characteristics of the network coupling matrix are relatively uninformative about the behavior of the constrained system.

Layer II stellate cells in the medial enthorinal cortex (MEC) express hyperpolarisation-activated cyclic-nucleotide-gated (HCN) channels that allow for rebound spiking via an [Formula: see text] current in response to hyperpolarising synaptic input. A computational modelling study by Hasselmo (Philos. Trans. R. Soc. Lond. B, Biol. Sci. 369:20120523, 2013) showed that an inhibitory network of such cells can support periodic travelling waves with a period that is controlled by the dynamics of the [Formula: see text] current. Hasselmo has suggested that these waves can underlie the generation of grid cells, and that the known difference in [Formula: see text] resonance frequency along the dorsal to ventral axis can explain the observed size and spacing between grid cell firing fields. Here we develop a biophysical spiking model within a framework that allows for analytical tractability. We combine the simplicity of integrate-and-fire neurons with a piecewise linear caricature of the gating dynamics for HCN channels to develop a spiking neural field model of MEC. Using techniques primarily drawn from the field of nonsmooth dynamical systems we show how to construct periodic travelling waves, and in particular the dispersion curve that determines how wave speed varies as a function of period. This exhibits a wide range of long wavelength solutions, reinforcing the idea that rebound spiking is a candidate mechanism for generating grid cell firing patterns. Importantly we develop a wave stability analysis to show how the maximum allowed period is controlled by the dynamical properties of the [Formula: see text] current. Our theoretical work is validated by numerical simulations of the spiking model in both one and two dimensions.

Neural mass models provide a useful framework for modelling mesoscopic neural dynamics and in this article we consider the Jansen and Rit neural mass model (JR-NMM). We formulate a stochastic version of it which arises by incorporating random input and has the structure of a damped stochastic Hamiltonian system with nonlinear displacement. We then investigate path properties and moment bounds of the model. Moreover, we study the asymptotic behaviour of the model and provide long-time stability results by establishing the geometric ergodicity of the system, which means that the system-independently of the initial values-always converges to an invariant measure. In the last part, we simulate the stochastic JR-NMM by an efficient numerical scheme based on a splitting approach which preserves the qualitative behaviour of the solution.

We show that point-neuron models with a Heaviside firing rate function can be ill posed. More specifically, the initial-condition-to-solution map might become discontinuous in finite time. Consequently, if finite precision arithmetic is used, then it is virtually impossible to guarantee the accurate numerical solution of such models. If a smooth firing rate function is employed, then standard ODE theory implies that point-neuron models are well posed. Nevertheless, in the steep firing rate regime, the problem may become close to ill posed, and the error amplification, in finite time, can be very large. This observation is illuminated by numerical experiments. We conclude that, if a steep firing rate function is employed, then minor round-off errors can have a devastating effect on simulations, unless proper error-control schemes are used.

Presented here is a model of neural tissue in a conductive medium stimulated by externally injected currents. The tissue is described as a conductively isotropic bidomain, i.e. comprised of intra and extracellular regions that occupy the same space, as well as the membrane that divides them, and the injection currents are described as a pair of source and sink points. The problem is solved in three spatial dimensions and defined in spherical coordinates [Formula: see text]. The system of coupled partial differential equations is solved by recasting the problem to be in terms of the membrane and a monodomain, interpreted as a weighted average of the intra and extracellular domains. The membrane and monodomain are defined by the scalar Helmholtz and Laplace equations, respectively, which are both separable in spherical coordinates. Product solutions are thus assumed and given through certain transcendental functions. From these electrical potentials, analytic expressions for current density are derived and from those fields the magnetic flux density is calculated. Numerical examples are considered wherein the interstitial conductivity is varied, as well as the limiting case of the problem simplifying to two dimensions due to azimuthal independence. Finally, future modeling work is discussed.