#### Abstract

A network of coupled limit cycle oscillators with delayed interactions is considered. The parameters characterizing the oscillator’s frequency and limit cycle are allowed to self-adapt. Adaptation is due to time-delayed state variables that mutually interact via a network. The self-adaptive mechanisms ultimately drive all coupled oscillators to a consensual cyclostationary state, where the values of the parameters are identical for all local systems. They are analytically expressible. The interplay between the spectral properties of the coupling matrix and the time delays determines the conditions for which convergence towards a consensual state takes place. Once reached, this consensual state subsists even if interactions are removed. In our class of models, the consensual values of the parameters depend neither on the delays nor on the network’s topologies.

*The farther back you can look, the farther forward you are likely to see*

Winston Churchill

#### 1. Introduction

The harmonic excitation of an elementary-damped harmonic oscillator
with , , and produces the well-known asymptotic response (c.f. [1]):
where and depend on the control parameters , * environment* here materialized by the input is totally insensitive to the oscillator , implying that in (1) is slaved by the external forcing. The next stage of complexity is to replace the harmonic oscillator in (1) by a Lienard system:
where is a nonlinear controller. In absence of external excitation in (3) (i.e., when ), we assume to asymptotically drive the orbits towards a stable limit cycle which is independent of the initial conditions—the paradigmatic illustration being here the Van der Pol oscillator. When and for a suitably selected range of parameters, the time asymptotic response of (3) can be qualitatively written as (c.f. [2, 3])
with being a* synchronized signal* with the same periodicity as the environment (i.e., ). By construction, the external forcing in (3) is, as before, insensitive to the Lienard oscillator. In the resulting synchronized regime, the oscillator is *caught* by the external excitation—in other words, the system * adjusts* itself to the environment but the environment remains insensible to the system. Observe that the dynamical response given by (4) only subsists as long as acts on the system. That is, as soon as the environment effect is removed (i.e., in (3)), the system (i.e., the limit cycle oscillator), after a transient time, recovers its original behavior—converges towards its limit cycle.

In our present paper, we will extend the previous classical * system-environment* relationship in order to allow more realistic situations where mutual interactions *permanently* affect both the * system* and the * environment*. Such * adaptive* mechanisms can modify individual dynamics on a *permanent* basis. To stylize this new situation, the basic dynamics given by (3) is modified as
where * mutual adaptation* of the * system-environment* dipole. Hence, the system of (5) has now to be considered globally—the * system* and its * environment*, are allowed to adaptively coevolve.

The particular case of (5) when (i.e., for all ) has been studied in [4, 5]. This type of dynamical system provides a cornerstone of bioinspired robotics where legs (or arms) of robots may be modeled by damped oscillators as (1) (i.e., ). To ensure a maximum leg stride, the damped oscillators must be excited by at the damped oscillator's resonant frequency (i.e.,

The case where neither nor is trivial has been covered in [7, 8], where the authors not only considered two adaptive coupled limit cycle oscillators, but mutually interacting through a network. Here, self-adapting oscillators can be applied to robot formation modeling. Each individual robot belonging to a swarm, circulating around a specific point, adapts its angular velocity in order to lower the amount of exchanged information to maintain the formation. Self-adaptation in networks is also considered in [9], where here the control signals (and not the local systems) of the variables of the CPG adapt, and this, to quickly react to new situations and produce several different behavioral patterns.

Building on what has been done in previous contributions [4, 5, 7, 8], we here consider a network of limit cycle oscillators interaction with time-delayed state variables. The general form of our dynamical system in the phase-radius coordinates (i.e., polar coordinates and ) is
where and govern the local dynamics, and are the state variables, is a parameter set determining the local characteristics and are coupling strengths. The coupling dynamics is the gradient of a potential

This paper is organized as follows: in Section 2 we define the three components that together form the global system. We then discuss the dynamics of our model in Section 3. An application is presented in Section 4 which is then followed by some numerical experiments in Section 5. Finally, we conclude in Section 6.

#### 2. Networks of Hopf Oscillators with Adaptive Mechanisms

We now present explicitly the local and coupling dynamics as well as the adaptive mechanisms on which we will focus.

##### 2.1. Local Dynamics

Each node of the network is equipped with a local dynamical system. In this contribution, a local system is a Hopf oscillator presented here in its polar coordinates: The state variables are and are, for the time being, fixed and constant parameters. The parameter controls the frequency of the th oscillator given by the phase dynamics . The radial dynamics produces a stable circular limit cycle with radius .

##### 2.2. Coupling Dynamics

Associated to an vertex * connected* and * undirected* network, denote by the weighted adjacency matrix with positive entries . Let be the corresponding Laplacian matrix (, where is the diagonal matrix with ). The coupling dynamics is given by the gradient of the positive semidefinite function
with and and where the matrix has entries . This matrix is positive semidefinite since all its eigenvalues are positive (i.e. nonnegative). This is a direct application of Gershgorin’s circle theorem [11]: for any eigenvalue of , there exists such that
and so . In particular, the coupling dynamics is
where are coupling strengths.

