Skip to main content
  • Submit
  • About
    • Editorial Board
    • PNAS Staff
    • FAQ
    • Accessibility Statement
    • Rights and Permissions
    • Site Map
  • Contact
  • Journal Club
  • Subscribe
    • Subscription Rates
    • Subscriptions FAQ
    • Open Access
    • Recommend PNAS to Your Librarian
  • Log in
  • My Cart

Main menu

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses
  • Submit
  • About
    • Editorial Board
    • PNAS Staff
    • FAQ
    • Accessibility Statement
    • Rights and Permissions
    • Site Map
  • Contact
  • Journal Club
  • Subscribe
    • Subscription Rates
    • Subscriptions FAQ
    • Open Access
    • Recommend PNAS to Your Librarian

User menu

  • Log in
  • My Cart

Search

  • Advanced search
Home
Home

Advanced Search

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses

New Research In

Physical Sciences

Featured Portals

  • Physics
  • Chemistry
  • Sustainability Science

Articles by Topic

  • Applied Mathematics
  • Applied Physical Sciences
  • Astronomy
  • Computer Sciences
  • Earth, Atmospheric, and Planetary Sciences
  • Engineering
  • Environmental Sciences
  • Mathematics
  • Statistics

Social Sciences

Featured Portals

  • Anthropology
  • Sustainability Science

Articles by Topic

  • Economic Sciences
  • Environmental Sciences
  • Political Sciences
  • Psychological and Cognitive Sciences
  • Social Sciences

Biological Sciences

Featured Portals

  • Sustainability Science

Articles by Topic

  • Agricultural Sciences
  • Anthropology
  • Applied Biological Sciences
  • Biochemistry
  • Biophysics and Computational Biology
  • Cell Biology
  • Developmental Biology
  • Ecology
  • Environmental Sciences
  • Evolution
  • Genetics
  • Immunology and Inflammation
  • Medical Sciences
  • Microbiology
  • Neuroscience
  • Pharmacology
  • Physiology
  • Plant Biology
  • Population Biology
  • Psychological and Cognitive Sciences
  • Sustainability Science
  • Systems Biology
Research Article

System crash as dynamics of complex networks

Yi Yu, Gaoxi Xiao, Jie Zhou, Yubo Wang, Zhen Wang, Jürgen Kurths, and Hans Joachim Schellnhuber
PNAS October 18, 2016 113 (42) 11726-11731; first published October 3, 2016; https://doi.org/10.1073/pnas.1612094113
Yi Yu
aSchool of Electrical and Electronic Engineering, Nanyang Technological University, Singapore 639798;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Gaoxi Xiao
aSchool of Electrical and Electronic Engineering, Nanyang Technological University, Singapore 639798;
bComplexity Institute, Nanyang Technological University, Singapore 639798;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: EGXXiao@ntu.edu.sg zhenwang0@gmail.com john@pik-potsdam.de
Jie Zhou
cDepartment of Physics, East China Normal University, Shanghai 200241, China;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Yubo Wang
aSchool of Electrical and Electronic Engineering, Nanyang Technological University, Singapore 639798;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Zhen Wang
dQingdao University, Qingdao, Shandong 266071, China;
eInterdisciplinary Graduate School of Engineering Sciences, Kyushu University, Fukuoka 816-8580, Japan;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: EGXXiao@ntu.edu.sg zhenwang0@gmail.com john@pik-potsdam.de
Jürgen Kurths
fPotsdam Institute for Climate Impact Research, 14473 Potsdam, Germany;
gDepartment of Physics, Humboldt University, 12489 Berlin, Germany;
hInstitute for Complex Systems and Mathematical Biology, University of Aberdeen, Aberdeen AB24 3UE, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Hans Joachim Schellnhuber
fPotsdam Institute for Climate Impact Research, 14473 Potsdam, Germany;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: EGXXiao@ntu.edu.sg zhenwang0@gmail.com john@pik-potsdam.de
  1. Contributed by Hans Joachim Schellnhuber, August 11, 2016 (sent for review June 8, 2016; reviewed by Ying-Cheng Lai and Matjaz Perc)

  • Article
  • Figures & SI
  • Info & Metrics
  • PDF
Loading

Significance

System crash, as an essential part of system evolution, sometimes happens in peculiar manners: Weakened systems may survive for a surprisingly long time before suddenly meeting their final ends, whereas seemingly unbeatable giants may drastically crash to virtual nonexistence. We propose a model that describes system crash as a consequence of some relatively simple local information-based individual behaviors: Individuals leave networks according to some most straightforward assessment of current and future benefits/risks. Of note, such a simple rule may enable a single push/mistake to cause multistage-style system crash. Our study helps to make sense of the process where complex systems go into unstoppable cascading declines and provides a viewpoint of predicting the fate of some social/natural systems.

Abstract

Complex systems, from animal herds to human nations, sometimes crash drastically. Although the growth and evolution of systems have been extensively studied, our understanding of how systems crash is still limited. It remains rather puzzling why some systems, appearing to be doomed to fail, manage to survive for a long time whereas some other systems, which seem to be too big or too strong to fail, crash rapidly. In this contribution, we propose a network-based system dynamics model, where individual actions based on the local information accessible in their respective system structures may lead to the “peculiar” dynamics of system crash mentioned above. Extensive simulations are carried out on synthetic and real-life networks, which further reveal the interesting system evolution leading to the final crash. Applications and possible extensions of the proposed model are discussed.

  • complex systems
  • system crash
  • pseudo-steady state
  • cascade behavior

