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

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

### [Effective networks: a model to predict network structure and critical  transitions from datasets](http://arxiv.org/abs/1907.02416v1) (1907.02416v1)
<i>Deniz Eroglu, Matteo Tanzi, Sebastian van Strien, Tiago Pereira</b>

<h10>2019-07-04</h10>
> Real-world complex systems such as ecological communities and neuron networks are essential parts of our everyday lives. These systems are composed of units which interact through intricate networks. The ability to predict sudden changes in network behaviour, known as critical transitions, from data is important to avert disastrous consequences of major disruptions. Predicting such changes is a major challenge as it requires forecasting the behaviour for parameter ranges for which no data on the system is available. In this paper, we address this issue for networks with weak individual interactions and chaotic local dynamics. We do this by building a model network, termed an effective network, consisting of the underlying local dynamics at each node and a statistical description of their interactions. We illustrate this approach by reconstructing the dynamics and structure of realistic neuronal interaction networks of the cat cerebral cortex. We reconstruct the community structure by analysing the stochastic fluctuations generated by the network and predict critical transitions for coupling parameters outside the observed range.

### [Phase reduction beyond the first order: the case of the mean-field  complex Ginzburg-Landau equation](http://arxiv.org/abs/1907.02276v1) (1907.02276v1)
<i>Iván León, Diego Pazó</b>

<h10>2019-07-04</h10>
> Phase reduction is a powerful technique that permits to describe the dynamics of a weakly perturbed limit-cycle oscillator in terms of its phase. For ensembles of oscillators, a classical example of phase reduction is the derivation of the Kuramoto model from the mean-field complex Ginzburg-Landau equation (MF-CGLE). Still, the Kuramoto model is a first-order phase approximation that displays either full synchronization or incoherence, but none of the nontrivial dynamics of the MF-CGLE. This fact calls for an expansion beyond the first order in the coupling constant. We develop an isochron-based scheme to obtain the second-order phase approximation, which reproduces the weak coupling dynamics of the MF-CGLE. The practicality of our method is evidenced by extending the calculation up to third order. Each new term of the power series expansion contributes with additional higher-order multi-body (i.e.non-pairwise) interactions. This points to intricate multi-body phase interactions as the source of pure collective chaos in the MF-CGLE at moderate coupling.

### [A dynamical systems model of unorganised segregation in two  neighbourhoods](http://arxiv.org/abs/1907.01941v1) (1907.01941v1)
<i>D. J. Haw, S. J. Hogan</b>

<h10>2019-07-03</h10>
> We present a complete analysis of the Schelling dynamical system [Haw2018] of two connected neighbourhoods, with or without population reservoirs, for different types of linear and nonlinear tolerance schedules. We show that stable integration is only possible when the minority is small and combined tolerance is large. Unlike the case of the single neighbourhood, limiting one population does not necessarily produce stable integration and may destroy it. We conclude that a growing minority can only remain integrated if the majority increases its own tolerance. Our results show that an integrated single neighbourhood may not remain so when a connecting neighbourhood is created.

### [Boltzmann distribution of sediment transport](http://arxiv.org/abs/1907.01880v1) (1907.01880v1)
<i>A. Abramian, O. Devauchelle, G. Seizilles, E. Lajeunesse</b>

<h10>2019-07-03</h10>
> The coupling of sediment transport with the flow that drives it allows rivers to shape their own bed. Cross-stream fluxes of sediment play a crucial, yet poorly understood, role in this process. Here, we track particles in a laboratory flume to relate their statistical behavior to the self organization of the granular bed they make up. As they travel downstream, the transported grains wander randomly across the bed's surface, thus inducing cross-stream diffusion. The balance of diffusion and gravity results in a peculiar Boltzmann distribution, in which the bed's roughness plays the role of thermal fluctuations, while its surface forms the potential well that confines the sediment flux.

