Welcome To Nonlinear Sciences

Adaptation And Self-Organizing Systems

---

Bus bunching as a synchronisation phenomenon (1812.00609v1)

Vee-Liem Saw, Ning Ning Chung, Wei Liang Quek, Yi En Ian Pang, Lock Yue Chew

2018-12-03

Bus bunching is a perennial phenomenon that not only diminishes the efficiency of a bus system, but also prevents transit authorities from keeping buses on schedule. We present a physical theory of buses serving a loop of bus stops as a ring of coupled self-oscillators, analogous to the Kuramoto model. Sustained bunching is a repercussion of the process of phase synchronisation whereby the phases of the oscillators are locked to each other. This emerges when demand exceeds a critical threshold. Buses also bunch at low demand, albeit temporarily, due to frequency detuning arising from different human drivers' distinct natural speeds. We calculate the critical transition when \emph{complete phase locking} (full synchronisation) occurs for the bus system, and posit the critical transition to \emph{completely no phase locking} (zero synchronisation). The intermediate regime comprising \emph{chimera states} is the phase where clusters of partially phase locked buses exist. Intriguingly, these theoretical results are in close correspondence to real buses in a university's shuttle bus system.

Surfactant and gravity dependent instability of two-layer channel flows: Linear theory covering all wave lengths (1801.09290v5)

Alexander F. Frenkel, David Halpern, Adam J. Schweiger

2018-01-28

A linear stability analysis of a two-layer plane Couette flow of two immiscible fluid layers with different densities, viscosities and thicknesses, bounded by two infinite parallel plates moving at a constant relative velocity to each other, with an insoluble surfactant along the interface and in the presence of gravity is carried out. The normal modes approach is applied to the equations governing flow disturbances. These equations, together with boundary conditions at the plates and the interface, yield a linear eigenvalue problem. When inertia is neglected velocity amplitudes are linear combinations of hyperbolic functions, and a quadratic dispersion equation for the complex growth rate is obtained where coefficients depend on the aspect ratio, the viscosity ratio, the basic velocity shear, the Marangoni number Ma that measures the effects of surfactant, and the Bond number Bo that measures the influence of gravity. An extensive investigation is carried out that examines the stabilizing or destabilizing influences of these parameters. There are two continuous branches of the normal modes: a robust branch that exists even with no surfactant, and a surfactant branch that vanishes when Ma . Due to the availability of the explicit forms for the growth rates, in many instances the numerical results are corroborated with analytical asymptotics. For the less unstable branch, a mid-wave interval of unstable wavenumbers (Halpern and Frenkel (2003)) sometimes co-exists with a long-wave one. We study the instability landscape, determined by the threshold curve of the long-wave instability and the critical curve of the mid-wave instability in the (Ma, Bo)-plane. The changes of the extremal points of the critical curves with the variation of the other parameters, such as the viscosity ratio, and the extrema bifurcation points are investigated.

The Inherent Instability of Disordered Systems (1812.00450v1)

Taeer Bar-Yam, Owen Lynch, Yaneer Bar-Yam

2018-12-02

The Multiscale Law of Requisite Variety is a scientific law relating, at each scale, the variation in an environment to the variation in internal state that is necessary for effective response by a system. While this law has been used to describe the effectiveness of systems in self-regulation, the consequences for failure have not been formalized. Here we use this law to consider the internal dynamics of an unstructured system, and its response to a structured environment. We find that, due to its inability to respond, a completely unstructured system is inherently unstable to the formation of structure. And in general, any system without structure above a certain scale is unable to withstand structure arising above that scale. To describe complicated internal dynamics, we develop a characterization of multiscale changes in a system. This characterization is motivated by Shannon information theoretic ideas of noise, but considers structured information. We then relate our findings to political anarchism showing that society requires some organizing processes, even if there is no traditional government or hierarchies. We also formulate our findings as an inverse second law of thermodynamics; while closed systems collapse into disorder, systems open to a structured environment spontaneously generate order.

Connecting empirical phenomena and theoretical models of biological coordination across scales (1812.00423v1)

Mengsen Zhang, Christopher Beetle, J. A. Scott Kelso, Emmanuelle Tognoli

2018-12-02

