Nonlinear Sciences

We report the diversity of scroll wave chimeras in the three-dimensional (3D) Kuramoto model with inertia for N^{3} identical phase oscillators placed in a unit 3D cube with periodic boundary conditions. In the considered model with inertia, we have found patterns which do not exist in this system without inertia. In particular, a scroll ring chimera is obtained from random initial conditions. In contrast to this system without inertia, where all chimera states have incoherent inner parts, these states can have partially coherent or fully coherent inner parts as exemplified by a scroll ring chimera. Solitary states exist in the considered model as separate states or can coexist with scroll wave chimeras in the oscillatory space. We also propose a method of construction of 3D images using solitary states as solutions of the 3D Kuramoto model with inertia.

The coexistence of coherent and incoherent domains, namely the appearance of chimera states, is being studied extensively in many contexts of science and technology since the past decade, though the previous studies are mostly built on the framework of one-dimensional and two-dimensional interaction topologies. Recently, the emergence of such fascinating phenomena has been studied in three-dimensional (3D) grid formation while considering only the nonlocal interaction. Here we study the emergence and existence of chimera patterns in a three dimensional network of coupled Stuart-Landau limit cycle oscillators and Hindmarsh-Rose neuronal oscillators with local (nearest neighbour) interaction topology. The emergence of different types of spatiotemporal chimera patterns is investigated by taking two distinct nonlinear interaction functions. We provide appropriate analytical explanations in the 3D grid of network formation and the corresponding numerical justifications are given. We extend our analysis on the basis of Ott-Antonsen reduction approach in the case of Stuart-Landau oscillators containing infinite number of oscillators. Particularly, in Hindmarsh-Rose neuronal network the existence of non-stationary chimera states are characterized by instantaneous strength of incoherence and instantaneous local order parameter. Besides, the condition for achieving exact neuronal synchrony is obtained analytically through a linear stability analysis. The different types of collective dynamics together with chimera states are mapped over a wide range of various parameter spaces.

We study the existence of chimera states, i.e. mixed states, in a globally coupled sine circle map lattice, with different strengths of inter-group and intra-group coupling. We find that at specific values of the parameters of the CML, a completely random initial condition evolves to chimera states, having a phase synchronised and a phase desynchronised group, where the space time variation of the phases of the maps in the desynchronised group shows structures similar to spatiotemporally intermittent regions. Using the complex order parameter we obtain a phase diagram that identifies the region in the parameter space which supports chimera states of this type, as well as other types of phase configurations such as globally phase synchronised states, two phase clustered states and fully phase desynchronised states. We estimate the volume of the basin of attraction of each kind of solution. The STI chimera region is studied in further detail via numerical and analytic stability analysis, and the Lyapunov spectrum is calculated. This state is identified to be hyperchaotic as the two largest Lyapunov exponents are found to be positive. The distributions of laminar and burst lengths in the incoherent region of the chimera show exponential behaviour. The average fraction of laminar/burst sites is identified to be the important quantity which governs the dynamics of the chimera. After an initial transient, these settle to steady values which can be used to reproduce the phase diagram in the chimera regime.

Chimera states, which consist of coexisting synchronous and asynchronous domains in networks of coupled oscillators, are in the focus of attention for over a decade. Although chimera morphology and properties have been investigated in a number of models, the mechanism responsible for their formation is still not well understood. To shed light in the chimera producing mechanism, in the present study we introduce an oscillatory model with variable frequency governed by a 3rd order equation. In this model single oscillators are constructed as bistable and depending on the initial conditions their frequency may result in one of the two stable fixed points, ω l and ω h (two-level synchronization). Numerical simulations demonstrate that these oscillators organize in domains with alternating frequencies, when they are nonlocally coupled in networks. In each domain the oscillators synchronize, sequential domains follow different modes of synchronization and the border elements between two consecutive domains form the asynchronous domains. We investigate the influence of the frequency coupling constant and of the coupling range on the chimera morphology and we show that the chimera multiplicity decreases as the coupling range increases. The frequency spectrum is calculated in the coherent and incoherent domains of this model. In the coherent domains single frequencies ( ω l or ω h ) are observed, while in the incoherent domains both ω l and ω h as well as their superpositions appear. This mechanism of creating domains of alternating frequencies offers a reasonable generic scenario for chimera state formation.

