Welcome To Nonlinear Sciences
Ivan A. Korneev, Andrei V. Slepnev, Vladimir V. Semenov, Tatiana Vadivasova
2019-07-15
The model of a non-autonomous memristor-based oscillator with a line of equilibria is studied. A numerical simulation of the system driven by a periodical force is combined with a theoretical analysis by means of the quasi-harmonic reduction. Both two mechanisms of synchronization are demonstrated: capture of the phase and frequency of oscillations and suppression by an external signal. Classification of undamped oscillations in an autonomous system with a line of equilibria as a special kind of self-sustained oscillations is concluded due to the possibility to observe the effect of frequency-phase locking in the same system in the presence of an external influence. It is established that the occurrence of phase locking in the considered system continuously depends both on parameter values and initial conditions. The simultaneous dependence of synchronization area boundaries on the initial conditions and the parameter values is also shown.
Ivan A. Korneev, Vladimir V. Semenov, Tatiana E. Vadivasova
2018-07-02
A model of two self-sustained oscillators interacting through memristive coupling is studied. Memristive coupling is realized by using a cubic memristor model. Numerical simulation is combined with theoretical analysis by means of quasi-harmonic reduction. It is shown that specifics of the memristor nonlinearity results in appearance of infinitely many equilibrium points, which form a line of equilibria in the phase space of the system under study. It is established that possibility to observe the effect of phase locking in the considered system depends both on parameter values and initial conditions. Consequently, boundaries of a synchronization area are determined by the initial conditions. It is demonstrated that addition of a small term into the memristor state equation gives rise to disappearance of the line of equilibria and eliminates the dependence of synchronization on the initial conditions.
Seth Frey, Curtis Atkisson
2019-07-15
Understanding the evolution of human social systems requires flexible formalisms for the emergence of institutions. Although game theory is normally used to model interactions individually, larger spaces of games can be helpful for modeling how interactions change. We introduce a framework for modeling "institutional evolution," how individuals change the games they are placed in. We contrast this with the more familiar within-game "behavioral evolution". Starting from an initial game, agents trace trajectories through game space by repeatedly navigating to more preferable games until they converge on attractor games that are preferred to all others. Agents choose between games on the basis of their "institutional preferences," which define between-game comparisons in terms of game-level features such as stability, fairness, and efficiency. Computing institutional change trajectories over the two-player space, we find that the attractors of self-interested economic agents over-represent fairness by 100% relative to baseline, even though those agents are indifferent to fairness. This seems to occur because fairness, as a game feature, co-occurs with the self-serving features these agents do prefer. We thus present institutional evolution as a mechanism for encouraging the spontaneous emergence of cooperation among inherently selfish agents. We then extend these findings beyond two players, and to two other types of evolutionary agent: the relative fitness maximizing agent of evolutionary game theory (who maximizes inequality), and the relative group fitness maximizing agent of multi-level/group selection theory (who minimizes inequality). This work provides a flexible, testable formalism for modeling the interdependencies of behavioral and institutional evolutionary processes.
Walid Ahmad, Mason A. Porter, Mariano Beguerisse-Díaz
2018-05-01
Network theory is a useful framework for studying interconnected systems of interacting agents. Many networked systems evolve continuously in time, but most existing methods for the analysis of time-dependent networks rely on discrete or discretized time. In this paper, we propose an approach for studying networks that evolve in continuous time by distinguishing between interactions, which we model as discrete contacts, and ties, which represent strengths of relationships as functions of time. To illustrate our tie-decay network formulation, we adapt the well-known PageRank centrality score to the tie-decay framework in a mathematically tractable and computationally efficient way. We demonstrate our framework on a synthetic example and then use it to study a network of retweets during the 2012 National Health Service controversy in the United Kingdom. Our work also provides guidance for similar generalizations of other tools from network theory to continuous-time networks with tie decay, including for applications to streaming data.
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.
Ivan A. Korneev, Vladimir V. Semenov, Tatiana E. Vadivasova
2018-07-02
A model of two self-sustained oscillators interacting through memristive coupling is studied. Memristive coupling is realized by using a cubic memristor model. Numerical simulation is combined with theoretical analysis by means of quasi-harmonic reduction. It is shown that specifics of the memristor nonlinearity results in appearance of infinitely many equilibrium points, which form a line of equilibria in the phase space of the system under study. It is established that possibility to observe the effect of phase locking in the considered system depends both on parameter values and initial conditions. Consequently, boundaries of a synchronization area are determined by the initial conditions. It is demonstrated that addition of a small term into the memristor state equation gives rise to disappearance of the line of equilibria and eliminates the dependence of synchronization on the initial conditions.
Rosario González-Férez, Manuel Iñarrea, J. Pablo Salas, Peter Schmelcher
2019-07-15
We investigate the structure and the nonlinear dynamics of two rigid polar rotors coupled through the dipole-dipole interaction in an external homogeneous electric field. In the field-free stable head-tail configuration, an excess energy is provided to one of the dipoles, and we explore the resulting three-dimensional classical dynamics. This dynamics is characterized in terms of the kinetic energy transfer between the dipoles, their orientation along the electric field, as well as their chaotic behavior. The field-free energy transfer mechanism shows an abrupt transition between equipartition and non-equipartition regimes, which is independent of the initial direction of rotation due to the existence of an infinite set of equivalent manifolds. The field-dressed dynamics is highly complex and strongly depends on the electric field strength and on the initial conditions. In the strong field regime, the energy equipartition and chaotic behavior dominate the dynamics.
Gabriel M. Lando, Alfredo M. Ozorio de Almeida
2019-07-15
An initial coherent state is propagated exactly by a kicked quantum Hamiltonian and its associated classical stroboscopic map. The classical trajectories within the initial state are regular for low kicking strengths, then bifurcate and become mainly chaotic as the kicking parameter is increased. Time-evolution is tracked using classical, quantum and semiclassical Wigner functions, obtained via the Herman-Kluk propagator. Quantitative comparisons are also included and carried out from probability marginals and autocorrelation functions. Sub-Planckian classical structure such as small stability islands and thin/folded classical filaments do impact semiclassical accuracy, but the approximation is seen to be accurate for multiple Ehrenfest times.
Stéphane Vannitsem, Roman Solé-Pomies, Lesley De Cruz
2019-01-18
The predictability of the atmosphere at short and long time scales, associated with the coupling to the ocean, is explored in a new version of the Modular Arbitrary-Order Ocean-Atmosphere Model (MAOOAM), based on a 2-layer quasi-geostrophic atmosphere and a 1-layer reduced-gravity quasi-geostrophic ocean. This version features a new ocean basin geometry with periodic boundary conditions in the zonal direction. The analysis presented in this paper considers a low-order version of the model with 40 dynamical variables. First the increase of surface friction (and the associated heat flux) with the ocean can either induce chaos when the aspect ratio between the meridional and zonal directions of the domain of integration is small, or suppress chaos when it is large. This reflects the potentially counter-intuitive role that the ocean can play in the coupled dynamics. Second, and perhaps more importantly, the emergence of long-term predictability within the atmosphere for specific values of the friction coefficient occurs through intermittent excursions in the vicinity of a (long-period) unstable periodic solution. Once close to this solution the system is predictable for long times, i.e. a few years. The intermittent transition close to this orbit is, however, erratic and probably hard to predict. This new route to long-term predictability contrasts with the one found in the closed ocean-basin low-order version of MAOOAM, in which the chaotic solution is permanently wandering in the vicinity of an unstable periodic orbit for specific values of the friction coefficient. The model solution is thus at any time influenced by the unstable periodic orbit and inherits from its long-term predictability.
Richard Janis Goldschmidt, Arkady Pikovsky, Antonio Politi
2019-07-14
In globally coupled ensembles of identical oscillators so-called chimera states can be observed. The chimera state is a symmetry-broken regime, where a subset of oscillators forms a cluster, a synchronized population, while the rest of the system remains a collection of non-synchronized, scattered units. We describe here a blinking chimera regime in an ensemble of seven globally coupled rotators (Kuramoto oscillators with inertia). It is characterized by a death-birth process, where a long-term stable cluster of four oscillators suddenly dissolves and is very quickly reborn with a new, reshuffled configuration. We identify three different kinds of rare blinking events and give a quantitative characterization by applying stability analysis to the long-lived chaotic state and to the short-lived regular regimes which arise when the cluster dissolves.
Markus Redeker
2019-07-13
This paper contains two methods to construct one-dimensional number-conserving cellular automata in terms of particle flows. One method is a sequence of increasingly stronger restrictions on the particle flow, which always ends with the specification of a number-conserving rule. The other is based on minimal flow functions, from which all others can be constructed. These constructions also provide a classification for number-conserving rules and a way to specify rules as a supremum of minimal flows. Other questions, like that about the nature of non-deterministic number-conserving rules, are treated briefly at the end.
Xinyu Du, Chao Wang, Tianze Wang, Zeyu Gao
2019-07-13
In this paper, we study reversibility of one-dimensional(1D) linear cellular automata(LCA) under null boundary condition, whose core problems have been divided into two main parts: calculating the period of reversibility and verifying the reversibility in a period. With existing methods, the time and space complexity of these two parts are still too expensive to be employed. So the process soon becomes totally incalculable with a slightly big size, which greatly limits its application. In this paper, we set out to solve these two problems using two efficient algorithms, which make it possible to solve reversible LCA of very large size. Furthermore, we provide an interesting perspective to conversely generate 1D LCA from a given period of reversibility. Due to our methods' efficiency, we can calculate the reversible LCA with large size, which has much potential to enhance security in cryptography system.
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.
@complexcity 