##### 2.3. Adaptive Mechanisms

In this section, we now allow the fixed and constant parameters and to (I)become time dependent, that is, , for , (II)and each of them have their own dynamics, depending solely on the state variables and , that is,

Among the numerous variants for changing the values of the parameter, we focus on those presented in [7, 8, 12], that is,
where are the entries of and

#### 3. Network's Dynamical System with Delay

We now discuss the resulting dynamics in the presence of a time delay affecting both the coupling dynamics and the adaptive mechanisms. We hence consider for . For (14), we have the following.

*Two Constants of Motion*. The functions

are constants of motion: in other words, if and are orbits of (14), then

*Existence of a Consensual Oscillatory State. *We can explicitly exhibit a consensual oscillatory state. Indeed, for given and ,

for all

Observe that in absence of the radial component and without adaptation, (14) yields the famous Kuramoto model with delays [13, 14]. In these contributions, the authors considered the coupling dynamics with delay of the form:

In the absence of time delay (i.e.,

*Convergence towards a Consensual Oscillatory State.* The Limit (18) raises two issues: (1) the existence itself and (2) the limit values and . For expository reasons, we first discuss the limit values and then the convergence conditions.

*Limit Values.* Thanks to the constant of motions in (15), we have
with and orbits of (14) with given initial conditions. Supposing that Limit (18) holds, we hence have
and so the asymptotic values are analytically expressed as
It is important to emphasize that the consensual values and do not depend on the network topology (i.e., not on ), nor on the initial conditions of the state variables (i.e., not on nor on the time delay (i.e., not on

*Convergence Conditions.* To this aim we study the first-order approximation of (14) in the vicinity of Solution (17) and assume that linear stability analysis is sufficient to infer convergence conditions for the nonlinear system. Accordingly, we study the asymptotic behavior of the small perturbations , and and write
Taking into account the constant of motions, we impose that

*First Order Approximation.* Rearranging the variables (i.e., the first are the , the second are the , the third are the , and finally the last are the ), the first-order approximation of (14) is
with the identity matrix , diagonal matrices

*Diagonalization.* Suppose now that

Changing the basis of System (24) with a bloc matrix (each bloc of size ) with

with

The case is worked out in the appendix. For , let us focus on the -dimensional systems and we rewrite (25a)-(25b) as linear second-order time delayed differential equations The convergence towards a consensual state is hence determined by the asymptotic stability of the zero solution of (26). Stability follows if, and only if, all roots of the corresponding characteristic equations have strictly negative real parts (c.f. [15] for details). For (25a), one can apply Theorem 3.3 in [16] which states that in this case the zero solution is asymptotically stable if and only if(I)(II)for . We emphasize that the consensual value influences the condition for convergence whereas it is not the case for . This leads to the idea that shaping the attractor is more delicate than tuning the angular velocity. This has been observed in [17].

*Adaptation only on *. If there is no adaptation on the radii (i.e.,

*No Time Delay in the Coupling Dynamics.* If there is no time delay in the coupling dynamics (i.e., the coupling dynamics is defined as in (11) with no delay), then (26) becomes, for , Invoking Theorem 3.5 in [16], the zero solution for (30a) and (30b) is asymptotically stable if and only if.

For (30a)

for (30b) for .

*Summary.* For a network (with arbitrary topology but with symmetric, positive entries adjacency matrix) of Hopf oscillators (as defined in Section 2.1 by (8)) interacting through time delayed Kuramoto type coupling (as defined in Section 2.2 by (11)) and with time-delay adaptive mechanisms (as defined in Section 2.3 by (13)) on the frequencies and amplitudes of the local systems—in other words, for (14), we have the following. (I)two constants of motions (c.f. (15)), (II)the existence of a consensual oscillatory state (c.f. (17)),(III)the consensual oscillatory state is linearly (i.e., locally) stable if all the roots of characteristic equations corresponding to (26) have strictly negative real parts, (IV)if there is no adaptation on the radii, the consensual oscillatory state is linearly (i.e., locally) stable if (27a), (27b), and (29) hold, (V)if there is no delay in the coupling dynamics, the consensual oscillatory state is linearly (i.e., locally) stable if (31), and (32) hold.

##### 3.1. Miscellaneous Remark: Delayed Stabilization Mechanism

In this section we discuss the particular case arising when the delay is introduced only in the stabilization mechanism (i.e., dissipative part) of the local dynamics. The dynamical system is

Equations (33) still admit the existence of a consensual oscillatory state and two constants of motion as in (15) and in (17), respectively. Linear stability analysis of the consensual state reduces to the study of for . The zero solution for (34a) is asymptotically stable. For (34b), we apply Theorem in [16] which states that in this case the zero solution is asymptotically stable if and only if(I)(II) for .

#### 4. Applications

Conceptually, the problem of reliably distributing time and frequency among several spatial remote locations is a * “leitmotiv”* in applications ranging from basic metrology, navigation and position determination, signal processing, computer communications, energy distribution networks, swarms robotics, bioengineering, multiagents systems, life sciences, acoustics, and musical art to give but only a highly nonexhaustive list. Presently, a strong research impetus is devoted to complex interacting oscillating systems able to exhibit self-adaptive capabilities leading to a resilient consensual dynamic. Whatever the configurations under study, communication delays between the collection of interacting subparts of the global systems are physically unavoidable. Depending on the underlying time scales, delays do strongly affect the resulting dynamics. Our class of models explicitly study the influence of delays and in particular their destabilizing effects, that modify the instantaneous behavior. By an appropriate tuning of control parameters (e.g., susceptibility constants), our class of models offer, via a unique formalism, the possibility to continuously explore interacting configurations ranging from slave-master (i.e., system-environment relationship) to fully decentralized regimes.

Alternatively, we may view this problematic in the context of * soft-controlled systems* which presently receive a sustain attention [19]. Here, a swarm of agents is infiltrated by a * lure agent* (sometimes called a shill in economy). While the lure exhibits all the features of any ordinary agent, it can be externally controlled by an operator. As the interactions between the lure and any agent of the swarm remain unaffected (i.e., the lure remains incognito to ordinary agents), the external control of the lure can ultimately drive the whole population to a specific configuration. In our class of dynamics, a suitable choice of the susceptibility constant of a given local system (i.e., oscillator) may convert it into a shill. Indeed, in view of (21), the ultimate consensual values and are weighted averages. Such weighted averages can be made to be strongly dependent on a very insensitive shill-stubborn to any external influence (i.e., with a very low susceptibility constant).

In absence of time delays, convergence towards a consensual state is observed even for large heterogeneities (widely dispersed initial frequencies and radii and large discrepancies among the susceptibility constants). Hence, a shill agent can be easily introduced. However, our present study shows, that time delays restrict the conditions for convergence towards a consensual state. As a consequence, the implementation of a lure is more delicate matter (i.e., here, the dynamics is far more sensitive to the value taken by the susceptibility constants).

#### 5. Numerical Simulations

*Adaptation on ** and **.* We consider four Hopf oscillators interacting on a network with topology as in Figure 1. We choose the coupling strengths and susceptibility constants as

The resulting dynamics is shown in Figure 2. With the same initial conditions, we carry out another numerical simulation with here

**(a)**

**(b)**

**(a)**

**(b)**

Note that in Figure 2(a), the converge close to , that is, close to the initial value . This because the first oscillator’s susceptibility is “small” and hence it is this local system that acts as a shill. It interacts with its neighbor in the same way as they act with it. It can control the behavior of the network and this without being connected to all other local systems.

In the extreme case, when *here adaptation enhances the syntonization capability* of the network.

*Adaptation only on *. Two Hopf oscillators, both having the same radius for the attractor (i.e.,

**(a)**

**(b)**

#### 6. Conclusion

It is more the rule than the exception that parameter adaptation in dynamical systems can be achieved via delayed mechanisms. This therefore converts ordinary differential equations arising in absence of delay to * functional differential equations*. Stability issues become more difficult to discuss since, dealing with functional differential equations, an infinite number of degrees of freedom is introduced into the dynamics. For the class of oscillatory networks with parametric adaptation we here considered, we are able to observe how the time delay affects the adaptation mechanisms. While it is intuitively expected, that large delays are likely to destabilize the dynamics, we are here able to analytically quantify the underlying critical delays. The analytical linear stability discussion is made possible since, for our class of dynamics, stability issues can be reduced to the study of two linear second order functional differential equations for which suitable theorems can be found. Finally, we emphasize that numerical simulations show that adaptation may enhance the emergence of common dynamical pattern whereas in classical synchronization, time delays may be too large for synchronous motion to be attained.

#### Appendix

Assuming that all perturbations and for converge to zero, let us now study the case for . Here, and so
and therefore and for all . Both of these constants and are zero. This is because the first orthonormal base vector (i.e., the normalized eigenvector for the eigenvalue ) is

#### Acknowledgments

M.-O. Hongler acknowledges partial support from the ESF (European Science Foundation) under the project entitled Exploring the Physics of Small Devices. J. Rodriguez acknowledges the support from the DFG-IRTG 1132 (Deutsche Forschungsgemeinschaft—International Research Training Group) under the project entitled Internationales Graduiertenkolleg—Stochastics and Real World Models.