Search

The paradigm of interdependent networks has recently been manifested in experimentally testable lab setup of interdependent superconducting networks. This system experiences an abrupt transition due to the thermal dissipation between the networks but its underlying mechanism remains elusive. Here we study the critical behavior and the underlying mechanism of the transition, unveiling its unique microscopic nature. The microscopic characteristics of the transition result in a macroscopic long-living plateau that lasts for thousands of seconds and increases with the size of the system. We characterize the critical behavior of the transition and find that the critical exponents are consistent with those predicted theoretically for percolation of abstract interdependent networks and interdependent ferromagnetic networks, supporting a common universal origin of interdependent systems.

The persistence of the global climate system is critical for assuring the sustainability of the natural ecosystem and the further development of the prosperity of socio-economics. In this paper, we develop a framework and analyze the time persistence of the yearly networks of climate and carbon flux, based on cross-correlations between sites, using daily data from China, the contiguous United States, and the Europe land region during 2000-2019. There are many studies on time persistence of single nodes, e.g., climate variables at a given location, however persistence at a network level has been rarely discussed. Here we develop a framework to study time persistence of network and we apply it to climate and carbon flux. Our framework for determining the persistence is based on analyzing the similarity between the network structures, i.e., the links of climate and carbon flux in different years of systems using the Jaccard index. Our Jaccard results reveal that the similarity of climate and carbon flux networks in different years are within the range of 0.51$\pm$ 0.09 (p-value<0.05), implying that the climate and carbon flux networks studied in the Earth's climate system are generally persistent and in a steady state. Our results suggest that close to 50% of the links appear regularly in different years. We find a very small decay in similarity when the gap between the years increases. However, we observe unique behavior of less similarity to other years in the carbon flux network of the Chinese region during the years 2004-2005 and 2015-2016. This seems to reflect China's carbon reduction policies in these specific years. Analyzing the persistence and evolution of the climate and carbon flux networks, enhance our understanding of the spatial and temporal evolution of the global climate system., Comment: 46 pages, 44 figures

The advancement of quantum communication technologies is calling for a better understanding of quantum network (QN) design from first principles, approached through network science. Pioneering studies have established a classical percolation mapping to model the task of entanglement transmission across QN. Yet, this mapping does not capture the stronger, yet not fully understood connectivity observed in QNs, which facilitates more efficient entanglement transmission than predicted by classical percolation. In this work, we explore the critical phenomena of the potential statistical theory underlying this enhanced connectivity, known as concurrence percolation. Compared to classical percolation, the concurrence percolation mapping employs a unique approach of "superposing" path connectivities, utilizing a different set of path connectivity rules, thereby boosting the overall network connectivity. Firstly, we analytically derive the percolation critical exponents for hierarchical, scale-free networks, particularly the UV flower model, characterized by two distinct network length scales, U$\leq$V. Our analysis confirms that classical and concurrence percolations, albeit both satisfying the hyperscaling relation, fall into separate universality classes. Most importantly, this separation stems from their different treatment of non-shortest path contributions to overall connectivity. Notably, as the longer path scale V increases, concurrence percolation retains unignorable dependence of both its critical threshold and critical exponents on V and thus, comparing with its classical counterpart, shows a higher resilience to the weakening of non-shortest paths. This higher resilience is also observed in real-world network topology, e.g., the Internet. Our findings reveal a first principle for QN design: longer paths still contribute significantly to QN connectivity -- as long as they are abundant., Comment: 16 pages, 4 figures

Earth system tipping elements have been predominantly investigated for their potential to trigger \textit{negative} ecological, climatic, and societal shifts. Yet, an overlooked but seminal avenue exists in the form of \textit{positive} tipping phenomena, whose underlying mechanisms and benefits remain largely underexplored. To bridge this gap, our research introduces a fundamental percolation-based framework to assess the criticality and resilience of planetary terrestrial vegetation systems. Leveraging high-resolution satellite data, we focus on greening-induced positive tipping dynamics driven by global warming. We feature the Qinghai-Tibetan Plateau (QTP) and the Sahel region as contrasting yet analogous case studies. Our analysis uncovers an intriguing phenomenon where vegetation fragmentation aligns with a percolation threshold, exhibiting a scale-invariant pattern characterized by nearly perfect power laws with three critical exponents. Remarkably, contrary to conventional destructive tipping elements, these regions act as favorable tipping elements, transitioning from fragmented to cohesive vegetation patterns due to anthropogenic climate change and afforestation efforts. Furthermore, we propose an \textit{optimal resilience enhancement model} to reinforce vegetation robustness while minimizing socio-economic costs. This study provides valuable insights into the favorable aspects of tipping elements under climate change and offers effective strategies for enhancing ecological resilience against environmental threats.

