Universal framework for edge controllability of complex networks

Dynamical processes occurring on the edges in complex networks are relevant to a variety of real-world situations. Despite recent advances, a framework for edge controllability is still required for complex networks of arbitrary structure and interaction strength. Generalizing a previously introduced class of processes for edge dynamics, the switchboard dynamics, and exploit- ing the exact controllability theory, we develop a universal framework in which the controllability of any node is exclusively determined by its local weighted structure. This framework enables us to identify a unique set of critical nodes for control, to derive analytic formulas and articulate efficient algorithms to determine the exact upper and lower controllability bounds, and to evaluate strongly structural controllability of any given network. Applying our framework to a large number of model and real-world networks, we find that the interaction strength plays a more significant role in edge controllability than the network structure does, due to a vast range between the bounds determined mainly by the interaction strength. Moreover, transcriptional regulatory networks and electronic circuits are much more strongly structurally controllable (SSC) than other types of real-world networks, directed networks are more SSC than undirected networks, and sparse networks are typically more SSC than dense networks.

switchboard dynamics as a general mathematical framework for edge dynamics and discovered that the structural controllability of edge dynamics can differ characteristically from that of nodal dynamics. Specifically, they considered directed complex networks G(V, E) with x denoting the vector specifying the state of each edge in the network, and with − y v and + y v being the state vectors corresponding to the incoming and outgoing edges of node v, respectively. Factors that can influence the evolution of the state vector + y v are the vector − y v , the vector of the damping terms τ v , and the external input vector u v . The edge dynamical process can then be described by the following switchboard dynamics: is switching matrix. Its row number equals the out-degree + k v and its column number is the in-degree − k v of node v. σ v is unity if node v is a driver node and is zero otherwise. ⊗ denotes the entry-wise product of the two vectors of the same size. This switchboard dynamics is suitable for modeling a variety of real world situations such as social communications and load-balancing or routing on the Internet 38,39 .
Reformulating the switchboard dynamics in terms of the edge variables yields a time invariant dynamical system: where  ∈ × W M M is the transpose of the adjacency matrix of the line graph L(G) (see Fig. 1(a,b) for an example of a line graph), in which w ij is nonzero if and only if the head of edge j is the tail of edge i. T is a diagonal matrix composed of the damping terms of each edge. H is a diagonal matrix where the i th element is σ v if node v is the tail of edge i. The controllability of system (2) can be assessed and quantified by employing the structural control theory through the assumption that W − T is a structural matrix, while omitting the effect of interaction strengths in the switching matrix 18,20 . Consequently, a minimum set of driven edges in the original network G can be identified by calculating the maximum matching of the line graph L(G). All the tail nodes of the driven edges are the driver nodes in G. One key result is that all divergent nodes ( > + − k k v v ) are driver nodes, and that one arbitrary node from each balanced component ( = > for all nodes in a connected component) is also a driver node 38 . The criterion for discerning driver nodes gives rise to several structural controllability properties of edge dynamics that differ markedly from those associated with nodal dynamics. Other findings include that most of (h) Driver node, driven edge and input signals for the structural adjacency matrix W. Panels (c-e) and (f-h) correspond to the lower and upper bounds of controllability, respectively. The linearly dependent rows in W in (d,g) stem from independent rows in the switching matrices in (c,f), respectively. The edges associated with linearly dependent rows in W are the driven edges that should be controlled. The external input signals u should be imposed on the tail nodes of the driven edges.
real-world networks are more controllable than their randomized counterparts, transcriptional regulatory networks are easy to be controlled, heterogeneous networks are more controllable than homogeneous networks, and a positive correlation between the in-and out-degrees can enhance the controllability 38 . Despite the interesting findings, since the structural controllability theory is valid only for directed networks, a number of open issues remain, such as the edge controllability of undirected networks, the effect of interaction strengths on the controllability, and the strong structural controllability associated with edge dynamics. It is worth of noting that strong structural controllability is an important index for measuring the robustness of controllability against uncertainties or variations in the interaction strengths among the edges.
In this paper, we generalize the existing framework of structural edge controllability 38 by developing a universal framework capable of characterizing the controllability of edge dynamics in arbitrary networks and interaction strengths. In particular, by bridging the exact controllability theory for nodal dynamics and the general switchboard dynamics, we find that, for an arbitrary network with any distributions of the interaction strengths among the edges, the role of a node in edge dynamics and strong structural controllability is exclusively determined by the local weighted structure of the node. We use this key result to uncover a number of phenomena associated uniquely with edge controllability, which have no counterparts in nodal controllability. Firstly, the set of driver nodes in edge controllability is unique in an arbitrary network, whereas for nodal controllability, there are many configurations of driver node sets in spite of the fixed number of driver nodes. The set of strongly structurally controllable (SSC) nodes is unique and is fully determined by the local topology of each node, whereas for nodal controllability, there exists no criterion to identify the SSC nodes. Secondly, interaction strengths among the edges play a more significant role in edge controllability than the network structure does. Particularly, there exist lower and upper bounds of edge controllability, which are determined by the interaction strengths. We prove rigorously that the lower and upper bounds are determined by the structural switching and the unweighted switching matrices, respectively. In fact, there is a vast range between the bounds, in which a broad spectrum of controllability can be achieved. Thirdly, applications of our framework to real-world networks show that transcriptional regulatory networks and electronic circuits have higher strong structural controllability. In addition, real directed networks are more SSC than undirected networks, and sparser networks are more SSC than denser networks. For all the results concerning general edge controllability and strongly structural controllability, we provide analytic formula and results from extensive numerical tests. We emphasize the universal nature of our edge controllability framework: it is applicable to arbitrary network structure (e.g., directed or undirected, weighted or unweighted) and interaction strengths among the edges. In fact, we demonstrate that for directed networks with structural switching matrices, a number of key results reduce to those of structural edge controllability.

