# Nonlinear Dynamics in a Chemical Reaction under an Amplitude-Modulated Excitation: Hysteresis, Vibrational Resonance, Multistability, and Chaos

Authors: A. V. Monwanou; A. A. Koukpémèdji; C. Ainamon; P. R. Nwagoum Tuwa; C. H. Miwadinou; J. B. Chabi Orou
Journal: Complexity
 (2020)

---

## Abstract

This paper deals with the effects of an amplitude-modulated (AM) excitation on the nonlinear dynamics of reactions between four molecules. The computation of the ﬁxed points of the autonomous nonlinear chemical system has been made in detail using the Cardan’s method. Hopf bifurcation has been also successfully checked. Routes to chaos have been investigated through bifurcations structures, Lyapunov exponent, phase portraits, and Poincaré section. The effects of the control force on chaotic motions have been strongly analyzed, and the control efficiency is found in the casesg=0(unmodulated case) and g≠0 with Ω=ω and Ω/w≠p/q; p and q are simple positive integers. Vibrational resonance (VR), hysteresis, and coexistence of several attractors have been studied in detail based on the relationship between the frequencies of the AM force. Results of analytical investigations are validated and complemented by numerical simulations.


---

## Body

## 1. Introduction

Nonlinear dynamics is a multidisciplinary ﬁeld that covers not only mathematics but also engineering, physics, chemistry, and biology. A monograph, authored by Strogatz [1], would be the ﬁrst place to start if one is interested in learning about this topic. In the past decade, the study of complex dynamics and chemistry of oscillating reaction under the inﬂuence of external perturbation received much attention and various results such as period-doubling bifurcation leading to chaotic motion, quasiperiodic route to chaos, coexistence of multiple attractors, hysteresis, and vibrational resonance have been obtained [2–12]. The chemical oscillating reactions in a continuously stirred tank reactor (CSTR) is one of the ﬁrst biochemical oscillations discovered. For considerable theoretical progress on the nature of chemical oscillation, the only known chemical oscillators were either biological in origin, like the glycolytic and oxidase-peroxidase systems; either discovered accidentally, like the Bray and BZ reactions; or variants of those reactions [2–12]. In these chemical oscillations, various dynamics behaviors are studied by many researchers. For instance, nonequilibrium phenomena such as oscillations, bistability, complex oscillations, and quasichaotic behavior of the reaction are revealed by these studies. One of the main challenges has been to predict and to control these phenomena in nonlinear chemical oscillations for potential applications (see [2–15]). The study of these oscillations has been made with a periodically external excitation or a parametric excitation [11, 12]. Moreover, the study of these oscillations is more advantageous for systems with the AM (amplitude modulation) [14]. However, some new dynamical phenomena including controllable frequency can be also presented. Binous and Bellagi [2] present the solution of four problems drawn from the chemical and biochemical engineering ﬁeld of study. These problems illustrate various important aspects of nonlinear dynamics such as limit cycles, quasiperiodic and chaotic behaviors, time series and phase portraits, power spectra, the time-delay reconstruction diagrams, Hopf bifurcation, bifurcation diagrams, and steady state multiplicity. Among a large number of examples, they have selected the following case studies because it illustrates the basic nonlinear dynamics concepts: the glycolytic oscillator model ﬁrst suggested by Selkov in 1968 in order to elucidate the mechanism that living cells use to obtain energy from sugar breakdown; the oregonator model derived by Field, Koros, and Noyes in the early 1970s, which elucidates the famous oscillatory behavior observed in the Belousov-Zhabotinsky (BZ) reaction; the steady state multiplicity in a biochemical reactor for both the Monod and substrate inhibition kinetics; and the three-variable autocatalator, ﬁrst proposed by Peng, Scott, and Showalter in 1990. Guruparan et al. [3] considered Brusselator chemical system driven by an amplitude-modulated (AM) force and studied numerically the dynamics of Brusselator chemical system driven by an amplitude-modulated force with widely different frequencies. They showed the occurrence of hysteresis and vibrational resonance and the coexistence of several period-T orbits, bifurcations of them, routes to chaos, and quasiperiodic and chaotic orbits. They have characterized periodic orbits, quasiperiodic orbits, chaotic orbits, hysteresis, and vibrational resonance using bifurcation diagram, maximal Lyapunov exponent, phase portrait, Poincaré map, and resonance plots. Shabunin et al. [5] introduced the lattice limit cycle (LLC) model as a minimal mean-ﬁeld scheme which can model reactive dynamics on lattices (low-dimensional supports) producing nonlinear limit cycle oscillations. They have found that, under the inﬂuence of an external periodic force, the dynamics of the LLC may be drastically modiﬁed. Synchronization phenomena, bifurcations, and transitions to chaos are also observed as a function of the excitation force. Taking advantage of the drastic change in the dynamics due to the periodic forcing, they found that it is possible to modify the output/product or the production rate of a chemical reaction at will, simply by applying a periodic force to it, without the need to change the support properties or the experimental conditions. Blekhman and Landa [6] considered resonances caused by a biharmonical external force with two different frequencies (the so-called vibrational resonances) using a bistable oscillator described by a Duffing equation as an example. It is shown that, in the case of a weakly damped oscillator, these resonances are conjugate; they occur as either the low and high frequency is varied. In addition, the resonances occur as the amplitude of the high-frequency excitation is varied. It is also shown that the high-frequency action induces the change in the number of stable steady states; these bifurcations are also conjugate and are the cause of the seeming resonance in an overdamped oscillator. Roy-Layinde et al. [7] examined the phenomenon of vibrational resonance (VR) and analyzed in a biharmonically driven two-ﬂuid plasma model with nonlinear dissipation. They derived analytically an equation for the slow oscillations of the system in themes of the parameters of the fast signal using the method of direct separation of motion. The presence of a high-frequency externally applied electric ﬁeld is found to signiﬁcantly modify the systems dynamics and, consequently, induce VR. They have identiﬁed the origin of the VR in the plasma model, not only from the effective plasma potential but also from the contributions of the effective nonlinear dissipation. Besides several dynamical changes, including multiple symmetry-breaking bifurcations, attractor escapes, and reversed period-doubling bifurcations, numerical simulations also revealed the occurrence of single and double resonances induced by symmetry-breaking bifurcations. Landa and McClintock [8] considered the effect of a high-frequency force on the response of a bistable system to a low-frequency signal for both the overdamped and weakly damped cases. They have shown that the response can be optimized by an appropriate choice of vibration amplitude. This vibrational resonance displays many analogies to the well-known phenomenon of stochastic resonance, but with the vibrational force ﬁlling the role usually played by noise. Jeevarathinam et al. [9] analyzed the vibrational resonance in the Duffing oscillator system in the presence of a gamma-distributed time-delayed feedback and an integrative time-delayed (uniformly distributed time delays over a ﬁnite interval) feedback. Particularly, applying a theoretical procedure, they obtain an expression for the response amplitude at the low frequency of the driving biharmonic force. For both double-well potential and single-well potential cases, they are able to identify the regions in parameter space where either two resonances, a single resonance, or no resonance occur. Theoretically predicted values of the response amplitude and the values of a control parameter at which resonance occurs are in good agreement with their numerical simulation. The analysis shows a strong inﬂuence of both types of time-delayed feedback on vibrational resonance. In the present paper, we seek hysteresis, vibrational resonance, and chaos in the system of reactions between four molecules when it is subjected to an external amplitude modulation excitation. More precisely, after an in-depth analysis of the ﬁxed points and of the Hopf bifurcation for the autonomous system, the effects of the modulated amplitude force on the dynamics of the chemical reaction considered have been studied in detail and the efficiency of the control force was analyzed. The paper is structured as follows: Section 2 gives the mathematical modeling of chemical model, while Section 3 analyses the ﬁxed points and their stability and the possibility to obtain the Hopf bifurcation. In Section 4, in-depth details of the bifurcation and the route to chaos when the system is under the periodically external excitation are presented. Section 5 deals with vibrational resonance, bifurcation, route to chaos, bistability, coexistence of attractors, and hysteresis for Ω≫ω. Section 6 analyses the effect of AM excitation when Ω=ω and Ω/w≠p/q, where p and q are simple positive integers. Finally, the conclusion of the research is given in Section 7.