K-core percolation is a fundamental dynamical process in complex networks with applications that span numerous real-world systems. Earlier studies focus primarily on random networks without spatial constraints and reveal intriguing mixed-order transitions. However, real-world systems, ranging from transportation and communication networks to complex brain networks, are not random but are spatially embedded. Here, we study k-core percolation on two-dimensional spatially embedded networks and show that, in contrast to regular percolation, the length of connections can control the transition type, leading to four different types of phase transitions associated with novel phenomena and a rich phase diagram. A key finding is the existence of a metastable phase in which microscopic localized damage, independent of system size, can cause a macroscopic phase transition, a result which cannot be achieved in traditional percolation. In this case, local failures can spontaneously propagate the damage radially until the system entirely collapses, a phenomenon analogous to the nucleation process. These findings suggest novel features and extreme vulnerabilities of spatially embedded k-core network systems, and highlight the necessity to take into account the characteristic length of links when designing robust spatial networks. Furthermore, our insight about the microscopic processes and their origin during the mixed order and first order abrupt transitions in k-core networks could shed light on the mechanisms of many systems where such transitions occur., Comment: 17 pages, 13 figures

Quantum networks have experienced rapid advancements in both theoretical and experimental domains over the last decade, making it increasingly important to understand their large-scale features from the viewpoint of statistical physics. This review paper discusses a fundamental question: how can entanglement be effectively and indirectly (e.g., through intermediate nodes) distributed between distant nodes in an imperfect quantum network, where the connections are only partially entangled and subject to quantum noise? We survey recent studies addressing this issue by drawing exact or approximate mappings to percolation theory, a branch of statistical physics centered on network connectivity. Notably, we show that the classical percolation frameworks do not uniquely define the network's indirect connectivity. This realization leads to the emergence of an alternative theory called ``concurrence percolation,'' which uncovers a previously unrecognized quantum advantage that emerges at large scales, suggesting that quantum networks are more resilient than initially assumed within classical percolation contexts, offering refreshing insights into future quantum network design.

Mainshocks are often followed by increased earthquake activity (aftershocks). According to the Omori-Utsu law, the rate of aftershocks decays as a power law over time. While aftershocks typically occur in the vicinity of the mainshock, previous studies have suggested that mainshocks can also trigger earthquakes in remote locations. Here we examine the earthquake rate in the days following mega-earthquakes (magnitude >= 7.5) and find that the rate is significantly lower beyond a certain distance from the epicenter compared to surrogate data. However, the remote earthquake rate after the strongest earthquakes (magnitude >= 8) can also be significantly higher than that of the rate based on surrogate data. Comparing our findings to the global ETAS model, we find that the model does not capture the earthquake rate found in the data, hinting at a potential missing mechanism. We suggest that the diminished earthquake rate is due the release of global energy/tension subsequent to substantial mainshock events. This conjecture holds the potential to enhance our comprehension of the intricacies governing post-seismic activity.

While network abrupt breakdowns due to overloads and cascading failures have been studied extensively, the critical exponents and the universality class of such phase transitions have not been discussed. Here, we study breakdowns triggered by failures of links and overloads in networks with a spatial characteristic link length $\zeta$. Our results indicate that this abrupt transition has features and critical exponents similar to those of interdependent networks, suggesting that both systems are in the same universality class. For weakly embedded systems (i.e., $\zeta$ of the order of the system size $L$) we observe a mixed-order transition, where the order parameter collapses following a long critical plateau. On the other hand, strongly embedded systems (i.e., $\zeta \ll L$) exhibit a pure first-order transition, involving nucleation and the growth of damage. The system's critical behavior in both limits is similar to that observed in interdependent networks., Comment: 6 pages, 6 figures

The dynamics of cascading failures in spatial interdependent networks significantly depend on the interaction range of dependency couplings between layers. In particular, for increasing range of dependency couplings, different types of phase transition accompanied by various cascade kinetics can be observed including mixed-order transition characterized by critical branching phenomena, first-order transition with nucleation cascades, and continuous second-order transition with weak cascades. We also describe the dynamics of cascades at the mutual mixed-order resistive transition in interdependent superconductors and show its similarity to that of percolation of interdependent abstract networks. Finally, we layout our perspectives for the experimental observation of these phenomena, their phase diagrams and the underlying kinetics, in the context of physical interdependent networks. Our studies of interdependent networks shed light on the possible mechanisms of three known types of phase transitions, second order, first order, and mixed order as well as predicting a novel fourth type where a microscopic intervention will yield a macroscopic phase transition.

To address the challenge of conflicting traffic flows that complete on opposing cycle times in a specific phase of the traffic light, we proposed a novel decentralized traffic light control methodology based on the identification of the congestion's bottleneck and its influence. By computing the cost of each congested tree, this methodology, called the Tree Method, prioritizes bottlenecks based on global influence, i.e., their impact on the entire road network. To test this methodology, we used SUMO (Simulation of Urban Mobility) as a traffic simulator to run various simulations with both realistic and abstract Origin-Destination (OD) matrices under different traffic loads. By prioritizing traffic flows based on their global cost, rather than merely local ones, the Tree Method accurately identified the root cause of traffic-tree congestion and its effect along the traffic upstream. It exhibited major improvements in both the throughput and the average travel time measurements. Further in-depth analyses of individual simulation runs and time-based performance demonstrated that the Tree Method is superior compared to other methods, not only in terms of average performance but also for the majority of drivers and throughout most of the time. Another significant advantage of this method is its simplicity of analysis, facilitating real-time adjustments in each cycle and aligning with the inherent feedback dynamics of traffic flows. This feature makes the Tree Method highly suitable for integration into actual traffic control systems.

