Abstract: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.

Moment neural network and an efficient numerical method for modeling irregular spiking activity

Dynamics of Moment Neuronal Networks

On Gaussian Random Neuronal Field Model: Moment Neuronal Network Approach

Moment Neuronal Networks: Stochastic Computation in Neuronal Systems

Numerical Methods for Solving Moment Equations in Kinetic Theory of Neuronal Network Dynamics

On a Gaussian Neuronal Field Model

Dynamic Modeling of Spike Count Data With Conway-Maxwell Poisson Variability

Estimation of Spike Train Statistics in Spontaneously Active Biological Neural Networks

Mean-field approximations of networks of spiking neurons with short-term synaptic plasticity

Generalized activity equations for spiking neural network dynamics

Hopf Bifurcation In Mean Field Explains Critical Avalanches In Excitation-Inhibition Balanced Neuronal Networks: A Mechanism For Multiscale Variability

Common-input models for multiple neural spike-train data

Dynamic Moment Analysis of the Extracellular Electric Field of a Biologically Realistic Spiking Neuron

Approximating Nonlinear Functions With Latent Boundaries in Low-Rank Excitatory-Inhibitory Spiking Networks

Toward stochastic neural computing

Oscillations and irregular persistent firing patterns in a homogeneous network of excitatory stochastic neurons with gap junctions in the mean-field limit

Rosetta stone for the population dynamics of spiking neuron networks

Next generation neural population models

Exact low-dimensional description for fast neural oscillations with low firing rates

Mean-field method for generic conductance-based integrate-and-fire neurons with finite timescales

Multiscale modeling via split-step methods in neural firing