Welcome To Nonlinear Sciences

Adaptation And Self-Organizing Systems

---

A method for the classification of chimera states of coupled oscillators and its application for creating a neural network information converter (1906.02680v1)

Andrei Velichko

2019-06-06

The paper presents a new method for the classification of chimera states, which characterizes the synchronization of two coupled oscillators more accurately. As an example of method application, a neural network information converter based on a network of pulsed oscillators is used, which can convert input information from digital to analogue type and perform information processing after training the network by selecting control parameters. In the proposed neural network scheme, the data arrives at the input layer in the form of current levels of the oscillators and is converted into a set of non-repeating states of the chimera synchronization of the output oscillator. By modelling a thermally coupled VO2-oscillator circuit, the network setup is demonstrated through the selection of coupling strength, power supply levels and the synchronization efficiency parameter. The distribution of solutions depending on the operating mode of the oscillators, prethreshold mode or generation mode are revealed. Technological approaches for the implementation of a neural network information converter are proposed, and examples of its application for image filtering are demonstrated. The proposed method for the classification of chimera states helps significantly expand the capabilities of neuromorphic and logical devices based on synchronization effects.

Onset of synchronization in networks of second-order Kuramoto oscillators with delayed coupling: Exact results and application to phase-locked loops (1906.02643v1)

David Métivier, Lucas Wetzel, Shamik Gupta

2019-06-06

We consider the inertial Kuramoto model of globally coupled oscillators characterized by both their phase and angular velocity, in which there is a time delay in the interaction between the oscillators. Besides the academic interest, we show that the model can be related to a network of phase-locked loops widely used in electronic circuits for generating a stable frequency at multiples of an input frequency. We study the model for a generic choice of the natural frequency distribution of the oscillators, to elucidate how a synchronized phase bifurcates from an incoherent phase as the coupling constant between the oscillators is tuned. We show that in contrast to the case with no delay, here the system in the stationary state may exhibit either a subcritical or a supercritical bifurcation between a synchronized and an incoherent phase, which is dictated by the value of the delay present in the interaction and the precise value of inertia of the oscillators. Our theoretical analysis, performed in the limit , is based on an unstable manifold expansion in the vicinity of the bifurcation, which we apply to the kinetic equation satisfied by the single-oscillator distribution function. We check our results by performing direct numerical integration of the dynamics for large , and highlight the subtleties arising from having a finite number of oscillators.

Universal and non-universal neural dynamics on small world connectomes: a finite size scaling analysis (1905.05280v2)

Mahdi Zarepour, Juan I. Perotti, Orlando V. Billoni, Dante R. Chialvo, Sergio A. Cannas

2019-05-13

Evidence of critical dynamics has been recently found in both experiments and models of large scale brain dynamics. The understanding of the nature and features of such critical regime is hampered by the relatively small size of the available connectome, which prevent among other things to determine its associated universality class. To circumvent that, here we study a neural model defined on a class of small-world network that share some topological features with the human connectome. We found that varying the topological parameters can give rise to a scale-invariant behavior belonging either to mean field percolation universality class or having non universal critical exponents. In addition, we found certain regions of the topological parameters space where the system presents a discontinuous (i.e., non critical) dynamical phase transition into a percolated state. Overall these results shed light on the interplay of dynamical and topological roots of the complex brain dynamics.

Co-Contagion Diffusion on Multilayer Networks (1903.06327v2)

Ho-Chun Herbert Chang, Feng Fu

2019-03-15

This study examines the interface of three elements during co-contagion diffusion: the \textbf{synergy} between contagions, the \textbf{dormancy} rate of each individual contagion, and the \textbf{multiplex network topology}. Dormancy is defined as a weaker form of "immunity," where dormant nodes no longer actively participate in diffusion, but are still susceptible to infection. The proposed model extends the literature on threshold models, and demonstrates intricate interdependencies between different graph structures. Our simulations show that first, the faster contagion induces branching on the slower contagion; second, shorter characteristic path lengths diminish the impact of dormancy in lowering diffusion. Third, when two long-range graphs are paired, the faster contagion depends on both dormancy rates, whereas the slower contagion depends only on its own; fourth, synergistic contagions are less sensitive to dormancy, and have a wider window to diffuse. Furthermore, when long-range and spatially constrained graphs are paired, ring vaccination occurs on the spatial graph and produces partial diffusion, due to dormant, surrounding nodes. The spatial contagion depends on both dormancy rates whereas the long-range contagion depends on only its own.

Bifurcation analysis of a TaO memristor model (1906.01377v1)

Y. V. Pershin, V. A. Slipko

2019-06-04

This paper presents a study of bifurcation in the time-averaged dynamics of TaO memristors driven by narrow pulses of alternating polarities. The analysis, based on a physics-inspired model, focuses on the stable fixed points and on how these are affected by the pulse parameters. Our main finding is the identification of a driving regime when two stable fixed points exist simultaneously. To the best of our knowledge, such bistability is identified in a single memristor for the first time. This result can be readily tested experimentally, and is expected to be useful in future memristor circuit designs.

