Higher-Order Components Dictate Higher-Order Contagion Dynamics in Hypergraphs.

The presence of the giant component is a necessary condition for the emergence of collective behavior in complex networked systems. Unlike networks, hypergraphs have an important native feature that components of hypergraphs might be of higher order, which could be defined in terms of the number of common nodes shared between hyperedges. Although the extensive higher-order component (HOC) could be witnessed ubiquitously in real-world hypergraphs, the role of the giant HOC in collective behavior on hypergraphs has yet to be elucidated. In this Letter, we demonstrate that the presence of the giant HOC fundamentally alters the outbreak patterns of higher-order contagion dynamics on real-world hypergraphs. Most crucially, the giant HOC is required for the higher-order contagion to invade globally from a single seed. We confirm it by using synthetic random hypergraphs containing adjustable and analytically calculable giant HOC.

K-selective percolation: A simple model leading to a rich repertoire of phase transitions.

We propose a K-selective percolation process as a model for iterative removals of nodes with a specific intermediate degree in complex networks. In the model, a random node with degree K is deactivated one by one until no more nodes with degree K remain. The non-monotonic response of the giant component size on various synthetic and real-world networks implies a conclusion that a network can be more robust against such a selective attack by removing further edges. From a theoretical perspective, the K-selective percolation process exhibits a rich repertoire of phase transitions, including double transitions of hybrid and continuous, as well as reentrant transitions. Notably, we observe a tricritical-like point on Erdos-Rényi networks. We also examine a discontinuous transition with unusual order parameter fluctuation and distribution on simple cubic lattices, which does not appear in other percolation models with cascade processes. Finally, we perform finite-size scaling analysis to obtain critical exponents on various transition points, including those exotic ones.

No-exclaves percolation.

Network robustness has been a pivotal issue in the study of system failure in network science since its inception. To shed light on this subject, we introduce and study a new percolation process based on a new cluster called an 'exclave' cluster. The entities comprising exclave clusters in a network are the sets of connected unfailed nodes that are completely surrounded by the failed (i.e., nonfunctional) nodes. The exclave clusters are thus detached from other unfailed parts of the network, thereby becoming effectively nonfunctional. This process defines a new class of clusters of nonfunctional nodes. We call it the no-exclave percolation cluster (NExP cluster), formed by the connected union of failed clusters and the exclave clusters they enclose. Here we showcase the effect of NExP cluster, suggesting a wide and disruptive collapse in two empirical infrastructure networks. We also study on two-dimensional Euclidean lattice to analyze the phase transition behavior using finite-size scaling. The NExP model considering the collective failure clusters uncovers new aspects of network collapse as a percolation process, such as quantitative change of transition point and qualitative change of transition type. Our study discloses hidden indirect damage added to the damage directly from attacks, and thus suggests a new useful way for finding nonfunctioning areas in complex systems under external perturbations as well as internal partial closures.

Critical behaviors of high-degree adaptive and collective-influence percolation.

How the giant component of a network disappears under attacking nodes or links addresses a key aspect of network robustness, which can be framed into percolation problems. Various strategies to select the node to be deactivated have been studied in the literature, for instance, a simple random failure or high-degree adaptive (HDA) percolation. Recently, a new attack strategy based on a quantity called collective-influence (CI) has been proposed from the perspective of optimal percolation. By successively deactivating the node having the largest CI-centrality value, it was shown to be able to dismantle a network more quickly and abruptly than many of the existing methods. In this paper, we focus on the critical behaviors of the percolation processes following degree-based attack and CI-based attack on random networks. Through extensive Monte Carlo simulations assisted by numerical solutions, we estimate various critical exponents of the HDA percolation and those of the CI percolations. Our results show that these attack-type percolation processes, despite displaying apparently more abrupt collapse, nevertheless exhibit standard mean-field critical behaviors at the percolation transition point. We further discover an extensive degeneracy in top-centrality nodes in both processes, which may provide a hint for understanding the observed results.

Coevolution and correlated multiplexity in multiplex networks.

