Synchronizability of multilayer star-ring networks with variable coupling strength

: We investigate the synchronizability of multilayer star-ring networks. Two types of multilayer networks, including aggregated coupling and divergent coupling, are established based on the connections between the hub node and the leaf nodes in the subnetwork. The eigenvalue spectrum of the two types of multilayer networks is strictly derived, and the correlation between topological parameters and synchronizability is analyzed by the master stability function framework. Moreover, the variable coupling strength has been investigated, revealing that it is significantly related to the synchronizability of the aggregated coupling while having no influence on the divergent coupling. Furthermore, the validity of the synchronizability analysis is obtained by implementing adaptive control on the multilayer star-ring networks previously mentioned. Calculations and comparisons show that the differences caused by the sizes of multilayer networks and interlayer coupling strength are not negligible. Finally, numerical examples are also provided to validate the effectiveness of the theoretical analysis.


Introduction
Complex networks are extensively explored in the field of nonlinear dynamics [1][2][3].Most of the previous studies of networks have focused on single-layer networks, while many real-world networks cannot be described specifically by single-layer networks alone.With advancements in theory and practice, it has gradually been found that networks do not isolate entities but rather interact with other networks in structure and function.For example, transportation networks [4] are composed of aviation networks, shipping networks and railway networks.The biological networks [5] consist of gene regulatory networks, metabolic networks and protein interaction networks.Therefore, it is imperative to study the multilayer networks formed by the interaction between each single-layer network in depth.It can be abstracted into a multilayer network model with intralayer and interlayer coupling to better interpret the association between sub-networks in the complex relationship.Recent research on multilayer networks has attracted some attention.Zhang et al. [6] studied the adaptive explosive synchronization of multilayer networks without considering the specific correlation between the oscillators and the natural frequency.Wang et al. [7] considered the nodes with Lur'e-type dynamics and investigated how to achieve synchronization from node to whole network in a directional topology.A series of studies have also been carried out in the field of topology identification for multilayer networks.Liu et al. [8] realized topology identification of the networks based on adaptive observers of chaotic external signals without complex linear independence assumptions.Different from previous research strategies, Mei et al. [9] studied the structure identification problem based on the compressed sensing method.In addition, in the study of multilayer network dynamics, Wang et al. [10] investigated the unique superdiffusion phenomenon caused by the directionality of connections.Turalska et al. [11] proved that the correlation existing in the networks would affect the dynamic cascade process and could prevent the spread of failures between layers.Therefore, it is imperative to undertake further exploration of multilayer networks.
Synchronization is a typical topic for the numerous dynamic behaviors observed in complex networks.It plays a significant role in various domains such as multi-agent cooperation, biological oscillators, signal transmission and numerous other applications [12][13][14][15].In 1990, Pecora and Carroll [16] first proposed synchronization control for two coupled chaotic oscillators.Subsequently, the synchronization problems of single-layer complex networks are studied one after another.Wu et al. [17] investigated the phenomenon of projective synchronization in drive-response networks involving complex-variable chaotic systems.Hong et al. [18] researched the phenomenon of collective synchronization in a system of coupled oscillators that were organized on small-world networks.Li et al. [19] employed linear matrix inequalities to examine the synchronization issue in complex dynamical networks with time-varying delays.Due to the complexity of multilayer networks, some methods of synchronization in single-layer networks are no longer applicable.Currently, some advances have been achieved in the field of synchronization research in multilayer networks.Zhao et al. [20] focused on the synchronization of multiplex delayed networks with stochastic perturbations via pinning control, which significantly improved the synchronization efficiency.Fan et al. [21] delved into the study of adaptive event-triggered prescribed performance learning synchronization for complex dynamical networks.The above studies focused on the synchronization of chaotic systems and networks with given topological parameters and coupling strengths.However, the relationship between topological parameters, the size of parameters and coupling strengths changing with time to synchronizability analysis is not clear, which is a challenging task in multilayer networks that deserves further research.
In large-scale complex systems, individual behavior tends to be consistent, which is an external manifestation of synchronizability, and it is very important to understand its influence from the perspective of topology.A supra-Laplacian matrix based on the diffusion dynamics equation of multilayer networks is proposed to analyze its eigenvalue spectrum for studying complex network synchronizability [22][23][24].So far, there have been some initial findings in the study of synchronizability in multilayer complex networks.Aguirre et al. [25] investigated the impact of node degree on the synchronization ability of two interconnected networks.They demonstrated that linking high-degree nodes in each network is the most effective approach for achieving synchronization.The two-layer star networks were derived strictly from the eigenvalue spectrum and analyzed for synchronizability by Xu et al. [26].Li [27] proposed a two-layer dumbbell network model to analyze the synchronizability and verified the synchronization through a simple numerical simulation.Deng et al. [28] researched the synchronizability of multilayer chain networks.However, in real-world complex networks, whether they are single-layer or multilayer, studying the synchronizability of networks with weighted, directed and variable coupling strength is more meaningful.Variable coupling strength plays an important role in information transmission and complex system control.Adjusting its size can effectively optimize traffic flow regulation and signal coordination in the intelligent transportation system [29].In the Internet of Things, its regulation can optimize the transmission and processing of information, enabling intelligent detection and control [30].Given its extensive application in engineering, it cannot be ignored in synchronization analysis.Therefore, apart from considering the influence of topological parameters, we also focus on analyzing the influence of variable coupling strengths on the synchronizability of the network.The common control methods in complex networks include adaptive control, sliding mode control, impulsive control and so on [31,32].Compared with other methods, the adaptive control strategy selected in this paper has the advantages of adjusting changes independently with network uncertainties, improving the efficiency of synchronization control, and having stronger robustness [33].Shi et al. [34] realized finite/fixed-time synchronization of complex networks with uncertain internal coupling based on the quantitative adaptive control method.Qin et al. [35] used an adaptive control strategy to study the robust H∞ synchronization problem of complex networks with multiple delays.
Inspired by the previous literature, we will examine two types of directed multilayer variable coupling strength star-ring networks: divergent coupling and aggregated coupling.In the case of divergent coupling, the hub node establishes connections directed toward the leaf nodes, thus forming a unidirectional communication flow, which plays an important role in improving the reliability and resource coordination of the sensor system [36].Conversely, aggregated coupling involves the leaf nodes directing connections toward the hub nodes, with advantages such as multi-source data aggregation and resource optimization [37,38].The major contributions of this paper can be summarized as follows: (i) The eigenvalue spectrum of the multilayer networks with divergent and aggregated coupling is strictly derived.
(ii) The correlation between topological parameters and synchronizability is analyzed with bounded and unbounded synchronous regions.
(iii) The variable coupling strength in divergent and aggregated coupling star-ring networks has been investigated, which is only related to the synchronizability of the aggregated coupling mode.
The remaining sections of the paper are organized as follows.In Section 2, several preliminaries are provided to establish a foundation for the subsequent discussions.The synchronizability of multilayer divergent and aggregated directionally coupled star-ring networks under varying coupling is illustrated in detail in Sections 3 and 4, respectively.To validate the accuracy of our synchronizability analysis, synchronization controls are provided in Section 5. Finally, in Section 6, we summarize our conclusions and provide a brief discussion of future work.

