<center><b> Welcome To Nonlinear Sciences</b></center>

# <center>Adaptation And Self-Organizing Systems</center> 
 <hr> 

### [Bus bunching as a synchronisation phenomenon](http://arxiv.org/abs/1812.00609v2) (1812.00609v2)
<i>Vee-Liem Saw, Ning Ning Chung, Wei Liang Quek, Yi En Ian Pang, Lock Yue Chew</b>

<h10>2018-12-03</h10>
> 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 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.

### [Directed Flow of Information in Chimera States](http://arxiv.org/abs/1904.11553v1) (1904.11553v1)
<i>Nicolás Deschle, Andreas Daffertshofer, Demian Battaglia, Erik A. Martens</b>

<h10>2019-04-25</h10>
> We investigated interactions within chimera states in a phase oscillator network with two coupled subpopulations. To quantify interactions within and between these subpopulations, we estimated the corresponding (delayed) mutual information that -- in general -- quantifies the capacity or the maximum rate at which information can be transferred to recover a sender's information at the receiver with a vanishingly low error probability. After verifying their equivalence with estimates based on the continuous phase data, we determined the mutual information using the time points at which the individual phases passed through their respective Poincar\'{e} sections. This stroboscopic view on the dynamics may resemble, e.g., neural spike times, that are common observables in the study of neuronal information transfer. This discretization also increased processing speed significantly, rendering it particularly suitable for a fine-grained analysis of the effects of experimental and model parameters. In our model, the delayed mutual information within each subpopulation peaked at zero delay, whereas between the subpopulations it was always maximal at non-zero delay, irrespective of parameter choices. We observed that the delayed mutual information of the desynchronized subpopulation preceded the synchronized subpopulation. Put differently, the oscillators of the desynchronized subpopulation were 'driving' the ones in the synchronized subpopulation. These findings were also observed when estimating mutual information of the full phase trajectories. We can thus conclude that the delayed mutual information of discrete time points allows for inferring a functional directed flow of information between subpopulations of coupled phase oscillators.

### [Gaussian statistics as an emergent symmetry of the stochastic Burgers  equation](http://arxiv.org/abs/1809.02158v2) (1809.02158v2)
<i>Enrique Rodriguez-Fernandez, Rodolfo Cuerno</b>

<h10>2018-09-06</h10>
> Symmetries play a conspicuous role in the large-scale behavior of critical systems. While in equilibrium they allow to classify asymptotics into different universality classes, out of equilibrium they can emerge, some times unexpectedly, as collective properties which are not explicit in the 'bare' interactions. Here we elucidate the emergence of an up-down symmetry in the asymptotic behavior of the stochastic scalar Burgers equation in one and two dimensions, manifested by the occurrence of Gaussian fluctuations for the physical field. This robustness of Gaussian behavior contradicts naive expectations, due to the detailed relation ---including the same set of symmetries--- between Burgers equation and the Kardar-Parisi-Zhang equation, which paradigmatically displays non-Gaussian fluctuations described by Tracy-Widom distributions. We reach our conclusions via a dynamic renormalization group study of the field statistics, confirmed by direct evaluation of the field probability distribution function from numerical simulations of the dynamical equation. All odd-order cumulants cancel exactly, while the excess kurtosis vanishes for large systems, indeed consistent with Gaussian behavior.

### [Concurrence of form and function in developing networks and its role in  synaptic pruning](http://arxiv.org/abs/1705.02773v2) (1705.02773v2)
<i>Ana P. Millán, J. J. Torres, S. Johnson, J. Marro</b>

<h10>2017-05-08</h10>
> A fundamental question in neuroscience is how structure and function of neural systems are related. We study this interplay by combining a familiar auto-associative neural network with an evolving mechanism for the birth and death of synapses. A feedback loop then arises leading to two qualitatively different types of behaviour. In one, the network structure becomes heterogeneous and dissasortative, and the system displays good memory performance; furthermore, the structure is optimised for the particular memory patterns stored during the process. In the other, the structure remains homogeneous and incapable of pattern retrieval. These findings provide an inspiring picture of brain structure and dynamics, are compatible with experimental results on early brain development, and may help to explain synaptic pruning. Other evolving networks -- such as those of protein interaction -- might share the basic ingredients for this feedback loop and other questions, and indeed many of their structural features are as predicted by our model.

### [Marginalized ordering and adaptive reaction time in bird flocks](http://arxiv.org/abs/1902.09288v2) (1902.09288v2)
<i>Narina Jung, Pilwon Kim</b>

