International Journal of Bifurcation and Chaos

Robust synchrony in network dynamics is governed by balanced colorings and the corresponding quotient network, also formalized in terms of graph fibrations. Dynamics and bifurcations are constrained — often in surprising ways — by the associated synchrony subspaces, which are invariant under all admissible ordinary differential equations (ODEs). The class of admissible ODEs is determined by a groupoid, whose objects are the input sets of nodes and whose morphisms are input isomorphisms between those sets. We define the coloring subgroupoid corresponding to a coloring, leading to groupoid interpretations of colorings and quotient networks. The first half of the paper is mainly tutorial. The second half, which is new, characterizes the structure of the network groupoid and proves that the groupoid of the quotient network is the quotient of the network groupoid by a normal subgroupoid of transition elements.

In this paper, a mathematical model of fluid flows in a convective thermal system is developed, and a five-dimensional dynamical system is developed for the investigation of the convective fluid dynamics. The analytical solutions of periodic motions to chaos of the convective fluid flows are developed for steady-state vortex flows, and the corresponding stability and bifurcations of periodic motions in the five-dimensional dynamical system are studied. The harmonic frequency-amplitude characteristics for periodic flows are obtained, which provide energy distribution in the parameter space. Analytical homoclinic orbits for the convective fluid flow systems are developed for the asymptotic convection through the infinite-many homoclinic orbits in the five-dimensional dynamical system. The dynamics of fluid flows in the convective thermal systems are revealed, and one can use such methodology to predict atmospheric and oceanic phenomena through thermal convections.

In this paper, we present a concise four-dimensional (4D) conservative chaotic system with a wide parameter range. Since there are no terms higher than first order, the circuit does not contain multipliers, resulting in a simple circuit implementation. The nonlinear dynamic characteristics, such as phase diagrams, equilibrium points, divergence, Poincaré cross-sections, Lyapunov exponents, bifurcation diagrams, and Lyapunov dimension, are analyzed in detail, which illustrates the conservativity. Besides, the system exhibits different offset boosting behaviors. Through offset boosting, the system can propagate along a line, convert signal polarity, control variable amplitude, generate coexisting attractors, and even induce changes in its state. Specially, we realize the transition from a single-vortex attractor to a multivortex one by some changes in the initial values. Furthermore, the Pearson correlation coefficient is used to demonstrate the higher initial value sensitivity of the system. Finally, the system is implemented through Multisim simulation and analog circuit separately, and their consistency validates the system effectively.

For the quadratic reversible systems of genus one, all of their periodic orbits are higher-order algebraic curves. When they are perturbed by polynomials of degree $n$, the numbers of zeros of their Abelian integrals will change and we study the upper bounds of these numbers by using the methods of Riccati equation and Picard–Fuchs equation. We consider both the highest and lowest degrees of polynomials, and more importantly, we consider the law of polynomials and the range of values for their variables. Consequently, some laws of the polynomials are discovered and many upper bounds are obtained, and these upper bounds are sharper than the results obtained by other techniques.

The rarity of species increases its market price, consequently leading to the overexploitation of the species and even the extinction of the species. We study how the harvest intensity and the additive Allee effect impact on the Gordon–Schaefer model. In addition, by Sotomayor’s theorem and Poincaré–Andronov theorem, we prove the existence of Hopf bifurcation, saddle-node bifurcation and transcritical bifurcation, respectively. Finally, we illustrate our results by numerical simulations. We find that both the cost per unit of harvest and the additive Allee effect have a significant impact on human exploitation of the population. As the additive Allee effect reduces to the weak Allee effect, the lower harvest cost encourages humans to increase the exploitation of species. This threshold is a switch that controls the strong Allee effect. If it exceeds its threshold, then the motivation of humans to exploit the species increases.

In this paper, the qualitative properties of a physically extended six-dimensional Lorenz system, with additional physical terms describing rotation and density, which was proposed in [Moon et al., 2019] have been investigated. The dissipation, invariance, Lyapunov exponents, Kaplan–Yorke dimension, ultimate boundedness and global attractivity of this six-dimensional Lorenz system have been discussed in detail according to the chaotic systems theory. We find that this system exhibits chaos phenomena for a new set of parameters. It is well known that the general method for studying the bounds of a chaotic system is to construct a suitable Lyapunov-like function (or the generalized positive definite and radically unbounded Lyapunov function). However, the higher the dimension of a chaotic system, the more difficult it is to construct the Lyapunov-like function. The innovation of this paper is that we first construct the suitable Lyapunov-like function for this six-dimensional Lorenz system, and then we prove that this system is not only globally bounded for varying parameters, but it also gives a collection of global absorbing sets for this system with respect to all parameters of this system according to Lyapunov’s direct method and the optimization method. Furthermore, we obtain the rate of the trajectories going from the exterior to the global absorbing set. Some numerical simulations are presented to validate our research results. Finally, we give a direct application of the results obtained in this paper. According to the results of this paper, we can conclude that the equilibrium point $O(0,0,0,0,0,0)$ of this system is globally exponentially stable.

