Showing 1–25 of 25 results for author: Sloot, P

Rate-Induced Transitions in Networked Complex Adaptive Systems: Exploring Dynamics and Management Implications Across Ecological, Social, and Socioecological Systems

Authors: Vítor V. Vasconcelos, Flávia M. D. Marquitti, Theresa Ong, Lisa C. McManus, Marcus Aguiar, Amanda B. Campos, Partha S. Dutta, Kristen Jovanelly, Victoria Junquera, Jude Kong, Elisabeth H. Krueger, Simon A. Levin, Wenying Liao, Mingzhen Lu, Dhruv Mittal, Mercedes Pascual, Flávio L. Pinheiro, Juan Rocha, Fernando P. Santos, Peter Sloot, Chenyang, Su, Benton Taylor, Eden Tekwa, Sjoerd Terpstra , et al. (5 additional authors not shown)

Abstract: Complex adaptive systems (CASs), from ecosystems to economies, are open systems and inherently dependent on external conditions. While a system can transition from one state to another based on the magnitude of change in external conditions, the rate of change -- irrespective of magnitude -- may also lead to system state changes due to a phenomenon known as a rate-induced transition (RIT). This st… ▽ More Complex adaptive systems (CASs), from ecosystems to economies, are open systems and inherently dependent on external conditions. While a system can transition from one state to another based on the magnitude of change in external conditions, the rate of change -- irrespective of magnitude -- may also lead to system state changes due to a phenomenon known as a rate-induced transition (RIT). This study presents a novel framework that captures RITs in CASs through a local model and a network extension where each node contributes to the structural adaptability of others. Our findings reveal how RITs occur at a critical environmental change rate, with lower-degree nodes tipping first due to fewer connections and reduced adaptive capacity. High-degree nodes tip later as their adaptability sources (lower-degree nodes) collapse. This pattern persists across various network structures. Our study calls for an extended perspective when managing CASs, emphasizing the need to focus not only on thresholds of external conditions but also the rate at which those conditions change, particularly in the context of the collapse of surrounding systems that contribute to the focal system's resilience. Our analytical method opens a path to designing management policies that mitigate RIT impacts and enhance resilience in ecological, social, and socioecological systems. These policies could include controlling environmental change rates, fostering system adaptability, implementing adaptive management strategies, and building capacity and knowledge exchange. Our study contributes to the understanding of RIT dynamics and informs effective management strategies for complex adaptive systems in the face of rapid environmental change. △ Less

Submitted 14 September, 2023; originally announced September 2023.

Comments: 25 pages, 4 figures, 1 box, supplementary information

MSC Class: 37G; 37N; 91B; 91C; 91D; 91E; 92D; 92D25; 92D40; 92F; 93A; 93A14; 93A16 ACM Class: I.6.3; I.6.m; J.3; J.4; J.m; K.4.2

Framework for developing quantitative agent based models based on qualitative expert knowledge: an organised crime use-case

Authors: Frederike Oetker, Vittorio Nespeca, Thijs Vis, Paul Duijn, Peter Sloot, Rick Quax

Abstract: In order to model criminal networks for law enforcement purposes, a limited supply of data needs to be translated into validated agent-based models. What is missing in current criminological modelling is a systematic and transparent framework for modelers and domain experts that establishes a modelling procedure for computational criminal modelling that includes translating qualitative data into q… ▽ More In order to model criminal networks for law enforcement purposes, a limited supply of data needs to be translated into validated agent-based models. What is missing in current criminological modelling is a systematic and transparent framework for modelers and domain experts that establishes a modelling procedure for computational criminal modelling that includes translating qualitative data into quantitative rules. For this, we propose FREIDA (Framework for Expert-Informed Data-driven Agent-based models). Throughout the paper, the criminal cocaine replacement model (CCRM) will be used as an example case to demonstrate the FREIDA methodology. For the CCRM, a criminal cocaine network in the Netherlands is being modelled where the kingpin node is being removed, the goal being for the remaining agents to reorganize after the disruption and return the network into a stable state. Qualitative data sources such as case files, literature and interviews are translated into empirical laws, and combined with the quantitative sources such as databases form the three dimensions (environment, agents, behaviour) of a networked ABM. Four case files are being modelled and scored both for training as well as for validation scores to transition to the computational model and application phase respectively. In the last phase, iterative sensitivity analysis, uncertainty quantification and scenario testing eventually lead to a robust model that can help law enforcement plan their intervention strategies. Results indicate the need for flexible parameters as well as additional case file simulations to be performed. △ Less

Submitted 21 July, 2023; originally announced August 2023.

Comments: 47 pages, 10 figures, 6 tables, appendix I-III

Cascades towards noise-induced transitions on networks revealed using information flows

Authors: Casper van Elteren, Rick Quax, Peter Sloot