<h10>2019-02-25</h10>
> We develop a model for a rich dynamics of a flock in a marginalized ordering state. The aim is to present an inter-individual coordination mechanism that keeps a flock constantly ready to respond to perturbations naturally present in biological systems. We extend the generalized Cucker-Smale model with the coupling of acceleration and introduce adaptive reaction times of each bird. We regard two key factors in the reaction times: (1) the local ordering state of each bird and (2) reaction sensitivity of a flock to the neighbor's momentum change with ![](http://latex2png.com/output//latex_e40752af2e0771d9e3b4727621bf6015.png). We show that our model displays innate fluctuations and rich dynamics as a reminiscent of natural flocks due to the adaptive reaction delay. This happens without relying on stochastic variables. We compute the correlation lengths of the fluctuations and find that the correlation of velocity and speed is scale-free, indicating some criticality of a flock. It is dynamically in a marginalized ordering states, rather than in either an ordered or a disordered state. Surprisingly, at a large value of ![](http://latex2png.com/output//latex_e40752af2e0771d9e3b4727621bf6015.png) (i.e., reaction sensitivity is high), the transition occurs from the standard diffusion to the super-diffusive Levy flights by increasing the strength of the velocity alignment. Our results indicate that the emergence of the long-term behaviors such as Levy flights can also be explained in terms of the inter-individual interaction that makes the system in a marginalized ordering state.

# <center>Chaotic Dynamics</center> 
 <hr> 

### [Role of acceleration in inducing chaotic fluctuations in particle  dynamics](http://arxiv.org/abs/1904.11760v1) (1904.11760v1)
<i>Surojit Dalui, Bibhas Ranjan Majhi, Pankaj Mishra</b>

<h10>2019-04-26</h10>
> The ongoing conjecture that the presence of horizon may induce chaos in an integrable system, is further investigated from the perspective of a uniformly accelerated frame. Particularly, we build up a model which consists of a particle trapped in harmonic oscillator and analyzed from uniformly accelerated frame (namely Rindler observer). Here the Rindler frame provides a Killing horizon without any intrinsic curvature to the system. This makes the present observations different from previous studies. We observe that for some particular values of parameters of the system (like acceleration, energy of the particle), the motion of the particle trapped in harmonic potential systematically goes from periodic state to the chaotic. This indicates that the existence of horizon alone, not the intrinsic curvature (i.e. the gravitational effect) in the background, is sufficient to induce the chaotic motion in the particle. We believe the present study further enlighten and balustrade the conjecture.

### [Arnold maps with noise: Differentiability and non-monotonicity of the  rotation number](http://arxiv.org/abs/1904.11744v1) (1904.11744v1)
<i>L. Marangio, J. Sedro, S. Galatolo, A. Di Garbo, M. Ghil</b>

