Chaotic behavior in Lotka-Volterra and May-Leonard models of biodiversity.

Quantification of chaos is a challenging issue in complex dynamical systems. In this paper, we discuss the chaotic properties of generalized Lotka-Volterra and May-Leonard models of biodiversity, via the Hamming distance density. We identified chaotic behavior for different scenarios via the specific features of the Hamming distance and the method of q-exponential fitting. We also investigated the spatial autocorrelation length to find the corresponding characteristic length in terms of the number of species in each system. In particular, the results concerning the characteristic length are in good accordance with the study of the chaotic behavior implemented in this work.

Environment driven oscillation in an off-lattice May-Leonard model.

Cyclic dominance of competing species is an intensively used working hypothesis to explain biodiversity in certain living systems, where the evolutionary selection principle would dictate a single victor otherwise. Technically the May-Leonard models offer a mathematical framework to describe the mentioned non-transitive interaction of competing species when individual movement is also considered in a spatial system. Emerging rotating spirals composed by the competing species are frequently observed character of the resulting patterns. But how do these spiraling patterns change when we vary the external environment which affects the general vitality of individuals? Motivated by this question we suggest an off-lattice version of the tradition May-Leonard model which allows us to change the actual state of the environment gradually. This can be done by introducing a local carrying capacity parameter which value can be varied gently in an off-lattice environment. Our results support a previous analysis obtained in a more intricate metapopulation model and we show that the well-known rotating spirals become evident in a benign environment when the general density of the population is high. The accompanying time-dependent oscillation of competing species can also be detected where the amplitude and the frequency show a scaling law of the parameter that characterizes the state of the environment. These observations highlight that the assumed non-transitive interaction alone is insufficient condition to maintain biodiversity safely, but the actual state of the environment, which characterizes the general living conditions, also plays a decisive role on the evolution of related systems.

Invasion-controlled pattern formation in a generalized multispecies predator-prey system.

Rock-scissors-paper game, as the simplest model of intransitive relation between competing agents, is a frequently quoted model to explain the stable diversity of competitors in the race of surviving. When increasing the number of competitors we may face a novel situation because beside the mentioned unidirectional predator-prey-like dominance a balanced or peer relation can emerge between some competitors. By utilizing this possibility in the present work we generalize a four-state predator-prey-type model where we establish two groups of species labeled by even and odd numbers. In particular, we introduce different invasion probabilities between and within these groups, which results in a tunable intensity of bidirectional invasion among peer species. Our study reveals an exceptional richness of pattern formations where five quantitatively different phases are observed by varying solely the strength of the mentioned inner invasion. The related transition points can be identified with the help of appropriate order parameters based on the spatial autocorrelation decay, on the fraction of empty sites, and on the variance of the species density. Furthermore, the application of diverse, alliance-specific inner invasion rates for different groups may result in the extinction of the pair of species where this inner invasion is moderate. These observations highlight that beyond the well-known and intensively studied cyclic dominance there is an additional source of complexity of pattern formation that has not been explored earlier.

How directional mobility affects coexistence in rock-paper-scissors models.

This work deals with a system of three distinct species that changes in time under the presence of mobility, selection, and reproduction, as in the popular rock-paper-scissors game. The novelty of the current study is the modification of the mobility rule to the case of directional mobility, in which the species move following the direction associated to a larger (averaged) number density of selection targets in the surrounding neighborhood. Directional mobility can be used to simulate eyes that see or a nose that smells, and we show how it may contribute to reduce the probability of coexistence.

Apex predator and the cyclic competition in a rock-paper-scissors game of three species.

This work deals with the effects of an apex predator on the cyclic competition among three distinct species that follow the rules of the rock-paper-scissors game. The investigation develops standard stochastic simulations but is motivated by a procedure which is explained in the work. We add the apex predator as the fourth species in a system that contains three species that evolve following the standard rules of migration, reproduction, and predation, and study how the system evolves in this new environment, in comparison with the case in the absence of the apex predator. The results show that the apex predator engenders the tendency to spread uniformly in the lattice, contributing to destroy the spiral patterns, keeping biodiversity but diminishing the average size of the clusters of the species that compete cyclically.

A novel procedure for the identification of chaos in complex biological systems.

