Geometrical properties of coupled oscillators at synchronization

We study the synchronization of N nearest neighbors coupled oscillators in a ring. We derive an analytic form for the phase difference among neighboring oscillators which shows the dependency on the periodic boundary conditions. At synchronization, we find two distinct quantities which characterize four of the oscillators, two pairs of nearest neighbors, which are at the border of the clusters before total synchronization occurs. These oscillators are responsible for the saddle node bifurcation, of which only two of them have a phase-lock of phase difference equals ± π/2. Using these properties we build a technique based on geometric properties and numerical observations to arrive to an exact analytic expression for the coupling strength at full synchronization and determine the two oscillators that have a phase-lock condition of ± π/2.     

Synchronization dynamics in a ring of four mutually coupled biological systems

This paper considers the synchronization dynamics in a ring of four mutually coupled biological systems described by coupled Van der Pol oscillators. The coupling parameter are non-identical between oscillators. The stability boundaries of the process are first evaluated without the influence of the local injection using the eigenvalues properties and the fourth-order Runge–Kutta algorithm. The effects of a locally injected trajectory on the stability boundaries of the synchronized states are performed using numerical simulations. In both cases, the stability boundaries and the main dynamical states are reported on the stability maps in the (K₁, K₂) plane.     

Fractional dynamics of systems with long-range interaction

We consider one-dimensional chain of coupled linear and nonlinear oscillators with long-range powerwise interaction defined by a term proportional to 1/∣n − m∣^α+1. Continuous medium equation for this system can be obtained in the so-called infrared limit when the wave number tends to zero. We construct a transform operator that maps the system of large number of ordinary differential equations of motion of the particles into a partial differential equation with the Riesz fractional derivative of order α, when 0 < α < 2. Few models of coupled oscillators are considered and their synchronized states and localized structures are discussed in details. Particularly, we discuss some solutions of time-dependent fractional Ginzburg–Landau (or nonlinear Schrodinger) equation.     

New exact solitary-wave special solutions for the nonlinear dispersive K(m, n) equations

In this paper, we study the nonlinear dispersive K(m, n) equations: u_t + (u^m)_x − (uⁿ)_xxx = 0 which exhibit solutions with solitary patterns. New exact solitary solutions are found. The two special cases, K(2, 2) and K(3, 3), are chosen to illustrate the concrete features of the decomposition method in K(m, n) equations. The nonlinear equations K(m, n) are studied for two different cases, namely when m = n being odd and even integers. General formulas for the solutions of K(m, n) equations are established.     

Hyperchaos in fractional order nonlinear systems

We numerically investigate hyperchaotic behavior in an autonomous nonlinear system of fractional order. It is demonstrated that hyperchaotic behavior of the integer order nonlinear system is preserved when the order becomes fractional. The system under study has been reported in the literature [Murali K, Tamasevicius A, Mykolaitis G, Namajunas A, Lindberg E. Hyperchaotic system with unstable oscillators. Nonlinear Phenom Complex Syst 3(1);2000:7–10], and consists of two nonlinearly coupled unstable oscillators, each consisting of an amplifier and an LC resonance loop. The fractional order model of this system is obtained by replacing one or both of its capacitors by fractional order capacitors. Hyperchaos is then assessed by studying the Lyapunov spectrum. The presence of multiple positive Lyapunov exponents in the spectrum is indicative of hyperchaos. Using the appropriate system control parameters, it is demonstrated that hyperchaotic attractors are obtained for a system order less than 4. Consequently, we present a conjecture that fourth-order hyperchaotic nonlinear systems can still produce hyperchaotic behavior with a total system order of 3 + ε, where 1 > ε > 0.     

Global chaos synchronization of new chaotic systems via nonlinear control

Nonlinear control is an effective method for making two identical chaotic systems or two different chaotic systems be synchronized. However, this method assumes that the Lyapunov function of error dynamic (e) of synchronization is always formed as V (e) = 1/2e^Te. In this paper, modification based on Lyapunov stability theory to design a controller is proposed in order to overcome this limitation. The method has been applied successfully to make two identical new systems and two different chaotic systems (new system and Lorenz system) globally asymptotically synchronized. Since the Lyapunov exponents are not required for the calculation, this method is effective and convenient to synchronize two identical systems and two different chaotic systems. Numerical simulations are also given to validate the proposed synchronization approach.     

Delayed transitions in non-linear replicator networks: About ghosts and hypercycles

In this paper we analyze delayed transition phenomena associated to extinction thresholds in a mean field model for hypercycles composed of three and four units, respectively. Hence, we extend a previous analysis carried out with the two-membered hypercycle [see Sardanyés J, Solé RV. Ghosts in the origins of life? Int J Bifurcation Chaos 2006;16(9), in press]. The models we analyze show that, after the tangent bifurcation, these hypercycles also leave a ghost in phase space. These ghosts, which actually conserve the dynamical properties of the coalesced coexistence fixed point, delay the flows before hypercycle extinction. In contrast with the two-component hypercycle, both ghosts show a plateau in the delay as ϕ → 0, thus displacing the power-law dependence to higher values of ϕ, in which the scaling law is now given by τ ∼ ϕ^β, with β = −1/3 (where τ is the delay and ϕ = ϵ − ϵ_c, the parametric distance above the extinction bifurcation point). These results suggest that the presence of the ghost is a general property of hypercycles. Such ghosts actually cause a memory effect which might increase hypercycle survival chances in fluctuating environments.     