Systems emerge and grow up, and, naturally, systems die. Although extensive studies have been carried out on the growth and evolution of various complex systems, in many cases adopting complex network-based approaches (e.g., refs. 1 and 2), hardly any similar studies exist on system crash. As a special type of system evolution, however, system crash defies straightforward description by any of the existing models of system evolution, to the best of our knowledge. A few important observations of system crash remain puzzling and require careful study: (i) Some “inferior” or “outdated” systems, although appearing to be doomed to fail, may survive for a very long time. Examples include “living fossils” such as coelacanths, which have survived for more than 80 million years (3), some depleted-but-not-eliminated species under competition or faced by seemingly unstoppable invasion (4), and social systems such as the final Qing dynasty in China, which went through 70 y of disastrous defeats and rebellions with surprising robustness until its sudden termination in 1911 (5). We term such systems as being in a “pseudo-steady state” before their final crashes. (ii) Systems that appear to be too big or too strong to fail sometimes crash rapidly. Well-known examples include the crashes of various ecological and biological systems, which, in some cases, occur much faster than expected (6, 7) or even for no obvious reason (8, 9). Similar phenomena exist in human society, e.g., the sudden crash of the Soviet Union (10) and the crash of the once biggest online social network (OSN), Friendster, in less than 1 y (11). To help understand the crashes of ecological (6, 7, 12), biological (8, 9), neurological (13), physical and cyberphysical (14), economic (15, 16), social (5, 10, 11), and many other complex systems, new modeling approaches to describe the process and dynamics of a system crash are needed.

Some relevant studies, although not focusing directly on system crash, may help in understanding a few aspects of system crash or system crash under some special cases, and are hence worth mention. The largest class of relevant studies is probably those on herd behavior, which refers to system behavior where individuals in a group act collectively without centralized control (17⇓⇓⇓⇓–22). For example, conformity theories (17), which focus on exploring the mechanism leading to uniform behavior, may help to explain a sudden crash when (and only when) most individuals are close to the borderline between alternatives; informational cascade theories (18, 19) point out that, if individuals follow the behavior of the preceding individuals without regarding their own information, a big shift in the system may be caused by a small shock; and studies on lock-in effects and switching cost (20, 21), which evaluate effects of the cost that individuals have to pay to change a choice, may help explain the existence of a pseudo-steady state. Such studies, however, have their respective limitations: Conformity theories in themselves cannot explain how individuals get close to the borderline; informational cascade theories may not work well for system crashes where individuals do not easily disregard their own information; and the critical question of how individuals overcome the lock-in effects in the system crash remains largely open.

Although most of the above studies have largely ignored effects of the underlying structures of connections (22), an important trend in research on complex systems is to study them in the content of various complex networks. The most notable studies that combine system crash dynamics and network analysis are probably those on k-core cascade theory (11, 23). Specifically, the k-core of a network is defined as the maximum subset of the network where each node is connected to at least k other nodes in the same subset. Assume that a network node with fewer than k connections may have a chance to leave the current system [e.g., an OSN user may need to have a large enough number of connections to justify the effort of staying on (11)]. Measuring the k-core of a network therefore estimates how many nodes may stay on when the value of k increases due to various reasons. This theory, however, can only explain a sudden system crash when there exists a big jump in the value k due to some special reasons. Other interesting developments include (i) system crash model based on node energy level and internodal energy transfer (24) and (ii) network-based ranking of the extinction risks of different species in food webs or mutualistic networks (12).

Altogether, although various theories have been developed in different areas of science to make sense of abrupt changes in complex systems, what is still lacking is a system dynamics model that could properly describe a system crash, including the pseudo-steady state in some cases, as evolution of the system dynamics in their respective system structures. Note that the system crash discussed in this contribution is different from the cascading failure of a complex system or multiple interdependence systems discussed in the literature (25, 26), where components of the system(s) are strongly coupled and the sequence of a cascading failure is mainly decided by such coupling effects, e.g., failure of a node A may lead to failure of its strongly coupled neighboring node B but not of another weakly coupled neighboring node C.

In this contribution, we propose a complex network-based model for describing the process and dynamics of a system crash, and the pseudo-steady state in some cases as well. Specifically, the proposed model is as follows: Given a complex network, each network node may leave the network at a certain probability either (i) when the node has fewer than ks connections, in which case it may not have enough support or benefit to stay on, or (ii) when the node has lost more than a certain proportion (denoted by q) of its neighbors, in which case it may become a more attractive option to leave the current system, either to lower/avoid the risk (e.g., to avoid becoming a victim of a sudden system/herd crash) or to join another system with a more promising future. As we see, in this model, the first part reflects a value/risk assessment of the present situation, whereas the second part measures the effects of a relatively simple counting-based “copying” action, which may be a result of certain calculations and predictions of the risk, benefit, and/or future developments of the system. We term this model the “KQ-cascade model.” Note that the KQ-cascade model is substantially different from the k-core-based model in ref. 11, and it is also different from the threshold model in refs. 27 and 28, which generally assumes that nodes leave when they lose q proportion of neighbors, with a main focus on analyzing the threshold value of q in causing a global crash after a trivial or random proportion of nodes are initially removed. The KQ-cascade model includes the latter two models as special cases (hereafter denoted the k-cascade and q-cascade models, respectively), but it leads to much more complex dynamics. For example, although it is difficult to provide a rigid proof, in our extensive numerical simulations on various synthetic and real-life networks, neither the k-cascade model nor the q-cascade model has ever allowed the emergence of the pseudo-steady state.

