Abstract
A dynamical system is controllable if by imposing appropriate external signals on a subset of its nodes, it can be driven from any initial state to any desired state in finite time. Here we study the impact of various network characteristics on the minimal number of driver nodes required to control a network. We find that clustering and modularity have no discernible impact, but the symmetries of the underlying matching problem can produce linear, quadratic or no dependence on degree correlation coefficients, depending on the nature of the underlying correlations. The results are supported by numerical simulations and help narrow the observed gap between the predicted and the observed number of driver nodes in real networks.
Similar content being viewed by others
Introduction
While during the past decade significant efforts have been devoted to understanding the structure, evolution and dynamics of complex networks1,2,3,4,5,6, only recently has attention turned to an equally important problem: our ability to control them. Given the problem's importance, recent work has extended the concept of pinning control7,8,9 and structural controllability10,11,12,13 to complex networks. Here we focus on the latter approach. A networked system is considered controllable if by imposing appropriate external signals on a subset of its components, called driver nodes, the system can be driven from any initial state to any final state in finite time14,15,16,17. As the control of a system requires a quantitative description of the governing dynamical rules, progress in this area was limited to small engineered systems. Yet, recently Liu et al.10 showed that the identification of the minimal number of driver nodes required to control a network, ND, can be derived from the network topology by mapping controllability16 to the maximum matching in directed networks18. The mapping indicated that ND is mainly determined by the degree distribution P(kin, kout). We know, however, that a series of characteristics, from degree correlations19,20,21 to local clustering22 and communities23,24,25,26, cannot be accounted for by P(kin, kout) alone, prompting us to ask: which network characteristics affect the system's controllability?
The three most commonly studied deviations from the random network configuration are (i) clustering, manifested as a higher clustering coefficient C than expected based on the degree distribution27; (ii) community structure, representing the agglomeration of nodes into distinct communities, captured by the modularity parameter Q25; (iii) degree correlations28. First, we motivate our work by showing that network characteristics other than the degree distribution also affect network control. Next, we use numerical simulations to identify the network characteristics that affect controllability, finding that only degree correlations have a discernible effect. We then analytically derive nD = ND/N for random networks with a given degree distribution and correlation profile. More detailed calculations are provided in the Supplementary Information Sec. III. Finally, we test our predictions on real networks.
Results
Prediction based on the degree distribution
To motivate our study we compared the observed ND to the prediction based on the degree sequence for several real networks. For this we randomize each network preserving its degree sequence and we calculate , the number of driver nodes for the randomized network. Plotting ND versus
on log-log scale indicates that the degree sequence correctly predicts the order of magnitude of ND despite known correlations19,20 (Fig. 1a). However, by plotting nD = ND/N versus
we observe clear deviations from the degree based prediction (Fig. 1b). Our goal is to understand the origin of these deviations and the degree to which network correlations can explain the observed nD.
(a) We compare ND for real systems to , representing the number of driver nodes needed to control their randomized counterparts. Randomization eliminates all local and global correlations, only preserving the degree sequence of the original system. We find that the degree sequence predicts the order of magnitude of ND correctly, however, small deviations are hidden by the log scale, needed to show the whole span of ND seen in real systems. (b) These deviations are more obvious if we compare the density of driver nodes nD = ND/N and
in linear scale, finding that for some systems (e.g. regulatory and p2p Internet networks) the degree sequence serves as a good predictor of nD, while for other systems (e.g. metabolic networks and food webs) nD deviates from the prediction based solely on the degree sequence.
Numerical simulations
We start from a directed network with Poisson29,30 or scale-free degree distribution31,32. The scale-free network is generated by the static model described in the Methods section. We use simulated annealing to add various network characteristics by link rewiring, while leaving the in- and out-degrees unchanged, tuning each measure to a desired value, for details see the Methods section. We computed nD using the Hopcroft-Karp algorithm33.
Clustering
We use the global clustering coefficient27 defined for directed networks as
![](http://media.springernature.com/lw326/springer-static/image/art%3A10.1038%2Fsrep01067/MediaObjects/41598_2013_Article_BFsrep01067_Equ1_HTML.jpg)
The simulations indicate that changes in C only slightly alter nD and that the effect is not systematic (Fig. 2a). Hence we conclude that C plays a negligible role in determining nD.
Modularity
We quantify the community structure using25,26:
![](http://media.springernature.com/lw315/springer-static/image/art%3A10.1038%2Fsrep01067/MediaObjects/41598_2013_Article_BFsrep01067_Equ2_HTML.jpg)
where Avw is the adjacency matrix, cv and cw are the communities the v and w nodes belong to, respectively. Specifying Q still leaves a great amount of freedom in the number and size of the communities. We therefore choose to randomly divide the nodes into NC equally sized groups and increase the edge density within these groups, elevating Q to the desired value.
The simulations indicate that this community structure has no effect on nD (Fig. 2b). While adding communities to networks can be achieved in many different ways and the effect of modularity can be explored in more detail (e.g. hierarchical organization of communities23,34,35, overlapping community structure24,36, etc), we have failed to detect systematic, modularity induced changes in nD, prompting us to conclude that Q does not play a leading role in nD.
Degree correlations
In directed networks each node has an in-degree (ki) and an out-degree (ko), thus we can define four correlation coefficients: correlations between the source node's in- and out-degree and the target node's in- and out-degree (Figs. 3, 4)28. We use the Pearson coefficient to quantify each correlation with a single parameter:
![](http://media.springernature.com/lw333/springer-static/image/art%3A10.1038%2Fsrep01067/MediaObjects/41598_2013_Article_BFsrep01067_Equ3_HTML.jpg)
where · sums over all edges, α, β ∈ {in, out} is the degree type, k(α) is the degree of the source node, j(β) is the degree of the target node. And
is the average degree of the nodes at the beginning of each link,
is the variance;
and σ(β) are defined similarly.
The impact of degree-degree correlations on the density of driver nodes (nD) for the Erdős-Rényi model (N = 10, 000) for average degrees 〈k〉 = 1 (red), 〈k〉 = 3 (green), 〈k〉 = 5 (blue), 〈k〉 = 7 (black) and 〈k〉 = 9 (orange).
The results are similar for the scale-free model (see Fig. 4). Each data point is an average of 100 independent runs.
The impact of degree-degree correlations on the density of driver nodes (nD) for the scale-free model (N = 10, 000, γ = 2.5) for average degrees 〈k〉 = 1 (red), 〈k〉 = 3 (green), 〈k〉 = 5 (blue), 〈k〉 = 7 (black) and 〈k〉 = 9 (orange).
The results are similar for the Erdős-Rényi model (see Fig. 3). Each data point is an average of 100 independent runs.
Simulations shown in Figs. 3 and 4 indicate that degree correlations systematically affect nD. We observe three distinct types of behavior:
-
i
nD depends monotonically on r(out-in), so that low (negative) correlations increase nD and high (positive) correlations lower nD (Figs. 3c, 4c);
-
ii
Both r(in-in) and r(out-out) increase nD, independent of the sign of the correlations (Figs. 3a, 3d, 4a, 4d);
- iii
The behavior is qualitatively the same for Erdős-Rényi (Fig. 3) and scale-free (Fig. 4) networks.
The diversity of these numerical results require a deeper explanation. Therefore in the remaining of the paper we focus on understanding analytically the role of degree correlations, which, by systematically altering nD, affect the system's controllability.
Analytical framework
The task of identifying the driver nodes can be mapped to the problem of finding a maximum matching of the network10. A matching is a subset of links that do not share start or end points. We call a node matched if a link in the matching points at it and we gain full control over a network if we control the unmatched nodes. The cavity method has been successfully used to calculate the size of the maximum matching for undirected37 and directed10 network ensembles with given degree distribution. Here we study network ensembles with a given degree correlation profile.
We calculate nD analytically for a given P(kin, kout) and selected degree-degree correlation e(jin, jout; kin, kout), representing the probability of a directed link pointing from a node with degrees jin and jout to a node with degrees kin and kout. In the absence of degree correlations (neutral case)
![](http://media.springernature.com/lw367/springer-static/image/art%3A10.1038%2Fsrep01067/MediaObjects/41598_2013_Article_BFsrep01067_Equ4_HTML.jpg)
where ,
and 〈k〉 is the average degree. To ensure analytical tractability we chose21
![](http://media.springernature.com/lw411/springer-static/image/art%3A10.1038%2Fsrep01067/MediaObjects/41598_2013_Article_BFsrep01067_Equ5_HTML.jpg)
![](http://media.springernature.com/lw411/springer-static/image/art%3A10.1038%2Fsrep01067/MediaObjects/41598_2013_Article_BFsrep01067_Equ6_HTML.jpg)
![](http://media.springernature.com/lw411/springer-static/image/art%3A10.1038%2Fsrep01067/MediaObjects/41598_2013_Article_BFsrep01067_Equ7_HTML.jpg)
![](http://media.springernature.com/lw411/springer-static/image/art%3A10.1038%2Fsrep01067/MediaObjects/41598_2013_Article_BFsrep01067_Equ8_HTML.jpg)
By fixing m(α–β)(j, k) (α, β ∈ {in, out}) we obtain a one parameter network ensemble characterized by r(α–β), where m(α–β)(j, k) satisfies the constraints
![](http://media.springernature.com/lw327/springer-static/image/art%3A10.1038%2Fsrep01067/MediaObjects/41598_2013_Article_BFsrep01067_Equ9_HTML.jpg)
![](http://media.springernature.com/lw303/springer-static/image/art%3A10.1038%2Fsrep01067/MediaObjects/41598_2013_Article_BFsrep01067_Equ10_HTML.jpg)
and all elements of e(α–β)(j, k) are between 0 and 1.
Our goal is to understand the relation between nD and the degree correlation coefficient r(α–β). Assuming that r(α–β) is small we treat the correlations as perturbations to the neutral case, discussing the impact of the four r(α–β) correlations separately.
Out-in correlations
Using equation (5c) and keeping the first nonzero correction we obtain (Supplementary Information Sec. III.):
![](http://media.springernature.com/lw414/springer-static/image/art%3A10.1038%2Fsrep01067/MediaObjects/41598_2013_Article_BFsrep01067_Equ11_HTML.jpg)
where is the fraction of driver nodes of the uncorrelated network; wi and
only depend on P(kin, kout)10 and
![](http://media.springernature.com/lw327/springer-static/image/art%3A10.1038%2Fsrep01067/MediaObjects/41598_2013_Article_BFsrep01067_Equ12_HTML.jpg)
Equation (8) predicts that depends linearly on r(out-in), a prediction supported by simulations for small r(out-in) (Figs. 3c and 4c). This behavior is also revealed by the equivalent problem of finding the maximum matching of graphs10. For a node A with out-degree k0, by definition only one edge can be in the matching. If the remainder k0 – 1 edges point to nodes with degree 1 (disassortative case), A inhibits them from being matched, so we have to control each of them individually, increasing nD. If the remainder k0 – 1 edges point to hubs (assortative case), these hubs are likely to be matched through another incoming edge, decreasing nD.
Out-out correlations
The cavity method indicates that for out-out correlations the first nonzero correction is of order (r(out-out))2:
![](http://media.springernature.com/lw411/springer-static/image/art%3A10.1038%2Fsrep01067/MediaObjects/41598_2013_Article_BFsrep01067_Equ13_HTML.jpg)
where (α ∈ {in, out}) only depends on P(kin, kout) and
![](http://media.springernature.com/lw388/springer-static/image/art%3A10.1038%2Fsrep01067/MediaObjects/41598_2013_Article_BFsrep01067_Equ14_HTML.jpg)
Equation (10) predicts that does not depend on the out-out correlation of the directly connected nodes, but only on the correlation between the second neighbors, hence its dependence is quadratic in r(out-out), a prediction supported by numerical simulations (Figs. 3d and 4d). Indeed, positive (negative) r(out-out) correlation between the immediate neighbors means that if node A has high out-degree, then node B is expected to have high (low) out-degree and therefore C is likely to have high out-degree (Fig. 5). That is, both positive and negative one-step out-out correlations induce positive two-step correlations, accounting for the symmetry of the effect observed in simulations (Figs. 3d and 4d).
In-in correlations
Switching the direction of each link does not change the matching, but turns out-out correlations into in-in correlations. So can be obtained by exchanging P(in)(kin) and P(out)(kout) in equation (10), predicting again a quadratic dependence on r(in-in), supported by the numerical simulations (Figs. 3a and 4a).
In-out correlations
The equations for do not depend on the in-degree of the source and the out-degree of the target of a link, hence we predict that r(in-out) does not play a role in network controllability, a prediction supported by the simulations (see Figs. 3b and 4b).
Taken together, we predict that the functional dependence of on degree correlations defines three classes of behaviors, depending on the matching problem's underlying symmetries:
has no dependence on r(in-out), linear dependence on r(out-in) and quadratic dependence on r(in-in) and r(out-out). These predictions are fully supported by numerical simulations Figs. 3 and( 4): for small r we see no dependence on r(in-out), an asymmetric, monotonic dependence on r(out-in) and a symmetric on r(in-in) and r(out-out).
To directly compare the analytical predictions to simulations we need to know the complete e(ji, jo; ki, ko) distribution, which is not explicitly set in our simulations. So to test the results we use a rewiring method that sets the e(ji, jo; ki, ko) distribution, not only the r correlation coefficient21. This method is not as robust as our original algorithm and the range of accessible r values is more restricted. However, since our results are based on perturbation scheme we only expect them to be correct for small r values. Indeed, we find that the predictions quantitatively reproduce the numerical results in a fair interval of r(α–β) (Fig. 6).
The analytic formulas are tested with simulations on an (a) Erdős-Rényi model and on a (b) scale-free model.
We used the algorithm proposed in21 to set e(α−β)(ji, jo; ki, ko). For (a) network we choose N = 1, 000 and 〈k〉 = 3; for (b) N = 1, 000, γ = 2.5 and 〈k〉 = 4. Each data point is an average over 100 independent runs; the errors represent by the standard deviation of the measurements.
Real networks
We test the predictions provided by the developed analytical and numerical tools on a set of publicly available network datasets. When complex systems are mapped to networks, the links connecting the nodes represent interactions between them. In this context self-loops represent self-interactions, with a strong, well understood impact on controllability10,38. While in some systems self-loops are obviously present (e.g. neural networks), in others they are manifestly absent (e.g. electric circuits39). Our purpose here is to test the effect of correlations, hence we rely on datasets that capture the wiring diagram of various complex systems with different correlation properties. Therefore, even if in a few of these maps self-loops are missing, it is beyond the scope of this work to complete these networks. However, when studying controllability of a particular system, careful thought has to be put into whether self-loops are present or not. We present a systematic study on the effect of self-loops in the Supplementary Information Sec. II.B.
To test the impact of our predictions on real networks we calculate
![](http://media.springernature.com/lw258/springer-static/image/art%3A10.1038%2Fsrep01067/MediaObjects/41598_2013_Article_BFsrep01067_Equ15_HTML.jpg)
where represent the number of driver nodes for the degree-preserved randomized version of the original network. Hence if Δ = 0 then P(kin, kout) accurately determines ND; if Δ ≠ 0 then the structural properties not captured by the degree sequence influence its controllability. We measure the correlations in several real networks and based on our numerical and analytical results we predict the sign of Δ (Fig. 7). We grouped the networks according to our predictions. We provide the details of each network dataset in the Supplementary Information Table SI.
Group A
The networks of p2p Internet (Gnutella filesharing clients) do not have strong correlations, therefore we expect nD to be correctly approximated by the prediction based on P(kin, kout) (i.e. Δ ≈ 0), in line with the empirical observations.
Group B
As in most networks the three relevant correlations coexist to some degree (Fig. 7), it is impossible to isolate their individual role. Yet, the networks in this group (electric circuits, metabolic networks, neural networks, power grids and food webs with exception of the Seagrass network) all have negative out-in and nonzero in-in and out-out correlations, each of which individually increase nD as we showed above. Therefore we predict Δ > 0, in line with the empirical observations.
Group C
Only the prison social-trust and the cell phone network feature significant positive out-in correlations. These networks also display nonzero in-in and out-out correlation, leading to the coexistence of two competing effects: out-in correlations decrease nD and the out-out and in-in correlations increase nD. Since the out-in correlation is a first order effect (equation (8)), while out-out and in-in correlations are only of second order (equation (10)), we expect a decrease in nD (i.e. Δ < 0), consistent with the empirical results.
Group D
The Seagrass food web and citation networks do not feature significant out-in correlations, only the secondary in-in and out-out correlations, hence we expect nD to increase (Δ > 0), consistent with the observations.
Group E
Only the transcriptional regulatory networks are somewhat puzzling in that they show degree correlations, yet the degree sequence still correctly gives nD. However, the simulations indicated that the effect of correlations is negligible for high nD. And our analytical results showed that the value of the correction depends on details of e(ji, jo; ki, ko), not captured by the Pearson coefficient r. These observations highlight that even though in most cases our qualitative predictions based on r are valid, in some cases further investigation is required.
Discussion
The goal of our paper was to clarify the higher order network characteristics that influence controllability. We studied the effect of three topological characteristics: clustering, modularity and degree correlations. We used numerical simulations to identify the role of the relevant characteristics, finding that changes in the clustering coefficient and the community structure have no systematic effect on the the minimum number of driver nodes nD. In contrast degree correlations showed a robust effect, whose magnitude and direction depends on the type of correlation. Using the cavity method we derived nD for networks with given degree distribution and correlation profiles, finding results that are consistent with our numerical simulations. For real networks these numerical and analytic results enabled us to qualitatively explain the deviation of the observed nD from the prediction based only on P(kin, kout).
Our results not only offer a new perspective on the role of topological properties on network controllability, but also raise several questions. Future research directions include determining the optimal network structure to minimize the number of necessary driver nodes and studying how different network characteristics influence the robustness of the control configuration.
Methods
Generating a scale-free network
We use the static model to generate directed scale-free networks40. We start from N disconnected nodes and assign a weight wi = (i + i0)–α to each node i (i = 1 … N). We randomly select two nodes i and j with probability proportional to wi and wj respectively and if they are yet not connected, we connect them. We allow self-loops, but avoid multi-edges. We repeat the process until L links have been placed. The resulting network has average degree 〈k〉 = 2L/N and P(in/out)(k) ~ k–γ for large k, where and maximum degree
.
To systematically study correlations, the starting network has to be uncorrelated. However, the presence of hubs may induce unwanted degree correlations41 and may also considerably limit the maximum and minimum correlations accessible via rewiring42. We overcome these difficulties by introducing a structural cutoff in the degrees, choosing i0 to ensure kmax < (〈k〉N)1/2 43. Note, that in the static model of Goh et al. i0 = 040.
As both in- and out-degree of node i is proportional to wi, the above procedure results in correlations between the in- and out-degrees of node i. To eliminate the correlations, we randomize the in-degree sequence while keeping the out-degree sequence unchanged.
Rewiring algorithm
We use degree preserving rewiring20 to add each network characteristic. Suppose that the chosen network characteristic is quantified by a metric X. To set its value to X*, we define the E(X) = |X − X*| energy, so E(X*) is a global minimum. We minimize this energy by simulated annealing44: (1) choose two links at random with uniform probability; (2) rewire the two links and calculate the energy E(X) of the resulted network; (3) accept the new configuration with probability
![](http://media.springernature.com/lw271/springer-static/image/art%3A10.1038%2Fsrep01067/MediaObjects/41598_2013_Article_BFsrep01067_Equ16_HTML.jpg)
where the β parameter is the inverse temperature; (4) repeat from step one and gradually increase β. Stop if |E(X) − E(X*)| is smaller than a predefined value.
Note, that keeping the degree sequence bounds the possible values of X that can be reached by rewiring. In all cases we study the full interval of accessible X values.
References
Ben-Naim, E., Frauenfelder, H. & Toroczkai, Z. (Eds.) Complex Networks (Springer, Berlin, 2004).
Newman, M., Barabási, A.-L. & Watts, D. J. The Structure and Dynamics of Networks (Princeton University Press, Princeton, 2006).
Caldarelli, G. Scale-Free Networks: Complex Webs in Nature and Technology (Oxford University Press, Oxford, 2007).
Barrat, A., Barthélemy, M. & Vespignani, A. Dynamical Processes on Complex Networks (Cambridge University Press, Cambridge, 2009).
Cohen, R. & Havlin, S. Complex Networks: Structure, Robustness and Function (Cambridge University Press, Cambridge, 2010).
Chen, G., Wang, X. & Li, X. Introduction to Complex Networks: Models, Structures and Dynamics (Higher Education Press, Beijing, 2012).
Wang, X. F. & Chen, G. Pinning control of scale-free dynamical networks. Physica A 310, 521–531 (2002).
Sorrentino, F., di Bernardo, M., Garofalo, F. & Chen, G. Controllability of complex networks via pinning. Phys. Rev. E 75, 1–6 (2007).
Gutiérrez, R., Sendiña-Nadal, I., Zanin, M., Papo, D. & Boccaletti, S. Targeting the dynamics of complex networks. Sci. Rep. 2 (2007).
Liu, Y.-Y., Slotine, J.-J. & Barabasi, A.-L. Controllability of complex networks. Nature 473, 167–173 (2011).
Wang, W.-X., Ni, X., Lai, Y.-C. & Grebogi, C. Optimizing controllability of complex networks by minimum structural perturbations. Phys. Rev. E 85, 026115 (2012).
Yan, G., Ren, J., Lai, Y.-C., Lai, C.-H. & Li, B. Controlling complex networks: How much energy is needed? Phys. Rev. Lett. 108, 218703 (2012).
Nepusz, T. & Vicsek, T. Controlling edge dynamics in complex networks. Nat. Phys. 8, 568–573 (2012). 10.1038/nphys2327.
Slotine, J.-J. & Li, W. Applied Nonlinear Control (Prentice-Hall, 1991).
Luenberger, D. G. Introduction to Dynamic Systems: Theory, Models, & Applications (John Wiley & Sons, New York, 1979).
Kalman, R. E. Mathematical description of linear dynamical systems. J. Soc. Indus. and Appl. Math. Ser. A 1, 152 (1963).
Chui, C. K. & Chen, G. Linear Systems and Optimal Control (Springer-Verlag, New York, 1989).
Lovász, L. & Plummer, M. D. Matching Theory (American Mathematical Society, Rhode Island, 2009).
Pastor-Satorras, R., Vazquez, A. & Vespignani, A. Dynamical and correlation properties of the internet. Phys. Rev. Lett. 87, 258701 (2001).
Maslov, S. & Sneppen, K. Specificity and stability in topology of protein networks. Science 296, 910–913 (2002).
Newman, M. E. J. Mixing patterns in networks. Phys. Rev. E 67, 026126 (2003).
Watts, D. J. & Strogatz, S. H. Collective dynamics of ‘small-world’ networks. Nature 393, 440–442 (1998).
Ravasz, E., Somera, A. L., Mongru, D. A., Oltvai, Z. N. & Barabási, A.-L. Hierarchical organization of modularity in metabolic networks. Science 297, 1551–1555 (2002).
Palla, G., Derényi, I., Farkas, I. & Vicsek, T. Uncovering the overlapping community structure of complex networks in nature and society. Nature 435, 814–818 (2005).
Leicht, E. A. & Newman, M. E. J. Community structure in directed networks. Phys. Rev. Lett. 100, 118703 (2008).
Fortunato, S. Community detection in graphs. Physics Reports 486, 75–174 (2010).
Barrat, A. & Weigt, M. On the properties of small-world network models. Eur. Phys. J. B 13, 547–560 (2000).
Foster, J. G., Foster, D. V., Grassberger, P. & Paczuski, M. Edge direction and the structure of networks. Proc. Natl. Acad. Sci. 107, 10815 (2010).
Erdös, P. & Rényi, A. On the evolution of random graphs. Publ. Math. Inst. Hung. Acad. Sci. 5, 17–60 (1960).
Bollobás, B. Random Graphs (Cambridge University Press, Cambridge, 2001).
Barabási, A.-L. & Albert, R. Emergence of scaling in random networks. Science 286, 509 (1999).
Albert, R. & Barabási, A.-L. Statistical mechanics of complex networks. Rev. Mod. Phys. 74, 47–97 (2002).
Hopcroft, J. E. & Karp, R. M. An n5/2 algorithm for maximum matchings in bipartite graphs. SIAM J. Comput. 2, 225 (1973).
Ravasz, E. & Barabási, A.-L. Hierarchical organization in complex networks. Phys. Rev. E 67, 026112 (2003).
Mones, E., Vicsek, L. & Vicsek, T. Hierarchy measure for complex networks. PLoS ONE 7(3): e33799 (2012).
Ahn, Y.-Y., Bagrow, J. P. & Lehmann, S. Link communities reveal multiscale complexity in networks. Nature 466, 761–764 (2010).
Zdeborová, L. & Mézard, M. The number of matchings in random graphs. J. Stat. Mech. 05, P05003 (2006).
Cowan, N. J., Chastain, E., Vilhena, D. A. & Bergstrom, C. T. Controllability of Real Networks. arXiv:1106.2573v3 (2011).
Lin, C.-T. Structural Controllability. IEEE Trans. Auto. Contr. 19, 201 (1974).
Goh, K.-I., Kahng, B. & Kim, D. Universal behavior of load distribution in scale-free networks. Phys. Rev. Lett. 87, 278701 (2001).
Boguñá, M., Pastor-Satorras, R. & Vespignani, A. Cut-offs and finite size effects in scale-free networks. The European Physical Journal B - Condensed Matter and Complex Systems 38, 205–209 (2004).
Menche, J., Valleriani, A. & Lipowsky, R. Asymptotic properties of degree-correlated scale-free networks. Phys. Rev. E 81, 046103 (2010).
Chung, F. & Lu, L. Connected components in random graphs with given expected degree sequences. Annals of Combinatorics 6, 125–145 (2002).
Press, W., Flannery, B., Teukolsky, S. & Vetterling, W. Numerical Recipes in C: The Art of Scientific Computing (Cambridge University Press, 1992).
Acknowledgements
This work was supported by the Network Science Collaborative Technology Alliance sponsored by the US Army Research Laboratory under Agreement Number W911NF-09-2-0053; the Defense Advanced Research Projects Agency under Agreement Number 11645021; the Defense Threat Reduction Agency award WMD BRBAA07-J-2-0035; FET IP project MULTIPLEX (3A532) and the generous support of Lockheed Martin. M. Pósfai has received funding from the European Union Seventh Framework Programme (FP7/2007–2013) under grant agreement No. 270833.
Author information
Authors and Affiliations
Contributions
All authors designed and did the research. M.P. analysed the empirical data and did the analytical and numerical calculations. A.-L.B. was the lead writer.
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Electronic supplementary material
Supplementary Information
Supplementary Information
Rights and permissions
This work is licensed under a Creative Commons Attribution-NonCommercial-NoDerivs 3.0 Unported License. To view a copy of this license, visit http://creativecommons.org/licenses/by-nc-nd/3.0/
About this article
Cite this article
Pósfai, M., Liu, YY., Slotine, JJ. et al. Effect of correlations on network controllability. Sci Rep 3, 1067 (2013). https://doi.org/10.1038/srep01067
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/srep01067