We demonstrate the presence of chaos in stochastic simulations that are widely used to study biodiversity in nature. The investigation deals with a set of three distinct species that evolve according to the standard rules of mobility, reproduction and predation, with predation following the cyclic rules of the popular rock, paper and scissors game. The study uncovers the possibility to distinguish between time evolutions that start from slightly different initial states, guided by the Hamming distance which heuristically unveils the chaotic behavior. The finding opens up a quantitative approach that relates the correlation length to the average density of maxima of a typical species, and an ensemble of stochastic simulations is implemented to support the procedure. The main result of the work shows how a single and simple experimental realization that counts the density of maxima associated with the chaotic evolution of the species serves to infer its correlation length. We use the result to investigate others distinct complex systems, one dealing with a set of differential equations that can be used to model a diversity of natural and artificial chaotic systems, and another one, focusing on the ocean water level.

Universal Braess paradox in open quantum dots.

We present analytical and numerical results that demonstrate the presence of the Braess paradox in chaotic quantum dots. The paradox that we identify, originally perceived in classical networks, shows that the addition of more capacity to the network can suppress the current flow in the universal regime. We investigate the weak localization term, showing that it presents the paradox encoded in a saturation minimum of the conductance, under the presence of hyperflow in the external leads. In addition, we demonstrate that the weak localization suffers a transition signal depending on the overcapacity lead and presents an echo on the magnetic crossover before going to zero due to the full time-reversal symmetry breaking. We also show that the quantum interference contribution can dominate the Ohm term in the presence of constrictions and that the corresponding Fano factor engenders an anomalous behavior.

Interfaces with internal structures in generalized rock-paper-scissors models.

In this work we investigate the development of stable dynamical structures along interfaces separating domains belonging to enemy partnerships in the context of cyclic predator-prey models with an even number of species N≥8. We use both stochastic and field theory simulations in one and two spatial dimensions, as well as analytical arguments, to describe the association at the interfaces of mutually neutral individuals belonging to enemy partnerships and to probe their role in the development of the dynamical structures at the interfaces. We identify an interesting behavior associated with the symmetric or asymmetric evolution of the interface profiles depending on whether N/2 is odd or even, respectively. We also show that the macroscopic evolution of the interface network is not very sensitive to the internal structure of the interfaces. Although this work focuses on cyclic predator-prey models with an even number of species, we argue that the results are expected to be quite generic in the context of spatial stochastic May-Leonard models.

Bright solitons from the nonpolynomial Schrödinger equation with inhomogeneous defocusing nonlinearities.

Extending the recent work on models with spatially nonuniform nonlinearities, we study bright solitons generated by the nonpolynomial self-defocusing (SDF) nonlinearity in the framework of the one-dimensional (1D) Muñoz-Mateo-Delgado (MM-D) equation (the 1D reduction of the Gross-Pitaevskii equation with the SDF nonlinearity), with the local strength of the nonlinearity growing at |x|→∞ faster than |x|. We produce numerical solutions and analytical ones, obtained by means of the Thomas-Fermi approximation, for nodeless ground states and for excited modes with one, two, three and four nodes, in two versions of the model, with steep (exponential) and mild (algebraic) nonlinear-modulation profiles. In both cases, the ground states and the single-node ones are completely stable, while the stability of the higher-order modes depends on their norm (in the case of the algebraic modulation, they are fully unstable). Unstable states spontaneously evolve into their stable lower-order counterparts.

von Neummann's and related scaling laws in rock-paper-scissors-type games.

We introduce a family of rock-paper-scissors-type models with Z(N) symmetry (N is the number of species), and we show that it has a very rich structure with many completely different phases. We study realizations that lead to the formation of domains, where individuals of one or more species coexist, separated by interfaces whose (average) dynamics is curvature driven. This type of behavior, which might be relevant for the development of biological complexity, leads to an interface network evolution and pattern formation similar to the ones of several other nonlinear systems in condensed matter and cosmology.

Junctions and spiral patterns in generalized rock-paper-scissors models.

