SpringerOpen Newsletter

Receive periodic news and updates relating to SpringerOpen.

Open Access Research

Criteria for robustness of heteroclinic cycles in neural microcircuits

Peter Ashwin1*, Özkan Karabacak2 and Thomas Nowotny3

Author Affiliations

1 Mathematics Research Institute, University of Exeter, Exeter, EX4 4QF, UK

2 Faculty of Electrical and Electronics Engineering, Electronics and Communication Department, Istanbul Technical University, TR-34469, Maslak-Istanbul, Turkey

3 Centre for Computational Neuroscience and Robotics, Informatics, University of Sussex, Falmer, Brighton, BN1 9QJ, UK

For all author emails, please log on.

The Journal of Mathematical Neuroscience 2011, 1:13  doi:10.1186/2190-8567-1-13


The electronic version of this article is the complete one and can be found online at: http://www.mathematical-neuroscience.com/content/1/1/13


Received:6 September 2011
Accepted:28 November 2011
Published:28 November 2011

© 2011 Ashwin et al.; licensee Springer

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

We introduce a test for robustness of heteroclinic cycles that appear in neural microcircuits modeled as coupled dynamical cells. Robust heteroclinic cycles (RHCs) can appear as robust attractors in Lotka-Volterra-type winnerless competition (WLC) models as well as in more general coupled and/or symmetric systems. It has been previously suggested that RHCs may be relevant to a range of neural activities, from encoding and binding to spatio-temporal sequence generation.

The robustness or otherwise of such cycles depends both on the coupling structure and the internal structure of the neurons. We verify that robust heteroclinic cycles can appear in systems of three identical cells, but only if we require perturbations to preserve some invariant subspaces for the individual cells. On the other hand, heteroclinic attractors can appear robustly in systems of four or more identical cells for some symmetric coupling patterns, without restriction on the internal dynamics of the cells.

1 Introduction

For some time, it has been recognized that robust heteroclinic cycles (RHCs) can be attractors in dynamical systems [1], and that RHCs can provide useful models for the dynamics in certain biological systems. Examples include Lotka-Volterra population models [2] in ecology and game dynamics [3]. Similar dynamics has been used to describe various neuronal microcircuits, in particular winnerless competition (WLC) dynamics [4] has been the subject of intense recent study. For example, [5] find conditions on the connectivity scheme of the generalised Lotka-Volterra model to guarantee the existence and structural robustness of a heteroclinic cycle in the system, [6] consider generalised “heteroclinic channels”, [7] use them as a model for sequential memory and [8] suggest that they may be used to describe binding problems. One question raised by these studies is whether Lotka-Volterra type dynamics is necessary to give robust heteroclinic cycles as attractors and how these cycles relate to those found in other models [9,10]. The purpose of this paper is to show that attracting heteroclinic cycles may be robust for a variety of reasons and appear in a variety of dynamical systems that model neural microcircuits. In doing so, we give a practical test for robustness of heteroclinic cycles within any particular context and demonstrate it in practice for several examples.

This paper was motivated by a recent paper on three synaptically coupled Hodgkin-Huxley type neurons in a ring that reported robust winnerless competition between neurons [11] without an explicit Lotka-Volterra type structure. This manifested as a cyclic progression between states where only one neuron is active (spiking) for a period of time. During this activity, the currently active neuron inhibits the activity of the next neuron in the ring while the third neuron recovers from previous inhibition.

One of the main observations of this paper is that the coupling structure and symmetries in this system are not sufficient to guarantee robustness of the heteroclinic behaviour observed in [11], but robustness can be demonstrated if we consider constraints in the system. For this case it is natural to investigate the invariance of a set of affine subspaces of the system’s phase space related to the type of synaptic coupling considered. More generally, we discuss cases of heteroclinic attractors that are robust, based purely on the coupling structure and the assumption that the cells are identical.

The paper is organized as follows: In Section 2 we consider the general problem of robustness of a heteroclinic cycle. We investigate a class of dynamical systems that have affine invariant subspaces and give a necessary and sufficient condition on the dimensionality of the invariant affine subspaces for the robustness of HCs in this class of systems. We translate these conditions into appropriate conditions for coupled systems. Section 3.1 reviews a simple example of winnerless competition and demonstrates robustness for Lotka-Volterra systems, while Section 3.2 discusses the three-cell problem of Nowotny et al.[11]. We demonstrate how the general results from Section 2 can be applied to show that the observed HC in the system (i) is not robust with respect to perturbations that only preserve its <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M1">View MathML</a> symmetry, but (ii) is robust with respect to perturbations that respect a specific set of invariant affine subspaces. Section 3.3 illustrates an example of a four-cell network of Hodgkin-Huxley type neurons where the coupling structure alone is sufficient for the robustness of HCs. We finish with a brief discussion in Section 4.

2 Robustness of heteroclinic cycles

Suppose we have a dynamical system given by a system of first order differential equations

<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M2">View MathML</a>

(1)

where <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M3">View MathML</a> and <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M4">View MathML</a>, the set of <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M5">View MathML</a> vector fields on <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M6">View MathML</a> with bounded global attractors.1 We say an invariant set Σ is a heteroclinic cycle (HC) if it consists of a union of hyperbolic equilibria <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M8">View MathML</a> and connecting orbits <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M9">View MathML</a>.2 We say that a heteroclinic cycle Σ is robust to perturbations in<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M10">View MathML</a> if <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M11">View MathML</a> and there is a <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M5">View MathML</a>-neighbourhood of f such that all <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M13">View MathML</a> within this neighbourhood have a heteroclinic cycle that is close to Σ.