<h10>2019-04-26</h10>
> Arnold's standard circle maps are widely used to study the quasi-periodic route to chaos and other phenomena associated with nonlinear dynamics in the presence of two rationally unrelated periodicities. In particular, the El Nino-Southern Oscillation (ENSO) phenomenon is a crucial component of climate variability on interannual time scales and it is dominated by the seasonal cycle, on the one hand, and an intrinsic oscillatory instability with a period of a few years, on the other. The role of meteorological phenomena on much shorter time scales, such as westerly wind bursts, has also been recognized and modeled as additive noise. We consider herein Arnold maps with additive, uniformly distributed noise. When the map's nonlinear term, scaled by the parameter ![](http://latex2png.com/output//latex_08ef962254d04558fe63cdc5bc67c170.png), is sufficiently small, i.e. {![](http://latex2png.com/output//latex_91cd9b4161ed82febb02ffeb26b1f34f.png)}, the map is known to be a diffeomorphism and the rotation number ![](http://latex2png.com/output//latex_0e1dd1a2da68033e2646e3d1cee354f3.png) is a differentiable function of the driving frequency ![](http://latex2png.com/output//latex_2d32baf7f1565d8a2542a3e420f3eef4.png). We concentrate on the rotation number's behavior as the nonlinearity becomes large, and show rigorously that ![](http://latex2png.com/output//latex_78723c0c29cdee2bf7d70b1175e2955f.png) is a differentiable function of ![](http://latex2png.com/output//latex_2d32baf7f1565d8a2542a3e420f3eef4.png), even for ![](http://latex2png.com/output//latex_1b471d942450ed6044d72967cdcb2b94.png), at every point at which the noise-perturbed map is mixing. We also provide a formula for the derivative of the rotation number. The reasoning relies on linear-response theory and a computer-aided proof. In the diffeomorphism case of {\ ![](http://latex2png.com/output//latex_4ba7048e8b727c284535a388f6f6ce8f.png),} the rotation number ![](http://latex2png.com/output//latex_78723c0c29cdee2bf7d70b1175e2955f.png) behaves monotonically with respect to ![](http://latex2png.com/output//latex_2d32baf7f1565d8a2542a3e420f3eef4.png). We show, using again a computer-aided proof, that this is not the case when ![](http://latex2png.com/output//latex_1b471d942450ed6044d72967cdcb2b94.png) and the map is not a diffeomorphism.

### [Probing Quantum Chaos in many-body quantum systems by the induced  dissipation](http://arxiv.org/abs/1904.11607v1) (1904.11607v1)
<i>Anna A . Bychek, Pavel S. Muraev, Andrey R. Kolovsky</b>

<h10>2019-04-25</h10>
> We theoretically analyze the depletion dynamics of an ensemble of cold atoms in a quasi one-dimensional optical lattice where atoms in one of the lattice sites are subject to decay. Unlike the previous studies of this problem in R. Labouvie, {\em et. al}, Phys. Rev. Lett. {\bf 116}, 235302 (2016) we focus on the case where the system is brought to the chaotic regime, that crucially modifies the depletion dynamics as compared to the regular case. It is shown that depletion of the affected site results in gradual depletion of the neighbouring sites according to ![](http://latex2png.com/output//latex_3f4eabb3ed9793fea7630278c801f042.png) scaling law. We also show that by measuring occupations of the lattice sites one can extract important information on chaotic dynamics of the original conservative system.

### [A practical method for estimating coupling functions in complex  dynamical systems](http://arxiv.org/abs/1904.11289v1) (1904.11289v1)
<i>Isao T. Tokuda, Zoran Levnajic, Kazuyoshi Ishimura</b>

<h10>2019-04-25</h10>
> A foremost challenge in modern network science is the inverse problem on reconstruction (inference) of coupling equations and network topology from the measurements of the network's dynamics. Of particular interest are the methods that can operate on real (empirical) data without interfering with the system. One of such earlier attempts [Tokuda et al., PRL, 2007] was a method suited for general limit-cycle oscillators, yielding both oscillators' natural frequencies and coupling functions between them (phase equations) from empirically measured time series. The present paper makes two advances along the same lines. First, it reviews the above method in a way comprehensive to domain-scientists other than physics. Second, it presents three new extensions: (i) algorithm for inference of the phase sensitivity function, (ii) coupling function to approximate interaction among phase-coherent chaotic oscillators, (iii) analysis of the experimental data from a forced Van der Pol electric circuit. This reaffirms the range of applicability of the method and makes it accessible to a wider scientific community.

### [Turbulent drag reduction: A universal perspective from energy fluxes](http://arxiv.org/abs/1904.11167v1) (1904.11167v1)
<i>Mahendra K. Verma</b>

<h10>2019-04-25</h10>
> Injection of dilute polymer in a turbulent flow suppresses frictional drag. This challenging and technologically important problem remains primarily unresolved due to the complex nature of the flow. An important factor in the drag reduction is the energy transfer from the velocity field to the polymers. In this paper we quantify this process using energy fluxes, as well as show its universality in diverse flows such as magnetohydrodynamics, quasi-static magnetohydrodynamics, and bubbly turbulence. We show that in such flows, the transfer from kinetic energy to elastic energy leads to a reduction in kinetic energy flux compared to the corresponding hydrodynamic turbulence. This leads to a reduction in nonlinearity of the velocity field that results into a more ordered flow and a suppression of turbulent drag.

# <center>Cellular Automata And Lattice Gases</center> 
 <hr> 

### [A Novel Antifragility Measure Based on Satisfaction and Its Application  to Random and Biological Boolean Networks](http://arxiv.org/abs/1812.06760v2) (1812.06760v2)
<i>Omar K. Pineda, Hyobin Kim, Carlos Gershenson</b>

<h10>2018-12-17</h10>
> Antifragility is a property that enhances the capability of a system in response to external perturbations. Although the concept has been applied in many areas, a practical measure of antifragility has not been developed yet. Here we propose a simply calculable measure of antifragility, based on the change of "satisfaction" before and after adding perturbations, and apply it to random Boolean networks (RBNs). Using the measure, we found that ordered RBNs are the most antifragile. Also, we demonstrated that seven biological systems are antifragile. Our measure and results can be used in various applications of Boolean networks (BNs) including creating antifragile engineering systems, identifying the genetic mechanism of antifragile biological systems, and developing new treatment strategies for various diseases.

### [Simply modified GKL density classifiers that reach consensus faster](http://arxiv.org/abs/1904.07411v1) (1904.07411v1)
<i>J. Ricardo G. Mendonça</b>

<h10>2019-04-16</h10>
> The two-state Gacs-Kurdyumov-Levin (GKL) cellular automaton has been a staple model in the study of complex systems due to its ability to classify binary arrays of symbols according to their initial density. We show that a class of modified GKL models over extended neighborhoods, but still involving only three cells at a time, achieves comparable density classification performance but in some cases reach consensus more than twice as fast. Our results suggest the time to consensus (relative to the length of the CA) as a complementary measure of density classification performance.

### [A Graph Theory Approach for Regional Controllability of Boolean Cellular  Automata](http://arxiv.org/abs/1904.07735v1) (1904.07735v1)
<i>Sara Dridi, Samira El Yacoubi, Franco Bagnoli, Allyx Fontaine</b>

<h10>2019-04-15</h10>
> Controllability is one of the central concepts of modern control theory that allows a good understanding of a system's behaviour. It consists in constraining a system to reach the desired state from an initial state within a given time interval. When the desired objective affects only a sub-region of the domain, the control is said to be regional. The purpose of this paper is to study a particular case of regional control using cellular automata models since they are spatially extended systems where spatial properties can be easily defined thanks to their intrinsic locality. We investigate the case of boundary controls on the target region using an original approach based on graph theory. Necessary and sufficient conditions are given based on the Hamiltonian Circuit and strongly connected component. The controls are obtained using a preimage approach.

### [Advanced Aspects of the Galactic Habitability](http://arxiv.org/abs/1904.01062v2) (1904.01062v2)
<i>Vladimir Đošović, Branislav Vukotić, Milan M. Ćirković</b>

<h10>2019-04-01</h10>
> Context. Astrobiological evolution of the Milky Way (or the shape of its "astrobiological landscape") has emerged as one of the key research topics in recent years. In order to build precise, quantitative models of the Galactic habitability, we need to account for two opposing tendencies of life and intelligence in the most general context: the tendency to spread to all available ecological niches (conventionally dubbed "colonization") and the tendency to succumb to various types of existential catastrophes ("catastrophism"). These evolutionary tendencies have become objects of study in fields such as ecology, macroevolution, risk analysis, and futures studies, while a serious astrobiological treatment has been lacking so far. Aims. Our aim is to numerically investigate the dynamics of opposed processes of expansion (panspermia, colonization) and extinction (catastrophic mechanisms) of life in the Galaxy. Methods. We employ a new type of numerical simulation based on 1D probabilistic cellular automaton with very high temporal resolution, in order to study astrobiological dynamics. Results. While the largest part of the examined parameter space shows very low habitability values, as expected, the remaining part has some observationally appealing features that imply, among other things, a reduction in the amount of fine-tuning necessary for resolving the Fermi paradox. Conclusions. Advanced aspects of Galactic habitability are amenable to precision studies using massive parallel computer simulations. There are regions of parameter space corresponding to a quasi-stationary state satisfying observable constraints and possessing viable SETI targets.

### [Models of infiltration into homogeneous and fractal porous media with  localized sources](http://arxiv.org/abs/1904.05231v1) (1904.05231v1)
<i>Fabio D. A. Aarao Reis, Vaughan R. Voller</b>

<h10>2019-04-10</h10>
> We study a random walk infiltration (RWI) model, in homogeneous and in fractal media, with localized sources at their boundaries. The particles released at a source, which is maintained at a constant density, execute unbiased random walks over a lattice; it represents solute infiltration by diffusion into a medium in contact with a reservoir. A scaling approach shows that the infiltrated length, area, or volume evolves in time as the number of distinct sites visited by a single random walker in the same medium. This is consistent with simulations of the lattice model and exact and numerical solutions of the corresponding diffusion equation. In a Sierpinski carpet, the infiltrated area is expected to evolve as t^{D_F/D_W} (Alexander-Orbach relation), where D_F is the fractal dimension of the medium and D_W is the random walk dimension; the numerical integration of the diffusion equation supports this relation and improves results of lattice random walk simulations. In a Menger sponge in which D_F>D_W (a fractal with a dimension close to 3), a linear time increase of the infiltrated volume is predicted and confirmed numerically. Thus, no evidence of fractality can be observed in infiltrated volumes or masses in media where random walks are not recurrent, although the tracer diffusion is anomalous. We compare our results with a fluid infiltration model in which the pressure head is constant at the source and the front displacement is driven by the local gradient of that head. Exact or numerical solutions in two and three dimensions and in a carpet show that this type of fluid infiltration is in the same universality class of RWI, with an equivalence between the head and the particle concentration. These results set a relation between different infiltration processes with localized sources and the recurrence properties of random walks in the same media.

 <br><hr> <center>Thank you for reading!<br> https://cdn.steemitimages.com/DQmbn3ovuKLM17k6aemZMrJj6iqKkYzXCYz5Qh1Fg7vPmRx/image.png <br> Don't forget to Follow and Resteem. @complexcity <br>Keeping everyone inform.</center>