Results
General switchboard dynamics. In the original switchboard dynamics 38 , the switching matrix S v must be a structural matrix, in which all nonzero elements are independent free parameters. Instead, we release the restriction of structural matrix S v and consider a general switchboard dynamics with any kind of switching matrices, in which the elements capture the interaction strengths among the edges. We exemplify two typical switching matrices, the weighted switching matrix and the unweighted switching matrix. In the former, all nonzero elements can be any values, and in the latter, all nonzero elements are one.
For directed networks, the general switchboard dynamics (GSBD) is described by Equation (2). The adjacency matrix W of the line graph is of the same type as S v of each node (Fig. 1). In contrast, for undirected networks, Equation (2) cannot be immediately adopted because, associated with each edge, the interaction and transmission are bidirectional. Two neighboring edges connecting with the same node can be the input and output of each other. To define GSBD for undirected networks, we split each undirected edge into two directed edges of opposite directions and use a pair of state variables, ′ x x ( , ) i i , to denote such an edge, where each variable corresponds to one of the directed edges. The state vector of the dynamical process occurring on undirected edges is then . The switching matrix can be written as For the whole network, the dynamical process can still be described by Equation (2), but the dimension of the matrix W is doubled: where M is the number of undirected edges in G. The properties of W are still determined by S v (see Supplementary Fig. S1 for an illustration). Our GSBD thus provides a more general characterization of the dynamics occurring on edges for arbitrary networks. The focus of our study is on the effect of the interaction strengths in the switching matrix on edge controllability. For a general switching matrix S v , the conventional structural control theory 20 is not applicable, due to the non-uniform interaction strengths among the edges and the undirected nature of the network structure. This calls for a more general control theory to determine/quantify the edge controllability.
Controllability framework of the GSBD. We make use of the exact controllability theory 21 developed recently to determine the controllability of the GSBD. By definition, the controllability of a network is measured by the minimum number N D of driver nodes. Prior to the identification of the driver nodes, we must ascertain the minimum number M D of the driven edges as determined by the matrix W − T in Equation (2). According to the exact controllability theory 35 , the damping matrix T with identical diagonal elements has no effect on the controllability of the network characterized by W − T, which can be proven rigorously. That is, the set of driven edges and driver nodes will not be affected by T, so it can be neglected. As a result, all self-loops of nodes stemming from −T in the line graph L(G) can be eliminated (see Fig. 1). We then determine M D of G characterized by the adjacency matrix W and N D of L(G), and present our key results. Figure 1 shows two representative cases of a simple directed network with structural and unweighted switching matrices. We can prove that the former and the latter cases generate the lower and upper bounds, respectively, of both M D and N D for any network (see Methods). The line graph L(G) is shown in Fig. 1(b). The driven edges in G correspond to the driver nodes in L(G). According to the exact controllability theory and the properties of line graphs, we can prove that the minimum number of driver nodes in L(G) (driven edges in G) is M − rank(W), where M is the number of edges in G (see Supplementary Note 1). That is, M D is the number of linearly dependent rows in W. For example, in the structural matrix W of Fig. 1(d), the row corresponding to x 5 is linearly dependent on the other rows, indicating that a control signal should be applied to x 5 . Figure 1(g) illustrates linearly dependent rows in the unweighted matrix W. Making use of a generic feature (see Methods and Supplementary Note 2) of line graphs, we obtain our first key result: the linearly dependent rows in W are exclusively determined by the linearly dependent rows in the switching matrices S i of all nodes. As shown in Fig. 1, the linearly dependent row (x 5 ) in W ( Fig. 1(d)) stems from S d (Fig. 1(c)) with a null set. The other rows in W are linearly independent because S a , S b and S c in Fig. 1(c) are row-full rank. Similarly, for the unweighted W in Fig. 1(g), the two linearly dependent rows x 1 and x 5 originate from S a and S d in Fig. 1(f), respectively. This key finding indicates that driven edges can be identified from the switching matrices S i by using the local information of nodes without relying on the line graph L(G), which gives where C is the number of connected components in G. β i = 1 if the switching matrices of all nodes in component i are square matrices with full rank, and β i = 0 otherwise. The second term from each connected component has little effect on M D . As a result, M D is determined by the rank of S i . After the driven edges are determined, we can immediately specify the driver nodes at the tail end of the driven edges, ensuring that external control signals applied to the driver nodes can directly pass on to the driven edges. As shown in Fig. 1(e,h), the driver nodes at the tail end of the driven edges can be identified based on local information contained in the switching matrices in Fig. 1(c,f), which does not depend on the global structure of L(G). In general, according to the characteristics of the driven edges and the relation between the driven edges and the driver nodes, we can prove that S i associated with driver node i satisfies This means that, if the switching matrix of a node is not full row-rank, the node must be a driver node. This is the general criterion for identifying the driver nodes based only on the local information of each node. From this result, the number N D of the driver nodes can be calculated through where β i is the same as that in Equation (3). Analogous to the expression of M D , the second term in Equation (5) has little effect on N D . We refer to full row-rank as full rank in the remaining paper for simplicity. It is worth noting that our edge controllability theory is not a trivial application of the exact controllability theory to the line graphs of a network. In particular, by exploiting the unique properties of a line graph, we prove that the controllability of a node is determined only by its local weight structure. As a result, the set of driver nodes in a general network with an arbitrary distribution of interaction strengths among the edges is unique. This is our key result here. Moreover, we also prove that, for the upper and lower controllability bounds, the driver nodes are determined only by the local topology of each node, enabling readily implementable algorithms to find the driver nodes associated with the bounds. As presented below, this key result can also address the issue of strong structural controllability, which is important for understanding the robustness of controlling networks in the presence of uncertainties or variations of interaction strengths among the edges.
Universal controllability bounds. We can prove that there exist universal upper and lower bounds for edge controllability for any network, and that any value of the controllability in between can be achieved by adjusting the interaction strengths among the edges. In particular, the upper and lower bounds are reached if S i of each node is an unweighted matrix and a structural matrix, respectively (see Methods for a detailed proof).
In general, it is necessary to calculate the rank(S i ) of all nodes to identify the driver nodes and to obtain N D and M D . However, for structural or unweighted S i , we are able to identify the driver nodes and the driven edges based solely on their in-degrees − k v and out-degrees + k v without having to calculate rank(S i ). As shown in Fig. 2, nodes in terms of their in-and out-degrees can be classified into three categories: For case (I), e.g., node a in Fig. 2(a), the structural S a is always full rank ( Fig. 2(b)), indicating that node a is non-essential in the sense that no external input signal is needed to control the outgoing edges of a, and that all the outgoing edges of a are non-essential as well. In contrast, for the unweighted S a , rank(S a ) is always unity and it satisfies the inequality < + S k rank( ) a a , as shown in Fig. 2(c). Thus, without any external input signal applied to the node, only one outgoing edge can be fully controlled, regardless of the degree of a. As a result, a in this case must be a driver node. Moreover, its one outgoing edge that can be arbitrarily selected is a non-essential edge, while the other outgoing edges are driven edges. For case (II), e.g., node b in Fig. 2(d), the structural S b is always full rank ( , indicating that node b and all of b's outgoing edges are non-essential. This is similar to the structural matrix S a in case (I), as shown in Fig. 2(e). In contrast, for the unweighted S b , despite the inequality > − + k k b b , the rank of S b is always unity and b is a driver node with a single non-essential edge among all b's outgoing edges, as shown in Fig. 2(f). This also implies that, for an unweighted switching matrix, the associated node must be a driver node if > For case (III), e.g., node c in Fig. 2(g), both the structural S c (Fig. 2(h)) and the unweighted S c (Fig. 2(i)) satisfy the inequality < + S k rank( ) c c . As a result, node c is a driver node for either structural or unweighted switching matrix, with different driven edges for the two scenarios.
Scientific RepoRts | 7: 4224 | DOI:10.1038/s41598-017-04463-5 From the above illustrations and arguments, we obtain our key result of identifying the driver nodes based on their + k v and − k v for the lower (structural switching matrix) and the upper (unweighted switching matrix) bounds. In particular, for a directed network, at the lower bound, a node with more outgoing than incoming edges must be a driver node. For the upper bound, a node with more than one outgoing edge must be a driver node, and a node without incoming edges and with a single outgoing edge must also be a driver node. For an undirected network, at the lower bound, a single driver node is required for each connected component, while for the upper bound, a node with more than one edge must be a driver node.
According to the criterion, for identifying the driver nodes and the driven edges, N D and M D associated with the lower and upper bounds can be calculated, as summarized in Table 1, where the contributions from isolated components are included (see Methods and Supplementary Note 3). Note that, for the lower bound in directed networks, our results reduce to those of structural edge controllability (see Supplementary Note 4). For the other scenarios, the results have not been reported prior to our work to our knowledge.
Next we will verify the universal controllability bounds for both directed and undirected model networks and offer further analytical results. According to Liu et al. 20 , the nodal controllability n D can be defined as the ratio of S a is row-full rank, so the two outgoing edges are non-essential (ordinary) edges and a is a non-essential node. (c) Unweighted switching matrix S a and the category of node a and its outgoing edges. In S a , there is a linearly dependent row corresponding to a driven outgoing edge. Node a at the tail end of the driven edge is a driver node. (d) Node b with three incoming edges and two outgoing edges. (e) Structural switching matrix S b with row-full rank and the category of node b and its outgoing edges. Node b and its outgoing edges are nonessential. (f) Unweighted switching matrix S b and the category of node b and its outgoing edges. The row-rank of S b is unity, leading to one driven edge. Node b becomes a driver node with one driven edge. (g) Node c with two incoming edges and three outgoing edges. (h) Structural switching matrix S c with deficient row-rank and the category of node b and its outgoing edges. The row-rank of S c is 2, so there are two non-essential edges. Node c is a driver node with one driven edge. (i) Unweighted switching matrix S c and the category of node c and its outgoing edges. The row-rank of S c is unity, indicating one non-essential edge and two driven edges. Node c with two driven outgoing edges is thus a driver node. The structural switching matrices in (b), (e,h) correspond to the lower bound, and the unweighted switching matrices in (c), (f,i) are associated with the upper bound.  Another remarkable result is that, except for very small values of the average degree 〈k〉, the range or "distance" between the upper and lower bounds is appreciable. Any value of n D and m D in between the bounds is achievable by properly setting the element values in the switching matrices S i . These results demonstrate that interaction strengths among the edges in the switching matrices play a more important role in edge controllability than the network structure does, in sharp contrast to the situation of controlling nodal dynamics (or nodal controllability). These findings provide a deeper understanding of the controllability of edge dynamics.

Directed Undirected
We also find a non-monotonic behavior in the upper bound of m D for both directed and undirected networks, and in the lower bound of n D for undirected networks as well. Such a behavior results from the combining effect of the first and the second terms in Equations (3) and (5), where the second term represents the contribution from each isolated component. Notice that the non-monotonic phenomenon occurs in the regime of relatively small values of 〈k〉, for which there are a number of isolated components. As 〈k〉 increased, the contribution from the second term diminishes because of the reduction in the number of isolated components, whereas the contribution from the first term begins to dominate, providing an explanation for the non-monotonic behavior.
We derive the analytical results of the bounds for different networks according to the classification of nodes and edges from local information, as shown in Fig. 2. The results are presented in Table 2 where, except for m D of directed networks associated with the lower bound, closed-form formulas can be obtained (see Supplementary Note 5 for detailed derivations). For directed networks, the analytical results of the lower bound reduce to the previous results of structural edge controllability, providing further validation of our theory.  Table S1. The transition between the upper and lower bounds has also been analyzed, as presented in Supplementary Note 6 and Fig. S4.

Strong structural controllability. Strong structural controllability is a critical notion of quantification if
the controllability of a network is robust against uncertainties or fluctuations in the interaction strengths. At present, to establish strong structural controllability even for nodal dynamics remains to be a challenging problem. Remarkably, we find that, for edge dynamics, strong structural controllability can be related to the controllability bounds in a straightforward manner, thus allowing us to obtain a straightforward but appealing analytic criterion to quantify strong structural controllability.
A network is SSC if, regardless of the values of the elements in the switching matrices, the quantities n D and m D do not change 20,40 . This means that the controllability of a fully SSC network cannot be affected by variations in the interaction strengths among the edges, exclusively determined by the network structure. Because the category of a node is determined only by its local information (Fig. 2), we can determine if a node is SSC based only on its switching matrix. In particular, a node is SSC if the nodal and edge categories do not change for any values of the elements in its switching matrix. We can prove that a node with ≤ is SSC for an arbitrarily directed network and a node with k v ≤ 1 is SSC for an arbitrary undirected network (see Methods).
The strong structural controllability of a network can be defined as the ratio of the number N SSC of the SSC nodes to the network size N, i.e., ssc ssc A network with a higher value of n ssc is more SSC, and a network with n ssc = 1 is fully SSC. As shown in Fig. 3(e,f), sparse networks with small values of 〈k〉 are nearly fully SSC. This is also reflected in Fig. 3(a-d) in the regime of small 〈k〉, where the upper and lower bounds are equal. Observe that n ssc decreases as 〈k〉 is increased, which means that sparser networks are generally more SSC than denser networks. The dependence of n ssc on 〈k〉 can be analytically predicted (Supplementary Note 7). Strong structural controllability is related with the controllability bounds, in the sense that the upper and lower bounds of a fully SSC network coincide with each other. This can be explained, as follows. The controllability bounds are determined by the interaction strengths. For a fully SSC network, its edge controllability is not affected by the interaction strengths. That is, interaction strengths do not induce any difference between the upper and lower bounds in an SSC network. As a result, the controllability bounds must be exactly the same in a fully SSC network.

Controllability properties of real-world networks.
Our theoretical framework and analytic predictions enable us to study the edge controllability of a variety of real directed and undirected networks. The upper and lower bounds of n D , m D and n ssc for different types of real networks are shown in Table 3. An interesting finding is that electronic circuits 41 and regulatory networks are more SSC, including the ownership network of US telecommunications and media corporations 42 and the transcriptional regulatory networks of Echerichia Coli 41 and Saccharomyces cerevisiae 41,43 . The values of n ssc of all the networks belonging to the two types are higher than 0.9, which are substantially larger than the mean value of n ssc of other networks. Consistent with the results of (1 ) Table 2. Analytical results of controllability for upper and lower bounds. Formulas for n D U and m D U (upper bounds) and n D L and m D L (lower bounds) for directed and undirected ER and SF (κ → ∞) networks. The average degree is 〈k〉 = 〈k in 〉 = 〈k out 〉 = M/N, I a (x) is the modified Bessel function of the first kind, ζ(x) is the Riemann zeta function, and n CC is the expected fraction of isolated components in the ER undirected networks. Note that the average degree is 〈k〉 = ζ(γ)/ζ(γ − 1) for SF networks with parameter κ → ∞ (see Supplementary Note 5 for details).
Scientific RepoRts | 7: 4224 | DOI:10.1038/s41598-017-04463-5 structural edge controllability, we find that regulatory networks are not only well-controllable under the edge dynamics but their controllability is also robust against uncertainties and variations in the interaction strengths. Except for the two types of networks, the effect of interaction strengths among the edges can be significant in other real world networks. For those networks, it is necessary to consider interaction strengths, in addition to the network structure, to gain insights into the edge controllability, especially when the networks are undirected.
Another finding for real-world networks is that the values of the lower bound n D L for undirected networks are much smaller than those for directed networks, in accordance with the results from the model networks. This indicates that real undirected networks have a much higher potential to be optimized for control through adjusting the interaction strengths. However, the values of n ssc for undirected networks are generally lower than those for directed networks, suggesting that undirected networks are more sensitive to perturbations of the interaction strengths. A trade-off thus exists between the controllability and sensitivity for undirected networks. The correspondence of the value of n ssc to the gap size between the upper and lower bounds in model networks also exists   Table 3. Controllability of edge dynamics in real networks. For each network, its type, class, name, number N of nodes, number M of edges, the upper bounds (n D U and m D U ), the lower bounds (n D L and m D L ), and strong structural controllability measure n ssc are shown. For data sources and references, see Supplementary Table S2.
We also study the dependence of the controllability bounds of the real-world networks on the average degree 〈k〉. As shown in Fig. 4(a-d), the values of n D U and n D L for real-world networks with different values of 〈k〉 indeed constitute the universal controllability bounds. Due to the small influence of network topology on controllability (Fig. 3), the values of n D U and n D L of real directed networks tend to spread in two relatively small regions, which is analogous to the behavior in model networks (Fig. 3). In contrast, the influence of the network topology on the bounds of undirected networks is much weaker than for directed networks, leading to much smaller variance about the theoretical estimation. Given the degree distribution, the controllability bounds can be predicted more precisely by eliminating the effect of the degree distribution. For each network, a good agreement is obtained between the analytical prediction and the simulation results, as shown in Fig. 4(e-h). Figure 5 shows the dependence of n scc on 〈k〉 for real networks. We find that the values of n scc are primarily determined by 〈k〉, especially for real undirected networks, in which the degree distribution has an insignificant effect on n scc . In fact, the role played by the distribution is marginal in our theoretical prediction of n scc for model networks with an exponential degree distribution. Nonetheless, for real-world networks, their degree distribution can be taken into account and we obtain a reasonable agreement between the theoretically predicted and numerical values of n scc for both directed and undirected networks (see Supplementary Fig. S5 for m D and Note 8 for the analytical results). We find that the effect of degree distribution on n ssc is quite similar to that on the controllability bounds.

Discussion
Most existing frameworks of controllability of complex networks focus on nodal dynamical processes. In real-world networks, however, edge dynamics can also be important, such as the Internet, transportation and modern social networks. Network controllability based on edge dynamics was first addressed by Nepusz and Vicsek 38 , who introduced a mathematical framework to describe generic edge processes through switchboard dynamics, and uncovered significant differences between network nodal and edge controllability. However, their framework, being fundamentally a structural controllability framework, was limited to directed networks. The goal of our work is to develop an edge controllability framework that is universally applicable to all types of complex networks, directed or undirected, weighted (with arbitrary distributions) or unweighted. To accomplish this goal, we propose a class of generalized switchboard dynamics and exploit the exact controllability theory 21 for complex networks.
Comprehensive mathematical analyses and extensive numerical tests with model and real-world networks have revealed a number of striking properties associated with edge controllability. For example, it is exclusively determined by the local topology and the interaction strengths, and this holds generally for arbitrary networks with any distributions of interaction strengths. This result provides a unique configuration of driver nodes, in stark contrast to what the nodal controllability theory can offer where it usually yields many configurations of the driver nodes. More strikingly, our framework is capable of providing a unique configuration of SSC nodes, a notion that has significant applied value for network control, whereas no structural controllability theory is able to offer a criterion for determining the SSC nodes. Our framework also allows us to address the challenging issue of upper and lower controllability bounds, and we prove that only given the local topology of each node, the controllability bounds of each node and those of the whole network can be completely determined. This result enables analytical predictions of the bounds and the SSC property based solely on the degree distribution. We also prove that the upper and lower bounds correspond to unweighted and structural switching matrices, respectively. A finding is that, for a fully SSC network, the gap between these bounds must vanish. In general, the interaction strengths play a more significant role in edge controllability than the network structure does, due to the typically large range between the controllability bounds, where an arbitrary degree of controllability in between the bounds can be achieved by adjusting the interaction strengths.
Applying our universal edge controllability framework to a large number of real-world networks, we find that transcriptional regulatory networks and electronic circuits possess the highest strong structural controllability. We also find that directed networks in the real world tend to be more SSC than undirected networks, and that sparser networks are usually more SSC than denser networks.
Our work raises a number of open questions, the answers to which would further deepen our understanding of the controllability of real world complex networks. For example, can a method be develop to implement target control for the general edge dynamics? What is the effect of the correlation between in-degrees and out-degrees on the edge controllability? Would it be possible to realize partial control of a subset of edges from a minimum number of driver nodes? Are there any approaches to optimize edge controllability through small perturbation? What is the energy cost in controlling general edge dynamics? Is it possible to treat nonlinear edge dynamics? Taken together, the framework developed here provides a base to address these questions, opening a new avenue towards fully controlling real networked systems in a wide range of fields.

Methods
Relationship between switching matrix S v and adjacency matrix W for line graphs. In GSBD, all nonzero elements in the adjacency matrix W of the line graph L(G) come from the switching matrices in the original directed or undirected network G, and the nonzero elements in the identical columns (rows) stem from the same switching matrix. An example is shown in Fig. 1, where the nonzero elements in W of L(G) (Fig. 1(d,g)) correspond to the elements with the same color in the switching matrices of the original network ( Fig. 1(c,f)). Furthermore, according to the generic property of directed line graphs, any two columns (rows) of the adjacency matrix W are either identical or orthogonal to each other. As a result, the contribution of a switching matrix S v to the rank of W is the rank of S v , leading to the following relation between the ranks of the switching matrices and the adjacency matrix: where N is the number of nodes in the original network G, and the relation is rigorous for any networks (see Supplementary Note 2). The relation in combination with the exact controllability theory 21 leads to our general formulas for M D (Eq. (3)) and N D (Eq. (5)).

Controllability bounds of nodes and edges.
The upper and lower bounds of controllability exist generally in a network. They can also be defined for each node in the network, as determined by the interaction strengths between the incoming and outgoing edges of the node. Because controllability is determined by the switching matrix of each node, we can ascertain the existence of the upper and lower bounds for any network by proving the existence of the bounds for each node. Specifically, for a node with > + k 0 v and > − k 0 v , the maximum rank of its switching matrix is equal to the smaller value of the numbers of rows ( + k v ) and columns ( − k v ), which corresponds to that of the structural switching matrix. In addition, the minimum rank of the switching matrix is S v = 1 if all rows (columns) are linearly dependent, which can be achieved in an unweighted switching matrix with identical element values. Moreover, the upper and lower bounds of a network can be reached if the corresponding bounds of every node are realized. We can thus conclude that the upper and lower bounds of M D and N D of a network correspond to unweighted and structural switching matrices, respectively. Driven edges and driver nodes associated with multiple network components. In a network with a single connected component, the driver nodes and the driven edges associated with the upper and lower bounds can be identified based on local topological information of nodes. An extreme case is where all nodes and edges are non-essential. In this case, we can randomly select a node to be a driver node on which an input signal is imposed, and the driver node is determined by the second term in Equation (5). In most cases, the contribution of the second term for a single connected component is negligible, compared with that of the first term. For a network consisting of a number of isolated components, the effect of each connected component on N D must be taken into account, especially for undirected networks. In particular, for an undirected network with structural switching matrices, each connected component requires a driver node because every node is non-essential in the component. In this case, the second term in Equation (5) cannot be neglected. Table 1 lists the complete formulas of N D and M D in terms of both local topology and multiple components.

Identification of SSC nodes.
A node is either SSC or weakly structurally controllable (WSC). For a WSC node, element values in its switching matrix S v will affect the rank of S v and the category that the node belongs to. For an SSC node, its lower and upper bounds coincide with each other so that its category is determined only by its local structure. As a result, whether a node is SSC can be discerned in terms of its in-and out-degrees.
More specifically, for a node with , any change in the value of the nonzero elements will not affect the rank of S v , because S v is always unity, i.e., the minimum rank. Thus, in this case, the category to which the node belongs will not be affected by the element values in S v . For a node with > Networks analyzed. The model networks used in this paper are the Erdös-Rényi random networks 44 and scale-free networks 45 . How to generate model networks is provided in Supplementary Note 9. The real networks are described in Supplementary Table S2.