Pulse coupled oscillators and the phase resetting curve

Limit cycle oscillators that are coupled in a pulsatile manner are referred to as pulse coupled oscillators. In these oscillators, the interactions take the form of brief pulses such that the effect of one input dies out before the next is received. A phase resetting curve (PRC) keeps track of how much an input advances or delays the next spike in an oscillatory neuron depending upon where in the cycle the input is applied. PRCs can be used to predict phase locking in networks of pulse coupled oscillators. In some studies of pulse coupled oscillators, a specific form is assumed for the interactions between oscillators, but a more general approach is to formulate the problem assuming a PRC that is generated using a perturbation that approximates the input received in the real biological network. In general, this approach requires that circuit architecture and a specific firing pattern be assumed. This allows the construction of discrete maps from one event to the next. The fixed points of these maps correspond to periodic firing modes and are easier to locate and analyze for stability compared to locating and analyzing periodic modes in the original network directly. Alternatively, maps based on the PRC have been constructed that do not presuppose a firing order. Specific circuits that have been analyzed under the assumption of pulsatile coupling include one to one lockings in a periodically forced oscillator or an oscillator forced at a fixed delay after a threshold event, two bidirectionally coupled oscillators with and without delays, a unidirectional N-ring of oscillators, and N all-to-all networks.     

Synchronization of delayed coupled neurons in presence of inhomogeneity

In principle, two directly coupled limit cycle oscillators can overcome mismatch in intrinsic rates and match their frequencies, but zero phase lag synchronization is just achievable in the limit of zero mismatch, i.e., with identical oscillators. Delay in communication, on the other hand, can exert phase shift in the activity of the coupled oscillators. In this study, we address the question of how phase locked, and in particular zero phase lag synchronization, can be achieved for a heterogeneous system of two delayed coupled neurons. We have analytically studied the possibility of inphase synchronization and near inphase synchronization when the neurons are not identical or the connections are not exactly symmetric. We have shown that while any single source of inhomogeneity can violate isochronous synchrony, multiple sources of inhomogeneity can compensate for each other and maintain synchrony. Numeric studies on biologically plausible models also support the analytic results.     

Mode analysis of a tubular structure of coupled non-linear oscillators for digestive-tract modelling

A commonly accepted mathematical model for the slow-wave electrical activity of the gastro-intestinal tract of humans and animals comprises a set of interconnected nonlinear oscillators. Using a van der Pol oscillator with third-power conductance characteristics as the unit oscillator a number of structures have been analysed using a matrix Krylov-Bogolioubov method linearisation. The mode analysis of one-dimensional chains and two-dimensional arrays has been reported. In this paper the method has been extended to consider a tubular structure which is relevant to modelling small-intestinal rhythms. It is shown that this structure is capable of producing stable single models, non-resonant double modes and degenerated modes. General expressions are obtained for anm×n structure and examples given of two special conditions of 3×4 (i.e. odd numbers of oscillators in a ring) and 4×3 cases. The analytical results obtained for these two cases have been vertified experimentally using an electronic implementation of coupled van der Pol oscillators. Results obtained using fifth-power non-linear oscillators are summarised.     

Phase Resetting Curves Allow for Simple and Accurate Prediction of Robust N:1 Phase Locking for Strongly Coupled Neural Oscillators

Existence and stability criteria for harmonic locking modes were derived for two reciprocally pulse coupled oscillators based on their first and second order phase resetting curves. Our theoretical methods are general in the sense that no assumptions about the strength of coupling, type of synaptic coupling, and model are made. These methods were then tested using two reciprocally inhibitory Wang and Buzsáki model neurons. The existence of bands of 2:1, 3:1, 4:1, and 5:1 phase locking in the relative frequency parameter space was predicted correctly, as was the phase of the slow neuron's spike within the cycle of the fast neuron in which it occurred. For weak coupling the bands are very narrow, but strong coupling broadens the bands. The predictions of the pulse coupled method agreed with weak coupling methods in the weak coupling regime, but extended predictability into the strong coupling regime. We show that our prediction method generalizes to pairs of neural oscillators coupled through excitatory synapses, and to networks of multiple oscillatory neurons. The main limitation of the method is the central assumption that the effect of each input dies out before the next input is received.     

A van der Pol coupled-oscillator model as a basis for developing a system for suppressing MHD instabilities in a tokamak

