Impulse-Induced Optimum Signal Amplification in Scale-Free Networks

Impulse-Induced Optimum Signal Amplification in Scale-Free Networks

Pedro J. Martínez and Ricardo Chacón Departamento de Física Aplicada, E.I.N.A., Universidad de Zaragoza, E-50018 Zaragoza, Spain, and Instituto de Ciencia de Materiales de Aragón, CSIC-Universidad de Zaragoza, E-50009 Zaragoza, Spain Departamento de Física Aplicada, E.I.I., Universidad de Extremadura, Apartado Postal 382, E-06006 Badajoz, Spain, and Instituto de Computación Científica Avanzada, Universidad de Extremadura, E-06006 Badajoz, Spain
August 20, 2019

Optimizing information transmission across a network is an essential task for controlling and manipulating generic information-processing systems. Here, we show how topological amplification effects in scale-free networks of signaling devices are optimally enhanced when the impulse transmitted by periodic external signals (time integral over two consecutive zeros) is maximum. This is demonstrated theoretically by means of a star-like network of overdamped bistable systems subjected to generic zero-mean periodic signals, and confirmed numerically by simulations of scale-free networks of such systems. Our results show that the enhancer effect of increasing values of the signal’s impulse is due to a correlative increase of the energy transmitted by the periodic signals, while it is found to be resonant-like with respect to the topology-induced amplification mechanism.

89.75.Hc, 05.45.Xt, 05.60.-k, 89.75.Fb

Introduction.Today, there is an emerging ”network perspective” approach to studying complex systems reflecting the ubiquitous presence of networks in nature and in human societies. In particular, there has been considerable interest in a class of networks known as scale-free networks due to their lack of a characteristic size [1-4]. They have the property that the degrees, , of the node follow a scale-free power-law distribution . Examples are diverse metabolic and cellular networks, computer networks such as the World Wide Web, and some social examples such as collaboration networks. Besides topological investigations [5,6], current interest in these (and other) networks has extended to their controllability [7,8], i.e., to the characterization and control of the dynamical properties of processes occurring in them, such as transport [9], synchronization of individual dynamical behaviour occurring at a network’s vertices [10,11], role of quenched spatial disorder in the optimal path problem in weighted networks [12], and dynamic pattern evolution [13]. One particular issue that has attracted much interest because of its importance in both biological and man-made information-processing systems is the propagation and enhancement of resonant collective behaviour across a network due to the application of weak external signals. In this regard, there has been recently studied the amplification of the response to weak harmonic signals in networks of bistable signaling devices [14-17]. In these works, however, the robustness of the signal amplification against diversity in the uniform distributions of periodic external signals was not studied. Clearly, the assumption of harmonic external signals means that all driving systemswhatever they might beare effectively taken as linear. This mathematically convenient choice is untenable for most natural and artificial information-processing systems due to their irreducible nonlinear nature. Thus, to approach signal amplification phenomena in real-world networks, it seems appropriate to consider distributions of periodic external signals which are the output of nonlinear systems, therefore being appropriately represented by generic Fourier series.

In this Letter, we study the interplay between heterogeneous connectivity, quenched spatial disorder, and generic zero-mean periodic signals in random scale-free networks of signaling devices through the instance of a simple deterministic overdamped bistable system. The system is given by


where is a unit-amplitude, zero-mean, -periodic signal, is the signal amplitude, is the coupling, is the Laplacian matrix of the network, is the degree of node , and is the adjacency matrix with entries 1 if is connected to , and 0 otherwise. Since there are infinitely many a priori independent parameters, , , and hence infinitely many different waveforms of , the relevant problem is how to characterize quantitatively the effect of the signal’s waveform on the topology-induced amplification scenario [14]. Here, we shall show that a relevant quantity properly characterizing the effectiveness of generic periodic signals (1) in the amplification scenario is the impulse transmitted by the signal over a half-period (hereafter referred to simply as the impulse, ) a quantity integrating the conjoint effects of the signal’s amplitude, period, and waveform (see Fig. 1 for an example). Remarkably, we found that the enhancer effect of increasing values of the impulse is due to a correlative increase of the energy transmitted by the periodic signals. Extensive numerical simulations of the system (1) were conducted for different network topologies to characterize the effect of the impulse on the amplification-synchronization scenario as the coupling strength is increased.