The study of interdependent networks has recently experienced a boost with the development of experimentally testable materials that physically realize their critical behaviors, calling for systematic studies that go beyond the percolation paradigm. Here we study the critical phase transition of interdependent spatial magnetic networks model where dependency couplings between networks are realized by a thermal interaction having a tunable spatial range. We show how the critical phenomena and the phase diagram of this realistic model are highly affected by the range of thermal dissipation and how the latter changes the transition from continuous to abrupt. Furthermore, we show that microscopic interventions of localized heating and localized magnetic field yield a macroscopic phase transition and novel phase diagrams. Our results provide novel and realistic insights about controlling the macroscopic phases of interdependent materials by means of localized microscopic interventions.

The models of $k$-core percolation and interdependent networks (IN) have been extensively studied in their respective fields. A recent study has revealed that they share several common critical exponents. However, several newly discovered exponents in IN have not been explored in $k$-core percolation, and the origin of the similarity still remains unclear. Here, we investigate k-core percolation in random networks. We find that for k-core percolation,the fractality of the giant component fluctuations is manifested by a fractal fluctuation dimension, $\widetilde d_f = 3/4$, within a correlation \emph{size} $N'$ that scales as $N' \propto (p-p_c)^{-\widetilde\nu}$, with $\widetilde\nu = 2$, same as found in IN. Indeed, here, $\widetilde\nu \equiv d\cdot \nu'$ and $\widetilde{d}_f \equiv d'_f/d$, where $\nu'$ and $d'_f$ are respectively the same as the correlation \emph{length} exponent and the fractal fluctuation dimension observed in $d$-dimensional IN spatial networks. These two new exponents found here for $k$-core percolation demonstrate the same scaling behaviors as found for IN with the same critical exponents, reinforcing the similarity between the two models. Furthermore, we suggest that these two models are similar since both have two types of interactions: short-range (SR) connectivity and long-range (LR) influences. In IN the LR are the influences of dependency links while in k-core we find here that for $k=1$ and $k=2$ the influences are short range while for $k\geq3$ the influence is long range. In addition, analytical arguments for a universal hyper-scaling relation for the fractal fluctuation dimension of the $k$-core giant component and for IN as well as for any mixed-order transition are established.Our analysis enhances the comprehension of k-core percolation and supports the generalization of the concept of fractal fluctuations in mixed-order phase transitions.

The impact and originality are two critical dimensions for evaluating scientific publications, measured by citation and disruption metrics respectively. Despite the extensive effort made to understand the statistical properties and evolution of each of these metrics, the relations between the two remain unclear. In this paper, we study the evolution during last 70 years of the correlation between scientific papers' citation and disruption, finding surprisingly a decreasing trend from positive to negative correlations over the years. Consequently, during the years, there are fewer and fewer disruptive works among the highly cited papers. These results suggest that highly disruptive studies nowadays attract less attention from the scientific community. The analysis on papers' references supports this trend, showing that papers citing older references, less popular references and diverse references become to have less citations. Possible explanations for the less attention phenomenon could be due to the increasing information overload in science, and citations become more and more prominent for impact. This is supported by the evidence that research fields with more papers have a more negative correlation between citation and disruption. Finally, we show the generality of our findings by analyzing and comparing six disciplines., Comment: 35 pages, 6+13 figures

Network motifs are the building blocks of complex networks and are significantly involved in the network dynamics such as information processing and local operations in the brain, biological marks for drug targets, identifying and predicting protein complexes in PPI networks, as well as echo chambers in social networks. Here we show that dense motifs such as cliques have different stable states than the network itself. These stable states enhance the dynamical stability of the network and can even turn local stable states into global ones. Moreover, we show how cliques create polarization phenomena and global opinion changes.

According to the Omori-Utsu law, the rate of aftershocks after a mainshock decays as a power law with an exponent close to 1. This well-established law was intensively used in the past to study and model the statistical properties of earthquakes. Moreover, according to the so-called inverse Omori law, the rate of earthquakes should also increase prior to a mainshock -- this law has received much less attention due to its large uncertainty. Here, we mainly study the inverse Omori law based on a highly detailed Southern California earthquake catalog, which is complete for magnitudes larger than M>0.3. First, we develop a technique to identify mainshocks, foreshocks, and aftershocks. We then find, based on a statistical procedure we developed, that the rate of earthquakes is higher a few days prior to a mainshock. We find that this increase is much smaller for a catalog with a magnitude threshold of m over 2.5 and for the Epidemic-Type Aftershocks Sequence (ETAS) model catalogs, even when used with a small magnitude threshold. We also analyze the rate of aftershocks after mainshocks and find that the Omori-Utsu law does not hold for many individual mainshocks and that it may be valid only statistically when considering many mainshocks together. Yet, the analysis of the ETAS model based on the Omori-Utsu law exhibits similar behavior as that of the real catalogs, indicating the validity of this law., Comment: 19 pages, 7 figures