The proposed model is tested on various synthetic and real-life networks. Considering the increasing importance of OSNs and the relative easiness in getting the data of OSN structures, the testing on real-life networks will mainly focus on OSNs, although the proposed model certainly applies to many different kinds of complex systems. For the well-known case of Friendster’s quick crash, where there exist relatively abundant data on the whole process (11), we perform individual-based simulations to mimic the actions of each individual. A good match is achieved with some interesting insights.

Results

The KQ-Cascade Model.

We start with an original network G(V,E), where V is the set of vertices and E is the set of edges. To keep the KQ-cascade model simple, we introduce only two key parameters: the critical degree ks and the loss-tolerance coefficient q. In the ith time step, any node with a degree k<ks or having lost more than q proportion of its original connections may leave the network with a probability fi(k).

An example of the KQ-cascade model is illustrated in Fig. 1. Note that the classic k-core cascading can be viewed as a special case of the KQ cascade where fi(k)=1 for k<ks and fi(k)=0 otherwise. As mentioned above, the value of ks quantifies the minimum support or benefit an individual needs to have to justify staying in the system: A less user-friendly OSN, for example, may lead to a higher ks value. The value of q, on the other hand, reveals the individuals’ risk tolerance level or the prospect on the future of the system they are staying in: A lower value of q reveals a lower risk tolerance or a less positive prospect; for those cases with competition between different systems, it may reveal a higher competition pressure from the competitor(s) as well. In real life, the value of the aggregated parameter q may be affected by many factors, e.g., risk tolerance level, switching cost, absolute/relative group size, and various environment factors.

Fig. 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 1.

Schematic illustration of KQ cascade with ks=3 and q=0.5. In a classic k-core cascade, the four nodes within the yellow square forming up a three-core would stay on while other nodes would leave. In the KQ cascade, however, because node 1 has lost more than 50% of its original neighbors, it will leave in the next time step, which leaves each of the three nodes within the blue triangle with fewer than three connections, leading to the final crash of the network.

Theoretical Analysis on Network Evolution.

Under the KQ cascade, a network may demonstrate a phase transition of cascade size, measured by the fraction of network nodes finally remaining in the network, under different values of ks and q: It may either crash into virtual nonexistence or have a nontrivial proportion of nodes remaining in the steady state.

A theoretical analysis on the network evolution under the KQ cascade, which allows an accurate prediction of such a phase transition in random networks by reproducing the whole procedure of the cascade, has been developed. The main idea of the analysis is presented in Methods, and further details can be found in SI Appendix, SI Appendix 1: Analytical Approach. To illustrate the accuracy of the analysis, we present analytical solutions versus simulation results averaged over 100 independent realizations (details for implementing these independent realizations are given in Methods) for the following paradigmatic random networks: an Erdös−Rényi (ER) network (29) with an average nodal degree z=20; an exponential (Exp) network with an average degree z=20 and a degree cutoff of 100; and a scale-free (SF) network with γ=2, a minimum degree of 3, and a degree cutoff of 100 (30). All three types of networks are generated using the configuration model (30) with size of N=104. For simplification, we assume that, in each time step, individuals fulfilling the conditions to leave may decide to do so with a constant probability f, which may be viewed as a special case of the leaving probability fi(k). As we observe in Fig. 2, the proposed theoretical analysis accurately predicts the evolution of the network size in all these cases. The accurate prediction of the network dynamics also allows us to calculate the threshold value of q leading to the network crash, denoted as qth, by adopting a simple trial-and-error approach, as illustrated in Fig. 3. The thresholds of some real-life networks, presented in SI Appendix, Fig. S1, and related discussions can be found in SI Appendix, SI Appendix 2: Cascade Threshold of Real-Life Networks. Further results of the proposed analysis in predicting the cascade process are given in SI Appendix, Figs. S2–S7. Simulation results showing how the degree distribution affects the resilience and cascade size (i.e., fraction of the remaining nodes) of random networks are presented in SI Appendix, Figs. S14–S20. Analytical and simulation results for the case where decline of a system shakes the individuals’ confidence, leading to an accelerated system crash (31), are reported in SI Appendix, SI Appendix 5: The Speedup Loss of Individuals in a System Crash.

Fig. 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 2.

KQ cascade in various networks. Comparison between analytical (lines) and simulation (symbols) results for KQ cascade in different networks with size N=104.

Fig. 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 3.

Threshold of KQ cascade. Comparison between analytical (lines) and simulation (symbols) results of cascade thresholds in different networks with size N=104. Both analytical and simulation results are obtained by adopting a trial-and-error approach, where, for each given ks, the value of q is increased by a step length of 0.01 until the threshold value is obtained.

Pseudo-Steady State and Sudden Crash.

We find that the proposed model enables the occurrence of the pseudo-steady state and a sudden crash of the systems. A few such cases in both random and real-life networks are illustrated in Fig. 4. As we can see, in the pseudo-steady state, the networks appear to be rather stable, with only a few nodes leaving in each time step. After a long period, however, the systems suddenly crash, sometimes within only a few steps. In some systems, e.g., LiveJournal, the pseudo-steady state may even appear more than once (more details will be discussed later). Note that, when q gets close to qth, the specific time when the crash starts can be rather sensitive to small fluctuations in the network structure. An example is given in Fig. 4A, where we show results in 10 random networks with the same parameters but different seeds for random number generation. Therefore, in this section, we shall only present our results for a single simulation case.

Fig. 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 4.