Distinct channels of interaction in a complex networked system define network layers, which coexist and cooperate for the system's function. Towards understanding such multiplex systems, we propose a modeling framework based on coevolution of network layers, with a class of minimalistic growing network models as working examples. We examine how the entangled growth of coevolving layers can shape the network structure and show analytically and numerically that the coevolution can induce strong degree correlations across layers, as well as modulate degree distributions. We further show that such a coevolution-induced correlated multiplexity can alter the system's response to the dynamical process, exemplified by the suppressed susceptibility to a social cascade process.

Contagion dynamics on hypergraphs with nested hyperedges.

In complex social systems encoded as hypergraphs, higher-order (i.e., group) interactions taking place among more than two individuals are represented by hyperedges. One of the higher-order correlation structures native to hypergraphs is the nestedness: Some hyperedges can be entirely contained (that is, nested) within another larger hyperedge, which itself can also be nested further in a hierarchical manner. Yet the effect of such hierarchical structure of hyperedges on the dynamics has remained unexplored. In this context, here we propose a random nested-hypergraph model with a tunable level of nestedness and investigate the effects of nestedness on a higher-order susceptible-infected-susceptible process. By developing an analytic framework called the facet approximation, we obtain the steady-state fraction of infected nodes on the random nested-hypergraph model more accurately than existing methods. Our results show that the hyperedge-nestedness affects the phase diagram significantly. Monte Carlo simulations support the analytical results.

Waiting time dynamics of priority-queue networks.

We study the dynamics of priority-queue networks, generalizations of the binary interacting priority-queue model introduced by Oliveira and Vazquez [Physica A 388, 187 (2009)]. We found that the original AND-type protocol for interacting tasks is not scalable for the queue networks with loops because the dynamics becomes frozen due to the priority conflicts. We then consider a scalable interaction protocol, an OR-type one, and examine the effects of the network topology and the number of queues on the waiting time distributions of the priority-queue networks, finding that they exhibit power-law tails in all cases considered, yet with model-dependent power-law exponents. We also show that the synchronicity in task executions, giving rise to priority conflicts in the priority-queue networks, is a relevant factor in the queue dynamics that can change the power-law exponent of the waiting time distribution.

Fractality in complex networks: critical and supercritical skeletons.

Fractal scaling--a power-law behavior of the number of boxes needed to tile a given network with respect to the lateral size of the box--is studied. We introduce a box-covering algorithm that is a modified version of the original algorithm introduced by Song [Nature (London) 433, 392 (2005)]; this algorithm enables easy implementation. Fractal networks are viewed as comprising a skeleton and shortcuts. The skeleton, embedded underneath the original network, is a special type of spanning tree based on the edge betweenness centrality; it provides a scaffold for the fractality of the network. When the skeleton is regarded as a branching tree, it exhibits a plateau in the mean branching number as a function of the distance from a root. For nonfractal networks, on the other hand, the mean branching number decays to zero without forming a plateau. Based on these observations, we construct a fractal network model by combining a random branching tree and local shortcuts. The scaffold branching tree can be either critical or supercritical, depending on the small worldness of a given network. For the network constructed from the critical (supercritical) branching tree, the average number of vertices within a given box grows with the lateral size of the box according to a power-law (an exponential) form in the cluster-growing method. The critical and supercritical skeletons are observed in protein interaction networks and the World Wide Web, respectively. The distribution of box masses, i.e., the number of vertices within each box, follows a power law Pm(M) approximately M(-eta). The exponent eta depends on the box lateral size l(B). For small values of l(B), eta is equal to the degree exponent gamma of a given scale-free network, whereas eta approaches the exponent tau=gamma/(gamma-1) as l(B) increases, which is the exponent of the cluster-size distribution of the random branching tree. Finally, we study the perimeter H(alpha) of a given box alpha, i.e., the number of edges connected to different boxes from a given box alpha as a function of the box mass M(B,alpha). It is obtained that the average perimeter over the boxes with box mass M(B) is likely to scale as approximately M(B), irrespective of the box size l(B).

Structure and evolution of online social relationships: Heterogeneity in unrestricted discussions.