Figure 1: (Color online) Top: Signal [Eq. (12)] vs where is the period and with , , , for four values of the shape parameter: (sinusoidal pulse), (nearly square-wave pulse), (double-humped pulse), and (sharp double-humped pulse). Bottom: The normalized impulse [Eq. (13)] vs .

To quantitatively describe this scenario, we used the average amplification over distinct initial conditions on the one hand, and the synchronization coefficient [18]


on the other, where the overlines indicate an average over nodes, while the angle brackets indicate a temporal average over a period .

Energy variation versus impulse.We start with a general argument showing the relationship between energy increases and impulse increases in isolated overdamped systems. Let us consider the family of dissipative nonlinear oscillators , with associated energy equation , where is a unit-amplitude, zero-mean, -periodic signal, is the energy function with being a generic potential. The energy equation can be recast into the form


Integration of Eq. (3) over any interval , , yields


Now, after applying the first mean value theorem for integration [19] to the last integral on the right-hand side of Eq. (4), one straightforwardly obtains


where and is the impulse. Let us consider an initial steady ( large) situation fixing the parameters and choosing a waveform such that the impulse is relatively small. Next, we only change the waveform. For sufficiently small values of and , one expects that both the dissipation work (integral in Eq. (5)) and will approximately maintain their initial values, while may increase from its initial value with the choice of a more convenient waveform, so that, in some cases depending upon the remaining parameters, the energy difference will increase with respect to the initial situation. Thus, the maximum probability of a maximal increase of the energy difference occurs when is also maximum, which establishes a clear correlation between impulse and energy transmitted over wide regions in parameter space. Specifically, one expects this impulse principle to be accurate when the signal period is the shortest significant timescale. Numerical simulations confirmed the accuracy and scope of this prediction. Figure 2 shows an illustrative example for an isolated overdamped bistable system (cf. Eq. (1)).

Figure 2: (Color online) Normalized average amplification versus shape parameter for an isolated overdamped bistable system [Eq. (1) with ] subjected to a signal given by Eq. (12) and seven values of the signal period: (circles), , (triangles), , (diamonds), (squares), and . The solid line shows the normalized impulse [Eq. (13)] versus shape parameter while the dashed vertical line indicates the value for which presents its single maximum as a function of . Notice that the period of the linearized motion around any of the potential minima, , corresponds to the shortest significant timescale.

Star-like network.Next, we consider a star-like network of overdamped bistable systems:


which describes the dynamics of a highly connected node (or hub), , and linked systems (or leaves), . We study the case of sufficiently small coupling, , and external signal amplitude, , such that the dynamics of the leaves may both be decoupled from that of the hub and be suitably described by linearizing their equations around one of the potential minima. Thus, one straightforwardly obtains


where while depending on the initial conditions. Since the initial conditions are randomly chosen, this means that the quantities behave as discrete random variables governed by Rademacher distributions. After inserting Eq. (7) into Eq. (6) and solving the resulting equation for the hub,



one straightforwardly obtains


where with being the equilibrium in the absence of any external signal, while is the hub potential with being the height of the potential barrier. For finite , the quantity behaves as a discrete random variable governed by a binomial distribution with zero mean and variance . One sees that the hub’s dynamics is affected by spatial quenched disorder through the term , while an estimate of its amplification [20] is given by


For sufficiently large , we may assume that the quantity behaves as a continuous random variable governed by a standard normal distribution, and hence


provides the final average amplification. Next, we demonstrate that the impulse is the relevant quantity controlling the effect of the external signal on the amplification by considering the illustrative example