Wave propagation properties in oscillatory chains with cubic nonlinearities via nonlinear map approach

Free wave propagation properties in one-dimensional chains of nonlinear oscillators are investigated by means of nonlinear maps. In this realm, the governing difference equations are regarded as symplectic nonlinear transformations relating the amplitudes in adjacent chain sites (n, n + 1) thereby considering a dynamical system where the location index n plays the role of the discrete time. Thus, wave propagation becomes synonymous of stability: finding regions of propagating wave solutions is equivalent to finding regions of linearly stable map solutions. Mechanical models of chains of linearly coupled nonlinear oscillators are investigated. Pass- and stop-band regions of the mono-coupled periodic system are analytically determined for period-q orbits as they are governed by the eigenvalues of the linearized 2D map arising from linear stability analysis of periodic orbits. Then, equivalent chains of nonlinear oscillators in complex domain are tackled. Also in this case, where a 4D real map governs the wave transmission, the nonlinear pass- and stop-bands for periodic orbits are analytically determined by extending the 2D map analysis. The analytical findings concerning the propagation properties are then compared with numerical results obtained through nonlinear map iteration.     

GNGA for general regions: Semilinear elliptic PDE and crossing eigenvalues

We consider the semilinear elliptic PDE Δu + f(λ, u) = 0 with the zero-Dirichlet boundary condition on a family of regions, namely stadions. Linear problems on such regions have been widely studied in the past. We seek to observe the corresponding phenomena in our nonlinear setting. Using the Gradient Newton Galerkin Algorithm (GNGA) of Neuberger and Swift, we document bifurcation, nodal structure, and symmetry of solutions. This paper provides the first published instance where the GNGA is applied to general regions. Our investigation involves both the dimension of the stadions and the value λ as parameters. We find that the so-called crossings and avoided crossings of eigenvalues as the dimension of the stadions vary influences the symmetry and variational structure of nonlinear solutions in a natural way.     

Exact special solitary solutions with compact support for the nonlinear dispersive K(m, n) equations

The nonlinear dispersive K(m, n) equations, u_t−(u^m)_x−(uⁿ)_xxx = 0 which exhibit compactons: solitons with compact support, are studied. New exact solitary solutions with compact support are found. The two special cases, K(2, 2) and K(3, 3), are chosen to illustrate the concrete features of the decomposition method in K(m, n) equations. General formulas for the solutions of K(m, n) equations are established.     

Synchronization of two low-dimensional Kerr chains

Synchronization of two chaotic low-dimensional chains (α₁, α₂, α₃) and (A₁, A₂, A₃) consisting of Kerr oscillators is studied. The synchronization has been achieved by the parallel coupling of α₁ with A₁, α₂ with A₂ and α₃ with A₃. We want to find whether and when the pairs (α₁, A₁), (α₂, A₂) and (α₃, A₃) synchronize non-simultaneously (three-time synchronism). The problem of synchronization is also studied for a number of couplings between the chains lower than the number of oscillators in a single chain. Both the ring and linear geometry of synchronization is investigated. The presented results suggest a possibility of multi-time synchronism in two coupled high-dimensional chains. It seems very promising for design of some devices for advanced signal processing.     

An algebraic equation for linear operators

Let T be a linear operator on a vector space V, possibly of infinite dimension, over a general field K. We solve the functional equation p(T) = F where p ∈ K[x] and F, an algebraic operator on V, are given. For nilpotent F we give an explicit linear system which determines the solutions by their similarity classes. The method is based on a canonical decomposition theorem.     

Pattern Frequency Sequences and Internal Zeros

Let q be a pattern and let S_n, q(c) be the number of n-permutations having exactly c copies of q. We investigate when the sequence (S_n, q(c))_c ≥ 0 has internal zeros. If q is a monotone pattern it turns out that, except for q = 12 or 21, the nontrivial sequences (those where n is at least the length of q) always have internal zeros. For the pattern q = 1(l + 1)l…2 there are infinitely many sequences which contain internal zeros and when l = 2 there are also infinitely many which do not. In the latter case, the only possible places for internal zeros are the next-to-last or the second-to-last positions. Note that by symmetry this completely determines the existence of internal zeros for all patterns of length at most 3.     

On the stagnation-point flow towards a stretching sheet with homogeneous–heterogeneous reactions effects