We study a family of one-dimensional quasi-periodically forced maps $F_{a,b}(x,𝜃)=(f_{a,b}(x,𝜃),𝜃+\omega )$, where $x$ is real, $𝜃$ is an angle, and $\omega$ is an irrational frequency, such that $f_{a,b}(x,𝜃)$ is a real piecewise-linear map with respect to $x$ of certain kind. The family depends on two real parameters, $a>0$ and $b>0$. For this family, we prove the existence of nonsmooth pitchfork bifurcations. For $a<1$ and any $b,$ there is only one continuous invariant curve. For $a>1,$ there exists a smooth map $b=b_0(a)$ such that: (a) For $b<b_0(a)$, $f_{a,b}$ has two continuous attracting invariant curves and one continuous repelling curve; (b) For $b=b_0(a),$ it has one continuous repelling invariant curve and two semi-continuous (noncontinuous) attracting invariant curves that intersect the unstable one in a zero-Lebesgue measure set of angles; (c) For $b>b_0(a),$ it has one continuous attracting invariant curve. The case $a=1$ is a degenerate case that is also discussed in the paper. It is interesting to note that this family is a simplified version of the smooth family $G_{a,b}(x,𝜃)=(arctan(ax)+bsin(𝜃),𝜃+\omega )$ for which there is numerical evidence of a nonsmooth pitchfork bifurcation. Finally, we also discuss the limit case when $a\to \infty$.

The impact of power-law memory on the dynamics of the Caputo standard fractional map is investigated in this paper. The definition of a complexity measure for the Caputo standard fractional map is introduced. This measure evaluates both the average algebraic complexity of a trajectory and the distribution of the trajectories of different types in the phase space of the system. The interplay between the small-scale spatial chaos and the large-scale spatial behavior is observed and measured during the transition of the Caputo standard fractional map from the circle map to the classical standard map. It is demonstrated that the impact of the fractional derivative on the complexity of the fractional system is not straightforward and is predetermined by the physical properties governing the dynamics of that system.

The slow–fast Leslie–Gower model with piecewise-smooth Holling type-I functional response and weak Allee effect is studied in this paper. It is shown that the model undergoes singular Hopf bifurcation and nonsmooth Hopf bifurcation as the parameters vary. The theoretical analysis implies that the predator’s food quality and Allee effect play an important role and lead to richer dynamical phenomena such as the globally stable equilibria, canard explosion phenomenon, a hyperbolically stable relaxation oscillation cycle enclosing almost two canard cycles with different stabilities and so on. Moreover, the predator and prey will coexist as multiple steady states or periodic oscillations for different positive initial populations and positive parameter values. Finally, we present some numerical simulations to illustrate the theoretical analysis such as the existence of one, two or three limit cycles.

For their good randomness and long iteration periods, chaotic maps have been widely used in cryptography. Recently, we have revealed the correlation between Lyapunov exponent and sequence randomness of multidimensional chaotic maps based on modular operation. Since the modular operation can realize the boundedness of chaotic state points, it is important to further reveal the deterministic correlation between Lyapunov exponent and modulus. First, we constructed an $n$-dimensional nondegenerate hyperchaotic map model with the desired Lyapunov exponents. Then, we gave the existence and uniqueness proof of quadrature rectangle decomposition theorem and revealed the correlation between Lyapunov exponent and modulus. The novelty lies in that (1) in order to realize the irreversibility of the iterative processes of chaotic maps, we constructed a chaotic map based on modular exponentiation, and its inverse function is the discrete logarithm problem; and (2) we reveal for the first time the correlation between Lyapunov exponent and modulus, and give the lower bound of the modulus of the nondegenerate chaotic map. In addition, to verify the effectiveness of the scheme, we constructed four-dimensional and five-dimensional chaotic maps, respectively, and analyzed their dynamical behaviors, and the results revealed that there exist linear or nonlinear correlation between Lyapunov exponent and modulus.