Pseudo-steady state and sudden crash. (A) Comparisons between theoretical and simulation results of 10 independent realizations in scale-free networks with the same parameters yet different seeds for random number generation: ks=4, q=0.29, f=0.2. (B) Simulation results showing the pseudo-steady state and sudden crash in different systems: Orkut (stars), Scale-free (triangles), LiveJournal (diamonds), Exponential (circles), and Erdős–Rényi (ER) (squares).

We reveal that the pseudo-steady state and sudden crash of the networks are caused by the emergence, growth, and final crash of a giant cluster composed of “vulnerable nodes,” where a vulnerable node is defined as the node that will fulfill the conditions of leaving the network if it loses one more neighbor. It is known that the vulnerable clusters composed by connected vulnerable nodes play an important role in a global cascade (27), as the departure/loss of a single node in such a cluster may trigger the cascading departure/loss of the whole cluster, which may further result in the crash of the entire system. It is, however, shown that a basic rule for leaving can lead to the emergence of a giant vulnerable cluster in complex networks.

Fig. 5 illustrates, in more detail, the growth and decline of vulnerable clusters. The simulations are carried out on a random exponential network, as defined in Theoretical Analysis on Network Evolution. To make comparisons, we choose two sets of parameters: ks = 14, q = 0.37, and f = 0.1 and ks = 14, q = 0.38, and f = 0.1. Although q=0.37 leads to a network crash, q=0.38 allows the network to survive; choosing such parameters thus enables us to closely observe a phase transition of the system. We find that, for the case where the network finally crashes, the number of vulnerable nodes slowly accumulates during the pseudo-steady state (Fig. 5B). The average relative size of all of the vulnerable clusters meanwhile remains between 10−4 and 2×10−4, meaning having only one or two nodes (Fig. 5D). This finding shows that most vulnerable clusters are tiny pieces scattered within the network. The size of the largest vulnerable cluster also remains rather small most of the time (Fig. 5C). Shortly before the sudden crash starts, however, there is a sharp increase in the size of the largest vulnerable cluster (Fig. 5C), when vulnerable clusters quickly connect together (Fig. 5D). At the moment when a sudden crash starts, almost all of the vulnerable nodes merge into a single cluster. For example, in the 69th time step, among 683 vulnerable nodes, 632 of them merge into the largest vulnerable cluster. The sudden emergence of the giant vulnerable cluster prepares a sufficient condition for the sudden crash to be easily triggered.

Fig. 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 5.

Evolution of the vulnerable clusters. Evolution of (A) network size, (B) overall size of vulnerable nodes, (C) the size of the largest vulnerable cluster, and (D) the average size of all of the vulnerable clusters in the random exponential network.

Measuring the Resilience of Some Real-Life Systems Against the KQ Cascade.

It is interesting to evaluate the resilience of a few real-life systems against the KQ cascade. In this section, we report numerical simulation results on three different OSNs, namely, the Orkut, LiveJournal, and YouTube networks (32). A summary of the basic information on these networks can be found in Methods.

Fig. 6 shows the dynamics of the cascade size in these networks with different values of ks and q. It is interesting to observe that Orkut demonstrates the strongest resilience against the KQ cascade among the three networks, whereas YouTube turns out to be the most fragile one. The conclusions, however, have to be taken with a pinch of salt because, as always, the data we have only reflect a fraction of the corresponding real-life networks, and it is not known how the network sampling was done in the first place. Some discussion of how much evaluating the resilience of a random sampling of a complex network may help to reveal the resilience of the whole system is presented in SI Appendix, SI Appendix 6: Can Random Sampling of a Network Reflect the Resilience.

Fig. 6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 6.

Cascade size of a few real-life networks. Cascade size of (A) Orkut online social network, (B) LiveJournal online social network, and (C) YouTube network. The cascade size is shown in color scale.

Another interesting observation is that, in real-life networks, there may be multistage phase transitions of the cascade size: The cascade size may go through multiple pseudo-steady states before the final crash. This, however, has never been observed in our extensive simulations on uncorrelated random networks. We believe that such phenomena are related to community structures (33) and degree correlations (34) existing in real-life networks. For example, the cascading departure/loss of a large number of individuals may have barely any impact on certain communities with dense intracommunity connections. It may be worth mentioning that such observations have been made in many OSNs: Friendster’s popularity was not significantly affected in Southeast Asia, especially the Philippines, throughout its fast decline; Orkut was especially popular in Brazil; and LiveJournal has 52% of customer visits from Russia (35). A related result in ref. 36 reports that, in a loosely coupled two-community network, system cascade may have one peak in each community, separated in time. Our simulation study shows that, in networks with fixed nodal degrees, multistage pseudo-steady state vanishes with the elimination of the community structures and degree correlations. Details are reported in SI Appendix, SI Appendix 7: Effects of Community Structures and Degree Correlations.

Cascade Decline of Friendster: A Possible Explanation.

The crash of Friendster offers an interesting case with relatively abundant data. In a recent study (11), it was found that, if we adopt the simple k-cascade model and let the threshold value of k increase continuously from 3 in July 2009 to 67 in June 2010 at a rate of about 6 per month, a good match between simulation results and historical records of Friendster’s crash could be achieved. A puzzling question remains, however: Although it is known that Friendster made a fatal mistake in 2009 when it changed its website interface, making it more difficult to use and hence increasing the threshold value ks, it is not clear how this threshold value could have gone up continuously when Friendster did not make a second mistake.

