Search Results
Ensemble analysis of complex network properties—an MCMC approach
2022, Pfeffer, Oskar, Molkenthin, Nora, Hellmann, Frank
What do generic networks that have certain properties look like? We use relative canonical network ensembles as the ensembles that realize a property R while being as indistinguishable as possible from a background network ensemble. This allows us to study the most generic features of the networks giving rise to the property under investigation. To test the approach we apply it to study properties thought to characterize ‘small-world networks’. We consider two different defining properties, the ‘small-world-ness’ of Humphries and Gurney, as well as a geometric variant. Studying them in the context of Erdős-Rényi and Watts-Strogatz ensembles we find that all ensembles studied exhibit phase transitions to systems with large hubs and in some cases cliques. Such features are not present in common examples of small-world networks, indicating that these properties do not robustly capture the notion of small-world networks. We expect the overall approach to have wide applicability for understanding network properties of real world interest, such as optimal ride-sharing designs, the vulnerability of networks to cascades, the performance of communication topologies in coordinating fluctuation response or the ability of social distancing measures to suppress disease spreading.
Efficient network immunization strategy based on generalized Herfindahl–Hirschman index
2021, Chen, Peng, Qi, Mingze, Lu, Xin, Duan, Xiaojun, Kurths, Jürgen
The topic of finding effective strategies to restrain epidemic spreading in complex networks is of current interest. A widely used approach for epidemic containment is the fragmentation of the contact networks through immunization. However, due to the limitation of immune resources, we cannot always fragment the contact network completely. In this study, based on the size distribution of connected components for the network, we designed a risk indicator of epidemic outbreaks, the generalized Herfindahl–Hirschman index (GHI), which measures the upper bound of the expected infection's prevalence (the fraction of infected nodes) in random outbreaks. An immunization approach based on minimizing GHI is developed to reduce the infection risk for individuals in the network. Experimental results show that our immunization strategy could effectively decrease the infection's prevalence as compared to other existing strategies, especially against infectious diseases with higher infection rates or lower recovery rates. The findings provide an efficient and practicable strategy for immunization against epidemic diseases.
Neural partial differential equations for chaotic systems
2021, Gelbrecht, Maximilian, Boers, Niklas, Kurths, Jürgen
When predicting complex systems one typically relies on differential equation which can often be incomplete, missing unknown influences or higher order effects. By augmenting the equations with artificial neural networks we can compensate these deficiencies. We show that this can be used to predict paradigmatic, high-dimensional chaotic partial differential equations even when only short and incomplete datasets are available. The forecast horizon for these high dimensional systems is about an order of magnitude larger than the length of the training data.
Particle dynamics and transport enhancement in a confined channel with position-dependent diffusivity
2020, Li, Yongge, Mei, Ruoxing, Xu, Yong, Kurths, Jürgen, Duan, Jinqiao, Metzler, Ralf
This work focuses on the dynamics of particles in a confined geometry with position-dependent diffusivity, where the confinement is modelled by a periodic channel consisting of unit cells connected by narrow passage ways. We consider three functional forms for the diffusivity, corresponding to the scenarios of a constant (D 0), as well as a low (D m) and a high (D d) mobility diffusion in cell centre of the longitudinally symmetric cells. Due to the interaction among the diffusivity, channel shape and external force, the system exhibits complex and interesting phenomena. By calculating the probability density function, mean velocity and mean first exit time with the Itô calculus form, we find that in the absence of external forces the diffusivity D d will redistribute particles near the channel wall, while the diffusivity D m will trap them near the cell centre. The superposition of external forces will break their static distributions. Besides, our results demonstrate that for the diffusivity D d, a high dependence on the x coordinate (parallel with the central channel line) will improve the mean velocity of the particles. In contrast, for the diffusivity D m, a weak dependence on the x coordinate will dramatically accelerate the moving speed. In addition, it shows that a large external force can weaken the influences of different diffusivities; inversely, for a small external force, the types of diffusivity affect significantly the particle dynamics. In practice, one can apply these results to achieve a prominent enhancement of the particle transport in two- or three-dimensional channels by modulating the local tracer diffusivity via an engineered gel of varying porosity or by adding a cold tube to cool down the diffusivity along the central line, which may be a relevant effect in engineering applications. Effects of different stochastic calculi in the evaluation of the underlying multiplicative stochastic equation for different physical scenarios are discussed. © 2020 The Author(s). Published by IOP Publishing Ltd on behalf of the Institute of Physics and Deutsche Physikalische Gesellschaft.
Improved earthquake aftershocks forecasting model based on long-term memory
2021, Zhang, Yongwen, Zhou, Dong, Fan, Jingfang, Marzocchi, Warner, Ashkenazy, Yosef, Havlin, Shlomo
A prominent feature of earthquakes is their empirical laws, including memory (clustering) in time and space. Several earthquake forecasting models, such as the epidemic-type aftershock sequence (ETAS) model, were developed based on these empirical laws. Yet, a recent study [1] showed that the ETAS model fails to reproduce the significant long-term memory characteristics found in real earthquake catalogs. Here we modify and generalize the ETAS model to include short- and long-term triggering mechanisms, to account for the short- and long-time memory (exponents) discovered in the data. Our generalized ETAS model accurately reproduces the short- and long-term/distance memory observed in the Italian and Southern Californian earthquake catalogs. The revised ETAS model is also found to improve earthquake forecasting after large shocks.
Evolution mechanism of principal modes in climate dynamics
2020, Zhang, Yongwen, Fan, Jingfang, Li, Xiaoteng, Liu, Wenqi, Chen, Xiaosong
Eigen analysis has been a powerful tool to distinguish multiple processes into different simple principal modes in complex systems. For a non-equilibrium system, the principal modes corresponding to the non-equilibrium processes are usually evolving with time. Here, we apply the eigen analysis into the complex climate systems. In particular, based on the daily surface air temperature in the tropics (30? S–30? N, 0? E–360? E) between 1979-01-01 and 2016-12-31, we uncover that the strength of two dominated intra-annual principal modes represented by the eigenvalues significantly changes with the El Niño/southern oscillation from year to year. Specifically, according to the ‘regional correlation’ introduced for the first intra-annual principal mode, we find that a sharp positive peak of the correlation between the El Niño region and the northern (southern) hemisphere usually signals the beginning (end) of the El Niño. We discuss the underlying physical mechanism and suppose that the evolution of the first intra-annual principal mode is related to the meridional circulations; the evolution of the second intra-annual principal mode responds positively to the Walker circulation. Our framework presented here not only facilitates the understanding of climate systems but also can potentially be used to study the dynamical evolution of other natural or engineering complex systems. © 2020 The Author(s).
Dynamical network size estimation from local observations
2020, Tang, Xiuchuan, Huo, Wei, Yuan, Ye, Li, Xiuting, Shi, Ling, Kurths, Jürgen
Here we present a method to estimate the total number of nodes of a network using locally observed response dynamics. The algorithm has the following advantages: (a) it is data-driven. Therefore it does not require any prior knowledge about the model; (b) it does not need to collect measurements from multiple stimulus; and (c) it is distributed as it uses local information only, without any prior information about the global network. Even if only a single node is measured, the exact network size can be correctly estimated using a single trajectory. The proposed algorithm has been applied to both linear and nonlinear networks in simulation, illustrating the applicability to real-world physical networks. © 2020 The Author(s). Published by IOP Publishing Ltd on behalf of the Institute of Physics and Deutsche Physikalische Gesellschaft.
Monte Carlo basin bifurcation analysis
2020, Gelbrecht, Maximilian, Kurths, Jürgen, Hellmann, Frank
Many high-dimensional complex systems exhibit an enormously complex landscape of possible asymptotic states. Here, we present a numerical approach geared towards analyzing such systems. It is situated between the classical analysis with macroscopic order parameters and a more thorough, detailed bifurcation analysis. With our machine learning method, based on random sampling and clustering methods, we are able to characterize the different asymptotic states or classes thereof and even their basins of attraction. In order to do this, suitable, easy to compute, statistics of trajectories with randomly generated initial conditions and parameters are clustered by an algorithm such as DBSCAN. Due to its modular and flexible nature, our method has a wide range of possible applications in many disciplines. While typical applications are oscillator networks, it is not limited only to ordinary differential equation systems, every complex system yielding trajectories, such as maps or agent-based models, can be analyzed, as we show by applying it the Dodds-Watts model, a generalized SIRS-model, modeling social and biological contagion. A second order Kuramoto model, used, e.g. to investigate power grid dynamics, and a Stuart-Landau oscillator network, each exhibiting a complex multistable regime, are shown as well. The method is available to use as a package for the Julia language. © 2020 The Author(s). Published by IOP Publishing Ltd on behalf of the Institute of Physics and Deutsche Physikalische Gesellschaft.
A recurrent plot based stochastic nonlinear ray propagation model for underwater signal propagation
2020, Haiyang, Yao, Haiyan, Wang, Yong, Xu, Kurths, Juergen
A stochastic nonlinear ray propagation model is proposed to carry out an exploration of the nonlinear ray theory in underwater signal propagation. The recurrence plot method is proposed to quantify the ray chaos and stochastics to optimize the model. Based on this method, the distribution function of the control parameter d is derived. Experiments and simulations indicate that this stochastic nonlinear ray propagation model provides a good explanation and description on the stochastic frequency shift in underwater signal propagation. © 2020 The Author(s). Published by IOP Publishing Ltd on behalf of the Institute of Physics and Deutsche Physikalische Gesellschaft.
Basin stability and limit cycles in a conceptual model for climate tipping cascades
2020, Wunderling, Nico, Gelbrecht, Maximilian, Winkelmann, Ricarda, Kurths, Jürgen, Donges, Jonathan F.
Tipping elements in the climate system are large-scale subregions of the Earth that might possess threshold behavior under global warming with large potential impacts on human societies. Here, we study a subset of five tipping elements and their interactions in a conceptual and easily extendable framework: the Greenland Ice Sheets (GIS) and West Antarctic Ice Sheets, the Atlantic meridional overturning circulation (AMOC), the El–Niño Southern Oscillation and the Amazon rainforest. In this nonlinear and multistable system, we perform a basin stability analysis to detect its stable states and their associated Earth system resilience. By combining these two methodologies with a large-scale Monte Carlo approach, we are able to propagate the many uncertainties associated with the critical temperature thresholds and the interaction strengths of the tipping elements. Using this approach, we perform a system-wide and comprehensive robustness analysis with more than 3.5 billion ensemble members. Further, we investigate dynamic regimes where some of the states lose stability and oscillations appear using a newly developed basin bifurcation analysis methodology. Our results reveal that the state of four or five tipped elements has the largest basin volume for large levels of global warming beyond 4 °C above pre-industrial climate conditions, representing a highly undesired state where a majority of the tipping elements reside in the transitioned regime. For lower levels of warming, states including disintegrated ice sheets on west Antarctica and Greenland have higher basin volume than other state configurations. Therefore in our model, we find that the large ice sheets are of particular importance for Earth system resilience. We also detect the emergence of limit cycles for 0.6% of all ensemble members at rare parameter combinations. Such limit cycle oscillations mainly occur between the GIS and AMOC (86%), due to their negative feedback coupling. These limit cycles point to possibly dangerous internal modes of variability in the climate system that could have played a role in paleoclimatic dynamics such as those unfolding during the Pleistocene ice age cycles.