<center><b> Welcome To Nonlinear Sciences</b></center>
# <center>Adaptation And Self-Organizing Systems</center>
<hr>
### [Modeling cochlear two-tone suppression using a system of nonlinear oscillators with feed-forward coupling](http://arxiv.org/abs/1812.08587v1) (1812.08587v1)
<i>Masanori Ouchi, Hiroya Nakao</b>
<h10>2018-12-20</h10>
> Mechanism of two-tone suppression is studied using a coupled-oscillator model of the cochlea with feed-forward coupling. Local amplification of sound signals is modeled by using Stuart-Landau oscillators near the Hopf bifurcation, and transmission of sound signals is described as feed-forward coupling between the oscillators. Effect of suppressor signals on the response to probe signals is analyzed by numerical simulations. It is found that the effect of suppression is qualitatively different depending on relative frequency between probe and suppressor signals. By analyzing a simplified two-oscillator model, we explain the mechanism of the suppression, where configuration of the oscillators plays an essential role.
### [When is sync globally stable in sparse networks of identical Kuramoto oscillators?](http://arxiv.org/abs/1812.08153v1) (1812.08153v1)
<i>Yury Sokolov, G. Bard Ermentrout</b>
<h10>2018-12-19</h10>
> Synchronization in systems of coupled Kuramoto oscillators may depend on their natural frequencies, coupling, and underlying networks. In this paper, we reduce the alternatives to only one by considering identical oscillators where the only parameter that is allowed to change is the underlying network. While such a model was analyzed over the past few decades by studying the size of the basin of attraction of the synchronized state on restricted families of graphs, here we address a qualitative question on general graphs. In an analogy to resistive networks with current sources, we describe an algorithm that produces initial conditions that are often outside of the basin of attraction of the synchronized state. In particular, if a graph allows a cyclic graph clustering with a sufficient number of clusters or contains a sufficiently long induced subpath without cut vertices of the graph then there is a non-synchronous stable phase-locked solution. Thus, we provide a partial answer to when the synchronized state is not globally stable.
### [Conformal Group Actions on Generalized Kuramoto Oscillators](http://arxiv.org/abs/1812.06539v2) (1812.06539v2)
<i>Max Lipton</b>
<h10>2018-12-16</h10>
> This paper unifies the recent results on generalized Kuramoto Model reductions. Lohe took a coupled system of  bodies on  governed by the Kuramoto equations  and used the method of Watanabe and Strogatz to reduce this system to  equations. Using a model of rigid rotations on a sphere as a guide, we show that the reduction is described by a smooth path in the Lie group of conformal transformations on the sphere, which is diffeomorphic to . Seeing the reduction this way allows us to apply geometric and topological reasoning in order to understand qualitative behavior of the Kuramoto Model.
### [Spectral properties of complex networks](http://arxiv.org/abs/1810.01254v2) (1810.01254v2)
<i>Camellia Sarkar, Sarika Jalan</b>
<h10>2018-10-01</h10>
> This review presents an account of the major works done on spectra of adjacency matrices drawn on networks and the basic understanding attained so far. We have divided the review under three sections: (a) extremal eigenvalues, (b) bulk part of the spectrum and (c) degenerate eigenvalues, based on the intrinsic properties of eigenvalues and the phenomena they capture. We have reviewed the works done for spectra of various popular model networks, such as the Erd\H{o}s-R\'enyi random networks, scale-free networks, 1-d lattice, small-world networks, and various different real-world networks. Additionally, potential applications of spectral properties for natural processes have been reviewed.
### [Antifragility of Random Boolean Networks](http://arxiv.org/abs/1812.06760v1) (1812.06760v1)
<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 demonstrate 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.
# <center>Chaotic Dynamics</center>
<hr>
### [Geodesic dynamics in Chazy-Curzon spacetimes](http://arxiv.org/abs/1812.08663v1) (1812.08663v1)
<i>F. L. Dubeibe, J. D. Arias H., J. E. Alfonso</b>
<h10>2018-12-20</h10>
> In the last decades, the dynamical studies around compact objects became a subject of active research, partially motivated by the observed differences in the profiles of the gravitational waves depending on the dynamics of the system. In this work, via the Poincar\'e section method, we conduct a thorough numerical analysis of the dynamical behavior of geodesics around Chazy-Curzon metrics. As the main result, we find only regular motions for the geodesics in all cases, which suggest the existence of the so-called Carter's constant in this kind of exact solutions. Moreover, our simulations indicate that in the two-particle Chazy-Curzon solution, some oscillatory motions take place as in the classical MacMillan problem.
### [A Universal Operator Growth Hypothesis](http://arxiv.org/abs/1812.08657v1) (1812.08657v1)
<i>Daniel E. Parker, Xiangyu Cao, Thomas Scaffidi, Ehud Altman</b>
<h10>2018-12-20</h10>
> We present a hypothesis for the universal properties of operators evolving under Hamiltonian dynamics in many-body systems. The hypothesis states that successive Lanczos coefficients in the continued fraction expansion of the Green's functions grow linearly with rate  in generic systems. The rate  --- an experimental observable --- governs the exponential growth of operator complexity in a sense we make precise. This exponential growth even prevails beyond semiclassical or large- limits. Moreover,  upper bounds a large class of operator complexity measures, including the out-of-time-order correlator. As a result, we conjecture a sharp bound on Lyapunov exponents , which generalizes the known universal low-temperature bound . We illustrate our results in paradigmatic examples such as non-integrable spin chains, the Sachdev-Ye-Kitaev model, and classical models. Finally we use the hypothesis in conjunction with the recursion method to develop a technique for computing diffusion constants.
### [Equilibrium points and basins of convergence in the triangular restricted four-body problem with a radiating body](http://arxiv.org/abs/1812.08641v1) (1812.08641v1)
<i>J. E. Osorio-Vargas, Guillermo A. González, F. L. Dubeibe</b>
<h10>2018-12-20</h10>
> The dynamics of the four-body problem have attracted increasing attention in recent years. In this paper, we extend the basic equilateral four-body problem by introducing the effect of radiation pressure, Poynting-Robertson drag, and solar wind drag. In our setup, three primaries lay at the vertices of an equilateral triangle and move in circular orbits around their common center of mass. Here, one of the primaries is a radiating body and the fourth body (whose mass is negligible) does not affect the motion of the primaries. We show that the existence and the number of equilibrium points of the problem depend on the mass parameters and radiation factor. Consequently, the allowed regions of motion, the regions of the basins of convergence for the equilibrium points, and the basin entropy will also depend on these parameters. The present dynamical model is analyzed for three combinations of mass for the primaries: equal masses, two equal masses, different masses. As the main results, we find that in all cases the libration points are unstable if the radiation factor is larger than 0.01 and hence able to destroy the stability of the libration points in the restricted four-body problem composed by Sun, Jupiter, Trojan asteroid and a test (dust) particle. Also, we conclude that the number of fixed points decreases with the increase of the radiation factor.
### [Stochastic resonance in a sinusoidal potential system: An analog simulation experiment](http://arxiv.org/abs/1608.06138v2) (1608.06138v2)
<i>Ivan Skhem Sawkmie, Mangal C. Mahato</b>
<h10>2016-08-22</h10>
> Recently, stochastic resonance was obtained numerically in an underdamped periodic potential system driven by a periodic force and a Gaussian white noise. In that numerical work, the occurrence of stochastic resonance was explained in terms of the existence of two dynamical states having different amplitude and phase lag. At zero temperature these two initial condition dependent dynamical states are stable. However, at elevated temperatures, these two states make transitions from one to the other at a mean rate. In the present work, we setup an analog simulation experiment to show the existence of the two dynamical states in a sinusoidal potential system as well as to verify the occurrence of stochastic resonance in the same system. The experimental procedure includes setting the initial conditions for the experiment.
### [Universal route to thermalization in weakly-nonlinear one-dimensional chains](http://arxiv.org/abs/1812.08279v1) (1812.08279v1)
<i>Lorenzo Pistone, Sergio Chibbaro, Miguel Bustamante, Yuri L'vov, Miguel Onorato</b>
<h10>2018-12-19</h10>
> We apply Wave Turbulence theory to describe the dynamics on nonlinear one-dimensional chains. We consider  and  Fermi-Pasta-Ulam-Tsingou (FPUT) systems, and the discrete nonlinear Klein-Gordon chain. We demonstrate that resonances are responsible for the irreversible transfer of energy among the Fourier modes. We predict that all the systems thermalize for large times, and that the equipartition time scales as a power-law of the strength of the nonlinearity. Our methodology is not limited to only these systems and can be applied to the case of a finite number of modes, such as in the original FPUT experiment, or to the thermodynamic limit, i.e. when the number of modes approach infinity. In the latter limit, we perform state of the art numerical simulations and show that the results are consistent with theoretical predictions. We suggest that the route to thermalization, based only on the presence of exact resonance, has universal features. Moreover, a by-product of our analysis is the asymptotic integrability, up to four wave interactions, of the discrete nonlinear Klein-Gordon chain.
# <center>Cellular Automata And Lattice Gases</center>
<hr>
### [Antifragility of Random Boolean Networks](http://arxiv.org/abs/1812.06760v1) (1812.06760v1)
<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 demonstrate 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.
### [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.
# <center>Pattern Formation And Solitons</center>
<hr>
### [Phase-space geometry of reaction--diffusion dynamics](http://arxiv.org/abs/1812.08684v1) (1812.08684v1)
<i>Fridtjof Brauns, Jacob Halatek, Erwin Frey</b>
<h10>2018-12-20</h10>
> Experimental studies of protein pattern formation (both in vivo and in vitro) have stimulated new interest in the dynamics of reaction--diffusion systems. However, a comprehensive theoretical understanding of the dynamics of such highly nonlinear, spatially extended systems is still missing. Here we show how a description in phase space, which has proven invaluable in shaping our intuition about the dynamics of nonlinear ordinary differential equations, can be generalized to mass-conserving reaction--diffusion (McRD) systems. We present a comprehensive theory for two-component McRD systems, which serve as paradigmatic minimal systems that encapsulate the core principles and concepts of our framework. The key insight is that shifting local chemical equilibria --- controlled by the local total density --- give rise to concentration gradients that drive diffusive redistribution of total density. This dynamic interplay can be embedded in the phase plane of the reaction kinetics in terms of simple geometric objects: the reactive nullcline (line of chemical equilibria) and the diffusive flux-balance subspace. On this phase-space level, physical insight can be gained from geometric criteria and graphical constructions --- the effects of nonlinearities on the global dynamics are simply encoded in the curved shape of the reactive nullcline. In particular, we show that the pattern-forming 'Turing instability' in McRD systems is a mass-redistribution instability, and that the features and bifurcations of patterns can be graphically determined by a flux-balance construction on the reactive nullcline. The fundamental elements of the theory presented suggest ways of experimentally characterizing pattern-forming systems on a mesoscopic level and are generalizable to a broad class of spatially extended non-equilibrium systems, and thereby pave the way toward an overarching theoretical framework.
### [Controlling Dispersive Hydrodynamic Wavebreaking in a Viscous Fluid Conduit](http://arxiv.org/abs/1812.08134v1) (1812.08134v1)
<i>Dalton V. Anderson, Michelle D. Maiden, Mark A. Hoefer</b>
<h10>2018-12-19</h10>
> The driven, cylindrical, free interface between two miscible, Stokes fluids with high viscosity contrast have been shown to exhibit dispersive hydrodynamics. A hallmark feature of dispersive hydrodynamic media is the dispersive resolution of wavebreaking that results in a dispersive shock wave. In the context of the viscous fluid conduit system, the present work introduces a simple, practical method to precisely control the location, time, and spatial profile of wavebreaking in dispersive hydrodynamic systems with only boundary control. The method is based on tracking the dispersionless characteristics backward from the desired wavebreaking profile to the boundary. In addition to the generation of step-like Riemann and box problems, the method is generalized to other, piecewise-linear dispersive hydrodynamic profiles including the triangle wave and N-wave. Quantitative agreement between the predicted location and time of dispersive wavebreaking and viscous fluid conduit experiment for a range of flow conditions is obtained. Observed space-time characteristics also agree with triangle and N-wave predictions. The characteristic boundary control method introduced here enables the experimental investigation of a variety of wavebreaking profiles and is expected to be useful in other dispersive hydrodynamic media.
### [Spontaneous center formation in Dictyostelium discoideum](http://arxiv.org/abs/1812.08061v1) (1812.08061v1)
<i>Estefania Vidal-Henriquez, Azam Gholami</b>
<h10>2018-12-19</h10>
> Dictyostelium discoideum (D.d.) is a widely studied amoeba due to its capabilities of development, survival, and self-organization. During aggregation it produces and relays a chemical signal (cAMP) which shows spirals and target centers. Nevertheless, the natural emergence of these structures is still not well understood. We present a mechanism for creation of centers and target waves of cAMP in D.d. by adding cell inhomogeneity to a well known reaction-diffusion model of cAMP waves and we characterize its properties. We show how stable activity centers appear spontaneously in areas of higher cell density with the oscillation frequency of these centers depending on their density. The cAMP waves have the characteristic dispersion relation of trigger waves and a velocity which increases with cell density. Chemotactically competent cells react to these waves and create aggregation streams even with very simple movement rules. Finally we argue in favor of the existence of bounded phosphodiesterase to maintain the wave properties once small cell clusters appear.
### [Stability of bubble-like fluxons in disk-shaped Josephson junctions in the presence of a coaxial dipole current](http://arxiv.org/abs/1812.07980v1) (1812.07980v1)
<i>Alicia G. Castro-Montes, Mónica A. García-Ñustes, Jorge A. González, Juan F. Marín, Diego Teca-Wellmann</b>
<h10>2018-12-19</h10>
> We investigate analytically and numerically the stability of bubble-like fluxons in disk-shaped heterogeneous Josephson junctions. Using ring solitons as a model of bubble fluxons in the two-dimensional sine-Gordon equation, we show that the insertion of coaxial dipole currents prevents their collapse. We characterize the onset of instability by introducing a single parameter that couples the radius of the bubble fluxon with the properties of the injected current. For different combination of parameters, we report the formation of stable oscillating bubbles, the emergence of internal modes, and bubble breakup due to internal mode instability. Our theoretical findings are in good agreement with numerical simulations. We discuss applications in quantum information technologies.
### [Embedded solitons in the double sine-Gordon lattice with the next neighbor interactions](http://arxiv.org/abs/1812.07856v1) (1812.07856v1)
<i>Yaroslav Zolotaryuk, Ivan O. Starodub</b>
<h10>2018-12-19</h10>
> Topological solitons can propagate without radiation in discrete media. These solutions are known as embedded solitons (ES). They come as isolated solutions and exist despite their resonance with the linear spectrum of the respective lattices. In this paper the properties of embedded solitons in the discrete double sine-Gordon equation with the next-neighbor and second-neighbor interactions are investigated. Depending on the sign of these interactions they can be either destructive or favorable for the ES creation. The ES existence area depends on the width of the linear spectrum: narrowing of the spectrum widens the ES existence range and vice versa. The application to the Josephson junction arrays is discussed.
# <center>Exactly Solvable And Integrable Systems</center>
<hr>
### [On symmetries of the Gibbons-Tsarev equation](http://arxiv.org/abs/1811.08199v2) (1811.08199v2)
<i>H. Baran, P. Blaschke, I. S. Krasil'shchik, M. Marvan</b>
<h10>2018-11-20</h10>
> We study the Gibbons-Tsarev equation  and, using the known Lax pair, we construct infinite series of conservation laws and the algebra of nonlocal symmetries in the covering associated with these conservation laws. We prove that the algebra is isomorphic to the Witt algebra. Finally, we show that the constructed symmetries are unique in the class of polynomial ones.
### [Cluster Toda chains and Nekrasov functions](http://arxiv.org/abs/1804.10145v2) (1804.10145v2)
<i>M. Bershtein, P. Gavrylenko, A. Marshakov</b>
<h10>2018-04-26</h10>
> In this paper the relation between the cluster integrable systems and -difference equations is extended beyond the Painlev\'e case. We consider the class of hyperelliptic curves when the Newton polygons contain only four boundary points. The corresponding cluster integrable Toda systems are presented, and their discrete automorphisms are identified with certain reductions of the Hirota difference equation. We also construct non-autonomous versions of these equations and find that their solutions are expressed in terms of 5d Nekrasov functions with the Chern-Simons contributions, while in the autonomous case these equations are solved in terms of the Riemann theta-functions.
### [On tau-functions for the KdV hierarchy](http://arxiv.org/abs/1812.08488v1) (1812.08488v1)
<i>Boris Dubrovin, Di Yang, Don Zagier</b>
<h10>2018-12-20</h10>
> For an  solution to the KdV hierarchy, the generating series of logarithmic derivatives of the tau-function of the solution can be expressed by the basic matrix resolvent via algebraic manipulations [5,7]. Based on this we develop in this paper two new formulae for the generating series by introducing the notion of  of wave functions of the solution. Applications to the Witten--Kontsevich tau-function, to the generalized Br\'ezin--Gross--Witten (BGW) tau-function, as well as to a modular deformation of the generalized BGW tau-function which we call the Lam\'e tau-function are also given.
### [Symbolic Computations of First Integrals for Polynomial Vector Fields](http://arxiv.org/abs/1710.08225v2) (1710.08225v2)
<i>Guillaume Chèze, Thierry Combot</b>
<h10>2017-10-23</h10>
> In this article we show how to generalize to the Darbouxian, Liouvillian and Riccati case the extactic curve introduced by J. Pereira. With this approach, we get new algorithms for computing, if it exists, a rational, Darbouxian, Liouvillian or Riccati first integral with bounded degree of a polynomial planar vector field. We give probabilistic and deterministic algorithms. The arithmetic complexity of our probabilistic algorithm is in , where  is the bound on the degree of a representation of the first integral and  is the exponent of linear algebra. This result improves previous algorithms. Our algorithms have been implemented in Maple and are available on authors' websites. In the last section, we give some examples showing the efficiency of these algorithms.
### [Quantum line operators from Lax pairs](http://arxiv.org/abs/1812.06701v1) (1812.06701v1)
<i>E. H Saidi</b>
<h10>2018-12-17</h10>
> Motivated by the realisation of Yang-Baxter equation of 2d Integrable models in the 4d gauge theory of Costello-Witten-Yamazaki (CWY), we study the embedding of integrable 2d Toda field models inside this construction. This is done by using the Lax formulation of 2d integrable systems and by thinking of the standard Lax pair  in terms of components of CWY gauge connection propagating along particular directions in the gauge bundle. We also use results of the CWY theory to build quantum line operators for 2d Toda models and compute the one loop contribution of two intersecting lines exchanging one gluon. Other features like local symmetries and comments on extension of our method to other 2d integrable models are also discussed.
<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