Rigorous justification of the Fokker-Planck equations of neural networks based on an iteration perspective

In this work, the primary goal is to establish rigorous connection between the Fokker-Planck equation of neural networks with its microscopic model: the diffusion-jump stochastic process that captures the mean field behavior of collections of neurons in the integrate-and-fire model. The proof is based on a novel iteration scheme: with an auxiliary random variable counting the firing events, both the density function of the stochastic process and the solution of the PDE problem admit series representations, and thus the difficulty in verifying the link between the density function and the PDE solution in each sub problem has been greatly mitigated. The iteration approach provides a generic frame in integrating the probability approach with PDE techniques, with which we prove that the density function of the diffusion-jump stochastic process is indeed the classical solution of the Fokker-Planck equation with a unique flux-shift structure.

[1]  Lawrence Sirovich,et al.  Dynamics of neural populations: Stability and synchrony , 2006, Network.

[2]  Aaditya V. Rangan,et al.  DYNAMICS OF CURRENT-BASED, POISSON DRIVEN, INTEGRATE-AND-FIRE NEURONAL NETWORKS " , 2010 .

[3]  Iddo Ben-Ari,et al.  Ergodic behavior of diffusions with random jumps from the boundary , 2009 .

[4]  Nicolas Brunel,et al.  Dynamics of Sparsely Connected Networks of Excitatory and Inhibitory Spiking Neurons , 2000, Journal of Computational Neuroscience.

[5]  Sylvain Rubenthaler,et al.  First hitting times for general non-homogeneous 1d diffusion processes: density estimates in small time , 2013 .

[6]  María J Cáceres,et al.  Beyond blow-up in excitatory integrate and fire neuronal networks: Refractory period and spontaneous activity. , 2014, Journal of theoretical biology.

[7]  A. Friedman Partial Differential Equations of Parabolic Type , 1983 .

[8]  Duane Q. Nykamp,et al.  A Population Density Approach That Facilitates Large-Scale Modeling of Neural Networks: Analysis and an Application to Orientation Tuning , 2004, Journal of Computational Neuroscience.

[9]  Leszek Słomiński,et al.  Stochastic differential equations with jump reflection at time-dependent barriers , 2010 .

[10]  Duane Q. Nykamp,et al.  A Population Density Approach That Facilitates Large-Scale Modeling of Neural Networks: Extension to Slow Inhibitory Synapses , 2001, Neural Computation.

[11]  W. Feller Diffusion processes in one dimension , 1954 .

[12]  D. Talay,et al.  Mean-Field Limit of a Stochastic Particle System Smoothly Interacting Through Threshold Hitting-Times and Applications to Neural Networks with Dendritic Component , 2014, SIAM J. Math. Anal..

[13]  Wenbo V. Li,et al.  Diffusions with holding and jumping boundary , 2013 .

[14]  Ricarda Schneider,et al.  Analysis and numerical solver for excitatory-inhibitory networks with delay and refractory periods , 2017, ESAIM: Mathematical Modelling and Numerical Analysis.

[15]  David Cai,et al.  Cascade-induced synchrony in stochastically driven neuronal networks. , 2010, Physical review. E, Statistical, nonlinear, and soft matter physics.

[16]  Maria Giovanna Garroni,et al.  Green Functions for Second Order Parabolic Integro-Differential Problems , 1993 .

[17]  Wulfram Gerstner,et al.  SPIKING NEURON MODELS Single Neurons , Populations , Plasticity , 2002 .

[18]  R. Pinsky,et al.  Spectral analysis of a family of second-order elliptic operators with nonlocal boundary condition indexed by a probability measure , 2007, 0707.0612.

[19]  Lawrence Sirovich,et al.  On the Simulation of Large Populations of Neurons , 2004, Journal of Computational Neuroscience.

[20]  Jingwei Hu,et al.  A structure preserving numerical scheme for Fokker-Planck equations of neuron networks: numerical analysis and exploration , 2019, ArXiv.

[21]  F. Delarue,et al.  Global solvability of a networked integrate-and-fire model of McKean–Vlasov type , 2012, 1211.0299.

[22]  F. Delarue,et al.  Particle systems with a singular mean-field self-excitation. Application to neuronal networks , 2014, 1406.1151.

[23]  Jun Peng A note on the first passage time of diffusions with holding and jumping boundary , 2014 .

[24]  Philip Protter Stochastic differential equations with jump reflection at the boundary , 1980, Advances in Applied Probability.

[25]  S. Shreve,et al.  Stochastic differential equations , 1955, Mathematical Proceedings of the Cambridge Philosophical Society.

[26]  W. Feller THE PARABOLIC DIFFERENTIAL EQUATIONS AND THE ASSOCIATED SEMI-GROUPS OF TRANSFORMATIONS , 1952 .

[27]  P. Roux,et al.  Global-in-time solutions and qualitative properties for the NNLIF neuron model with synaptic delay , 2018, Communications in Partial Differential Equations.

[28]  Xiao-Jing Wang,et al.  Mean-Field Theory of Irregularly Spiking Neuronal Populations and Working Memory in Recurrent Cortical Networks , 2003 .

[29]  P. Goldman-Rakic,et al.  Synaptic mechanisms and network dynamics underlying spatial working memory in a cortical network model. , 2000, Cerebral cortex.

[30]  Maria E. Schonbek,et al.  Classical Solutions for a Nonlinear Fokker-Planck Equation Arising in Computational Neuroscience , 2011, 1109.1298.

[31]  Nicolas Brunel,et al.  Fast Global Oscillations in Networks of Integrate-and-Fire Neurons with Low Firing Rates , 1999, Neural Computation.

[32]  Louis Tao,et al.  A numerical solver for a nonlinear Fokker-Planck equation representation of neuronal network dynamics , 2011, J. Comput. Phys..

[33]  María J. Cáceres,et al.  Analysis of nonlinear noisy integrate & fire neuron models: blow-up and steady states , 2010, Journal of mathematical neuroscience.

[34]  M. Mattia,et al.  Population dynamics of interacting spiking neurons. , 2002, Physical review. E, Statistical, nonlinear, and soft matter physics.

[35]  Olivier D. Faugeras,et al.  Noise-Induced Behaviors in Neural Mean Field Dynamics , 2011, SIAM J. Appl. Dyn. Syst..