The main parameters of tokamak discharges are known to be limited by large-scale MHD instabilities. Sometimes, the instabilities lead to a rapid (on time scales of tens of microseconds) disruption of the discharge current and to the release of all the energy stored in the plasma column at the discharge chamber wall. This process, which is called the disruptive instability, may have irreversible catastrophic consequences for the operation of a fusion reactor. In the present paper, a study is made of the dynamics of self-oscillations in systems of two and six van der Pol coupled oscillators. A van der Pol coupled-oscillator model is used to develop a multivariable feedback controller based on the combined principle of compensating for internal cross feedbacks within the object and introducing damping feedbacks in each control channel. By using mathematical simulation methods, it is shown that the controller designed guarantees the suppression of self-oscillations in a system of van der Pol oscillators over a fairly broad range of parameters of the object under control (and thereby provides the structural stability of the object). The nonlinear control system model makes it possible to suppress coupled MHD perturbations developing in a tokamak plasma.     

On partial contraction analysis for coupled nonlinear oscillators

We describe a simple yet general method to analyze networks of coupled identical nonlinear oscillators and study applications to fast synchronization, locomotion, and schooling. Specifically, we use nonlinear contraction theory to derive exact and global (rather than linearized) results on synchronization, antisynchronization, and oscillator death. The method can be applied to coupled networks of various structures and arbitrary size. For oscillators with positive definite diffusion coupling, it can be shown that synchronization always occurs globally for strong enough coupling strengths, and an explicit upper bound on the corresponding threshold can be computed through eigenvalue analysis. The discussion also extends to the case when network structure varies abruptly and asynchronously, as in flocks of oscillators or dynamic elements.     

On self synchronization--aspects of flagellar and ciliary beating

Mathematical models involving hydrodynamically coupled oscillators are proposed to explain some self-synchronization phenomena observed in flagella and cilia. It is shown using elementary mathematical analysis and computer simulations that the behaviors of the models proposed are similar to those observed in experiment.     

Two coupled neural oscillators as a model of the circadian pacemaker

Pittendrigh first found that the circadian rhythm of locomotor activity in nocturnal rodents split into two components. Hoffman then reported that the splitting phenomenon was even more reproducible in the small diurnal primate Tupaia. These “splitting” experiments and many other experiments suggest that two coupled oscillators may constitute the circadian pacemaker system. Pittendrigh proposed a phenomenological two-oscillator model. Daan and Berde developed a quantitative model assuming that the interaction between the two constituent oscillators is by instantaneous resets. Their model system can simulate several qualitative features in the experimental data. As the assumption of instantaneous resets seems to be unnatural, we study two limit cycle oscillators, which are coupled continuously to each other, as a model of the circadian pacemaker. We assume the following points, (i) One oscillator in a resting state does not affect another oscillator, (ii) Two oscillators are identical, (iii) The coupling is symmetrical. By the theory of Hopf bifurcation it is found that the general two-oscillator system has two stable periodic solutions. One is the in-phase solution where the two constituent oscillators oscillate in phase synchrony. Another is the anti-phase solution where the two oscillators oscillate 180 ° out of phase. The former corresponds to a single pattern of locomotor activity and the latter corresponds to a splitting pattern. Furthermore, we study specific two-neural oscillators, which are linearly coupled to each other. By the method of secondary bifurcation we find that the model shows simultaneous stability of the two alternative phase relationships and the hysteresis phenomena found in Tupaia. A natural period of the uncoupled constituent oscillator is longer than that of the in-phase solution but it is shorter than that of the anti-phase solution. This is in agreement with the data of Tupaia.     

Rapid synchronization through fast threshold modulation

Synchronization properties of locally coupled neural oscillators were investigated analytically and by computer simulation. When coupled in a manner that mimics excitatory chemical synapses, oscillators having more than one time scale (relaxation oscillators) are shown to approach synchrony using mechanisms very different from that of oscillators with a more sinusoidal waveform. The relaxation oscillators make critical use of fast modulations of their thresholds, leading to a rate of synchronization relatively independent of coupling strength within some basin of attraction; this rate is faster for oscillators that have conductance-based features than for neural caricatures such as the FitzHugh-Nagumo equations that lack such features. Computer simulations of one-dimensional arrays show that oscillators in the relaxation regime synchronize much more rapidly than oscillators with the same equations whose parameters have been modulated to yield a more sinusoidal waveform. We present a heuristic explanation of this effect based on properties of the coupling mechanisms that can affect the way the synchronization scales with array length. These results suggest that the emergent synchronization behavior of oscillating neural networks can be dramatically influenced by the intrinsic properties of the network components. Possible implications for perceptual feature binding and attention are discussed.Supported in part by NASA (NGT-50497)Supported in part by NSF (DMS-8901913), and NIMH-47150 Present address and address for correspondence: Department of Brain and Cognitive Sciences, Massachusetts Institute of Technology, E25-618, Cambridge, MA 02139, USA     

Integer-ratio entrainment in mutually-coupled non-linear oscillators

