Interdependent networks as an important structure of the real system not only include one-to-one dependency relationship but also include more realistic undirected multiple interdependent relationship. The study on interdependent networks plays an important role in designing more resilient real systems. Here, we mainly focus on the case of interdependent networks with a multiple-to-multiple correspondence of interdependent nodes by generalizing feedback and nonfeedback conditions. We develop a new mathematical framework and study numerically and analytically the percolation of interdependent networks with partial multiple-to-multiple dependency links by using percolation theory. By analyzing the process of cascading failure, the size of the giant component and the critical threshold are analytically obtained and testified by simulation results for couple Erdös-Rnyi and scale-free networks. The results imply that the system shows a discontinuous phase transition for different coupling strengths. We find that the system becomes more resilient and easy to defend by increasing the coupling strength and the connectivity density. Our model has the potential to shed light on designing more resilient real-world dependent systems.
The interdependent network as an important research field is attracting a lot of attention, but most of the previous works are limited to a one-to-one dependency relationship. The phenomenon of the dependent network with a multiple-to-multiple dependency relationship exists widely in the real world. The existing conclusions are not suitable to explain the resilience among more complicated multiple dependency relationships and also challenge the existing theoretical methods. To fill this gap, we develop a framework with two partially dependent networks with multiple-to-multiple dependency relationship to study the robustness of the system. We investigated the phase transition behaviors and testified the simulation results by using the percolation theory. Interestingly, the system shows a discontinuous phase transition behavior for strong coupling strength, which is a remarkable contrast to two interacting network models. The results help us to obtain a deeper understanding of the dependent system with multiple-to-multiple relationship, and the developed theoretical method could also be applied for designing more resilient real systems.
I. INTRODUCTION
A complex network as a momentous research field is attracting together researchers from different scientific communities working on these areas. Recently, interdependent networks with different types of dependent links have drawn a lot of attention.1–11 The initial investigations of network properties are mainly on single networks, which are viewed as independent and do not exchange any information with other networks.12–21 However, the real systems are always coupled or are interdependent on each other to function. More recently, some analytical frameworks for exploring the percolation properties of two interdependent networks or multilayered interdependent networks have been put forward because many real systems interact with each other virtually.22–27 Previous studies have shown two strongly interdependent networks consisting of two networks with one-to-one correspondence, where distinct phase transition and discontinuous phase transition occur, from a single network, which shows a continuous transition.22 Buldyrev et al. developed a framework for understanding the robustness of two dependent networks with a one-to-one correspondence of dependent nodes subject to such cascading failures.27 Gao et al. introduced a more general case of percolation of interdependent networks, Network of Networks, by considering undirected one-to-one correspondence with feedback and nonfeedback conditions between networks, respectively.28 Shao et al. presented a model of cascading failure on interdependent networks by considering directed support-dependence relations between two coupled network systems.3 Dong et al. introduced a coupling network model with similar types of connectivity-link within and between communities’ structure. The results highlight that the coupling community structure can significantly affect the resilience of the system in that it removes the phase transition present in a single module, and the network remains resilient at this transition.29 However, a node within the network not only depends on one node in another network but most often mutually depends on multiple nodes in the real system. Examples of this case include global trade coupled networks, power grid networks, and communication networks. Due to a shortage of resources in various countries, there exist trade dependencies between countries. Since each national demand is different and necessary, multiple-to-multiple dependent relationship can be described between countries. In addition, based on economic considerations, one power station often provides power to several communication stations, and one communication station functions by depending on multiple power stations. The undirected one here means that there is a demand for each system in real systems. This real undirected multiple dependent relationship is different from directed multiple support-dependence relations, where active nodes have at least one directed support node in the system. For an interdependent network with two networks and , a single node in network depends on multiple nodes in network and functions since at least one of the dependent nodes of another network is functional. Similarly, a node in network also depends on several nodes in network . In an undirected multiple dependence relationship, node failures in one network will cause the failure of interdependent nodes in the other network and vice versa. Furthermore, this recursive process leads to a cascading failure until a stable state is reached, in which no node fails. Based on this motivation, we develop a theoretical framework for understanding the robustness of the interdependent networks with a multiple-to-multiple correspondence of interdependent nodes. Furthermore, the system threshold value, which plays an important role in determining the system robustness, is also numerically predicted by using the percolation theory.
II. MODEL DESCRIPTION
We simply assume that two independent networks and with the same number of nodes , which have the given degree distributions and . Here, and are the degrees of a node in single network and , respectively. The topological structure of undirected multiple-to-multiple dependency links between networks and is described by degree distributions and , which means that a node in network has dependency links from network and a node in network has support links from network . The process of the cascading of failures for small networks is described in Fig. 1. Without loss of generality, we assume that the attack on network occurs before network . At the initial state, let every node in the interdependent system be active. At every stage, the nodes without any dependency link from the other network or separated from the giant component within the network are considered to be inactive and removed from the interdependent system. This process will be continuous until there is no further node removed in the system. We first consider the process of failure in network . Node becomes inactive because of the random attack, and nodes and are removed due to the absence of one interdependency link and not belonging to the set of the giant component, respectively. Since there are failures of nodes in network , some nodes in network are triggered to fail. Correspondingly, nodes and become inactive due to loss of all interdependent links after fails in the interdependent system. After the cascading failure, the system achieves a stable state and no more nodes fail.
Demonstration of the process of cascading failures of dependent networks and with undirected multiple-to-multiple dependency links. Curves represent connectivity-links within networks and dashed lines represent different multiple dependency links between networks (blue) and (orange).
Demonstration of the process of cascading failures of dependent networks and with undirected multiple-to-multiple dependency links. Curves represent connectivity-links within networks and dashed lines represent different multiple dependency links between networks (blue) and (orange).
III. ANALYTICAL FRAMEWORK
Here, we develop an analytical framework to study the system robustness by applying the method of generating functions19,30 for partially interdependent networks. Let and be the possibility that a randomly chosen connectivity-link in network is connected to the giant component of networks and , respectively. The symbol denotes the probability of a node in network with in-degree and connecting by dependency nodes within the giant component of the network . Analogously, we define probability for the dependency networks and . For partially interdependent networks and , satisfies
where we assumed that only a fraction of nodes within networks and depends on each other. and are the generation function and branch generation function of network respectively. Similarly, we can get the expression of and . represents the probability of a randomly selected link attaching a node with in-degree in network excluded in the largest component of network . represents the probability of the node with in-degree and has interdependent nodes in network , which does not belong to the giant component of network . can be explained similarly as . Similarly, we obtain the expression of with the same method as follows:
The probability of a node in network with in-degree connecting via an interdependency link to a node, which included in the giant component from network , using the following formula can be obtained:
The probability can be similarly written as
Furthermore, and , which describe the probability that a chosen node belongs to the giant component of the steady networks and at the stable state, can be obtained as
IV. COMPARING THEORETICAL RESULTS WITH SIMULATION RESULTS
In this section, we study two interdependent ER networks, whose degree distribution follow the Poisson distribution () with the equal average degree . Equations (1)–(5) become
From Fig. 2, one can observe that the simulation results agree well with the theoretical predictions and the system shows a discontinuous phase transition. As increases, and gradually decrease to zero at the critical threshold . As shown in Fig. 2(a), the system becomes more and more vulnerable and easy to defend as the average degree increases. Also, for strong coupling strength , the system becomes more robust, but for weak coupling strength , the system becomes more vulnerable. This follows the fact that fraction of nodes will fail due to the absence of dependency links for weak coupling strength, whereas the network connectivity mainly depends on the nodes with interdependency links in the system. Now, we will give an analysis of how to obtain the critical value . Figure 3 gives graphic solutions of and by using Eq. (8). From Fig. 3(a), we find that there exists a trivial solution at . As , forms a tangent point, at which starts to appear at a finite value. It means that can be obtained from the condition , as shown in Fig. 3(b). Furthermore, Fig. 3(c) shows that has two intersection points for . For this case, the bigger value is chosen to , since it denotes the size of the biggest component of the system. In order to explore the influence of coupling strength on system robustness, we show the critical value as a function of in Fig. 4. From Fig. 4, we notice as a function of with different parameters for . The results imply that the system becomes more robust as increases, as shown in Fig. 4. As the network density increases, the system becomes more resilient and easy to defend. This will prompt a better understanding and further realization in designing more resilient dependent systems.
and as functions of with parameters . Analytic predictions in comparison with simulation results. Symbols and lines represent simulation results and theoretical results. (a) and . (b) and (c) Different coupling strength with 4 and 8, respectively. Simulation results are averaged over 100 realizations with the size .
and as functions of with parameters . Analytic predictions in comparison with simulation results. Symbols and lines represent simulation results and theoretical results. (a) and . (b) and (c) Different coupling strength with 4 and 8, respectively. Simulation results are averaged over 100 realizations with the size .
Graphic solutions of and with parameters and for different [(a) , (b) , (c) ] from Eq. (8) and the simulation results are averaged over 100 realizations with size .
Graphic solutions of and with parameters and for different [(a) , (b) , (c) ] from Eq. (8) and the simulation results are averaged over 100 realizations with size .
Next, we will compare our theoretical results with the numerical results in partially interdependent networks with multiple-to-multiple relationship, whose degree distributions follow power law distribution 20 and dependency links follow Poisson degree distribution. Thus, Eqs. (1)–(5) become
By comparing analytical and simulation results, we find that the simulation results agree well with the analytical results, as shown in Fig. 5. Figure 5(a) shows that becomes smaller if the value of becomes bigger, which means that one can get a robust system by increasing the dependency density. From Fig. 5(b), the results imply that the system becomes more robust as the fraction of dependency nodes increases for our model. For different power law exponents , as a function of is shown in Fig. 6(a) for = = 4. We find that the value of is getting smaller with a smaller , which means that the interdependent system is more robust for smaller power law exponents. For fixed , values get smaller as the average dependency degree increases, and the system here becomes more steady and robust, as shown in Fig. 6(b).
Analytic predictions (line) in comparison with simulation results (symbol) with different parameters. (a) The parameters , with different . (b) The parameters , with different . Simulation results are averaged over 100 realizations with size .
Analytic predictions (line) in comparison with simulation results (symbol) with different parameters. (a) The parameters , with different . (b) The parameters , with different . Simulation results are averaged over 100 realizations with size .
as a function of . (a) For different with . (b) For different with .
V. CONCLUSIONS
In this paper, the robustness of two multiple-to-multiple dependent networks is numerically and analytically studied by extending previous works.28,31 By defining the mechanism of cascading failures of the system, we developed a framework and obtained the theoretical prediction for the size of the giant component at the stable state. The analytical results agree well with the simulation results for coupled networks. Additionally, the graphical solution of the critical threshold is numerically analyzed for the above models. The results suggest that the system behaves as a discontinuous phase transition for the above two systems. The robustness of a single network increases as the coupling strength increases. We expect that the increased network density can also strengthen the robustness of the network for both Erdös-Rnyi and scale-free networks. Moreover, the studies imply that one can obtain a more robust system by increasing the power exponent index of the system. The case of robustness of the general case of coupled networks with multiple-to-multiple dependent relationship has been studied very recently.
ACKNOWLEDGMENTS
We acknowledge the Major Program of National Natural Science Foundation of China (NNSFC) (Grant Nos. 71690242, 91546118, and 61803184), the Israel-Italian Collaborative Project NECST, the Israel Science Foundation, ONR, the Japan Science Foundation, BSFNSF, and DTRA (Grant No. HDTRA-1-10-1-0014) for financial support and Key Research Program of Frontier Sciences, CAS, Grant No. QYZDJ-SSWSYS019. F.W. acknowledges the fellowship program funded by the Planning and Budgeting Committee of the Council for Higher Education of Israel.