State-space kinetic Ising model reveals task-dependent entropy flow in sparsely active nonequilibrium neuronal dynamics

Abstract: Neuronal ensemble activity, including coordinated and oscillatory patterns, exhibits hallmarks of nonequilibrium systems with time-asymmetric trajectories to maintain their organization. However, assessing time asymmetry from neuronal spiking activity remains challenging. The kinetic Ising model provides a framework for studying the causal, nonequilibrium dynamics in spiking recurrent neural networks. Recent theoretical advances in this model have enabled time-asymmetry estimation from large-scale steady-state data. Yet, neuronal activity often exhibits time-varying firing rates and coupling strengths, violating steady-state assumption. To overcome these limitations, we developed a state-space kinetic Ising model that accounts for non-stationary and nonequilibrium properties of neural systems. This approach incorporates a mean-field method for estimating time-varying entropy flow, a key measure for maintaining the system's organization by dissipation. Applying this method to mouse visual cortex data revealed greater variability in causal couplings during task engagement despite the reduced neuronal activity with increased sparsity. Moreover, higher-performing mice exhibited increased entropy flow in higher-firing neurons during task engagement, suggesting that stronger directed activity emerged in a fewer neurons within sparsely active populations. These findings underscore the model's utility in uncovering intricate asymmetric causal dynamics in neuronal ensembles and linking them to behavior through the thermodynamic underpinnings of neural computation.