Chaotic Dynamics

---

The route to chaos in routing games: Population increase drives period-doubling instability, chaos & inefficiency with Price of Anarchy equal to one (1906.02486v1)

Thiparat Chotibut, Fryderyk Falniowski, Michał Misiurewicz, Georgios Piliouras

2019-06-06

We study a learning dynamic model of routing (congestion) games to explore how an increase in the total demand influences system performance. We focus on non-atomic routing games with two parallel edges of linear cost, where all agents evolve using Multiplicative Weights Updates with a fixed learning rate. Previous game-theoretic equilibrium analysis suggests that system performance is close to optimal in the large population limit, as seen by the Price of Anarchy reduction. In this work, however, we reveal a rather undesirable consequence of non-equilibrium phenomena driven by population increase. As the total demand rises, we prove that the learning dynamics unavoidably become non-equilibrating, typically chaotic. The Price of Anarchy predictions of near-optimal performance no longer apply. To the contrary, the time-average social cost may converge to its worst possible value in the large population limit. Every system has a carrying capacity, above which the dynamics is non-equilibrating. If the equilibrium flow is a symmetric split, the system exhibits one period-doubling bifurcation. A single periodic attractor of period two replaces the attracting fixed point when the demand exceeds the carrying capacity. In general, for asymmetric equilibrium flows, increasing the demand destabilizes the system, so that the system eventually becomes Li-Yorke chaotic with positive topological entropy. This demand-driven instability emerges from any pair of linear cost functions. Remarkably, in any non-equilibrating regime, the time-average flows on the edges converge {\it exactly} to the equilibrium flows, a property akin to no-regret learning in zero-sum games. Our results extend to any sequence of shrinking learning rates, e.g., , by allowing for a dynamically increasing population size.

Bifurcation without parameters in a chaotic system with a memristive element (1906.02445v1)

Tom Birkoben, Moritz Drangmeister, Finn Zahari, Serhiy Yanchuk, Philipp Hövel, Hermann Kohlstedt

2019-06-06

We investigate the effect of memory on a chaotic system experimentally and theoretically. For this purpose, we use Chua's oscillator as an electrical model system showing chaotic dynamics extended by a memory element in form of a double-barrier memristive device. The device consists of Au/NbO/AlO/Al/Nb layers and exhibits strong analog-type resistive changes depending on the history of the charge flow. In the extended system strong changes in the dynamics of chaotic oscillations are observable. The otherwise fluctuating amplitudes of the Chua system are disrupted by transient silent states. After developing a model for Chua's oscillator with a memristive device, the numerical treatment reveals the underling dynamics as driven by the slow-fast dynamics of the memory element. Furthermore, the stabilizing and destabilizing dynamic bifurcations are identified that are passed by the system during its chaotic behavior.

Building a Maxey--Riley framework for surface ocean inertial particle dynamics (1906.01080v2)

F. J. Beron-Vera, M. J. Olascoaga, P. Miron

2019-05-29

A Maxey-Riley set for surface ocean inertial (i.e., buoyant, finite-size) particle dynamics is derived by vertically integrating the original Maxey-Riley set, adapted to account for Earth's rotation and sphericity effects, across a sufficiently small spherical particle, floating at an unperturbed air-sea interface with unsteady nonuniform winds and ocean currents above and below, respectively. The inertial particle velocity is shown to exponentially decay in time to a velocity that lies close to an average of seawater and air velocities, weighted by a function of the seawater-to-particle density ratio. Such a weighted average velocity turns out to fortuitously be of the type commonly discussed in the search-and-rescue literature, which alone cannot explain the observed role of anticyclonic mesoscale eddies as traps for marine debris or the formation of great garbage patches in the subtropical gyres, phenomena dominated by finite-size effects. A heuristic extension of the theory is proposed to describe the motion of nonspherical particles by means of a simple shape factor correction, and recommendations are made for incorporating wave-induced Stokes drift, consistently accounting for memory effects in the presence of recurrent motions, and allowing for inhomogeneities of the carrying fluid density. The new Maxey--Riley set outperforms an ocean adaptation that ignored wind drag effects and the first reported adaption that attempted to incorporate them.

-body chaos and the continuum limit in numerical simulations of self-gravitating systems, revisited (1901.08981v2)

Pierfrancesco Di Cintio, Lapo Casetti

2019-01-25

We revisit the r^{o}le of discreteness and chaos in the dynamics of self-gravitating systems by means of -body simulations with active and frozen potentials, starting from spherically symmetric stationary states and considering the orbits of single particles in a frozen -body potential as well as the orbits of the system in the full -dimensional phase space. We also consider the intermediate case where a test particle moves in the field generated by non-interacting particles, which in turn move in a static smooth potential. We investigate the dependence on and on the softening length of the largest Lyapunov exponent both of single particle orbits and of the full -body system. For single orbits we also study the dependence on the angular momentum and on the energy. Our results confirm the expectation that orbital properties of single orbits in finite- systems approach those of orbits in smooth potentials in the continuum limit and that the largest Lyapunov exponent of the full -body system does decrease with , for sufficiently large systems. However, single orbits in frozen models and active self-consistent models have different largest Lyapunov exponents and the -dependence of the values in non-trivial, so that the use of frozen -body potentials to gain information on large- systems or on the continuum limit may be misleading in certain cases.