## 2. Chemical Model and Its Mathematical Equations

In this work, we consider the autonomous system of reactions between four molecules as in [5, 16]:(1)2X+2Y⟶k13Y+S,X+S⟶k22Y,Y+S⟶k32S.X, Y, and S are molecules of the types X and Y and empty lattice sites, respectively, while k1,k2, and k3 are kinetic constants of the corresponding reactions. The complex oscillatory dynamics and the formation of spatial patterns observed experimentally have been predicted with great success by reactive multimolecular patterns [5, 12, 16–21]. Indeed, at the beginning of the 80 years, experimentally, “vacation models” were written by reactive quadrimolecular stages. In these models, increasing the number of vacant sites results in an increase in autocatalytic behavior, which results in oscillatory behavior. The reactions between nitrogen monoxide and carbon monoxide (NO + CO), nitrogen monoxide and ammonia (NO + NH3), and nitrogen monoxide and dihydrogen (NO + H2), which all use the platinum surface Pt as catalyst, are well-known models of vacancies [18]. In each of these reactions, an autocatalytic behavior linked to vacant sites is observed, and this explains, amongst other nonlinear phenomena, the ”surface explosion” phenomena observed frequently in heterogeneous catalysis representing narrow peaks of the product concentrations at regular temporal intervals [5, 12, 16–21].In the mean-ﬁeld approach, the kinetic equations for the evolution of the relative concentrations of the molecules (x, y) and of empty sites (S) can be written as follows [5, 16]:(2)dxdt=−2k1x2y2+k2xs,(3)dydt=k1x2y2−k3ys,(4)dsdt=2k1x2y2−k2xs+k3ys.It is easily shown that a condition of conservation is fulfilled:(5)dxdt+dydt+dsdt=0,and hence, x+y+s=cste. By choosing this constant equal to 1, we have x+y+s=1, and putting it in equation (3), the system is reduced to the following system:(6)dxdt=−2k1x2y2+k2x1−x−y,(7)dydt=k1x2y2−k3y1−x−y.The complex dynamic behaviors of this system under the action of a sinusoidal force have been studied in-depth [5, 16]. Recently, the importance of the amplitude modulated force has been proven for the control of complex dynamics of certain systems such as mechanical and electrical systems [14, 15] and biochemical and chemical systems [3, 11–13]. In this work, we consider the dynamics of molecules in the chemical reaction under the external amplitude-modulated excitation, which can be described by a system of two differential equations of order 1. So we modify the original model (7) of the LLC oscillator by adding an amplitude modulation force to the right themes of the ﬁrst equation. So, we have(8)dxdt=−2k1x2y2+k2x1−x−y+F+2GcosΩ0tsinw0t,(9)dydt=k1x2y2−k3y1−x−y,where F is the amplitude of unmodulated force, G is the degree of modulation, and w0 and Ω0 are the frequencies of the AMF.Using the time rescalingτ=k1t and notation x˙=dx/dτ, we reduce the system (8) as follows:(10)x˙=−2x2y2+αx1−x−y+f+2gcosΩτsinωτ(11)y˙=x2y2−βy1−x−y,with α=k2/k1>0;β=k3/k1>0,f=F/k1,g=G/k1,Ω=Ω0/k1, andw=w0/k1.