Let us suppose that <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M4">View MathML</a> has a HC Σ between equilibria <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M15">View MathML</a>. As the connection <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M16">View MathML</a> is contained within <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M17','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M17">View MathML</a>, this implies that <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M18','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M18">View MathML</a>. In order for the connection from <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M15">View MathML</a> to <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M20','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M20">View MathML</a> to be robust with respect to arbitrary <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M5">View MathML</a> perturbations it is necessary that the intersection is transverse [12], meaning that

<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M22','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M22">View MathML</a>

(2)

Using the fact that <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M23','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M23">View MathML</a> for any hyperbolic equilibrium and adding these for all equilibria along the cycle, we find that

<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M24','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M24">View MathML</a>

(3)

This implies that it is not possible for Equation 2 to be satisfied for all connections. Hence our first statement is the following (which can be thought of a special case of the Kupka-Smale Theorem [12], see also [13]).

Proposition 1A heteroclinic cycle between<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M25','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M25">View MathML</a>hyperbolic equilibria is never robust to general<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M5">View MathML</a>perturbations in<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M27','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M27">View MathML</a>.

The heteroclinic cycle may however be robust to a constrained set of perturbations. We explore this in the following sections.

2.1 Conditions for robustness of heteroclinic cycles with constraints

A subset <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M28','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M28">View MathML</a> is an affine subspace if it can be written as <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M29','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M29">View MathML</a> for some real-valued <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M30','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M30">View MathML</a> matrix A and vector <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M31','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M31">View MathML</a> (this is a linear subspace if b can be chosen to be zero). For a given phase space <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M6">View MathML</a>, suppose that we have a (finite) set of non-empty affine subspaces

<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M33','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M33">View MathML</a>

(4)

that are closed under intersection; i.e. the intersection <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M34','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M34">View MathML</a> of any two subspaces <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M35','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M35">View MathML</a> is an element of <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M36','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M36">View MathML</a> unless it is empty. We include <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M37','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M37">View MathML</a>, which is trivially invariant, so <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M36','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M36">View MathML</a> is always non-empty. For a given <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M36','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M36">View MathML</a>, we define the set of vector fields (in<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M27','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M27">View MathML</a>) respecting<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M36','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M36">View MathML</a> to be

<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M42','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M42">View MathML</a>

(5)

and call the subspaces in <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M36','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M36">View MathML</a>invariant subspaces in the phase space of the dynamical systems described by <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M44','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M44">View MathML</a>.

A set of invariant affine subspaces <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M36','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M36">View MathML</a> may arise from a variety of modelling assumptions; for example,

• If f is a Lotka-Volterra type population model that leaves some subspaces corresponding to the absence of one or more “species” invariant then <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M46','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M46">View MathML</a> where <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M36','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M36">View MathML</a> is the set of the invariant subspaces forced by the absence of these species.

• If f is symmetric (equivariant) for some group action G and <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M36','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M36">View MathML</a> is the set of fixed point subspaces of G then <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M46','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M46">View MathML</a> because fixed point subspaces are invariant under the dynamics of equivariant systems [14], Theorem 1.17]. Note that for an orthogonal group action, the fixed point subspaces are linear subspaces. It is known that symmetries impose further constraints on the dynamics such as repeated eigenvalues or missing terms in Taylor expansions [14] but we focus here only on the invariant subspaces.

• If f is a realization of a particular coupled cell system with a given coupling structure then <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M46','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M46">View MathML</a> where <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M36','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M36">View MathML</a> corresponds to the set of possible cluster states (also called synchrony subspaces or polydiagonals in the literature [15-17]).

Note that <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M52','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M52">View MathML</a> inherits a subset topology from <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M27','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M27">View MathML</a>; for a discussion of homoclinic and heteroclinic phenomena in general and their associated bifurcations in particular, we refer to the review [13].

