Multivariable Harmonic Balance for Central Pattern Generators

Detecting coexisting oscillatory patterns in delay coupled Lur'e systems

This work addresses the problem of pattern analysis in networks consisting of delay-coupled identical Lur'e systems. We study a class of nonlinear systems, which, being isolated, are globally asymptotically stable. Assembling such systems into a network via time-delayed coupling may result in the change of network equilibrium stability under parameter variation in the coupling. In this work, we focus on cases where a Hopf bifurcation causes the change of stability of the network equilibrium and leads to the occurrence of oscillatory modes (patterns). Moreover, some of these patterns can co-exist for the same set of coupling parameters, which makes the analysis by means of common methods, such as the Lyapunov-Krasovskii method or the analysis of Poincaré maps, cumbersome. A numerically efficient algorithm, aiming at the computation of the oscillatory patterns occurring in such networks, is presented. Moreover, we show that our approach is able to deal with co-existing patterns, and both stable and unstable regimes can be simultaneously computed, which gives deep insight into the network dynamics. In order to illustrate the efficiency of the method, we present two examples in which the instability of the network equilibria is caused by a subcritical and a supercritical Hopf bifurcation. In addition, a bifurcation analysis of the subcritical case is performed in order to further explain the occurrence of the detected coexisting modes.

Intercellular delay regulates the collective period of repressively coupled gene regulatory oscillator networks

Most biological rhythms are generated by a population of cellular oscillators coupled through intercellular signaling. Recent experimental evidence shows that the collective period may differ significantly from the autonomous period in the presence of intercellular delays. The phenomenon has been investigated using delay-coupled phase oscillators, but the proposed phase model contains no direct biological mechanism, which may weaken the model's reliability in unraveling biophysical principles. Based on a published gene regulatory oscillator model, we analyze the collective period of delay-coupled biological oscillators using the multivariable harmonic balance technique. We prove that, in contradiction to the common intuition that the collective period increases linearly with the coupling delay, the collective period turns out to be a periodic function of the intercellular delay. More surprisingly, the collective period may even decrease with the intercellular delay when the delay resides in certain regions. The collective period is given in a closed-form in terms of biochemical reaction constants and thus provides biological insights as well as guidance in synthetic-biological-oscillator design. Simulation results are given based on a segmentation clock model to confirm the theoretical predictions.

Analysis of a neural oscillator

Although the Matsuoka neural oscillator, which was originally proposed as a model of central pattern generators, has widely been used for various robots performing rhythmic movements, its characteristics are not clearly explained even now. This article shows two closed-form relations that express the frequency and amplitude of the generated oscillation as functions of the parameters of the model. Although they are derived based on a rough linear approximation, they accord with the result obtained by a simulation considerably. The obtained relations also give us some nontrivial predictions about the properties of the oscillator.

Multivariable harmonic balance analysis of the neuronal oscillator for leech swimming

Biological systems, and particularly neuronal circuits, embody a very high level of complexity. Mathematical modeling is therefore essential for understanding how large sets of neurons with complex multiple interconnections work as a functional system. With the increase in computing power, it is now possible to numerically integrate a model with many variables to simulate behavior. However, such analysis can be time-consuming and may not reveal the mechanisms underlying the observed phenomena. An alternative, complementary approach is mathematical analysis, which can demonstrate direct and explicit relationships between a property of interest and system parameters. This paper introduces a mathematical tool for analyzing neuronal oscillator circuits based on multivariable harmonic balance (MHB). The tool is applied to a model of the central pattern generator (CPG) for leech swimming, which comprises a chain of weakly coupled segmental oscillators. The results demonstrate the effectiveness of the MHB method and provide analytical explanations for some CPG properties. In particular, the intersegmental phase lag is estimated to be the sum of a nominal value and a perturbation, where the former depends on the structure and span of the neuronal connections and the latter is roughly proportional to the period gradient, communication delay, and the reciprocal of the intersegmental coupling strength.