With the advancement in the information age, people are using electronic media more frequently for communications, and social relationships are also increasingly resorting to online channels. While extensive studies on traditional social networks have been carried out, little has been done on online social networks. Here we analyze the structure and evolution of online social relationships by examining the temporal records of a bulletin board system (BBS) in a university. The BBS dataset comprises of 1908 boards, in which a total of 7446 students participate. An edge is assigned to each dialogue between two students, and it is defined as the appearance of the name of a student in the from- and to-field in each message. This yields a weighted network between the communicating students with an unambiguous group association of individuals. In contrast to a typical community network, where intracommunities (intercommunities) are strongly (weakly) tied, the BBS network contains hub members who participate in many boards simultaneously but are strongly tied, that is, they have a large degree and betweenness centrality and provide communication channels between communities. On the other hand, intracommunities are rather homogeneously and weakly connected. Such a structure, which has never been empirically characterized in the past, might provide a new perspective on the social opinion formation in this digital era.

Strength of weak layers in cascading failures on multiplex networks: case of the international trade network.

Many real-world complex systems across natural, social, and economical domains consist of manifold layers to form multiplex networks. The multiple network layers give rise to nonlinear effect for the emergent dynamics of systems. Especially, weak layers that can potentially play significant role in amplifying the vulnerability of multiplex networks might be shadowed in the aggregated single-layer network framework which indiscriminately accumulates all layers. Here we present a simple model of cascading failure on multiplex networks of weight-heterogeneous layers. By simulating the model on the multiplex network of international trades, we found that the multiplex model produces more catastrophic cascading failures which are the result of emergent collective effect of coupling layers, rather than the simple sum thereof. Therefore risks can be systematically underestimated in single-layer network analyses because the impact of weak layers can be overlooked. We anticipate that our simple theoretical study can contribute to further investigation and design of optimal risk-averse real-world complex systems.

Layer-switching cost and optimality in information spreading on multiplex networks.

We study a model of information spreading on multiplex networks, in which agents interact through multiple interaction channels (layers), say online vs. offline communication layers, subject to layer-switching cost for transmissions across different interaction layers. The model is characterized by the layer-wise path-dependent transmissibility over a contact, that is dynamically determined dependently on both incoming and outgoing transmission layers. We formulate an analytical framework to deal with such path-dependent transmissibility and demonstrate the nontrivial interplay between the multiplexity and spreading dynamics, including optimality. It is shown that the epidemic threshold and prevalence respond to the layer-switching cost non-monotonically and that the optimal conditions can change in abrupt non-analytic ways, depending also on the densities of network layers and the type of seed infections. Our results elucidate the essential role of multiplexity that its explicit consideration should be crucial for realistic modeling and prediction of spreading phenomena on multiplex social networks in an era of ever-diversifying social interaction layers.

Extremal dynamics on complex networks: analytic solutions.

The Bak-Sneppen model displaying punctuated equilibria in biological evolution is studied on random complex networks. By using the rate equation and the random walk approaches, we obtain the analytic solution of the fitness threshold xc to be 1/((k)f+1), where (k)f=(k2)/(k) (=(k)) in the quenched (annealed) updating case, where kn is the nth moment of the degree distribution. Thus, the threshold is zero (finite) for the degree exponent gamma<3 (gamma>3) for the quenched case in the thermodynamic limit. The theoretical value xc fits well to the numerical simulation data in the annealed case only. Avalanche size, defined as the duration of successive mutations below the threshold, exhibits a critical behavior as its distribution follows a power law, Pa(s) approximately s(-3/2).

Nonlocal evolution of weighted scale-free networks.

We introduce the notion of globally updating evolution for a class of weighted networks, in which the weight of a link is characterized by the amount of data packet transport flowing through it. By noting that the packet transport over the network is determined nonlocally, this approach can explain the generic nonlinear scaling between the strength and the degree of a node. We demonstrate by a simple model that the strength-driven evolution scheme recently introduced can be generalized to a nonlinear preferential attachment rule, generating the power-law behaviors in degree and in strength simultaneously.

Robustness of the avalanche dynamics in data-packet transport on scale-free networks.