Suppose that for a vector field <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M46','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M46">View MathML</a> we have a heteroclinic cycle Σ between hyperbolic equilibria <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M55','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M55">View MathML</a> (<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M56','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M56">View MathML</a>) with connections <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M16">View MathML</a> from <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M15">View MathML</a> to <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M20','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M20">View MathML</a>. We define

<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M60','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M60">View MathML</a>

(6)

i.e. the smallest subspace in <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M36','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M36">View MathML</a> that contains <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M16">View MathML</a>. The invariant set <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M63','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M63">View MathML</a> is clearly well defined because <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M36','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M36">View MathML</a> is closed under intersections. We define the connection scheme of the heteroclinic cycle to be the sequence

<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M65','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M65">View MathML</a>

(7)

The following theorem gives necessary and sufficient conditions for such a heteroclinic cycle to be robust to perturbations in <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M52','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M52">View MathML</a>, depending on its connection scheme (we will require robustness to preserve the connection scheme). More precisely it depends on the following equation being satisfied:

<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M67','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M67">View MathML</a>

(8)

for each i. Note that there is a slight complication for the sufficient condition - it may be necessary to perturb the system slightly within <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M52','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M52">View MathML</a> to unfold the intersection to general position and remove a tangency between <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M69','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M69">View MathML</a> and <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M70','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M70">View MathML</a>. This complication has the benefit that it allows us to make statements about particular connections without needing to verify that the intersection of manifolds is transverse.

Theorem 1Let Σ be a heteroclinic cycle for<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M46','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M46">View MathML</a>between hyperbolic equilibria<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M72','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M72">View MathML</a>with connection scheme Equation 7.

1. If the cycle Σ is robust to perturbations in<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M52','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M52">View MathML</a>then Equation 8is satisfied for<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M56','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M56">View MathML</a>.

2. Conversely, if Equation 8is satisfied for<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M56','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M56">View MathML</a>then there is a nearby<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M76','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M76">View MathML</a> (with<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M77','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M77">View MathML</a>arbitrarily close tof) such that Σ is a heteroclinic cycle for<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M77','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M77">View MathML</a>that is robust to perturbations in<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M52','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M52">View MathML</a>.

Proof We will abbreviate <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M80','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M80">View MathML</a>. Because <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M16">View MathML</a> is a connection from <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M15">View MathML</a> to <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M20','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M20">View MathML</a>, there is a non-trivial intersection of <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M17','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M17">View MathML</a> within <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M85','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M85">View MathML</a>. As <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M85','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M85">View MathML</a> is the smallest invariant subspace containing <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M16">View MathML</a>, typical points <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M88','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M88">View MathML</a> will have a neighbourhood in <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M85','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M85">View MathML</a> that contain no points in any other <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M90','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M90">View MathML</a>. In a neighbourhood of this y, perturbations of f in <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M52','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M52">View MathML</a> have no restriction other than they should leave <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M85','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M85">View MathML</a> invariant.

The stability of the intersection of the unstable and stable manifolds depends on the dimension of the unstable manifolds (also called the Morse index[13]) for these equilibria for the vector field restricted to <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M85','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M85">View MathML</a>. Pick any codimension one section <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M94','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M94">View MathML</a> transverse to the connection at y. We have

<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M95','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M95">View MathML</a>

(9)

and within P, the invariant manifolds have dimensions

<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M96','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M96">View MathML</a>

(10)

The intersection of these invariant manifolds may not be transverse within P, but it will be for a dense set of nearby vector fields. In particular, if

<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M97','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M97">View MathML</a>

(11)

then there will be an open dense set of perturbations of f that remove the intersection, giving lack of robustness of <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M16">View MathML</a> and hence we obtain a proof for case 1. On the other hand, if Equation 11 is not satisfied, we can choose a vector field <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M77','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M77">View MathML</a> that is identical to f except on a small neighbourhood ofy- there it is chosen to preserve the connection but to perturb the manifolds so that the intersection is transverse. Transversality of the intersection then implies robustness of the connection and hence we obtain a proof for case 2. □

Note that caution is necessary in interpreting this result for a number of reasons:

1. Just because a given heteroclinic connection is not robust due to this result does not necessarily imply that there is no robust connection from <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M15">View MathML</a> to <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M20','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M20">View MathML</a> at all. Indeed, it may be [18] that there are several connections from <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M15">View MathML</a> to <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M20','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M20">View MathML</a> and that perturbations will break some but not all of them. In this sense, it may be that at the same time, one heteroclinic cycle is not robust, but another heteroclinic cycle between the same equilibria is robust.

2. We consider robustness to perturbations that preserve the connection scheme - there are situations where a typical perturbation may break a connection but preserve a nearby connection in a larger invariant subspace. This situation will typically only occur in exceptional cases.

3. The structure of general robust heteroclinic cycles may be very complex even if we only examine cases forced by symmetries - they easily form networks with multiple cycles. There may be multiple or even a continuum of connections between two equilibria, and they may be embedded in more general “heteroclinic networks” where there may be connections to “heteroclinic subcycles” [16,19,20].

4. Theorem 1 does not consider any dynamical stability (attraction) properties of the heteroclinic cycles.

5. In what follows, we slightly abuse notation by saying that a heteroclinic cycle is robust if the cycle for an arbitrarily small perturbation of the vector field is robust.

If <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M69','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M69">View MathML</a> is not contained in <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M70','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M70">View MathML</a> then the heteroclinic cycle Σ cannot be asymptotically stable. We say that an invariant set Σ is a regular heteroclinic cycle if it consists of a union of equilibria and a set of connecting orbits <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M106','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M106">View MathML</a> with <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M107','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M107">View MathML</a>. The following result is stated in [13] for the case of symmetric systems.

Theorem 2Suppose that Σ is a regular heteroclinic cycle for<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M46','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M46">View MathML</a>between hyperbolic equilibria<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M72','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M72">View MathML</a>. Suppose that<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M20','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M20">View MathML</a>is a sink for the dynamics reduced to<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M63','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M63">View MathML</a>, i.e.

<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M112','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M112">View MathML</a>

(12)

for alli. Then the heteroclinic cycle is robust to perturbations within<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M52','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M52">View MathML</a>.

Proof Suppose that <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M114','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M114">View MathML</a>. Since <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M69','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M69">View MathML</a> is contained in <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M70','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M70">View MathML</a> by regularity of the HC, and because <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M117','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M117">View MathML</a>, we find <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M118','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M118">View MathML</a>. Hence, Equation 8 follows and we could apply Theorem 1 case 2. In fact this is a simpler case in that because <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M119','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M119">View MathML</a> the intersection must already be transverse - one does not need to consider any perturbations to force transversality of the intersection. □

2.2 Cluster states for coupled systems

RHCs may appear in coupled systems due to a variety of constraints from the coupling structure - these are associated with cluster states (also called synchrony subspaces [15] or polydiagonals for the network [21]). Consider a network of N systems each with phase space <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M120','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M120">View MathML</a> and coupled to each other to give a set of differential equations on <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M6">View MathML</a>, with <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M122','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M122">View MathML</a>, of the form

<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M123','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M123">View MathML</a>

(13)

for <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M124','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M124">View MathML</a><a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M125','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M125">View MathML</a>. We write <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M126','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M126">View MathML</a> with <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M127','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M127">View MathML</a>. We define a cluster state for a class of ODEs to be a partition <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M128','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M128">View MathML</a> of <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M129','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M129">View MathML</a> such that the linear subspace

<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M130','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M130">View MathML</a>

is dynamically invariant for all ODEs in that class. For a given symmetry or coupling structure, we identify a list of possible cluster states and use these to test for robustness of any given heteroclinic cycle using Theorem 1.

We remark that the simplest (and indeed only, up to relabelling) coupling structure for a network of three identical cells found by [15] to admit heteroclinic cycles can be represented as a system of the form

<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M131','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M131">View MathML</a>

(14)

For an open set of choices of <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M132','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M132">View MathML</a>, the heteroclinic cycle involves two saddles within the subspace <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M133','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M133">View MathML</a> and connections that are contained within <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M134','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M134">View MathML</a> in one direction and <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M135','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M135">View MathML</a> in the other. This represents a system of three identical units coupled in a specific way, where each unit has two different input types; we refer to [15] for details. It can be quite difficult to find a suitable function f that gives a robust heteroclinic cycle in this case. Nevertheless, once one has found a heteroclinic cycle, it can be shown to be robust using Theorem 1 (case 2).

Other examples of robust heteroclinic cycles between equilibria for systems of coupled phase oscillators are given in [22,23]. For such systems the final state equations are obtained by reducing the dynamics to phase difference variables. In this case, each equilibrium represents the oscillatory motion of oscillators with some fixed phase difference.

2.3 Robust heteroclinic cycles between periodic orbits

In cases where a phase difference reduction is not possible, one may need to study heteroclinic cycles between periodic orbits in order to explain heteroclinic behaviour. Unlike heteroclinic cycles between equilibria, heteroclinic cycles between periodic orbits can be robust under general perturbations since for a hyperbolic periodic orbit p, <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M136','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M136">View MathML</a>. Hence, the condition Equation 2 can be satisfied. For instance, consider a system on <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M137','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M137">View MathML</a> with two hyperbolic periodic orbits p and q for which the stable and unstable manifolds <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M138','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M138">View MathML</a>, <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M139','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M139">View MathML</a>, <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M140','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M140">View MathML</a>, and <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M141','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M141">View MathML</a> are two-dimensional. In this case, <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M139','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M139">View MathML</a> and <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M140','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M140">View MathML</a> (and similarly, <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M141','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M141">View MathML</a> and <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M138','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M138">View MathML</a>) intersect transversely, and therefore, a heteroclinic cycle between p and q can exist robustly. However, for this heteroclinic cycle only one orbit connects p to q, whereas infinitely many orbits which are backward asymptotic to p move away from the heteroclinic cycle. As a result, such a robust heteroclinic cycle cannot be asymptotically stable.

To overcome this difficulty we assume that the connections of a heteroclinic cycle between periodic orbits consist of unstable manifolds of periodic orbits and these are contained in the stable manifold of the next periodic orbit. Namely, we say an invariant set Σ is a heteroclinic cycle that contains all unstable manifolds if it consists of a union of periodic orbits and/or equilibria <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M8">View MathML</a> and a set of connecting manifolds <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M147','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M147">View MathML</a> with <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M107','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M107">View MathML</a>.

Theorem 3Suppose that Σ is a heteroclinic cycle that contains all unstable manifolds for<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M46','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M46">View MathML</a>between hyperbolic equilibria or periodic orbits<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M72','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M72">View MathML</a>. If there exists a finite sequence<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M151','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M151">View MathML</a>of elements in<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M36','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M36">View MathML</a>such that<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M153','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M153">View MathML</a>and

<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M154','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M154">View MathML</a>

(15)

(in other words, <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M20','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M20">View MathML</a>is a sink for the dynamics reduced to<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M63','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M63">View MathML</a>) for all<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M56','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M56">View MathML</a>then Σ is robust to perturbations within<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M52','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M52">View MathML</a>.

Proof Consider a unique orbit <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M159','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M159">View MathML</a>. Since <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M70','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M70">View MathML</a>contains a neighbourhood of<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M20','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M20">View MathML</a>in<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M63','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M63">View MathML</a>, <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M16">View MathML</a> is robust by the same reasoning as in the proof of Theorem 2. This implies that the manifold of connections <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M164','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M164">View MathML</a> is robust for all i. □

Note that a heteroclinic cycle may contain all unstable manifolds but not be attracting even in a very weak sense (essentially asymptotically stable [24]). Conversely, a heteroclinic cycle may not contain all unstable manifolds but may be essentially asymptotically stable.

3 Robust heteroclinic behaviour in neural models

We discuss three examples of cases where robust heteroclinic behaviour can be found in simple neural microcircuits.

3.1 Winnerless competition in Lotka-Volterra rate models

The review [25] includes a discussion of winnerless competition and related phenomena. This has focused on the dynamics of Lotka-Volterra type models for firing rates, justified by an approximation of Fukai and Tanaka [26]. In their most general form, these are written as

<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M165','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M165">View MathML</a>

(16)

where <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M15">View MathML</a> for <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M125','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M125">View MathML</a> is the firing rate of some neuron (or neural assembly) and <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M168','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M168">View MathML</a> is a nonlinear function that represents both the intrinsic firing and that due to interaction with the other cells in the network. These systems have a very rich set of invariant subspaces because of the invariance of all subspaces where <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M169','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M169">View MathML</a>. More precisely, given any subset <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M170','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M170">View MathML</a> there is an invariant subspace corresponding to

<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M171','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M171">View MathML</a>

for example <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M172','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M172">View MathML</a>. This gives a total of <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M173','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M173">View MathML</a> invariant subspaces for the dynamics of Equation 16. Using these one can find a connection scheme involving these <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M174','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M174">View MathML</a> such that Theorem 1 can be applied to check robustness of a specific heteroclinic cycle to perturbations that preserve the form Equation 16. For example, the following rate model for the pyloric CPG of the lobster stomatogastric ganglion is discussed in [25]:

<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M175','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M175">View MathML</a>

(17)

with <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M164','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M164">View MathML</a> representing the stimulus and <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M177','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M177">View MathML</a> the firing rate of the ith neuron. In the absence of stimulus <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M178','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M178">View MathML</a> this exhibits heteroclinic cycles for <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M179','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M179">View MathML</a>

<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M180','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M180">View MathML</a>

and <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M181','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M181">View MathML</a>. These heteroclinic cycles connect three equilibria of Equation 17, namely <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M182','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M182">View MathML</a>. Calculating linearizations of Equation 17 at these equilibria one can show that, for the equilibrium <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M15">View MathML</a>, three linearly independent eigenvectors are contained in <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M184','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M184">View MathML</a><a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M90','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M90">View MathML</a> and <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M186','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M186">View MathML</a> with eigenvalues <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M187','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M187">View MathML</a><a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M188','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M188">View MathML</a><a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M189','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M189">View MathML</a>, respectively, where <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M190','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M190">View MathML</a> are different indices. Hence, when ρ is chosen as above, it follows that

<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M191','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M191">View MathML</a>

Finally, from Theorem 1 case 2, we can conclude that the heteroclinic cycle between saddle equilibria <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M192','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M192">View MathML</a><a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M193','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M193">View MathML</a> and <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M194','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M194">View MathML</a> is robust for the robust connection scheme

<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M195','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M195">View MathML</a>

3.2 Robustness of a heteroclinic cycle in a rate model with synaptic coupling

We now turn to the robustness of heteroclinic cycles in a specific model of <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M179','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M179">View MathML</a> coupled neurons derived from a Hodgkin-Huxley type model with synaptic coupling [11], a case where we do not have the Lotka-Volterra structure Equation 16. If the synaptic time scales are slow compared to the time scale of the individual spikes, then the full conductance based model can be reduced systematically to an approximate rate model [11], equations (13,14)]:

<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M197','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M197">View MathML</a>

(18)

where time variable t is in ms. The unitless dynamical variables <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M198','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M198">View MathML</a> represent the fraction of presynaptically released and <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M16">View MathML</a> the fraction of postsynaptically bound neurotransmitter for the ith neuron (<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M125','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M125">View MathML</a>), and

<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M201','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M201">View MathML</a>

characterises the rate response of the neurons to input current. We have introduced a smoothing factor <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M202','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M202">View MathML</a>, with small <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M203','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M203">View MathML</a> to ensure that F is <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M5">View MathML</a> without affecting the overall structure of the model. We use parameters as in Table 1 and couple <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M179','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M179">View MathML</a> cells in a ring using different coupling strength in each direction:

<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M206','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M206">View MathML</a>

(19)

A typical timeseries showing an attracting heteroclinic cycle for this system is shown in Figure 1.

thumbnailFig. 1. A trajectory approaching a heteroclinic cycle for the rate model (18) with <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M207','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M207">View MathML</a> and parameters as in Table 1. Observe that the trajectory cycles between the neighbourhoods of three saddle equilibria where two of the three <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M16">View MathML</a> are close to saturated at <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M209','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M209">View MathML</a>. The switching dynamics between equilibria continues to slow down indicating that the trajectory is approaching the actual heteroclinic orbit. The bottom panel illustrates this further by showing that the <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M16">View MathML</a> continue to approach the equilibria over the whole duration of the simulation. The simulation is performed using a numerical scheme that carefully resolves the behaviour near the equilibria.

Table 1. Numerical values of parameters used for simulation of the rate model (18); see [11] for a discussion of the derivation of the model and the meaning of the parameters.

The heteroclinic cycle <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M215','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M215">View MathML</a> connects the saddle equilibria <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M216','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M216">View MathML</a>, <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M217','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M217">View MathML</a>, <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M218','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M218">View MathML</a> listed in Table 2, all of which have one-dimensional unstable manifolds (unstable eigenvalue 0.0062) and five-dimensional stable manifolds (stable eigenvalues −0.0066, −0.01, −0.02, −0.02, −0.02). Adjusting any of these parameters appears to preserve the heteroclinic attractor. This raises the question whether the symmetry in the system is necessary or sufficient to ensure robustness of a heteroclinic cycle. We investigate the robustness of this cycle in the light of Theorem 1 to show that in fact the presence of this symmetry is neither necessary nor sufficient to ensure robustness.

Table 2. Equilibria of (18) involved in the heteroclinic cycle for <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M207','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M207">View MathML</a> and parameters as in Table 1.

Theorem 4There are heteroclinic cycles in the system Equation 18with parameters in Table 1. These cycles:

are not robust to perturbations that preserve the<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M239','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M239">View MathML</a>symmetry of cyclic permutation of the cells.

are robust to perturbations that preserve the affine subspaces associated with<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M240','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M240">View MathML</a>.

Proof(We do not rigorously prove that the heteroclinic cycles exists; this should in principle be possible via rigorous methods with an error bounded integrator - see for example[27].) To show the first part, note that the only invariant subspaces in <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M241','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M241">View MathML</a> permitted by <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M239','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M239">View MathML</a> permutation symmetries are

<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M243','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M243">View MathML</a>

Since <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M244','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M244">View MathML</a> in all cases, Theorem 1 case 1 implies that typical symmetry-preserving perturbations of the system destroy the heteroclinic cycle.

To see the second part, let us consider the set of vector fields on <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M245','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M245">View MathML</a> that preserve the property that <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M246','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M246">View MathML</a> is invariant: this means that we assume that the following set of subspaces are invariant:

<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M247','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M247">View MathML</a>

(20)

Examining the equilibria in Table 2 we note that the <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M15">View MathML</a> are connected in the following connection scheme:

<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M249','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M249">View MathML</a>

For the particular choice of parameters in Table 1, there is a heteroclinic cycle between three equilibria <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M250','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M250">View MathML</a>, <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M251','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M251">View MathML</a>, <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M252','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M252">View MathML</a>. These equilibria have unstable/stable manifolds that intersect to form a heteroclinic loop and satisfy

<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M253','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M253">View MathML</a>

Hence the criteria of Theorem 1 (case 2) are satisfied and the heteroclinic cycle is robust with respect to <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M5">View MathML</a>-perturbations that preserve the subspaces Equation 20. □

3.3 Robustness of heteroclinic cycles for a delay-coupled Hodgkin-Huxley type model

One might suspect that Theorem 4 can be generalized to show that internal constraints might be needed to give robustness of HCs for larger numbers of cells, but this is not the case as long as the cells are assumed identical. For example [28-30] find robust cycles in systems of four or more identical, globally coupled phase oscillators with no further constraints.

To illustrate this, we give an example of a robust heteroclinic attractor for a model system of four synaptically coupled neurons. We use a modification of Rinzel’s neuron model [31] presented by Rubin [32] with synaptic coupling [32]. Due to the global coupling of the system, the invariant subspaces are all nontrivial cluster states.

Consider N all-to-all synaptically coupled neurons with delay coupling (using units of mV for voltages, ms for time, mS/cm2 for conductances, μA/cm2 for currents, and μF/cm2 for capacitance):

<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M255','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M255">View MathML</a>

(21)

for <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M256','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M256">View MathML</a>, where

<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M257','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M257">View MathML</a>

We consider the parameters <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M258','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M258">View MathML</a>, <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M259','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M259">View MathML</a>, <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M260','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M260">View MathML</a>, <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M261','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M261">View MathML</a>, <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M262','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M262">View MathML</a>, <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M263','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M263">View MathML</a>, <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M264','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M264">View MathML</a>, <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M265','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M265">View MathML</a> and synaptic coupling parameters <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M266','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M266">View MathML</a>, <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M267','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M267">View MathML</a>, <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M268','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M268">View MathML</a>.

The dynamics of this model is oscillatory for these parameter values. For the purpose of visualizing the dynamics, we define an approximate phase using a projection of the oscillation signal onto the <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M269','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M269">View MathML</a> plane (see Figure 2). A reference point <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M270','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M270">View MathML</a> is chosen and the approximate phase is given as

<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M271','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M271">View MathML</a>

thumbnailFig. 2. Oscillation of an uncoupled Hodgkin-Huxley type neuron (21) in the <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M272','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M272">View MathML</a>-plane. The phase variable is estimated with respect to the reference point <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M270','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M270">View MathML</a>.

For two different neurons we use the synchronization index

<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M274','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M274">View MathML</a>

as a measure of their phase synchronization. The neurons i and j are completely phase synchronized when <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M275','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M275">View MathML</a>.

For <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M276','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M276">View MathML</a>, a heteroclinic cycle exists as shown in Figure 3. This is a heteroclinic cycle between two saddle periodic orbits with the same clustering, that is <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M277','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M277">View MathML</a>. These saddle periodic orbits <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M216','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M216">View MathML</a> and <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M217','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M217">View MathML</a> form a heteroclinic cycle with a scheme

<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M280','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M280">View MathML</a>

where

<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M281','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M281">View MathML</a>

and <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M282','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M282">View MathML</a>. Theorem 3 implies that as long as (a) the periodic orbits <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M216','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M216">View MathML</a> (resp. <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M217','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M217">View MathML</a>) are hyperbolic, and (b) they are sinks when considered within the subspaces <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M285','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M285">View MathML</a> (resp. <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M286','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M286">View MathML</a>) then the connection is robust. Condition (a) is generically satisfied. We have checked (b) using simulations by choosing different initial conditions for the dynamics reduced to <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M285','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M285">View MathML</a> and <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M286','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M286">View MathML</a>.

thumbnailFig. 3. A solution of (21) for <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M276','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M276">View MathML</a> and <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M290','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M290">View MathML</a> approaching to a heteroclinic cycle between two cluster states with the same clustering <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M277','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M277">View MathML</a> but a different effective phase difference. In the upper graph, synchronization indices for neuron pairs <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M292','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M292">View MathML</a> and <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M293','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M293">View MathML</a> are plotted (after block averaging of size 1000), whereas in the lower graph a shorter time series of membrane voltages shows transitions between different synchronized clusters.

Coupled phase oscillators are used as simplified models for weakly coupled limit cycle oscillators, and one can find one-to-one correspondence between solutions if the coupling is weak enough [33]. In particular, the heteroclinic cycle depicted in Figure 3 corresponds in clustering type to a heteroclinic cycle found in [22] (see Figure 4) for a system of four globally coupled phase oscillators. The two saddle periodic orbits mentioned above correspond to the two saddle equilibria in Figure 4.

thumbnailFig. 4. A robust heteroclinic cycle for four all-to-all coupled phase oscillator system analogous to the cycle found in Figure 3 for the Hodgkin-Huxley type system. The heteroclinic cycle consists of two saddle equilibria <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M216','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M216">View MathML</a> and <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M217','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M217">View MathML</a> and connections <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M221','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M221">View MathML</a> and <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M223','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M223">View MathML</a> on invariant subspaces. The invariant subspaces are embedded in a cube that represents a unit cell for the torus of phase difference space- in this representation the vertices represent in-phase solutions where all oscillators are synchronized. (Adapted from [22].)

For <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M298','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M298">View MathML</a>, more complex heteroclinic cycles can appear as seen in Figure 5. This is a heteroclinic cycle that connects different cluster states of type<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M299','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M299">View MathML</a>. Note that the transition times between clusters are fixed but the duration of stay at each cluster gets longer and longer - a feature of attracting heteroclinic cycles. In the case of noisy systems, the dynamics switches from one cluster to another randomly around a graph of connections between symmetric cluster states [10].

thumbnailFig. 5. A trajectory of the system of <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M298','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M298">View MathML</a> neural oscillators (21) for <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M301','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M301">View MathML</a> approaching to a heteroclinic cycle between two clusters. In the first five graphs, five synchronization indices are plotted (after block averaging of size 1000) for different synchronized pairs, whereas in the last graph transitions between synchronized clusters are illustrated.

For globally coupled networks of <a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M302','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M302">View MathML</a> phase oscillators, robust heteroclinic cycles between cluster states have been found in [22,28,30]. Such robust heteroclinic cycles of coupled phase oscillators involve robust connections between saddle-type cluster states, where the robustness of the connections relies on them being contained within another nontrivial cluster state that corresponds to partially breaking the clusters and reforming them in a different way.

4 Discussion

In this paper we have introduced a testable criterion for robustness for a given cycle of heteroclinic connections within constrained settings - this test involves finding the connection scheme and then applying Theorem 1. We have attempted to clarify the similarity between winnerless competition dynamics in Lotka-Volterra systems as a special case of robust heteroclinic dynamics that respect some set of invariant subspaces in a connection scheme.

Winnerless competition has previously been used to describe the competition of modes where at each mode a different neuron or neuron ensemble is active and other neurons or neuron ensembles remain inactive [8,34]. This type of competition relies on a stable robust heteroclinic cycle where robustness is due to the constraints on the individual dynamics of neurons. However, models where constraints are only on the coupling structure can admit a general phenomenon, namely robust heteroclinic cycles between cluster states. The model analyzed in Section 3.3 is an example with RHCs between cluster states. This dynamics relies on a stable robust heteroclinic cycle where robustness is due to the invariant subspaces forced by the coupling structure. In this case, the heteroclinic cycle connects saddle equilibria or saddle periodic orbits that represent different cluster states.

We have not discussed the robustness of attraction properties of RHCs - mere existence of a RHC is not enough to guarantee that it will be an attractor, but we mention that as attraction properties are determined by open conditions on eigenvalues of the saddles (e.g. [1,24,35]), continuity of variation of the eigenvalues will guarantee that attractivity is also a robust property.

For larger numbers of cells in symmetric or asymmetric arrays there may be very many such invariant subspaces, giving a wide range of possible robust heteroclinic cycles. Some of these are constructed in [15] for small numbers of coupled cells, but up to now there does not seem to be an easy way to explore which cycles are possible and which are not within any particular system. On the other hand, verifying that a particular heteroclinic cycle is, or is not, robust is a more tractable question that we address here. Note that which cycles exist may depend not just on having a valid connection scheme for some constrained set of vector fields, but also on the constraints not preventing the existence of the appropriate saddles or connections between them.

Finally, we remark that there is evidence of metastable states in neural systems (e.g. [36-38]) that are supportive of the presence of approximate robust heteroclinic cycles. There are also suggestions that heteroclinic cycles may facilitate certain computational properties of neural systems - see for example [7,39,40].

Competing interests

The authors declare that they have no competing interests.

Footnotes

1We work within the class of continuously differentiable vector fields (<a onClick="popup('http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.mathematical-neuroscience.com/content/1/1/13/mathml/M5">View MathML</a>) to ensure, by the Hartman Grobman theorem [12], that hyperbolic equilibria are robust - this is a minimal requirement to discuss robustness of heteroclinic cycles.

2We take the subscripts modulo p.

References

  1. Krupa, M: Robust heteroclinic cycles. J. Nonlinear Sci.. 7(2), 129–176 (1997). Publisher Full Text OpenURL

  2. May, RM, Leonard, WJ: Nonlinear aspects of competition between three species. SIAM J. Appl. Math.. 29(2), 243–253 (1975). Publisher Full Text OpenURL

  3. Hofbauer, J, Sigmund, K: Evolutionary Games and Population Dynamics, Cambridge University Press, Cambridge (1998)

  4. Rabinovich, MI, Volkovskii, A, Lecanda, P, Huerta, R, Abarbanel, HDI, Laurent, G: Dynamical encoding by networks of competing neuron groups: Winnerless competition. Phys. Rev. Lett.. 87(6), (2001)

  5. Afraimovich, VS, Rabinovich, MI, Varona, P: Heteroclinic contours in neural ensembles and the winnerless competition principle. Int. J. Bifurcat. Chaos. 14(4), 1195–1208 (2004). Publisher Full Text OpenURL

  6. Bick, C, Rabinovich, MI: On the occurrence of stable heteroclinic channels in Lotka-Volterra models. Dyn. Syst.. 25, 97–110 (2010). Publisher Full Text OpenURL

  7. Seliger, P, Tsimring, LS, Rabinovich, MI: Dynamics-based sequential memory: winnerless competition of patterns. Phys. Rev. E (3). 67, (2003)

  8. Rabinovich, MI, Afraimovich, VS, Varona, P: Heteroclinic binding. Dynamical Systems. 25(3), 433–442 (2010). Publisher Full Text OpenURL

  9. Wordsworth, J, Ashwin, P: Spatiotemporal coding of inputs for a system of globally coupled phase oscillators. Phys. Rev. E. 78, (2008)

  10. Ashwin, P, Orosz, G, Wordsworth, J, Townley, S: Dynamics on networks of clustered states for globally coupled phase oscillators. SIAM J. Appl. Dyn. Syst.. 6(4), 728–758 (2007). Publisher Full Text OpenURL

  11. Nowotny, T, Rabinovich, MI: Dynamical origin of independent spiking and bursting activity in neural microcircuits. Phys. Rev. Lett.. 98, (2007)

  12. Palis, J, de Melo, W: Geometric Theory of Dynamical Systems, Springer-Verlag, New York (1982) An introduction, translated from the Portuguese by A.K. Manning

  13. Homburg, AJ, Sandstede, B: Homoclinic and heteroclinic bifurcations in vector fields. Handbook of Dynamical Systems III. 379–524 (2010)

  14. Golubitsky, M, Stewart, I: The Symmetry Perspective, Birkhäuser Verlag, Basel (2003)

  15. Aguiar, MAD, Ashwin, P, Dias, APS, Field, M: Dynamics of coupled cell networks: synchrony, heteroclinic cycles and inflation. J. Nonlinear Sci.. 21, 271–323 (2011). Publisher Full Text OpenURL

  16. Field, M: Combinatorial dynamics. Dynamical Systems. 19, 217–243 (2004). Publisher Full Text OpenURL

  17. Golubitsky, M, Stewart, I: Nonlinear dynamics of networks: the groupoid formalism. Bull. Am. Math. Soc.. 43, 305–364 (2006). Publisher Full Text OpenURL

  18. Field, M: Lectures on Bifurcations, Dynamics and Symmetry, Longman, Harlow (1996)

  19. Ashwin, P, Field, M: Heteroclinic networks in coupled cell systems. Arch. Ration. Mech. Anal.. 148(2), 107–143 (1999). Publisher Full Text OpenURL

  20. Chawanya, T, Ashwin, P: A minimal system with a depth-two heteroclinic network. Dynamical Systems. 25, 397–412 (2010). Publisher Full Text OpenURL

  21. Stewart, I, Golubitsky, M, Pivato, M: Symmetry groupoids and patterns of synchrony in coupled cell networks. SIAM J. Appl. Dyn. Sys.. 2(4), 609–646 (2003). Publisher Full Text OpenURL

  22. Ashwin, P, Burylko, O, Maistrenko, Y: Bifurcation to heteroclinic cycles and sensitivity in three and four coupled phase oscillators. Physica D. 237, 454–466 (2008). Publisher Full Text OpenURL

  23. Karabacak, O, Ashwin, P: Heteroclinic ratchets in networks of coupled oscillators. J. Nonlinear Sci.. 20, 105–129 (2010). Publisher Full Text OpenURL

  24. Podvigina, O, Ashwin, P: On local attraction properties and a stability index for heteroclinic connections. Nonlinearity. 24(3), 887–929 (2011). Publisher Full Text OpenURL

  25. Rabinovich, MI, Varona, P, Selverston, AI, Abarbanel, HDI: Dynamical principles in neuroscience. Rev. Mod. Phys.. 95, 519–536 (2006)

  26. Fukai, T, Tanaka, S: A simple neural network exhibiting selective activation of neuronal ensembles: from winner-take-all to winners-share-all. Neural Comput.. 9, 77–97 (1997). PubMed Abstract | Publisher Full Text OpenURL

  27. Proceedings of the International Conference on Interval Methods and Computer Aided Proofs in Science and Engineering - Interval’96. Springer, Dordrecht (1997) [Held in Würzburg, September 30-October 2, 1996, edited by J. Wolff von Gudenberg, Reliab. Comput. 3(3) (1997)]

  28. Ashwin, P, Swift, JW: The dynamics of n weakly coupled identical oscillators. J. Nonlinear Sci.. 2, 69–108 (1992). Publisher Full Text OpenURL

  29. Ashwin, P, Burylko, O, Maistrenko, Y, Popovych, O: Extreme sensitivity to detuning for globally coupled phase oscillators. Phys. Rev. Lett.. 96(5), (2006)

  30. Hansel, D, Mato, G, Meunier, C: Clustering and slow switching in globally coupled phase oscillators. Phys. Rev. E. 48(5), 3470–3477 (1993). Publisher Full Text OpenURL

  31. Rinzel, J: Excitation dynamics: insights from simplified membrane models. Fed. Proc.. 44, 2944–2946 (1985). PubMed Abstract OpenURL

  32. Rubin, JI: Surprising effects of synaptic excitation. J. Comput. Neurosci.. 18(3), 333–342 (2005). PubMed Abstract | Publisher Full Text OpenURL

  33. Hansel, D, Mato, G, Meunier, C: Phase dynamics for weakly coupled Hodgkin-Huxley neurons. Europhys. Lett.. 23, 367–372 (1993). Publisher Full Text OpenURL

  34. Komarov, MA, Osipov, GV, Suykens, JAK: Sequentially activated groups in neural networks. Europhys. Lett.. 86(6), (2009)

  35. Krupa, M, Melbourne, I: Asymptotic stability of heteroclinic cycles in systems with symmetry. Ergod. Theory Dyn. Syst.. 15, 121–147 (1995)

  36. Beggs, J, Plenz, D: Neuronal avalanches are diverse and precise activity patterns that are stable for many hours in cortical slice cultures. J. Neurosci.. 24, 5216–5229 (2004). PubMed Abstract | Publisher Full Text OpenURL

  37. Mazor, O, Laurent, G: Transient dynamics versus fixed points in odor representations by locust antennal lobe projection neurons. Neuron. 48, 661–673 (2005). PubMed Abstract | Publisher Full Text OpenURL

  38. Perez-Orive, J, Mazor, O, Turner, GC, Cassenaer, S, Wilson, RI, Laurent, G: Oscillations and sparsening of odor representations in the mushroom body. Science. 297, 359–365 (2002). PubMed Abstract | Publisher Full Text OpenURL

  39. Ashwin, P, Borresen, J: Encoding via conjugate symmetries of slow oscillations for globally coupled oscillators. Phys. Rev. E. 70(2), (2004)

  40. Neves, FS, Timme, M: Controlled perturbation-induced switching in pulse-coupled oscillator networks. J. Phys. A. 42, (2009)