Persistent neural activity associated with working memory (WM) lasts for a limited time duration. Current theories suggest that its termination is \textit{actively} obtained via inhibitory currents, and there is currently no theory regarding the possibility of a \textit{passive} memory-loss mechanism that terminates memory persistent activity. Here, we develop an analytical-framework, based on synaptic strength, and show via simulations and fitting to wet-lab experiments, that passive memory-loss might be a result of an ionic-current long-term plateau, i.e., very slow reduction of memory followed by abrupt loss. We describe analytically the plateau, when the memory state is just below criticality. These results, including the plateau, are supported by experiments performed on rats. Moreover, we show that even just above criticality, forgetfulness can occur due to neuronal noise with ionic-current fluctuations, yielding a plateau, representing memory with very slow decay, and eventually a fast memory decay. Our results could have implications for developing new medications, targeted against memory impairments, through modifying neuronal noise.

The COVID-19 pandemic has evolved over time through multiple spatial and temporal dynamics. The varying extent of interactions among different geographical areas can result to a complex pattern of spreading so that influences between these areas can be hard to discern. Here, we use cross-correlation analysis to detect synchronous evolution and potential inter-influences in the time evolution of new COVID-19 cases at the county level in the USA. Our analysis identified two main time periods with distinguishable features in the behavior of correlations. In the first phase, there were few strong correlations which only emerged between urban areas. In the second phase of the epidemic, strong correlations became widespread and there was a clear directionality of influence from urban to rural areas. In general, the effect of distance between two counties was much weaker than that of the counties population. Such analysis can provide possible clues on the evolution of the disease and may identify parts of the country where intervention may be more efficient in limiting the disease spread., Comment: 18 pages, 7 figures

Tipping elements of the Earth system may shift abruptly and irreversibly from one state to another at tipping points, resulting in a growing threat to our society. Yet, it is not fully clear how to assess and quantify the influence of a tipping element and how to explore the teleconnections between different tipping elements. To fill this knowledge gap, we propose a climate network approach to quantitatively analyze the global impacts of a prominent tipping element, the Amazon Rainforest Area (ARA). We find that regions, such as, the Tibetan Plateau (TP) and West Antarctic ice sheet, are characterized by higher network weighted links and exhibit strong correlations with the ARA. We then identify a teleconnection propagation path between the ARA and the TP. This path is robust under climate change as simulated by various climate models of CMIP5 and CMIP6. In addition, we detect early warning signals for critical transition in the snow cover extent on the Tibetan Plateau by applying critical slowing down indicators, lag-1 autocorrelation and detrended fluctuation analysis. We find that the snow cover of the TP has been losing stability since 2008, revealing that the TP is operating like a tipping element and approaching a potential tipping point. We further uncover that various climate extremes between the ARA and the TP are significantly synchronized under climate change. Our framework provides new insights into how tipping elements are linked to each other and into the potential predictability of cascading tipping dynamics.

In scientific research, collaboration is one of the most effective ways to take advantage of new ideas, skills, resources, and for performing interdisciplinary research. Although collaboration networks have been intensively studied, the question of how individual scientists choose collaborators to study a new research topic remains almost unexplored. Here, we investigate the statistics and mechanisms of collaborations of individual scientists along their careers, revealing that, in general, collaborators are involved in significantly fewer topics than expected from controlled surrogate. In particular, we find that highly productive scientists tend to have higher fraction of single-topic collaborators, while highly cited, i.e., impactful, scientists have higher fraction of multi-topic collaborators. We also suggest a plausible mechanism for this distinction. Moreover, we investigate the cases where scientists involve existing collaborators into a new topic. We find that compared to productive scientists, impactful scientists show strong preference of collaboration with high impact scientists on a new topic. Finally, we validate our findings by investigating active scientists in different years and across different disciplines., Comment: 26+23 pages, 6+18 figures

Dynamical processes on complex networks, ranging from biological, technological and social systems, show phase transitions between distinct global states of the system. Often, such transitions rely upon the interplay between the structure and dynamics that takes place on it, such that weak connectivity, either sparse network or frail interactions, might lead to global activity collapse, while strong connectivity leads to high activity. Here, we show that controlling dynamics of a fraction of the nodes in such systems acts as an external field in a continuous phase transition. As such, it defines corresponding critical exponents, both at equilibrium and of the transient time. We find the critical exponents for a general class of dynamics using the leading orders of the dynamic functions. By applying this framework to three examples, we reveal distinct universality classes.

We study the geometrical features of the order parameter's fluctuations near the critical point of mixed-order phase transitions in randomly interdependent spatial networks. In contrast to continuous transitions, where the structure of the order parameter at criticality is fractal, in mixed-order transitions the structure of the order parameter is known to be compact. Remarkably, we find that although being compact, the fluctuations of the order parameter close to mixed-order transitions are fractal up to a well-defined correlation length $\xi'$, which diverges when approaching the critical threshold. We characterize the self-similar nature of these critical fluctuations through their fractal dimension, $d_f'=3d/4$, and correlation length exponent, $\nu'=2/d$, where $d$ is the dimension of the system. By means of percolation and magnetization, we demonstrate that $d_f'$ and $\nu'$ are independent on the symmetry of the underlying process for any $d$ of the underlying networks.

