## Xt 4t t t

Typically \$ is taken to be a constant, i.e., x(t) = x0, —t < t < 0 , which is reasonable for most experimental systems. It should be noted that only solutions which are asymptotically stable can be accurately approximated using numerical integration.

There are two main programs available for the numerical integration of delay differential equations. The widely-used (and free) package XPPAUT (Ermentrout, 2005) can perform numerical integration using a variety of fixed step numerical methods, including Runge-Kutta. It has a good graphical user interface for visualizing the results. Perhaps the most useful aspect of this program is the ease with which parameters and initial conditions can be changed. The recent book of Ermentrout (2002) gives a overview of the package including many examples. Information on how to download the package as well as documentation and tutorials are available at www.math.pitt.edu/~bard/xpp/xpp.html. Within MATLAB there is the function DDE23 (Shampine and Thompson, 2001) which is a variable step size numerical integration routine for delay differential equations. A tutorial is on this routine available at www.mathworks.com/dde_tutorial. Results maybe visualized using the extensive graphing tools of MATLAB.

Numerical bifurcation analysis consists of two parts, the approximation of a solution and the calculation of the stability of this solution. The approximation of a solution in a numerical bifurcation package is not done using numerical integration, but rather using numerical continuation. Numerical continuation uses a given solution for a particular parameter value to find a solution for a different (but close) parameter value. This is only easily implemented for equilibrium and periodic solutions. Both stable and unstable solutions can be found. Once an equilibrium solution is found to a desired accuracy, approximations for a finite set of the eigenvalues with the largest real part can be determined, which will determine the stability of the equilibrium point. The stability of periodic orbits can be numerically determined in a similar way. Numerical bifurcation packages generally track the stability of equilibrium points and periodic orbits, indicating where bifurcations occur.

There is one package available that does numerical bifurcation analysis for delay differential equations, DDE-BIFTOOL (Engelborghs et al., 2001). This package runs on MATLAB. An overview of the numerical methods used in this package and some examples applications can be found in the paper of Engelborghs et al. (2002). The user manual and information on how to download the package are available at

www.cs.kuleuven.ac.be/cwis/research/twr/research/software/delay/

3 Effects of delay

In this section I will outline some of the effects of delay that have been documented in the literature.

### 3.1 Creation of Oscillations

Time delays are commonly associated with type II oscillations, i.e. oscillations created by a Hopf bifurcation (Breakspear and Jirsa, 2006), for the following reason. There are many examples of systems that have a stable equilibrium point if the time delay is zero (or sufficiently small), but have oscillatory behaviour if the delay is large enough. In these systems, the oscillation is created via a Hopf bifurcation at a critical value of the delay. This is sometimes referred to as a delay-induced oscillation. One of the simplest examples of this is the following model for recurrent inhibition due to Plant (1981):

v(t) = v(t) — 1 v3(t) — w(t) + c(v(t — t ) — vo) W(t) = p(v(t) + a — bw(t)) .

This is a Fitzhugh-Nagumo model neuron with a delayed term which represents recurrent feedback. Plant considered parameters such that the system with no feedback has a stable equilibrium point and showed that this stability is maintained for the system with feedback and sufficiently small delay. He then showed that when c < 0 (i.e. the recurrent feedback is inhibitory), there is a Hopf bifurcation at a critical value of the delay, leading to oscillations.

### 3.2 Oscillator Death

One of the most publicized (Strogatz, 1998) effects of time delays is the fact that the presence of time delays in the coupling between oscillators can destroy the oscillations. This phenomenon, usually called oscillator death or amplitude death was first noted by Ramana Reddy et al. (1998), in their analysis of a simple model of type II oscillators with gap junctional coupling. Subsequently Ramana Reddy et al. (2000) observed this phenomenon experimentally in a system of two intrinsically oscillating circuits with the same type of coupling. There are many papers related to delay induced oscillator death in the coupled oscillator literature, which I will not attempt to review here. Instead I will focus the discussion on results relevant to neural models.