Preliminaries
Consider a complex multilayer network with M layers and N interconnected nodes in each layer.The evolution of a full network dynamical equation is described as follows [39]: where xi K ∈ N denotes the state of the ith node in the Kth subnetwork, and f(xi K ): where wij K is defined as follows: if there exists a connection between node i and j(i ≠ j) within the Kth layer, then wij K = wji K = 1; otherwise, wij K = wji K = 0.In addition, the diagonal elements of W K satisfy the following condition: .
Therefore, L K = -c K W K is a Laplacian matrix of the Kth layer.Similarly, the interlayer weight di KZ = di ZK = 1 if node i in the Kth layer is connected to node j in the Zth layer.Otherwise, di KZ , .
Then, D = (di KZ )∈ M×M is the negative Laplacian matrix of the interlayer.Let  L be the supra-Laplacian matrix that represents the interlayer topology and  I be the supra-Laplacian matrix that represents the intralayer topology.Then, the supra-Laplacian matrix of the above system is described as [40]: Let L I denote the Laplacian matrix of the interlayer networks.Specifically, the interlayer supra-Laplacian matrix is where ⊗ is the Kronecker product and IN represents the N × N identity matrix.The intralayer supra-Laplacian matrix of a multiplex network is defined as: where K = 1, 2, …, M, and L K is the Laplacian matrix of the Kth layer.
In accordance with the master stability function framework [41], we can analyze whether the eigenvalues of the multilayer network fall in the corresponding synchronous regions to judge the synchronizability of the network.Both bounded and unbounded regions are usually the two main categories for discussion.If the network has an unbounded synchronous region, the minimum nonzero eigenvalue λ2 of is key to research the synchronizability.Moreover, synchronizability is positively correlated with λ2.On the other hand, if the network has a bounded synchronous region, the synchronizability depends on the ratio R = λmax /λ2, which is the maximum eigenvalue λmax to the minimum nonzero eigenvalue λ2 of .At this time, synchronizability is negatively correlated with R.
For the necessity of subsequent theoretical deduction, a lemma and two hypotheses are presented here: Lemma 2.1.[42] A and B are N × N matrices, M is the number of layers, then Hypothesis 2.1.[43] Suppose that there exists α > 0 where, ei K (t) and s(t) represent the synchronization errors and target system of the network, respectively.
Hypothesis 2.2.[44] For any p, q n , there exists a positive definite matrix Q n×n , and the following always holds:

