Welcome To Nonlinear Sciences

Adaptation And Self-Organizing Systems

---

Monte Carlo gPC methods for diffusive kinetic flocking models with uncertainties (1902.04518v1)

Jose Antonio Carrillo, Mattia Zanella

2019-02-12

In this paper we introduce and discuss numerical schemes for the approximation of kinetic equations for flocking behavior with phase transitions that incorporate uncertain quantities. This class of schemes here considered make use of a Monte Carlo approach in the phase space coupled with a stochastic Galerkin expansion in the random space. The proposed methods naturally preserve the positivity of the statistical moments of the solution and are capable to achieve high accuracy in the random space. Several tests on a kinetic alignment model with self propulsion validate the proposed methods both in the homogeneous and inhomogeneous setting, shading light on the influence of uncertainties in phase transition phenomena driven by noise such as their smoothing and confidence bands.

Classifying attention deficit hyperactivity disorder in children with non-linearities in actigraphy (1902.03530v1)

Jeremi K. Ochab, Krzysztof Gerc, Magdalena Fafrowicz, Ewa Gudowska-Nowak, Tadeusz Marek, Maciej A. Nowak, Halszka Oginska, Katarzyna Oles, Małgorzata Rams, Jerzy Szwed, Dante R. Chialvo

2019-02-10

Objective This study provides an objective measure based on actigraphy for Attention Deficit Hyperactivity Disorder (ADHD) diagnosis in children. We search for motor activity features that could allow further investigation into their association with other neurophysiological disordered traits. Method The study involved (48 eligible) male participants aged years (8 controls, and 7 in each group: ADHD combined subtype, ADHD hyperactive-impulsive subtype, and autism spectrum disorder, ASD) wearing a wristwatch actigraph continuously for a week ( losses in daily records) in two acquisition modes. We analyzed 47 quantities: from sleep duration or movement intensity to theory-driven scaling exponents or non-linear prediction errors of both diurnal and nocturnal activity. We used them in supervised classification to obtain cross-validated diagnostic performance. Results We report the best performing measures, including a nearest neighbors 4-feature classifier providing accuracy, sensitivity and specificity in a binary ADHD vs control classification and accuracy (against baseline), sensitivity and specificity in 4-class task (two ADHD subtypes, ASD, and control). The most informative feature is skewness of the shape of Zero Crossing Mode (ZCM) activity. Mean and standard deviation of nocturnal activity are among the least informative. Conclusion Actigraphy causes only minor discomfort to the subjects and is inexpensive. The range of existing mathematical and machine learning tools also allow it to be a useful add-on test for ADHD or differential diagnosis between ADHD subtypes. The study was limited to a small, male sample without the inattentive ADHD subtype.

Topological Control of Synchronization Patterns: Trading Symmetry for Stability (1902.03255v1)

Joseph D. Hart, Yuanzhao Zhang, Rajarshi Roy, Adilson E. Motter

2019-02-08

Symmetries are ubiquitous in network systems and have profound impacts on the observable dynamics. At the most fundamental level, many synchronization patterns are induced by underlying network symmetry, and a high degree of symmetry is believed to enhance the stability of identical synchronization. Yet, here we show that the synchronizability of almost any symmetry cluster in a network of identical nodes can be enhanced precisely by breaking its structural symmetry. This counterintuitive effect holds for generic node dynamics and arbitrary network structure and is, moreover, robust against noise and imperfections typical of real systems, which we demonstrate by implementing a state-of-the-art optoelectronic experiment. These results lead to new possibilities for the topological control of synchronization patterns, which we substantiate by presenting an algorithm that optimizes the structure of individual clusters under various constraints.

Spatial eco-evolutionary feedbacks mediate coexistence in prey-predator systems (1902.03016v1)

Eduardo H. Colombo, Ricardo Martínez-García, Cristóbal, López, Emilio Hernández-García

2019-02-08

