Phase-resetting curves determine synchronization, phase locking, and clustering in networks of neural oscillators

Networks of model neurons were constructed and their activity was predicted using an iterated map based solely on the phase-resetting curves (PRCs). The predictions were quite accurate provided that the resetting to simultaneous inputs was calculated using the sum of the simultaneously active conductances, obviating the need for weak coupling assumptions. Fully synchronous activity was observed only when the slope of the PRC at a phase of zero, corresponding to spike initiation, was positive. A novel stability criterion was developed and tested for all-to-all networks of identical, identically connected neurons. When the PRC generated using N-1 simultaneously active inputs becomes too steep, the fully synchronous mode loses stability in a network of N model neurons. Therefore, the stability of synchrony can be lost by increasing the slope of this PRC either by increasing the network size or the strength of the individual synapses. Existence and stability criteria were also developed and tested for the splay mode in which neurons fire sequentially. Finally, N/M synchronous subclusters of M neurons were predicted using the intersection of parameters that supported both between-cluster splay and within-cluster synchrony. Surprisingly, the splay mode between clusters could enforce synchrony on subclusters that were incapable of synchronizing themselves. These results can be used to gain insights into the activity of networks of biological neurons whose PRCs can be measured.