Many real-world networks are embedded in space, and their resilience in the presence of reinforced nodes has not been studied. Here we model such networks using a spatial network model that have an exponential distribution of link length $r$ having a characteristic length $\zeta$. We find that reinforced nodes can significantly increase the resilience of the networks which varies with strength of spatial embedding. We also study different reinforced node distribution strategies for improving the network resilience. Interestingly, we find that the best strategy is highly dependent on the stage of the percolation process, i.e., the expected fraction of failures. Finally, we show that the reinforced nodes are analogous to an external field in percolation phase transition i.e., having the same critical exponents and that the critical exponents satisfy Widom's relation., Comment: 7 pages, 10 figures

Interdisciplinary research increasingly fuels innovation, and is considered to be a key to tomorrow breakthrough. Yet little is known about whether interdisciplinary research manifests delayed impact. Here, we use the time to reach the citation peak to quantify the highest impact time and citation dynamics, and examine its relationship with interdisciplinarity. Using large scale publication datasets, our results suggest that interdisciplinary papers show significant delayed impact both microscopically per paper and macroscopically collectively, as it takes longer time for interdisciplinary papers to reach their citation peak. Furthermore, we study the underlying forces of such delayed impact, finding that the effect goes beyond the Matthew effect (i.e., the rich-get-richer effect). Finally, we find that team size and content conventionality only partly account for this effect. Overall, our results suggest that governments, research administrators, funding agencies should be aware of this general feature of interdisciplinary science, which may have broad policy implications., Comment: 21 pages, 4 figures

Multi-stability is a widely observed phenomenon in real complex networked systems, such as technological infrastructures, ecological systems, gene regulation, transportation and more. When a system functions normally but there exists also a potential state with abnormal low activity, although the system is at equilibrium it might make a transition into the low activity undesired state due to external disturbances and perturbations. Thus, such a system can be regarded as unsustainable, due to the danger of falling into the potential inactive state. Here we explore, analytically and by simulations, how supporting the activity of a fraction $\rho$ of nodes can turn an unsustainable system to be sustainable by eliminating the inactive potential stable state. We thus unveil a new sustainability phase diagram in the presence of a fraction of controlled nodes $\rho$. This phase diagram could provide guidelines to sustain a network by external intervention and/or by strengthening the connectivity of the network.

Epidemics on complex networks is a widely investigated topic in the last few years, mainly due to the last pandemic events. Usually, real contact networks are dynamic, hence much effort has been invested in studying epidemics on evolving networks. Here we propose and study a model for evolving networks based on varying degrees, where at each time step a node might get, with probability $r$, a new degree and new neighbors according to a given degree distribution, instead of its former neighbors. We find analytically, using the generating functions framework, the epidemic threshold and the probability for a macroscopic spread of disease depending on the rewiring rate $r$. Our analytical results are supported by numerical simulations. We find surprisingly that the impact of the rewiring rate $r$ has qualitative different trends for networks having different degree distributions. That is, in some structures, such as random regular networks the dynamics enhances the epidemic spreading while in others such as scale free the dynamics reduces the spreading. In addition, for scale-free networks, we reveal that fast dynamics of the network, $r=1$, changes the epidemic threshold to nonzero rather than zero found for $r<1$, which is similar to the known case of $r=0$, i.e., a static network. Finally, we find the epidemic threshold also for a general distribution of the recovery time.

Many realistic systems such as infrastructures are characterized by spatial structure and anisotropic alignment. Here we propose and study a model for dealing with such characteristics by introducing a parameter that controls the strength of the anisotropy in the spatial network. This parameter is added to an existing isotropic model used to describe networks under spatial constraints, thus generalizing the spatial model to take into account both spatial and anisotropic features. We study the resilience of such networks by using a percolation process and find that anisotropy has a negative impact on a network's robustness. In addition, our results suggest that the anisotropy in this model does not affect the critical exponent of the correlation length, $\nu$, which remains the same as the known $\nu$ in 2D isotropic lattices., Comment: 9 pages, 5 figures

Many real-world complex systems including human interactions can be represented by temporal (or evolving) networks, where links activate or deactivate over time. Characterizing temporal networks is crucial to compare such systems and to study the dynamical processes unfolding on them. A systematic method to characterize simultaneously the temporal and topological relations of active links (also called contacts or events), in order to compare different real-world networks and to detect their common patterns or differences is still missing. In this paper, we propose a method to characterize to what extent contacts that happen close in time occur also close in topology. Specifically, we study the interrelation between temporal and topological properties of contacts from three perspectives: (1) the autocorrelation of the time series recording the total number of contacts happened at each time step in a network; (2) the interplay between the topological distance and interevent time of two contacts; (3) the temporal correlation of contacts within local neighborhoods beyond a node pair. By applying our method on 13 real-world temporal networks, we found that temporal-topological correlation of contacts is more evident in virtual contact networks than in physical contact ones. This could be due to the lower cost and easier access of online communications than physical interactions, allowing and possibly facilitating social contagion, i.e., interactions of one individual may influence the activity of its neighbors. We also identify different patterns between virtual and physical networks and among physical contact networks at, e.g., school and workplace, in the formation of correlation in local neighborhoods. Detected patterns and differences may further inspire the development of more realistic temporal network models, that could reproduce jointly temporal and topological properties of contacts., Comment: Manuscript accepted by Journal of Complex Networks