Abstract: Abrupt, system-wide transitions can be endogenously generated by seemingly stable networks of interacting dynamical units, such as mode switching in neuronal networks or public opinion changes in social systems. However, it remains poorly understood how such `noise-induced transitions' emerge from the interplay of network structure and dynamics on the network. Here, we report on two key roles that… ▽ More Abrupt, system-wide transitions can be endogenously generated by seemingly stable networks of interacting dynamical units, such as mode switching in neuronal networks or public opinion changes in social systems. However, it remains poorly understood how such `noise-induced transitions' emerge from the interplay of network structure and dynamics on the network. Here, we report on two key roles that nodes can play in the progression towards noise-induced tipping points. The models used are dynamical networks where the nodes are governed by the Boltzmann-Gibbs distribution, but the concept is easily generalized. First, so-called `initiator nodes' absorb and then transmit short-lived fluctuations to neighboring nodes, making them temporarily more dynamic. These neighbor nodes can then in turn transmit fluctuations to their neighbors, and so on, leading to a domino-effect where the more stable a node is (i.e., high average free energy barrier), the more neighbors are needed that have become temporarily dynamic. Interestingly, towards the tipping point we identify so-called `stabilizer nodes' whose state information becomes part of the long-term memory of the system, after which the domino-effect is reversed and settles the node in their new stable attractor. We validate these roles by targeted interventions that make tipping points more (or less) likely to begin or lead to systemic change. This opens up possibilities for understanding and controlling endogenously generated metastable behavior. △ Less

Submitted 29 May, 2024; v1 submitted 28 July, 2022; originally announced July 2022.

Comments: Should contain color

Inferring temporal dynamics from cross-sectional data using Langevin dynamics

Authors: Pritha Dutta, Rick Quax, Loes Crielaard, Peter M. A. Sloot

Abstract: Cross-sectional studies are widely prevalent since they are more feasible to conduct compared to longitudinal studies. However, cross-sectional data lack the temporal information required to study the evolution of the underlying processes. Nevertheless, this is essential to develop predictive computational models which is the first step towards causal modelling. We propose a method for inferring c… ▽ More Cross-sectional studies are widely prevalent since they are more feasible to conduct compared to longitudinal studies. However, cross-sectional data lack the temporal information required to study the evolution of the underlying processes. Nevertheless, this is essential to develop predictive computational models which is the first step towards causal modelling. We propose a method for inferring computational models from cross-sectional data using Langevin dynamics. This method can be applied to any system that can be described as effectively following a free energy landscape, such as protein folding, stem cell differentiation and reprogramming, and social systems involving human interaction and social norms. A crucial assumption in our method is that the data-points are gathered from a system in (local) equilibrium. The result is a set of stochastic differential equations which capture the temporal dynamics, by assuming that groups of data-points are subject to the same free energy landscape and amount of noise. Our method is a 'baseline' method which initiates the development of computational models which can be iteratively enhanced through the inclusion of expert knowledge. We validate the proposed method against two population-based longitudinal datasets and observe significant predictive power in comparison with random choice algorithms. We also show how the predictive power of our 'baseline' model can be enhanced by incorporating domain expert knowledge. Our method addresses an important obstacle for model development in fields dominated by cross-sectional datasets. △ Less

Submitted 23 February, 2021; originally announced February 2021.

Comments: 17 pages, 3 figures, "The code for the proposed method is written in Mathematica programming language and is available at https://github.com/Pritha17/langevin-crosssectional"

The dynamic importance of nodes is poorly predicted by static network features

Authors: Casper van Elteren, Rick Quax, Peter Sloot

Abstract: One of the most central questions in network science is: which nodes are most important? Often this question is answered using structural properties such as high connectedness or centrality in the network. However, static structural connectedness does not necessarily translate to dynamical importance. To demonstrate this, we simulate the kinetic Ising spin model on generated networks and one real-… ▽ More One of the most central questions in network science is: which nodes are most important? Often this question is answered using structural properties such as high connectedness or centrality in the network. However, static structural connectedness does not necessarily translate to dynamical importance. To demonstrate this, we simulate the kinetic Ising spin model on generated networks and one real-world weighted network. The dynamic impact of nodes is assessed by causally intervening on node state probabilities and measuring the effect on the systemic dynamics. The results show that structural features such as network centrality or connectedness are actually poor predictors of the dynamical impact of a node on the rest of the network. A solution is offered in the form of an information theoretical measure named integrated mutual information. The metric is able to accurately predict the dynamically most important node ('driver' node) in networks based on observational data of non-intervened dynamics. We conclude that the driver node(s) in networks are not necessarily the most well-connected or central nodes. Indeed, the common assumption of network structural features being proportional to dynamical importance is false. Consequently, great care should be taken when deriving dynamical importance from network data alone. These results highlight the need for novel inference methods that take both structure and dynamics into account. △ Less

Submitted 1 July, 2021; v1 submitted 14 April, 2019; originally announced April 2019.

A New Approach to Privacy-Preserving Clinical Decision Support Systems

Authors: Thomas Attema, Emiliano Mancini, Gabriele Spini, Mark Abspoel, Jan de Gier, Serge Fehr, Thijs Veugen, Maran van Heesch, Daniël Worm, Andrea De Luca, Ronald Cramer, Peter M. A. Sloot

Abstract: Background: Clinical decision support systems (CDSS) are a category of health information technologies that can assist clinicians to choose optimal treatments. These support systems are based on clinical trials and expert knowledge; however, the amount of data available to these systems is limited. For this reason, CDSSs could be significantly improved by using the knowledge obtained by treating p… ▽ More Background: Clinical decision support systems (CDSS) are a category of health information technologies that can assist clinicians to choose optimal treatments. These support systems are based on clinical trials and expert knowledge; however, the amount of data available to these systems is limited. For this reason, CDSSs could be significantly improved by using the knowledge obtained by treating patients. This knowledge is mainly contained in patient records, whose usage is restricted due to privacy and confidentiality constraints. Methods: A treatment effectiveness measure, containing valuable information for treatment prescription, was defined and a method to extract this measure from patient records was developed. This method uses an advanced cryptographic technology, known as secure Multiparty Computation (henceforth referred to as MPC), to preserve the privacy of the patient records and the confidentiality of the clinicians' decisions. Results: Our solution enables to compute the effectiveness measure of a treatment based on patient records, while preserving privacy. Moreover, clinicians are not burdened with the computational and communication costs introduced by the privacy-preserving techniques that are used. Our system is able to compute the effectiveness of 100 treatments for a specific patient in less than 24 minutes, querying a database containing 20,000 patient records. Conclusion: This paper presents a novel and efficient clinical decision support system, that harnesses the potential and insights acquired from treatment data, while preserving the privacy of patient records and the confidentiality of clinician decisions. △ Less

Submitted 3 December, 2018; v1 submitted 2 October, 2018; originally announced October 2018.

Comments: 15 pages, 4 figures

Information processing features can detect behavioral regimes of dynamical systems

Authors: Rick Quax, Gregor Chliamovitch, Alexandre Dupuis, Jean-Luc Falcone, Bastien Chopard, Alfons G. Hoekstra, Peter M. A. Sloot

Abstract: In dynamical systems, local interactions between dynamical units generate correlations which are stored and transmitted throughout the system, generating the macroscopic behavior. However a framework to quantify and study this at the microscopic scale is missing. Here we propose an 'information processing' framework based on Shannon mutual information quantities between the initial and future stat… ▽ More In dynamical systems, local interactions between dynamical units generate correlations which are stored and transmitted throughout the system, generating the macroscopic behavior. However a framework to quantify and study this at the microscopic scale is missing. Here we propose an 'information processing' framework based on Shannon mutual information quantities between the initial and future states. We apply it to the 256 elementary cellular automata (ECA), which are the simplest possible dynamical systems exhibiting behaviors ranging from simple to complex. Our main finding for ECA is that only a few features are needed for full predictability and that the 'information integration' (synergy) feature is always most predictive. Finally we apply the formalism to foreign exchange (FX) and interest-rate swap (IRS) time series data and find that the 2008 financial crisis marks a sudden and sustained regime shift (FX and EUR IRS) resembling tipping point behavior. The USD IRS market exhibits instead a slow and steady progression which appears consistent with the hypothesis that this market is (part of) the driving force behind the crisis. Our work suggests that the proposed framework is a promising way of predicting emergent complex systemic behaviors in terms of the local information processing of units. △ Less

Submitted 27 September, 2017; originally announced September 2017.

Quantifying synergistic information using intermediate stochastic variables

Authors: Rick Quax, Omri Har-Shemesh, Peter M. A. Sloot

Abstract: Quantifying synergy among stochastic variables is an important open problem in information theory. Information synergy occurs when multiple sources together predict an outcome variable better than the sum of single-source predictions. It is an essential phenomenon in biology such as in neuronal networks and cellular regulatory processes, where different information flows integrate to produce a sin… ▽ More Quantifying synergy among stochastic variables is an important open problem in information theory. Information synergy occurs when multiple sources together predict an outcome variable better than the sum of single-source predictions. It is an essential phenomenon in biology such as in neuronal networks and cellular regulatory processes, where different information flows integrate to produce a single response, but also in social cooperation processes as well as in statistical inference tasks in machine learning. Here we propose a metric of synergistic entropy and synergistic information from first principles. The proposed measure relies on so-called synergistic random variables (SRVs) which are constructed to have zero mutual information about individual source variables but non-zero mutual information about the complete set of source variables. We prove several basic and desired properties of our measure, including bounds and additivity properties. In addition, we prove several important consequences of our measure, including the fact that different types of synergistic information may co-exist between the same sets of variables. A numerical implementation is provided, which we use to demonstrate that synergy is associated with resilience to noise. Our measure may be a marked step forward in the study of multivariate information theory and its numerous applications. △ Less

Submitted 19 October, 2016; v1 submitted 3 February, 2016; originally announced February 2016.

ACM Class: E.4; H.1.1; G.3

Information geometric analysis of phase transitions in complex patterns: the case of the Gray-Scott reaction-diffusion model

Authors: Omri Har Shemesh, Rick Quax, Alfons G. Hoekstra, Peter M. A. Sloot

Abstract: The Fisher-Rao metric from Information Geometry is related to phase transition phenomena in classical statistical mechanics. Several studies propose to extend the use of Information Geometry to study more general phase transitions in complex systems. However, it is unclear whether the Fisher-Rao metric does indeed detect these more general transitions, especially in the absence of a statistical mo… ▽ More The Fisher-Rao metric from Information Geometry is related to phase transition phenomena in classical statistical mechanics. Several studies propose to extend the use of Information Geometry to study more general phase transitions in complex systems. However, it is unclear whether the Fisher-Rao metric does indeed detect these more general transitions, especially in the absence of a statistical model. In this paper we study the transitions between patterns in the Gray-Scott reaction-diffusion model using Fisher information. We describe the system by a probability density function that represents the size distribution of blobs in the patterns and compute its Fisher information with respect to changing the two rate parameters of the underlying model. We estimate the distribution non-parametrically so that we do not assume any statistical model. The resulting Fisher map can be interpreted as a phase-map of the different patterns. Lines with high Fisher information can be considered as boundaries between regions of parameter space where patterns with similar characteristics appear. These lines of high Fisher information can be interpreted as phase transitions between complex patterns. △ Less

Submitted 4 January, 2016; v1 submitted 7 December, 2015; originally announced December 2015.

Comments: 13 pages, 8 figures

Topological characterization of S[B] systems: From data to models of complexity

Authors: Emanuela Merelli, Matteo Rucco, Peter Sloot, Luca Tesei

Abstract: In this paper we propose a methodology for deriving a model of a complex system by exploiting the information extracted from Topological Data Analysis. Central to our approach is the S[B] paradigm in which a complex system is represented by a two-level model. One level, the structural S one, is derived using the newly introduced quantitative concept of Persistent Entropy. The other level, the beha… ▽ More In this paper we propose a methodology for deriving a model of a complex system by exploiting the information extracted from Topological Data Analysis. Central to our approach is the S[B] paradigm in which a complex system is represented by a two-level model. One level, the structural S one, is derived using the newly introduced quantitative concept of Persistent Entropy. The other level, the behavioral B one, is characterized by a network of interacting computational agents described by a Higher Dimensional Automaton. The methodology yields also a representation of the evolution of the derived two-level model as a Persistent Entropy Automaton. The presented methodology is applied to a real case study, the Idiotypic Network of the mammal immune system. △ Less

Submitted 25 May, 2015; v1 submitted 21 May, 2015; originally announced May 2015.

Comments: 26 pages, 10 figures

Slope Instability of the Earthen Levee in Boston, UK: Numerical Simulation and Sensor Data Analysis

Authors: N. B. Melnikova, D. Jordan, V. V. Krzhizhanovskaya, P. M. A. Sloot

Abstract: The paper presents a slope stability analysis for a heterogeneous earthen levee in Boston, UK, which is prone to occasional slope failures under tidal loads. Dynamic behavior of the levee under tidal fluctuations was simulated using a finite element model of variably saturated linear elastic perfectly plastic soil. Hydraulic conductivities of the soil strata have been calibrated according to piezo… ▽ More The paper presents a slope stability analysis for a heterogeneous earthen levee in Boston, UK, which is prone to occasional slope failures under tidal loads. Dynamic behavior of the levee under tidal fluctuations was simulated using a finite element model of variably saturated linear elastic perfectly plastic soil. Hydraulic conductivities of the soil strata have been calibrated according to piezometers readings, in order to obtain correct range of hydraulic loads in tidal mode. Finite element simulation was complemented with series of limit equilibrium analyses. Stability analyses have shown that slope failure occurs with the development of a circular slip surface located in the soft clay layer. Both models (FEM and LEM) confirm that the least stable hydraulic condition is the combination of the minimum river levels at low tide with the maximal saturation of soil layers. FEM results indicate that in winter time the levee is almost at its limit state, at the margin of safety (strength reduction factor values are 1.03 and 1.04 for the low-tide and high-tide phases, respectively); these results agree with real-life observations. The stability analyses have been implemented as real-time components integrated into the UrbanFlood early warning system for flood protection. △ Less

Submitted 29 January, 2014; originally announced January 2014.

Quantitative Comparison Between Crowd Models for Evacuation Planning and Evaluation

Authors: Vaisagh Viswanathan, Chong Eu Lee, Michael Harold Lees, Siew Ann Cheong, Peter M. A. Sloot

Abstract: Crowd simulation is rapidly becoming a standard tool for evacuation planning and evaluation. However, the many crowd models in the literature are structurally different, and few have been rigorously calibrated against real-world egress data, especially in emergency situations. In this paper we describe a procedure to quantitatively compare different crowd models or between models and real-world da… ▽ More Crowd simulation is rapidly becoming a standard tool for evacuation planning and evaluation. However, the many crowd models in the literature are structurally different, and few have been rigorously calibrated against real-world egress data, especially in emergency situations. In this paper we describe a procedure to quantitatively compare different crowd models or between models and real-world data. We simulated three models: (1) the lattice gas model, (2) the social force model, and (3) the RVO2 model, and obtained the distributions of six observables: (1) evacuation time, (2) zoned evacuation time, (3) passage density, (4) total distance traveled, (5) inconvenience, and (6) flow rate. We then used the DISTATIS procedure to compute the compromise matrix of statistical distances between the three models. Projecting the three models onto the first two principal components of the compromise matrix, we find the lattice gas and RVO2 models are similar in terms of the evacuation time, passage density, and flow rates, whereas the social force and RVO2 models are similar in terms of the total distance traveled. Most importantly, we find that the zoned evacuation times of the three models to be very different from each other. Thus we propose to use this variable, if it can be measured, as the key test between different models, and also between models and the real world. Finally, we compared the model flow rates against the flow rate of an emergency evacuation during the May 2008 Sichuan earthquake, and found the social force model agrees best with this real data. △ Less

Submitted 10 January, 2014; v1 submitted 1 January, 2014; originally announced January 2014.

Comments: 12 pages, 25 figures, accepted in EPJ B

ACM Class: I.2.11; I.6.4; I.6.5; I.6.6; J.2

Evolutionary Design of Numerical Methods: Generating Finite Difference and Integration Schemes by Differential Evolution

Authors: C. D. Erdbrink, V. V. Krzhizhanovskaya, P. M. A. Sloot

Abstract: Classical and new numerical schemes are generated using evolutionary computing. Differential Evolution is used to find the coefficients of finite difference approximations of function derivatives, and of single and multi-step integration methods. The coefficients are reverse engineered based on samples from a target function and its derivative used for training. The Runge-Kutta schemes are trained… ▽ More Classical and new numerical schemes are generated using evolutionary computing. Differential Evolution is used to find the coefficients of finite difference approximations of function derivatives, and of single and multi-step integration methods. The coefficients are reverse engineered based on samples from a target function and its derivative used for training. The Runge-Kutta schemes are trained using the order condition equations. An appealing feature of the evolutionary method is the low number of model parameters. The population size, termination criterion and number of training points are determined in a sensitivity analysis. Computational results show good agreement between evolved and analytical coefficients. In particular, a new fifth-order Runge-Kutta scheme is computed which adheres to the order conditions with a sum of absolute errors of order 10^-14. Execution of the evolved schemes proved the intended orders of accuracy. The outcome of this study is valuable for future developments in the design of complex numerical methods that are out of reach by conventional means. △ Less

Submitted 30 December, 2013; originally announced December 2013.

Comments: 19 pages, 7 figures, 10 tables, 4 appendices

Categorical and Geographical Separation in Science

Authors: Julian Sienkiewicz, Krzysztof Soja, Janusz A. Holyst, Peter M. A. Sloot

Abstract: We perform the analysis of scientific collaboration at the level of universities. The scope of this study is to answer two fundamental questions: (i) can one indicate a category (i.e., a scientific discipline) that has the greatest impact on the rank of the university and (ii) do the best universities collaborate with the best ones only? Using two university ranking lists (ARWU and QS) as well as… ▽ More We perform the analysis of scientific collaboration at the level of universities. The scope of this study is to answer two fundamental questions: (i) can one indicate a category (i.e., a scientific discipline) that has the greatest impact on the rank of the university and (ii) do the best universities collaborate with the best ones only? Using two university ranking lists (ARWU and QS) as well as data from the Science Citation Index we show how the number of publications in certain categories correlates with the university rank. Moreover, using complex networks analysis, we give hints that the scientific collaboration is highly embedded in the physical space and the number of common papers decays with the distance between them. We also show the strength of the ties between universities is proportional to product of their total number of publications. △ Less

Submitted 15 December, 2013; v1 submitted 2 July, 2013; originally announced July 2013.

Comments: 17 pages, 12 figures, 3 tables

Distributed simulation of city inundation by coupled surface and subsurface porous flow for urban flood decision support system

Authors: V. V. Krzhizhanovskaya, N. B. Melnikova, A. M. Chirkin, S. V. Ivanov, A. V. Boukhanovsky, P. M. A. Sloot

Abstract: We present a decision support system for flood early warning and disaster management. It includes the models for data-driven meteorological predictions, for simulation of atmospheric pressure, wind, long sea waves and seiches; a module for optimization of flood barrier gates operation; models for stability assessment of levees and embankments, for simulation of city inundation dynamics and citizen… ▽ More We present a decision support system for flood early warning and disaster management. It includes the models for data-driven meteorological predictions, for simulation of atmospheric pressure, wind, long sea waves and seiches; a module for optimization of flood barrier gates operation; models for stability assessment of levees and embankments, for simulation of city inundation dynamics and citizens evacuation scenarios. The novelty of this paper is a coupled distributed simulation of surface and subsurface flows that can predict inundation of low-lying inland zones far from the submerged waterfront areas, as observed in St. Petersburg city during the floods. All the models are wrapped as software services in the CLAVIRE platform for urgent computing, which provides workflow management and resource orchestration. △ Less

Submitted 1 February, 2013; originally announced February 2013.

Comments: Pre-print submitted to the 2013 International Conference on Computational Science

Journal ref: Procedia Computer Science, Volume 18, 2013, Pages 1046-1056

Inference of the Russian drug community from one of the largest social networks in the Russian Federation

Authors: L. J. Dijkstra, A. V. Yakushev, P. A. C. Duijn, A. V. Boukhanovsky, P. M. A. Sloot

Abstract: The criminal nature of narcotics complicates the direct assessment of a drug community, while having a good understanding of the type of people drawn or currently using drugs is vital for finding effective intervening strategies. Especially for the Russian Federation this is of immediate concern given the dramatic increase it has seen in drug abuse since the fall of the Soviet Union in the early n… ▽ More The criminal nature of narcotics complicates the direct assessment of a drug community, while having a good understanding of the type of people drawn or currently using drugs is vital for finding effective intervening strategies. Especially for the Russian Federation this is of immediate concern given the dramatic increase it has seen in drug abuse since the fall of the Soviet Union in the early nineties. Using unique data from the Russian social network 'LiveJournal' with over 39 million registered users worldwide, we were able for the first time to identify the on-line drug community by context sensitive text mining of the users' blogs using a dictionary of known drug-related official and 'slang' terminology. By comparing the interests of the users that most actively spread information on narcotics over the network with the interests of the individuals outside the on-line drug community, we found that the 'average' drug user in the Russian Federation is generally mostly interested in topics such as Russian rock, non-traditional medicine, UFOs, Buddhism, yoga and the occult. We identify three distinct scale-free sub-networks of users which can be uniquely classified as being either 'infectious', 'susceptible' or 'immune'. △ Less

Submitted 13 May, 2013; v1 submitted 20 November, 2012; originally announced November 2012.

Comments: 12 pages, 11 figures

Free-surface flow simulations for discharge-based operation of hydraulic structure gates

Authors: C. D. Erdbrink, V. V. Krzhizhanovskaya, P. M. A. Sloot

Abstract: We combine non-hydrostatic flow simulations of the free surface with a discharge model based on elementary gate flow equations for decision support in operation of hydraulic structure gates. A water level-based gate control used in most of today's general practice does not take into account the fact that gate operation scenarios producing similar total discharged volumes and similar water levels m… ▽ More We combine non-hydrostatic flow simulations of the free surface with a discharge model based on elementary gate flow equations for decision support in operation of hydraulic structure gates. A water level-based gate control used in most of today's general practice does not take into account the fact that gate operation scenarios producing similar total discharged volumes and similar water levels may have different local flow characteristics. Accurate and timely prediction of local flow conditions around hydraulic gates is important for several aspects of structure management: ecology, scour, flow-induced gate vibrations and waterway navigation. The modelling approach is described and tested for a multi-gate sluice structure regulating discharge from a river to the sea. The number of opened gates is varied and the discharge is stabilized with automated control by varying gate openings. The free-surface model was validated for discharge showing a correlation coefficient of 0.994 compared to experimental data. Additionally, we show the analysis of CFD results for evaluating bed stability and gate vibrations. △ Less

Submitted 19 November, 2012; originally announced November 2012.

Comments: 25 pages, 16 figures

MSC Class: 76D05

Journal ref: Journal of Hydroinformatics, V. 16, N 1, pp. 189-206, 2014

Continuous Models of Epidemic Spreading in Heterogeneous Dynamically Changing Random Networks

Authors: S. V. Ivanov, A. V. Boukhanovsky, P. M. A. Sloot

Abstract: Modeling spreading processes in complex random networks plays an essential role in understanding and prediction of many real phenomena like epidemics or rumor spreading. The dynamics of such systems may be represented algorithmically by Monte-Carlo simulations on graphs or by ordinary differential equations (ODEs). Despite many results in the area of network modeling the selection of the best comp… ▽ More Modeling spreading processes in complex random networks plays an essential role in understanding and prediction of many real phenomena like epidemics or rumor spreading. The dynamics of such systems may be represented algorithmically by Monte-Carlo simulations on graphs or by ordinary differential equations (ODEs). Despite many results in the area of network modeling the selection of the best computational representation of the model dynamics remains a challenge. While a closed form description is often straightforward to derive, it generally cannot be solved analytically; as a consequence the network dynamics requires a numerical solution of the ODEs or a direct Monte-Carlo simulation on the networks. Moreover, Monte-Carlo simulations and ODE solutions are not equivalent since ODEs produce a deterministic solution while Monte-Carlo simulations are stochastic by nature. Despite some recent advantages in Monte-Carlo simulations, particularly in the flexibility of implementation, the computational cost of an ODE solution is much lower and supports accurate and detailed output analysis such as uncertainty or sensitivity analyses, parameter identification etc. In this paper we propose a novel approach to model spreading processes in complex random heterogeneous networks using systems of nonlinear ordinary differential equations. We successfully apply this approach to predict the dynamics of HIV-AIDS spreading in sexual networks, and compare it to historical data. △ Less

Submitted 19 November, 2012; originally announced November 2012.

Modeling Earthen Dike Stability: Sensitivity Analysis and Automatic Calibration of Diffusivities Based on Live Sensor Data

Authors: N. B. Melnikova, V. V. Krzhizhanovskaya, P. M. A. Sloot

Abstract: The paper describes concept and implementation details of integrating a finite element module for dike stability analysis Virtual Dike into an early warning system for flood protection. The module operates in real-time mode and includes fluid and structural sub-models for simulation of porous flow through the dike and for dike stability analysis. Real-time measurements obtained from pore pressure… ▽ More The paper describes concept and implementation details of integrating a finite element module for dike stability analysis Virtual Dike into an early warning system for flood protection. The module operates in real-time mode and includes fluid and structural sub-models for simulation of porous flow through the dike and for dike stability analysis. Real-time measurements obtained from pore pressure sensors are fed into the simulation module, to be compared with simulated pore pressure dynamics. Implementation of the module has been performed for a real-world test case - an earthen levee protecting a sea-port in Groningen, the Netherlands. Sensitivity analysis and calibration of diffusivities have been performed for tidal fluctuations. An algorithm for automatic diffusivities calibration for a heterogeneous dike is proposed and studied. Analytical solutions describing tidal propagation in one-dimensional saturated aquifer are employed in the algorithm to generate initial estimates of diffusivities. △ Less

Submitted 21 November, 2012; v1 submitted 18 November, 2012; originally announced November 2012.

Journal ref: Journal of Hydrology 496 (2013), pp. 154-165

Are motorways rational from slime mould's point of view?

Authors: Andrew Adamatzky, Selim Akl, Ramon Alonso-Sanz, Wesley van Dessel, Zuwairie Ibrahim, Andrew Ilachinski, Jeff Jones, Anne V. D. M. Kayem, Genaro J. Martinez, Pedro de Oliveira, Mikhail Prokopenko, Theresa Schubert, Peter Sloot, Emanuele Strano, Xin-She Yang

Abstract: We analyse the results of our experimental laboratory approximation of motorways networks with slime mould Physarum polycephalum. Motorway networks of fourteen geographical areas are considered: Australia, Africa, Belgium, Brazil, Canada, China, Germany, Iberia, Italy, Malaysia, Mexico, The Netherlands, UK, USA. For each geographical entity we represented major urban areas by oat flakes and inocul… ▽ More We analyse the results of our experimental laboratory approximation of motorways networks with slime mould Physarum polycephalum. Motorway networks of fourteen geographical areas are considered: Australia, Africa, Belgium, Brazil, Canada, China, Germany, Iberia, Italy, Malaysia, Mexico, The Netherlands, UK, USA. For each geographical entity we represented major urban areas by oat flakes and inoculated the slime mould in a capital. After slime mould spanned all urban areas with a network of its protoplasmic tubes we extracted a generalised Physarum graph from the network and compared the graphs with an abstract motorway graph using most common measures. The measures employed are the number of independent cycles, cohesion, shortest paths lengths, diameter, the Harary index and the Randic index. We obtained a series of intriguing results, and found that the slime mould approximates best of all the motorway graphs of Belgium, Canada and China, and that for all entities studied the best match between Physarum and motorway graphs is detected by the Randic index (molecular branching index). △ Less

Submitted 13 March, 2012; originally announced March 2012.

Simulating City-level Airborne Infectious Diseases

Authors: Mei Shan, Zhou Xuan, Zhu Yifan, Zu Zhenghu, Zheng Tao, A. V. Boukhanovsky, P. M. A Sloot

Abstract: With the exponential growth in the world population and the constant increase in human mobility, the danger of outbreaks of epidemics is rising. Especially in high density urban areas such as public transport and transfer points, where people come in close proximity of each other, we observe a dramatic increase in the transmission of airborne viruses and related pathogens. It is essential to have… ▽ More With the exponential growth in the world population and the constant increase in human mobility, the danger of outbreaks of epidemics is rising. Especially in high density urban areas such as public transport and transfer points, where people come in close proximity of each other, we observe a dramatic increase in the transmission of airborne viruses and related pathogens. It is essential to have a good understanding of the `transmission highways' in such areas, in order to prevent or to predict the spreading of infectious diseases. The approach we take is to combine as much information as is possible, from all relevant sources and integrate this in a simulation environment that allows for scenario testing and decision support. In this paper we lay out a novel approach to study Urban Airborne Disease spreading by combining traffic information, with geo-spatial data, infection dynamics and spreading characteristics. △ Less

