<center><b> Welcome To Nonlinear Sciences</b></center>
# <center>Adaptation And Self-Organizing Systems</center>
<hr>
### [Emergence of Exploitation as Symmetry Breaking in Iterated Prisoner's Dilemma](http://arxiv.org/abs/1905.06602v1) (1905.06602v1)
<i>Yuma Fujimoto, Kunihiko Kaneko</b>
<h10>2019-05-16</h10>
> In society, mutual cooperation, defection, and asymmetric exploitative relationships are common. Whereas cooperation and defection are studied extensively in the literature on game theory, asymmetric exploitative relationships between players are little explored. In a recent study, Press and Dyson demonstrate that if only one player can learn about the other, asymmetric exploitation is achieved in the prisoner's dilemma game. In contrast, however, it is unknown whether such one-way exploitation is stably established when both players learn about each other symmetrically and try to optimize their payoffs. Here, we first formulate a dynamical system that describes the change in a player's probabilistic strategy with reinforcement learning to obtain greater payoffs, based on the recognition of the other player. By applying this formulation to the standard prisoner's dilemma game, we numerically and analytically demonstrate that an exploitative relationship can be achieved despite symmetric strategy dynamics and symmetric rule of games. This exploitative relationship is stable, even though the exploited player, who receives a lower payoff than the exploiting player, has optimized the own strategy. Whether the final equilibrium state is mutual cooperation, defection, or exploitation, crucially depends on the initial conditions: Punishment against a defector oscillates between the players, and thus a complicated basin structure to the final equilibrium appears. In other words, slight differences in the initial state may lead to drastic changes in the final state. Considering the generality of the result, this study provides a new perspective on the origin of exploitation in society.
### [Phase-locked states in oscillating neural networks and their role in neural communication](http://arxiv.org/abs/1905.06038v1) (1905.06038v1)
<i>Gemma Huguet, Alberto Pérez-Cervera, Tere M-Seara</b>
<h10>2019-05-15</h10>
> The theory of communication through coherence (CTC) proposes that brain oscillations reflect changes in the excitability of neurons, and therefore the successful communication between two oscillating neural populations depends not only on the strength of the signal emitted but also on the relative phases between them. More precisely, effective communication occurs when the emitting and receiving populations are properly phase locked so the inputs sent by the emitting population arrive at the phases of maximal excitability of the receiving population. To study this setting, we consider a population rate model consisting of excitatory and inhibitory cells modelling the receiving population, and we perturb it with a time-dependent periodic function modelling the input from the emitting population. We consider the stroboscopic map for this system and compute numerically the fixed and periodic points of this map and their bifurcations as the amplitude and the frequency of the perturbation are varied. From the bifurcation diagram, we identify the phase-locked states as well as different regions of bistability. We explore carefully the dynamics emphasizing its implications for the CTC theory. In particular, we study how the input gain depends on the timing between the input and the inhibitory action of the receiving population. Our results show that naturally an optimal phase locking for CTC emerges, and provide a mechanism by which the receiving population can implement selective communication. Moreover, the presence of bistable regions, suggests a mechanism by which different communication regimes between brain areas can be established without changing the structure of the network
### [Semiclassical Phase Reduction Theory for Quantum Synchronization](http://arxiv.org/abs/1905.05949v1) (1905.05949v1)
<i>Yuzuru Kato, Naoki Yamamoto, Hiroya Nakao</b>
<h10>2019-05-15</h10>
> We develop a general theoretical framework of semiclassical phase reduction for analyzing synchronization of quantum limit-cycle oscillators. The dynamics of quantum dissipative systems exhibiting limit-cycle oscillations are reduced to a simple, one-dimensional classical stochastic differential equation approximately describing the phase dynamics of the system under the semiclassical approximation. The density matrix and power spectrum of the original quantum system can be approximately reconstructed from the reduced phase equation. The developed framework enables us to analyze synchronization dynamics of quantum limit-cycle oscillators using the standard methods for classical limit-cycle oscillators in a quantitative way. As an example, we analyze synchronization of a quantum van der Pol oscillator under harmonic driving and squeezing, including the case that the squeezing is strong and the oscillation is asymmetric. The developed framework provides insights into the relation between quantum and classical synchronization and will facilitate systematic analysis and control of quantum nonlinear oscillators.
### [Hidden stochastic, quantum and dynamic information of Markov diffusion process and its evaluation by an entropy integral measure under the impulse controls actions, applied to information observer](http://arxiv.org/abs/1207.3091v4) (1207.3091v4)
<i>Vladimir S. Lerner</b>
<h10>2012-07-12</h10>
> Hidden information emerges under impulse interactions with Markov diffusion process modeling interactive random environment. Impulse yes no action cuts Markov correlations revealing Bit of hidden information connected correlated states. Information appears phenomenon of interaction cutting correlations carrying entropy. Each inter action models Kronicker impulse, delta impulse models interaction between the Kronicker impulses. Each impulse step down action cuts maximum of impulse minimal entropy and impulse step up action transits cutting minimal entropy to each step up action of merging delta function. Delta step down action kills delivering entropy producing equivalent minimax information. The merging action initiates quantum microprocess. Multiple cutting entropy is converting to information micro macroprocess. Cutting impulse entropy integrates entropy functional EF along trajectories of multidimensional diffusion process. Information which delivers ending states of each impulse integrates information path functional IPF along process trajectories. Hidden information evaluates Feller kernel whose minimal path transforms Markov transition probability to probability of Brownian diffusion. Each transitive transformation virtually observes origin of hidden information probabilities correlated states. IPF integrates observing Bits along minimal path assembling information Observer. Minimax imposes variation principle on EF and IPF whose extreme equations describe observing micro and macroprocess which describes irreversible thermodynamics. Hidden information curries free information frozen from correlated connections. Free information binds observing micro macro processes in information macrodynamics. Each dynamic three free information composes triplet structures. Three structural triplets assemble information network. Triple networks free information cooperate information Observer.
### [Coevolutionary dynamics of a variant of the cyclic Lotka-Volterra model with three-agent interactions](http://arxiv.org/abs/1905.05591v1) (1905.05591v1)
<i>Filippo Palombi, Stefano Ferriani, Simona Toti</b>
<h10>2019-05-14</h10>
> We study a variant of the cyclic Lotka-Volterra model with three-agent interactions. Inspired by a multiplayer variation of the Rock-Paper-Scissors game, the model describes an ideal ecosystem in which cyclic competition among three species develops through cooperative predation. Its rate equations in a well-mixed environment display a degenerate Hopf bifurcation, occurring as reactions involving two predators plus one prey have the same rate as reactions involving two preys plus one predator. We estimate the magnitude of the stochastic noise at the bifurcation point, where finite size effects turn neutrally stable orbits into erratically diverging trajectories. In particular, we compare analytic predictions for the extinction probability, derived in the Fokker-Planck approximation, with numerical simulations based on the Gillespie stochastic algorithm. We then extend the analysis of the phase portrait to heterogeneous rates. In a well-mixed environment, we observe a continuum of degenerate Hopf bifurcations, generalizing the above one. Neutral stability ensues from a complex equilibrium between different reactions. Remarkably, on a two-dimensional lattice, all bifurcations disappear as a consequence of the spatial locality of the interactions. In the second part of the paper, we investigate the effects of mobility in a lattice metapopulation model with patches hosting several agents. We find that strategies propagate along the arms of rotating spirals, as they usually do in models of cyclic dominance. We observe propagation instabilities in the regime of large wavelengths. We also examine three-agent interactions inducing nonlinear diffusion.
# <center>Chaotic Dynamics</center>
<hr>
### [On the Automatic Parameter Selection for Permutation Entropy](http://arxiv.org/abs/1905.06443v1) (1905.06443v1)
<i>Audun Myers, Firas Khasawneh</b>
<h10>2019-05-15</h10>
> Permutation Entropy (PE) has been shown to be a useful tool for time series analysis due to its low computational cost and noise robustness. This has drawn for its successful application in many fields. Some of these include damage detection, disease forecasting, and financial volatility analysis. However, to successfully use PE, an accurate selection of two parameters is needed: the permutation dimension  and embedding delay . These parameters are often suggested by experts based on a heuristic or by a trial and error approach. unfortunately, both of these methods can be time-consuming and lead to inaccurate results. To help combat this issue, in this paper we investigate multiple schemes for automatically selecting these parameters with only the corresponding time series as the input. Specifically, we develop a frequency-domain approach based on the least median of squares and the Fourier spectrum, as well as extend two existing methods: Permutation Auto-Mutual Information (PAMI) and Multi-scale Permutation Entropy (MPE) for determining . We then compare our methods as well as current methods in the literature for obtaining both  and  against expert-suggested values in published works. We show that the success of any method in automatically generating the correct PE parameters depends on the category of the studied system. Specifically, for the delay parameter , we show that our frequency approach provides accurate suggestions for periodic systems, nonlinear difference equations, and ECG/EEG data, while the mutual information function computed using adaptive partitions provides the most accurate results for chaotic differential equations. For the permutation dimension , both False Nearest Neighbors and MPE provide accurate values for  for most of the systems with  being suitable in most cases.
### [Impact of Network Topology on the Stability of DC Power Grids](http://arxiv.org/abs/1905.06429v1) (1905.06429v1)
<i>Julian F. Wienand, David Eidmann, Julian Kremers, Jobst Heitzig, Frank Hellmann, Jürgen Kurths</b>
<h10>2019-05-15</h10>
> We probe the stability of Watts-Strogatz DC power grids, in which droop-controlled producers, constant power load consumers and power lines obey Kirchhoff's circuit laws. The concept of survivability is employed to evaluate the system's response to voltage perturbations in dependence on the network topology. Following a fixed point analysis of the power grid model, we extract three main indicators of stability through numerical studies: the share of producers in the network, the node degree and the magnitude of the perturbation. Based on our findings, we investigate the local dynamics of the perturbed system and derive explicit guidelines for the design of resilient DC power grids. Depending on the imposed voltage and current limits, the stability is optimized for low node degrees or a specific share of producers.
### [Chaos in time delay systems, an educational review](http://arxiv.org/abs/1901.04826v3) (1901.04826v3)
<i>Hendrik Wernecke, Bulcsú Sándor, Claudius Gros</b>
<h10>2019-01-15</h10>
> The time needed to exchange information in the physical world induces a delay term when the respective system is modeled by differential equations. Time delays are hence ubiquitous, being furthermore likely to induce instabilities and with it various kinds of chaotic phases. Which are then the possible types of time delays, induced chaotic states, and methods suitable to characterize the resulting dynamics? This review presents an overview of the field that includes an in-depth discussion of the most important results, of the standard numerical approaches and of several novel tests for identifying chaos. Special emphasis is placed on a structured representation that is straightforward to follow. Several educational examples are included in addition as entry points to the rapidly developing field of time delay systems.
### [High harmonic generation with nearly circular polarized pulses](http://arxiv.org/abs/1905.05989v1) (1905.05989v1)
<i>Jonathan Dubois, C. Chandre, T. Uzer</b>
<h10>2019-05-15</h10>
> According to conventional wisdom, increasing ellipticity reduces high harmonic generation by several orders because the recollision probability decreases. This is the obvious conclusion drawn from the motion of an electron in a laser field without an envelope. In contrast, we report on a recollision channel with large return energy and a substantial probability, regardless the ellipticity, in which the laser envelope plays a dominant role in the energy gained by the electron, and in the conditions for the electron to come back to the core.
### [Classical three rotor problem: periodic solutions, stability and chaos](http://arxiv.org/abs/1811.05807v2) (1811.05807v2)
<i>Govind S. Krishnaswami, Himalaya Senapati</b>
<h10>2018-11-14</h10>
> This paper concerns the classical dynamics of three coupled rotors: equal masses moving on a circle subject to attractive cosine inter-particle potentials. It is a simpler variant of the gravitational three-body problem and also arises as the classical limit of a model of coupled Josephson junctions. Unlike in the gravitational problem, there are no singularities (neither collisional nor non-collisional), leading to global existence and uniqueness of solutions. In appropriate units, the non-negative energy E of the relative motion is the only free parameter. We find analogues of the Euler-Lagrange family of periodic solutions: pendulum and isosceles solutions at all energies and choreographies up to moderate energies. The model displays order-chaos-order behavior: it is integrable at zero and infinitely high energies but displays a fairly sharp transition from regular to chaotic behavior as E is increased beyond  and a more gradual return to regularity. The transition to chaos is manifested in a dramatic rise of the fraction of the area of the Hill region of Poincar\'e surfaces occupied by chaotic sections and also in the spontaneous breaking of discrete symmetries of Poincar\'e sections present at lower energies. Interestingly, the above pendulum solutions alternate between being stable and unstable, with the transition energies cascading geometrically from either side at E = 4. The transition to chaos is also reflected in the curvature of the Jacobi-Maupertuis metric that ceases to be everywhere positive when E exceeds four. Examination of Poincar\'e sections also indicates global chaos in a band of energies  slightly above this transition.
# <center>Cellular Automata And Lattice Gases</center>
<hr>
### [Lenia - Biology of Artificial Life](http://arxiv.org/abs/1812.05433v3) (1812.05433v3)
<i>Bert Wang-Chak Chan</b>
<h10>2018-12-13</h10>
> We report a new system 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. They differ from other cellular automata patterns in being geometric, metameric, fuzzy, resilient, adaptive, and rule-generic. We present basic observations of the system regarding the properties of space-time and basic settings. We provide a broad 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.
### [Escape dynamics based on bounded rationality](http://arxiv.org/abs/1902.06541v2) (1902.06541v2)
<i>Lingxiao Wang, Yin Jiang</b>
<h10>2019-02-18</h10>
> The bounded rationality plays a vital role in the collective behavior of the evacuation process. Also investigating human behavior in such an extreme situation is a continuing concern within social psychology. In this paper, we construct a cellular automaton (CA) model for the escape dynamics, and the bounded rational behavior induced by heterogeneous information is introduced. The non-trivial behavior shows in the replicator dynamics method with mean field approximation, where people's perception of the distribution of population and velocity is reduced to an average value in a certain direction. Analyzing the escape efficiency shows that under the premise of rationality, the bounded rational strategy can get higher performance. Interestingly, a quantifiable meta-stable state appears in the escape process, and the escape time is power-law dependent on system size.
### [Sandpile monomorphisms and harmonic functions](http://arxiv.org/abs/1904.12209v1) (1904.12209v1)
<i>Moritz Lang, Mikhail Shkolnikov</b>
<h10>2019-04-27</h10>
> The abelian sandpile model is a cellular automaton defined on a finite convex domain  of the standard square lattice . Its recurrent configurations form an abelian group, the sandpile group. Little is known about the structure of this group, or the relationships between sandpile groups defined on different domains. In this article, we show that the sandpile group is isomorphic to the rational-valued discrete harmonic functions which take integer-values on the boundary  of the complement of the domain, modulo the integer-valued harmonic functions on the same domain. We use this isomorphism to derive unique coordinates for every recurrent sandpile configuration as well as an alternative formula for the order of the sandpile group. Furthermore, we derive several lemmata on the existence of cyclic subgroups for sandpile groups on  square domains. Finally, we show that there exist families of monomorphisms from each sandpile group on an  domain to all sandpile groups on  domains, with . The latter result can likely be extended to sandpile groups defined on certain other, non-rectangular domains of , or on domains of higher dimensional lattices , . Our work proposes that injective limits of the sandpile group with respect to domains of certain shapes might exist.
### [A Novel Antifragility Measure Based on Satisfaction and Its Application to Random and Biological Boolean Networks](http://arxiv.org/abs/1812.06760v2) (1812.06760v2)
<i>Omar K. Pineda, Hyobin Kim, Carlos Gershenson</b>
<h10>2018-12-17</h10>
> Antifragility is a property that enhances the capability of a system in response to external perturbations. Although the concept has been applied in many areas, a practical measure of antifragility has not been developed yet. Here we propose a simply calculable measure of antifragility, based on the change of "satisfaction" before and after adding perturbations, and apply it to random Boolean networks (RBNs). Using the measure, we found that ordered RBNs are the most antifragile. Also, we demonstrated that seven biological systems are antifragile. Our measure and results can be used in various applications of Boolean networks (BNs) including creating antifragile engineering systems, identifying the genetic mechanism of antifragile biological systems, and developing new treatment strategies for various diseases.
### [Simply modified GKL density classifiers that reach consensus faster](http://arxiv.org/abs/1904.07411v1) (1904.07411v1)
<i>J. Ricardo G. Mendonça</b>
<h10>2019-04-16</h10>
> The two-state Gacs-Kurdyumov-Levin (GKL) cellular automaton has been a staple model in the study of complex systems due to its ability to classify binary arrays of symbols according to their initial density. We show that a class of modified GKL models over extended neighborhoods, but still involving only three cells at a time, achieves comparable density classification performance but in some cases reach consensus more than twice as fast. Our results suggest the time to consensus (relative to the length of the CA) as a complementary measure of density classification performance.
<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>