A decision-making Fokker–Planck model in computational neuroscience

In computational neuroscience, decision-making may be explained analyzing models based on the evolution of the average firing rates of two interacting neuron populations, e.g., in bistable visual perception problems. These models typically lead to a multi-stable scenario for the concerned dynamical systems. Nevertheless, noise is an important feature of the model taking into account both the finite-size effects and the decision’s robustness. These stochastic dynamical systems can be analyzed by studying carefully their associated Fokker–Planck partial differential equation. In particular, in the Fokker–Planck setting, we analytically discuss the asymptotic behavior for large times towards a unique probability distribution, and we propose a numerical scheme capturing this convergence. These simulations are used to validate deterministic moment methods recently applied to the stochastic differential system. Further, proving the existence, positivity and uniqueness of the probability density solution for the stationary equation, as well as for the time evolving problem, we show that this stabilization does happen. Finally, we discuss the convergence of the solution for large times to the stationary state. Our approach leads to a more detailed analytical and numerical study of decision-making models applied in computational neuroscience.