Systems of globally coupled phase oscillators can have robust attractors that are heteroclinic networks. Such a heteroclinic network is generated, where the phases cluster into three groups, within a specific regime of parameters when the phase oscillators are globally coupled using the function $g(\varphi) = -\sin(\varphi + \alpha) + r \sin(2\varphi + \beta)$. The resulting network switches between 30 partially synchronised states for a system of $N=5$ oscillators. Considering the states that are visited and the time spent at those states a spatio-temporal code can be generated for a given navigation around the network. We explore this phenomenon further by investigating the effect that noise has on the system, how this system can be used to generate a spatio-temporal code derived from specific inputs and how observation of a spatio-temporal code can be used to determine the inputs that were presented to the system to generate a given coding. We show that it is possible to find chaotic attractors for certain parameters and that it is possible to detail a genetic algorithm that can find the parameters required to generate a specific spatio-temporal code, even in the presence of noise. In closing we briefly explore the dynamics where $N>5$ and discuss this work in relation to winnerless competition.

Thesis advisor: Jan R. Engelbrecht / Thesis advisor: Renato E. Mirollo / Many phenomena in nature that involve ordering in time can be understood as collective behavior of coupled oscillators. One paradigm for studying a population of self-sustained oscillators is the Kuramoto model, where each oscillator is described by a phase variable, and interacts with other oscillators through trigonometric functions of phase differences. This dissertation studies $N$ identical Kuramoto oscillators in a general form \[ \dot{\theta}_{j}=A+B\cos\theta_{j}+C\sin\theta_{j}\qquad j=1,\dots,N, \] where coefficients $A$, $B$, and $C$ are symmetric functions of all oscillators $(\theta_{1},\dots,\theta_{N})$. Dynamics of this model live in group orbits of M\"obius transformations, which are low-dimensional manifolds in the full state space. When the system is a phase model (invariant under a global phase shift), trajectories in a group orbit can be identified as flows in the unit disk with an intrinsic hyperbolic metric. A simple criterion for such system to be a gradient flow is found, which leads to new classes of models that can be described by potential or Hamiltonian functions while exhibiting a large number of constants of motions. A generalization to extended phase models with non-identical couplings gives rise to richer structures of fixed points and bifurcations. When the coupling weights sum to zero, the system is simultaneously gradient and Hamiltonian. The flows mimic field lines of a two-dimensional electrostatic system consisting of equal amounts of positive and negative charges. Bifurcations on a partially synchronized subspace are discussed as well. / Thesis (PhD) — Boston College, 2017. / Submitted to: Boston College. Graduate School of Arts and Sciences. / Discipline: Physics.