Eco-evolutionary frameworks can explain certain features of communities in which ecological and evolutionary processes occur over comparable timescales. In the particular case of prey-predator systems, a combination of empirical and theoretical studies have explored this possibility, showing that the evolution of prey traits, predator traits or the coevolution of both can contribute to the stability of the community, as well as to the emergence of various types of population cycles. However, these studies overlook that interactions are spatially constrained, a crucial ingredient known to foster species coexistence per se. Here, we investigate whether evolutionary dynamics interacts with the spatial structure of a prey-predator community in which both species show limited mobility and predators perceptual ranges are subject to natural selection. In these conditions, our results unveil an eco-evolutionary feedback between species spatial mixing and predators perceptual range: different levels of species mixing select for different perceptual ranges, which in turn reshape the spatial distribution of preys and their interaction with predators. This emergent pattern of interspecific interactions feeds back to the efficiency of the various perceptual ranges, thus selecting for new ones. Finally, since prey-predator mixing is the key factor that regulates the intensity of predation, we explore the community-level implications of such feedback and show that it controls both coexistence times and species extinction probabilities.

A biologically motivated three-species exclusion model: effects of leaky scanning and overlapping genes on initiation of protein synthesis (1812.10143v2)

Bhavya Mishra, Debashish Chowdhury

2018-12-25