We investigate the population dynamics in generalized rock-paper-scissors models with an arbitrary number of species N. We show that spiral patterns with N arms may develop both for odd and even N, in particular in models where a bidirectional predation interaction of equal strength between all species is modified to include one N-cyclic predator-prey rule. While the former case gives rise to an interface network with Y-type junctions obeying the scaling law L∝t1/2, where L is the characteristic length of the network and t is the time, the latter can lead to a population network with N-armed spiral patterns, having a roughly constant characteristic length scale. We explicitly demonstrate the connection between interface junctions and spiral patterns in these models and compute the corresponding scaling laws. This work significantly extends the results of previous studies of population dynamics and could have profound implications for the understanding of biological complexity in systems with a large number of species.

Modulation of localized solutions in a system of two coupled nonlinear Schrödinger equations.

In this work we study localized solutions of a system of two coupled nonlinear Schrödinger equations, with the linear (potential) and nonlinear coefficients engendering spatial and temporal dependencies. Similarity transformations are used to convert the nonautonomous coupled equations into autonomous ones and we use the trial orbit method to help us solving them, presenting solutions in a general way. Numerical experiments are then used to verify the stability of the localized solutions.

Conductance peaks in open quantum dots.

We present a simple measure of the conductance fluctuations in open ballistic chaotic quantum dots, extending the number of maxima method originally proposed for the statistical analysis of compound nuclear reactions. The average number of extreme points (maxima and minima) in the dimensionless conductance T as a function of an arbitrary external parameter Z is directly related to the autocorrelation function of T(Z). The parameter Z can be associated with an applied gate voltage causing shape deformation in quantum dot, an external magnetic field, the Fermi energy, etc. The average density of maxima is found to be <ρ(Z)>=α(Z)/Z(c), where α(Z) is a universal constant and Z(c) is the conductance autocorrelation length, which is system specific. The analysis of <ρ(Z)> does not require large statistic samples, providing a quite amenable way to access information about parametric correlations, such as Z(c).

One-dimensional reduction of the three-dimenstional Gross-Pitaevskii equation with two- and three-body interactions.

We deal with the three-dimensional Gross-Pitaevskii equation which is used to describe a cloud of dilute bosonic atoms that interact under competing two- and three-body scattering potentials. We study the case where the cloud of atoms is strongly confined in two spatial dimensions, allowing us to build an unidimensional nonlinear equation,controlled by the nonlinearities and the confining potentials that trap the system along the longitudinal coordinate. We focus attention on specific limits dictated by the cubic and quintic coefficients, and we implement numerical simulations to help us to quantify the validity of the procedure.

Modulation of breathers in the three-dimensional nonlinear Gross-Pitaevskii equation.

In this paper we present analytical breather solutions of the three-dimensional nonlinear generalized Gross-Pitaevskii equation. We use an Ansatz to reduce the three-dimensional equation with space- and time-dependent coefficients into a one-dimensional equation with constant coefficients. The key point is to show that both the space- and time-dependent coefficients of the nonlinear equation can contribute to modulate the breather excitations. We briefly discuss the experimental feasibility of the results in Bose-Einstein condensates.

Solitons with cubic and quintic nonlinearities modulated in space and time.

This work deals with soliton solutions of the nonlinear Schrödinger equation with cubic and quintic nonlinearities. We extend the procedure put forward in a recent paper [J. Belmonte-Beitia, Phys. Rev. Lett. 100, 164102 (2008)], and we solve the equation in the presence of a linear background and cubic and quintic interactions which are modulated in space and time. As a result, we show how a simple parameter can be used to generate brightlike or darklike localized nonlinear waves which oscillate in several distinct ways, driven by the space and time dependence of the parameters that control the trapping potential and the cubic and quintic nonlinearities.

New global defect structures.

We investigate the presence of defects in systems described by real scalar field in (D,1) spacetime dimensions. We show that when the potential assumes specific form, there are models which support stable global defects for D arbitrary. We also show how to find first-order differential equations that solve the equations of motion, and how to solve models in D dimensions via soluble problems in D=1. We illustrate the procedure examining specific models and finding explicit solutions.

Tiling the plane without supersymmetry

We present a way of tiling the plane with a regular hexagonal network of defects. The network is stable and follows in consequence of the three-junctions that appear in a model of two real scalar fields that presents Z3 symmetry. The Z3 symmetry is effective in both the vacuum and defect sectors, and no supersymmetry is required to build the network.