The increasing urbanization process we have been witnessing in the last decades is accompanied by a significant increase in traffic congestion in cities around the world. The effect of the congestion is represented in the enormous time people spent on roads leading to significant money waste and air pollution. Here, we present a new methodology for identification, cost evaluation, and thus, prioritization of congestion sources, the jam bottlenecks. It extends existing methods as it is based on network analysis of the entire road network and can be applied to different traffic models. Our results show that the macro-stability, presented by scaling characteristics of the traffic bottlenecks, overshadows the existence of meso-dynamics, where the bottlenecks change their location in time and space. This means that to plan and manage traffic jams in different locations and at different times, it is required to implement a framework, as developed here, that tracks traffic and evaluates the relative effect of each evolving bottleneck on the entire road network.

Interdependence is a fundamental ingredient to analyze the stability of many real-world complex systems featuring functional liasons. Yet, physical realizations of this coupling are still unknown, due to the lack of a theoretical framework for their study. To address this gap, we develop an interdependent magnetization framework and show that dependency links between $K-1$ pairwise networks of Ising spins can be rigorously mapped to directed $K$-spin interactions or to adaptive thermal couplings. We adopt the thermal portrait to determine analytically the phase diagram of the model under different structural configurations and we corroborate our results by extensive simulations. We find that interdependence acts like an entropic force that amplifies site-to-site thermal fluctuations, yielding unusual forms of vulnerability and making the system's functioning often unrecoverable. Finally, we discover an isomorphism between the ground state of random multi-spin models and interdependent percolation on randomly coupled networks. This connection raises new perspectives of cross-fertilization, providing unfamiliar methods with relevant implications in the study of constraint satisfaction as well as to the functional robustness of interdependent systems., Comment: 11 pages, 4 figures

The importance of understanding human mobility patterns has led many studies to examine their spatial-temporal scaling laws. These studies mainly reveal that human travel can be highly non-homogeneous with power-law scaling distributions of distances and times. However, investigating and quantifying the extent of variability in time and space when traveling the same air distance has not been addressed so far. Using taxi data from five large cities, we focus on several novel measures of distance and time to explore the spatio-temporal variations of taxi travel routes relative to their typical routes during peak and nonpeak periods. To compare all trips using a single measure, we calculate the distributions of the ratios between actual travel distances and the average travel distance as well as between actual travel times and the average travel time for all origin destinations (OD) during peak and nonpeak periods. In this way, we measure the scaling of the distribution of all single trip paths with respect to their mean trip path. Our results surprisingly demonstrate very broad distributions for both the distance ratio and time ratio, characterized by a long-tail power-law distribution. Moreover, all analyzed cities have larger exponents in peak hours than in nonpeak hours. We suggest that the interesting results of shorter trip lengths and times, characterized by larger exponents during rush hours, are due to the higher availability of travelers in rush hours compared to non-rush hours...

The performance of distributing entanglement between two distant nodes in a large-scale quantum network (QN) of partially entangled bipartite pure states is generally benchmarked against the classical entanglement percolation (CEP) scheme. Improvements beyond CEP were only achieved by nonscalable strategies for restricted QN topologies. This paper explores and amplifies a new and more effective mapping of a QN, referred to as concurrence percolation theory (ConPT), that suggests using deterministic rather than probabilistic protocols for scalably improving on CEP across arbitrary QN topologies. More precisely, we implement ConPT via a deterministic entanglement transmission (DET) scheme that is fully analogous to resistor network analysis, with the corresponding series and parallel rules represented by deterministic entanglement swapping and concentration protocols, respectively. The main contribution of this paper is to establish a powerful mathematical framework, which is applicable to arbitrary d-dimensional information carriers (qudits), that provides different natural optimality metrics in terms of generalized k-concurrences (a family of fundamental entanglement measures) for different QN topology. In particular, we conclude that the introduced DET scheme (a) is optimal over the well-known nested repeater protocol for distilling entanglement from partially entangled qubits and (b) leads to higher success probabilities of obtaining a maximally entangled state than using CEP. The implementation of the DET scheme is experimentally feasible as tested on IBM's quantum computation platform., Comment: 18 pages, 6 figures. Final version

Local Field potential (LFP) in the basal ganglia (BG) nuclei in the brain have attracted much research and clinical interest. However, the origin of this signal is still under debate throughout the last decades. The question is whether it is a local subthreshold phenomenon, synaptic input to neurons or it is a flow of electrical signals merged as volume conduction which are generated from simultaneous firing neurons in the cerebral cortex and obeys the Maxwell equations. In this study, we recorded in a monkey brain simultaneously LFP's from the cerebral cortex, in the frontal lobe and primary motor cortex (M1) and in sites in all BG nuclei: the striatum, globus pallidus, and subthalamic nucleus. All the records were taken from human primate model (vervet monkey), during spontaneous activity. Developing and applying a novel method to identify significant cross-correlations (potential links) while removing "spurious" correlations, we found a tool that may discriminate between the two major phenomena of synaptic inputs (as we define as information flow) and volume conduction. We find mainly two major paths flows of field potential, that propagates with two different time delays, from the primary motor cortex, and from the frontal cortex. Our results indicate that the two path flows may represent the two mechanisms of volume conduction and information flow.

