Search

The voltage-conductance kinetic model for the collective behavior of neurons has been studied by scientists and mathematicians for two decades, but the rigorous analysis of its solution structure has been only partially obtained in spite of plenty of numerical evidence in various scenarios. In this work, we consider a simplified voltage-conductance model in which the velocity field in the voltage variable is in a separable form. The long time behavior of the simplified model is fully investigated, leading to the following dichotomy: either the density function converges to the global equilibrium, or the firing rate diverges as time goes to infinity. Besides, the fast conductance asymptotic limit is justified and analyzed, where the solution to the limit model either blows up in finite time or globally exists leading to time periodic solutions. An important implication of these results is that the nonseparable velocity field, or physically the leaky mechanism, is a key element for the emergence of periodic solutions in the original model. [ABSTRACT FROM AUTHOR]

The noisy leaky integrate-and-fire (NLIF) model describes the voltage configurations of neuron networks with an interacting many-particles system at a microscopic level. When simulating neuron networks of large sizes, computing a coarse-grained mean-field Fokker--Planck equation solving the voltage densities of the networks at a macroscopic level practically serves as a feasible alternative in its high efficiency and credible accuracy when the interaction within the network remains relatively low. However, the macroscopic model fails to yield valid results of the networks when simulating considerably synchronous networks with active firing events. In this paper, we propose a multiscale solver for the NLIF networks, inheriting the macroscopic solver's low cost and the microscopic solver's high reliability. For each temporal step, the multiscale solver uses the macroscopic solver when the firing rate of the simulated network is low, while it switches to the microscopic solver when the firing rate tends to blow up. Moreover, the macroscopic and microscopic solvers are integrated with a high-precision switching algorithm to ensure the accuracy of the multiscale solver. The validity of the multiscale solver is analyzed from two perspectives: first, we provide practically sufficient conditions that guarantee the mean-field approximation of the macroscopic model and present rigorous numerical analysis on simulation errors when coupling the two solvers; second, the numerical performance of the multiscale solver is validated through simulating several large neuron networks, including networks with either instantaneous or periodic input currents which prompt active firing events over some time. [ABSTRACT FROM AUTHOR]