in which and are Jacobian elliptic functions of parameter [ is the complete elliptic integral of the first kind] [21] and is a normalization function (see Fig. 1, top) which is introduced for the elliptic signal to have the same amplitude 1 and period for any wave form (i.e., ). When , then , i.e., one recovers the previously studied case of a harmonic signal [14], whereas, for the limiting value , the signal vanishes. Note that, as a function of , the impulse per unit of amplitude,


presents a single maximum at (see Fig. 1, bottom). In this case, Eq. (11) predicts that and that the signal amplification increases on average as the impulse is increased, i.e., as the shape parameter (see Fig. 3, left), which is accurately confirmed by numerical simulations (see Fig. 3, right). One also has from Eqs. (11) and (12) that , as a function of only , presents a sharp single maximum at for all , which indicates that the topology-induced amplification mechanism is robust against diversity in the uniform distributions of periodic signals in star-like networks.

Figure 3: (Color online) Theoretical average amplification in the parameter plane (left panel, Eqs. (11) and (12)) and corresponding numerical results (right panel) for a starlike network [Eq. (6)] with and . The dashed lines indicate the values , for which presents its single maximum in the parameter plane.

Scale-free network.Next, we discuss the possibility of extending the results obtained for a star-like network to Barabási-Albert (BA) networks [2] of the same overdamped bistable systems. Indeed, a highly connected node in the BA network can be thought of as a hub of a local star-like network with a certain degree picked up from the degree distribution. Thus, one can expect that the enhancer effect of the impulse will act at any scale to yield a significant enhancement of the signal amplification over the whole scale-free network in the weak coupling regime. Figure 4 shows an illustrative example where the averaged amplification is plotted against the coupling (top panel) and the shape parameter (bottom panel). One sees that becomes ever larger as approaches over the complete range of values of , confirming the predictions of the above theoretical analysis. Remarkably, the dependence of on follows in detail the respective dependence of the impulse (see Figs. 1 bottom and 4 bottom). We found that this scenario remains the same in any random realization of the network connectivity and for any value of (see Fig. 5, top panel). The synchronization decreases monotonically as is increased for any value of , with the decrease being ever faster as (Fig. 5, bottom panel). This can be understood as the result of two conjoint mechanisms: the impulse-induced enhancement of amplification and the amplification-induced lowering of synchronization in the weak coupling regime. Finally, we tested the robustness of the present results against possible finite-size effects: we found that the features of the impulse-induced amplification scenario remain the same for much larger networks [22].

Figure 4: (Color online) Top panel: Average amplification versus coupling for a BA scale-free network and four values of the shape parameter: , , (stars), and (squares). Note that the first relative maximum of occurs around for the four values of while the network has a maximal active hub having leaves. For this effective star-like network the theoretically predicted maximum occurs at . Bottom panel: Average amplification versus shape parameter for a BA scale-free network and three values of the coupling: , , and (squares). Averaged degree , , and the remaining fixed parameters are as in Fig. 3.
Figure 5: (Color online) Average of the average amplification over random realizations of the network connectivity for two values of the shape parameter ((circles) and (squares)) (top panel) and corresponding synchronization coefficient [Eq. (2)] for and (bottom panel) versus coupling for a BA scale-free network with . Other fixed parameters are as in Fig. 3.

Conclusion.We have shown through the example of a network of overdamped bistable systems that maximizing the impulse transmitted by the periodic external signals strongly enhances topology-induced signal amplification in scale-free networks. We have analytically demonstrated that this resonant-like effect of the impulse is due to a correlative increase of the energy transmitted by the periodic signals as the impulse is increased, while it may be completely characterized in the simple model of a starlike network. Remarkably, our results indicate that varying the impulse does not significantly change the values of the coupling strength for which amplification is maximum, which means that the topology-induced amplification mechanism is robust against diversity in the uniform distributions of periodic external signals. Thus, the present findings provide a reliable criterion for the optimization of topology-induced amplification processes in scale-free networks by controlling the periodic external signals according to the impulse principle. This novel principle opens up new avenues for studying external-signal-induced amplification processes in complex networks, including, for example, control of chaos by weak periodic excitations [23], or optimization of signal transmission in neuronal networks [24,25] in which a cooperative effect with the underlying noise of the neural medium is expected under certain conditions.