The work of Ramana Reddy et al. (1998, 1999) shows that when two or more intrinsically oscillating elements are connected with gap junctional coupling of sufficient strength with a sufficiently large delay then the oscillations may be destroyed. Their work focused on systems where the elements were identical except for the frequency of the intrinsic oscillations and the coupling was all-to-all and symmetric (all the coupling coefficients were the same). Their model oscillator was just the normal form for the Hopf bifurcation. This behaviour has also been seen for a delayed, linearly coupled (i.e. (3) with no xj(t) term) pair of van der Pol oscillators (Wirkus and Rand, 2002), and for a pair of Fitzhugh-Nagumo oscillators with delayed gap junctional coupling (Campbell and Smith, 2007). To my knowledge this has yet to be observed for other biophysical models of neural oscillators, however, it may be expected to occur for most type II oscillators. Atay (2003b) obtained results for a network of weakly nonlinear oscillators with a symmetric connection matrix and gap junctional coupling. He showed that if the intrinsic frequency of the oscillators is sufficiently similar then oscillator death can occur.

Several studies have shown that the type of oscillator death described above does not occur for type II oscillators with sigmoidal coupling (Buric and Todorovic, 2003; Campbell et al., 2004; Shayer and Campbell, 2000). However, a different type of oscillator death can occur (Buric and Todorovic, 2003; Buric et al., 2005; Campbell et al., 2004; Shayer and Campbell, 2000): for elements which are intrinsically excitable (i.e. not oscillating when decoupled), oscillations induced by instantaneous coupling may be lost if a time delay is introduced.

The work of Buric et al. (2005) has shown that for the type I oscillator of Terman and Wang (1995), there is no oscillator death of this latter type with either gap junctional or sigmoidal coupling. Their work also suggests that delay induced oscillator death of the first type is not possible.

The study of type II oscillator death in coupled neural systems combines various techniques of Sect. 2. Oscillator death can occur when increasing the time delay causes the stabilization of an equilibrium point. Values of the delay where this occurs will correspond to places where the characteristic (19) has an eigenvalue with zero real part and dR^X) < 0. To have oscillator death, however, one must also show that the periodic orbit is eliminated. This means that at the value of t where the equilibrium point stabilizes, there is a "reverse" Hopf bifurcation destroying the stable limit cycle. This may be checked via numerical simulations or numerical continuations (see subsection 2.5), or by showing, as outlined in subsection 2.2, that the Hopf bifurcation is subcritical. Buric et al. (2005) and Buric and Todorovic (2003, 2005) have shown that for excitable Fitzhugh-Nagumo neurons, the restabilization of the equilibrium point is not always accompanied by oscillator death. In the case that the Hopf bifurcation is subcritical, the stable oscillator may persist with the stable equilibrium point giving a region of bistability. In their model, for larger values of t the periodic orbit is eliminated in a saddle-node bifurcation of limit cycles, leading to oscillator death.

The results of Buric et al. on type I oscillator death are primarily based on numerical simulations. To my knowledge there has been virtually no mathematical study of this situation. Recall that type I oscillators are those where the oscillation is created by an infinite period bifurcation (Breakspear and Jirsa, 2006). If such a bifurcation takes place in the coupled system with no time delay, introducing a time delay will not change the presence of the saddle-node bifurcation, however, it may affect whether this bifurcation occurs on an invariant circle. Continuity arguments would suggest that for sufficiently small delay, the saddle-node bifurcation will still occur on the invariant circle, leading to the creation of a periodic orbit at exactly the same bifurcation point as for the undelayed system. What happens for large delay remains to be investigated.

### 3.3 Attractor Switching and Multistability

A significant observation about ANNs of the form (14), is that many intersections between different Hopf bifurcation curves and between Hopf bifurcation curves and pitchfork bifurcation curves can occur (Belair et al., 1996; Shayer and Campbell, 2000; Yuan and Campbell, 2004). Figure 1 shows that this occurs in our coupled Fitzhugh-Nagumo model as well. These intersection points are called codimension two bifurcation points. Such points can lead to more complicated dynamics including: the existence of solutions with multiple frequencies (quasiperiodicity), the coexistence of more than one stable solution (multistability) or the switching of the system from one type of solution to another as a parameter is varied (Guckenheimer and Holmes, 1983, Chap. 7), (Kuznetsov, 1995, Chap. 8). In ordinary differential equations, such points are quite rare. In delay differential equations, however, such points are more common as the time delay forces there to be multiple branches of Hopf bifurcation.