## 3. Equilibrium Points and Its Stability Analysis

We consider the system (10) in the absence of the excitation AM and we seek the ﬁxed points. Indeed, we have(12)−2x2y2+αx1−x−y=0,(13)x2y2−βy1−x−y=0.The resolution of the system formed by these two equations comes down to that of the following equation:(14)4β2α2y4−βy1−2βαy−y=0.The ﬁxed points of the autonomous system areE00,0 and E∗x∗,y∗ whose coordinates x∗,y∗ are obtained after solving the following equation:(15)4β2α2y3−β1−2βαy−y=0.To solve equation (15), we set p=α/2+α2/4β and q=−α2/4β.So, we have(16)y3+py+q=0.By taking the Cardan technique [11], the resolution of equation (16) amounts to looking for the roots of a trinome. If ∆ is the discriminant of this trinomial, we have (17)D=27Δ=4p3+27q2.Here,α > 0 and β > 0 so p3>0 and q2>0 and then D > 0.Equation (16) has only one real solution, which is T=−q±Δ/2. So, we have E3=x3,y3, where(18)x3=β2α1−δ3+β2α1+δ3,y3=α28β1−δ3+α28β1−δ3,with(19)δ=α+2β327αβ+1On the contrary, by combining the equations of the system (12), we have x+y=1.Thus, we haveE10,1 and E2=1,0.In total, the autonomous system admits exactly four ﬁxed pointsE00,0 (trivial), E10,1 (semitrivial), E21,0 (semitrivial), and E3x3,y3 (nontrivial).The study of the stability of each ﬁxed point is made by searching the Jacobian of the autonomous system associated. The Jacobian of the autonomous system associated with each ﬁxed pointE00,0 or E∗x∗,y∗ is as follows:(20)J=−4x∗y∗2+α1−x∗−y∗−αx∗−4x∗2y∗−αx∗2x∗y∗2+βy∗2x∗2y∗−β1−x∗−y∗+βy∗.ForE00,0, the eigenvalues of the Jacobian J are λ1=α and λ2=−β, which are reals and have opposite signs because α>0 and β>0. Hence, E0 is a saddle point. Then, the two eigenvalues of the Jacobian J associated with E10,1 are λ1=0 and λ2=β>0, and therefore, E1 is the unrobust point. On the contrary, for E21,0, the eigenvalues are λ1=0 and λ2=−α<0, and therefore, E2 is a robust point. Finally, for E3x3,y3, the characteristic equation giving the eigenvalues of J is(21)λ2+σ1λ+σ2=0,with σ1=−a1+a2 and σ2=a1a4−a2a3, where(22)a1=−4x3y32+α1−2x3−y3,a2=−4x32y3−αx3,a3=2x3y32+βy3,a4=2x32y3−β1−x3−2y3.Ifσ12=4σ2, then equation (21) has only one solution, and the ﬁxed point E3 is a center. For σ12>4σ2, the eigenvalues are (23)λ1=−σ1−σ12−4σ22,λ2=−σ1+σ12−4σ22.In this case, the equilibrium pointE3 is a node. If σ12<4σ2, then (24)λ1=−σ1−iσ12+4σ22,λ2=−σ1+i−σ12+4σ22,and E3 is either a focus or a center. By virtue of the Routh–Hurwitz criterion [22, 23], the ﬁxed point E3 is stable if and only if σ1>0and σ2>0 if not the ﬁxed point is unstable. We now look for the Hopf bifurcation assuming that λ=iww>0. Let us insert λ in equation (21), and we have (25)−w2+iσ1w+σ2=0,with w=σ2 and σ1=a1+a4=0.Assumingβ ﬁxed and taking α as a bifurcation parameter, differentiating both sides of the characteristic equation (21) associated to E3 with respect to αH, we obtain(26)dλdα=−12dσ1dα−12λdσ2dα,(27)Redλdα=−12dσ1dα,where α=αH and λ=iw. Thus, if for α=αH, σ1=0, and Redλ/dα≠0, then α=αH is a Hopf bifurcation value for the equilibrium point system E3. In short, we have the following Theorem 1.Theorem 1.
Ifβ is ﬁxed, system (8) in absence of the AM force has a Hopf bifurcation at equilibrium point E3x3,y3 when α passes through the critical value αH, which veriﬁes simultaneously the following relations (16 and 17): (28)α−β−12αα21−δβ1/3+βα21−δβ1/3−2αβ21−δα1/3+ββ21−δα1/3−α21−δβ2/3β21−δα1/3+α21−δβ1/3β21−δα2/3−12αα21+δβ1/3+βα21+δβ1/3−2α21−δβ1/3β21−δα1/3α21+δβ1/3+β21−δα2/3α21+δβ1/3−β21−δα1/3α21+δβ2/3−2αβ21+δα1/3+ββ21+δα1/3−α21−δβ2/3β21+δα1/3+2α21−δβ1/3β21−δα1/3β21+δα1/3−2α21−δβ1/3α21−δβ1/3β21+δα1/3+2α21+δβ1/3β21−δα1/3β21+δα1/3−α21+δβ2/3β21+δα1/3+β21+δα2/3α21−δβ1/3+β21+δα2/3α21+δβ1/3=0,(29)−4y32+4x3y3−2α+βdx3dα+2x32−8x3y3−α+2βdy3dα−2x3−y3+1≠0.To support this theorem, we represent the real part of the eigenvalues associated with the equilibrium pointE3x3,y3 in the case where the eigenvalues are complex as a function of α for fixed values of β. Figure 1 shows the critical values αH≈0.00845,αH≈0.0237,αH≈0.059,αH≈0.0875 of α corresponding to the Hopf bifurcation points H1,H2,H3 and H4 respectively. By a simple calculation, we show that each of these values verifies simultaneously (28) and (29). It should also be noted that the critical value αH for the Hopf bifurcation increases with β. At the moment of its birth, through Hopf bifurcation of the point E3x3,y3, the oscillations have infinitesimally small amplitude and near-harmonic shape (see Figure 2).Figure 1
Real part of the eigenvalues versusα for different values of β. H1,H2,H3 and H4are the points where Re (λ) = 0.Figure 2
Phase portraits of the autonomous system for parameter values corresponding to Hopf bifurcation at fixed pointE3=x3,y3 (a): αH′=0.00845; β=0.01 and (b): αH′=0.0237; β=0.027.
(a)(b)


