Search

The metapopulation network model is effectively used to study the spatial spread of epidemics with individuals mobility. Considering the time-varying nature of individual activity and the preferences for attractive destinations in population mobility, this paper develops a time-varying network model in which activity of a population is correlated with its attractiveness. Based on the model, the spreading processes of the SIR disease on different correlated networks are studied, and global migration thresholds are derived. It is observed that increasing the correlation between activity and attractiveness results in a reduced outbreak threshold but suppresses the disease outbreak size and introduces greater heterogeneity in the spatial distribution of infected individuals. We also investigate the impact of non-pharmacological interventions (self-isolation and self-protection) on the spread of epidemics in different correlation networks. The results show that the simultaneous implementation of these measures is more effective in negatively correlated networks than in positively correlated or non-correlated networks, and the prevalence is reduced significantly. In addition, both self-isolation and self-protection strategies increase the migration threshold of the spreading and thus slow the spread of the epidemic. However, the effectiveness of each strategy in reducing the density of infected populations varies depending on different correlated networks. Self-protection is more effective in positively correlated networks, whereas self-isolation is more effective in negatively correlated networks. These findings contribute to a better understanding of epidemic spreading in large-scale time-varying metapopulation networks and provide insights for epidemic prevention and control.

In an ecosystem, environmental changes as a result of natural and human processes can cause some key parameters of the system to change with time. Depending on how fast such a parameter changes, a tipping point can occur. Existing works on rate-induced tipping, or R-tipping, offered a theoretical way to study this phenomenon but from a local dynamical point of view, revealing, e.g., the existence of a critical rate for some specific initial condition above which a tipping point will occur. As ecosystems are subject to constant disturbances and can drift away from their equilibrium point, it is necessary to study R-tipping from a global perspective in terms of the initial conditions in the entire relevant phase space region. In particular, we introduce the notion of the probability of R-tipping defined for initial conditions taken from the whole relevant phase space. Using a number of real-world, complex mutualistic networks as a paradigm, we discover a scaling law between this probability and the rate of parameter change and provide a geometric theory to explain the law. The real-world implication is that even a slow parameter change can lead to a system collapse with catastrophic consequences. In fact, to mitigate the environmental changes by merely slowing down the parameter drift may not always be effective: only when the rate of parameter change is reduced to practically zero would the tipping be avoided. Our global dynamics approach offers a more complete and physically meaningful way to understand the important phenomenon of R-tipping., Comment: 8 pages, 5 figures

We propose an efficient and accurate measure for ranking spreaders and identifying the influential ones in spreading processes in networks. While the edges determine the connections among the nodes, their specific role in spreading should be considered explicitly. An edge connecting nodes i and j may differ in its importance for spreading from i to j and from j to i. The key issue is whether node j, after infected by i through the edge, would reach out to other nodes that i itself could not reach directly. It becomes necessary to invoke two unequal weights wij and wji characterizing the importance of an edge according to the neighborhoods of nodes i and j. The total asymmetric directional weights originating from a node leads to a novel measure si which quantifies the impact of the node in spreading processes. A s-shell decomposition scheme further assigns a s-shell index or weighted coreness to the nodes. The effectiveness and accuracy of rankings based on si and the weighted coreness are demonstrated by applying them to nine real-world networks. Results show that they generally outperform rankings based on the nodes' degree and k-shell index, while maintaining a low computational complexity. Our work represents a crucial step towards understanding and controlling the spread of diseases, rumors, information, trends, and innovations in networks., Comment: 9 pages, 8 figures

We investigate fundamental nonlinear dynamics of ferrofluidic Taylor-Couette flow - flow confined between two concentric independently rotating cylinders - consider small aspect ratio by solving the ferrohydrodynamical equations, carrying out systematic bifurcation analysis. Without magnetic field, we find steady flow patterns, previously observed with a simple fluid, such as those containing normal one- or two vortex cells, as well as anomalous one-cell and twin-cell flow states. However, when a symmetry-breaking transverse magnetic field is present, all flow states exhibit stimulated, finite two-fold mode. Various bifurcations between steady and unsteady states can occur, corresponding to the transitions between the two-cell and one-cell states. While unsteady, axially oscillating flow states can arise, we also detect the emergence of new unsteady flow states. In particular, we uncover two new states: one contains only the azimuthally oscillating solution in the configuration of the twin-cell flow state, and another a rotating flow state. Topologically, these flow states are a limit cycle and a quasiperiodic solution on a two-torus, respectively. Emergence of new flow states in addition to observed ones with classical fluid, indicates that richer but potentially more controllable dynamics in ferrofluidic flows, as such flow states depend on the external magnetic field., Comment: 29 pages, 16 figures

