Numerical methods for solving moment equations in kinetic theory of neuronal network dynamics

Recently developed kinetic theory and related closures for neuronal network dynamics have been demonstrated to be a powerful theoretical framework for investigating coarse-grained dynamical properties of neuronal networks. The moment equations arising from the kinetic theory are a system of (1+1)-dimensional nonlinear partial differential equations (PDE) on a bounded domain with nonlinear boundary conditions. The PDEs themselves are self-consistently specified by parameters which are functions of the boundary values of the solution. The moment equations can be stiff in space and time. Numerical methods are presented here for efficiently and accurately solving these moment equations. The essential ingredients in our numerical methods include: (i) the system is discretized in time with an implicit Euler method within a spectral deferred correction framework, therefore, the PDEs of the kinetic theory are reduced to a sequence, in time, of boundary value problems (BVPs) with nonlinear boundary conditions; (ii) a set of auxiliary parameters is introduced to recast the original BVP with nonlinear boundary conditions as BVPs with linear boundary conditions - with additional algebraic constraints on the auxiliary parameters; (iii) a careful combination of two Newton's iterates for the nonlinear BVP with linear boundary condition, interlaced with a Newton's iterate for solving the associated algebraic constraints is constructed to achieve quadratic convergence for obtaining the solutions with self-consistent parameters. It is shown that a simple fixed-point iteration can only achieve a linear convergence for the self-consistent parameters. The practicability and efficiency of our numerical methods for solving the moment equations of the kinetic theory are illustrated with numerical examples. It is further demonstrated that the moment equations derived from the kinetic theory of neuronal network dynamics can very well capture the coarse-grained dynamical properties of integrate-and-fire neuronal networks.

[1]  E. Çinlar.  On the superposition of m-dimensional point processes , 1968, Journal of Applied Probability.

[2]  G. Toscani,et al.  Diiusive Relaxation Schemes for Discrete-velocity Kinetic Equations , 2007 .

[3]  Shi Jin.  Runge-Kutta Methods for Hyperbolic Conservation Laws with Stiff Relaxation Terms , 1995 .

[4]  A. Bourlioux,et al.  High-order multi-implicit spectral deferred correction methods for problems of reactive flow , 2003 .

[5]  C. D. Levermore,et al.  Moment closure hierarchies for kinetic theories , 1996 .

[6]  Louis Tao,et al.  KINETIC THEORY FOR NEURONAL NETWORK DYNAMICS , 2006 .

[7]  A. Treves.  Mean-field analysis of neuronal spike dynamics , 1993 .

[8]  Michael L. Minion,et al.  Conservative multi-implicit spectral deferred correction methods for reacting gas dynamics , 2004 .

[9]  L. Greengard,et al.  Spectral Deferred Correction Methods for Ordinary Differential Equations , 2000 .

[10]  C. D. Levermore,et al.  Hyperbolic conservation laws with stiff relaxation terms and entropy , 1994 .

[11]  Aaditya V. Rangan,et al.  Fast numerical methods for simulating large-scale integrate-and-fire neuronal networks , 2007, Journal of Computational Neuroscience.

[12]  Bengt Fornberg,et al.  A practical guide to pseudospectral methods: Introduction , 1996 .

[13]  David Cai,et al.  Maximum-entropy closures for kinetic theories of neuronal network dynamics. , 2006, Physical review letters.

[14]  G. Russo,et al.  Implicit-explicit runge-kutta schemes and applications to hyperbolic systems with relaxation , 2005 .

[15]  Lorenzo Pareschi,et al.  Diffusive Relaxation Schemes for Multiscale Discrete-Velocity Kinetic Equations , 1998 .

[16]  Vladimir Rokhlin,et al.  On the numerical solution of two-point boundary value problems II , 1994 .

[17]  Aaditya V. Rangan.  Ph.D. thesis: Adaptive solvers for partial differential and differential-algebraic equations , 2003 .

[18]  R. Shapley,et al.  A neuronal network model of macaque primary visual cortex (V1): orientation selectivity and dynamics in the input layer 4Calpha. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

[19]  H. Grad.  On the kinetic theory of rarefied gases , 1949 .

[20]  J. Damasceno,et al.  Numerical simulation of two-dimensional flows over a circular cylinder using the immersed boundary method , 2003 .

[21]  J. Cowan,et al.  A mathematical theory of the functional dynamics of cortical and thalamic nervous tissue , 1973, Kybernetik.

[22]  Christof Koch,et al.  Biophysics of Computation: Information Processing in Single Neurons (Computational Neuroscience Series) , 1998 .

[23]  Peter A. W. Lewis,et al.  Stochastic point processes : statistical analysis, theory, and applications , 1973 .

[24]  Z. Xin,et al.  The relaxation schemes for systems of conservation laws in arbitrary space dimensions , 1995 .

[25]  Giovanni Russo,et al.  Uniformly Accurate Schemes for Hyperbolic Systems with Relaxation , 1997 .

[26]  Leslie Greengard,et al.  A Fast Adaptive Numerical Method for Stiff Two-Point Boundary Value Problems , 1997, SIAM J. Sci. Comput..

[27]  Louis Tao,et al.  An embedded network approach for scale-up of fluctuation-driven systems with preservation of spike information. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[28]  M. Shelley,et al.  An effective kinetic representation of fluctuation-driven neuronal networks with application to simple and complex cells in visual cortex. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[29]  David McLaughlin,et al.  Coarse-Grained Reduction and Analysis of a Network Model of Cortical Response: I. Drifting Grating Stimuli , 2002, Journal of Computational Neuroscience.

[30]  R. Natalini.  Convergence to equilibrium for the relaxation approximations of conservation laws , 1996 .