## 4. Bifurcation and Route to Chaos forg = 0

Before studying the influence of the AM excitation on the dynamics of the system, let us analyze the effect of each parameterα, β, and f on the dynamics of the unmodulated system. Indeed, we represent the bifurcation diagram, the Lyapunov exponent, the phase space, and the Poincaré section taking as control parameters f and β. In the remains of this work, we take ω=0.0481, α>0, ß>0, x (0) = 0.5, and y (0) = 0.5. From all results of the simulations carried out in the case of the unmodulated system (g= 0), we note that the chemical reaction studied can present a very rich dynamic such as periodic, multiperiodic, quasiperiodic, and chaotic oscillations. In addition, we observe phenomena such as hysteresis, multistability, coexistence of periodic, multiperiodic, quasiperiodic, and chaotic attractors (see Figures 3–6). All these phenomena prove that very complex dynamics can be observed in the nonlinear chemical reaction studied. Thus, Figure 3 represents the bifurcation diagram and corresponding Lyapunov exponents for the system when f varies in the domain 0 ≤ f ≤ 0.003, the other parameters being fixed. From this figure, it emerges that the chemical oscillations are quasiperiodic, multiperiodic, and chaotic. When we compare the dynamics of the system when f increases from 0 to 0.003 (blue color) to that is obtained when f decreases from 0.003 to 0 (red color), we note apart from chemical oscillations of the same nature which are observed on a large domain that the system presents the coexistence of the attractors of period 6T with quasiperiodic attractors for 0.0011162 < f < 0.001131, attractors of period 5T with quasiperiodic attractors for 0.001521 < f < 0.00156103 and 0.001598 ≤ f < 0.0017485, attractors of period 5T with attractors of period 10T for 0.00156103 < f < 0.001598, and attractors of period 5T with chaotic attractors for 0.00298 < f ≤ 0.003. The chaotic behavior of chemical oscillations predicted by the bifurcation diagram and its corresponding Lyapunov exponent (Figure 3) is confirmed by the phase space (Figure 4(a)) and the corresponding Poincaré section (Figure 4(b)) obtained for f  = 0.0025. For f = 0.00165 and with two different initial conditions, we have plotted the phase space and the corresponding Poincaré section in Figure 5. We clearly note through this figure the coexistence of the quasiperiodic attractors (Figures 5(a) and 5(b)) with the attractors of period 5T (Figures 5(c) and 5(d)) thus justifying this phenomenon. To understand the influence of the nonlinear parameters α and β on the nature of the oscillations of the chemical reaction considered, we consider now β as the bifurcation parameter. Thus, Figure 6 shows the effect of the parameter α on the bifurcation diagram when β varies between 0.03 and 0.1. As it can be seen in this figure, the chaotic, periodic, multiperiodic, and quasiperiodic behaviors also depend on β. We take note that chaotic oscillations can only exist when 0.02 ≤ α < 0.05, and outside this domain, the chemical oscillations are periodic. It should be noted that the numerical simulations have shown that α and β have almost the same effects on the chemical dynamics of the studied reaction.Figure 3
Bifurcation diagram and its corresponding Lyapunov exponent of nonlinear chemical reaction withω  = 0.0481, α = 0.0413, β = 0.0667, and g = 0. Bifurcation diagrams and its corresponding Lyapunov exponent are obtained by scanning the parameter f upwards (blue) and downwards (red).
(a)(b)Figure 4
Chaotic phase portrait and its Poincaré section of nonlinear chemical reaction with parameters of Figure3 and f = 0.0025.
(a)(b)Figure 5
Phase space and its corresponding Poincaré section of nonlinear chemical reaction withf = 0.00165 and other parameters of Figure 3. (a, b) Quasiperiodic attractor for initial conditions (0.5, 0.5) and (c, d) attractor of period 5T for initial conditions (0.1, 0.1).
(a)(b)(c)(d)Figure 6
Effect ofα on bifurcation diagram of nonlinear chemical reaction with ω = 0.0481, f = 0.00165, and g = 0: (a) α = 0.02; (b) α = 0.03; (c) α = 0.05.
(a)(b)(c)


