Pairwise maximum-entropy models and their Glauber dynamics: bimodality, bistability, non-ergodicity problems, and their elimination via inhibition Vahid Rostami 1*, PierGianLuca Porta Mana 1, Moritz Helias 1,2 16 May 2016 1 Institute of Neuroscience and Medicine (INM-6) and Institute for Advanced Simulation (IAS- 6) and JARA BRAIN Institute I, Jülich Research Centre, Germany 2 Department of Physics, Faculty 1, RWTH Aachen University, Germany *
[email protected] Abstract Pairwise maximum-entropy models have been used in recent neuroscientific literature to predict the activity of neuronal populations, given only the time-averaged correlations of the neuron activities. This paper provides evidence that the pairwise model, applied to experimental recordings, predicts a bimodal distribution for the population-averaged activity, and for some population sizes the second mode peaks at high activities, with 90% of the neuron population active within time-windows of few milliseconds. This bimodality has several undesirable consequences: 1. The presence of two modes is unrealistic in view of observed neuronal activity. 2. The prediction of a high-activity mode is unrealistic on arXiv:1605.04740v1 [q-bio.NC] 16 May 2016 neurobiological grounds. 3. Boltzmann learning becomes non-ergodic, hence the pairwise model found by this method is not the maximum entropy distribution; similarly, solving the inverse problem by common variants of mean-field approximations has the same problem. 4. The Glauber dynamics associated with the model is either unrealistically bistable, or does not reflect the distribution of the pairwise model. This bimodality is first demonstrated for an experimental dataset comprising 159 neuron activities recorded from the motor cortex of macaque monkey.