### [Policies for allocation of information in task-oriented groups: elitism  and egalitarianism outperform welfarism](http://arxiv.org/abs/1902.08183v2) (1902.08183v2)
<i>Sandro M. Reia, Paulo F. Gomes, José F. Fontanari</b>

<h10>2019-02-21</h10>
> Communication or influence networks are probably the most controllable of all factors that are known to impact on the problem-solving capability of task-forces. In the case connections are costly, it is necessary to implement a policy to allocate them to the individuals. Here we use an agent-based model to study how distinct allocation policies affect the performance of a group of agents whose task is to find the global maxima of NK fitness landscapes. Agents cooperate by broadcasting messages informing on their fitness and use this information to imitate the fittest agent in their influence neighborhoods. The larger the influence neighborhood of an agent, the more links, and hence information, the agent receives. We find that the elitist policy in which agents with above-average fitness have their influence neighborhoods amplified, whereas agents with below-average fitness have theirs deflated, is optimal for smooth landscapes, provided the group size is not too small. For rugged landscapes, however, the elitist policy can perform very poorly for certain group sizes. In addition, we find that the egalitarian policy, in which the size of the influence neighborhood is the same for all agents, is optimal for both smooth and rugged landscapes in the case of small groups. The welfarist policy, in which the actions of the elitist policy are reversed, is always suboptimal, i.e., depending on the group size it is outperformed by either the elitist or the egalitarian policies.

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

### [Exact localized and ballistic eigenstates in disordered chaotic spin  ladders and the Fermi-Hubbard model](http://arxiv.org/abs/1811.07903v2) (1811.07903v2)
<i>Thomas Iadecola, Marko Znidaric</b>

<h10>2018-11-19</h10>
> We demonstrate the existence of exact atypical many-body eigenstates in a class of disordered, interacting one-dimensional quantum systems that includes the Fermi-Hubbard model as a special case. These atypical eigenstates, which generically have finite energy density and are exponentially many in number, are populated by noninteracting excitations. They can exhibit Anderson localization with area-law eigenstate entanglement or, surprisingly, ballistic transport at any disorder strength. These properties differ strikingly from those of typical eigenstates nearby in energy, which we show give rise to diffusive transport as expected in a chaotic quantum system. We discuss how to observe these atypical eigenstates in cold-atom experiments realizing the Fermi-Hubbard model, and comment on the robustness of their properties.

### [Phase reduction beyond the first order: the case of the mean-field  complex Ginzburg-Landau equation](http://arxiv.org/abs/1907.02276v1) (1907.02276v1)
<i>Iván León, Diego Pazó</b>

<h10>2019-07-04</h10>
> Phase reduction is a powerful technique that permits to describe the dynamics of a weakly perturbed limit-cycle oscillator in terms of its phase. For ensembles of oscillators, a classical example of phase reduction is the derivation of the Kuramoto model from the mean-field complex Ginzburg-Landau equation (MF-CGLE). Still, the Kuramoto model is a first-order phase approximation that displays either full synchronization or incoherence, but none of the nontrivial dynamics of the MF-CGLE. This fact calls for an expansion beyond the first order in the coupling constant. We develop an isochron-based scheme to obtain the second-order phase approximation, which reproduces the weak coupling dynamics of the MF-CGLE. The practicality of our method is evidenced by extending the calculation up to third order. Each new term of the power series expansion contributes with additional higher-order multi-body (i.e.non-pairwise) interactions. This points to intricate multi-body phase interactions as the source of pure collective chaos in the MF-CGLE at moderate coupling.

### [Multifaceted nonlinear dynamics in ![](http://latex2png.com/output//latex_7b56e43ee24b86f2108ea09eb3096ca6.png)-symmetric coupled  Liénard oscillators](http://arxiv.org/abs/1812.10126v2) (1812.10126v2)
<i>Jyoti Prasad Deka, Arjunan Govindarajan, Kulkarni, Amarendra K. Sarma</b>