Identifying the most influential spreaders is an important issue in controlling the spreading processes in complex networks. Centrality measures are used to rank node influence in a spreading dynamics. Here we propose a node influence measure based on the centrality of a node and its neighbors' centrality, which we call the neighborhood centrality. By simulating the spreading processes in six real-world networks, we find that the neighborhood centrality greatly outperforms the basic centrality of a node such as the degree and coreness in ranking node influence and identifying the most influential spreaders. Interestingly, we discover a saturation effect in considering the neighborhood of a node, which is not the case of the larger the better. Specifically speaking, considering the 2-step neighborhood of nodes is a good choice that balances the cost and performance. If further step of neighborhood is taken into consideration, there is no obvious improvement and even decrease in the ranking performance. The saturation effect may be informative for studies that make use of the local structure of a node to determine its importance in the network., Comment: 17 pages, 8 figures

Recent study shows that the accuracy of the k-shell method in determining node coreness in a spreading process is largely impacted due to the existence of core-like group, which has a large k-shell index but a low spreading efficiency. Based on analysis of the structure of core-like groups in real-world networks, we discover that nodes in the core-like group are mutually densely connected with very few out-leaving links from the group. By defining a measure of diffusion importance for each edge based on the number of out-leaving links of its both ends, we are able to identify redundant links in the spreading process, which have a relatively low diffusion importance but lead to form the locally densely connected core-like group. After filtering out the redundant links and applying the k-shell method to the residual network, we obtain a renewed coreness for each node which is a more accurate index to indicate its location importance and spreading influence in the original network. Moreover, we find that the performance of the ranking algorithms based on the renewed coreness are also greatly enhanced. Our findings help to more accurately decompose the network core structure and identify influential nodes in spreading processes., Comment: 18 pages, 14 figures

It is known that in classical fluids turbulence typically occurs at high Reynolds numbers. But can turbulence occur at low Reynolds numbers? Here we investigate the transition to turbulence in the classic Taylor-Couette system in which the rotating fluids are manufactured ferrofluids with magnetized nanoparticles embedded in liquid carriers. We find that, in the presence of a magnetic field turbulence can occur at Reynolds numbers that are at least one order of magnitude smaller than those in conventional fluids. This is established by extensive computational ferrohydrodynamics through a detailed bifurcation analysis and characterization of behaviors of physical quantities such as the energy, the wave number, and the angular momentum through the bifurcations. A striking finding is that, as the magnetic field is increased, the onset of turbulence can be determined accurately and reliably. Our results imply that experimental investigation of turbulence can be greatly facilitated by using ferrofluids, opening up a new avenue to probe into the fundamentals of turbulence and the challenging problem of turbulence control.

We investigate the quasi two-dimensional Taylor-Couette system in the regime where the radius ratio is close to unity - a transitional regime between three and two dimensions. By systematically increasing the Reynolds number we observe a number of standard transitions, such as one from the classical Taylor vortex flow (TVF) to wavy vortex flow (WVF), as well as the transition to fully developed turbulence. Prior to the onset of turbulence we observe intermittent burst patterns of localized turbulent patches, confirming the experimentally observed pattern of very short wavelength bursts (VSWBs). A striking finding is that, for Reynolds number larger than the onset of VSWBs, a new type of intermittently bursting behaviors emerge: burst patterns of azimuthally closed rings of various orders. We call them ring-burst patterns, which surround the cylinder completely but remain localized and separated by non-turbulent mostly wavy structures in the axial direction. We use a number of quantitative measures, including the cross-flow energy, to characterize the ring-burst patterns and to distinguish them from the background flow. The ring-burst patterns are interesting because it does not occur in either three- or two-dimensional Taylor-Couette flow: it occurs only in the transition, quasi two-dimensional regime of the system, a regime that is less studied but certainly deserves further attention so as to obtain deeper insights into turbulence.

Epidemic threshold is one of the most important features of the epidemic dynamics. Through a lot of numerical simulations in classic Susceptible-Infected-Recovered (SIR) and Susceptible-Infected-Susceptible (SIS) models on various types of networks, we study the simulated identification of epidemic thresholds on finite-size networks. We confirm that the susceptibility measure goes awry for the SIR model due to the bimodal distribution of outbreak sizes near the critical point, while the simulated thresholds of the SIS and SIR models can be accurately determined by analyzing the peak of the epidemic variability. We further verify the accuracy of theoretical predictions derived by the heterogeneous mean-field theory (HMF) and the quenched mean-field theory (QMF), by comparing them with the simulated threshold of the SIR model obtained from the variability measure. The results show that the HMF prediction agrees very well with the simulated threshold, except the case that the networks are disassortive, in which the QMF prediction is more close to the simulated threshold., Comment: 14 pages, 5 figures