In the ANN models, the following behaviour associated with the codimen-sion two points has been observed (Belair et al., 1996; Campbell et al., 2005; Shayer and Campbell, 2000; Yuan and Campbell, 2004): (i) multistability between a periodic solution and one or more equilibrium points; (ii) bistability between two periodic solutions (both synchronous or one synchronous and one asynchronous); and (iii) switching from one stable solutions to another as the delay is changed for a fixed coupling strength or as the coupling strength is changed for a fixed delay. The switching in (iii) may take place through a region of bistability or a region where the trivial solution is stable. Note that situation (i) leads to a different type of oscillator death than that discussed in the previous subsection: a slight perturbation can cause the system to switch from the stable oscillatory solution to the stable equilibrium solution, with no change in the parameter values.

Most of this behaviour has been confirmed in systems with biophysically relevant models for the neurons. In their studies of rings of Fitzhugh-Nagumo oscillators with time delayed gap-junctional or sigmoidal coupling, BuriC and TodoroviC (2003) and Buric et al. (2005) have documented almost all the behaviour observed in the ANN models including switching between different oscillation patterns and bistability between different oscillation patterns. For a system of two van der Pol oscillators with linear delayed coupling (i.e. (3) with no Xj(t) term), Sen and Rand (2003) have numerically observed and Wirkus and Rand (2002) have analytically proven the following sequence as the time delay is increased: in-phase oscillations ^ bistability between in-phase and anti-phase oscillations ^ anti-phase oscillations. They also observed the reverse sequence for different values of the coupling strength. Delay-induced bistability between in-phase oscillations and suppression oscillations (i.e. one cell oscillates and the other is quiescent) has been observed in models of hip-pocampal interneurons (Skinner et al., 2005a,b). Here the delay was synaptic and modelled via an extra equation representing the chemical kinetics of the synapse. Bistability between different types of travelling pulses has been observed in certain integrate-and-fire networks with delayed excitatory synaptic connections (Golomb and Ermentrout, 1999, 2000). In particular, they observe a switch from continuous travelling pulses to lurching travelling pulses as the time delay is increased with a transition region where there is bista-bility between the two types. This behaviour seems to be associated with a subcritical Hopf bifurcation.

Foss et al. (1996) and Milton and Foss (1997) have studied multistability in models for a delayed recurrent neural loop. Their model consists of a single excitatory neuron with delayed inhibitory feedback. They showed that up to three stable oscillatory patterns can coexist and that switching between the attractors can be induced by small perturbations in the neuron voltage (Foss et al., 1996) or by noise (Foss et al., 1997). These results have been replicated in experimental studies of a hybrid neural computer device consisting of an Aplysia motorneuron dynamically clamped to a computer which provides the delayed feedback (Foss and Milton, 2000, 2002). A possible cause of the multistability in these delayed feedback systems maybe period doubling bifurcations (Ikeda and Matsumoto, 1987). Bistability between different oscillation patterns was also observed in preparations of small Aplysia neural circuits (Kleinfeld et al., 1990).

Bifurcation induced transitions between different attractors have been observed in several experiments. In an experimental electrical circuit system, Ramana Reddy et al. (2000) have observed the sequence: in-phase oscillations ^ no oscillations ^ anti-phase oscillations as the time delay in the (gap-junctional) coupling is increased. Transitions from in-phase to anti-phase oscillations have been observed in human bimanual coordination experiments (Kelso et al., 1981; Kelso, 1984; Carson et al., 1994); see also the review article of Jantzen and Kelso (2006). One model which explains these experiments incorporates time delays in the coupling (Haken et al., 1985).

### 3.4 Synchronization

There are several approaches to studying synchronization. I will not review the details here, but give some indication which of these have been extended to delay differential equations and what the results are.

There is a very large literature on synchronization in artificial neural networks, some of which addresses systems with time delays (Campbell et al., 2006; Wu et al., 1999; Yuan and Campbell, 2004; Zhou et al., 2004a,b). Most of these papers use Lyapunov functionals to show that the all solutions synchronize as t ^ to, for appropriate parameter values. Although the equations of the individual elements are not relevant for modelling biophysical neurons, the techniques of analysis may be carried over to neural systems. A common conclusion in many of these papers is that if the strength of the coupling is small enough, one can achieve synchronization for all t > 0. However, synchronization may mean that all elements asymptotically approach the same equilibrium point.

