A numerical solver for a nonlinear Fokker-Planck equation representation of neuronal network dynamics

JOURNAL OF COMPUTATIONAL PHYSICS

卷:230

Article

Caceres, Maria J.¹  Carrillo, Jose A.^2,3  Tao, Louis⁴ 

[1] Univ Granada, Dept Matemat Aplicada, E-18071 Granada, Spain

[2] Univ Autonoma Barcelona, Dept Matemat, E-08193 Bellaterra, Spain

[3] Univ Autonoma Barcelona, ICREA, E-08193 Bellaterra, Spain

[4] Peking Univ, Coll Life Sci, Ctr Bioinformat, Natl Lab Prot Engn & Plant Genet Engn, Beijing 100871, Peoples R China

关键词: Neuronal network;    Kinetic equations;    Fokker-Planck equation;    Direct simulation Monte Carlo;    Deterministic simulations;    WENO methods;    Chang-Cooper method;   

DOI  :  10.1016/j.jcp.2010.10.027

来源: Elsevier

【 摘 要 】

To describe the collective behavior of large ensembles of neurons in neuronal network, a kinetic theory description was developed in [15,14], where a macroscopic representation of the network dynamics was directly derived from the microscopic dynamics of individual neurons, which are modeled by conductance-based, linear, integrate-and-fire point neurons. A diffusion approximation then led to a nonlinear Fokker-Planck equation for the probability density function of neuronal membrane potentials and synaptic conductances. In this work, we propose a deterministic numerical scheme for a Fokker-Planck model of an excitatory-only network. Our numerical solver allows us to obtain the time evolution of probability distribution functions, and thus, the evolution of all possible macroscopic quantities that are given by suitable moments of the probability density function. We show that this deterministic scheme is capable of capturing the bistability of stationary states observed in Monte Carlo simulations. Moreover, the transient behavior of the firing rates computed from the Fokker-Planck equation is analyzed in this bistable situation, where a bifurcation scenario, of asynchronous convergence towards stationary states, periodic synchronous solutions or damped oscillatory convergence towards stationary states, can be uncovered by increasing the strength of the excitatory coupling. Finally, the computation of moments of the probability distribution allows us to validate the applicability of a moment closure assumption used in [15] to further simplify the kinetic theory. (C) 2010 Elsevier Inc. All rights reserved.

【 授权许可】

Free