Submitted 30 December, 2011; originally announced January 2012.

Comments: 20 pages, 5 figures

The diminishing role of hubs in dynamical processes on complex networks

Authors: Rick Quax, Andrea Apolloni, Peter M. A. Sloot

Abstract: It is notoriously difficult to predict the behaviour of a complex self-organizing system, where the interactions among dynamical units form a heterogeneous topology. Even if the dynamics of each microscopic unit is known, a real understanding of their contributions to the macroscopic system behaviour is still lacking. Here we develop information-theoretical methods to distinguish the contribution… ▽ More It is notoriously difficult to predict the behaviour of a complex self-organizing system, where the interactions among dynamical units form a heterogeneous topology. Even if the dynamics of each microscopic unit is known, a real understanding of their contributions to the macroscopic system behaviour is still lacking. Here we develop information-theoretical methods to distinguish the contribution of each individual unit to the collective out-of-equilibrium dynamics. We show that for a system of units connected by a network of interaction potentials with an arbitrary degree distribution, highly connected units have less impact on the system dynamics as compared to intermediately connected units. In an equilibrium setting, the hubs are often found to dictate the long-term behaviour. However, we find both analytically and experimentally that the instantaneous states of these units have a short-lasting effect on the state trajectory of the entire system. We present qualitative evidence of this phenomenon from empirical findings about a social network of product recommendations, a protein-protein interaction network, and a neural network, suggesting that it might indeed be a widespread property in nature. △ Less

