^{1}

^{1}

^{*}

A transient delay-coupled network was proposed by modifying the standard delay-coupled network with a transient coupling technique to enlarge the synchronization domain of the network, where the synchronization domain is the interval of the coupling strength for which the network gets synchronized. The coupling of the transient delay-coupled network is activated when the systems are in a particular region (coupling region) of the phase space and inactivated otherwise, which is different from the standard coupling. The specific synchronization performance of the transient delay-coupled network was investigated through case studies. The relationships between the synchronization domain and the coupling region were obtained by gauging the synchronization index. It is understood that the synchronization domain changes in a non-smooth manner with the variation of the coupling region. In particular, the synchronization domain of a transient delay-coupled network is much larger than that of the standard delay-coupled network when the coupling region is appropriately determined.

Synchronization is a typical dynamic phenomenon in coupled nonlinear systems [

Synchronization of networks with time-varying topologies has drawn considerable attention, as researchers have sought to understand the complex collective behaviors of real networks with time-varying structures [

The transient coupling is one of the simple types of switching structures proposed to enhance the synchronization of driving-driven systems [

Connection delays are generally incorporated into coupled systems, particularly in large-scale networks and long-distance secure communication networks. Long connection delays result from finite information transmission as well as processing speeds [

It is fundamentally necessary to enhance the synchronization performance of delay-coupled networks with applicable methods. To the best of the author’s knowledge, transient delay-coupling has not been introduced into delay-coupled networks, and the effects of the transient delay-coupling on the synchronization of delay-coupled networks have not been investigated. Inspired by this, the transient delay-coupling is introduced into the delay-coupled network to enhance the synchronization performance in this paper, which leads to the transient delay-coupled network. Its specific synchronization features are studied through case studies. The rest of the article is organized as follows: In Section 2, transient delay-coupling is introduced into the delay-coupled network. In Section 3, the specific synchronization performances of the transient delay-coupled network are studied. In the last section, conclusions are drawn from the above investigation.

A fully coupled network with n identical nodes and a connection delay is described by

X ˙ i = F ( X i ) + g n − 1 ⋅ ( a i i H ( X i ) + ∑ j = 1 , j ≠ 1 n H ( X j ( t − τ ) ) ) , i = 1 , 2 , ⋯ , n , (1)

where X ˙ i = F ( X i ) describes the dynamics of the ith node, X i = ( x i 1 , x i 2 , ⋯ , x i m ) T ∈ R m is a vector of the state variables, F ( ⋅ ) = ( f 1 ( ⋅ ) , f 2 ( ⋅ ) , ⋯ , f m ( ⋅ ) ) T ∈ R m is a vector of nonlinear functions, g is the total coupling strength of each node, τ represents the connection delay, H ( ⋅ ) = ( h 1 ( ⋅ ) , h 2 ( ⋅ ) , ⋯ , h m ( ⋅ ) ) T ∈ R m is the vector of inner coupling functions, and a i j ( i , j = 1 , 2 , ⋯ , n ) are the coupling coefficients with and a i j = − ∑ j = 1 , j ≠ i n a i j .

In the complete synchronization manifold, the dynamics of the synchronized state X 1 = X 2 = ⋯ = X n = X is determined by

X ˙ = F ( X ) + g ( H ( x − τ ) − H ( x ) ) . (2)

For a fully coupled network, all the perturbed equations on the transverse directions are the same, described by the master stability function [

η ˙ = D F η − g ( D H η + 1 n − 1 D H τ η ( t − τ ) ) , (3)

where η represents the deviation in the transverse directions of the synchronization manifold, DF is the Jacobian matrix of F ( ⋅ ) , evaluated on the synchronized state X ( t ) , and DH and D H τ are the Jacobian matrix of H ( ⋅ ) evaluated at the synchronized state X ( t ) and delayed state X ( t − τ ) , respectively. The synchronization manifold is stable, and synchronization is attained when the maximum conditional Lyapunov exponent determined by the master stability function (3) is negative [

The evolution of the master stability function (3) depends on both the present and past states [

λ η t 0 ( θ ) = lim t → ∞ 1 t ⋅ ln ( ‖ η t 0 + t ( θ ) ‖ C ‖ η t 0 ( θ ) ‖ C ) , (4)

where ‖ ϕ ( θ ) ‖ C = max θ ∈ [ − T , 0 ] ‖ ϕ ( θ ) ‖ L 2 is the norm of the function ϕ ( θ ) defined in Banach space C, and ‖ ⋅ ‖ L 2 denotes the Euclidean norm defined in R m .

Let t i = t 0 + i h ( i = 0 , 1 , ⋯ , + ∞ ) denote the time series with the step size h, the conditional Lyapunov exponent λ η t 0 ( θ ) is then described as

λ η t 0 ( θ ) = lim i → ∞ 1 i h ⋅ ln ( ‖ η t i ( θ ) ‖ C ‖ η t 0 ( θ ) ‖ C ) = lim i → ∞ 1 i h ⋅ ln ( ∏ k = 1 i ‖ η t k ( θ ) ‖ C ‖ η t k − 1 ( θ ) ‖ C ) = lim i → ∞ 1 i h ∑ k = 1 i ln ( ‖ η t k ( θ ) ‖ C ‖ η t k − 1 ( θ ) ‖ C ) = lim i → ∞ 1 i h ∑ k = 1 i l k

l k = ln ( ‖ η t k ( θ ) ‖ C ‖ η t k − 1 ( θ ) ‖ C ) is regarded as the local conditional Lyapunov exponent corresponding to the point X ( t k ) of the synchronized state [

The coupling without a connection delay is not invasive, while the coupling force vanishes when synchronization is attained. However, coupling with a connection delay is invasive, and the coupling force is still acting on the synchronized system when the synchronization is achieved. For delay-coupled systems, the conditions for complete synchronization are strict. One of the conditions is that all the coupled nodes have the same input degree. Thus, the transient coupling scheme proposed for coupled systems without connection delays [

In this article, the transient delay-coupling is introduced into the network given by Equation (1) by letting

a i j ( i ≠ j ) s = { a i j ( i ≠ j ) X i ∈ Ω and X i ∈ Ω 0 X i ∉ Ω and X i ∉ Ω , (5)

where Ω ⊆ R m is the coupling region. Equation (1) becomes

#Math_44# (6)

where, a i i s = − ∑ j = 1 , j ≠ i n a i j s . The ith node will cut off all the connections with other nodes when X i ∉ Ω , and all the other coupled nodes still have the same input degree. When Ω = R m , the transient delay-coupling in Equation (5) degrades into static delay-coupling with as a i i s = a i j .

The synchronization performance of the network given by Equation (6) is closely associated with the coupling region Ω . It is challenging to determine the proper coupling regions analytically such that the network can get synchronization in a large synchronization domain because of the strong nonlinear relationship between the synchronization and the coupling regions [

This section examines the specific synchronization features of the transient delay-coupled networks through case studies. To show the advantage of the transient delay-coupling over the standard delay-coupling, we compare the synchronization performance of the transient delay-coupled network and that of the standard delay-coupled network.

Case 1: A fully delay-coupled network with nodes of the Hindmarsh-Rose (HR) neuron coupled via the x component only is considered, as follows:

{ x ˙ i = y i − a x i 3 + b x i 2 − z i + I + g n − 1 ⋅ ∑ j = 1 , j ≠ i n a i j s ( x j ( t − τ ) − x i ) y ˙ i = c − d x i 2 − y i z ˙ i = r [ s ( x i + x 0 ) − z i ] (7)

where x i , y i and z i are the state variables of the ith node, representing the membrane potential, exchange of ions across the neuron membrane, and a slowly changing adaptation current, respectively; I is the membrane input current for biological neurons, r is usually a small parameter related to the slow variable z i , and x 0 represents the resting potential of the system [

The HR system is the minimal model that can reproduce several electrical firing activities observed in real biological neuron activities. HR neuron networks have been studied intensively to understand the collective dynamics of real neural networks [

The collective dynamics are evaluated by the synchronization index defined as

S d e x = ( n − 1 ) − 1 ∑ i = 1 n − 1 〈 ‖ X i + 1 ( t ) − X i ( t ) ‖ L 2 〉 . (8)

where, X i ( t ) = ( x i ( t ) , y i ( t ) , z i ( t ) ) T is the vector of state variables of ith node, 〈 ⋅ 〉 represents the average over a sufficiently long time after the transient process. Smaller synchronization index S d e x corresponds to better collective dynamics, and S d e x = 0 indicates complete synchronization.

As the synchronization performance of the transient delay-coupled network (7) is closely related to the coupling region Ω , a proper coupling region Ω should first be established. There are different ways to change the coupling signal with different types of coupling regions. For the sake of simplifying the study, consider one of the simplest types of coupling regions:

Ω S l = { ( x , y , z ) | x ≤ S l } . (9)

As shown in

Without loss of generality, the synchronization performances of the network provided by Equation (7) with network scale n = 10 and different connection delays τ were studied for g ∈ [ 0 , 1.5 ] . The synchronization domain is the interval of the coupling strength g for which the network gets synchronized. The relationships between the synchronization domain and the coupling region ( S l ) are obtained in

To demonstrate the advantage of the transient delay-coupling over the standard delay-coupling, comparisons between the control performance of the transient delay-coupled network with a proper controlling region Ω S l and that of the standard delay-coupled network are given in

Case 2: The transient delay-coupled Rössler systems connected linearly via the x component only is considered, as follows:

{ x ˙ i = − ( y i + z i ) + g n − 1 ⋅ ∑ j = 1 , j ≠ i n a i j s ( x j ( t − τ ) − x i ) y ˙ i = x i + a y i z ˙ i = b + z i ( x i − c ) (10)

The dynamics of coupled Rössler systems have been studied intensively to understand the collective phenomena of coupled chaotic systems [

The collective dynamics of the transient delay-coupled network given by Equation (10) are again evaluated by the synchronization index defined in Equation (8), and the simplest type of coupling region Ω S l described in Equation (9) is used to construct the transient delay-coupling in the network. The synchronization performances of the network with network scale n = 10 and different

connection delays τ are examined for g ∈ [ 0 , 10 ] . The relationships between the synchronization domain and the controlling region are obtained in

To instantiate the advantage of the transient delay-coupling over the standard delay-coupling, comparisons between the control performance of the transient delay-coupled network with a proper controlling region Ω S l and that of the standard delay-coupled network are given in

Summarizing the conclusions obtained from the above case studies, the synchronization performance of the transient delay-coupled network is closely related to the coupling region. The synchronization domain changes in a non-smooth fashion with the variation of the controlling region. In particular, the transient delay-coupling proposed in this paper is beneficial to the synchronization performance of the delay-coupled network; the synchronization domain of the transient delay-coupled network is much larger than that of the standard delay-coupled network when the coupling region is appropriately determined.

Delay-coupling is invasive as the delay-coupling force still acts on the synchronized system when the synchronization is achieved. The connection delay destabilized the synchronization manifold frequently and deteriorated the synchronization performance of the delay-coupled network. It is of fundamental importance to improve the synchronization performance of the delay-coupled network. The transient delay-coupled network was proposed, for the first time, by applying the transient coupling technique to the standard delay-coupled network. Case studies have indicated that the synchronization performance of the delay-coupled network is much better than that of the standard delay-coupled network when the coupling region is appropriately chosen. From the practical application point of view in engineering, the transient delay-coupled network has some advantages over the standard delay-coupled network. For example, in some real networks, the coupling strength is restricted to a particular limited domain of the parameter space. In this case, the synchronization domain of the transient delay-coupled network can be changed and enlarged to make the coupling strength fall within the synchronization domain by choosing an appropriate coupling region, and thus, the synchronization is attained. Different than the continuous coupling force of the standard delay-coupled network, the coupling force of the transient delay-coupled network is only activated in the coupling region. Consequently, the transient delay-coupled network is low-cost and conserves energy, especially when the coupling region is small. Besides, the results obtained in the present paper can help to understand the synchronization of real delay-coupled networks with time-varying structures.

This paper was supported by the NSF of China under Grants 11872197, 12072140, and the JXNSF of China under Grant 20192BAB202002.

The authors declare no conflicts of interest regarding the publication of this paper.

Zhang, Y.Y. and Zheng, Y.G. (2021) Synchronization of Transient Delay-Coupled Network. Journal of Computer and Communications, 9, 25-37. https://doi.org/10.4236/jcc.2021.99002