Phase models and clustering in networks of oscillators with delayed coupling

We consider a general model for a network of oscillators with time delayed, circulant coupling. We use the theory of weakly coupled oscillators to reduce the system of delay differential equations to a phase model where the time delay enters as a phase shift. We use the phase model to study the existence and stability of cluster solutions. Cluster solutions are phase locked solutions where the oscillators separate into groups. Oscillators within a group are synchronized while those in different groups are phase-locked. We give model independent existence and stability results for symmetric cluster solutions. We show that the presence of the time delay can lead to the coexistence of multiple stable clustering solutions. We apply our analytical results to a network of Morris Lecar neurons and compare these results with numerical continuation and simulation studies

A New Approach for Determining Phase Response Curves Reveals that Purkinje Cells Can Act as Perfect Integrators

Cerebellar Purkinje cells display complex intrinsic dynamics. They fire spontaneously, exhibit bistability, and via mutual network interactions are involved in the generation of high frequency oscillations and travelling waves of activity. To probe the dynamical properties of Purkinje cells we measured their phase response curves (PRCs). PRCs quantify the change in spike phase caused by a stimulus as a function of its temporal position within the interspike interval, and are widely used to predict neuronal responses to more complex stimulus patterns. Significant variability in the interspike interval during spontaneous firing can lead to PRCs with a low signal-to-noise ratio, requiring averaging over thousands of trials. We show using electrophysiological experiments and simulations that the PRC calculated in the traditional way by sampling the interspike interval with brief current pulses is biased. We introduce a corrected approach for calculating PRCs which eliminates this bias. Using our new approach, we show that Purkinje cell PRCs change qualitatively depending on the firing frequency of the cell. At high firing rates, Purkinje cells exhibit single-peaked, or monophasic PRCs. Surprisingly, at low firing rates, Purkinje cell PRCs are largely independent of phase, resembling PRCs of ideal non-leaky integrate-and-fire neurons. These results indicate that Purkinje cells can act as perfect integrators at low firing rates, and that the integration mode of Purkinje cells depends on their firing rate

Analysis and Modeling of Ensemble Recordings from Respiratory Pre-Motor Neurons Indicate Changes in Functional Network Architecture after Acute Hypoxia

We have combined neurophysiologic recording, statistical analysis, and computational modeling to investigate the dynamics of the respiratory network in the brainstem. Using a multielectrode array, we recorded ensembles of respiratory neurons in perfused in situ rat preparations that produce spontaneous breathing patterns, focusing on inspiratory pre-motor neurons. We compared firing rates and neuronal synchronization among these neurons before and after a brief hypoxic stimulus. We observed a significant decrease in the number of spikes after stimulation, in part due to a transient slowing of the respiratory pattern. However, the median interspike interval did not change, suggesting that the firing threshold of the neurons was not affected but rather the synaptic input was. A bootstrap analysis of synchrony between spike trains revealed that both before and after brief hypoxia, up to 45% (but typically less than 5%) of coincident spikes across neuronal pairs was not explained by chance. Most likely, this synchrony resulted from common synaptic input to the pre-motor population, an example of stochastic synchronization. After brief hypoxia most pairs were less synchronized, although some were more, suggesting that the respiratory network was transiently “rewired” after the stimulus. To investigate this hypothesis, we created a simple computational model with feed-forward divergent connections along the inspiratory pathway. Assuming that (1) the number of divergent projections was not the same for all presynaptic cells, but rather spanned a wide range and (2) that the stimulus increased inhibition at the top of the network; this model reproduced the reduction in firing rate and bootstrap-corrected synchrony subsequent to hypoxic stimulation observed in our experimental data

Clustering behaviour in networks with time delayed all-to-all coupling

Networks of coupled oscillators arise in a variety of areas. Clustering is a type of oscillatory network behavior where elements of a network segregate into groups. Elements within a group oscillate synchronously, while elements in different groups oscillate with a fixed phase difference. In this thesis, we study networks of N identical oscillators with time delayed, global circulant coupling with two approaches. We first use the theory of weakly coupled oscillators to reduce the system of delay differential equations to a phase model where the time delay enters as a phase shift. We use the phase model to determine model independent existence and stability results for symmetric cluster solutions. We show that the presence of the time delay can lead to the coexistence of multiple stable clustering solutions. We then perform stability and bifurcation analysis to the original system of delay differential equations with symmetry. We first study the existence of Hopf bifurcations induced by coupling time delay, and then use symmetric Hopf bifurcation theory to determine how these bifurcations lead to different patterns of symmetric cluster oscillations. We apply our results to two specfi c examples: a network of FitzHugh-Nagumo neurons with diffusive coupling and a network of Morris-Lecar neurons with synaptic coupling. In the case studies, we show how time delays in the coupling between neurons can give rise to switching between different stable cluster solutions, coexistence of multiple stable cluster solutions and solutions with multiple frequencies