Submitted 30 September, 2013; v1 submitted 23 November, 2011; originally announced November 2011.

Comments: Published version

ACM Class: F.1.1

Journal ref: Journal of the Royal Society Interface, vol. 10, nr 88 2013. ISSN: 1742-5662

Increasing risk behavior can outweigh the benefits of anti-retroviral drug treatment on the HIV incidence among men-having-sex-with-men in Amsterdam

Authors: Shan Mei, Rick Quax, David van de Vijver, Yifan Zhu, A. V. Boukhanovsky, P. M. A. Sloot

Abstract: The transmission through contacts among MSM (men who have sex with men) is one of the dominating contributors to HIV prevalence in industrialized countries. In Amsterdam, the capital of the Netherlands, the MSM risk group has been traced for decades. This has motivated studies which provide detailed information about MSM's risk behavior statistically, psychologically and sociologically. Despite th… ▽ More The transmission through contacts among MSM (men who have sex with men) is one of the dominating contributors to HIV prevalence in industrialized countries. In Amsterdam, the capital of the Netherlands, the MSM risk group has been traced for decades. This has motivated studies which provide detailed information about MSM's risk behavior statistically, psychologically and sociologically. Despite the era of potent antiretroviral therapy, the incidence of HIV among MSM increases. In the long term the contradictory effects of risk behavior and effective therapy are still poorly understood. Using a previously presented Complex Agent Network model, we describe steady and casual partnerships to predict the HIV spreading among MSM. Behavior-related parameters and values, inferred from studies on Amsterdam MSM, are fed into the model; we validate the model using historical yearly incidence data. Subsequently, we study scenarios to assess the contradictory effects of risk behavior and effective therapy, by varying corresponding values of parameters. Finally, we conduct quantitative analysis based on the resulting incidence data. The simulated incidence reproduces the ACS historical incidence well and helps to predict the HIV epidemic among MSM in Amsterdam. Our results show that in the long run the positive influence of effective therapy can be outweighed by an increase in risk behavior of at least 30% for MSM. Conclusion: We recommend, based on the model predictions, that lowering risk behavior is the prominent control mechanism of HIV incidence even in the presence of effective therapy. △ Less

