**Welcome To Nonlinear Sciences**

Adaptation And Self-Organizing Systems

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

### Partial Phase Cohesiveness in Networks of Kuramoto Oscillator Networks (1906.01065v1)

Yuzhen Qin, Yu Kawano, Oscar Portoles, Ming Cao

2019-06-03

Partial, instead of complete, synchronization has been widely observed in various networks including, in particular, brain networks. Motivated by data from human brain functional networks, in this technical note, we analytically show that partial synchronization can be induced by strong regional connections in coupled subnetworks of Kuramoto oscillators. To quantify the required strength of regional connections, we first obtain a critical value for the algebraic connectivity of the corresponding subnetwork using the incremental 2 norm. We then introduce the concept of the generalized complement graph, and obtain another condition on the weighted nodal degree by using the incremental infinity norm. Under these two conditions, regions of attraction for partial phase cohesiveness are estimated in the forms of the incremental 2 and infinity norms, respectively. Our result based on the incremental infinity norm is the first known criterion that is applicable to non-complete graphs. Numerical simulations are performed on a two-level network to illustrate our theoretical results; more importantly, we use real anatomical brain network data to show how our results may reveal the interplay between anatomical structure and empirical patterns of synchrony.

### Neural networks grown and self-organized by noise (1906.01039v1)

Guruprasad Raghavan, Matt Thomson

2019-06-03

