Network controllability

Controlling a simple network.

Network Controllability is concerned about the structural controllability of a network. Controllability describes our ability to guide a dynamical system from any initial state to any desired final state in finite time, with a suitable choice of inputs. This definition agrees well with our intuitive notion of control. The controllability of general directed and weighted complex networks has recently been the subject of intense study by a number of groups in wide variety of networks, worldwide. Recent studies by Sharma et al on multi-type biological networks (gene- gene, miRNA- gene , Protein-protein interaction networks) identified control targets in phenotypically characterized Osteosarcoma showing important role of genes and proteins responsible for maintaining tumor microenvironment.

Background

Consider the canonical linear time-invariant dynamics on a complex network ${\dot {\mathbf {X} }}(t)=\mathbf {A} \cdot \mathbf {X} (t)+\mathbf {B} \cdot \mathbf {u} (t)$ where the vector $\mathbf {X} (t)=(x_{1}(t),\cdots ,x_{N}(t))^{\mathrm {T} }$ captures the state of a system of $N$ nodes at time $t$. The $N\times N$ matrix $\mathbf {A}$ describes the system's wiring diagram and the interaction strength between the components. The $N\times M$ matrix $\mathbf {B}$ identifies the nodes controlled by an outside controller. The system is controlled through the time dependent input vector $\mathbf {u} (t)=(u_{1}(t),\cdots ,u_{M}(t))^{\mathrm {T} }$ that the controller imposes on the system. To identify the minimum number of driver nodes, denoted by $N_{\mathrm {D} }$, whose control is sufficient to fully control the system's dynamics, Liu et al. attempted to combine the tools from structural control theory, graph theory and statistical physics. They showed that the minimum number of inputs or driver nodes needed to maintain full control of the network is determined by the 'maximum matching’ in the network, that is, the maximum set of links that do not share start or end nodes. From this result, an analytical framework, based on the in-out degree distribution, was developed to predict $n_{\mathrm {D} }=N_{\mathrm {D} }/N$ for scale-free and Erdős–Rényi Graphs. However, more recently it has been demonstrated that network controllability (and other structure-only methods which use exclusively the connectivity of a graph, $\mathbf {A}$, to simplify the underlying dynamics), both undershoot and overshoot the number and which sets of driver nodes best control network dynamics, highlighting the importance of redundancy (e.g. canalization) and non-linear dynamics in determining control.

It is also notable, that Liu's et al formulation would predict same values of ${n_{\mathrm {D} }}$ for a chain graph and for a weak densely connected graph. Obviously, both these graphs have very different in and out degree distributions. A recent unpublished work, questions whether degree, which is a purely local measure in networks, would completely describe controllability and whether even slightly distant nodes would have no role in deciding network controllability. Indeed, for many real-word networks, namely, food webs, neuronal and metabolic networks, the mismatch in values of ${n_{\mathrm {D} }}^{real}$ and ${n_{\mathrm {D} }}^{\mathrm {rand\_degree} }$ calculated by Liu et al. is notable. If controllability is decided mainly by degree, why are ${n_{\mathrm {D} }}^{real}$ and ${n_{\mathrm {D} }}^{\mathrm {rand\_degree} }$ so different for many real world networks? They argued (arXiv:1203.5161v1), that this might be due to the effect of degree correlations. However, it has been shown that network controllability can be altered only by using betweenness centrality and closeness centrality, without using degree (graph theory) or degree correlations at all.

A schematic diagram shows the control of a directed network. For a given directed network (Fig. a), one calculates its maximum matching: a largest set of edges without common heads or tails. The maximum matching will compose of a set of vertex-disjoint directed paths and directed cycles (see red edges in Fig.b). If a node is a head of a matching edge, then this node is matched (green nodes in Fig.b). Otherwise, it is unmatched (white nodes in Fig.b). Those unmatched nodes are the nodes one needs to control, i.e. the driver nodes. By injecting signals to those driver nodes, one gets a set of directed path with starting points being the inputs (see Fig.c). Those paths are called "stems". The resulting digraph is called U-rooted factorial connection. By "grafting" the directed cycles to those "stems", one gets "buds". The resulting digraph is called the cacti (see Fig.d). According to the structural controllability theorem, since there is a cacti structure spanning the controlled network (see Fig.e), the system is controllable. The cacti structure (Fig.d) underlying the controlled network (Fig.e) is the "skeleton" for maintaining controllability.

Structural Controllability

The concept of the structural properties was first introduced by Lin (1974) and then extended by Shields and Pearson (1976) and alternatively derived by Glover and Silverman (1976). The main question is whether the lack of controllability or observability are generic with respect to the variable system parameters. In the framework of structural control the system parameters are either independent free variables or fixed zeros. This is consistent for models of physical systems since parameter values are never known exactly, with the exception of zero values which express the absence of interactions or connections.

Maximum Matching

In graph theory, a matching is a set of edges without common vertices. Liu et al. extended this definition to directed graph, where a matching is a set of directed edges that do not share start or end vertices. It is easy to check that a matching of a directed graph composes of a set of vertex-disjoint simple paths and cycles. The maximum matching of a directed network can be efficiently calculated by working in the bipartite representation using the classical Hopcroft–Karp algorithm, which runs in O(E√N) time in the worst case. For undirected graph, analytical solutions of the size and number of maximum matchings have been studied using the cavity method developed in statistical physics. Liu et al. extended the calculations for directed graph.

By calculating the maximum matchings of a wide range of real networks, Liu et al. asserted that the number of driver nodes is determined mainly by the networks degree distribution $P(k_{\mathrm {in} },k_{\mathrm {out} })$. They also calculated the average number of driver nodes for a network ensemble with arbitrary degree distribution using the cavity method. It is interesting that for a chain graph and a weak densely connected graph, both of which have very different in and out degree distributions; the formulation of Liu et al. would predict same values of ${n_{\mathrm {D} }}$. Also, for many real-word networks, namely, food webs, neuronal and metabolic networks, the mismatch in values of ${n_{\mathrm {D} }}^{real}$ and ${n_{\mathrm {D} }}^{\mathrm {rand\_degree} }$ calculated by Liu et al. is notable. If controllability is decided purely by degree, why are ${n_{\mathrm {D} }}^{real}$ and ${n_{\mathrm {D} }}^{\mathrm {rand\_degree} }$ so different for many real world networks? It remains open to scrutiny whether control robustness" in networks is influenced more by using betweenness centrality and closeness centrality over degree (graph theory) based metrics.

While sparser graphs are more difficult to control, it would obviously be interesting to find whether betweenness centrality and closeness centrality or degree heterogeneity plays a more important role in deciding controllability of sparse graphs with almost similar degree distributions.

Control of composite quantum systems and algebraic graph theory

A control theory of networks has also been developed in the context of universal control for composite quantum systems, where subsystems and their interactions are associated to nodes and links, respectively. This framework permits to formulate Kalman's criterion with tools from algebraic graph theory via the minimum rank of a graph and related notions.

See also

• Controllability Gramian

External links

• The network controllability project website
• The video showing network controllability