Submitted 28 December, 2010; originally announced December 2010.

Complex Agent Networks explaining the HIV epidemic among homosexual men in Amsterdam

Authors: Shan Mei, P. M. A Sloot, Rick Quax, Yifan Zhu, Weiping Wang

Abstract: Simulating the evolution of the Human Immunodeficiency Virus (HIV) epidemic requires a detailed description of the population network, especially for small populations in which individuals can be represented in detail and accuracy. In this paper, we introduce the concept of a Complex Agent Network(CAN) to model the HIV epidemics by combining agent-based modelling and complex networks, in which a… ▽ More Simulating the evolution of the Human Immunodeficiency Virus (HIV) epidemic requires a detailed description of the population network, especially for small populations in which individuals can be represented in detail and accuracy. In this paper, we introduce the concept of a Complex Agent Network(CAN) to model the HIV epidemics by combining agent-based modelling and complex networks, in which agents represent individuals that have sexual interactions. The applicability of CANs is demonstrated by constructing and executing a detailed HIV epidemic model for men who have sex with men (MSM) in Amsterdam, including a distinction between steady and casual relationships. We focus on MSM contacts because they play an important role in HIV epidemics and have been tracked in Amsterdam for a long time. Our experiments show good correspondence between the historical data of the Amsterdam cohort and the simulation results. △ Less