We study the avalanche dynamics in the data-packet transport on scale-free networks through a simple model. In the model, each vertex is assigned a capacity proportional to the load with the proportionality constant 1+a . When the system is perturbed by a single vertex removal, the load of each vertex is redistributed, followed by subsequent failures of overloaded vertices. The avalanche size depends on the parameter a as well as which vertex triggers it. We find that there exists a critical value a(c) at which the avalanche size distribution follows a power law. The critical exponent associated with it appears to be robust as long as the degree exponent is between 2 and 3 and is close in value to that of the distribution of the diameter changes by single vertex removal.

Load distribution in weighted complex networks.

We study the load distribution in weighted networks by measuring the effective number of optimal paths passing through a given vertex. The optimal path, along which the total cost is minimum, crucially depends on the cost distribution function p(c) (c) . In the strong disorder limit, where p(c) (c) approximately c(-1) , the load distribution follows a power law both in the Erdös-Rényi (ER) random graphs and in the scale-free (SF) networks, and its characteristics are determined by the structure of the minimum spanning tree. The distribution of loads at vertices with a given vertex degree also follows the SF nature similar to the whole load distribution, implying that the global transport property is not correlated to the local structural information. Finally, we measure the effect of disorder by the correlation coefficient between vertex degree and load, finding that it is larger for ER networks than for SF networks.

Spectra and eigenvectors of scale-free networks.

We study the spectra and eigenvectors of the adjacency matrices of scale-free networks when bidirectional interaction is allowed, so that the adjacency matrix is real and symmetric. The spectral density shows an exponential decay around the center, followed by power-law long tails at both spectrum edges. The largest eigenvalue lambda1 depends on system size N as lambda1 approximately N1/4 for large N, and the corresponding eigenfunction is strongly localized at the hub, the vertex with largest degree. The component of the normalized eigenfunction at the hub is of order unity. We also find that the mass gap scales as N(-0.68).

Betweenness centrality correlation in social networks.

Scale-free (SF) networks exhibiting a power-law degree distribution can be grouped into the assortative, dissortative, and neutral networks according to the behavior of the degree-degree correlation coefficient. Here we investigate the betweenness centrality (BC) correlation for each type of SF networks. While the BC-BC correlation coefficients behave similarly to the degree-degree correlation coefficients for the dissortative and neutral networks, the BC correlation is nontrivial for the assortative ones found mainly in social networks. The mean BC of neighbors of a vertex with BC g(i) is almost independent of g(i), implying that each person is surrounded by almost the same influential environments of people no matter how influential the person may be.

An evaluation of bond strengths of denture repair resins by a torsional method.

Chemical cure resin materials are generally used in the repair of dentures. Different repair resins used may yield different results. The bond strength of three autopolymerizing resins were evaluated using a torsional test method. The results showed that Palapress and Caulk resins had a higher repair strength than Rapid Repair resin.

Shear bond strength of denture teeth to base: a comparative study.

Fracture and debonding of plastic teeth from the denture base are common clinical problems. In this study, a shear, or peeling, load was applied to the tooth-denture base junction. The bond strength of a high-impact, heat-cured denture base to three commercial brands of denture teeth, Bioform, Dentacryl, and TNR, was assessed. Ninety-three percent of the specimens exhibited cohesive failure within the body of the tooth and not adhesive failure at the tooth-denture base junction. Some of the teeth exhibited brittle fracture while others underwent distortion on loading. Dentacryl displayed the highest cohesive bond strength, followed by Bioform and TNR.

Multiple resource demands and viability in multiplex networks.

Many complex systems demand manifold resources to be supplied from distinct channels to function properly, e.g., water, gas, and electricity for a city. Here, we study a model for viability of such systems demanding more than one type of vital resource be produced and distributed by resource nodes in multiplex networks. We found a rich variety of behaviors such as discontinuity, bistability, and hysteresis in the fraction of viable nodes with respect to the density of networks and the fraction of resource nodes. Our result suggests that viability in multiplex networks is not only exposed to the risk of abrupt collapse but also suffers excessive complication in recovery.