We study classical and wave chaotic dynamics of light both on a curved surface and on a table billiard with nonuniform distribution of refractive index. Inspired by the concept of transformation optics, we demonstrate that these two systems are fundamentally equivalent in terms of both light rays and waves, under a conformal coordinate transformation. This total equivalency allows us to study simultaneously a typical family of curved surfaces and their corresponding transformed billiards. Here we consider the truncated Tannery's pear and its equivalent flat billiard with non-uniform distribution of refractive index, when an off-centered hole is pierced to introduce chaos. We find that the degree of chaos is fully controlled by the single geometric parameter of Tannery's pears. This is proved by exploring in the transformed billiard the dependence with this geometric parameter of the Poincaré surface of section, the Lyapunov exponent and the the statistics of eigenmodes and eigenfrequency spectrum. Finally, a simple interpretation of our findings naturally emerges when considering transformed billiards, which allows to extend our prediction to other class of curved surfaces. The powerful analogy we reveal here between two a priori unrelated systems not only brings forward a novel approach to control the degree of chaos, but also provides potentialities for further studies and applications in various fields, such as billiards design, optical fibers, or laser microcavities.

Synchronization among rhythmic elements is modeled by coupled phase-oscillators each of which has the so-called natural frequency. A symmetric natural frequency distribution induces a continuous or discontinuous synchronization transition from the nonsynchronized state, for instance. It has been numerically reported that asymmetry in the natural frequency distribution brings new types of bifurcation diagram having, in the order parameter, oscillation or a discontinuous jump which emerges from a partially synchronized state. We propose a theoretical classification method of five types of bifurcation diagrams including the new ones, paying attention to the generality of the theory. The oscillation and the jump from partially synchronized states are discussed respectively by the linear analysis around the nonsynchronized state and by extending the amplitude equation up to the third leading term. The theoretical classification is examined by comparing with numerically obtained one.

By a small-size complex network of coupled chaotic Hindmarsh-Rose circuits, we study experimentally the stability of network synchronization to the removal of shortcut links. It is shown that the removal of a single shortcut link may destroy either completely or partially the network synchronization. Interestingly, when the network is partially desynchronized, it is found that the oscillators can be organized into different groups, with oscillators within each group being highly synchronized but are not for oscillators from different groups, showing the intriguing phenomenon of cluster synchronization. The experimental results are analyzed by the method of eigenvalue analysis, which implies that the formation of cluster synchronization is crucially dependent on the network symmetries. Our study demonstrates the observability of cluster synchronization in realistic systems, and indicates the feasibility of controlling network synchronization by adjusting network topology.

We explore the coherent dynamics in a small network of three coupled parametric oscillators and demonstrate the effect of frustration on the persistent beating between them. Since a single-mode parametric oscillator represents an analog of a classical Ising spin, networks of coupled parametric oscillators are considered as simulators of Ising spin models, aiming to efficiently calculate the ground state of an Ising network - a computationally hard problem. However, the coherent dynamics of coupled parametric oscillators can be considerably richer than that of Ising spins, depending on the nature of the coupling between them (energy preserving or dissipative), as was recently shown for two coupled parametric oscillators. In particular, when the energy-preserving coupling is dominant, the system displays everlasting coherent beats, transcending the Ising description. Here, we extend these findings to three coupled parametric oscillators, focusing in particular on the effect of frustration of the dissipative coupling. We theoretically analyze the dynamics using coupled nonlinear Mathieu's equations, and corroborate our theoretical findings by a numerical simulation that closely mimics the dynamics of the system in an actual experiment. Our main finding is that frustration drastically modifies the dynamics. While in the absence of frustration the system is analogous to the two-oscillator case, frustration reverses the role of the coupling completely, and beats are found for small energy-preserving couplings.

In this study, we propose a universal scenario explaining the 1/f fluctuation, including pink noises, in Hamiltonian dynamical systems with many degrees of freedom under long-range interaction. In the thermodynamic limit, the dynamics of such systems can be described by the Vlasov equation, which has an infinite number of Casimir invariants. In a finite system, they become pseudoinvariants, which yield quasistationary states. The dynamics then exhibit slow motion over them, up to the timescale where the pseudo-Casimir-invariants are effective. Such long-time correlation leads to 1/f fluctuations of collective variables, as is confirmed by direct numerical simulations. The universality of this collective 1/f fluctuation is demonstrated by taking a variety of Hamiltonians and changing the range of interaction and number of particles.

We show how to couple phase-oscillators on a graph so that collective dynamics "searches" for the coloring of that graph as it relaxes toward the dynamical equilibrium. This translates a combinatorial optimization problem (graph coloring) into a functional optimization problem (finding and evaluating the global minimum of dynamical non-equilibrium potential, done by the natural system's evolution). Using a sample of graphs, we show that our method can serve as a viable alternative to the traditional combinatorial algorithms. Moreover, we show that, with the same computational cost, our method efficiently solves the harder problem of improper coloring of weighed graphs.