Here we focus on a basic statistical measure of earthquake catalogs that has not been studied before, the asymmetry of interevent time series (e.g., reflecting the tendency to have more aftershocks than spontaneous earthquakes). We define the asymmetry metric as the ratio between the number of positive interevent time increments minus negative increments and the total (positive plus negative) number of increments. Such asymmetry commonly exists in time series data for non-linear geophysical systems like river flow which decays slowly and increases rapidly. We find that earthquake interevent time series are significantly asymmetric, where the asymmetry function exhibits a significant crossover to weak asymmetry at large lag-index. We suggest that the Omori law can be associated with the large asymmetry at short time intervals below the crossover whereas overlapping aftershock sequences and the spontaneous events can be associated with a fast decay of asymmetry above the crossover. We show that the asymmetry is better reproduced by a recently modified ETAS model with two triggering processes in comparison to the standard ETAS model which only has one.

To reduce the spread and the effect of the COVID-19 global pandemic, non-pharmaceutical interventions (NPIs) have been adopted on multiple occasions by governments. In particular lockdown policies, i.e., generalized mobility restrictions, have been employed to fight the first wave of the pandemic. We analyze data reflecting mobility levels over time in Italy before, during and after the national lockdown, in order to assess some direct and indirect effects. By applying methodologies based on percolation and network science approaches, we find that the typical network characteristics, while very revealing, do not tell the whole story. In particular, the Italian mobility network during lockdown has been damaged much more than node- and edge-level metrics indicate. Additionally, many of the main Provinces of Italy are affected by the lockdown in a surprisingly similar fashion, despite their geographical and economic dissimilarity. Based on our findings we offer an approach to estimate unavailable high-resolution economic dimensions, such as real time Province-level GDP, based on easily measurable mobility information., Comment: 15 pages, 8 figures

Many systems such as critical infrastructure exhibit a modular structure with many links within the modules and few links between them. One approach to increase the robustness of these systems is to reinforce a fraction of the nodes in each module, so that the reinforced nodes provide additional needed sources for themselves as well as for their nearby neighborhood. Since reinforcing a node can be an expensive task, the efficiency of the decentralization process by reinforced nodes is vital. In our study we analyze a new model which combines both above mentioned features of real complex systems - modularity and reinforced nodes. Using tools from percolation theory, we derived an analytical solution for any partition of reinforced nodes; between nodes which have links that connect them to other modules ("inter-nodes") and nodes which have connections only within their modules ("intra-nodes"). Among our results, we find that near the critical percolation point ($p\approx p_c$) the robustness is greatly affected by the distribution. In particular, we find a partition of reinforced nodes which yields an optimal robustness and we show that the optimal partition remains constant for high average degrees.

The structure of real-world multilayer infrastructure systems usually exhibits anisotropy due to constraints of the embedding space. For example, geographical features like mountains, rivers and shores influence the architecture of critical infrastructure networks. Moreover, such spatial networks are often non-homogeneous but rather have a modular structure with dense connections within communities and sparse connections between neighboring communities. When the networks of the different layers are interdependent, local failures and attacks may propagate throughout the system. Here we study the robustness of spatial interdependent networks which are both anisotropic and heterogeneous. We also evaluate the effect of localized attacks having different geometrical shapes. We find that anisotropic networks are more robust against localized attacks and that anisotropic attacks, surprisingly, even on isotropic structures, are more effective than isotropic attacks.

The global spread of the COVID-19 pandemic has followed complex pathways, largely attributed to the high virus infectivity, human travel patterns, and the implementation of multiple mitigation measures. The resulting geographic patterns describe the evolution of the epidemic and can indicate areas which are at risk of an outbreak. Here, we analyze the spatial correlations of new active cases in USA at the county level and characterize the extent of these correlations at different times. We show that the epidemic did not progress uniformly and we identify various stages which are distinguished by significant differences in the correlation length. Our results indicate that the correlation length may be large even during periods when the number of cases declines. We find that correlations between urban centers were much more significant than between rural areas and this finding indicates that long-range spreading was mainly facilitated by travel between cities, especially at the first months of the epidemic. We also show the existence of a percolation transition in November 2020, when the largest part of the country was connected to a spanning cluster, and a smaller-scale transition in January 2021, with both times corresponding to the peak of the epidemic in the country., Comment: 16 pages, 5 figures