In recent years entrainment conditions for mutually-coupled, non-linear oscillators have been studied for a number of biomedical applications and using different analytical methods. The emphasis has been on entrainment between oscillators of similar frequencies. In this paper entrainment conditions are considered for oscillators having intrinsic frequency ratio of about 3:1 and which exhibit integer-ratio synchronization. This condition has application in the study of blood pressure regulation particularly in relation to respiratory effects. Coupling has been observed between respiration and the vasomotor activity associated with the baroreceptor reflex, which has an intrinsic 0·1 Hz component. At normal breathing frequencies the frequency ratio of the respiratory and vasomotor components is in the region of 3:1 hence integer-ratio entrainment is feasible. Using a coupled van der Pol model the entrainment zones for different parameters are described. The parameters considered allow for varying amounts of output, output rate and delay in the intercoupling structure. In particular, it is shown that the entrainment regions are strongly affected by the nature of the coupling. Within these zones the harmonic balance method is developed to provide an analytical solution to frequency, amplitudes and phase conditions. The assumed solution is valid only for certain regions of the stability zones and the reason for this is demonstrated and the means whereby this can be overcome are indicated.     

Synchronization of coupled nonidentical genetic oscillators

The study of the collective dynamics of synchronization among genetic oscillators is essential for the understanding of the rhythmic phenomena of living organisms at both molecular and cellular levels. Genetic oscillators are biochemical networks, which can generally be modelled as nonlinear dynamic systems. We show in this paper that many genetic oscillators can be transformed into Lur'e form by exploiting the special structure of biological systems. By using a control theory approach, we provide a theoretical method for analysing the synchronization of coupled nonidentical genetic oscillators. Sufficient conditions for the synchronization as well as the estimation of the bound of the synchronization error are also obtained. To demonstrate the effectiveness of our theoretical results, a population of genetic oscillators based on the Goodwin model are adopted as numerical examples.     

Hydrocarbon-nitrosamine synergism as a possible amplifying factor in lung tumorigenesis by tobacco smoke.

Twenty-five coupled relaxation sawtooth oscillators have been investigated for the occurrence of mutual synchronization, using digital computer simulation. It is shown that mutual synchronization occurs only for relatively strong coupling. Synchronization depends heavily upon the precise waveform of the oscillators in the uncoupled state. The results are compared with a number of biological phenomena.     

Stabilization of metapopulation cycles: toward a classification scheme

The stability of population oscillations in ecological systems is considered. Experiments suggest that in many cases the single patch dynamics of predator-prey or host-parasite systems is extinction prone, and stability is achieved only when the spatial structure of the population is expressed via desynchronization between patches. A few mechanisms have been suggested so far to explain the inability of dispersal to synchronize the system. Here we compare a recently discovered mechanism, based on the dependence of the angular velocity on the oscillation amplitude, with other, already known conditions for desynchronization. Using a toy model composed of diffusively coupled oscillators we suggest a classification scheme for stability mechanisms, a scheme that allows for either a priori (based on the system parameters) or a posteriori (based on local measurements) identification of the dominant process that yields desynchronization.     

Analysis of clustered firing patterns in synaptically coupled networks of oscillators

Oscillators in networks may display a variety of activity patterns. This paper presents a geometric singular perturbation analysis of clustering, or alternate firing of synchronized subgroups, among synaptically coupled oscillators. We consider oscillators in two types of networks: mutually coupled, with all-to-all inhibitory connections, and globally inhibitory, with one excitatory and one inhibitory population of oscillators, each of arbitrary size. Our analysis yields existence and stability conditions for clustered states, along with formulas for the periods of such firing patterns. By using two different approaches, we derive complementary conditions, the first set stated in terms of time lengths determined by intrinsic and synaptic properties of the oscillators and their coupling and the second set stated in terms of model parameters and phase space structures directly linked to parameters. These results suggest how biological components may interact to produce the spindle sleep rhythm in thalamocortical networks. Received: 9 September 1999 / Revised version: 7 July 2000 / Published online: 24 November 2000     

Multiple pulse interactions and averaging in systems of coupled neural oscillators

Oscillators coupled strongly are capable of complicated behavior which may be pathological for biological control systems. Nevertheless, strong coupling may be needed to prevent asynchrony. We discuss how some neural networks may be designed to achieve only simple locking behavior when the coupling is strong. The design is based on the fact that the method of averaging produces equations that are capable only of locking or drift, not pathological complexity. Furthermore, it is shown that oscillators that interact by means of multiple pulses per cycle, dispersed around the cycle, behave like averaged equations, even if the number of pulses is small. We discuss the biological intuition behind this scheme, and show numerically that it works when the oscillators are taken to be composites, each unit of which is governed by a well-known model of a neural oscillator. Finally, we describe numerical methods for computing from equations for coupled limit cycle oscillators the averaged coupling functions of our theory.Research partially supported by the National Science Foundation under grants DMS 8796235 and DMS 8701405 and the Air Force Office of Scientific Research under University Research Contract F 49620-C-0131 to Northeastern University     

