<center><b> Welcome To Nonlinear Sciences</b></center>
# <center>Adaptation And Self-Organizing Systems</center>
<hr>
### [Coexisting Ordered States, Local Equilibrium Points, and Broken Ergodicity in a Non-turbulent Rayleigh-B{é}nard Convection at Steady-state](http://arxiv.org/abs/1812.06002v1) (1812.06002v1)
<i>Atanu Chatterjee, Yash Yadati, Nicholas Mears, Germano Iannacchione</b>
<h10>2018-12-14</h10>
> The branch of thermodynamics permeates all physical processes in the universe. Historically, its foundations lie in the equilibrium description of matter, where states are stationary, and forces and potentials within a system are completely balanced. Systems in nature however operate under conditions that are far-from-equilibrium. One of the major reasons that prevents developing a general understanding of out-of-equilibrium thermodynamics is the incomplete definition of temperature, a key thermodynamic parameter. In this study, we revisit the notion of temperature and thermal fluctuations in a Rayleigh-B{\'e}nard convection that has been driven out-of-equilibrium. We perform a series of analysis, both temporally and spatially on the thermal images that were obtained from real-time Infra Red calorimetry of the Rayleigh-B{\'e}nard system as it evolved to an out-of-equilibrium steady-state, and gradually relaxes back to room temperature equilibrium. Our study uncovers exciting results about how temperature should be interpreted far-from-equilibrium, while providing insights on how local equilibrium points can not only coexist, but thrive in a macroscopically driven system.
### [Repulsively Coupled Kuramoto-Sakaguchi Phase Oscillators Ensemble Subject to Common Noise](http://arxiv.org/abs/1812.05883v1) (1812.05883v1)
<i>Chen Chris Gong, Chunming Zheng, Ralf Toenjes, Arkady Pikovsky</b>
<h10>2018-12-14</h10>
> Phase oscillators evolving by a M\"obius transformation of the unit circle according to the Watanabe-Strogatz theory in principle cannot contract in phase space to form more than one cluster. In numerical simulation, however, Kuramoto-Sakaguchi phase oscillators under common multiplicative noise, which alone leads to the synchronization of the oscillators, and under repulsive coupling, which alone leads to the desynchronization of the oscillators, have been shown to form stable multicluster states. Using Fokker-Planck formulation, we show that two-cluster states under such a model are non-attractive. Integrating with various time steps reveals that clustering is a numerical artifact, explained by the existence of higher order Fourier terms in the errors of the employed numerical integration schemes. Monitoring the induced change in certain integrals of motion we quantify these errors. On the other hand, in ensembles of general limit cycle oscillators, such as Van der Pol oscillators, due to an anharmonic phase response function, as well as additional amplitude dynamics, multiclusters can occur naturally. Thus that common noise can induce clustering in general oscillator systems under repulsive coupling remains valid.
### [Local Symmetry and Global Structure in Adaptive Voter Models](http://arxiv.org/abs/1812.05464v1) (1812.05464v1)
<i>Philip S. Chodrow, Peter J. Mucha</b>
<h10>2018-12-13</h10>
> "Coevolving" or "adaptive" voter models (AVMs) are natural systems for modeling the emergence of mesoscopic structure from local networked processes driven by conflict and homophily. Because of this, many methods for approximating the long-run behavior of AVMs have been proposed over the last decade. However, most such methods are either restricted in scope, expensive in computation, or inaccurate in predicting important statistics. In this work, we develop a novel, second-order moment closure approximation method for studying the equilibrium mesoscopic structure of AVMs and apply it to binary-state rewire-to-random and rewire-to-same model variants with random state-switching. This framework exploits an asymmetry in voting events that enables us to derive analytic approximations for the fast-timescale dynamics. The resulting numerical approximations enable the computation of key properties of the model behavior, such as the location of the fragmentation transition and the equilibrium active edge density, across the entire range of state densities. Numerically, they are nearly exact for the rewire-to-random model, and competitive with other current approaches for the rewire-to-same model. We conclude with suggestions for model refinement and extensions to more complex models.
### [Graph Comparison via the Non-backtracking Spectrum](http://arxiv.org/abs/1812.05457v1) (1812.05457v1)
<i>Andrew Mellor, Angelica Grusovin</b>
<h10>2018-12-13</h10>
> The comparison of graphs is a vitally important, yet difficult task which arises across a number of diverse research areas including biological and social networks. There have been a number of approaches to define graph distance however often these are not metrics (rendering standard data-mining techniques infeasible), or are computationally infeasible for large graphs. In this work we define a new metric based on the spectrum of the non-backtracking graph operator and show that it can not only be used to compare graphs generated through different mechanisms, but can reliably compare graphs of varying size. We observe that the family of Watts-Strogatz graphs lie on a manifold in the non-backtracking spectral embedding and show how this metric can be used in a standard classification problem of empirical graphs.
### [On the potential for open-endedness in neural networks](http://arxiv.org/abs/1812.04907v1) (1812.04907v1)
<i>Nicholas Guttenberg, Nathaniel Virgo, Alexandra Penn</b>
<h10>2018-12-12</h10>
> Natural evolution gives the impression of leading to an open-ended process of increasing diversity and complexity. If our goal is to produce such open-endedness artificially, this suggests an approach driven by evolutionary metaphor. On the other hand, techniques from machine learning and artificial intelligence are often considered too narrow to provide the sort of exploratory dynamics associated with evolution. In this paper, we hope to bridge that gap by reviewing common barriers to open-endedness in the evolution-inspired approach and how they are dealt with in the evolutionary case - collapse of diversity, saturation of complexity, and failure to form new kinds of individuality. We then show how these problems map onto similar issues in the machine learning approach, and discuss how the same insights and solutions which alleviated those barriers in evolutionary approaches can be ported over. At the same time, the form these issues take in the machine learning formulation suggests new ways to analyze and resolve barriers to open-endedness. Ultimately, we hope to inspire researchers to be able to interchangeably use evolutionary and gradient-descent-based machine learning methods to approach the design and creation of open-ended systems.
# <center>Chaotic Dynamics</center>
<hr>
### [Diffusion on hierarchical systems of weakly-coupled networks](http://arxiv.org/abs/1303.2650v4) (1303.2650v4)
<i>Grzegorz Siudem, Janusz A. Hołyst</b>
<h10>2013-03-10</h10>
> We analyse diffusion dynamics on weakly-coupled networks (interconnected networks) by means of separation of time scales. Using an adiabatic approximation we reduced the system dynamics to a Markov chain with aggregated variables and derived a transport equation that is analogous to Fick's First Law and includes a driving force. Entropy production is a sum of microscopic entropy transport, which results from the particle's migration between networks of different topologies and macroscopic entropy production of the Markov chain. Equilibrium particles partition between different sub-networks depends only on internal sub-network parameters. Our framework, confirmed by numerical simulations, is also useful for considering diffusion in nested systems corresponding to hierarchical networks with several different time scales thus it can serve to uncover hidden hierarchy levels from observations of diffusion processes.
### [Symmetry broken states in an ensemble of globally coupled pendulums](http://arxiv.org/abs/1812.05141v1) (1812.05141v1)
<i>M. I. Bolotov, V. O. Munyaev, L. A. Smirnov, A. E. Hramov</b>
<h10>2018-12-12</h10>
> We consider the rotational dynamics in an ensemble of identical globally coupled pendulums. This model is essentially a generalization of the standard Kuramoto model, which takes into account the inertia and the intrinsic nonlinearity of the community elements. There exist the wide variety of in-phase and out-of-phase regimes. Many of these states appear due to instability of in-phase rotational mode that leads to partially broken symmetry in the considered model. Our theoretical analysis allows one to find the boundaries of the instability domain for ensembles with arbitrary number of pendulums. In the case of three elements, for system control parameter sets corresponding to the unstable in-phase rotations we numerically find a number of out-of-phase rotation modes and study in detail their stability and bifurcations. We demonstrate that chaotic behavior arises throw the cascade period doubling bifurcations and torus destruction. As a result, we obtain a sufficiently detailed picture of the symmetry breaking and existence of various regular and chaotic states in an ensemble of identical globally coupled pendulums. In particular, we show that the family of symmetry-broken states includes chimera attractors. In these regimes, one part of the oscillators moving irregularly is synchronized, and the other pendulums rotate separately. We directly demonstrate the example of chimera states in the system under discussion for the cases of three and four number of elements.
### [Entanglement spreading in a minimal model of maximal many-body quantum chaos](http://arxiv.org/abs/1812.05090v1) (1812.05090v1)
<i>Bruno Bertini, Pavel Kos, Tomaz Prosen</b>
<h10>2018-12-12</h10>
> We study the entanglement dynamics in the self-dual kicked Ising model, which we recently identified as a minimal model for maximal many-body quantum chaos. We consider the time-evolution of the entanglement entropies of finite blocks of spins after the system is initialised in a class of factorised initial states. Generalising the recently introduced "duality method" we determine the time evolution of all R\'enyi entropies exactly. This further confirms the special status of the self-dual kicked Ising model. Even though the model is non-integrable and displays spectral correlations described by the circular orthogonal ensemble of random matrices, it allows for exact calculations. Our results saturate the so called "minimal cut" bound and are in agreement with those found in the contexts of random unitary circuits with infinite-dimensional local Hilbert space and conformal field theory.
### [Permutation entropy revisited](http://arxiv.org/abs/1812.05075v1) (1812.05075v1)
<i>Stuart J Watt, Antonio Politi</b>
<h10>2018-12-12</h10>
> Time-series analysis in terms of ordinal patterns is revisited by introducing a generalized permutation entropy , which depends on two different window lengths: , implicitly defining the resolution of the underlying partition; , playing the role of an embedding dimension, analogously to standard nonlinear time-series analysis. The -dependence provides information on the structure of the corresponding invariant measure, while the -dependence helps determining the Kolmogorov-Sinai entropy. We finally investigate the structure of the partition with the help of principal component analysis, finding that, upon increasing , the single atoms become increasingly elongated.
### [Chaotic few-body vortex dynamics in rotating Bose--Einstein condensates](http://arxiv.org/abs/1812.04759v1) (1812.04759v1)
<i>Tiantian Zhang, James Schloss, Andreas Thomasen, Lee James O'Riordan, Thomas Busch, Angela White</b>
<h10>2018-12-12</h10>
> We investigate a small vortex-lattice system of four co-rotating vortices in an atomic Bose--Einstein condensate and find that the vortex dynamics display chaotic behaviour after a system quench introduced by reversing the direction of circulation of a single vortex through a phase-imprinting process. By tracking the vortex trajectories and Lyapunov exponent, we show the onset of chaotic dynamics is not immediate, but occurs at later times and is accelerated by the close-approach and separation of all vortices in a scattering event. The techniques we develop could potentially be applied to create locally induced chaotic dynamics in larger lattice systems as a stepping stone to study the role of chaotic events in turbulent vortex dynamics.
# <center>Cellular Automata And Lattice Gases</center>
<hr>
### [Lenia - Biology of Artificial Life](http://arxiv.org/abs/1812.05433v1) (1812.05433v1)
<i>Bert Wang-Chak Chan</b>
<h10>2018-12-13</h10>
> We report a new model of artificial life called Lenia (from Latin lenis "smooth"), a two-dimensional cellular automaton with continuous space-time-state and generalized local rule. Computer simulations show that Lenia supports a great diversity of complex autonomous patterns or "lifeforms" bearing resemblance to real-world microscopic organisms. More than 400 species in 18 families have been identified, many discovered via interactive evolutionary computation. We present basic observations of the model regarding the properties of space-time and basic settings. We provide a board survey of the lifeforms, categorize them into a hierarchical taxonomy, and map their distribution in the parameter hyperspace. We describe their morphological structures and behavioral dynamics, propose possible mechanisms of their self-propulsion, self-organization and plasticity. Finally, we discuss how the study of Lenia would be related to biology, artificial life, and artificial intelligence.
### [A decidability result for the halting problem of cellular automata in the pentagrid](http://arxiv.org/abs/1812.03882v1) (1812.03882v1)
<i>Maurice Margenstern</b>
<h10>2018-12-10</h10>
> In this paper, we investigate the halting problem for deterministic cellula automata in the pentagrid. We prove that the problem is decidable when the cellular automaton starts its computation from a finite configuration and when it has at most two states, one of them being a quiescent state.
### [One-sided scaling limit of multicolor box-ball system](http://arxiv.org/abs/1808.08074v4) (1808.08074v4)
<i>Atsuo Kuniba, Hanbaek Lyu</b>
<h10>2018-08-24</h10>
> 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. Using a Markov chain method, we show that if we scale the rows of each of the invariant Young diagrams by , it converges to some limiting shape as  at an exponential rate. Furthermore, we determine the limiting shape by ratios of Schur polynomials with initial ball densities as parameters. We also derive similar results through an alternative method using the Fermionic form and Thermodynamic Bethe Ansatz, which apply once we condition the initial measure on the set of highest states. By a large deviations principle, we identify the limiting shapes of invariant Young diagrams corresponding to the unconditioned and conditioned initial measures.
### [The cellular automaton pulsing model, experiments with DDLab](http://arxiv.org/abs/1811.12307v1) (1811.12307v1)
<i>Andrew Wuensche, Edward Coxon</b>
<h10>2018-11-29</h10>
> 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](http://arxiv.org/abs/1811.08310v1) (1811.08310v1)
<i>Bastien Fernandez</b>
<h10>2018-11-20</h10>
> 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.
# <center>Pattern Formation And Solitons</center>
<hr>
### [Coexisting Ordered States, Local Equilibrium Points, and Broken Ergodicity in a Non-turbulent Rayleigh-B{é}nard Convection at Steady-state](http://arxiv.org/abs/1812.06002v1) (1812.06002v1)
<i>Atanu Chatterjee, Yash Yadati, Nicholas Mears, Germano Iannacchione</b>
<h10>2018-12-14</h10>
> The branch of thermodynamics permeates all physical processes in the universe. Historically, its foundations lie in the equilibrium description of matter, where states are stationary, and forces and potentials within a system are completely balanced. Systems in nature however operate under conditions that are far-from-equilibrium. One of the major reasons that prevents developing a general understanding of out-of-equilibrium thermodynamics is the incomplete definition of temperature, a key thermodynamic parameter. In this study, we revisit the notion of temperature and thermal fluctuations in a Rayleigh-B{\'e}nard convection that has been driven out-of-equilibrium. We perform a series of analysis, both temporally and spatially on the thermal images that were obtained from real-time Infra Red calorimetry of the Rayleigh-B{\'e}nard system as it evolved to an out-of-equilibrium steady-state, and gradually relaxes back to room temperature equilibrium. Our study uncovers exciting results about how temperature should be interpreted far-from-equilibrium, while providing insights on how local equilibrium points can not only coexist, but thrive in a macroscopically driven system.
### [Solitary waves in the Ablowitz-Ladik equation with power-law nonlinearity](http://arxiv.org/abs/1806.10898v3) (1806.10898v3)
<i>J. Cuevas-Maraver, P. G. Kevrekidis, B. A. Malomed, L. Guo</b>
<h10>2018-06-28</h10>
> We introduce a generalized version of the Ablowitz-Ladik model with a power-law nonlinearity, as a discretization of the continuum nonlinear Schr\"{o}dinger equation with the same type of the nonlinearity. The model opens a way to study the interplay of discreteness and nonlinearity features. We identify stationary discrete-soliton states for different values of nonlinearity power , and address changes of their stability as frequency  of the standing wave varies for given . Along with numerical methods, a variational approximation is used to predict the form of the discrete solitons, their stability changes, and bistability features by means of the Vakhitov-Kolokolov criterion (developed from the first principles). Development of instabilities and the resulting asymptotic dynamics are explored by means of direct simulations.
### [Lenia - Biology of Artificial Life](http://arxiv.org/abs/1812.05433v1) (1812.05433v1)
<i>Bert Wang-Chak Chan</b>
<h10>2018-12-13</h10>
> We report a new model of artificial life called Lenia (from Latin lenis "smooth"), a two-dimensional cellular automaton with continuous space-time-state and generalized local rule. Computer simulations show that Lenia supports a great diversity of complex autonomous patterns or "lifeforms" bearing resemblance to real-world microscopic organisms. More than 400 species in 18 families have been identified, many discovered via interactive evolutionary computation. We present basic observations of the model regarding the properties of space-time and basic settings. We provide a board survey of the lifeforms, categorize them into a hierarchical taxonomy, and map their distribution in the parameter hyperspace. We describe their morphological structures and behavioral dynamics, propose possible mechanisms of their self-propulsion, self-organization and plasticity. Finally, we discuss how the study of Lenia would be related to biology, artificial life, and artificial intelligence.
### [Solitary Waves of the Camassa-Holm Derivative Nonlinear Schrödinger equation](http://arxiv.org/abs/1812.05149v1) (1812.05149v1)
<i>L. J. Guo, C. B. Ward, I. K. Mylonas, P. G. Kevrekidis</b>
<h10>2018-12-12</h10>
> In this paper we examine a deformation of the derivative nonlinear Schr\"odinger (DNLS) equation, the so-called Camassa-Holm DNLS (CH-DNLS) equation. We use two asymptotic multiscale expansion methods to reduce this model to both the modified Korteweg-de Vries (MKdV) equation and the Korteweg-de Vries(KdV) equation. Using their exact soliton solutions, we construct approximate solutions of the original CH-DNLS equation, namely dark and anti-dark solitary waves. Direct numerical simulations demonstrate the validity of these approximate solutions and illustrate their dynamical evolution, including their potential for nearly elastic head-on collisions in the case examples considered.
### [Desynchronization in heterogeneous networks due to uniform exposure](http://arxiv.org/abs/1812.05122v1) (1812.05122v1)
<i>V. A. Kostin, G. V. Osipov</b>
<h10>2018-12-12</h10>
> We study synchronization and rhythmic patterns generated in the heterogeneous cluster of FitzHughNagumo oscillators with transition between self-oscillating and excitable elements. Such cluster models the sinoatrial node of the heart, and the particular point of focus is the role of different physiological responses to mechanical stress: stretch-activated current and conductivity change. The comparison between them finds the difference which can serve to identify the dominant mechanism in specific experimental situations.
# <center>Exactly Solvable And Integrable Systems</center>
<hr>
### [Splitting of surface defect partition functions and integrable systems](http://arxiv.org/abs/1709.04926v2) (1709.04926v2)
<i>Saebyeok Jeong</b>
<h10>2017-09-14</h10>
> We study Bethe/gauge correspondence at the special locus of Coulomb moduli where the integrable system exhibits the splitting of degenerate levels. For this investigation, we consider the four-dimensional pure  supersymmetric  gauge theory, with a half-BPS surface defect constructed with the help of an orbifold or a degenerate gauge vertex. We show that the non-perturbative Dyson-Schwinger equations imply the Schr\"odinger-type and the Baxter-type differential equations satisfied by the respective surface defect partition functions. At the special locus of Coulomb moduli the surface defect partition function splits into parts. We recover the Bethe/gauge dictionary for each summand.
### [Repulsively Coupled Kuramoto-Sakaguchi Phase Oscillators Ensemble Subject to Common Noise](http://arxiv.org/abs/1812.05883v1) (1812.05883v1)
<i>Chen Chris Gong, Chunming Zheng, Ralf Toenjes, Arkady Pikovsky</b>
<h10>2018-12-14</h10>
> Phase oscillators evolving by a M\"obius transformation of the unit circle according to the Watanabe-Strogatz theory in principle cannot contract in phase space to form more than one cluster. In numerical simulation, however, Kuramoto-Sakaguchi phase oscillators under common multiplicative noise, which alone leads to the synchronization of the oscillators, and under repulsive coupling, which alone leads to the desynchronization of the oscillators, have been shown to form stable multicluster states. Using Fokker-Planck formulation, we show that two-cluster states under such a model are non-attractive. Integrating with various time steps reveals that clustering is a numerical artifact, explained by the existence of higher order Fourier terms in the errors of the employed numerical integration schemes. Monitoring the induced change in certain integrals of motion we quantify these errors. On the other hand, in ensembles of general limit cycle oscillators, such as Van der Pol oscillators, due to an anharmonic phase response function, as well as additional amplitude dynamics, multiclusters can occur naturally. Thus that common noise can induce clustering in general oscillator systems under repulsive coupling remains valid.
### [Periodic parabola solitons for the nonautonomous KP equation](http://arxiv.org/abs/1812.05321v1) (1812.05321v1)
<i>Yingyou Ma, Zhiqiang Chen, Xin Yu</b>
<h10>2018-12-13</h10>
> Kadomtsev-Petviashvili (KP) equation, who can describe different models in fluids and plasmas, has drawn investigation for its solitonic solutions with various methods. In this paper, we focus on the periodic parabola solitons for the (2+1) dimensional nonautonomous KP equations where the necessary constraints of the parameters are figured out. With Painleve analysis and Hirota bilinear method, we find that the solution has six undetermined parameters as well as analyze the features of some typical cases of the solutions. Based on the constructed solutions, the conditions of their convergence are also discussed.
### [Noncommutative Painlevé equations and systems of Calogero type](http://arxiv.org/abs/1710.00736v3) (1710.00736v3)
<i>Marco Bertola, Mattia Cafasso, Vladimir Roubtsov</b>
<h10>2017-10-02</h10>
> All Painlev\'e equations can be written as a time-dependent Hamiltonian system, and as such they admit a natural generalization to the case of several particles with an interaction of Calogero type (rational, trigonometric or elliptic). Recently, these systems of interacting particles have been proved to be relevant in the study of -models. An almost two decade old open question by Takasaki asks whether these multi-particle systems can be understood as isomonodromic equations, thus extending the Painlev\'e correspondence. In this paper we answer in the affirmative by displaying explicitly suitable isomonodromic Lax pair formulations. As an application of the isomonodromic representation we provide a construction based on discrete Schlesinger transforms, to produce solutions for these systems for special values of the coupling constants, starting from uncoupled ones; the method is illustrated for the case of the second Painlev\'e equation.
### [Conservation laws of the system of equations of one-dimensional shallow water over uneven bottom in Lagrange's variables](http://arxiv.org/abs/1812.05096v1) (1812.05096v1)
<i>Alexander V. Aksenov, Konstantin P. Druzhkov</b>
<h10>2018-12-12</h10>
> The system of equations of one-dimensional shallow water over uneven bottom in Euler's and Lagrange's variables is considered. Intermediate system of equations is introduced. Hydrodynamic conservation laws of intermediate system of equations is used to find all first order conservation laws of shallow water equations in Lagrange's variable for all bottom profiles. The obtained conservation laws are compared with the hydrodynamic conservation laws of the system of equations of one-dimensional shallow water over uneven bottom in Euler's variables. Bottom profiles are given for which there are additional conservation laws.
<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>
hiveblocks