## 5. Effect of AM Force on the System withΩ ≫ ω

### 5.1. Vibrational Resonance and Amplitude Response

In a nonlinear dynamical system driven by a biharmonic force consisting of a low-frequencyω and a high-frequency Ω with Ωω, when the amplitude g of the high-frequency force is varied, the amplitude response at the low-frequency ω exhibits a resonance. This high-frequency force-induced resonance is called vibrational resonance [3, 6–9]. To determine the VR, we use the amplitude of the response at the frequency ω of the signal. Indeed, using the fourth-order Runge–Kutta algorithm, with time step size, we numerically integrate the system (10) of the chemical reaction under consideration. Thus, the numerical solution xτ allows to calculate the amplitude response Q through the following formula:(30)Q=Qs2+Qc2f,where(31)Qs=2nπ∫0nTxτsinwτdτ,Qc=2nπ∫0nTxτcoswτdτ,with T=2π/ω the response period and n = 500.We computeQ with a low-frequency force only, a high-frequency only, and with both forces.Forf=0, Q is determined as Qs2+Qc2, Qs and Qc and representing the Fourier coefficients of the output signal at the frequency 2πT and Q the amplitude of the response to this same frequency. Indeed, for Ω = 80ω with ω being small, the different results are shown in Figures 7–9. Figures 7(a) and 7(b) represent the variation in Q depending, respectively, on f for g = 0 and g for f = 0 when ω = 0.0481. From Figure 7(a), we note the presence of several peaks, and while f = 0, we have a single maximum for which g = 0.003301 and Q = 0.00092802. When g varies from 0 to 0.05 with f 6 = 0 and Ω = 80ω, the VR appears and we observe that its persistence, its form, its maximum amplitude, and the value of g at the resonance depends not only on f and ω but also on nonlinear parameters α and β of the nonlinear chemical action considered. Indeed, a multiresonance appears for ω < 0.06 (see Figure 8(a)) and disappears when 0.06 ≤ ω < 0.09. For 0.09 ≤ ω ≤ 0.1, the VR becomes a double resonance with two different Qmax (see Figure 8(b)) and disappears again when 0.1 < ω ≤ 0.14 and then reappears in the form of a monoresonance when ω > 0.14 (see Figure 8(c)). Finally, we note that when g varies from 0 to 0.05, the maximum amplitude Qmax decreases when the frequency ω increases.Figure 7
(a) Variation in numerically computedQ against the control parameter f with ω  = 0.0481, α = 0.0413, β = 0.0667, and g = 0; (b) variation in numerically computed Q against the control parameter g with ω  = 0.0481, Ω=80ω, α = 0.0413, β = 0.0667, and f = 0.
(a)(b)Figure 8
Effect of frequencyω on response amplitude Q versus g with ω = 0.0481, Ω = 80ω, α = 0.0413, β = 0.0667, and f = 0.001.
(a)(b)(c)Figure 9
(a) Effect of amplitudef on response amplitude Q versus g with ω = 0.016, Ω = 80ω, α = 0.0413, and β = 0.0667; (b) effect of nonlinear parameter α on response amplitude Q versus g with ω = 0.016, Ω=80ω, β = 0.0667, and f = 0.001; (c) effect of nonlinear parameter β on response amplitude Q versus g with ω = 0.016, Ω = 80 ω, α = 0.0413, and f = 0.001.
(a)(b)(c)Figure9 shows when g varies from 0 to 0.05, Ω  = 80ω with ω=0.16, and the influence of the parameters f, α, are β (Figures 9(a)–9(c), respectively) is shown on the VR. From Figure 9(a), we note that the maximum value of the amplitude of the response decreases when f increases and we obtain a double resonance when f = 0.0011. Figure 9(b) shows that the parameter α also influences the vibrational resonance. Finally, we notice through Figure 9(c) that VR exists when 0.06 < β < 0.08 and that the parameter β also has the same effects on VR as in the case of α. In conclusion, the birth and the disappearance of the vibrational resonance in the chemical reaction studied can be strongly controlled not only by the AM force but also by the parameters α and β of the chemical reaction.