Living neural networks emerge through a process of growth and self-organization that begins with a single cell and results in a brain, an organized and functional computational device. Artificial neural networks, however, rely on human-designed, hand-programmed architectures for their remarkable performance. Can we develop artificial computational devices that can grow and self-organize without human intervention? In this paper, we propose a biologically inspired developmental algorithm that can 'grow' a functional, layered neural network from a single initial cell. The algorithm organizes inter-layer connections to construct a convolutional pooling layer, a key constituent of convolutional neural networks (CNN's). Our approach is inspired by the mechanisms employed by the early visual system to wire the retina to the lateral geniculate nucleus (LGN), days before animals open their eyes. The key ingredients for robust self-organization are an emergent spontaneous spatiotemporal activity wave in the first layer and a local learning rule in the second layer that 'learns' the underlying activity pattern in the first layer. The algorithm is adaptable to a wide-range of input-layer geometries, robust to malfunctioning units in the first layer, and so can be used to successfully grow and self-organize pooling architectures of different pool-sizes and shapes. The algorithm provides a primitive procedure for constructing layered neural networks through growth and self-organization. Broadly, our work shows that biologically inspired developmental algorithms can be applied to autonomously grow functional 'brains' in-silico.

### How does active participation effect consensus: Adaptive network model of opinion dynamics and influence maximizing rewiring (1906.00868v1)

Markus Brede

2019-06-03

In this paper we study the impact of active participation -- or deliberately seeking out other agents with an aim to convince them -- on the dynamics of consensus formation. For this purpose, we propose an adaptive network model in which two processes shape opinion dynamics at interwoven time-scales: (i) agents adapt their opinions subject to influence from social network neighbours who hold opinions within a tolerance interval and (ii) agents rewire network connections with an aim to maximize their own influence on overall system opinion. We study this system in both an endogenous setting in which all agents are subject to influence and also attempt to maximize influence, and in a setting of exogenous control, in which external agents not subject to influence adaptively attempt to maximize their influence. In both settings we find three regimes of stationary opinion configurations: (i) for low a regime of two evenly balanced radicalized opinion clusters at the extremes of the opinion space, (ii) for intermediate a 'winner-takes-most' regime of two unevenly sized radicalized opinion clusters, and (iii) for large a regime in which very low spread compromise consensus states can be reached. Comparing to adaptive processes of random and deliberately spread-reducing rewiring, we demonstrate that in regime (iii) competitive influence maximization can achieve near-minimal opinion spread within near-optimal times. Further, we also show that competitive influence maximizing rewiring can reduce the impact of small influential minorities on consensus states.

Chaotic Dynamics

### Computation of kinematic and magnetic -effect and eddy diffusivity tensors by Padé approximation (1906.01423v1)

Sílvio M. A. Gama, Roman Chertovskih, Vladislav Zheligovsky

2019-06-04

We present examples of Pad'e approximation of the -effect and eddy viscosity/diffusivity tensors in various flows. Expressions for the tensors derived in the framework of the standard multiscale formalism are employed. Algebraically the simplest case is that of a two-dimensional parity-invariant six-fold rotation-symmetric flow, where eddy viscosity is negative, indicating intervals of large-scale instability of the flow. Turning to the kinematic dynamo problem for three-dimensional flows of an incompressible fluid, we explore application of Pad'e approximants for computation of tensors of magnetic -effect and, for parity-invariant flows, of magnetic eddy diffusivity. We construct Pad'e approximants of the tensors expanded in power series in the inverse molecular diffusivity around . This yields the values of the dominant growth rate due to the action of the -effect or eddy diffusivity to satisfactory accuracy for , several dozen times smaller than the threshold, above which the power series is convergent. For one sample flow, we observe eddy diffusivity tending to negative infinity when tends from above to the point of the onset of small-scale dynamo action in a symmetry-invariant subspace where a neutral small-scale magnetic mode resides. However, 49 first coefficients in the power series in prove insufficient for Pad'e approximants to reproduce this behaviour. We do computations in Fortran in the standard

`double' (real*8) and extended`

quadruple' (real*16) precision, as well as perform symbolic calculations in Mathematica.

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

### Dimensional scaling of flame propagation in discrete particulate clouds (1906.01138v1)

Fredric Y. K. Lam, XiaoCheng Mi, Andrew J. Higgins

2019-06-04

The critical dimension necessary for a flame to propagate in suspensions of fuel particles in oxidizer is studied analytically and numerically. Two types of models are considered: First, a continuum model, wherein the individual particulate sources are not resolved and the heat release is assumed spatially uniform, is solved via conventional finite difference techniques. Second, a discrete source model, wherein the heat diffusion from individual sources is modeled via superposition of the Green's function of each source, is employed to examine the influence of the random, discrete nature of the media. Heat transfer to cold, isothermal walls and to a layer of inert gas surrounding the reactive medium are considered as the loss mechanisms. Both cylindrical and rectangular (slab) geometries of the reactive medium are considered, and the flame speed is measured as a function of the diameter and thickness of the domains, respectively. In the continuum model with inert gas confinement, a universal scaling of critical diameter to critical thickness near 2:1 is found. In the discrete source model, as the time scale of heat release of the sources is made small compared to the interparticle diffusion time, the geometric scaling between cylinders and slabs exhibits values greater than 2:1. The ability of the flame in the discrete regime to propagate in thinner slabs than predicted by continuum scaling is attributed to the flame being able to exploit local fluctuations in concentration across the slab to sustain propagation. As the heat release time of the sources is increased, the discrete source model reverts back to results consistent with the continuum model. Implications of these results for experiments are discussed.

### Eigenvalue Statistics for Generalized Symmetric and Hermitian Matrices (1906.01023v1)

Adway Kumar Das, Anandamohan Ghosh

2019-06-03

The Nearest Neighbour Spacing (NNS) distribution can be computed for generalized symmetric 2x2 matrices having different variances in the diagonal and in the off-diagonal elements. Tuning the relative value of the variances we show that the distributions of the level spacings exhibit a crossover from clustering to repulsion as in GOE. The analysis is extended to 3x3 matrices where distributions of NNS as well as Ratio of Nearest Neighbour Spacing (RNNS) show similar crossovers. We show that it is possible to calculate NNS distributions for Hermitian matrices (N=2, 3) where also crossovers take place between clustering and repulsion as in GUE. For large symmetric and Hermitian matrices we use interpolation between clustered and repulsive regimes and identify phase diagrams with respect to the variances.

### Many-body Chaos in a Thermalised Fluid (1906.00016v1)

Dheeraj Kumar, Subhro Bhattacharjee, Samriddhi Sankar Ray

2019-05-31

We use a new measure of many-body chaos for classical systems---cross-correlators---to show that in a thermalised fluid (obtained from a non-linear, prototypical equation of hydrodynamics sharing formal similarities with models of turbulence) characterised by a temperature and degrees of freedom, the Lyapunov exponent scales as . This bound, obtained from detailed numerical simulations and theoretical estimates, provides compelling evidence not only for recent conjectures for chaotic, equilibrium, classical many-body systems, as well as, numerical results from frustrated spin systems, but also, remarkably, show that scales linearly with the degrees of freedom in a finite-dimensional, classical, chaotic system.

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.

Don't forget to Follow and Resteem. @complexcity

Keeping everyone inform.