We apply the proposed cascade model to the Friendster network, which consists of 65,608,366 nodes and 1,806,067,135 connections (11). As obtaining the exact number of Friendster users over time is difficult, following the work in ref. 11, we use the Google search volume to approximate the popularity evolution of Friendster. Specifically, the curve is still figured by obtaining the search volume of “www.friendster.com.” Two reference points are set, one in June 2009, when Friendster began to decline as users were not happy with the changed interface (probably also due to the fast growth of Facebook) (11), and the other in July 2010, when Friendster was reported to have only about 10 million active users left, less than 10% of its peak size. A trial-and-error approach shows that, when we set ks=20, q=0.2, and f=0.15, simulation results based on the proposed model match well with the real-life data (Fig. 7). Note that, as mentioned above, for those cases where systems crash mainly due to strong competition, a smaller value of q may imply a higher competition pressure from a stronger competitor. Although q=0.2 may appear to be quite a low value, adopting such a value is not without basis: Even when Friendster was still at its peak user size in 2009, the Google search volume of “Facebook” was already more than 20 times higher than that of “Friendster” (37), as illustrated in Fig. 7A. People thus may have had good reason to believe, back in 2009–2010, that although Friendster was larger, Facebook would certainly boom (and such belief helped Facebook actually to boom). Also note that, at that time, many users may have registered their accounts in both Friendster and Facebook; leaving Friendster at q=0.2 therefore did not necessarily mean that they had to lose 80% of their online social connections; instead, it might only have meant getting rid of some inconveniences once and for all.

Fig. 7.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 7.

Cascade decline of the Friendster. (A) Comparison between the Google search volumes of “Friendster” and “Facebook.” In January 2009, the search index of “Facebook” was 27, whereas, for “Friendster,” it was 1. (B) The monthly Google search volume for “www.friendster.com” (triangles) and simulation results adopting the proposed model (circles). The size of 1 corresponds to the user size at the peak time (excluding those users with zero degree), which was about 68 millions. Decline started at a size of 0.94 (i.e., 64 million) in June 2009 and stopped at a size of 0.15 (about 10 million) in June 2010.

With only a snapshot showing the aggregated connections of all of the ever-existing users until the moment the snapshot was taken, it is not a surprise that we have to adopt a trial-and-error approach to estimate the ks and q values. Nevertheless, it is encouraging to see that, without making the assumption that the ks value increases over time, a good match between simulation results and real-life data can be achieved. When detailed data showing network topology in different stages of system decline are available, our model may allow a more accurate estimation of parameter values and, consequently, a more accurate reflection/prediction of system dynamics.

Discussion

In this work, we have introduced a network-based system dynamics model for describing the crash of complex systems and the pseudo-steady state in some cases. In the proposed model, a network node may choose to leave either when the number of connections it has becomes too low or when it has lost more than a certain proportion of its neighbors. We have derived a theoretical analysis for the crash process on a random graph with an arbitrary nodal degree distribution and an arbitrary leaving probability. Based on the proposed model, a pseudo-steady state and sudden crash phenomenon could be steadily observed in certain ranges of parameters and be easily explained. Further, the resilience of some real-life networks has been evaluated, and a possible explanation for the sudden crash of Friendster has been presented.

The proposed model may find wide applications in helping understand and predict the declines of various complex systems, especially complex social systems. Studies on such applications would be of future research interest, in particular:

  • i) Research areas may heat up and cool down. Although “early movers” may leave a research area when important work has been done or low-hanging fruit has been collected, many others may only make up their minds to leave when their colleagues are leaving (similar to q cascade) or when they have lost their collaborators (similar to k cascade). Dynamics of the decline of a research area needs to be studied.

  • ii) Collective intelligence systems such as Wikipedia may have an increasing coordination cost when growing in scope (38). Whether and how participants of a collective intelligence system may decide to leave due to increasing coordination cost are surely worth careful studies.

  • iii) Decentralized adoptions of new technologies, such as voluntary installations of solar panels on house roofs, may be subject to certain constraints. For example, there may be an upper bound to the penetration level of grid-tied photovoltaic power (39). Certain policies therefore may have to be installed to restrict the adoption of the technology. Predicting whether adopting such policies might result in decline or even crash of the technology adoption is of significant importance.

It would also be of future research interest to figure out how much a model of system crash dynamics might help to identify tipping points of a system (40), and to model system dynamics when individuals’ decisions to leave involve certain kinds of more sophisticated Bayesian reasoning (41).

Methods

Theoretical Analysis: Main Idea.

The main idea of the theoretical analysis is to model the degree transition in every step as a Markov process: By taking each pair of original−current nodal degrees as a state and calculating the state transition in each time step, we can reproduce the evolution process of the network. Specifically, we construct (i) a degree transition matrix Di where the element Djki denotes the probability that a node with an original degree j becomes a degree-k node at the beginning of the ith time step and (ii) a matrix Ui where the element Ujki denotes the probability that a node with an original degree j has a degree k at the beginning of the ith time step and it does not leave the network in this time step. We haveUjki=Djkiμ(k,j),[1]