### 5.2. Hysteresis, Coexistence of Attractors, and Multistability

In this section, we analyze the effect ofg on chemical dynamics when Ω≫ω by looking the hysteresis and coexistence of attractors phenomena. For this, we represent the bifurcation diagram by varying g with Ω=80ω (Figure 10).Figure 10
Bifurcation diagram and its corresponding Lyapunov exponent of nonlinear chemical reaction withω  = 0.0481, Ω  = 80ω, α = 0.0413, f = 0.002, and β = 0.0667. Bifurcation diagrams and its corresponding Lyapunov exponent are obtained by scanning the parameter g upwards (a) and downwards (b).
(a)(b)When we compare the dynamics of the system wheng increases from 0 to 0.1 (blue color, Figure 10(a)) to that obtained when g decreases from 0.1 to 0 (red color, Figure 10(b)), we note that the chemical oscillations of the same kind are observed but with amplitudes that outside x different especially for multiperiodical oscillations. This difference in amplitude noted shows that the chemical oscillations do not follow the same path as going back and forth when we increase g from 0 to 0.1 and when we decrease it from 0.1 to 0: this phenomenon is called hysteresis. There also needs to be a domain where multiperiodic attractors coexist with attractors of period 2T and attractors of period 4T on the one hand and also the coexistence on the other hand of the chaotic attractors with the attractors of period 4T and multiperiodic attractors. The phenomena of hysteresis and coexistence of attractors are very visible in Figure 10 (confirmed by Figure 11) and show multistable behaviors and hysteresis in the system. Thus, for Ω=80ω, the chemical oscillations become more complex for the reaction considered. Figure 12 shows the effect of β on the chemical oscillations and the two phenomena obtained. We deduce that when the parameter β ∈ [0, 0.06], the chemical oscillations are not chaotic and there is coexistence of periodic attractors and the presence of hysteresis.Figure 11
Poincaré section of nonlinear chemical reaction withg = 0.09 and other parameters of Figure 10: (a) chaotic attractor for initial conditions (0.1, 0.1); (b) attractor of period 4T for initial conditions (0.5, 0.5).
(a)(b)Figure 12
Effect ofβ on bifurcation diagram and its corresponding Lyapunov exponent of nonlinear chemical reaction with ω = 0.0481, Ω = 80 ω, α = 0.0413, and f = 0.002: (a) β = 0.04; (b) β = 0.08. Bifurcation diagrams and its corresponding Lyapunov exponent are obtained by scanning the parameter g upwards (blue) and downwards (red).
(a)(b)The chemical oscillations become very varied, and we note periodic, multiperiodic, quasiperiodic, and chaotic behaviors with the presence and persistence of hysteresis and the coexistence of multiple attractors whenβ ∈ [0.06, 0.7]. Finally, for β > 0.7, we observe a total dominance of chaotic oscillations and the disappearance of the hysteresis phenomenon. The parameters α and f have similar effects to that of β reason for which we decided not to represent the figures translating these effects. We notice that the parameter g and the frequencies of the AM force have a very important effect on the dynamics of the system as well as on the phenomena of hysteresis and coexistence of attractors.More precisely,g can be used to reduce the domains where the oscillations are quasiperiodic, the domains of coexistence of attractors, and to make disappear completely hysteresis.


## 5.1. Vibrational Resonance and Amplitude Response