Synchronizability of multilayer network with divergent coupling
Based on the theory of Section 2, a multilayer star-ring network with divergent coupling is constructed, and its characteristic spectrum is strictly derived.Furthermore, when the synchronous region is unbounded or bounded, the influence of different topological parameters on synchronizability is studied in this section.

The structure of the multilayer network with divergent coupling
In this section, the synchronizability analysis of the divergent directionally coupled star-ring network is considered, and its topology is shown in Figure 1.The left panel describes the intralayer connection of each subnetwork in the multilayer network, and the right panel shows the one-to-one connection of the subnetwork nodes of each layer.L (1) , L (2) , …, L (M) represent each layer.Assume that the topologies of each subnetwork are the same.There are N nodes in each subnetwork, and M represents the number of layers.The blue nodes and the red nodes represent the leaf nodes and the hub nodes, respectively.The nodes at each layer are connected one to one.It is assumed that the connections between nodes within a layer exhibit directionality, while the interlayer connections of each subnetwork are undirected.The red nodes all point one way to the blue nodes and the blue node points to the next blue node.In the aggregated coupling network, d0 represents the corresponding coupling strength between each leaf node in each subnetwork.The remaining d represents the coupling strength between each leaf node.The intralayer coupling strength of blue nodes is given by a0.Suppose that the intralayer coupling strength of the red node points to the blue nodes is denoted as a and the total intralayer coupling strength is (N-1) a at the initial time.At some time, γ (γ is even) of coupling strength becomes σa, then the coupling strength of the remaining (N-γ-1) nodes becomes (N-σγ-1) a/(N-γ-1).We have denoted (N-σγ-1) a/(N-γ-1) as δa for convenience.

The eigenvalue spectrum and synchronizability
The supra-Laplacian matrix of the network can be written as: .
, and According to Lemma 2.1, the corresponding characteristic polynomial of  is expressed as follows:


From the characteristic polynomial above, whether the topological parameter N is odd or even affects the eigenvalues, which will be discussed below.
Table 1.λ2 and R change with N, a, a0, d, d0, and M when N is odd or even.
Table 2. Synchronizability changes with N, a, a0, d, d0, and M when N is odd or even.
1) N is odd, the eigenvalues  are: .
The relationships between λ2, R, and the topological parameters are depicted in Table 1.Similarly, Table 2 presents the correlations between the network's synchronizability and its topological parameters.

Comparisons and discussions
Numerical experiments between topological parameters and λ2, R of multilayer networks will be investigated.That is, including N, a, a0, d, d0, and M to discuss the synchronizability.
1) Influence of the number of subnetwork nodes N Take a0 = 0.6, d = 1, M = 15, d0 = 0.1, and a = 1 to examine the synchronizability with N (N varies from 9 to 399 if N is odd; N varies from 10 to 400 if N is even).Figure 2(a) shows that λ2 doesn't change with N in the subnetwork, that is, the synchronizability is not affected.The relationship between R and N is shown in Figure 2(b), which indicates that synchronizability is negatively correlated with N.
2) Influence of the intralayer coupling strength a Take a0 = 0.6, M = 15, d = 1, d0 = 0.1, N (odd) = 399, and N (even) = 400 to examine the synchronizability with a.As is shown in Figure 3(a), λ2 is unchanged with a, that is, the synchronizability is also unchanged with a. Then Figure 3(b) displays the relationships between R and a, which means that synchronizability is negatively correlated with a. 3) Influence of the intralayer coupling strength a0 Take a = 1, M = 15, d = 1, d0 = 0.1, N (odd) = 399, and N (odd) = 400 to explore the synchronizability with a0.Notice from Figure 4(a) that when N is even, λ2 unchanged with a0, that is, the synchronizability is not affected.However, when N is odd, λ2 exhibits a positive correlation with a0, indicating that a higher value of a0 is associated with increased synchronizability.Notice from Figure 4(b) that when N is even, R remains unchanged, so the synchronizability is not affected.However, R decreases with a0, so the synchronizability strengthens.