where μ(k,j) reflects the probability that a node with an original degree j and a degree k at the beginning of the ith time step does not leave the system in this time step,μ(k,j)={1k≥(1−q)j and k≥ks,1−fi(k)otherwise,[2]

where, as mentioned in Results, fi(k) denotes the probability that a node fulfilling the conditions to leave may actually leave the network in the ith time step. As some nodes will lose a proportion of their neighbors in the leaving process, their degrees need to be recalculated. Use the matrix Ti to keep record of the transition within this time step, where Tk′ki denotes the probability that a node with a degree k′ at the beginning of the ith time step ends up with a degree k at the end of this step. The degree transition matrix Di+1 hence can be calculated asDi+1=UiTi.[3]

The whole KQ-cascade process can be reproduced by iterative calculation of Eqs. 1 and 3. The detailed calculations of these matrixes are discussed in SI Appendix, SI Appendix 1: Analytical Approach.

This analysis can be easily extended to a more general case where network nodes leave the network following an arbitrary criterion ϕ as long as the new criterion can be reflected by properly changing μ(k,j) in Eq. 2 accordingly. Such extensions would allow an easy coverage of a wide range of existing work, e.g., the classic k-core problem.

Numerical Simulations.

For an easy reference, Table 1 summarizes the real-life networks adopted in the numerical simulations. Note that, although the data of Friendster contain 117 million identifications (IDs), only 65 million of them have a connection record, as reflected in Table 1; the other 52 million IDs are for private users with confidential connection information.

View this table:
  • View inline
  • View popup
Table 1.

Basic information of the networks evaluated in this work

In both the synthetic and real-life networks, numerical simulations are carried out in the following steps:

  • i) Initialization: record the degree of each node as its original degree. Set values of the parameters ks and q.

  • ii) First time step: for any node with its original degree lower than ks, remove it at a constant probability f. Update the record of each node’s degree as its current degree.

  • iii) Iterations in the following time steps: check all nodes’ current degrees as recorded at the beginning of the time step. If a node has a current degree lower than ks or has lost more than q proportion of its original degree, it has a probability f to leave the network. The records of the nodes’ current degrees are only updated when all of the nodes have decided whether to stay on or to leave, at the end of each time step. Repeat the iterations until the network diminishes or no more nodes can leave the network.

For each synthetic network with f<1, we carry out 100 independent realizations: First we generate 10 random networks and then carry out 10 rounds of simulations on each of these networks. Unless otherwise specified, we present the average of these 100 independent realizations. When f=1, in each of the 10 random networks, only a single round of simulation is needed. For the real-life networks, when f<1, we also carry out 10 rounds of simulation and present the average results. Note that, in all of the figures, the error bars are rather small; hence they are omitted.

Also note that different values of f only speed up or slow down the network evolution process; they have no influence on the crash threshold and the cascade size at the steady state. For the synthetic networks, we typically adopt small values such as f=0.1 or f=0.2 to generate smoother curves. For real-life networks, we let f=0.5 to speed up the computation.

Acknowledgments

This work is partially supported by the Ministry of Education (MOE), Singapore, under Research Grants RG 69/12, RG 28/14, MOE2013-T2-2-006, and MOE2014-T2-1-028. Part of this work is an outcome of the Future Resilient Systems project at the Singapore-ETH (Swiss Federal Institute of Technology in Zurich) Centre, which is funded by the National Research Foundation of Singapore under its Campus for Research Excellence and Technological Enterprise program.

Footnotes

  • ↵1Y.Y., G.X., and Z.W. contributed equally to this work.

  • ↵2To whom correspondence may be addressed. Email: EGXXiao{at}ntu.edu.sg, zhenwang0{at}gmail.com, or john{at}pik-potsdam.de.
  • Author contributions: Y.Y., G.X., J.Z., Y.W., Z.W., J.K., and H.J.S. designed research; Y.Y., G.X., J.Z., Y.W., and Z.W. performed research; and Y.Y., G.X., Z.W., J.K., and H.J.S. wrote the paper.

  • Reviewers: Y.-C.L., Arizona State University; and M.P., University of Maribor.

  • The authors declare no conflict of interest.

  • This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1612094113/-/DCSupplemental.

View Abstract