As I have mentioned elsewhere in this chapter, a basic principle of delay differential equations such as (13) is that the behaviour of the system for small delay is often qualitatively similar to that for zero delay. Thus if the neurons are synchronized for a given value of the coupling with zero delay they should remain synchronized for small enough delays in the coupling. Unfortunately, quantifying "small enough" may be difficult and will generally depend on the particular neural model involved. Recall the example illustrated in Fig. 1. We showed that for c > 0 large enough (sufficiently large excitatory coupling) the undelayed system exhibits synchronized oscillations. We expect these oscillations to persist for t > 0 at least until one reaches the first thick Hopf bifurcation curve where synchronous oscillations are destroyed. (If the Hopf bifurcation is subcritical, the oscillations may persist above the curve). Thus, for this particular example, the Hopf bifurcation curve gives a lower bound on "how small" the delay must be to preserve the synchronization found for zero delay. Note that this does not preclude synchronization occurring for larger values of the delay, which is the case in this example. A similar situation is seen for coupled van der Pol oscillators in (Wirkus and Rand, 2002). Another example is the work of Fox et al. (2001) who studied relaxation oscillators with excitatory time delayed coupling. They showed that synchrony achieved for zero delay is preserved for delays up to about 10% of the period of the oscillation, for a variety of different models. The one exception is when the right hand side of the equation is not a differentiable function, in which case synchronization is lost for t > 0. Crook et al. (1997) observed a similar phenomenon for a continuum model of the cortex, with excitatory coupling and distance dependent delays. Namely, they found for small enough delay the synchronous oscillation is stable, but for larger delays this oscillation loses stability to a travelling wave.

More complicated situations occur when both excitatory and inhibitory connections exist. Ermentrout and Kopell (1998); Kopell et al. (2000); Kar-bowski and Kopell (2000) have studied a model for hippocampal networks of excitatory and inhibitory neurons where two types of synchronous oscillation are possible. They show that persistence of the synchronous oscillations with delays depends subtly on the currents present in the cells and the connections present between cells.

So far I have discussed synchronization in spite of delays. I now move on to the more interesting case of synchronization because of delays. This situation can occur when there are inhibitory synaptic connections in the network. This has been extensively documented and studied when the delay is modelled by slow kinetics of the synaptic gating variable (van Vreeswijk et al., 1994; Wang and Buzsaki, 1998; Wang and Rinzel, 1992, 1993; White et al., 1998). Further, Maex and De Schutter (2003) suggest that the type of delay is not important, just the fact that it leads to a separation in time between when the pre-synaptic neuron generates an action potential and the post-synaptic neuron receives it. They confirm this for a network of multi-compartment model neurons with fast synaptic kinetics and a discrete conduction delay. This idea is further supported by the observation of synchronization via discrete delayed inhibition in a number of artificial neural network models (Campbell et al., 2004, 2005; Shayer and Campbell, 2000). Finally we illustrate this with our coupled Fitzhugh-Nagumo model. Consider the part of Fig. 1 with c < 0 (inhibitory coupling). For sufficiently large coupling strength and zero delay the system tends to an asynchronous phase-locked state. This state persists for t > 0 sufficiently small, however, for t large enough a stable synchronous state may be created in the Hopf bifurcation corresponding to the thin curve.

Only a few studies have looked at synchronization with time delayed gap-junctional coupling. One example is the work of Dhamala et al. (2004) which shows that for two gap junctional coupled Hindmarsh-Rose neurons synchronization is achieved for smaller coupling strengths if there is a nonzero time delay in the coupling. Another is the work of Buric et al. (2005).

### 4 Distributed Delays

There are very few results concerning neural systems with distributed delays, thus I will review some general results, mostly from the population biology literature, which should carry over to neural systems. What has emerged from this literature is a general principle that a system with a distribution of delays is inherently more stable than the same system with a discrete delay. Some specific results to support this are described below.

Bernard et al. (2001) analyzed the linear stability of a scalar system with one and two delays in terms of generic properties of the distribution g, such as the mean, variance and skewness. For the uniform and continuous distributions, they have shown that stability regions are larger than those with a discrete delay.

