Patterns of synchronization in 2D networks of inhibitory neurons

Neural firing in many inhibitory networks displays synchronous assembly or clustered firing, in which subsets of neurons fire synchronously, and these subsets may vary with different inputs to, or states of, the network. Most prior analytical and computational modeling of such networks has focused on 1D networks or 2D networks with symmetry (often circular symmetry). Here, we consider a 2D discrete network model on a general torus, where neurons are coupled to two or more nearest neighbors in three directions (horizontal, vertical, and diagonal), and allow different coupling strengths in all directions. Using phase model analysis, we establish conditions for the stability of different patterns of clustered firing behavior in the network. We then apply our results to study how variation of network connectivity and the presence of heterogeneous coupling strengths influence which patterns are stable. We confirm and supplement our results with numerical simulations of biophysical inhibitory neural network models. Our work shows that 2D networks may exhibit clustered firing behavior that cannot be predicted as a simple generalization of a 1D network, and that heterogeneity of coupling can be an important factor in determining which patterns are stable.