The effects of homogeneous–heterogeneous reactions on the steady boundary layer flow near the stagnation point on a stretching surface is studied. The possible steady-states of this system are analyzed in the case when the diffusion coefficients of both reactant and auto catalyst are equal. The strength of this effect is represented by the dimensionless parameter K and K_s. It is shown that for a fluid of small kinematic viscosity, a boundary layer is formed when the stretching velocity is less than the free stream velocity and an inverted boundary layer is formed when the stretching velocity exceeds the free stream velocity. The uniqueness of this problem lies on the fact that the solutions are possible for all values of λ > 0 (stretching surface), while for λ < 0 (shrinking surface), solutions are possible only for its limited range.     

On the positive definite solutions of nonlinear matrix equation X + A⋆ X−δA = Q

In this paper we consider the positive definite solutions of nonlinear matrix equation X + A^⋆X^−δA = Q, where δ ∈ (0, 1], which appears for the first time in [S.M. El-Sayed, A.C.M. Ran, On an iteration methods for solving a class of nonlinear matrix equations, SIAM J. Matrix Anal. Appl. 23 (2001) 632–645]. The necessary and sufficient conditions for the existence of a solution are derived. An iterative algorithm for obtaining the positive definite solutions of the equation is discussed. The error estimations are found.     

Bifurcations of travelling wave solutions in the (N + 1)-dimensional sine–cosine-Gordon equations

In this paper, the (N + 1)-dimensional sine–cosine-Gordon equations are studied. The existence of solitary wave, kink and anti-kink wave, and periodic wave solutions are proved, by using the method of bifurcation theory of dynamical systems. All possible bounded exact explicit parametric representations of the above travelling solutions are obtained.     

The effect of a small background inhomogeneity on the asymptotic properties of linear perturbations

General regularities in the evolution of one-dimensional unstable linear perturbations on a weakly inhomogeneous background are studied when, at the initial instant, the perturbations are concentrated in the δ-neighbourhood of a certain point. Times are considered when these perturbations do not fall outside the limits of a certain domain of size l such that δ ? l ? L, where L is the large characteristic size of the background inhomogeneity. With contain assumptions, the effect of the background inhomogeneity on the asymptotic behaviour of the perturbations at long times is taken into account in a general form. The first corrections to the well known asymptotic relation for the evolution of perturbations on a homogeneous background, that arise because of background inhomogeneity, are obtained using Hamilton's method. An example of the use of the proposed approximate method is considered and the error in the approximation is estimated.     

Complex dynamics in Duffing–Van der Pol equation

Duffing–Van der Pol equation with fifth nonlinear-restoring force and two external forcing terms is investigated. The threshold values of existence of chaotic motion are obtained under the periodic perturbation. By second-order averaging method and Melnikov method, we prove the criterion of existence of chaos in averaged system under quasi-periodic perturbation for ω₂ = nω₁ + εσ, n = 1, 3, 5, and cannot prove the criterion of existence of chaos in second-order averaged system under quasi-periodic perturbation for ω₂ = nω₁ + εσ, n = 2, 4, 6, 7, 8, 9, 10, where σ is not rational to ω₁, but can show the occurrence of chaos in original system by numerical simulation. Numerical simulations including heteroclinic and homoclinic bifurcation surfaces, bifurcation diagrams, Lyapunov exponent, phase portraits and Poincaré map, not only show the consistence with the theoretical analysis but also exhibit the more new complex dynamical behaviors. We show that cascades of interlocking period-doubling and reverse period-doubling bifurcations from period-2 to -4 and -6 orbits, interleaving occurrence of chaotic behaviors and quasi-periodic orbits, transient chaos with a great abundance of period windows, symmetry-breaking of periodic orbits in chaotic regions, onset of chaos which occurs more than one, chaos suddenly disappearing to period orbits, interior crisis, strange non-chaotic attractor, non-attracting chaotic set and nice chaotic attractors. Our results show many dynamical behaviors and some of them are strictly departure from the behaviors of Duffing–Van der Pol equation with a cubic nonlinear-restoring force and one external forcing.     

The stochastic soliton-like solutions of stochastic KdV equations

By means of a generalized method and symbolic computation, we consider a stochastic KdV equation U_t + f(t)U ♢ U_x + g(t)U_xxx = W(t) ♢ R^♢(t, U, U_x, U_xxx). We construct new and more general formal solutions. At the same time, we recover all the solutions found by Xie [Phys. Lett. A 310 (2003) 161]. The solutions obtained include the nontravelling wave and coefficient function’s stochastic soliton-like solutions, singular stochastic soliton-like solutions, stochastic triangular functions solutions.     

Comment on the 3+1 dimensional Kadomtsev–Petviashvili equations

We comment on traveling wave solutions and rational solutions to the 3+1 dimensional Kadomtsev–Petviashvili (KP) equations: (u_t + 6uu_x + u_xxx)_x ± 3u_yy ± 3u_zz = 0. We also show that both of the 3+1 dimensional KP equations do not possess the three-soliton solution. This suggests that none of the 3+1 dimensional KP equations should be integrable, and partially explains why they do not pass the Painlevé test. As by-products, the one-soliton and two-soliton solutions and four classes of specific three-soliton solutions are explicitly presented.