Jirsa and Ding (2004) have analyzed an n x n linear system with linear decay and arbitrary connections with a common delay. They have shown, under some mild assumptions, that the stability region of the trivial solution for any distribution of delays is larger than and contains the stability region for a discrete delay.

Campbell and Ncube (2006) have shown that it is more difficult to get delay induced oscillations with distributions of delays of the form (6) with Tm = 0. For large variance (m =1) delay induced instability is impossible and for smaller variance (m > 1) the mean delay needed for instability is much larger than the discrete delay value. They have also shown that sufficiently small variance in the distribution is needed to get the bifurcation interactions which may lead to multistability, oscillator death and attractor switching discussed above.

Atay (2003a, 2006) has studied the same model as Ramana Reddy et al. (1998) only with distributed delays of the form (6) with g given by (7). He shows it is easier to destroy oscillations with a distribution of delays than with a discrete delay, in the sense that there is a larger region of oscillator death in the parameter space consisting of the mean delay and the strength of the coupling. As the variance of the distribution increases the size of this region increases.

Thiel et al. (2003) studied a scalar equation representing a mean field approximation for a population pyramidal cells with recurrent feedback, first formulated by Mackey and an der Heiden (1984). They show that having a uniform distribution of delays simplifies the dynamics of the system. The size of the stability region of the equilibrium point is larger and larger mean delays are needed to induce oscillations. Complex phenomena such as chaos are less likely to occur, or totally precluded if the variance of the distribution is sufficiently large. The model with a distribution of delays better explains the appearance of periodic bursts of activity when penicillin is added to a hippocampal slice preparation (which reduces the coupling strength).

### 5 Summary and Future Directions

In this chapter I showed how time delays due to conduction along the axon or dendrite or due to transmission across the synapse could be modelled with delay differential equations. I outlined some of the tools available for analyzing such equations and reviewed some of the literature about such models. Some key observations are:

- Time delays can lead to the creation of type II oscillations, especially in systems with delayed inhibitory coupling.

- Time delays can destroy type II oscillations in a network of intrinsically oscillatory neurons with gap junctional coupling.

- If a system has a stable synchronous oscillation when there is no delay in the coupling, the solution remains stable for small enough delay, but may lose stability for larger delay.

- A system with inhibitory coupling which does not have a stable synchronous oscillation for zero delay, may have one if the delay is large enough.

- Time delays may lead to bistability between different type II oscillatory solutions (e.g. synchronous and anti-phase) or switching between different type II oscillatory solutions.

There are a number of problems which still require further study. These include: determining the effect of delay on the generation and destruction of type I oscillations (infinite period bifurcations), applying and/or extending the methods used to study synchronization in artificial neural networks to biophysical neural networks, and studying the effect of distributions of delays on biophysical neural networks.

References

Atay FM (2003a) Distributed delays facilitate amplitude death of coupled oscillators.

Phys. Rev. Lett. 91 (9), 094101 Atay FM (2003b) Total and partial amplitude death in networks of diffusively coupled oscillators. Physica D 183, 1-18 Atay FM (2006) Oscillator death in coupled functional differential equations near

Hopf bifurcation. J. Differential Equations 221 (1), 190-209 Belair J, Campbell SA, van den Driessche P (1996) Frustration, stability and delay-induced oscillations in a neural network model. SIAM J. Applied Mathematics 56 (1), 245-255

Bernard S, Belair J, Mackey MC (2001) Sufficient conditions for stability of linear differential equations with distributed delay. Discrete and Continuous Dynamical Systems 1B, 233-256 Breakspear M, Jirsa VK (2006) Neuronal dynamics and brain connectivity. In: Mcintosh R, Jirsa VK (eds), Handbook of Brain Connectivity. Springer-Verlag, New York