In a nonlinear dynamical system driven by a biharmonic force consisting of a low-frequencyω and a high-frequency Ω with Ωω, when the amplitude g of the high-frequency force is varied, the amplitude response at the low-frequency ω exhibits a resonance. This high-frequency force-induced resonance is called vibrational resonance [3, 6–9]. To determine the VR, we use the amplitude of the response at the frequency ω of the signal. Indeed, using the fourth-order Runge–Kutta algorithm, with time step size, we numerically integrate the system (10) of the chemical reaction under consideration. Thus, the numerical solution xτ allows to calculate the amplitude response Q through the following formula:(30)Q=Qs2+Qc2f,where(31)Qs=2nπ∫0nTxτsinwτdτ,Qc=2nπ∫0nTxτcoswτdτ,with T=2π/ω the response period and n = 500.We computeQ with a low-frequency force only, a high-frequency only, and with both forces.Forf=0, Q is determined as Qs2+Qc2, Qs and Qc and representing the Fourier coefficients of the output signal at the frequency 2πT and Q the amplitude of the response to this same frequency. Indeed, for Ω = 80ω with ω being small, the different results are shown in Figures 7–9. Figures 7(a) and 7(b) represent the variation in Q depending, respectively, on f for g = 0 and g for f = 0 when ω = 0.0481. From Figure 7(a), we note the presence of several peaks, and while f = 0, we have a single maximum for which g = 0.003301 and Q = 0.00092802. When g varies from 0 to 0.05 with f 6 = 0 and Ω = 80ω, the VR appears and we observe that its persistence, its form, its maximum amplitude, and the value of g at the resonance depends not only on f and ω but also on nonlinear parameters α and β of the nonlinear chemical action considered. Indeed, a multiresonance appears for ω < 0.06 (see Figure 8(a)) and disappears when 0.06 ≤ ω < 0.09. For 0.09 ≤ ω ≤ 0.1, the VR becomes a double resonance with two different Qmax (see Figure 8(b)) and disappears again when 0.1 < ω ≤ 0.14 and then reappears in the form of a monoresonance when ω > 0.14 (see Figure 8(c)). Finally, we note that when g varies from 0 to 0.05, the maximum amplitude Qmax decreases when the frequency ω increases.Figure 7
(a) Variation in numerically computedQ against the control parameter f with ω  = 0.0481, α = 0.0413, β = 0.0667, and g = 0; (b) variation in numerically computed Q against the control parameter g with ω  = 0.0481, Ω=80ω, α = 0.0413, β = 0.0667, and f = 0.
(a)(b)Figure 8
Effect of frequencyω on response amplitude Q versus g with ω = 0.0481, Ω = 80ω, α = 0.0413, β = 0.0667, and f = 0.001.
(a)(b)(c)Figure 9
(a) Effect of amplitudef on response amplitude Q versus g with ω = 0.016, Ω = 80ω, α = 0.0413, and β = 0.0667; (b) effect of nonlinear parameter α on response amplitude Q versus g with ω = 0.016, Ω=80ω, β = 0.0667, and f = 0.001; (c) effect of nonlinear parameter β on response amplitude Q versus g with ω = 0.016, Ω = 80 ω, α = 0.0413, and f = 0.001.
(a)(b)(c)Figure9 shows when g varies from 0 to 0.05, Ω  = 80ω with ω=0.16, and the influence of the parameters f, α, are β (Figures 9(a)–9(c), respectively) is shown on the VR. From Figure 9(a), we note that the maximum value of the amplitude of the response decreases when f increases and we obtain a double resonance when f = 0.0011. Figure 9(b) shows that the parameter α also influences the vibrational resonance. Finally, we notice through Figure 9(c) that VR exists when 0.06 < β < 0.08 and that the parameter β also has the same effects on VR as in the case of α. In conclusion, the birth and the disappearance of the vibrational resonance in the chemical reaction studied can be strongly controlled not only by the AM force but also by the parameters α and β of the chemical reaction.


## 5.2. Hysteresis, Coexistence of Attractors, and Multistability

In this section, we analyze the effect ofg on chemical dynamics when Ω≫ω by looking the hysteresis and coexistence of attractors phenomena. For this, we represent the bifurcation diagram by varying g with Ω=80ω (Figure 10).Figure 10
Bifurcation diagram and its corresponding Lyapunov exponent of nonlinear chemical reaction withω  = 0.0481, Ω  = 80ω, α = 0.0413, f = 0.002, and β = 0.0667. Bifurcation diagrams and its corresponding Lyapunov exponent are obtained by scanning the parameter g upwards (a) and downwards (b).
(a)(b)When we compare the dynamics of the system wheng increases from 0 to 0.1 (blue color, Figure 10(a)) to that obtained when g decreases from 0.1 to 0 (red color, Figure 10(b)), we note that the chemical oscillations of the same kind are observed but with amplitudes that outside x different especially for multiperiodical oscillations. This difference in amplitude noted shows that the chemical oscillations do not follow the same path as going back and forth when we increase g from 0 to 0.1 and when we decrease it from 0.1 to 0: this phenomenon is called hysteresis. There also needs to be a domain where multiperiodic attractors coexist with attractors of period 2T and attractors of period 4T on the one hand and also the coexistence on the other hand of the chaotic attractors with the attractors of period 4T and multiperiodic attractors. The phenomena of hysteresis and coexistence of attractors are very visible in Figure 10 (confirmed by Figure 11) and show multistable behaviors and hysteresis in the system. Thus, for Ω=80ω, the chemical oscillations become more complex for the reaction considered. Figure 12 shows the effect of β on the chemical oscillations and the two phenomena obtained. We deduce that when the parameter β ∈ [0, 0.06], the chemical oscillations are not chaotic and there is coexistence of periodic attractors and the presence of hysteresis.Figure 11
Poincaré section of nonlinear chemical reaction withg = 0.09 and other parameters of Figure 10: (a) chaotic attractor for initial conditions (0.1, 0.1); (b) attractor of period 4T for initial conditions (0.5, 0.5).
(a)(b)Figure 12
Effect ofβ on bifurcation diagram and its corresponding Lyapunov exponent of nonlinear chemical reaction with ω = 0.0481, Ω = 80 ω, α = 0.0413, and f = 0.002: (a) β = 0.04; (b) β = 0.08. Bifurcation diagrams and its corresponding Lyapunov exponent are obtained by scanning the parameter g upwards (blue) and downwards (red).
(a)(b)The chemical oscillations become very varied, and we note periodic, multiperiodic, quasiperiodic, and chaotic behaviors with the presence and persistence of hysteresis and the coexistence of multiple attractors whenβ ∈ [0.06, 0.7]. Finally, for β > 0.7, we observe a total dominance of chaotic oscillations and the disappearance of the hysteresis phenomenon. The parameters α and f have similar effects to that of β reason for which we decided not to represent the figures translating these effects. We notice that the parameter g and the frequencies of the AM force have a very important effect on the dynamics of the system as well as on the phenomena of hysteresis and coexistence of attractors.More precisely,g can be used to reduce the domains where the oscillations are quasiperiodic, the domains of coexistence of attractors, and to make disappear completely hysteresis.