Establishing long-distance quantum entanglement, i.e., entanglement transmission, in quantum networks (QN) is a key and timely challenge for developing efficient quantum communication. Traditional comprehension based on classical percolation assumes a necessary condition for successful entanglement transmission between any two infinitely distant nodes: they must be connected by at least a path of perfectly entangled states (singlets). Here, we relax this condition by explicitly showing that one can focus not on optimally converting singlets but on establishing concurrence -- a key measure of bipartite entanglement. We thereby introduce a new statistical theory, concurrence percolation theory (ConPT), remotely analogous to classical percolation but fundamentally different, built by generalizing bond percolation in terms of "sponge-crossing" paths instead of clusters. Inspired by resistance network analysis, we determine the path connectivity by series/parallel rules and approximate higher-order rules via star-mesh transforms. Interestingly, we find that the entanglement transmission threshold predicted by ConPT is lower than the known classical-percolation-based results and is readily achievable on any series-parallel networks such as the Bethe lattice. ConPT promotes our understanding of how well quantum communication can be further systematically improved versus classical statistical predictions under the limitation of QN locality -- a "quantum advantage" that is more general and efficient than expected. ConPT also shows a percolation-like universal critical behavior derived by finite-size analysis on the Bethe lattice and regular two-dimensional lattices, offering new perspectives for a theory of criticality in entanglement statistics., Comment: Phys. Rev. Lett. 126, 170501 (2021)

From mass extinction to cell death, complex networked systems often exhibit abrupt dynamic transitions between desirable and undesirable states. Such transitions are often caused by topological perturbations, such as node or link removal, or decreasing link strengths. The problem is that reversing the topological damage, namely retrieving the lost nodes or links, or reinforcing the weakened interactions, does not guarantee the spontaneous recovery to the desired functional state. Indeed, many of the relevant systems exhibit a hysteresis phenomenon, remaining in the dysfunctional state, despite reconstructing their damaged topology. To address this challenge, we develop a two-step recovery scheme: first - topological reconstruction to the point where the system can be revived, then dynamic interventions, to reignite the system's lost functionality. Applying this method to a range of nonlinear network dynamics, we identify the recoverable phase of a complex system, a state in which the system can be reignited by microscopic interventions, for instance, controlling just a single node. Mapping the boundaries of this dynamical phase, we obtain guidelines for our two-step recovery.

We studied the evolution of the European Patent Office (EPO) patents applicants' collaborations network, within a 35 years span of data (1978-2013). Focusing on the Giant Component (GC) formation process over many time-windows, distributed throughout the data timeline, we found that the features governing this phenomenon are indicative of emerging globalization in the applicants' collaborations. The timeline appears to be divided into three regimes, corresponding to three states of the network's evolution. In the early years state, the GC takes long to form and the instant of its creation is easily pinpointed, while it features geographically segregated groups of applicants with technologically similar activities. In contrast, in the late years state, the GC forms quickly, the exact point of its creation is harder to spot, the applicants' activities are more disparate technologically, while their inter-regional collaborations are significantly increased. The middle years are an intermediate state between the two extreme of early and late years. Moreover, we concluded that the critical patents, which induce the GC's formation, are typically introduced by large-sized applicants and also that top patent-producing applicants are likely to submit critical patents, albeit at a lower rate than their overall patent submission. Lastly, we uncovered the crucial role that Japan plays in the network's coherence, through its prominent participation in the GC and the critical patents., Comment: 20 pages, 19 figures

Epidemic spread on networks is one of the most studied dynamics in network science and has important implications in real epidemic scenarios. Nonetheless, the dynamics of real epidemics and how it is affected by the underline structure of the infection channels are still not fully understood. Here we apply the SIR model and study analytically and numerically the epidemic spread on a recently developed spatial modular model imitating the structure of cities in a country. The model assumes that inside a city the infection channels connect many different locations, while the infection channels between cities are less and usually directly connect only a few nearest neighbor cities in a two-dimensional plane. We find that the model experience two epidemic transitions. The first lower threshold represents a local epidemic spread within a city but not to the entire country and the second higher threshold represents a global epidemic in the entire country. Based on our analytical solution we proposed several control strategies and how to optimize them. We also show that while control strategies can successfully control the disease, early actions are essentials to prevent the disease global spread., Comment: 8 pages, 7 figures

Global climate change, extreme climate events, earthquakes and their accompanying natural disasters pose significant risks to humanity. Yet due to the nonlinear feedbacks, strategic interactions and complex structure of the Earth system, the understanding and in particular the predicting of such disruptive events represent formidable challenges for both scientific and policy communities. During the past years, the emergence and evolution of Earth system science has attracted much attention and produced new concepts and frameworks. Especially, novel statistical physics and complex networks-based techniques have been developed and implemented to substantially advance our knowledge for a better understanding of the Earth system, including climate extreme events, earthquakes and Earth geometric relief features, leading to substantially improved predictive performances. We present here a comprehensive review on the recent scientific progress in the development and application of how combined statistical physics and complex systems science approaches such as, critical phenomena, network theory, percolation, tipping points analysis, as well as entropy can be applied to complex Earth systems (climate, earthquakes, etc.). Notably, these integrating tools and approaches provide new insights and perspectives for understanding the dynamics of the Earth systems. The overall aim of this review is to offer readers the knowledge on how statistical physics approaches can be useful in the field of Earth system science.