Synchronizability of multilayer network with aggregated coupling
In contrast to the divergent coupling, a multilayer star-ring network with aggregated coupling is built.In addition, the characteristic spectrum of the network is rigorously derived.Through theoretical deduction, we find that variable coupling strength has a significant influence on synchronizability, except for topological parameters.Furthermore, the influence of variable coupling strength on synchronizability under different function variation rules is discussed in this section.

The structure of the multilayer network with aggregated coupling
We consider an aggregated directionally coupled star-ring network as shown in Figure 8. Only, unlike the topology of Figure 1, the blue nodes all point to the red nodes within each layer.

The eigenvalue spectrum and synchronizability
The supra-Laplacian matrix of the network can be written as: .
, and According to Lemma 2.1, the corresponding characteristic polynomial  is expressed as follows:
It is obvious that the variable coupling strengths σa and δa will affect the synchronizability from the above eigenvalues.Different from Section 3, the impact of varying coupling strengths σa and δa on synchronizability will be investigated in this section.
Remark 1.Compared with the unidirectional connection between the hub node and the leaf nodes in [45], we discuss the influence of topological parameters and variable coupling strength on synchronizability through divergent and aggregated coupling.Moreover, the characteristic spectrum of multilayer divergent and aggregated coupling networks is also accurately obtained.And it is concluded that the topological parameters influencing the synchronizability of the network are different when the directionality of the node connections is changed.

Comparisons and discussions
According to the above theoretical derivation, the number of nodes in the subnetwork has no effect on synchronizability.Moreover, λ2 and R in the two cases are basically similar, so N is selected as odd for further study.To investigate the network synchronizability under variable coupling strengths σa and δa, the partial coupling strength is varied.Let the variable coupling strength follow the principles of the proportional function, the inverse proportional function, the exponential function, the logarithmic function, and the trigonometric function, respectively.The principles of variable coupling strength are as follows: σ1 = t, σ2 = 1/t, σ3 = 0.9 t , σ4 = lg(t), σ5 = sint+1.5.

1) Variable coupling strength follows the proportional function
We take σ1 = t (then δ1 = (N-tγ-1) /(N-γ-1)) for the numerical simulation.Then, λ2 = min {σ1a, δ1a}, R = λmax/λ2 = max {Md+σ1a+2a0, Md+δ1a+2a0}/ min {σ1a, δ1a}.Figure 9(a) shows the coupling strength a versus the proportional function, and the color depth corresponds to the value of λ2.If t is fixed, λ2 increases with a, and synchronizability is improved.If a is fixed, λ2 first increases and then decreases with t.Obviously, when t = 1 and a = 5, λ2 reaches its maximum and the synchronizability is the strongest.Figure 9(b) shows the coupling strength a versus the proportional function, and the color depth corresponds to the value of R. If t is fixed, R decreases with a, and synchronizability is improved.Additionally, if a is fixed, R first decreases and then increases as t increases.Obviously, when t = 1 and a = 5, R reaches its minimum, and the synchronizability is the strongest.
2) Variable coupling strength follows the inverse proportional function We take σ2 = 1/t (then δ2 = (N-(1/t) γ-1) /(N-γ-1)) for the numerical simulation.Then, λ2 = min {σ2a, δ2a}, R = λmax/λ2 = max {Md+σ2a+2a0, Md+δ2a+2a0}/ min {σ2a, δ2a}.Figure 10(a) shows the coupling strength a versus the inverse proportional function, and the color depth corresponds to the value of λ2.If t is fixed, λ2 increases with a, and synchronizability is improved.Additionally, if a is fixed, λ2 first increases and then decreases with t.Obviously, when t = 1 and a = 5, λ2 reaches its maximum, and the synchronizability is the strongest.Figure 10(b) shows the coupling strength a versus the inverse proportional function, and the color depth corresponds to the value of R. If t is fixed, R decreases with a, and synchronizability is improved.Additionally, if a is fixed, R first decreases and then increases as t increases.Obviously, when t = 1 and a = 5, R reaches its minimum, and the synchronizability is the strongest.