Totally asymmetric simple exclusion process (TASEP) was originally introduced as a model for the traffic-like collective movement of ribosomes on a messenger RNA (mRNA) that serves as the track for the motor-like forward stepping of individual ribosomes. In each step, a ribosome elongates a protein by a single unit using the track also as a template for protein synthesis. But, pre-fabricated, functionally competent, ribosomes are not available to begin synthesis of protein; a subunit directionally scans the mRNA in search of the pre-designated site where it is supposed to bind with the other subunit and begin the synthesis of the corresponding protein. However, because of `leaky' scanning, a fraction of the scanning subunits miss the target site and continue their search beyond the first target. Sometimes such scanners successfully identify the site that marks the site for initiation of the synthesis of a different protein. In this paper, we develop an exclusion model, with three interconvertible species of hard rods, to capture some of the key features of these biological phenomena and study the effects of the interference of the flow of the different species of rods on the same lattice. More specifically, we identify the meantime for the initiation of protein synthesis as appropriate mean {\it first-passage} time that we calculate analytically using the formalism of backward master equations. In spite of the approximations made, our analytical predictions are in reasonably good agreement with the numerical data that we obtain by performing Monte Carlo simulations. We also compare our results with a few experimental facts reported in the literature and propose new experiments for testing some of our new quantitative predictions.

Chaotic Dynamics

---

Negative magnetic eddy diffusivity due to oscillogenic -effect (1711.02390v3)

Alexander Andrievsky, Roman Chertovskih, Vladislav Zheligovsky

2017-11-07

We study large-scale kinematic dynamo action of steady mirror-antisymmetric flows of incompressible fluid, that involve small spatial scales only, by asymptotic methods of the multiscale stability theory. It turns out that, due to the magnetic -effect in such flows, the large-scale mean field experiences harmonic oscillations in time on the scale O() without growth or decay. Here is the spatial scale ratio and is the fast time of the order of the flow turnover time. The interaction of the accompanying fluctuating magnetic field with the flow gives rise to an anisotropic magnetic eddy diffusivity, whose dependence on the direction of the large-scale wave vector generically exhibits a singular behaviour, and thus to negative eddy diffusivity for whichever molecular magnetic diffusivity. Consequently, such flows always act as kinematic dynamos on the time scale O(); for the directions at which eddy diffusivity is infinite, the large-scale mean-field growth rate is finite on the scale O(). We investigate numerically this dynamo mechanism for two sample flows.

Electro-Oxidation of p-Silicon in Fluoride-Containing Electrolyte: A Physical Model for the Regime of Negative Differential Resistance (1902.04328v1)

Munir Salman, Maximilian Patzauer, Dominique Koster, Fabio La Mantia, Katharina Krischer

2019-02-12

When Si is anodically oxidized in a fluoride containing electrolyte, an oxide layer is grown. Simultaneously, the layer is etched by the fluoride containing electrolyte. The resulting stationary state exhibits a negative slope of the current-voltage characteristics in a certain range of applied voltage. We propose a physical model that reproduces this negative slope. In particular, our model assumes that the oxide layer consists of both partially and fully oxidized Si and that the etch rate depends on the effective degree of oxidation. Finally, we show that our simulations are in good agreement with measurements of the current-voltage characteristics, the oxide layer thickness, the dissolution valence, and the impedance spectra of the electrochemical system.

Surface waves enhance particle dispersion (1902.04034v1)

Mohammad Farazmand, Themistoklis Sapsis

2019-02-11

We study the horizontal dispersion of passive tracer particles on the free surface of gravity waves in deep water. For random linear waves with the JONSWAP spectrum, the Lagrangian particle trajectories are computed using an exact nonlinear model known as the John--Sclavounos equation. We show that the single-particle dispersion exhibits an unusual super-diffusive behavior. In particular, for large times , the variance of the tracer increases as a quadratic function of time, i.e., . This dispersion is markedly faster than Taylor's single-particle dispersion theory which predicts that the variance of passive tracers grows linearly with time for large . Our results imply that the wave motion significantly enhances the dispersion of fluid particles. We show that this super-diffusive behavior is a result of the long-term correlation of the Lagrangian velocities of fluid parcels on the free surface.

Stochastic bistable systems, and competing hysteresis and phase coexistence (1902.03991v1)

Mahendra K. Verma, Abhishek Kumar, Adhip Pattanayak

2019-02-11

In this paper we describe the solution of a stochastic bistable system from a dynamical perspective. We show how a single framework with variable noise can explain hysteresis at zero temperature and two-state coexistence in the presence of noise. This feature is similar to the phase transition of thermodynamics. Our mathematical model for bistable systems also explains how the width of a hysteresis loop shrinks in the presence of noise, and how variation in initial conditions can take such systems to different final states.

Theory of coupled parametric oscillators beyond coupled Ising spins (1901.07372v2)

Marcello Calvanese Strinati, Leon Bello, Avi Pe'er, Emanuele G. Dalla Torre

2019-01-18

Periodically driven parametric oscillators offer a convenient way to simulate classical Ising spins. When many parametric oscillators are coupled dissipatively, they can be analogous to networks of Ising spins, forming an effective coherent Ising machine (CIM) that efficiently solves computationally hard optimization problems. In the companion letter, we studied experimentally the minimal realization of a CIM, i.e. two coupled parametric oscillators. We found that the presence of an energy-conserving coupling between the oscillators can dramatically change the dynamics, leading to everlasting beats, which transcend the Ising description. Here, we analyze this effect theoretically by solving numerically and, when possible, analytically the equations of motion of two parametric oscillators. Our main tools include: (i) a Floquet analysis of the linear equations, (ii) a multi-scale analysis based on a separation of scales between the parametric oscillations and the beats, and (iii) the numerical identification of limit cycles and attractors. Using these tools, we fully determine the phase boundaries and critical exponents of the model, as a function of the intensity and the phase of the coupling and of the pump. Our study highlights the universal character of the phase diagram and its independence on the specific type of nonlinearity present in the system. Furthermore, we identify new phases of the model with more than two attractors, possibly describing a larger spin algebra.

Cellular Automata And Lattice Gases

---

A Survey of the Individual-Based Model applied in Biomedical and Epidemiology (1902.02784v1)

E. G. Nepomuceno, D. F. Resende, M. J. Lacerda

2019-02-07

Individual-based model (IBM) has been used to simulate and to design control strategies for dynamic systems that are subject to stochasticity and heterogeneity, such as infectious diseases. In the IBM, an individual is represented by a set of specific characteristics that may change dynamically over time. This feature allows a more realistic analysis of the spread of an epidemic. This paper presents a literature survey of IBM applied to biomedical and epidemiology research. The main goal is to present existing techniques, advantages and future perspectives in the development of the model. We evaluated 89 articles, which mostly analyze interventions aimed at endemic infections. In addition to the review, an overview of IBM is presented as an alternative to complement or replace compartmental models, such as the SIR (Susceptible-Infected-Recovered) model. Numerical simulations also illustrate the capabilities of IBM, as well as some limitations regarding the effects of discretization. We show that similar side-effects of discretization scheme for compartmental models may also occur in IBM, which requires careful attention.

Ergodicity versus non-ergodicity for Probabilistic Cellular Automata on rooted trees (1710.00084v2)

Bruno Kimura, Wioletta Ruszel, Cristian Spitoni

2017-09-29

In this article we study a class of shift-invariant and positive rate probabilistic cellular automata (PCA) on rooted d-regular trees . In a first result we extend the results of [10] on trees, namely we prove that to every stationary measure of the PCA we can associate a space-time Gibbs measure on . Under certain assumptions on the dynamics the converse is also true. A second result concerns proving sufficient conditions for ergodicity and non-ergodicity of our PCA on d-ary trees for and characterizing the invariant product Bernoulli measures.

Throughput reduction by the simultaneous effect of two traveling routes equipped with parking sites in a two-dimensional stochastic lattice model (1901.10390v1)

Satori Tsuzuki, Daichi Yanagisawa, Katsuhiro Nishinari

2019-01-29

This paper proposes a two-dimensional stochastic lattice model that comprises a junction of two traveling routes: the domestic route and the international route, each of which has parking sites. In our model, the system distributes the arrived particles to either of the two routes and selects one of the parking sites in the route for each particle, which stops at the parking site once during its travel. Because each particle has antennas in the back and front directions to detect other approaching particles, the effect of the volume exclusion of each particle extends in the moving direction. The system displays interesting behavior; remarkably, the dependence of the throughput on the distribution ratio of particles to the domestic route reduces after reaching the maximum parking capacity of the domestic route. Our simulations and analysis with the queueing model describe this phenomenon and suggest the following fact: As the distribution ratio of particles to the international route decreases, the throughput of the international route reduces, and simultaneously, that of the domestic route saturates. The simultaneous effect of the decrease and saturation causes reduction in the throughput of the entire system.

Dissipation of traffic congestion using agent-based car-following model with modified optimal velocity (1901.06058v1)

Manit Klawtanong, Surachate Limkumnerd

2019-01-18

We investigate dynamical properties of traffic flow using the stochastic car-following model with modified optimal velocity on circular road. The safety distance following the two-second rule and autonomous vehicles, acting as agents, obeying simple requirements are incorporated into the model. The dynamic safety distance increases in a light traffic condition where the average driving velocity is high, while decreases in a dense traffic condition in anticipation of slower traffic motion. The results show that the presence of the agents can enhance overall velocity and traffic current of the system, and postpone the traffic congestion. In a particular phase region, imposing a speed limit enables the system to leave the congested flow phase. The density-dependent speed limit in agent-free condition is obtained to achieve the optimal traffic flow.

Search Space Reduction of Asynchrony Immune Cellular Automata by Center Permutivity (1901.01534v1)

Luca Mariot, Luca Manzoni, Alberto Dennunzio

2019-01-06

We continue the study of asynchrony immunity in cellular automata (CA), which can be considered as a weaker version of correlation immunity in the context of vectorial Boolean functions. The property could have applications as a countermeasure for side-channel attacks in CA-based cryptographic primitives, such as S-boxes and pseudorandom number generators. We first give some theoretical results on the necessary conditions that a CA rule must satisfy in order to meet asynchrony immunity, the most important one being center permutivity. Next, we perform an exhaustive search of all asynchrony immune CA rules of neighborhood size up to , leveraging on the discovered theoretical properties to greatly reduce the size of the search space.

Pattern Formation And Solitons

---

Electro-Oxidation of p-Silicon in Fluoride-Containing Electrolyte: A Physical Model for the Regime of Negative Differential Resistance (1902.04328v1)

Munir Salman, Maximilian Patzauer, Dominique Koster, Fabio La Mantia, Katharina Krischer

2019-02-12

When Si is anodically oxidized in a fluoride containing electrolyte, an oxide layer is grown. Simultaneously, the layer is etched by the fluoride containing electrolyte. The resulting stationary state exhibits a negative slope of the current-voltage characteristics in a certain range of applied voltage. We propose a physical model that reproduces this negative slope. In particular, our model assumes that the oxide layer consists of both partially and fully oxidized Si and that the etch rate depends on the effective degree of oxidation. Finally, we show that our simulations are in good agreement with measurements of the current-voltage characteristics, the oxide layer thickness, the dissolution valence, and the impedance spectra of the electrochemical system.

Dark solitons for an extended quintic nonlinear Schrödinger equation: Application to water waves at (1810.12500v2)

F. Tsitoura, T. P. Horikis, D. J. Frantzeskakis

2018-10-30

We study the existence, formation and dynamics of gray solitons for an extended quintic nonlinear Schr"odinger (NLS) equation. The considered model finds applications to water waves, when the characteristic parameter - where is the wavenumber and is the undistorted water's depth - takes the critical value . It is shown that this model admits approximate dark soliton solutions emerging from an effective Korteweg-de Vries equation and that two types of gray solitons exist: fast and slow, with the latter being almost stationary objects. Analytical results are corroborated by direct numerical simulations.

Drift of spectrally stable shifted states on star graphs (1902.03612v1)

Adilbek Kairzhan, Dmitry E. Pelinovsky, Roy H. Goodman

2019-02-10

When the coefficients of the cubic terms match the coefficients in the boundary conditions at a vertex of a star graph and satisfy a certain constraint, the nonlinear Schr"{o}dinger (NLS) equation on the star graph can be transformed to the NLS equation on a real line. Such balanced star graphs have appeared in the context of reflectionless transmission of solitary waves. Steady states on such balanced star graphs can be translated along the edges with a translational parameter and are referred to as the shifted states. When the star graph has exactly one incoming edge and several outgoing edges, the steady states are spectrally stable if their monotonic tails are located on the outgoing edges. These spectrally stable states are degenerate minimizers of the action functional with the degeneracy due to the translational symmetry. Nonlinear stability of these spectrally stable states has been an open problem up to now. In this paper, we prove that these spectrally stable states are nonlinearly unstable because of the irreversible drift along the incoming edge towards the vertex of the star graph. When the shifted states reach the vertex as a result of the drift, they become saddle points of the action functional, in which case the nonlinear instability leads to their destruction. In addition to rigorous mathematical results, we use numerical simulations to illustrate the drift instability and destruction of the shifted states on the balanced star graph.

Salt Polygons are Caused by Convection (1902.03600v1)

Jana Lasser, Joanna M. Nield, Marcel Ernst, Volker Karius, Lucas Goehring

2019-02-10

From fairy circles to patterned ground and columnar joints, natural patterns spontaneously appear in many complex geophysical settings. Here, we consider the origins of polygons in the crusts of salt playa and salt pans. These beautifully regular features, a meter or so in diameter, have a worldwide distribution and are important to the transport of salt and dust in arid regions, yet there has been no convincing mechanism known for their formation. We present the first evidence that they are the surface expression of buoyancy-driven convection in the porous soil beneath a salt crust. By combining consistent results from direct field observations, analogue experiments, linear stability theory, and numerical simulations, we further explain the conditions under which salt polygons will form, as well as how their characteristic size emerges.

Cesium bright matter-wave solitons and soliton trains (1902.03144v1)

Tadej Mežnaršič, Tina Arh, Jure Brence, Jaka Pišljar, Katja Gosar, Žiga Gosar, Rok Žitko, Erik Zupanič, Peter Jeglič

2019-02-08

A setup for studying bright matter-wave solitons of a cesium Bose-Einstein condensate (BEC) is presented. Production of a single soliton is demonstrated and dependence of soliton atom number on the interatomic interaction is studied. Formation of soliton trains in the quasi one-dimensional confinement is shown. Additionally, trains forming from a BEC have been observed outside confinement, in free space. In the end a double BEC production setup for studying soliton collisions is described.

Exactly Solvable And Integrable Systems

---

New determinants in the 8VSOS model with domain-wall boundaries (1812.10781v2)

W. Galleas

2018-12-27

In this letter we show the partition function of the 8VSOS model with domain-wall boundaries satisfies the same type of functional equations as its six-vertex model counterpart. We then use these refined functional equations to obtain novel determinantal representations for the aforementioned partition function.

Global description of action-angle duality for a Poisson-Lie deformation of the trigonometric Sutherland system (1710.08760v2)

L. Feher, I. Marshall

2017-10-24

Integrable many-body systems of Ruijsenaars--Schneider--van Diejen type displaying action-angle duality are derived by Hamiltonian reduction of the Heisenberg double of the Poisson-Lie group . New global models of the reduced phase space are described, revealing non-trivial features of the two systems in duality with one another. For example, after establishing that the symplectic vector space underlies both global models, it is seen that for both systems the action variables generate the standard torus action on , and the fixed point of this action corresponds to the unique equilibrium positions of the pertinent systems. The systems in duality are found to be non-degenerate in the sense that the functional dimension of the Poisson algebra of their conserved quantities is equal to half the dimension of the phase space. The dual of the deformed Sutherland system is shown to be a limiting case of a van Diejen system.

Scaling laws for harmonically trapped two-species mixtures at thermal equilibrium (1902.04007v1)

Francisco Jauffred, Roberto Onofrio, Bala Sundaram

2019-02-11

We discuss the scaling of the interaction energy with particle numbers for a harmonically trapped two-species mixture at thermal equilibrium experiencing interactions of arbitrary strength and range. In the limit of long-range interactions and weak coupling, we recover known results for the integrable Caldeira-Leggett model in the classical limit. In the case of short-range interactions and for a balanced mixture, numerical simulations show scaling laws with exponents that depend on the interaction strength, its attractive or repulsive nature, and the dimensionality of the system. Simple analytic considerations based on equilibrium statistical mechanics and small interspecies coupling quantitatively recover the numerical results. The dependence of the scaling on interaction strength helps to identify a threshold between two distinct regimes. Our thermalization model covers both local and extended interactions allowing for interpolation between different systems such as fully ionized gases and neutral atoms, as well as parameters describing integrable and chaotic dynamics.

Discretization and superintegrability all rolled into one (1902.03884v1)

A. V. Tsiganov

2019-02-11

Abelian integrals arise in the mathematical description of various physical processes. According to Abel's theorem these integrals are related to motion of a set of points along a plane curve around fixed points, which are relatively little used in physics applications. We propose to interpret coordinates of the fixed points either as parameters of exact discretization or as additional first integrals for equations of motion reduced to Abelian quadratures.

Recurrence Relations of the Multi-Indexed Orthogonal Polynomials V : Racah and -Racah types (1804.10352v2)

Satoru Odake

2018-04-27

In previous papers, we discussed the recurrence relations of the multi-indexed orthogonal polynomials of the Laguerre, Jacobi, Wilson and Askey-Wilson types. In this paper we explore those of the Racah and -Racah types. For the -indexed (-)Racah polynomials, we derive term recurrence relations with variable dependent coefficients and term () recurrence relations with constant coefficients. Based on the latter, the generalized closure relations and the creation and annihilation operators of the quantum mechanical systems described by the multi-indexed (-)Racah polynomials are obtained. In appendix we present a proof and some data of the recurrence relations with constant coefficients for the multi-indexed Wilson and Askey-Wilson polynomials.

---

Thank you for reading!

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