Buric N, Grozdanovic I, Vasovic N (2005) Type I vs. type II excitable systems with delayed coupling. Chaos, Solitons and Fractals 23, 1221-1233 Buric N, Todorovic D (2003) Dynamics of Fitzhugh-Nagumo excitable systems with delayed coupling. Physical Review E 67, 066222 Buric N, Todorovic D (2005) Bifurcations due to small time-lag in coupled excitable systems. International Journal of Bifurcations and Chaos 15 (5), 1775-1785 Campbell SA, Edwards R, van den Dreissche P (2004) Delayed coupling between two neural network loops. SIAM J. Applied Mathematics 65 (1), 316-335 Campbell SA, Ncube I (2006) Some effects of gamma distribution on the dynamics of a scalar delay differential equation. Preprint Campbell SA, Ncube I, Wu J (2006) Multistability and stable asynchronous periodic oscillations in a multiple-delayed neural system. Physica D 214 (2), 101-119 Campbell SA, Smith A (2007) Phase models and delayed coupled Fitzhugh-Nagumo oscillators. Preprint

Campbell SA, Yuan Y, Bungay S (2005) Equivariant Hopf bifurcation in a ring of identical cells with delayed coupling. Nonlinearity 18, 2827-2846 Carson RG, Byblow WD, Goodman D (1994) The dynamical substructure of bimanual coordination. In: Swinnen S, Heuer H, Massion J, Casaer P (eds), Interlimb coordination: Neural, dynamical and cognitive constraints. Academic Press, San Diego, pp 319-337

Crook SM, Ermentrout GB, Vanier MC, Bower JM (1997) The role of axonal delay in the synchronization of networks of coupled cortical oscillators. J. Computational Neuroscience 4, 161-172 Cushing JM (1977) Integrodifferential Equations and Delay Models in Popluation Dynamics. Vol. 20 of Lecture Notes in Biomathematics. Springer-Verlag, Berlin; New York

Desmedt JE, Cheron G (1980) Central somatosensory conduction in man: Neural generators and interpeak latencies of the far-field components recorded from neck and right or left scalp and earlobes. Electro. Clin. Electroencephalog. 50, 382-403

Dhamala M, Jirsa VK, Ding M (2004) Enhancement of neural synchrony by time delay. Physical Review E 92 (7), 074104 Diekmann O, van Gils SA, Verduyn Lunel SM, Walther H-O (1995) Delay Equations.

Springer-Verlag, New York Engelborghs K, Luzyanina T, Roose D (2002) Numerical bifurcation analysis of delay differential equations using DDE-BIFTOOL. ACM Transactions on Mathematical Software 28 (1), 1-21 Engelborghs K, Luzyanina T, Samaey G (2001) DDE-BIFTOOL v. 2.00: a MAT-LABpackage for bifurcation analysis of delay differential equations. Tech. Rep. TW-330, Department of Computer Science, K.U. Leuven, Leuven, Belgium Ermentrout GB (1994) An introduction to neural oscillators. In: Ventriglia F (ed),

Neural Modeling and Neural Networks. Pergamon, Oxford, UK, pp 79-110 Ermentrout GB (2002) Similating, Analyzing and Animating Dynamical Systems:

A Guide to XPPAUT for Researcher and Students. SIAM, Philadelphia, PA Ermentrout GB (2005) XPPAUT5.91 - the differential equations tool. Department of Mathematics, University of Pittsburgh, Pittsburgh, PA Ermentrout GB, Kopell N (1998) Fine structure of neural spiking and synchronization in the presence of conduction delays. PNAS 95 (3), 1259-64 Fitzhugh R (1960) Thresholds and plateaus in the Hodgkin-Huxley nerve equations.

J. General Physiology 43, 867 Foss J, Longtin A, Mensour B, Milton JG (1996) Multistability and delayed recurrent loops. Phys. Rev. Letters 76, 708-711 Foss J, Milton JG (2000) Multistability in recurrent loops arising from delay. J.

Neurophysiol. 84, 975-985 Foss J, Milton JG (2002) Multistability in delayed recurrent loops. In: Milton J, Jung P (eds), Epilepsy as a Dynamic Disease. Springer-Verlag, New York, pp 283-295

Foss J, Moss F, Milton JG (1997) Noise, multistability and delayed recurrent loops.

Phys. Rev. E 55, 4536-4543 Fox JJ, Jayaprakash C, Wang DL, Campbell SR (2001) Synchronization in relaxation oscillator networks with conduction delays. Neural Computation 13, 1003-1021

Golomb D, Ermentrout GB (1999) Continuous and lurching travelling pulses in neuronal networks with delay and spatially decaying connectivity. PNAS 96, 13480-13485