Synchronization control
A detailed and rigorous theoretical derivation is presented for the eigenvalue spectrum and synchronizability of the star-ring network under different coupling cases in Sections 3 and 4. Based on the theoretical analysis, simulations are conducted to investigate the effect of different topological parameters on synchronizability.In this section, synchronization control experiments on multilayer star-ring networks will be performed to verify the above analytical results.For simplicity, consider a multilayer network with M layers, each containing N nodes, which is described by the relevant descriptions for Eq (8) are similar to Eq (1), u K i (t) is an adaptive controller particularly.Assume the Hindmarsh-Rose [46,47] system is selected as the dynamic system described by:     where, x, y, and z are defined as the membrane potential of the neuron, the recovery variable, and the slow-varying adaptive current of the neuron in the system, respectively.(10) then synchronization errors of the networks are defined as: Subsequently, the error systems can be derived as follows: If limt→∞ ||ei K (t)|| = 0 is satisfied, it can be inferred that the multilayer complex networks and the target system have achieved synchronization.Designing an adaptive controller as follows: e H e (13) where h K i >0 and d K i (t) is the adaptive feedback control gains.Theorem 5.1.Based on Hypothesises 2.1 and 2.2, the error system is asymptotically synchronization under the controller, if there exists d * > 0 satisfying where, , Consider the following Lyapunov-krasovskii [48] candidate as:   e e (15) one gains Where, β = ((e (1) (t)) T , (e (2) (t)) T , …, (e (M) (t)) T ) T .The error systems ( 16) are shown to achieve asymptotic synchronization by Theorem 5.1.Therefore, the proof is now completed.Remark 2. In addition to the theoretical derivation of synchronizability in [49], the synchronization control of the constructed multilayer star-ring network model is added.Precise synchronization control can be achieved by changing the size of the topological parameters correlated to synchronizability.
In accordance with the master stability function theory, if the inner coupling continuous matrix Γ is an identity matrix, the synchronous region is unbounded.Therefore, the synchronizability of the network is determined by the λ2 of its Laplacian matrix [50].We select a divergent directionally coupled star-ring network for simulation, and N is even.From Table 1, λ2 = Md.The following topological parameters of the multilayer star-ring network are selected for numerical simulation: N = 50, d0 = 1, a = 1, a0 = 1, d = 1 (when considering the variation of M), and M = 3 (when considering the variation of d). Figure 15 shows the synchronization errors ||e|| versus time.Figure 15 reveals that the synchronization speed increases with M, indicating greater synchronizability at higher M. From Figure 15(b), the speed of synchronization exhibits a gradual increase with d.Therefore, the synchronizability becomes stronger as d increases.Figure 16 shows the evolution of the error in each dimension of the node over time for different numbers of layers M. As the number of layers M increases, it can clearly be seen that synchronizability improves.Figure 17 shows the variation of the error of each dimension node over time with d, and the results show that a larger d leads to faster synchronization.
Conclusive evidence is provided by the synchronization control experiments in this section, confirming the validity of the synchronizability analysis results.As the network structure becomes larger, leading to an increase in synchronizability.Additionally, the strength of interlayer coupling also enhances synchronizability.Consequently, it is possible to effectively implement control by determining the key topology parameters that have a correlation with synchronizability.