P.J.M. and R.C. acknowledge financial support from the Ministerio de Economía y Competitividad (MINECO, Spain) through FIS2011-25167 and FIS2012-34902 projects, respectively. P.J.M. acknowledges financial support from the Comunidad de Aragón (DGA, Spain. Grupo FENOL) and European Social Funds. R.C. acknowledges financial support from the Junta de Extremadura (JEx, Spain) through project GR15146.


  • (1) S. N. Dorogovtsev and J. F. F. Mendes, Evolution of networks: From biological nets to the Internet and WWW (Oxford Univ. Press, Oxford, 2003).
  • (2) A.-L. Barabási and R. Albert, Science 286, 509 (1999).
  • (3) R. Albert and A.-L. Barabási, Rev. Mod. Phys. 74, 47 (2002).
  • (4) G. Caldarelli, Scale-free Networks: Complex Webs in Nature and Technology (Oxford Univ. Press, Oxford, 2007).
  • (5) M. A. Serrano, D. Krioukov, and M. Boguñá, Phys. Rev. Lett. 100, 078701 (2008); Phys. Rev. Lett. 100, 199002(E) (2008).
  • (6) V. Colizza, A. Flammini, M. A. Serrano, and A. Vespignani, Nat. Phys. 2, 110 (2006).
  • (7) Y.-Y. Liu, J.-J. Slotine, and A.-L. Barabási, Nature 473, 167 (2011).
  • (8) T. Nepusz and T. Vicsek, Nat. Phys. 8, 568 (2012).
  • (9) E. López, S. V. Buldyrev, S. Havlin, and H. E. Stanley, Phys. Rev. Lett. 94, 248701 (2005).
  • (10) L. Donetti, P. I. Hurtado, and M. A. Muñoz, Phys. Rev. Lett. 95, 188701 (2005).
  • (11) C. Zhou and J. Kurths, Phys. Rev. Lett. 96, 164102 (2006).
  • (12) Y. Chen, E. López, S. Havlin, and H. E. Stanley, Phys. Rev. Lett. 96, 068702 (2006).
  • (13) H. Zhou and R. Lipowsky, Proc. Natl. Acad. Sci. U.S.A 102, 10052 (2005).
  • (14) J. A. Acebrón, S. Lozano, and A. Arenas, Phys. Rev. Lett. 99, 128701 (2007); Phys. Rev. Lett. 99, 229902(E) (2007).
  • (15) X. Liang, Z. Liu, and B. Li, Phys. Rev. E 80, 046102 (2009).
  • (16) T. Kondo, Z. Liu, and T. Munakata, Phys. Rev. E 81, 041115 (2010).
  • (17) J. Zhou, Y. Zhou, and Z. Liu, Phys. Rev. E 83, 046107 (2011).
  • (18) P. Balenzuela, P. Rué, S. Boccaletti, and J. García-Ojalvo, New J. Phys. 16, 013036 (2014).
  • (19) I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series, and Products (Academic Press, San Diego, 1980).
  • (20) This measure is based on the root-mean-square of a -periodic function.
  • (21) P. F. Byrd and M. D. Friedman, Handbook of Elliptic Integrals for Engineers and Scientists (Springer-Verlag, Berlin, 1971).
  • (22) See Supplemental Material at xxx for an additional figure.
  • (23) R. Chacón, Control of Homoclinic Chaos by Weak Periodic Perturbations (World Scientific, Singapore, 2005).
  • (24) J. J. Torres, I. Elices, and J. Marro, PLoS ONE 10(3), e0121156 (2015).
  • (25) S. J. Schiff, Neural Control Engineering (MIT Press, Cambridge, 2012).
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description