Identifying the most influential spreaders is an important issue in understanding and controlling spreading processes on complex networks. Recent studies showed that nodes located in the core of a network as identified by the k-shell decomposition are the most influential spreaders. However, through a great deal of numerical simulations, we observe that not in all real networks do nodes in high shells are very influential: in some networks the core nodes are the most influential which we call true core, while in others nodes in high shells, even the innermost core, are not good spreaders which we call core-like group. By analyzing the k-core structure of the networks, we find that the true core of a network links diversely to the shells of the network, while the core-like group links very locally within the group. For nodes in the core-like group, the k-shell index cannot reflect their location importance in the network. We further introduce a measure based on the link diversity of shells to effectively distinguish the true core and core-like group, and identify core-like groups throughout the networks. Our findings help to better understand the structural features of real networks and influential nodes., Comment: 18 pages, 16 figures

A two-state epidemic model in networks with links mimicking two kinds of relationships between connected nodes is introduced. Links of weights w1 and w0 occur with probabilities p and 1-p, respectively. The fraction of infected nodes rho(p) shows a non-monotonic behavior, with rho drops with p for small p and increases for large p. For small to moderate w1/w0 ratios, rho(p) exhibits a minimum that signifies an optimal suppression. For large w1/w0 ratios, the suppression leads to an absorbing phase consisting only of healthy nodes within a range p_L =< p =< p_R, and an active phase with mixed infected and healthy nodes for p < p_L and p>p_R. A mean field theory that ignores spatial correlation is shown to give qualitative agreement and capture all the key features. A physical picture that emphasizes the intricate interplay between infections via w0 links and within clusters formed by nodes carrying the w1 links is presented. The absorbing state at large w1/w0 ratios results when the clusters are big enough to disrupt the spread via w0 links and yet small enough to avoid an epidemic within the clusters. A theory that uses the possible local environments of a node as variables is formulated. The theory gives results in good agreement with simulation results, thereby showing the necessity of including longer spatial correlations., Comment: 9 pages, 6 figures

The spread of disease on complex networks has attracted widely attention in the physics community. Recent works have demonstrated that heterogeneous degree and weight distributions have a significant influence on the epidemic dynamics. In this study, a novel edge-weight based compartmental approach is developed to estimate the epidemic threshold and epidemic size (final infected density) on networks with general degree and weight distributions, and a remarkable agreement with numerics is obtained. Even in complex network with the strong heterogeneous degree and weight distributions, this approach is worked. We then propose an edge-weight based removal strategy with different biases, and find that such a strategy can effectively control the spread of epidemic when the highly weighted edges are preferentially removed, especially when the weight distribution of a network is extremely heterogenous. The theoretical results from the suggested method can accurately predict the above removal effectiveness., Comment: 9 pages, 5 figures

The spread of disease through a physical-contact network and the spread of information about the disease on a communication network are two intimately related dynamical processes. We investigate the asymmetrical interplay between the two types of spreading dynamics, each occurring on its own layer, by focusing on the two fundamental quantities underlying any spreading process: epidemic threshold and the final infection ratio. We find that an epidemic outbreak on the contact layer can induce an outbreak on the communication layer, and information spreading can effectively raise the epidemic threshold. When structural correlation exists between the two layers, the information threshold remains unchanged but the epidemic threshold can be enhanced, making the contact layer more resilient to epidemic outbreak. We develop a physical theory to understand the intricate interplay between the two types of spreading dynamics., Comment: 29 pages, 14 figures

Recently, the impacts of spatiotemporal heterogeneities of human activities on spreading dynamics have attracted extensive attention. In this paper, to study heterogeneous response times on information spreading, we focus on the susceptible-infected spreading dynamics with adjustable power-law response time distribution based on uncorrelated scale-free networks. We find that the stronger the heterogeneity of response times is, the faster the information spreading is in the early and middle stages. Following a given heterogeneity, the procedure of reducing the correlation between the response times and degrees of individuals can also accelerate the spreading dynamics in the early and middle stages. However, the dynamics in the late stage is slightly more complicated, and there is an optimal value of the full prevalence time changing with the heterogeneity of response times and the response time-degree correlation, respectively. The optimal phenomena results from the efficient allocation of heterogeneous response times., Comment: 8 pages, 6 figures