In this paper, we study the dynamic behaviors and chaos synchronization of the Chen system described by Caputo–Hadamard fractional derivative. First, the existence and uniqueness of a solution to the Chen system with Caputo–Hadamard derivative are proved by qualitative analysis. Further, the stability of equilibria of the considered system is analyzed with the aid of Routh–Hurwitz criteria. Meanwhile, the bifurcation condition of the Caputo–Hadamard Chen system is compared with the integer-order Chen system, where the differences between the two systems are demonstrated numerically. In the study of chaos synchronization of the drive–response Chen systems with Caputo–Hadamard derivative, two control schemes are developed: three nonlinear controllers and single linear controller. The feasibility of two control schemes is verified, and the synchronization performances of these two schemes are compared by numerical simulations. Based on this, the influence of the fractional-order on chaos synchronization performance is illustrated as well.

In this work, we present a method for the Global Asymptotic Stabilization (GAS) of an affine control chaotic Lorenz system, via admissible (bounded and regular) feedback controls, where the control bounds are given by a class of (convex) polytopes. The proposed control design method is based on the control Lyapunov function (CLF) theory introduced in [Artstein, 1983; Sontag, 1998]. Hence, we first recall, with parameters including those in [Lorenz, 1963], that these equations are point-dissipative, i.e. there is an explicit absorbing ball ${ℬ}$ given by the level set of a certain Lyapunov function, $V_{\infty }(x)$. However, since the minimum point of $V_{\infty }(x)$ does not coincide with any rest point of Lorenz system, we apply a modified solution to the “uniting CLF problem” (to unify local (possibly optimal) controls with global ones, proposed in [Andrieu & Prieur, 2010]) in order to obtain a CLF $V(x)$ for the affine system with minimum at a desired equilibrium point. Finally, we achieve the GAS of “any” rest point of this system via bounded and regular feedback controls by using the proposed CLF method, which also contains the following controllers: (i) damping controls outside ${ℬ}$, so they collaborate with the beneficial stable free dynamics, and (ii) (possibly optimal) feedback controls inside ${ℬ}$ that stabilize the control system at “any” desired rest point of the (unforced) Lorenz system.

In this paper, a diffusive zooplankton–phytoplankton model with time delay and hunting cooperation is established. First, the existence of all positive equilibria and their local stability are proved when the system does not include time delay and diffusion. Then, the existence of Hopf bifurcation at the positive equilibrium is proved by taking time delay as the bifurcation parameter, and the direction of Hopf bifurcation and the stability of bifurcating periodic solutions are investigated by using the center manifold theorem and the normal form theory in partial differential equations. Next, according to the theory of Turing bifurcation, the conditions for the occurrence of Turing bifurcation are obtained by taking the intraspecific competition rate of the prey as the bifurcation parameter. Furthermore, the corresponding amplitude equations are discussed by using the standard multi-scale analysis method. Finally, some numerical simulations are given to verify the theoretical results.

We study canard cycles and their cyclicity of a fast–slow Leslie–Gower predator–prey system with Allee effect. More specifically, we find necessary and sufficient conditions of the exact number (zero, one or two) of positive equilibria of the slow–fast system and its location (or their locations), and then we further completely determine its (or their) dynamics under explicit conditions. Besides, by geometric singular perturbation theory and the slow–fast normal form, we find explicit sufficient conditions to characterize singular Hopf bifurcation and canard explosion of the system. Additionally, the cyclicity of canard cycles is completely solved, and of particular interest is that we show the existence and uniqueness of a canard cycle, whose cyclicity is at most two, under corresponding precise explicit conditions.

In this paper, the impact of the strong Allee effect and ratio-dependent Holling–Tanner functional response on the dynamical behaviors of a predator–prey system is investigated. First, the positivity and boundedness of solutions of the system are proved. Then, stability and bifurcation analysis on equilibria is provided, with explicit conditions obtained for Hopf bifurcation. Moreover, global dynamics of the system is discussed. In particular, the degenerate singular point at the origin is proved to be globally asymptotically stable under various conditions. Further, a detailed bifurcation analysis is presented to show that the system undergoes a codimension-$1$ Hopf bifurcation and a codimension-$2$ cusp Bogdanov–Takens bifurcation. Simulations are given to illustrate the theoretical predictions. The results obtained in this paper indicate that the strong Allee effect and proportional dependence coefficient have significant impact on the fundamental change of predator–prey dynamics and the species persistence.