Golomb D, Ermentrout GB (2000) Effects of delay on the type and velocity of travelling pulses in neuronal networks with spatially decaying connectivity. Network: Comput. Neural Syst. 11, 221-246 Gopalsamy K, Leung I (1996) Delay induced periodicity in a neural netlet of excitation and inhibition. Physica D 89, 395-426 Guckenheimer J, Holmes PJ (1983) Nonlinear Oscillations, Dynamical Systems and

Bifurcations of Vector Fields. Springer-Verlag, New York Haken H, Kelso JAS, Bunz H (1985) A theoretical model of phase transitions in human hand movements. Biological Cybernetics 51, 347-356 Hale JK, Verduyn Lunel SM (1993) Introduction to Functional Differential Equations. Springer-Verlag, New York Hoppensteadt FC, Izhikevich EM (1997) Weakly connected neural networks.

Springer-Verlag, New York Ikeda I, Matsumoto K (1987) High dimensional chaotic behaviour in systems with time-delayed feedback. Physica D 29, 223-235 Izhikevich EM (1998) Phase models with explicit time delays. Physical Review E 58 (1), 905-908

Jantzen KJ, Kelso JAS (2006) Neural coordination dynamics of human sensorimo-tor behaviour: A review. In: McIntosh R, Jirsa VK (eds), Handbook of Brain Connectivity. Springer-Verlag, New York Jirsa VK, Ding M (2004) Will a large complex system with delays be stable? Physical

Review Letters 93, 070602 Karbowski J, Kopell N (2000) Multispikes and synchronization in a large neural network with temporal delays. Neural Computation 12, 1573-1606 Keener J, Sneyd J (1998) Mathematical Physiology. Springer-Velag, New York Kelso JAS (1984) Phase transitions and critical behaviour in human bimanual coordination. American J. Physiology: Regulatory, Integrative and Comparative Physiology 15, R1000-R1004 Kelso JAS, Holt KG, Rubin P, Kugler PN (1981) Patterns of human interlimb coordination emerge from nonlinear limit cycle oscillatory processes: theory and data. J. Motor Behaviour 13, 226-261 Kleinfeld D, Raccuia-Behling F, Chiel HJ (1990) Circuits constructed from identified Aplysia neurons exhibit multiple patterns of activity. Biophysical J. 57 (4), 697-715

Koch C (1999) Biophysics of Computation. Oxford University Press, New York Kolmanovskii VB, Nosov VR (1986) Stability of functional differential equations.

Vol. 180 of Mathematics in Science and Engineering. Academic Press Kopell N, Ermentrout GB (2002) Mechanisms of phase-locking and frequency control in pairs of coupled neural oscillators. In: Fiedler B (ed), Handbook of Dynamical Systems, vol 2: Toward Applications. Elsevier, Amsterdam Kopell N, Ermentrout GB, Whittington MA, Traub R (2000) Gamma rhythms and beta rhythms have different synchronization properties. PNAS 97 (4), 1867-1872 Kuznetsov YA (1995) Elements of Applied Bifurcation Theory. Vol. 112 of Applied

Mathematical Sciences. Springer-Verlag, Berlin; New York MacDonald N (1978) Time lags in biological models. Vol. 27 of Lecture notes in biomathematics. Springer-Verlag, Berlin; New York

Mackey MC, an der Heiden U (1984) The dynamics of recurrent inhibition. J. Mathematical Biology 19, 211-225 Maex R, De Schutter E (2003) Resonant synchronization in heterogeneous networks of inhibitory neurons. J. Neuroscience 23 (33), 10503-10514 Milton JG, Foss J (1997) Oscillations and multistability in delayed feedback control. In: Othmer HG, Adler FR, Lewis MA, Dallon JC (eds), The Art of Mathematical Modeling: Case Studies in Ecology, Physiology and Cell Biology. Prentice Hall, New York, pp 179-198 Nagumo J, Arimoto S, Yoshizawa S (1962) An active pulse transmission line simulating nerve axon. Proc. IRE 50, 2061-2070 Olgac N, Sipahi R (2002) An exact method for the stability analysis of time-delayed linear time-invariant (LTI) systems. IEEE Transactions on Automatic Control 47 (5), 793-797