## 6. Effect of AM Force in the System withω = Ω

Here, the effect of the AM force is analyzed in the caseω = Ω. For this reason, we chose g as the bifurcation parameter, and the results obtained are shown in Figure 13. From the analysis of Figure 13, it appears that the variation in α can make disappear chaotic and quasiperiodic oscillations but preserves the phenomena of hysteresis and coexistence of attractors. In addition, the simulations (the figures of which are not shown here) have shown that the parameters β and f reduce the two phenomena and also act on the chemical oscillations of the model studied. Figure 14 represents the Poincaré section of nonlinear chemical reaction and shows the 3T-periodic, 4T-periodic, and chaotic oscillations as also obtained in Figure 13. Finally, we look for the effect of the AM force when the frequencies of the amplitude modulation force are not resonant and such that their relationship to each other is irrational.Figure 13
Effect of parameterα on bifurcation diagram of nonlinear chemical reaction with Ω = ω = 0.0481, f = 0.002, and β = 0.04: (a) α = 0.03; (b) α = 0.028; (c) α = 0.02. Bifurcation diagrams and its corresponding Lyapunov exponent are obtained by scanning the parameter g upwards (blue) and downwards (red).
(a)(b)(c)Figure 14
Effect ofα on Poincaré section of nonlinear chemical reaction with g = 0.002 and other parameters of Figure 13: (a) α = 0.03: attractor of period 3T; (b) α = 0.028: attractor of period 4T; (c) α = 0.02: chaotic attractor.
(a)(b)(c)For this, we takeω = 0.0481 and Ω = (5 −1)/2; we also seek the dynamics of the system by constructing the bifurcation diagram considering g as control parameter. Figures 15 and 16, respectively, represent the effects of β and f on the bifurcation diagram. When we analyze these figures, we note the same remarks as the case where ω = Ω except that here each of the parameters β and f can be used to reduce or even eliminate the chaotic chemical oscillations.Figure 15
Effect of parameterβ on bifurcation diagram of nonlinear chemical reaction with ω = 0.0481, Ω = (5−1)/2, f = 0.002, and α = 0.28: (a) β = 0.04; (b) β = 0.06; (c) β = 0.07. Bifurcation diagrams and its corresponding Lyapunov exponent are obtained by scanning the parameter g upwards (blue) and downwards (red).
(a)(b)(c)Figure 16
Effect of parameterf on bifurcation diagram of nonlinear chemical reaction with ω = 0.0481, Ω = (5−1)/2, α = 0.28, and β = 0.04: (a) f = 0; (b) f = 0.002; (c) f = 0.0025. Bifurcation diagrams and its corresponding Lyapunov exponent are obtained by scanning the parameter g upwards (blue) and downwards (red).
(a)(b)(c)


## 7. Conclusions

In this work, we have studied the influence of an amplitude-modulated force on the chemical oscillations of reactions between four molecules. We have looked for fixed points and their natures when the system is autonomous. It appears that the autonomous system admits exactly four fixed points. The nature of the nontrivial fixed point depends of parametersα and β, and Hopf bifurcation is obtained for these two parameters for the nontrivial equilibrium point. The dynamics of the chemical system considered has been widely analyzed in the unmodulated case, and periodic, multiperiodic, quasiperiodic, and chaotic chemical oscillations have been obtained. Domains of existence of hysteresis, of coexistence of quasiperiodic attractors with attractors of period 5T and attractors of period 6T, and of coexistence of multiperiodic attractors with chaotic attractors have been obtained. From the simulations made, it should be noted that, for certain values of the nonlinear parameters α and β, the complex dynamics can be controlled and even reduced to a periodic oscillation. For g 6 = 0 and for Ωω, complex phenomena such as vibrational resonance (VR), hysteresis, and coexistence of attractors appear and are very remarkable and can be well controlled either by the AM force or by the parameters α and β. The domains of existence of these different phenomena are less important when Ω = ω. When the frequencies of the amplitude modulation force are not resonant and are such that their relationship to each other is irrational, the multiperiodic attractors when they exist are of period ≥ 4T and their domains of existence are more important. Finally, our work confirmed the chaotic oscillation and the Hopf bifurcation obtained in [5] and proved that the appearance of VR, hysteresis, multistability, and coexistence of various attractors in the chemical reaction is considered. The presence of VR shows that the amplitude of oscillations of concentrations can be explored, while the hysteresis phenomenon indicates that the amplitude of oscillations of concentrations will not be able to follow the same path of evolution when the control parameter increases or decreases following the same path. The multistability and the coexistence of attractors obtained reveal that the stable or unstable nature of the chemical oscillations and of the chemical attractors in the same domain of the parameters is considered. The use of the amplitude modulated force has shown that the various complex phenomena obtained can be controlled.



---

*Source: 8823458-2020-10-05.xml*