Welcome To Nonlinear Sciences
Robert West, Mauro Mobilia
2019-07-11
Rock-paper-scissors games metaphorically model cyclic dominance in ecology and microbiology. In a static environment, these models are characterized by fixation probabilities obeying two different "laws" in large and small well-mixed populations. Here, we investigate the evolution of these three-species models subject to a randomly switching carrying capacity modeling the endless change between states of resources scarcity and abundance. Focusing mainly on the zero-sum rock-paper-scissors game, equivalent to the cyclic Lotka-Volterra model, we study how the
of demographic and environmental noise influences the fixation properties. More specifically, we investigate which species is the most likely to prevail in a population of fluctuating size and how the outcome depends on the environmental variability. We show that demographic noise coupled with environmental randomness "levels the field" of cyclic competition by balancing the effect of selection. In particular, we show that fast switching effectively reduces the selection intensity proportionally to the variance of the carrying capacity. We determine the conditions under which new fixation scenarios arise, where the most likely species to prevail changes with the rate of switching and the variance of the carrying capacity. Random switching has a limited effect on the mean fixation time that scales linearly with the average population size. Hence, environmental randomness makes the cyclic competition more egalitarian, but does not prolong the species coexistence. We also show how the fixation probabilities of close-to-zero-sum rock-paper-scissors games can be obtained from those of the zero-sum model by rescaling the selection intensity.
Marya Bazzi, Lucas G. S. Jeub, Alex Arenas, Sam D. Howison, Mason A. Porter
2016-08-22
Multilayer networks allow one to represent diverse and coupled connectivity patterns --- e.g., time-dependence, multiple subsystems, or both --- that arise in many applications and which are difficult or awkward to incorporate into standard network representations. In the study of multilayer networks, it is important to investigate mesoscale (i.e., intermediate-scale) structures, such as dense sets of nodes known as communities, to discover network features that are not apparent at the microscale or the macroscale. In this paper, we introduce a generative model for mesoscale structure in multilayer networks. Our model is very general, with the ability to produce many features of empirical multilayer networks, and it explicitly incorporates a user-specified dependency structure between layers. Our results provide a standardized set of null models, together with an associated set of principles from which they are derived, for studies of mesoscale structures in multilayer networks. We discuss the parameters and properties of our generative model, and we illustrate examples of its use with benchmark models for community-detection methods and algorithms in multilayer networks.
Marcos Cardinot, Colm O'Riordan, Josephine Griffith, Attila Szolnoki
2019-07-11
It is generally believed that in a situation where individual and collective interests are in conflict, the availability of optional participation is a key mechanism to maintain cooperation. Surprisingly, this effect is sensitive to the use of microscopic dynamics and can easily be broken when agents make a fully rational decision during their strategy updates. In the framework of the celebrated prisoner's dilemma game, we show that this discrepancy can be fixed automatically if we leave the strict and frequently artifact condition of a fully occupied interaction graph, and allow agents to change not just their strategies but also their positions according to their success. In this way, a diluted graph where agents may move offers a natural and alternative way to handle artifacts arising from the application of specific and sometimes awkward microscopic rules.
Leonardo Rydin Gorjão, Jan Heysel, Klaus Lehnertz, M. Reza Rahimi Tabar
2019-07-11
We introduce the bivariate jump-diffusion process, comprising two-dimensional diffusion and two-dimensional jumps, that can be coupled to one another. We present a data-driven, non-parametric estimation procedure of higher-order Kramers--Moyal coefficients that allows one to reconstruct relevant aspects of the underlying jump-diffusion processes and recover the underlying parameters of a jump-diffusion process. The procedure is validated with numerically integrated data using synthetic bivariate continuous and discontinuous time series. We further evaluate the possibility of estimating the parameters of the jump-diffusion model via a data driven analyses of the higher-order Kramers--Moyal coefficients, and the limitations arising from the scarcity of points in the data or disproportionate parameters in the system.
Benjamin Schäfer, G. Cigdem Yalcin
2019-07-11
A reliable supply of electricity is critical for our modern society and any large scale disturbances of the electrical system causes substantial costs. In 2015, one overloaded transmission line caused a cascading failure in the Turkish power grid, affecting about 75 million people. We focus on the dynamical and statistical properties of the Turkish power grid as a real-world system example that can be modelled based on complex networks structures and we propose for the first time a model that incorporates the dynamical properties of the Turkish power grid as a complex network topology by investigating these cascading failures. We find that the network damage depends on the load and generation distribution in the network with centralized generation being more susceptible to failures than a decentralized one. Furthermore, economic considerations on transmission line capacity are shown to conflict with stability.
Alexander Haluszczynski, Christoph Räth
2019-07-12
The prediction of complex nonlinear dynamical systems with the help of machine learning techniques has become increasingly popular. In particular, reservoir computing turned out to be a very promising approach especially for the reproduction of the long-term properties of a nonlinear system. Yet, a thorough statistical analysis of the forecast results is missing. Using the Lorenz and R"ossler system we statistically analyze the quality of prediction for different parametrizations - both the exact short-term prediction as well as the reproduction of the long-term properties (the "climate") of the system as estimated by the correlation dimension and largest Lyapunov exponent. We find that both short and longterm predictions vary significantly among the realizations. Thus special care must be taken in selecting the good predictions as predictions which deliver better short-term prediction also tend to better resemble the long-term climate of the system. Instead of only using purely random Erd"os-Renyi networks we also investigate the benefit of alternative network topologies such as small world or scale-free networks and show which effect they have on the prediction quality. Our results suggest that the overall performance with respect to the reproduction of the climate of both the Lorenz and R"ossler system is worst for scale-free networks. For the Lorenz system there seems to be a slight benefit of using small world networks while for the R"ossler system small world and Erd"os -Renyi networks performed equivalently well. In general the observation is that reservoir computing works for all network topologies investigated here.
Balachandra Suri, Ravi Kumar Pallantla, Michael F. Schatz, Roman O. Grigoriev
2019-07-11
Recent studies suggest that unstable recurrent solutions of the Navier-Stokes equation provide new insights into dynamics of turbulent flows. In this study, we compute an extensive network of dynamical connections between such solutions in a weakly turbulent quasi-two-dimensional Kolmogorov flow that lies in the inversion-symmetric subspace. In particular, we find numerous isolated heteroclinic connections between different types of solutions -- equilibria, periodic, and quasi-periodic orbits -- as well as continua of connections forming higher-dimensional connecting manifolds. We also compute a homoclinic connection of a periodic orbit and provide strong evidence that the associated homoclinic tangle forms the chaotic repeller that underpins transient turbulence in the symmetric subspace.
Giovanni Conti, Gualtiero Badin
2019-07-11
We investigate the statistical mechanics of a family of two dimensional (2D) fluid flows, described by the generalized Euler equations, or {\alpha}-models. We aim to study the equilibrium mechanics, using initially a point-vortex approximation and then exploiting the full continuous equations, invoking the maximization of appropriate entropy functionals. The point-vortex approximation highlights an important difference between the 2D turbulence and local dynamics models. In the latter, it is in fact possible to derive a statistical measure only considering two conserved quantities as constraints for the maximization problem, the Hamiltonian and the angular impulse. This result does not hold for 2D turbulence. Both the continuous and the point vortex approximation allow for the derivation of mean field equations that act as constraints for the functional relation between the streamfunction and the active scalar of the model considered. Further, the analysis of the continuous equations suggests the existence of a selective decay principle for the whole family of models. To test these ideas we use numerical simulations of the partial differential equations of the {\alpha}-models starting from different sets of initial conditions (i.c.s). For random i.c.s, all the solutions tend to a dipolar structure. The functional relation between the active scalar and the streamfunction shows an increase of nonlinearity with a decrease of the locality of the dynamics. We then test the evolution of the specific case of SQG for i.c.s in the form of a hyperbolic saddle, that is a candidate for the possible formation of singularity through a self-similar cascade though secondary instabilities. Results show the presence of a scale dependent selective decay associated to the breaking of the frontal structures emerging from the flow, suggesting a relation with the change of topology of the flow.
Jianjin Wang, Sergey V. Dmitriev, Daxing Xiong
2019-06-26
The power-law length (
) divergence of thermal conductivity ( ) in one-dimensional (1D) systems, i.e., , has been predicted by theories and also corroborated by experiments. The theoretical predictions of the exponent are usually ranging from to ; however sometimes, the experimental observations can be higher, e.g., - . This dispute has not yet been settled. Here we show the first convincing evidence that an exponent of that falls within experimental observations, can occur in a theoretical model of 1D long-range interacting Fermi-Pasta-Ulam chain. This, for the first time, theoretically supports the possibility of a higher divergent exponent and thus sheds new light on understanding of extremely high thermal conductivity in 1D materials at macroscopic scales.
Xihaier Luo, Ahsan Kareem
2019-03-26
Aerodynamic pressure field over bluff bodies immersed in boundary layer flows is correlated both in space and time. Conventional approaches for the analysis of distributed aerodynamic pressures, e.g., the proper orthogonal decomposition (POD), can only offer relevant spatial patterns in a set of coherent structures. This study provides an operator-theoretic approach that describes dynamic pressure fields in a functional space rather than conventional phase space via the Koopman operator. Subsequently, spectral analysis of the Koopman operator provides a spatiotemporal characterization of the pressure field. An augmented dynamic mode decomposition (DMD) method is proposed to perform the spectral decomposition. The augmentation is achieved by the use of the Takens's embedding theorem, where time delay coordinates are considered. Consequently, the identified eigen-tuples (eigenvalues, eigenvectors, and time evolution) can capture not only dominant spatial structures but also identify each structure with a specific frequency and a corresponding temporal growth/decay. This study encompasses learning the evolution dynamics of the random aerodynamic pressure field over a scaled model of a finite height prism using limited wind tunnel data. The POD analysis of the experimental data was also carried out. To demonstrate the unique feature of the proposed approach, the DMD and POD based learning results including algorithm convergence, data sufficiency, and modal analysis are examined. The ensuing observations offer a glimpse of the complex dynamics of the surface pressure field over bluff bodies that lends insights to features previously masked by conventional analysis approaches.
Lionel Levine, Hanbaek Lyu, John Pike
2017-06-18
In this paper, we consider the soliton cellular automaton introduced in [Takahashi 1990] with a random initial configuration. We give multiple constructions of a Young diagram describing various statistics of the system in terms of familiar objects like birth-and-death chains and Galton-Watson forests. Using these ideas, we establish limit theorems showing that if the first
boxes are occupied independently with probability , then the number of solitons is of order for all , and the length of the longest soliton is of order for , order for , and order for . Additionally, we uncover a condensation phenomenon in the supercritical regime: For each fixed , the top soliton lengths have the same order as the longest for , whereas all but the longest have order at most for . As an application, we obtain scaling limits for the lengths of the longest increasing and decreasing subsequences in a random stack-sortable permutation of length in terms of random walks and Brownian excursions.
Sergio J. Martinez, Ivan M. Mendoza, Genaro J. Martinez, Shigeru Ninagawa
2019-07-06
Universality in cellular automata theory is a central problem studied and developed from their origins by John von Neumann. In this paper, we present an algorithm where any Turing machine can be converted to one-dimensional cellular automaton with a 2-linear time and display its spatial dynamics. Three particular Turing machines are converted in three universal one-dimensional cellular automata, they are: binary sum, rule 110 and a universal reversible Turing machine.
Jan de Gier, Andreas Schadschneider, Johannes Schmidt, Gunter M. Schütz
2019-07-01
Dynamical universality classes are distinguished by their dynamical exponent
and unique scaling functions encoding space-time asymmetry for, e.g. slow-relaxation modes or the distribution of time-integrated currents. So far the universality class of the Nagel-Schreckenberg (NaSch) model, which is a paradigmatic model for traffic flow on highways, was not known except for the special case . Here the model corresponds to the TASEP (totally asymmetric simple exclusion process) that is known to belong to the superdiffusive Kardar-Parisi-Zhang (KPZ) class with . In this paper, we show that the NaSch model also belongs to the KPZ class \cite{KPZ} for general maximum velocities . Using nonlinear fluctuating hydrodynamics theory we calculate the nonuniversal coefficients, fixing the exact asymptotic solutions for the dynamical structure function and the distribution of time-integrated currents. Performing large-scale Monte-Carlo simulations we show that the simulation results match the exact asymptotic KPZ solutions without any fitting parameter left. Additionally, we find that nonuniversal early-time effects or the choice of initial conditions might have a strong impact on the numerical determination of the dynamical exponent and therefore lead to inconclusive results. We also show that the universality class is not changed by extending the model to a two-lane NaSch model with dynamical lane changing rules.
Kazushige Endo
2019-06-29
We propose some conjectures for asymptotic distribution of probabilistic Burgers cellular automaton (PBCA) which is defined by a simple motion rule of particles including a probabilistic parameter. Asymptotic distribution of configurations converges to a unique steady state for PBCA. We assume some conjecture on the distribution and derive the asymptotic probability expressed by GKZ hypergeometric function. If we take a limit of space size to infinity, a relation between density and flux of particles for infinite space size can be evaluated. Moreover, we propose two extended systems of PBCA of which asymptotic behavior can be analyzed as PBCA.
Peter Banda, John Caughman, Martin Cenek, Christof Teuscher
2017-03-27
The search for symmetry as an unusual yet profoundly appealing phenomenon, and the origin of regular, repeating configuration patterns have long been a central focus of complexity science and physics. To better grasp and understand symmetry of configurations in decentralized toroidal architectures, we employ group-theoretic methods, which allow us to identify and enumerate these inputs, and argue about irreversible system behaviors with undesired effects on many computational problems. The concept of so-called configuration shift-symmetry is applied to two-dimensional cellular automata as an ideal model of computation. Regardless of the transition function, the results show the universal insolvability of crucial distributed tasks, such as leader election, pattern recognition, hashing, and encryption. By using compact enumeration formulas and bounding the number of shift-symmetric configurations for a given lattice size, we efficiently calculate the probability of a configuration being shift-symmetric for a uniform or density-uniform distribution. Further, we devise an algorithm detecting the presence of shift-symmetry in a configuration. Given the resource constraints, the enumeration and probability formulas can directly help to lower the minimal expected error and provide recommendations for system's size and initialization. Besides cellular automata, the shift-symmetry analysis can be used to study the non-linear behavior in various synchronous rule-based systems that include inference engines, Boolean networks, neural networks, and systolic arrays.
@complexcity 