References

  1. ↵
    1. Boccaletti S,
    2. Latora V,
    3. Moreno Y,
    4. Chavez M,
    5. Hwang D-U
    (2006) Complex networks: Structure and dynamics. Phys Rep 424(4-5):175–308.
    .
    OpenUrlCrossRef
  2. ↵
    1. Song C,
    2. Havlin S,
    3. Makse HA
    (2006) Origins of fractality in the growth of complex networks. Nat Phys 2:275–281.
    .
    OpenUrlCrossRef
  3. ↵
    1. Forey PL
    (1998) History of the Coelacanth Fishes (Chapman & Hall, London).
    .
  4. ↵
    1. Zhi Y, et al.
    (2007) Inter-specific competition, Spartina alterniflora is replacing spartina anglica in costal China. Estuar Coast Shelf Sci 74(3):437–448.
    .
    OpenUrlCrossRef
  5. ↵
    1. Fairbank JK,
    2. Liu KC
    (1980) Late Ch’ing, 1800−1911, Part 2, The Cambridge History of China (Cambridge Univ Press, Cambridge, UK), Vol 11.
    .
  6. ↵
    1. Sahney S,
    2. Benton MJ,
    3. Falcon-Lang HJ
    (2010) Rainforest collapse triggered Pennsylvanian tetrapod diversification in Euramerica. Geology 38(12):1079–1082.
    .
    OpenUrlAbstract/FREE Full Text
  7. ↵
    1. Angus I,
    2. Butler S
    (2011) Too Many People? Population, Immigration, and the Environmental Crisis (Haymarket, Chicago).
    .
  8. ↵
    1. Houlahan JE,
    2. Findlay CS,
    3. Schmidt BR,
    4. Meyer AH,
    5. Kuzmin SL
    (2000) Quantitative evidence for global amphibian population declines. Nature 404(6779):752–755.
    .
    OpenUrlCrossRefPubMed
  9. ↵
    1. Stuart SN, et al.
    (2004) Status and trends of amphibian declines and extinctions worldwide. Science 306(5702):1783–1786.
    .
    OpenUrlAbstract/FREE Full Text
  10. ↵
    1. Gaddis JL
    (2006) The Cold War: A New History (Penguin, New York).
    .
  11. ↵
    1. Garcia D,
    2. Mavrodiev P,
    3. Schweitzer F
    (2013) Social resilience in online communities: The autopsy of Friendster. Proceedings of the First ACM Conference on Online Social Networks (Assoc Comput Machinery, Boston), pp 39–50.
    .
  12. ↵
    1. Domínguez-García V,
    2. Muñoz MA
    (2015) Ranking species in mutualistic networks. Sci Rep 5:8182.
    .
    OpenUrlCrossRefPubMed
  13. ↵
    1. Wenning GK,
    2. Colosimo C,
    3. Geser F,
    4. Poewe W
    (2004) Multiple system atrophy. Lancet Neurol 3(2):93–103.
    .
    OpenUrlCrossRefPubMed
  14. ↵
    1. Krishna PV,
    2. Saritha V,
    3. Sultana HP
    (2014) Challenges, Opportunities, and Dimensions of Cyber-Physical Systems (IGI Global, Hershey, PA).
    .
  15. ↵
    1. Sieczka P,
    2. Sornette D,
    3. Holyst JA
    (2011) The Lehman Brothers effect and bankruptcy cascades. Eur Phys J B 82:257–259.
    .
    OpenUrlCrossRef
  16. ↵
    1. Podobnik B, et al.
    (2015) The cost of attack in competing networks. J R Soc Interface 12(112):20150770.
    .
    OpenUrlAbstract/FREE Full Text
  17. ↵
    1. Bernheim BD
    (1994) A theory of conformity. J Polit Econ 102(5):841–877.
    .
    OpenUrlCrossRef
  18. ↵
    1. Bikhchandani S,
    2. Hirshleifer D,
    3. Welch I
    (1992) A theory of fads, fashion, custom and cultural change as informational cascades. J Polit Econ 100(5):992–1026.
    .
    OpenUrlCrossRef
  19. ↵
    1. Chen Y-F
    (2008) Herd behavior in purchasing books online. Comput Human Behav 24(5):1977–1992.
    .
    OpenUrlCrossRef
  20. ↵
    1. Wilson CM,
    2. Price CW
    (2006) Do consumers switch to the best supplier? Oxf Econ Pap 62(4):647–668.
    .
    OpenUrl
  21. ↵
    1. Villas-Boas JM
    (2015) A short survey on switching costs and dynamic competition. Int J Res Mark 32(2):219–222.
    .
    OpenUrlCrossRef
  22. ↵
    1. Raafat RM,
    2. Chater N,
    3. Frith C
    (2009) Herding in humans. Trends Cogn Sci 13(10):420–428.
    .
    OpenUrlCrossRefPubMed
  23. ↵
    1. Dorogovtsev SN,
    2. Goltsev AV,
    3. Mendes JFF
    (2006) k-Core organization of complex networks. Phys Rev Lett 96(4):040601.
    .
    OpenUrlCrossRefPubMed
  24. ↵
    1. Dasgupta K, et al.
    (2008) Social ties and their relevance to churn in mobile telecom networks. Proceedings of the 11th International Conference on Extending Database Technology: Advances in Database Technology (Assoc Comput Machinery, Boston), pp 668–677.
    .
  25. ↵
    1. Buldyrev SV,
    2. Parshani R,
    3. Paul G,
    4. Stanley HE,
    5. Havlin S
    (2010) Catastrophic cascade of failures in interdependent networks. Nature 464(7291):1025–1028.
    .
    OpenUrlCrossRefPubMed
  26. ↵
    1. Saavedra S,
    2. Stouffer DB,
    3. Uzzi B,
    4. Bascompte J
    (2011) Strong contributors to network persistence are the most vulnerable to extinction. Nature 478(7368):233–235.
    .
    OpenUrlCrossRefPubMed
  27. ↵
    1. Watts DJ
    (2002) A simple model of global cascades on random networks. Proc Natl Acad Sci USA 99(9):5766–5771.
    .
    OpenUrlAbstract/FREE Full Text
  28. ↵
    1. Hackett A,
    2. Melnik S,
    3. Gleeson JP
    (2011) Cascades on a class of clustered random networks. Phys Rev E Stat Nonlin Soft Matter Phys 83(5 Pt 2):056107.
    .
    OpenUrlCrossRefPubMed
  29. ↵
    1. Erdös P,
    2. Rényi A
    (1959) On random graphs, I. Publ Math (Debrecen) 6:290–297.
    .
    OpenUrl
  30. ↵
    1. Catanzaro M,
    2. Boguñá M,
    3. Pastor-Satorras R
    (2005) Generation of uncorrelated random scale-free networks. Phys Rev E Stat Nonlin Soft Matter Phys 71(2 Pt 2):027103.
    .
    OpenUrlCrossRefPubMed
  31. ↵
    1. Spanjers W
    (2008) Loss of confidence and currency crises. Int J Econ Res 5:219–237.
    .
    OpenUrl
  32. ↵
    1. Yang J,
    2. Leskovec J
    (2012) Defining and evaluating network communities based on ground-truth. Proceedings of the ACM SIGKDD Workshop on Mining Data Semantics (Assoc Comput Machinery, Boston), pp 1–8.
    .
  33. ↵
    1. Newman MEJ
    (2006) Modularity and community structure in networks. Proc Natl Acad Sci USA 103(23):8577–8582.
    .
    OpenUrlAbstract/FREE Full Text
  34. ↵
    1. Newman MEJ
    (2002) Assortative mixing in networks. Phys Rev Lett 89(20):208701.
    .
    OpenUrlCrossRefPubMed
  35. ↵
    Alexa (2014) Traffic Rank for Livejournal.com. Available at www.alexa.com/siteinfo/livejournal.com. Accessed November 30, 2014.
    .
  36. ↵
    1. Galstyan A,
    2. Cohen P
    (2007) Cascading dynamics in modular networks. Phys Rev E Stat Nonlin Soft Matter Phys 75(3 Pt 2):036109.
    .
    OpenUrlCrossRefPubMed
  37. ↵
    Google (2014) Google Trend: Search terms Friendster vs. Facebook. Available at https://www.google.com/trends/explore?date=all&q=friendster,facebook. Accessed November 30, 2014.
    .
  38. ↵
    1. Kittur A,
    2. Suh B,
    3. Pendleton A,
    4. Chi EH
    (2007) He says, she says: Conflict and coordination in Wikipedia. Proceedings of the SIGCHI Conference on Human Factors in Computing Systems (Assoc Comput Machinery, Boston), pp 453–462.
    .
  39. ↵
    1. Eltawil MA,
    2. Zhao Z
    (2010) Gird-connected photovoltaic power systems: Technical and potential problems—A review. Renew Sustain Energy Rev 14(1):112–129.
    .
    OpenUrlCrossRef
  40. ↵
    1. Dakos V,
    2. Bascompte J
    (2014) Critical slowing down as early warning for the onset of collapse in mutualistic communities. Proc Natl Acad Sci USA 111(49):17546–17551.
    .
    OpenUrlAbstract/FREE Full Text
  41. ↵
    1. Baltag A,
    2. Christoff Z,
    3. Hansen JU,
    4. Smets S
    (2013) Logical models of informational cascades. Working paper (Vrije Universiteit Brussel, Brussels). Available at www.vub.ac.be/CLWF/SS/CASCADES-COLLEGE.pdf.
    .
    1. Cho E,
    2. Myers SA,
    3. Leskovec J
    (2011) Friendship and mobility: User movement in location-based social networks. Proceedings of the 17th ACM SIGKDD International Conference on Knowledge Discovery and Data Mining (Assoc Comput Machinery, Boston), pp 1082–1090.
    .