Conclusions and discussion
The eigenvalue spectrum of the supra-Laplacian matrix of multilayer star-ring networks is rigorously derived in this paper, considering divergent and aggregated coupling within the subnetwork.Compared with Sections 3 and 4, when the synchronized regions are bounded and unbounded, the variable coupled strengths σa and δa are closely connected to the synchronizability of the aggregated coupling but have no correlation with the divergent coupling.Synchronizability is firmly related to the topological parameters, including the number of nodes N, intralayer coupling strengths a and a0, interlayer coupling strengths d and d0, and the number of layers M in the divergent coupling model.In the case of an unbounded synchronous region, when N is odd, it is primarily influenced by the topological parameter a0.A bounded synchronous region depends on several factors, including N, a, a0, d0, and M. In addition, when N is even, M and d are relevant to the synchronizability with the unbounded synchronous region.When N is even, it is determined by N, a, d0, d, and M with the bounded synchronous region.Particularly, the synchronizability is significantly related to the intralayer variable coupling strengths σa and δa with aggregated coupling.When σ changes with different functions, the tendency toward synchronizability is found.Finally, a numerical simulation of a star-ring network composed of the Hindmarsh-Rose systems is given to verify the validity of the synchronizability analysis.
Based on the research on multilayer star-ring networks in this paper, coupling mode, coupling strength and topological parameters play an important role in the analysis of synchronizability.Through the analysis of synchronizability and the synchronization process, it has some theoretical reference value for understanding and exploring practical problems in engineering.In particular, adjusting traffic network signals, optimizing control in the Internet of Things and improving

Figure 1 .
Figure 1.Schematic diagram of the multilayer star-ring network with divergent coupling.

Figure 2 .Figure 3 .
Figure 2. λ2 and R of the divergent coupled star-ring network change with N.

4 )
Influence of the intralayer coupling strength a0 Take a = 1, a0 = 1, M = 15 (N = odd), M = 15, 16, 17 (N = even), d0 = 0.1, N (odd) = 399, and N(even) = 400) to explore the synchronizability with d.From Figure 5(a), λ2 is unchanged with d when N is odd, which means that synchronizability is unchanged.Additionally, λ2 increases with d by selecting three group data when N is even, and the synchronizability becomes stronger.From Figure 5(b), R is unchanged with d when N is odd, which means that synchronizability is unchanged.And R decreases with d by selecting three group data when N is even, the synchronizability becomes stronger.

Figure 4 .
Figure 4. λ2 and R of the divergent coupled star-ring network change with a0.

Figure 5 .
Figure 5. λ2 and R of the divergent coupled star-ring network change with d.
Figure 6(a) clearly indicates that λ2 is not affected by the change of d0, that is, the synchronizability is stable.As is shown in Figure 6(b), R increases with d0.It is revealed that the greater the coupling strength d0, the less favorable the synchronizability.6) Influence of the number of layers M Take a = 1, a0 = 1, d = 1 (N = odd), d = 1, 1.3, 1.6 (N = even), N (odd) = 399, and N (even) = 400 to analyze the synchronizability with M.

Figure 7 (
a) shows that λ2 is unchanged with M when N is odd, which means the synchronizability is unchanged.Additionally, λ2 increases with the M by selecting three group data, when N is even, the synchronizability becomes stronger with the M.
Figure 7(b)   shows that R increases slowly with M when N is odd, and then the synchronizability decreases.And R decreases with M when N is even; the synchronizability strengthens with M.

Figure 6 .
Figure 6.λ2 and R of the divergent coupled star-ring network change with d0.

Figure 7 .
Figure 7. λ2 and R of the divergent coupled star-ring network change with M.

Figure 8 .
Figure 8. Schematic diagram of the multilayer star-ring network with aggregated coupling.
Figure 13(a),(b) indicates the λ2 and R for varying coupling strength a.We could find λ2 and R vary periodically.

Figure 9 .
Figure 9. λ2 and R when the coupling strength varies as the proportional function.

Figure 10 .
Figure 10.λ2 and R when the coupling strength varies as inverse proportional function.

Figure 11 .
Figure 11.λ2 and R when the coupling strength varies as exponential function.

Figure 12 .
Figure 12. λ2 and R when the coupling strength varies as logarithmic function.

Figure 13 .
Figure 13.λ2 and R when the coupling strength varies as trigonometric function.

Figure 14 .
Figure 14.Hindmarsh-Rose system and time evolution of x, y, z.
The target system s(t) of the network is defined as follows:

Figure 15 .
Figure 15.The synchronization errors ||e|| of the multilayer divergent coupled star-ring networks evolve with time.
0 is the intralayer coupling strength, and H(xj K ):  N →  N is the inner coupling continuous function of the intralayer.p > 0 is the interlayer coupling strength, and Γ(xi Z ):  N →  N is the inner continuous coupling function of the interlayer.W K = (wij K )∈ NN represents the coupling weight configuration matrix of the Kth layer,