Submitted 16 February, 2009; v1 submitted 5 December, 2008; originally announced December 2008.

Comments: 21 pages, 4 figures, Mathematics and Computers in Simulation, added references

Towards Distributed Petascale Computing

Authors: A. G. Hoekstra, S. F. Portegies Zwart, M. Bubak, P. M. A. Sloot

Abstract: In this chapter we will argue that studying such multi-scale multi-science systems gives rise to inherently hybrid models containing many different algorithms best serviced by different types of computing environments (ranging from massively parallel computers, via large-scale special purpose machines to clusters of PC's) whose total integrated computing capacity can easily reach the PFlop/s sca… ▽ More In this chapter we will argue that studying such multi-scale multi-science systems gives rise to inherently hybrid models containing many different algorithms best serviced by different types of computing environments (ranging from massively parallel computers, via large-scale special purpose machines to clusters of PC's) whose total integrated computing capacity can easily reach the PFlop/s scale. Such hybrid models, in combination with the by now inherently distributed nature of the data on which the models `feed' suggest a distributed computing model, where parts of the multi-scale multi-science model are executed on the most suitable computing environment, and/or where the computations are carried out close to the required data (i.e. bring the computations to the data instead of the other way around). We presents an estimate for the compute requirements to simulate the Galaxy as a typical example of a multi-scale multi-physics application, requiring distributed Petaflop/s computational power. △ Less

Submitted 19 March, 2007; originally announced March 2007.

Comments: To appear in D. Bader (Ed.) Petascale, Computing: Algorithms and Applications, Chapman & Hall / CRC Press, Taylor and Francis Group