PreviousNext
Back to top
Article Alerts
Email Article

Thank you for your interest in spreading the word on PNAS.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
System crash as dynamics of complex networks
(Your Name) has sent you a message from PNAS
(Your Name) thought you would like to see the PNAS web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
System crash as dynamics of complex networks
Yi Yu, Gaoxi Xiao, Jie Zhou, Yubo Wang, Zhen Wang, Jürgen Kurths, Hans Joachim Schellnhuber
Proceedings of the National Academy of Sciences Oct 2016, 113 (42) 11726-11731; DOI: 10.1073/pnas.1612094113

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Request Permissions
Share
System crash as dynamics of complex networks
Yi Yu, Gaoxi Xiao, Jie Zhou, Yubo Wang, Zhen Wang, Jürgen Kurths, Hans Joachim Schellnhuber
Proceedings of the National Academy of Sciences Oct 2016, 113 (42) 11726-11731; DOI: 10.1073/pnas.1612094113
Digg logo Reddit logo Twitter logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Mendeley logo Mendeley
Proceedings of the National Academy of Sciences: 113 (42)
Table of Contents

Submit

Sign up for Article Alerts

Article Classifications

  • Physical Sciences
  • Applied Mathematics
  • Social Sciences
  • Social Sciences

Jump to section

  • Article
    • Abstract
    • Results
    • Discussion
    • Methods
    • Acknowledgments
    • Footnotes
    • References
  • Figures & SI
  • Info & Metrics
  • PDF

You May Also be Interested in

Abstract depiction of a guitar and musical note
Science & Culture: At the nexus of music and medicine, some see disease treatments
Although the evidence is still limited, a growing body of research suggests music may have beneficial effects for diseases such as Parkinson’s.
Image credit: Shutterstock/agsandrew.
Scientist looking at an electronic tablet
Opinion: Standardizing gene product nomenclature—a call to action
Biomedical communities and journals need to standardize nomenclature of gene products to enhance accuracy in scientific and public communication.
Image credit: Shutterstock/greenbutterfly.
One red and one yellow modeled protein structures
Journal Club: Study reveals evolutionary origins of fold-switching protein
Shapeshifting designs could have wide-ranging pharmaceutical and biomedical applications in coming years.
Image credit: Acacia Dishman/Medical College of Wisconsin.
White and blue bird
Hazards of ozone pollution to birds
Amanda Rodewald, Ivan Rudik, and Catherine Kling talk about the hazards of ozone pollution to birds.
Listen
Past PodcastsSubscribe
Goats standing in a pin
Transplantation of sperm-producing stem cells
CRISPR-Cas9 gene editing can improve the effectiveness of spermatogonial stem cell transplantation in mice and livestock, a study finds.
Image credit: Jon M. Oatley.

Similar Articles

Site Logo
Powered by HighWire
  • Submit Manuscript
  • Twitter
  • Facebook
  • RSS Feeds
  • Email Alerts

Articles

  • Current Issue
  • Latest Articles
  • Archive

PNAS Portals

  • Anthropology
  • Chemistry
  • Classics
  • Front Matter
  • Physics
  • Sustainability Science
  • Teaching Resources

Information

  • Authors
  • Editorial Board
  • Reviewers
  • Librarians
  • Press
  • Site Map
  • PNAS Updates

Feedback    Privacy/Legal

Copyright © 2021 National Academy of Sciences. Online ISSN 1091-6490