Olgac N, Sipahi R (2005) Complete stability robustness of third-order LTI multiple time-delay systems. Automatica 41, 1413-1422 Plant RE (1981) A Fitzhugh differential-difference equation modeling recurrent neural feedback. SIAM J. Applied Mathematics 40 (1), 150-162 Ramana Reddy DV, Sen A, Johnston GL (1998) Time delay induced death in coupled limit cycle oscillators. Physical Review Letters 80, 5109-5112 Ramana Reddy DV, Sen A, Johnston GL (1999) Time delay effects on coupled limit cycle oscillators at Hopf bifurcation. Physica D 129, 15-34 Ramana Reddy DV, Sen A, Johnston GL (2000) Experimental evidence of time-delay-induced death in coupled limit-cycle oscillators. Physical Review Letters 85 (16), 3381-3384

Sen AK, Rand R (2003) A numerical investigation of the dynamics of a system of two time-delayed coupled relaxation oscillators. Communications on Pure and Applied Mathematics 2 (4), 567-577 Shampine LF, Thompson S (2001) Solving DDEs in MATLAB. Applied Numerical

Mathematics 37, 441-458 Shayer LP, Campbell SA (2000) Stability, bifurcation and multistability in a system of two coupled neurons with multiple time delays. SIAM J. Applied Mathematics 61 (2), 673-700

Shepherd G (1994) Neurobiology. Oxford University Press, New York Skinner FK, Bazzazi H, Campbell SA (2005a) Two-cell to n-cell heterogeneous, inhibitory networks: precise linking of multistable and coherent properties. J. Computational Neuroscience 18 (3), 343-352 Skinner FK, Chung JYJ, Ncube I, Murray PA, Campbell SA (2005b) Using heterogeneity to predict inhibitory model characteristics. J. Neurophysiology 93, 1898-1907

Stepan G (1989) Retarded Dynamical Systems. Vol. 210 of Pitman Research Notes in Mathematics. Longman Group, Essex Strogatz SH (1998) Death by delay. Nature 394, 316-317

Terman D, Wang DL (1995) Global competition and local cooperation in a network of neural oscillators. Physica D 81, 148-176 Thiel A, Schwegler H, Eurich CW (2003) Complex dynamics is abolished in delayed recurrent systems with distributed feedback times. Complexity 8 (4), 102-108 van Vreeswijk C, Abbott LF, Ermentrout GB (1994) When inhibition not excitation synchronizes neural firing. J. Computational Neuroscience 1, 313-321

Wang X-J, Buzsaki G (1998) Gamma oscillation by synaptic inhibition in a hip-

pocampal interneuronal network model. J. Neuroscience 16 (20), 6-16 Wang X-J, Rinzel J (1992) Alternating and synchronous rhythms in reciprocally inhibitory model neurons. Neural Computation 4, 84-97 Wang X-J, Rinzel J (1993) Spindle rhythmicity in the reticularis thalami nucleus: synchronization among mutually inhibitory neurons. Neuroscience 53, 899-904 White JA, Chow CC, Ritt J, Soto-Trevino C, Kopell N (1998) Synchronization and oscillatory dynamics in heterogeneous, mutually inhibited neurons. J. Computational Neuroscience 5, 5-16 Wirkus S, Rand R (2002) The dynamics of two coupled van der Pol oscillators with delay coupling. Nonlinear Dynamics 30, 205-221 Wischert W, Wunderlin A, Pelster A, Olivier M, Groslambert J (1994) Delay-induced instabilities in nonlinear feedback systems. Physical Review E 49 (1), 203-219

Wu J, Faria T, Huang YS (1999) Synchronization and stable phase-locking in a network of neurons with memory. Math. Comp. Modelling 30 (1-2), 117-138 Yuan Y, Campbell SA (2004) Stability and synchronization of a ring of identical cells with delayed coupling. J. Dynamics and Differential Equations 16 (1), 709-744 Zhou J, Chen T, Lan Z (2004a) Robust synchronization of coupled delayed recurrent neural networks. In: Advances in Neural Networks - ISNN 2004. Vol. 3173 of Lecture Notes in Computer Science. Springer-Verlag, Berlin, pp 144-149 Zhou S, Liao Z, Yu J, Wong K (2004b) Chaos and its synchronization in two-neuron systems with discrete delays. Chaos, Solitons and Fractals 21, 133-142