Abstract:Recently developed kinetic theory for neuronal network dynamics has been demonstrated to be a powerful theoretical framework for investigating coarse-grained dynamical properties of neuronal networks. The kinetic equations are a system of (1+1)dimensional nonlinear partial differential equations (PDE) on a bounded domain with the following features: (i) the boundary conditions are nonlinear and they are themselves a functional of the present solution; (ii) the PDEs themselves are self-consistently specified by parameters which are functions of the boundary values of the solution; and (iii) the PDEs can be stiff in space and time. Numerical methods are presented here for efficiently and accurately solving these kinetic 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 fixedpoint iteration can only achieve a linear convergence for the self-consistent parameters. The practicability and efficiency of our numerical methods for solving the kinetic equations are illustrated with numerical examples. It is further demonstrated that the kinetic equations can very well capture the coarse-grained dynamical properties of integrate-and-fire neuronal networks. Corresponding author: Aaditya V. Rangan, Courant Institute of Mathematical Sciences, New York University, New York, NY 10012. Phone: 212-998-3303, Fax: 212-995-4121, email: rangan@cims.nyu.edu.

Library-based Numerical Reduction of the Hodgkin–Huxley Neuron for Network Simulation

Library-based Fast Algorithm for Simulating the Hodgkin-Huxley Neuronal Networks

Application of Efficient Numerical Methods in Solution of Ordinary Differential Equations for Modeling Electrical Activity in Cardiac Cells

Design Fast Algorithms For Hodgkin-Huxley Neuronal Networks

A Combined Offline–Online Algorithm for Hodgkin–Huxley Neural Networks

Coarse-grained Event Tree Analysis for Quantifying Hodgkin-Huxley Neuronal Network Dynamics

Numerical Methods for Solving Kinetic Equations of Neuronal Network Dynamics

Fast Algorithms for Simulation of Neuronal Dynamics Based on the Bilinear Dendritic Integration Rule

Exponential Time Differencing Algorithm for Pulse-Coupled Hodgkin-Huxley Neural Networks.

Automatic Fitting Of Neuron Parameters

Dynamics of the Exponential Integrate-and-fire Model with Slow Currents and Adaptation

On a variable step size modification of Hines' method in computational neuroscience

Channel Based Generating Function Approach to the Stochastic Hodgkin-Huxley Neuronal System

Structure-preserving numerical integrators for Hodgkin-Huxley-type systems

Compact hybrid type electronic neuron and computational model of its dynamics

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

Langevin Approach with Rescaled Noise for Stochastic Channel Dynamics in Hodgkin–Huxley Neurons

Channel-Based Langevin Approach for the Stochastic Hodgkin-Huxley Neuron

Pseudo-Lyapunov Exponents and Predictability of Hodgkin-Huxley Neuronal Network Dynamics.

Analog-Digital Simulations of Full Conductance-Based Networks of Spiking Neurons with Spike Timing Dependent Plasticity

Hypergraph-Based Numerical Spiking Neural Membrane Systems with Novel Repartition Protocols