Directed momentum current induced by the PT-symmetric driving (1901.01699v2)

Wen-Lei Zhao, Jiaozi Wang, Xiaohui Wang, Peiqing Tong

2019-01-07

We investigate the directed momentum current in the quantum kicked rotor model with symmetric deriving potential. For the quantum non-resonance case, the values of quasi-energy become to be complex when the strength of imaginary part of the kicking potential exceeds \textbf{a} threshold value, which demonstrates the appearance of the spontaneous symmetry breaking. In the vicinity of the phase-transition point, the momentum current exhibits a staircase growth with time. Each platform of the momentum current corresponds to the mean momentum of some eigenstates of the Floquet operator whose imaginary parts of the quasi-energy are significantly large. Above the phase-transition point, the momentum current increases linearly with time. Interestingly, its acceleration rate exhibits a kind of "quantized" increment with the kicking strength. We propose a modified classical acceleration mode of the kicked rotor model to explain such an intriguing phenomenon. Our theoretical prediction is in good agreement with numerical results.

Cellular Automata And Lattice Gases

---

Paradox of integration---Cellular automata approach (1906.00393v1)

Krzysztof Malarz, Krzysztof Kułakowski

2019-06-02

We discuss the self-deprecating strategy introduced by Peter Blau as one of stages of the process of social integration. Recently we have introduced a two-dimensional space of status, real and surface one ( and ), and we have demonstrated that with this setup, the self-deprecating strategy efficiently prevents the rejection (Malarz and Ku{\l}akowski, arXiv:1903.04291). There, the process of reducing the conflict was described by master equations, i.e. a set of differential equations describing evolution of density of actors of status . Here we reformulate the problem in terms of probabilistic asynchronous cellular automata. The obtained results for number of actors of status are qualitatively the same as in the previous approach. Namely, an enhancement of the surface status compensates a deficiency of the real one. Yet, the results depend on the order of application of the automaton rules; the above results are obtained for the backward typewriter order.

Coarse Graining of Partitioned Cellular Automata (1905.10391v1)

Pedro C. S. Costa, Fernando de Melo

2019-05-24

Partitioned cellular automata are known to be an useful tool to simulate linear and nonlinear problems in physics, specially because they allow for a straightforward way to define conserved quantities and reversible dynamics. Here we show how to construct a local coarse graining description of partitioned cellular automata. By making use of this tool we investigate the effective dynamics in this model of computation. All examples explored are in the scenario of lattice gases, so that the information lost after the coarse graining is related to the number of particles. It becomes apparent how difficult it is to remain with a deterministic dynamics after coarse graining. Several examples are shown where an effective stochastic dynamics is obtained after a deterministic dynamics is coarse grained. These results suggest why random processes are so common in nature. Although all the cases presented assume one-dimensional lattices, we show how our approach can be extended to higher dimensions.

Simply modified GKL density classifiers that reach consensus faster (1904.07411v2)

J. Ricardo G. Mendonça

2019-04-16

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.

A trust model for spreading gossip in social networks (1905.11204v1)

Rinni Bhansali, Laura P. Schaposnik

2019-05-23

We introduce here a multi-type bootstrap percolation model, which we call T-Bootstrap Percolation (T-BP), and apply it to study information propagation in social networks. In this model, a social network is represented by a graph G whose vertices have different labels corresponding to the type of role the person plays in the network (e.g. a student, an educator, etc.). Once an initial set of vertices of G is randomly selected to be carrying a gossip (e.g. to be infected), the gossip propagates to a new vertex provided it is transmitted by a minimum threshold of vertices with different labels. By considering random graphs, which have been shown to closely represent social networks, we study different properties of the T-BP model through numerical simulations, and describe its implications when applied to rumour spread, fake news, and marketing strategies.

Two-species hardcore reversible cellular automaton: matrix ansatz for dynamics and nonequilibrium stationary state (1903.10590v2)

Marko Medenjak, Vladislav Popkov, Tomaž Prosen, Eric Ragoucy, Matthieu Vanicat

2019-03-25

In this paper we study the statistical properties of a reversible cellular automaton in two out-of-equilibrium settings. In the first part we consider two instances of the initial value problem, corresponding to the inhomogeneous quench and the local quench. Our main result is an exact matrix product expression of the time evolution of the probability distribution, which we use to determine the time evolution of the density profiles analytically. In the second part we study the model on a finite lattice coupled with stochastic boundaries. Once again we derive an exact matrix product expression of the stationary distribution, as well as the particle current and density profiles in the stationary state. The exact expressions reveal the existence of different phases with either ballistic or diffusive transport depending on the boundary parameters.

---

Thank you for reading!

Don't forget to Follow and Resteem. @complexcity
Keeping everyone inform.