Enhancement of Network Synchronizability via Two Oscillatory System

The loss of synchronizability at large coupling strength is of major concern especially in the fields of secure communication and complex systems. Because theoretically, the coupling mode that can surely stabilize the chaotic/hyperchaotic synchronized state is vector coupling (using all the coordinates) which is in contrast to the practical demand of information exchange using lesser number of coordinates (commonly via a single coordinate). In the present work, we propose that if the node dynamics are given by a pair of oscillators (say, {\it two oscillatory system} TOS) rather than by a conventional way of single oscillator (say, {\it single oscillatory system} SOS), then the information exchange via a single coordinate could be sufficient to stabilize the chaotic/hyperchaotic synchronization manifold at large coupling strength. The frameworks of drive-response system and Master Stability Function (MSF) have been used to study the TOS effect by varying TOS parameters with and without feedback (feedback means quorum sensing conditions). The TOS effect has been found numerically both in the chaotic (R{\"o}ssler, Chua and Lorenz) and hyperchaotic (electrical circuit) systems. However, since threshold also increases as a side effect of TOS, the extent of $\beta$ enhancement depends on the choice of oscillator model like larger for R{\"o}ssler, intermediate for Chua and smaller for Lorenz.


I. INTRODUCTION
Synchronization is not always an obvious emergent behavior of the interacting dynamical systems even in the case when they are identical limit cycle oscillators. For an example, N identical Rössler oscillators (periodic/chaotic) coupled via a single coordinate (scalar coupling) on a ring network in the nearest neighborhood configuration exhibit desynchronization at large coupling strength [1,2]. This desynchronization behavior (known as Short Wavelength Bifurcation SWB [1]) arises because the minimum coupling strength (threshold) requires for the synchronization increases with the increase in number of oscillators whereas the synchronization manifold in case of scalar coupling mode, losses its stability (riddle basin behavior [3,4]) as the coupling strength (γ) augments towards a critical value (overload-tolerance). Thus, the small values of overload-tolerance (say γ dsyn ) limit the size of a synchronizable network, e.g. the synchronized state of a ring network having x 1 -coupled chaotic Rössler becomes unstable with the increase in N from 18 to 19 at γ dsyn =1.5. Moreover, since the mechanism which could increment γ dsyn has not been developed so far, the network modification methods are generally used to tackle the issue of strong coupling such as (1) adding the additional edges between the nodes deterministically (Pristine World) or/and stochastically (Small World) [5], (2) modifying the ring network topology to a synchronizable topology like a unidirectional tree network or a star network (a hub of Scale Free net- * harpartap@mail.saitama-u.ac.jp work) [6][7][8]. Practically these modifications could be considered as the distribution of overload among the nodes and theoretically these alterations imply the minimization of eigen-ratio (R) of the Laplacian/coupling matrix (largest eigenvalue to the smallest nonzero eigenvalue). The concept of R minimization comes from the theory of Master Stability Function (MSF) [9] which says that a complex network of size N is synchronizable if R < β where β = γ dsyn /γ syn (γ syn means threshold) [5]. Therefore, to study the chaotic complete synchronization behavior in the complex systems, finding the different pathways that could reduce R remain a primary goal for the researchers [5][6][7][8][9][10][11][12][13][14][15][16][17][18][19]. However, modification (1) increases the coupling cost of synchronization whereas the modification (2) decreases the robustness by increasing the centralization in a network.
In summary, the primary question which has not been answered yet is "how to achieve the finitely large overload-tolerance in case of scalar coupling as similar to the vector coupling scenario (wherein all coordinates are used), i.e. maximization of γ dsyn via scalar coupling"? The importance of this question lies on the fact that getting stability by using all the coordinates is neither useful (as in case of secure communication) nor realistic (as in case of complex system). Moreover, this theoretical problem also appears in many different forms such as: is it possible to stabilize the chaotic/hyperchaotic synchronization manifold? or is it possible to make a large non-centralized chaotic/hyperchaotic synchronizable network? or is it possible to surely stabilize the chaotic/hyperchaotic sub system under local/global parameter fluctuations?, etc. The answer to this question could be considered as a possible solution to the real problems such as overload failure in the real networks like arXiv:1706.02426v1 [nlin.CD] 8 Jun 2017 Internet system and power grid system [20,21], the stability issue of chaotic/hyperchaotic transmitter-receiver system in the field of secure communication [4,22,23], etc. It should be noted that since large coupling strength means infinite coupling (γ = ∞) in drive-response system [4,22,23], MSF complements the drive-response formulation by explicitly incorporating γ dependence which results into an elegant relation between the node property (β) of a just two node system (coupled bidirectionally) with the structural property (R) of an arbitrary network having N nodes (provided coupling matrix has zero row sum). In other words, if a system is stable in the drive-response framework then at the large coupling strength MSF shows negative values and vice-versa (discussed later).
In the present work, an attempt has been made to answer the primary question of maximization of γ dsyn via scalar coupling in case of linear interactions under small perturbations and global parameter fluctuations (each node experiences same fluctuations, i.e. identical node scenario). We argue that γ dsyn could be maximized, if the node dynamics are given by a pair of oscillators (say, two oscillatory system TOS) rather than by a conventional way of single oscillator (say, single oscillatory system SOS) as schematic shown in Fig. 1 on a ring network (i.e. a basic non-centralized network). In addition, since the stabilization of synchronization manifold happens only due to the emergence of dissipative factors via TOS (explained later by drive-response framework), to maximize γ dsyn the following two conditions should be met: 1. TOS should be in mTOS configuration as dTOS configuration behaves same as SOS (discussed later), 2. only that coordinate can be employed whose self dissipation could stabilize the unstable fixed point by adding its linear dissipative term in the autonomous system (discussed later). According to condition 2, all the three coordinates of Chua [24] and Lorenz [25] oscillators fulfill this criteria whereas only two out of the three coordinates of Rössler [26] (i.e. x 1 and x 2 ) and also only two out of the four coordinates of piecewise electronic hyperchaotic system [27] (i.e. x 1 and x 2 ), satisfy this condition. However, there also exist a possibility when none of the coordinates of an oscillator model like hyperchaotic Rössler [28], could meet the condition and hence can not be employed as the oscillatory dynamics for TOS. Furthermore, since threshold increases with the increase in number of oscillators, γ syn also increases as the side effect of mTOS along with the incrementation of γ dsyn . This makes β enhancement dependent on the choice of oscillator model like larger for Rössler, intermediate for Chua and smaller for Lorenz.
The paper is organized as follows. In Sec. II, the stability issue of scalar coupling and the theory of TOS are given in terms of MSF and drive-response frameworks along with the equations of the employed oscillator models. The stabilization via TOS for the scenarios of chaotic (Rössler, Chua and Lorenz) as well as hyperchaotic (electrical circuit) oscillators are presented in Sec. III. Finally, FIG. 1. Schematics of a ring network (nearest neighborhood configuration bidirectionally coupled) for two scenarios: each node behaves as (a) conventional single oscillatory system (SOS), (b) proposed two oscillatory system TOS. The subplots (b1) and (b2) depict the implementation of TOS, i.e. by using both the oscillators (dTOS) and one of the oscillator (mTOS), respectively, where 'd' means 'di' and 'm' means 'mono'. Oscillators x and y are coupled via all the coordinates (vector coupling) whereas oscillators x and x (y and y) are coupled via a single coordinate (scalar coupling).

FIG. 2.
Schematics of a two node network (unidirectional/bidirectional) wherein oscillators x and y are coupled via all the coordinates (vector coupling) whereas oscillators x and x (y and y) are coupled via a single coordinate (scalar coupling). the paper is concluded in Sec. IV.

II. THEORY
Under small perturbations (linear analysis), Master Stability Function (MSF) serves as the necessary and sufficient condition to judge the basin stability at large coupling strength. This is analogous to the driveresponse system wherein the conditional Lyapunov exponents surely determine the stability of a sub-system (response) by running the system from the different initial conditions. Therefore, the sub-Jacobian (Jacobian of response) method has been employed to investigate the reason behind the emergence of TOS effect.

A. MSF
To understand the formulism of MSF [9], consider the equations of motion for a complex network having N identical nodes for the conventional single oscillatory system (SOS): Where x i is a m-dimensional vector of i th node whose autonomous behavior is described by Here, H(x) (R m → R m ) represents m × m matrix (linear coupling function) which gives the information about the coordinates of x involved in the coupling, i.e. its all other elements are zero except a diagonal element corresponding to the employed coordinate (for more details see ref. [9]). In Eq. 1, G is a N × N coupling matrix (captures network's architect) which could be symmetric or not but it must be real with zero row sum [8] so that synchronous state (x i = s, i=1,..,N; s = F (s)) could become a solution of Eq. 1. For a symmetric case (Laplacian), G ii = −g i (g i is the degree/connections of i th node) and G ij = 1 if i th node is connected to j th node otherwise G ij = 0. The parameter γ(> 0) represents the coupling strength.
Corresponding to the N eigenvalues, there are N maximum Lyapunov exponents, λ(γδ k ), wherein λ(γδ 1 ) (λ(0) or λ 0 ) describes the state of the synchronized regime as λ 0 = 0 implies periodic and λ 0 > 0 means chaotic. The linear stability of this synchronized state is decided by the remaining N − 1 maximum Lyapunov exponents (transverse), i.e. the given synchronous state is stable if λ(γδ k ) < 0 (k > 1), and these N − 1 exponents can be found from a single variational equation by using scaling relation (γδ k ) [2]. This means that at fixed γ, the variational equation in the presence of δ max may experience stronger coupling and hence its stabilization or destabilization may occur earlier than the other modes. In other words, λ(γδ max ) plays a vital role in the stability condition wherein all the modes should be simultaneously stabilized or one could say that λ(γδ max ) can solely provides the stability of a network which is the case and hence it is called MSF for given DF (s) and DH(s). This is because, as long as G is diagonalizable its eigenvalues corresponding to different network topologies as well as sizes can be found from any topology using scaling relation (for more details see ref. [4,9]). Therefore, MSF of two bidirectionally coupled nodes, i.e. λ(2γ), yields same β (=γ dsyn /γ syn ) as an arbitrary network having N nodes does,i.e. β (=σ dsyn /σ syn ) where σ = γδ max . Now we can define γ syn and γ dsyn (similarly σ syn , σ dsyn ), i.e. 2γ syn is the minimum coupling strength at which λ(2γ)=0 or λ(2γ) → 0 + (emergence of synchrony) and 2γ dsyn is the maximum coupling strength at which λ(2γ)=0 or λ(2γ) → 0 − (emergence of desynchrony).

Problem of scalar coupling
In contrast to the vector coupling scenario (λ(σ → ∞) < 0), in case of scalar coupling wherein λ(σ) > 0 as σ > σ dsyn (σ dsyn is small), it becomes a challenge to synchronize a ring network topology ( Fig. 1(a)) because its R (i.e. 1/sin 2 (π/N )) grows faster with the increase in N than any other topology, e.g. for a star network R = N . Therefore, the previous works on the enhancement of network's synchronizability using network modification methods could be considered as the work on the minimization of R = 1/sin 2 (π/N ).

B. Drive-Response
To understand the framework of drive-response [22,23], consider the equations of motion for two unidirectionally coupled identical nodes for the conventional single oscillatory system SOS ( Fig. 2(a)): or diag(0,0,1) 3×3 (scalar coupling). Now we choose x 1 coordinate for the interactions and assume that the coupling term vanishes, i.e. x 2 1 → x 1 1 , as γ → ∞ which is only possible ifẋ 2 1 =ẋ 1 1 . This means that x 2 node instead of generating its own x 2 1 signal must use x 1 1 signal from x 1 node in order to ensure x 2 1 = x 1 1 (γ = ∞) for all the time, i.e. even when x 2 j = x 1 j (j = 2, .., m). Hence, Eq. 3 becomesẋ 2 = F (x 2 , x 1 1 ) andẋ 2 =(ẋ 2 2 ,..,ẋ 2 m ). In this scenario, Eq. 3 depicts x 1 -driving system wherein x 1 is drive and x 2 is a response sub system. The maximum conditional Lyapunov exponent (λ c ) of sub system (sub Jacobian) decides the stability of synchronization manifold, i.e. λ c < 0 implies stabilization (for more details see ref. [4]). Thus, this becomes analogous to the MSF behavior of two bidirectionally x 1 -coupled nodes wherein the condition of stabilized synchronization manifold is λ(2γ) < 0. Therefore, by using scaling relation one could also relate with the N nodes scenario, i.e. λ(γδ max ) < 0 (as discussed above).
(depending upon the choice of oscillator) which sometimes can not provide the stabilization against the parameter fluctuation as the scenario of x 2 -driving chaotic Rössler (shown in ref [23]). Secondly, to stabilize a hyperchaotic drive-response usually a scalar signal is generated by using BK method [32] which requires many parameter (2m). Hence, this necessitates some simpler mechanism which could work for both chaotic as well as hyperchaotic scenarios.

C. TOS
In contrast to SOS (ẋ = F (x)), the proposed two oscillatory system (TOS) is made up of two bidirectionally coupled oscillators (identical/nonidentical) by employing their all coordinates (vector coupling) so that the stable synchronization behavior (complete/generalized) could be ensured. Here, one should not confuse the used terminology of two oscillators with two nodes, i.e. in TOS scenario, the dynamics of a single node are generated by a pair of oscillators as: Similar to x, y is also a m-dimensional vector whose autonomous behavior is described by same function F (y) (R m → R m ) where µ and µ are the intrinsic parameters of the models. The parameters γ 1 and γ 2 represent the intra-coupling constants of TOS and Γ v shows the vector coupling scenario between x and y, i.e. Γ v =diag(1,..,1) m×m . For the scenario of γ 2 =θγ 1 (feedback), TOS could be considered as a two oscillatory representation of the quorum sensing network of x oscillators interacting with each other via a y oscillator (medium) [31]. The parameter θ represents the population density of x at each node of a network whereas y may have periodic/chaotic dynamics, i.e. other than the conventional steady state behavior used in quorum sensing [29][30][31].
Further, the nodes having TOS dynamics may interact with each other in two ways, i.e. (i) each node uses either its x or y (mTOS, also quorum sensing type) and (ii) each node uses its both x and y (dTOS), as shown in Fig. 1(b) and Fig. 2(b). It has been found that only former way (mTOS) could lead to TOS effect since the latter way (dTOS) behaves same as SOS (found analytically as well as numerically). Moreover, since threshold increases with the increase in number of oscillators, γ syn also increases as the side effect of mTOS along with the maximization of γ dsyn . This situation is more clear in terms of a two node network (Fig. 2) wherein SOS is the case of two oscillators ( Fig. 2(a)) whereas mTOS depicts the scenario of four oscillators (Fig. 2(b2)).
To demonstrate the TOS effect in terms of MSF (λ(σ → ∞) < 0) and drive-response system (λ c < 0) for scalar coupling, we use a ring network of size N (Fig. 1) and a two node network (Fig. 2), respectively. Furthermore, it should be noted that TOS could be ap-plied to any arbitrary network topology since TOS alters only the node dynamics not the network structure, e.g. a ring network of TOS has same R as for SOS, i.e. R = 1/sin 2 (π/N ) where δ ring min = 4sin 2 (π/N ) and δ ring max = 4 (N is even [2]).

MSF in case of TOS
The equations of motion for the nodes of a ring network having TOS dynamics ( Fig. 1(b)), are: Fig. 1(b1) and (b2), Eq. 4 and Eq. 5 respectively, represent the scenarios of N identical diffusively coupled dTOS and mTOS. The bold terms of Eq. 4-5 show the nearest neighborhood interactions To understand the different behavior of dTOS and mTOS, consider the block diagonalized linear variational equations of Eq. 4-5 by using technique of spatial Fourier modes (given in ref. [34]): Where η k Moreover, DF (s, µ) and DF (y s , µ ) are the Jacobian matrices evaluated at the synchronization manifold: The bold terms of Eq. 6-7 depict the form of γδ k DH(s) (Eq. 2) in case of ring network topology. Eq. 6 and 7 are corresponding to Eq. 4 (dTOS) and 5 (mTOS), respectively. By using η k = γ 2 η k x + γ 1 η k y , Eq. 6 and Eq. 7 become: To simplify, we assume DF (s, µ) ≈ DF (y s , µ ) = D since TOS due to vector coupling (Γ v ) and strong values of intra-coupling constants (γ 1 ,γ 2 ), could behave as a system whose divergence rates from the synchronization manifold (complete/generalized) for both the oscillators could be equal even when δµ = 0 (δµ = µ − µ). This assumption is valid for every oscillator model for a given range of δµ (similar type of assumption had also been used previously [35]). Therefore using this assumption, Eq. 8 (dTOS) becomes: η k = (D−4γsin 2 (πk/N)Γ)η k (10) and Eq. 9 (mTOS) becomes: The form of Eq. 10 evidently shows that dTOS behaves same as SOS (Eq. 2). On the contrary, Eq. 11 depicts that in case of mTOS, the perturbation η also depends upon its x-component (η x ) in addition to coupling strength (γ) and eigenmodes (k). Moreover, since the evolution of η x does not depend on γ and k (Eq. 11b), it could act as a stability factor. However, this analysis (without numerics) does not provide any hint that γ-independent Eq. 11b could lead to the maximization of γ dsyn . Thus, we need to further investigate by using drive-response framework.

Drive-Response in case of TOS
The equations of motion for two unidirectionally coupled identical nodes having TOS dynamics are: Similar to Eq. 3, the node x 1 − y 1 is drive whereas the node x 2 − y 2 is response. Eq. 12 and Eq. 13 depict the scenarios of dTOS and mTOS, respectively. Now we choose x 1 , y 1 coordinates for dTOS whereas y 1 coordinate in case of mTOS for the interactions and assume that the coupling term vanishes at γ → ∞. Hence, Eq. 12 and 13 become Eq. 14-15 and Eq. 16-17, respectively (same as discussed for Eq. 3): Eq. 14-15 depict that for dTOS, x 2 and y 2 depend only on their respective missing components as x 1 1 only drives x 2 and y 1 1 only drives y 2 which is similar to Eq. 3. On the other hand, Eq. 16 shows that in case of mTOS both x 2 and y 2 are driven by same y 1 1 . Now we need to analyze y 1 1 driving effect on x 2 (which is missing in dTOS). For this we decompose only Eq. 16 (because Eq. 17 is same as Eq. 15): Where i = 2, .., m. The relevance of y 1 1 -drive would clearly emerge when one takes partial derivative of Eq. 18 w.r.t. x 2 1 which results into the addition of −γ 1 in the first element of sub-Jacobian whose counter balance term, i.e. +γ 1 , is missing in the non diagonal element (i.e. 0) because Eq. 17 does not haveẏ 2 1 . Thus, this uncompensated −γ 1 becomes the source of TOS effect, i.e. λ c < 0 (λ(σ → ∞) < 0). This is in contrast to the sub-Jacobian for all other cases (Eq. 14-15, 17 and 19) wherein the negative term (−γ 1 /γ 2 ) of the diagonal elements gets balanced by the positive terms (γ 1 /γ 2 ) of non-diagonal elements. Therefore, we can say that Eq. 18 demonstrates the role of γ-independent Eq. 11b towards the maximization of γ dsyn . In addition, Eq. 18 also explains the need of condition 2 (Sec. I), i.e. x 1 (y 1 ) should stabilize the unstable fixed point of x (y) since the TOS effect emerges only due to the additional linear dissipative term (−γ 1 x 1 ). Moreover, this also shows the relevance of γ 1 , i.e. to induce the TOS effect γ 101 ≤ γ 1 ≤ γ 102 where γ 101 and γ 102 depend on the choice of oscillator model. In the present work, the lower bound γ 101 has been found by exploiting the quorum sensing conditions [33] (the upper bound γ 102 has been located just by scanning the intra-coupling parameter space).

Chaotic
Rössler oscillator [26]: where d and µ R are the bifurcation parameters.

III. RESULTS AND DISCUSSIONS
In the present work, we employ 5 different scenarios: x 1 and x 2 -coupled/driving chaotic Rössler, x 3 -coupled/driving chaotic Chua, x 3 -coupled/driving chaotic Lorenz and x 1 -coupled/driving hyperchaotic oscillator, to generalize the TOS effect, i.e. the stabilization of chaotic/hyperchaotic synchronization manifold (λ c < 0) or the maximization of γ dsyn (λ(4γ → ∞) < 0) in case of scalar coupling. To reiterate, since these different scenarios implies different DF and/or DH (Eq. 2), each scenario has its own MSF. In addition, 'coupled' term depicts explicit γ dependence, i.e. the stability is shown by MSF using λ(4γ) (ring network) whereas 'driving' term shows the drive-response framework wherein the stabilization is given by maximum conditional Lyapunov exponent λ c of sub-Jacobian. Furthermore, the robustness of TOS effect has also been shown by varying the TOS parameters (γ 1 , γ 2 ) and the model parameters of Rössler (µ R ), Chua (µ C ) and Lorenz (µ L ). It has been found that the intra-coupling feedback situation (discussed in Sec. IIC) of the TOS parameters provides more robustness but lesser β enhancement, in comparison to no feedback situation.

A. Rössler
The emergent effect of TOS for Rössler dynamics is given in Fig. 3-7 wherein Fig. 3 depicts the stabilization of x 2 -driving Rössler using λ c against the variation of its bifurcation parameter d for both SOS and TOS configurations. Fig. 3(a) shows the scenario (which had also been previously reported in figure 10 of ref [23]) wherein the subsystem (x 2 -driven Rössler) loses its stability close to the bifurcation points especially at the periodic windows in case of SOS. In contrast to Fig. 3(a), Fig. 3(b) evidently demonstrates that the mTOS configuration stabilizes the subsystem at all the domains  whereas the dTOS configuration exhibits qualitatively similar behavior as of SOS (consistent with the theory). Moreover, Fig. 4 shows the MSF behavior for x 2coupled chaotic Rössler (for SOS and TOS) corresponding to d = 6.8 (chaos) of Fig. 3 at which synchronization manifold is stable for mTOS (λ c < 0) and unstable for SOS (λ c > 0). Now it is evident from Fig. 4 that λ(4γ → ∞) < 0 only happens for mTOS configuration (where '∞' means finitely large value). This also proves that at the large coupling strength MSF exhibits same information as the corresponding drive-response system does, i.e. λ(4γ → ∞) ≈ λ c . Analogous to Fig. 4, MSF behavior shown in Fig. 5 for x 1 -coupled chaotic Rössler also illustrates that only mTOS stabilizes the chaotic synchronization manifold or maximizes γ dsyn . Moreover, the insets of Fig. 4 and 5 show the incrementation of γ syn as the side effect of mTOS configuration. Furthermore, the mTOS effect shown in Fig. 5 is studied under parameter fluctuations, i.e. by varying TOS's critical parameter γ 1 (discussed in Sec. IIC2) and the model parameter µ R (nonidentical TOS), for no-feedback (Fig. 6) and feedback (Fig. 7) scenarios. In contrast to Fig. 6(a) wherein weak-mTOS effect observed at both small and large values of γ 1 , Fig. 7(a) evidently demonstrates the robustness of mTOS effect (λ(4γ → ∞) < 0) over the broader range of parameter space in case of feedback scenario. Moreover, it also suggests that for the induction of mTOS effect at large γ 1 , γ 1 should be less than γ 2 (one can say that the upper bound of γ 1 , γ 102 , also depends on γ 2 along with the model parameters). However, the behavior of incrementation of γ syn increases with the increase in γ 2 which reduces β enhancement as depicted by Fig. 7(b)-(d) in comparison to Fig. 6(b)-(d) wherein γ 2 is fixed. Moreover, it should be noted that the extra strong-enhancement observed only in no-feedback scenario ( Fig. 6(b)) is the consequence of decrement of γ syn due to the stable periodic synchronization manifold, λ 0 = 0 (Fig. 6(e)), which emerge at γ 1 → γ 2 for µ R < 0 (Rössler's steady state domain). But this behavior is missing in case of feedback scenario (Fig. 7(b)) because γ 2 increases with the increase in γ 1 which results into chaotic dynamics, λ 0 > 0 (Fig. 7(e)), even at µ R < 0.

B. Chua
Analogous to Fig. 5, Fig. 8 demonstrates the mTOS effect for x 3 -coupled chaotic Chua and the robustness of this effect is given in Fig. 9-10 (similar to Fig. 6-7). Here also, the feedback scenario (Fig. 10) provides more robustness but lesser β enhancement than the case of no-feedback (Fig. 9). Moreover, by comparing Fig. 9(a) with Fig. 6(a), one could realize the dependence of γ 102 on the oscillator model whereas Fig. 10(a) depicts that this dependence gets reduced in the presence of feedback (similar to Fig. 7(a)). But the major difference arises due to change in oscillator model from Rössler to Chua, is the extent of β enhancement as shown in Fig. 9(b)-10(b) which is lesser in contrast to Fig. 6(b)-7(b). This is because of more augmentation of γ syn in case of Chua ( Fig. 9(c)-(d) and 10(c)-(d)) than the case of Rössler (Fig. 6(c)-(d) and 7(c)-(d)). Instead of this discrepancy, Fig. 9(a)-10(a) in conjunction with Fig. 9(e)-10(e) clearly demonstrate the independence of γ dsyn maximization from the state of synchronization manifold which is similar to the case of Rössler.

C. Lorenz
Similar to Rössler and Chua, the case of x 3 -coupled chaotic Lorenz given in Fig. 11-13 also depict more robustness but lesser β enhancement in case of feedback (Fig. 13) than the case of no-feedback (Fig. 12). In addition, consistent with Rössler and Chua, Fig. 12(a)-13(a) in conjunction with Fig. 12(e)-13(e), evidently demonstrate the independence of γ dsyn maximization from the state of synchronization manifold. Again the difference emerge in the extent of β enhancement (Fig. 12(b)-13(b)) which is much smaller in contrast to Rössler because γ syn augments much more in case of Lorenz.

D. Hyperchaotic
Fig. 14 demonstrates the mTOS effect, λ(4γ → ∞) < 0, in case of x 1 -coupled hyperchaotic electronic system. Moreover, it is worth noticing that since λ(4γ → ∞) < 0 implies λ c < 0 (as shown for Rössler), we can say that a hyperchaotic drive-response system could also be stabilized without employing BK method [32]. Furthermore, similar to the above chaotic oscillators, in case of hyperchaotic oscillator feedback scenario yields more robustness but lesser β enhancement than the scenario of nofeedback and the state of synchronization manifold does not effect the maximization of γ dsyn (results not shown).

IV. CONCLUSIONS
In the present work, without altering the network topology and without employing multi parameter BK method (previous approaches), an effort has been made to stabilize the chaotic as well as hyperchaotic synchronization manifold at large coupling strength in case of scalar coupling, just by modifying the node dynamics from a single oscillator (SOS) to a pair of oscillators (TOS). The impact of modified node configuration for both chaotic and hyperchaotic dynamics, has been studied using Master Stability Function MSF (for a fundamental non-centralized ring network) along with the drive-response framework (for a two node network). The presented results evidently dictates the potency of TOS in mTOS configuration, i.e. maximization of γ dsyn or λ(4γ → ∞) < 0 or λ c < 0, under the broad range of parametric fluctuations. Furthermore, results also demonstrate the enhancement of β even when the threshold (γ syn ) augments to large values. Therefore, we believe that the proposed TOS would play a vital role both in the fields of complex system as well as secure communication wherein the network synchronizability at large coupling strength in case of scalar coupling, remains a primary concern especially for the hyperchaotic systems.

V. ACKNOWLEDGEMENT
Author would like to thank Vandana H. Singh for her valuable suggestions to improve the manuscript presentation. Financial support (501T015013, 502T016522 and TJC23) from Saitama University Japan, is acknowledged.