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

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

### [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.

### [Kerr Frequency combs and dissipative Kerr solitons in non-linear  coupled-cavity waveguides](http://arxiv.org/abs/1907.01492v1) (1907.01492v1)
<i>J. P. Vasco, V. Savona</b>

<h10>2019-07-02</h10>
> We derive a non-linear formalism from first principles to study Kerr frequency combs and dissipative Kerr solitons (DKS) in coupled-cavity waveguides (CCW), with Kerr non-linearity and two-photon absorption. We arrive to a system of non-linear coupled-mode equations describing the dynamics of the slowly-varying field amplitudes of the CCW Bloch modes, where the corresponding power threshold for comb generation depends explicitly on the main CCW figures of merit, namely, non-linear mode volume, normal mode quality factor and slow-light group index. We apply our formalism to a simple tight binding model and to a globally optimized silicon photonic crystal CCW at telecom wavelengths, and find clear signatures of low-threshold frequency combs and DKS. Our result sets the CCW as a new paradigm for comb generation via dispersion engineering and slow-light non-linear enhancement.

### [Coexisting Ordered States, Local Equilibrium-like Domains, and Broken  Ergodicity in a Non-turbulent Rayleigh-Bénard Convection at Steady-state](http://arxiv.org/abs/1812.06002v3) (1812.06002v3)
<i>Atanu Chatterjee, Yash Yadati, Nicholas Mears, Germano Iannacchione</b>

<h10>2018-12-14</h10>
> A challenge in fundamental physics and especially in thermodynamics is to understand emergent order in far-from-equilibrium systems. While at equilibrium, temperature plays the role of a key thermodynamic variable whose uniformity in space and time defines the equilibrium state the system is in, this is not the case in a far-from-equilibrium driven system. When energy flows through a finite system at steady-state, temperature takes on a time-independent but spatially varying character. In this study, the convection patterns of a Rayleigh-B{\'e}nard fluid cell at steady-state is used as a prototype system where the temperature profile and fluctuations are measured spatio-temporally. The thermal data is obtained by performing high-resolution real-time infrared calorimetry on the convection system as it is first driven out-of-equilibrium when the power is applied, achieves steady-state, and then as it gradually relaxes back to room temperature equilibrium when the power is removed. Our study provides new experimental data on the non-trivial nature of thermal fluctuations when stable complex convective structures emerge. The thermal analysis of these convective cells at steady-state further yield local equilibrium-like statistics. In conclusion, these results correlate the spatial ordering of the convective cells with the evolution of the system's temperature manifold.

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

### [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.

### [Chaos and Complexity in Quantum Mechanics](http://arxiv.org/abs/1905.13534v3) (1905.13534v3)
<i>Tibra Ali, Arpan Bhattacharyya, S. Shajidul Haque, Eugene H. Kim, Nathan Moynihan, Jeff Murugan</b>

<h10>2019-05-31</h10>
> We propose a new diagnostic for quantum chaos. We show that time evolution of complexity for a particular type of target state can provide equivalent information about the classical Lyapunov exponent and scrambling time as out-of-time-order correlators. Moreover, for systems that can be switched from a regular to unstable (chaotic) regime by a tuning of the coupling constant of the interaction Hamiltonian, we find that the complexity defines a new time scale. We interpret this time scale as recording when the system makes the transition from regular to chaotic behaviour.

### [Data-driven prediction of a multi-scale Lorenz 96 chaotic system using a  hierarchy of deep learning methods: Reservoir computing, ANN, and RNN-LSTM](http://arxiv.org/abs/1906.08829v2) (1906.08829v2)
<i>Ashesh Chattopadhyay, Pedram Hassanzadeh, Krishna Palem, Devika Subramanian</b>

<h10>2019-06-20</h10>
> In this paper, the performance of three deep learning methods for predicting short-term evolution and reproducing the long-term statistics of a multi-scale spatio-temporal Lorenz 96 system is examined. The methods are: echo state network (a type of reservoir computing, RC-ESN), deep feed-forward artificial neural network (ANN), and recurrent neural network with long short-term memory (RNN-LSTM). This Lorenz system has three tiers of nonlinearly interacting variables representing slow/large-scale (![](http://latex2png.com/output//latex_4296370263158f14b42066659babdfc1.png)), intermediate (![](http://latex2png.com/output//latex_533af54fb84df0c659a7fe0c0d009d64.png)), and fast/small-scale (![](http://latex2png.com/output//latex_81a0d8d0563f6028a7dd52e0a8620057.png)) processes. For training or testing, only ![](http://latex2png.com/output//latex_4296370263158f14b42066659babdfc1.png) is available; ![](http://latex2png.com/output//latex_533af54fb84df0c659a7fe0c0d009d64.png) and ![](http://latex2png.com/output//latex_81a0d8d0563f6028a7dd52e0a8620057.png) are never known/used. It is shown that RC-ESN substantially outperforms ANN and RNN-LSTM for short-term prediction, e.g., accurately forecasting the chaotic trajectories for hundreds of numerical solver's time steps, equivalent to several Lyapunov timescales. RNN-LSTM and ANN show some prediction skills as well; RNN-LSTM bests ANN. Furthermore, even after losing the trajectory, data predicted by RC-ESN and RNN-LSTM have probability density functions (PDFs) that closely match the true PDF, even at the tails. PDF of the ANN data deviates from the true PDF. Implications, caveats, and applications to data-driven and inexact, data-assisted surrogate modeling of complex dynamical systems such as weather/climate are discussed.

### [Can quantum many-body systems behave as strongly chaotic, being  completely integrable ?](http://arxiv.org/abs/1907.01893v1) (1907.01893v1)
<i>Samy Mailoud, Fausto Borgonovi, Felix Izrailev</b>

<h10>2019-06-29</h10>
> We study the paradigmatic Lieb-Liniger (LL) model belonging to the class of integrable quantum many-body systems, by considering its statistical properties in the many-body Hilbert space. We demonstrate that, for a fixed total momentum, the properties of both energy spectra fluctuations and many-body eigenstates follow the predictions of the random matrix theory. Specifically, for a finite number of bosons in a finite number of momentum levels and sufficiently strong interaction, the level spacing distribution manifests the crossover from Poisson to Wigner-Dyson statistics. In the latter situation, the many-body eigenstates can be treated as fully random, in spite of the deterministic nature of the Hamiltonian matrix elements. By studying the quench dynamics of an initially excited state of the non interacting Hamiltonian, we have discovered a remarkable relation between the thermodynamic entropy emerging after the relaxation of the system to equilibrium, and the diagonal entropy related to the initial state. Our semi-analytical predictions of the onset of statistical relaxation in the LL model are fully confirmed by numerical data.

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