Coordination is ubiquitous in living systems. Existing theoretical models of coordination -- from bacteria to brains -- focus on either gross statistics in large-scale systems () or detailed dynamics in small-scale systems (mostly ). Both approaches have proceeded largely independent of each other. The present work bridges this gap with a theoretical model of biological coordination that captures key experimental observations of mid-scale social coordination at multiple levels of description. It also reconciles in a single formulation two well-studied models of large- and small-scale biological coordination (Kuramoto and extended Haken-Kelso-Bunz). The model adds second-order coupling (from extended Haken-Kelso-Bunz) to the Kuramoto model. We show that second-order coupling is indispensable for reproducing empirically observed phenomena and gives rise to a phase transition from mono- to multi-stable coordination across scales. This mono-to-multistable transition connects the emergence and growth of behavioral complexity in small and large systems.

The continuum limit of the Kuramoto model on sparse random graphs (1802.03787v3)

Georgi S. Medvedev

2018-02-11

In this paper, we study convergence of coupled dynamical systems on convergent sequences of graphs to a continuum limit. We show that the solutions of the initial value problem for the dynamical system on a convergent graph sequence tend to that for the nonlocal diffusion equation on a unit interval, as the graph size tends to infinity. We improve our earlier results in [Arch. Ration. Mech. Anal., 21 (2014), pp. 781--803] and extend them to a larger class of graphs, which includes directed and undirected, sparse and dense, random and deterministic graphs. There are three main ingredients of our approach. First, we employ a flexible framework for incorporating random graphs into the models of interacting dynamical systems, which fits seamlessly with the derivation of the continuum limit. Next, we prove the averaging principle for approximating a dynamical system on a random graph by its deterministic (averaged) counterpart. The proof covers systems on sparse graphs and yields almost sure convergence on time intervals of order where is the number of vertices. Finally, a Galerkin scheme is developed to show convergence of the averaged model to the continuum limit. The analysis of this paper covers the Kuramoto model of coupled phase oscillators on a variety of graphs including sparse Erd\H{o}s-R{' e}nyi, small-world, and power law graphs.

Chaotic Dynamics

---

Chaos signatures in the short and long time behavior of the out-of-time ordered correlator (1806.04281v3)

Ignacio García-Mata, Marcos Saraceno, Rodolfo A. Jalabert, Augusto J. Roncaglia, Diego A. Wisniacki

2018-06-12

Two properties are needed for a classical system to be chaotic: exponential stretching and mixing. Recently, out-of-time order correlators were proposed as a measure of chaos in a wide range of physical systems. While most of the attention has previously been devoted to the short time stretching aspect of chaos, characterized by the Lyapunov exponent, we show for quantum maps that the out-of-time correlator approaches its stationary value exponentially with a rate determined by the Ruelle-Pollicot resonances. This property constitutes clear evidence of the dual role of the underlying classical chaos dictating the behavior of the correlator at different timescales.

Inhomogeneities and caustics in the sedimentation of noninertial particles in incompressible flows (1801.10116v2)

Gábor Drótos, Pedro Monroy, Emilio Hernández-García, Cristóbal López

2018-01-30

In an incompressible flow, fluid density remains invariant along fluid element trajectories. This implies that the spatial distribution of non-interacting noninertial particles in such flows cannot develop density inhomogeneities beyond those that are already introduced in the initial condition. However, in certain practical situations, density is measured or accumulated on (hyper-) surfaces of dimensionality lower than the full dimensionality of the flow in which the particles move. An example is the observation of particle distributions sedimented on the floor of the ocean. In such cases, even if the initial distribution of noninertial particles is uniform within a finite support in an incompressible flow, advection in the flow will give rise to inhomogeneities in the observed density. In this paper we analytically derive, in the framework of an initially homogeneous particle sheet sedimenting towards a bottom surface, the relationship between the geometry of the flow and the emerging distribution. From a physical point of view, we identify the two processes that generate inhomogeneities to be the stretching within the sheet, and the projection of the deformed sheet onto the target surface. We point out that an extreme form of inhomogeneity, caustics, can develop for sheets. We exemplify our geometrical results with simulations of particle advection in a simple kinematic flow, study the dependence on various parameters involved, and illustrate that the basic mechanisms work similarly if the initial (homogeneous) distribution occupies a more general region of finite extension rather than a sheet.

Memory effects in the Fermi-Pasta-Ulam Model (1806.06210v3)

Graziano Amati, Hugues Meyer, Tanja Schilling

2018-06-16

We study the Intermediate Scattering Function (ISF) of the strongly-nonlinear Fermi-Pasta Ulam Model at thermal equilibrium, using both numerical and analytical methods. From the molecular dynamics simulations we distinguish two limit regimes, as the system behaves as an ideal gas at high temperature and as a harmonic chain for low excitations. At intermediate temperatures the ISF relaxes to equilibrium in a nontrivial fashion. We then calculate analytically the Taylor coefficients of the ISF to arbitrarily high orders (the specific, simple shape of the two-body interaction allows us to derive an iterative scheme for these.) The results of the recursion are in good agreement with the numerical ones. Via an estimate of the complete series expansion of the scattering function, we can reconstruct within a certain temperature range its coarse-grained dynamics. This is governed by a memory-dependent Generalized Langevin Equation (GLE), which can be derived via projection operator techniques. Moreover, by analyzing the first series coefficients of the ISF, we can extract a parameter associated to the strength of the memory effects in the dynamics.

Integrability and Chaos in Figure Skating (1809.11146v3)

Vaughn Gzenda, Vakhtang Putkaradze

2018-09-28

We derive and analyze a three dimensional model of a figure skater. We model the skater as a three-dimensional body moving in space subject to a non-holonomic constraint enforcing movement along the skate's direction and holonomic constraints of continuous contact with ice and pitch constancy of the skate. For a static (non-articulated) skater, we show that the system is integrable if and only if the projection of the center of mass on skate's direction coincides with the contact point with ice and some mild (and realistic) assumptions on the directions of inertia's axes. The integrability is proved by showing the existence of two new constants of motion linear in momenta, providing a new and highly nontrivial example of an integrable non-holonomic mechanical system. We also consider the case when the projection of the center of mass on skate's direction does not coincide with the contact point and show that this non-integrable case exhibits apparent chaotic behavior, by studying the divergence of nearby trajectories We also demonstrate the intricate behavior during the transition from the integrable to chaotic case. Our model shows many features of real-life skating, especially figure skating, and we conjecture that real-life skaters may intuitively use the discovered mechanical properties of the system for the control of the performance on ice.

Periodicity of quantum correlations in the quantum kicked top (1806.06184v2)

Udaysinh T. Bhosale, M. S. Santhanam

2018-06-16

Quantum kicked top is a fundamental model for time-dependent, chaotic Hamiltonian system and has been realized in experiments as well. As the quantum kicked top can be represented as a system of qubits, it is also popular as a testbed for the study of measures of quantum correlations such as entanglement, quantum discord and other multipartite entanglement measures. Further, earlier studies on kicked top have led to a broad understanding of how these measures are affected by the classical dynamical features. In this work, relying on the invariance of quantum correlation measures under local unitary transformations, it is shown exactly these measures display periodic behaviour either as a function of time or as a function of the chaos parameter in this system. As the kicked top has been experimentally realised using cold atoms as well as superconducting qubits, it is pointed out that these periodicities must be factored in while choosing of experimental parameters so that repetitions can be avoided.

Cellular Automata And Lattice Gases

---

The cellular automaton pulsing model, experiments with DDLab (1811.12307v1)

Andrew Wuensche, Edward Coxon

2018-11-29

The cellular automaton (CA) pulsing model (arXiv:1806.06416) described the surprising phenomenon of spontaneous, sustained and robust rhythmic oscillations, pulsing dynamics, when random wiring is applied to a 2D `glider' rule running in a 3-value totalistic CA. Case studies, pulsing measures, possible mechanisms, and implications for oscillatory networks in biology were presented. In this paper we summarise the results, extend the entropy-density and density-return map plots to include a linked history, look at totalistic glider rules with neighborhoods of 3, 4 and 5, as well as 6 and 7 studied previously, introduce methods to automatically recognise the wavelength, and extend results for randomly asynchronous updating. We show how the model is implemented in DDLab to validate results, output data, and allow experiments and research by others.

Selective chaos of travelling waves in feedforward chains of bistable maps (1811.08310v1)

Bastien Fernandez

2018-11-20

We study the chaos of travelling waves (TW) in unidirectional chains of bistable maps. Previous numerical results suggested that this property is selective, {\sl viz.}\ given the parameters, there is at most a single (non-trivial) velocity for which the corresponding set of wave profiles has positive topological entropy. However, mathematical proofs have remained elusive, in particular because the related symbolic dynamics involves entire past sequences. Here, we consider instead inite (short) rank approximations for which the symbolic dynamics has finite memory. For every possible velocity, we compute the existence domains of all possible finite type subshifts of TW with positive entropy. In all examples, chaos of TW turns out to be selective, indeed.

2D Hexagonal Cellular Automata: The Complexity of the Forms (1811.12387v1)

Vural Erdogan

2018-11-09

We created two dimensional hexagonal cellular automata to obtain complexity. Considering the game of life rules, Wolfram's works about life-like structures and John von Neumann's self-replication, self-maintenance, self-reproduction problems, we developed 2-states and 3-states hexagonal growing algorithms that reach large populations through random initial states. Unlike the game of life, we used six neighbourhoods cellular automata instead of eight or four neighbourhoods. First simulations explained that whether we are able to obtain sort of oscillators, blinkers, and gliders. Inspired by Wolfram's 1D cellular automata complexity and life-like structures, we simulated 2D synchronous, discrete, deterministic cellular automata to reach life-like forms with 2-states cells. The life-like formations and the oscillators have been explained how they contribute to initiating self-maintenance together with self-reproduction and self-replication. After comparing the simulation results, we decided to develop the algorithm for another step. Appending a new state to the same algorithm, that we used for reaching life-like structures, led us to experiment new branching and fractal forms. All these studies tried to demonstrate that complex life forms might come from uncomplicated rules.

Theoretical investigation, simulation and empirical analysis of the growth pattern of traffic oscillations in the Euler coordinates (1806.04287v2)

Junfang Tian, Martin Treiber, Rui Jiang, Xiaopeng Li, Bin Jia, Ziyou Gao

2018-06-12

The formation and development of oscillations is an important traffic flow phenomenon. Recent studies found that along a vehicle platoon described in the Lagrangian specification, traffic oscillations grow in a concave way. Since stationary bottlenecks are more intuitively described in the Eulerian framework, this paper investigates whether the concave growth pattern of traffic oscillations in the Lagrangian coordinates can be transferred to the Euler coordinates (i.e. the concave increase in standard deviation is no longer measured across the vehicle indices but as a function of the road location). To this end, we theoretically unify these two ways of measuring oscillations by revealing their mapping relationship. We show that the growth pattern measured in the Lagrangian coordinates can be transferred to the Euler coordinates. We believe this finding is nontrivial since the scenarios are significantly different: while in vehicle platoons (Lagrangian view, non-penetrable moving bottleneck), the speed variance for a given vehicle is ideally constant, the drivers in the Eulerian setting (penetrable stationary bottleneck triggering the waves) experience all amplitudes, first the big ones and then the small ones. To test this proposition, we performed simulation using two different kinds of car-following models. Simulation results validate the theoretical analysis. Finally, we performed empirical analysis using the NGSIM data, which also validates the theoretical analysis.

One-sided scaling limit of multicolor box-ball system (1808.08074v3)

Atsuo Kuniba, Hanbaek Lyu

2018-08-24

A random box-ball system starts with occupying each of the first boxes independently with a ball of random color from , where balls of color 0 are considered as empty boxes. The time evolution is defined by a successive application of the combinatorial , and possesses a -tuple of Young diagrams as the complete set of conserved quantities. We show that the induced carrier process is an irreducible Markov chain and that the associated energy function scaled by converges as at an exponential rate. This implies that each of the invariant Young diagrams converges to some limiting shape as if we scale its rows by . By establishing a large deviations principle, we show that this limiting shape is unchanged under conditioning the initial measure on the highest states. Using these results, we give an explicit formula for the limiting shape of the invariant Young diagrams by Thermodynamic Bethe Ansatz.

Pattern Formation And Solitons

---

Hydrodynamically-bound states of a pair of microrollers: a dynamical system insight (1811.05168v4)

Blaise Delmotte

2018-11-13

Recent work has identified persistent cluster states which were shown to be assembled and held together by hydrodynamic interactions alone [Driscoll \textit{et al.} (2017) Nature Physics, 13(4), 375]. These states were seen in systems of colloidal microrollers; microrollers are colloidal particles which rotate about an axis parallel to the floor and generate strong, slowly decaying, advective flows. To understand these bound states, we study a simple, yet rich, model system of two microrollers. Here we show that pairs of microrollers can exhibit hydrodynamic bound states whose nature depends on a dimensionless number, denoted , that compares the relative strength of gravitational forces and external torques. Using a dynamical system framework, we characterize these various states in phase space and analyze the bifurcations of the system as varies. In particular, we show that there is a critical value, , above which active flows can beat gravity and lead to stable motile orbiting, or "leapfrog", trajectories, reminiscent of the self-assembled motile structures, called "critters", observed by Driscoll \textit{et al}. We identify the conditions for the emergence of these trajectories and study their basin of attraction. This work shows that a wide variety of stable bound states can be obtained with only two particles. Our results aid in understanding the mechanisms that lead to spontaneous self-assembly in hydrodynamic systems, such as microroller suspensions, as well as how to optimize these systems for particle transport.

Pattern Universes (1811.10546v1)

Alan C. Newell, Shankar C. Venkataramani

2018-11-26

Our aim in this essay is to introduce a new paradigm for thinking about the evolution of the universe. Building on previous work on patterns and pattern formation, we suggest a new theoretical construct modeling aspects of the universe on a wide range of scales, from the very small (Planck scales) to the very large (galactic scales). Our "toy models" display parallels to many of the challenges to and mysteries of current thinking in cosmology. We certainly do not claim that what we suggest resolves all the challenges to current theories. Rather, it presents an alternative viewpoint with which to approach these challenges. While it may or may not have anything to do with the real universe, it nevertheless has several intriguing parallels.

Dynamics of Nonlinear Random Walks on Complex Networks (1807.09514v3)

Per Sebastian Skardal, Sabina Adhikari

2018-07-25

In this paper we study the dynamics of nonlinear random walks. While typical random walks on networks consist of standard Markov chains whose static transition probabilities dictate the flow of random walkers through the network, nonlinear random walks consist of nonlinear Markov chains whose transition probabilities change in time depending on the current state of the system. This framework allows us to model more complex flows through networks that may depend on the current system state. For instance, under humanitarian or capitalistic direction, resource flow between institutions may be diverted preferentially to poorer or wealthier institutions, respectively. Importantly, the nonlinearity in this framework gives rise to richer dynamical behavior than occurs in linear random walks. Here we study these dynamics that arise in weakly and strongly nonlinear regimes in a family of nonlinear random walks where random walkers are biased either towards (positive bias) or away from (negative bias) nodes that currently have more random walkers. In the weakly nonlinear regime we prove the existence and uniqueness of a stable stationary state fixed point provided that the network structure is primitive that is analogous to the stationary distribution of a typical (linear) random walk. We also present an asymptotic analysis that allows us to approximate the stationary state fixed point in the weakly nonlinear regime. We then turn our attention to the strongly nonlinear regime. For negative bias we characterize a period-doubling bifurcation where the stationary state fixed point loses stability and gives rise to a periodic orbit below a critical value. For positive bias we investigate the emergence of multistability of several stable stationary state fixed points.

Solitary Wave Propagation in Elastic Bars with Multiple Sections and Layers (1811.09887v1)

M. R. Tranter

2018-11-24

In this paper we present a numerical scheme for solving a system of Boussinesq-type equations. This can correspond to longitudinal displacements in a multi-layered elastic bar with delamination, with conditions on the interface between the sections of the bar. The method is initially presented for two coupled equations in each section and multiple sections in the bar, and later extended to any number of layers. Previous works have presented a similar method constructed using finite-difference methods, however these only solved for two sections of the bar at a time which limited the scope of studies using these methods. The new method presented here solves for all sections at a given time step and therefore the transmitted and reflected waves in each section of the bar can be studied. The new results are shown to be in excellent agreement with previously obtained results and a further study is performed showing that the delamination width can be inferred from the changes to the incident soliton. The generalised form of this method, for any number of sections and layers with coupling terms independent of time derivatives, can be used to study the behaviour of longitudinal waves in more complicated waveguides in future studies.

Excitation of Peregrine-type waveforms from vanishing initial conditions in the presence of periodic forcing (1811.09812v1)

Nikos I. Karachalios, Paris Kyriazopoulos, Konstantinos Vetas

2018-11-24

We show by direct numerical simulations that spatiotemporally localized wave forms, strongly reminiscent of the Peregrine rogue wave, can be excited by vanishing initial conditions for the periodically driven nonlinear Schr"odinger equation. The emergence of the Peregrine-type waveforms can be potentially justified, in terms of the existence and modulational instability of spatially homogeneous solutions of the model, and the continuous dependence of the localized initial data for small time intervals. We also comment on the persistence of the above dynamics, under the presence of small damping effects, and justify, that this behavior should be considered as far from approximations of the corresponding integrable limit.

Exactly Solvable And Integrable Systems

---

The geometry of nonholonomic Chaplygin systems revisited (1812.01422v1)

Luis C. García-Naranjo, Juan C. Marrero

2018-12-04

We consider nonholonomic Chaplygin systems and associate to them a tensor field on the shape space, that we term the gyroscopic tensor, and that measures the interplay between the constraint distribution and the kinetic energy metric. We show how this tensor may be naturally used to derive an almost symplectic description of the reduced dynamics. Moreover, we express sufficient conditions for measure preservation and Hamiltonisation via Chaplygin's reducing multiplier method in terms of the properties of this tensor. The theory is used to give a new proof of the remarkable Hamiltonisation of the multi-dimensional Veselova system obtained by Fedorov and Jovanovi'c in J. Nonlinear Sci. 2004, and Reg. Chaot. Dyn. 2009.

-Rosochatius system, superintegrability, supersymmetry (1812.00930v1)

Evgeny Ivanov, Armen Nersessian, Hovhannes Shmavonyan

2018-12-03

We propose new superintegrable mechanical system on the complex projective space involving a potential term together with coupling to a constant magnetic fields. This system can be viewed as a -analog of both the flat singular oscillator and its spherical analog known as "Rosochatius system". We find its constants of motion and calculate their (highly nonlinear) algebra. We also present its classical and quantum solutions. The system belongs to the class of "K"ahler oscillators" admitting supersymmetric extension. We show that, in the absence of magnetic field and with the special choice of the characteristic parameters, one can construct Poinacar'e supersymmetric extension of the system considered.

Domain-wall boundaries through non-diagonal twists in the six-vertex model (1807.03659v3)

W. Galleas

2018-07-10

In this work we elaborate on a previous result relating the partition function of the six-vertex model with domain-wall boundary conditions to eigenvalues of a transfer matrix. More precisely, we express the aforementioned partition function as a determinant of a matrix with entries being eigenvalues of the anti-periodic six-vertex model's transfer matrix.

Restoring discrete Painlevé equations from an E-associated one (1812.00712v1)

Basil Grammaticos, Alfred Ramani, Ralph Willox

2018-12-03

We present a systematic method for the construction of discrete Painlev'e equations. The method, dubbed `restoration', allows one to obtain all discrete Painlev'e equations that share a common autonomous limit, up to homographic transformations, starting from any one of those limits. As the restoration process crucially depends on the classification of canonical forms for the mappings in the QRT family, it can in principle only be applied to mappings that belong to that family. However, as we show in this paper, it is still possible to obtain the results of the restoration even when the initial mapping is not of QRT type (at least for the system at hand, but we believe our approach to be of much wider applicability). For the equations derived in this paper we also show how, starting from a form where the independent variable advances one step at a time, one can obtain versions corresponding to multistep evolutions.

Projective loops generate rational loop groups (1812.01456v1)

Gang Wang, Oliver Goertsches, Erxiao Wang

2018-12-03

Rational loops played a central role in Uhlenbeck's construction of harmonic maps into U(n) (chiral model in physics), and they are generated by simple elements with one pole and one zero constructed from Hermitian projections. It has been believed for long time that nilpotent loops should be added to generate rational loop groups with noncompact reality conditions. We prove a somewhat unexpected theorem that projective loops are enough to generate the rational loop groups of GL(n,C), GL(n,R), and U(p, q).

---

Thank you for reading!

Don't forget to Follow and Resteem. @complexcity
Keeping everyone inform.