Interaction of coupled nonlinear oscillators having different intrinsic resting levels

The interaction among coupled oscillators is governed by oscillator properties (intrinsic frequency and amplitude) and coupling mechanisms. This study considers another oscillator property, the intrinsic resting level, and evaluates its role in governing oscillator interactions. The results of computer experiments on a chain of either three or five bidirectionally coupled nonlinear oscillators, suggest that an intrinsic resting level gradient, if present, is one of the factors governing the interaction between coupled oscillators. If there is no intrinsic frequency gradient, then an intrinsic resting level gradient is sufficient to produce many features of interaction among coupled oscillators. If both intrinsic frequency and intrinsic resting level gradients are present, then both of them determine the manner in which the coupled oscillators interact with each other.     

Numerical Method Using Cubic B-Spline for a Strongly Coupled Reaction-Diffusion System

In this paper, a numerical method for the solution of a strongly coupled reaction-diffusion system, with suitable initial and Neumann boundary conditions, by using cubic B-spline collocation scheme on a uniform grid is presented. The scheme is based on the usual finite difference scheme to discretize the time derivative while cubic B-spline is used as an interpolation function in the space dimension. The scheme is shown to be unconditionally stable using the von Neumann method. The accuracy of the proposed scheme is demonstrated by applying it on a test problem. The performance of this scheme is shown by computing and error norms for different time levels. The numerical results are found to be in good agreement with known exact solutions.     

Coupled Van der Pol oscillators — A model of excitatory and inhibitory neural interactions

A system of mutually coupled Van der Pol equations is derived from an extended version of the Wilson and Cowan model for the dynamics of a number of excitatory and inhibitory neural subsets. In the lowest order of approximation, interactions between excitatory and inhibitory subsets appear as linear elastic coupling, while those within and between excitatory and excitatory subsets appear as nonlinear frictional coupling. The case of two coupled oscillators is investigated by the method of averaging and the stability conditions for two mode oscillations are obtained. Internal resonance is also discussed briefly in the case of identical oscillators.     

Phase response characteristics of model neurons determine which patterns are expressed in a ring circuit model of gait generation

In order to assess the relative contributions to pattern-generation of the intrinsic properties of individual neurons and of their connectivity, we examined a ring circuit composed of four complex physiologically based oscillators. This circuit produced patterns that correspond to several quadrupedal gaits, including the walk, the bound, and the gallop. An analysis using the phase response curve (PRC) of an uncoupled oscillator accurately predicted all modes exhibited by this circuit and their phasic relationships – with the caveat that in certain parameter ranges, bistability in the individual oscillators added nongait patterns that were not amenable to PRC analysis, but further enriched the pattern-generating repertoire of the circuit. The key insights in the PRC analysis were that in a gait pattern, since all oscillators are entrained at the same frequency, the phase advance or delay caused by the action of each oscillator on its postsynaptic oscillator must be the same, and the sum of the normalized phase differences around the ring must equal to an integer. As suggested by several previous studies, our analysis showed that the capacity to exhibit a large number of patterns is inherent in the ring circuit configuration. In addition, our analysis revealed that the shape of the PRC for the individual oscillators determines which of the theoretically possible modes can be generated using these oscillators as circuit elements. PRCs that have a complex shape enable a circuit to produce a wider variety of patterns, and since complex neurons tend to have complex PRCs, enriching the repertoire of patterns exhibited by a circuit may be the function of some intrinsic neuronal complexity. Our analysis showed that gait transitions, or more generally, pattern transitions, in a ring circuit do not require rewiring the circuit or any changes in the strength of the connections. Instead, transitions can be achieved by using a control parameter, such as stimulus intensity, to sculpt the PRC so that it has the appropriate shape for the desired pattern(s). A transition can then be achieved simply by changing the value of the control parameter so that the first pattern either ceases to exist or loses stability, while a second pattern either comes into existence or gains stability. Our analysis illustrates the predictive value of PRCs in circuit analysis and can be extended to provide a design method for pattern-generating circuits. Received: 20 November 1996 / Accepted: 29 July 1997     

Power-rate synchronization of coupled genetic oscillators with unbounded time-varying delay

In this paper, a new synchronization problem for the collective dynamics among genetic oscillators with unbounded time-varying delay is investigated. The dynamical system under consideration consists of an array of linearly coupled identical genetic oscillators with each oscillators having unbounded time-delays. A new concept called power-rate synchronization, which is different from both the asymptotical synchronization and the exponential synchronization, is put forward to facilitate handling the unbounded time-varying delays. By using a combination of the Lyapunov functional method, matrix inequality techniques and properties of Kronecker product, we derive several sufficient conditions that ensure the coupled genetic oscillators to be power-rate synchronized. The criteria obtained in this paper are in the form of matrix inequalities. Illustrative example is presented to show the effectiveness of the obtained results.