<h10>2018-12-25</h10>
> We propose a generalized parity-time (![](http://latex2png.com/output//latex_7b56e43ee24b86f2108ea09eb3096ca6.png)) -symmetric Li\'enard oscillator with two different orders of nonlinear position-dependent dissipation. We study the stability of the stationary states by using the eigenvalues of Jacobian and evaluate the stability threshold thereafter. In the first order nonlinear damping model, we discover that the temporal evolution of both gain and lossy oscillators attains a complete convergence towards the stable stationary state leading to simultaneous emergence of oscillation and amplitude deaths. Also, the system displays a remarkable manifestation of transient chaos in the lossy oscillator while the gain counterpart exhibits blow-up dynamics for certain choice of initial conditions and control parameters. Employing an external driving force on the loss oscillator, we find that the blow-up dynamics can be controlled and a pure aperiodic state is achievable. On the other hand, the second order nonlinear damping model yields a completely different dynamics on contrary to the first order where the former reveals a conventional quasi-periodic route to chaos upon decreasing the natural frequency of both gain and loss oscillators. An electronic circuit scheme for the experimental realization of the proposed system has also been put forward.

### [Extremely High Length-Divergent Thermal Conductivity in Long-Range  Interacting Fermi-Pasta-Ulam Chains](http://arxiv.org/abs/1906.11086v2) (1906.11086v2)
<i>Jianjin Wang, Sergey V. Dmitriev, Daxing Xiong</b>

<h10>2019-06-26</h10>
> The power-law length (![](http://latex2png.com/output//latex_0ee8776a295f6d06670ea84906a0e45a.png)) divergence of thermal conductivity (![](http://latex2png.com/output//latex_39a1acc0a1fe620307fc63cb0815e3a7.png)) in one-dimensional (1D) systems, i.e., ![](http://latex2png.com/output//latex_b20866e36c2c81c0136fd88a226ed4ee.png), has been predicted by theories and also corroborated by experiments. The theoretical predictions of the exponent ![](http://latex2png.com/output//latex_879bee8d9bc4848dee78650baeff05fd.png) are usually ranging from ![](http://latex2png.com/output//latex_90e13cc294761693602fb52e9bfa10a3.png) to ![](http://latex2png.com/output//latex_1dc26b2661fa7239c25eb86fa8572e08.png); however sometimes, the experimental observations can be higher, e.g., ![](http://latex2png.com/output//latex_39328b34a56faaf05d072b0ec10cb72b.png)-![](http://latex2png.com/output//latex_3868dd6a68cc6c62a53407e92b0150c9.png). This dispute has not yet been settled. Here we show the first convincing evidence that an exponent of ![](http://latex2png.com/output//latex_3e09866b5ba9179843d4b951ca50c1f1.png) 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.

### [Sensitivity computation of statistically stationary quantities in  turbulent flows](http://arxiv.org/abs/1905.09362v2) (1905.09362v2)
<i>Nisha Chandramoorthy, Qiqi Wang</b>

<h10>2019-05-22</h10>
> It is well-known that linearized perturbation methods for sensitivity analysis, such as tangent or adjoint equation-based, finite difference and automatic differentiation are not suitable for turbulent flows. The reason is that turbulent flows exhibit chaotic dynamics, leading to the norm of an infinitesimal perturbation to the state growing exponentially in time. As a result, these conventional methods cannot be used to compute the derivatives of long-time averaged quantities to control or design inputs. The ensemble-based approaches and shadowing-based approaches to circumvent the problems of the conventional methods in chaotic systems, also suffer from computational impracticality and lack of consistency guarantees, respectively. We introduce the space-split sensitivity, or the S3 algorithm, which is a Monte-Carlo approach to the chaotic sensitivity computation problem. In this work, we derive the S3 algorithm under simplifying assumptions on the dynamics and present a numerical validation on a low-dimensional example of chaos.

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

### [Kardar-Parisi-Zhang Universality of the Nagel-Schreckenberg Model](http://arxiv.org/abs/1907.00636v1) (1907.00636v1)
<i>Jan de Gier, Andreas Schadschneider, Johannes Schmidt, Gunter M. Schütz</b>

<h10>2019-07-01</h10>
> Dynamical universality classes are distinguished by their dynamical exponent ![](http://latex2png.com/output//latex_bd61b29a408a64d11f1a55bca65fd0c4.png) 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 ![](http://latex2png.com/output//latex_2fdb6000f3667b6e6887c72c14280977.png). 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 ![](http://latex2png.com/output//latex_612eb39efa60bedf469ce3141e6664ff.png). In this paper, we show that the NaSch model also belongs to the KPZ class \cite{KPZ} for general maximum velocities ![](http://latex2png.com/output//latex_d470cde12bfb573e851ddd17b54b2b41.png). 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.

### [Evaluation on asymptotic distribution of particle systems expressed by  probabilistic cellular automata](http://arxiv.org/abs/1907.01635v1) (1907.01635v1)
<i>Kazushige Endo</b>

<h10>2019-06-29</h10>
> 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.

### [Shift-Symmetric Configurations in Two-Dimensional Cellular Automata:  Irreversibility, Insolvability, and Enumeration](http://arxiv.org/abs/1703.09030v2) (1703.09030v2)
<i>Peter Banda, John Caughman, Martin Cenek, Christof Teuscher</b>

<h10>2017-03-27</h10>
> 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.

### [Improving the bus flow in a Bus Rapid Transit system: an approach based  on cellular automata simulations](http://arxiv.org/abs/1906.11916v1) (1906.11916v1)
<i>M. A. Uribe-Laverde, W. F. Oquendo-Patiño</b>

<h10>2019-06-27</h10>
> We studied the bus flow in a Bus Rapid Transit (BRT) system using a novel approach based on a cellular automata (CA) that properly accounts for bus interactions. The model quantitatively reproduces the bus queuing behaviour for both fixed and random dwell times, with one bus service. For more bus services, our results show that the bus flow strongly depends not only on how the buses are distributed among the services, but also on how the bus stops are physically arranged at the stations. The latter dependence has hitherto not been considered in common BRT optimization methods. By using a genetic algorithm, we could find the bus distribution that maximizes the bus flow for a given stops arrangement at the stations, and explore its evolution with the bus density. Our results suggest a set of guidelines which could be applied on working BRT systems to improve the bus flow.

### [Auto-generation of a centerline graph from a geometrically complex  roadmap of real-world traffic systems using a hierarchical quadtree for  cellular automata simulations](http://arxiv.org/abs/1903.09419v2) (1903.09419v2)
<i>Satori Tsuzuki, Daichi Yanagisawa, Katsuhiro Nishinari</b>

<h10>2019-03-22</h10>
> This paper proposes a method of auto-generation of a centerline graph from a geometrically complex roadmap of real-world traffic systems by using a hierarchical quadtree for cellular automata simulations. Our method is summarized as follows. First, we store the binary values of the monochrome image of target roadmap (one and zero represent the road and the other areas, respectively) in the two-dimensional square map. Second, we recursively divide the square map into sub-leafs by a quadtree until the summed-up value of pixels included inside the leaf becomes equal to or less than one. Third, we gradually remove the distal leaves that are adjacent to the leaves whose depths are shallower than the distal leaf. After that, we trace the remaining distal leaves of the tree using Morton's space-filling curve, while selecting the leaves that keep a certain distance among the previously selected leaves as the nodes of the graph. Finally, each selected node searches the neighboring nodes and stores them as the edges of the graph. We demonstrate our method by generating a centerline graph from a complex roadmap of a real-world airport and by carrying out a typical network analysis using Dijkstra's method.

 <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>