删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

Superperiodicity, chaos and coexisting orbits of ion-acoustic waves in a four-component nonextensive

本站小编 Free考研考试/2022-01-02

Asit Saha,1, Jharna Tamang,1, Guo-Cheng Wu2, Santo Banerjee,3,41 Department of Mathematics, Sikkim Manipal Institute of Technology, Sikkim Manipal University, Majitar, Rangpo, East-Sikkim 737136, India
2 Data Recovery Key Laboratory of Sichuan Province, College of Mathematics and Information Science, Neijiang Normal University, Neijiang 641100, China
3 Department of Mathematical Sciences, "Giuseppe Luigi Lagrange" Politecnico di Torino, Corso Duca degli Abruzzi 24, Torino, Italy

First author contact: 4 Author to whom any correspondence should be addressed.
Received:2020-03-9Revised:2020-06-9Accepted:2020-06-11Online:2020-10-15


Abstract
Superperiodicity, chaos and coexisting orbits of ion-acoustic waves (IAWs) are studied in a multi-component plasma consisting of fluid ions, q -nonextensive cold and hot electrons and Maxwellian hot positrons. The significant impacts of the system parameters on superperiodic and nonlinear periodic IAWs are presented. Considering an external periodic perturbation various types of quasiperiodic and chaotic features for IAWs are studied in different parametric ranges through time series’ plots, phase spaces and Lyapunov exponents. It has been observed that there exist some coexisting orbits for IAWs. Coexisting orbits for IAWs in a classical electron–positron–ion plasma system are reported.
Keywords: dynamical system;supernonlinear periodic wave;quasiperiodicity;multistability


PDF (1255KB)MetadataMetricsRelated articlesExportEndNote|Ris|BibtexFavorite
Cite this article
Asit Saha, Jharna Tamang, Guo-Cheng Wu, Santo Banerjee. Superperiodicity, chaos and coexisting orbits of ion-acoustic waves in a four-component nonextensive plasma. Communications in Theoretical Physics[J], 2020, 72(11): 115501- doi:10.1088/1572-9494/aba256

1. Introduction

Our early universe is composed of completely ionized gases, which comprised of opposite charge particles of same masses and charge densities [1]. Negatively charged particles (electrons) and positively charged particles (positrons) are considered to be the fundamental composites of ionized gases. Electron–positron particles are main constituents of many astrophysical surroundings, such as, pulsar magnetospheres [2], active galactic nuclei [3], solar flares [4], fireballs producing γ -ray bursts [5], etc. Plasmas having electron and positron are detected in laboratory examinations [68]. Electron–positron pair in plasmas may be generated during intense short laser pulse transmission [9]. Many astrophysical [4, 5] and laboratory [1013] observation shows that plasmas are composed of electrons, ions and positrons. Surko and Murphy [10] studied that for a large scale of parameters, desolation of electrons and positrons occur when 1×1012 cm−3 is considered as density of electrons with temperature of 1 eV and the desolation of positron takes more than 1 s. Hence, the examination of electron–positron–ion (EPI) plasmas is beneficial to perceive features of astrophysical and laboratory environments. Lately, transmissions of different waves in plasmas bearing composition of EPI have gained immense attention of researchers [1419].

Recently, investigations of supernonlinear waves (SNWs) [20, 21] were initiated in plasmas. Dubinov et al [22] examined the properties of SNWs in multi-component plasma environments. The authors reported that SNWs can exist in plasma environments when there are at least three different plasma components. The multi-component plasma environment shows complexity which prompts existence of SNW features. In the form of plasma waves, SNWs exist in different physical environments such as, electrostatic and MHD waves [23]. All SNWs are greatly anharmonic waves with a substantially large amplitude [24]. The studies of solitary solution of SNWs known as supersolitons were reported earlier [2528], but were not specified as supersolitons. Verheest et al [2931] studied the existence of supersolitons in different three-component plasma systems. Singh et al [32] studied the presence of positive potential supersolitons with nonlinear potential solitons. Rufai et al [33] investigated a plasma limitation suggesting that the plasma system should have at least three-component for the transmission of SNWs. Lakhina et al [34] examined plasma system with three-components namely, two kinds of ions and Maxwellian electrons, observed absence of SNWs. Therefore, the necessary condition for SNWs to restrain in a plasma system is composition of three or more plasma components. Recently, El-Wakil et al [35] reported transmission of SNWs in four-constituent nonextensive dusty plasmas.

The kinetic motions of a charged particle are governed by Maxwellian for a plasma system exhibiting low range of interactions and collisions. However, due to existence of stationary states in some system, charged particles of the system tend to diverge from the standard Maxwellian distribution and converge to non-Maxwellian distributions. Some of non-Maxwellian distribution are q -nonextensive, kappa, nonthermal distributions, etc. Tsallis [36] proposed the Boltzmann–Gibbs–Shannon q -nonextensive entropy. In plasmas, q- nonextensive distribution has huge implication on plasma constituents that are widely used in astrophysical and cosmological environments such as, collisionless thermal plasma, hadronic matter, stellar polytropes, quark-gluon plasma, etc. The strength of nonextensivity q varies as −1<q <1 and q >1. It is not normalizable for q <−1 and tends to Maxwellian distribution for $q\to 1$ . Tsallis [37, 38] established the nonextensive distribution for systems that show extensive interactivity and collisions, and are not definable under the traditional Boltzmann–Gibbs statistics. The nonextensitivity of combined (A +B) is projected by ${S}_{q}^{(A+B)}={S}_{q}^{(A)}+{S}_{q}^{(B)}+(1-q){S}_{q}^{(A)}{S}_{q}^{(B)}$, where A and B are two distinct systems. This formulation classify the basis of system dynamics generalized by Tsallis. The entropy generalized by Tsallis is associated to plasma dynamics [39, 40], Hamiltonian systems [41] with extensive interaction, and nonlinear gravitational model [42]. Moreover, Liu et al [43] and Lima et al [39] studied plasma systems that follow velocity distribution other than Maxwellian. Many works [4448] were examined to study the Tsallis q -entropy and to enhance the theory of generalized statistics.

Nonlinear dynamics and chaotic phenomenon have been widely applied in various fields of engineering and science such as, biomedical engineering [49, 50], communication and cryptography [5155]. Of late, it has been investigated that dynamical systems can produce a series of different stable features for a particular set of parameters. The corresponding phenomena is defined as the existence of coexisting attractors in the system. Coexisting attractors have been observed in several real world phenomenon [5661]. However, the study of coexisting attractors in plasmas, as a new area of research in a chaotic system, is still in its beginning. The multistability is not reported much in plasma systems. In the work [62], it has been investigated that under the evidence of external forcing, the butterfly pattern can be distorted in a plasma perturbation model. In this article, we have proposed that the nonautonomous system can produce chaos and various coexisting attractors, which are quantified using phase space and Lyapunov exponents.

Lately, the application of bifurcation of planar dynamical systems [63] has enriched the study of underlying dynamics of plasma systems. Many physical systems are constantly in dynamic motions, application of bifurcation concept helps to analyze transitions and instabilities occurring in the system. The best possible way to describe the bifurcation analysis is through phase plane analysis which can be controlled by system parameters [63]. The bifurcation theory was employed in plasma for the first to describe nonlinear waves by Samanta et al [64]. Employing bifurcation theory, many researchers [65, 66] investigated various acoustic waves in plasmas. The features of supersolitons [6770] were examined using bifurcation concept in multi-component plasmas. Very recently, few works [7173] were published on supernonlinear periodic waves in plasmas. Tamang and Saha [73] reported positron-acoustic superperiodic waves and its dynamic motions in plasmas. Recently, Saha et al [74] reported bifurcation and coexisting attractors under the dynamics of Schrödinger equation in quantum plasmas. However, coexisting orbits of ion-acoustic waves (IAWs) using the planar dynamical systems are not reported in classical EPI plasma as far as we have surveyed. Therefore, the main objective of this study is to present superperiodicity, chaos and coexisting orbits of IAWs in a multi-component nonextensive plasma containing cold fluid ions, cold electrons, hot electrons and positrons.

The article is systematized as follows: in section 2, plasma model is considered. In section 3, conservative dynamical system (CDS) and superperiodic wave feature are obtained. In section 4, quasiperiodic and chaotic orbits under external perturbation are examined. In section 5, coexisting orbits are shown. Section 6 provides conclusions of the work.

2. Model equations

We inspect the propagation of IAWs in a multi-component plasma system containing fluid ions, Maxwellian hot positrons and q -nonextensive cold and hot electrons. The following normalized model equations [75] represent the dynamics of IAWs:$ \begin{eqnarray}\displaystyle \frac{\partial n}{\partial t}+\displaystyle \frac{\partial }{\partial x}({nv})=0,\end{eqnarray}$$ \begin{eqnarray}\displaystyle \frac{\partial v}{\partial t}+v\displaystyle \frac{\partial v}{\partial x}=-\displaystyle \frac{\partial \phi }{\partial x},\end{eqnarray}$$ \begin{eqnarray}\displaystyle \frac{{\partial }^{2}\phi }{\partial {x}^{2}}={n}_{h}+{n}_{c}-\alpha {n}_{p}-(1-\alpha )n,\end{eqnarray}$ where ${n}_{p}={{\rm{e}}}^{-\gamma \phi }$ is number density of hot positrons. Here, n and nh (nc ) are number densities of ions and hot electrons (cold electrons), respectively. Here, ion velocity and electrostatic potential are denoted by v and φ, respectively. γ is given by ${T}_{\mathrm{eff}}/{T}_{p}$ with effective temperature ${T}_{\mathrm{eff}}={T}_{c}/({\mu }_{c}+{\mu }_{h}\beta )$, where Th, Tc and Tp denote temperatures of hot electrons, cold electrons and hot positrons, respectively. An equilibrium density ratio α is specified by $\left(\tfrac{{n}_{{\rm{p}}0}}{{n}_{{\rm{e}}0}}\right)$, where ‘0’ is used for unperturbed quantity. Here, nh, nc and np are normalized by ne0 . v is normalized by ${C}_{s}={({T}_{\mathrm{eff}}/{m}_{i})}^{1/2}$ . Here, x and t are normalized by ${\lambda }_{D}={({T}_{\mathrm{eff}}/4\pi {n}_{{\rm{e}}0}{e}^{2})}^{1/2}$ and ${\omega }_{\mathrm{pi}}^{-1}={(4\pi {n}_{{\rm{e}}0}{e}^{2}/{m}_{i})}^{-1/2}$, respectively, where mi is mass of ion and electronic charge e .

The nonextensive velocity distribution of electrons is achieved by considering the following function [76]:
$ \begin{eqnarray*}{f}_{e}(v)={C}_{q}{\left\{1+(q-1)\left[\displaystyle \frac{{m}_{e}{v}^{2}}{2{k}_{{\rm{B}}}{T}_{e}}-\tfrac{e\phi }{{k}_{{\rm{B}}}{T}_{e}}\right]\right\}}^{\tfrac{1}{(q-1)}},\end{eqnarray*}$
where the distribution fe (v) satisfies the laws of thermodynamics and prosper the Tsallis entropy (q). Here, kB is the Boltzmann constant, Te, me and v are temperature, mass and velocity of electrons, respectively and Cq is the normalizing constant expressed as
$ \begin{eqnarray*}{C}_{q}={n}_{{\rm{e}}0}\displaystyle \frac{{\rm{\Gamma }}\left(\tfrac{1}{1-q}\right)}{{\rm{\Gamma }}\left(\tfrac{1}{1-q}-\tfrac{1}{2}\right)}\sqrt{\displaystyle \frac{{m}_{e}(1-q)}{2\pi {k}_{{\rm{B}}}{T}_{e}}}\ \mathrm{for}\ -1\lt q\lt 1,\end{eqnarray*}$
and
$ \begin{eqnarray*}{C}_{q}={n}_{{\rm{e}}0}\displaystyle \frac{1+q}{2}\tfrac{{\rm{\Gamma }}\left(\tfrac{1}{q-1}+\tfrac{1}{2}\right)}{{\rm{\Gamma }}\left(\tfrac{1}{q-1}\right)}\sqrt{\displaystyle \frac{{m}_{e}(q-1)}{2\pi {k}_{{\rm{B}}}{T}_{e}}}\ \mathrm{for}\ q\gt 1.\end{eqnarray*}$
Performing integration on fe (v), electron number density with nonextensive distribution is acquired as:
$ \begin{eqnarray*}{n}_{e}={n}_{{\rm{e}}0}{\left\{1+(q-1)\tfrac{e\phi }{{k}_{{\rm{B}}}{T}_{e}}\right\}}^{1/(q-1)+1/2}.\end{eqnarray*}$
Thus, the normalized cold and hot electron number densities can be expressed as$ \begin{eqnarray}{n}_{c}={\mu }_{c}{\left\{1+\displaystyle \frac{1}{{\mu }_{c}+{\mu }_{h}\beta }(q-1)\phi \right\}}^{\tfrac{1}{q-1}+\tfrac{1}{2}},\end{eqnarray}$$ \begin{eqnarray}{n}_{h}={\mu }_{h}{\left\{1+\displaystyle \frac{\beta }{{\mu }_{c}+{\mu }_{h}\beta }(q-1)\phi \right\}}^{\tfrac{1}{q-1}+\tfrac{1}{2}},\end{eqnarray}$ where q >−1, and it measures strength of nonextensivity. Here, $\beta =\tfrac{{T}_{c}}{{T}_{h}}$, ${\mu }_{c}=\tfrac{{n}_{{\rm{c}}0}}{{n}_{{\rm{e}}0}}$ and ${\mu }_{h}=\tfrac{{n}_{{\rm{h}}0}}{{n}_{{\rm{e}}0}}$ . The charge-neutrality condition is given by μc +μh =1.

3. CDS and superperiodic waves

It is effective to obtain all possible nonlinear waves of the considered plasma system. To study this, the model equations are changed into a CDS [7779] employing the wave alteration $\xi =x-{Vt}$, where the speed of wave is denoted by V . Substituting this alteration to equation (1 ) and performing integration subjected to conditions n =1, v =0, as $\xi \to \pm \infty $, one can deduce$ \begin{eqnarray}n=\displaystyle \frac{V}{V-v}.\end{eqnarray}$ Similarly, from equation (2 ) we obtain the following expression subjected to conditions v =0, φ =0, as $\xi \to \pm \infty $ ,$ \begin{eqnarray}V-v=\sqrt{{V}^{2}-2\phi }.\end{eqnarray}$ From equations (6 ) and (7 ), we get$ \begin{eqnarray}n=\displaystyle \frac{V}{\sqrt{{V}^{2}-2\phi }}.\end{eqnarray}$ Substituting equations (4 ), (5 ) and (8 ) in equation (3 ), we deduce$ \begin{eqnarray}\begin{array}{rcl}\displaystyle \frac{{{\rm{d}}}^{2}\phi }{{\rm{d}}{\xi }^{2}} & = & {\mu }_{h}{\left\{1+\displaystyle \frac{\beta }{{\mu }_{c}+{\mu }_{h}\beta }(q-1)\phi \right\}}^{\tfrac{1}{q-1}+\tfrac{1}{2}}\\ & & +{\mu }_{c}{\left\{1+\displaystyle \frac{1}{{\mu }_{c}+{\mu }_{h}\beta }(q-1)\phi \right\}}^{\tfrac{1}{q-1}+\tfrac{1}{2}}\\ & & -\alpha {{\rm{e}}}^{-\gamma \phi }-(1-\alpha )\displaystyle \frac{V}{\sqrt{{V}^{2}-2\phi }}.\end{array}\end{eqnarray}$ We rewrite the above equation considering up to 4th degree term as$ \begin{eqnarray}\displaystyle \frac{{{\rm{d}}}^{2}\phi }{{\rm{d}}{\xi }^{2}}=a\phi +b{\phi }^{2}+c{\phi }^{3}+d{\phi }^{4},\end{eqnarray}$ where
$ \begin{eqnarray*}\begin{array}{rcl}a & = & \displaystyle \frac{q+1}{2}+\alpha \gamma -\displaystyle \frac{1}{{V}^{2}}(1-\alpha ),\\ b & = & \displaystyle \frac{{\mu }_{h}{\beta }^{2}+{{mu}}_{c}}{{\left({\mu }_{h}\beta +{{mu}}_{c}\right)}^{2}}\displaystyle \frac{(q+1)(3-q)}{8}-\displaystyle \frac{\alpha {\gamma }^{2}}{2}-\displaystyle \frac{3}{2{V}^{4}}(1-\alpha ),\end{array}\end{eqnarray*}$
$ \begin{eqnarray*}\begin{array}{rcl}c & = & \displaystyle \frac{{\mu }_{h}{\beta }^{3}+{{mu}}_{c}}{{\left({\mu }_{h}\beta +{{mu}}_{c}\right)}^{3}}\displaystyle \frac{(q+1)(3-q)(5-3q)}{48}\\ & & +\displaystyle \frac{\alpha {\gamma }^{3}}{6}-\displaystyle \frac{5}{2{V}^{6}}(1-\alpha ),\end{array}\,\end{eqnarray*}$
$ \begin{eqnarray*}\begin{array}{rcl}d & = & \displaystyle \frac{{\mu }_{h}{\beta }^{4}+{{mu}}_{c}}{{\left({\mu }_{h}\beta +{{mu}}_{c}\right)}^{4}}\displaystyle \frac{(q+1)(3-q)(5-3q)(7-5q)}{384}\\ & & -\displaystyle \frac{\alpha {\gamma }^{4}}{24}-\displaystyle \frac{35}{8{V}^{8}}(1-\alpha ).\end{array}\,\end{eqnarray*}$
Equation (10 ) can be framed into a dynamical system as:$ \begin{eqnarray}\left\{\begin{array}{l}\tfrac{{\rm{d}}\phi }{{\rm{d}}\xi }=z,\\ \tfrac{{\rm{d}}z}{{\rm{d}}\xi }=a\phi +b{\phi }^{2}+c{\phi }^{3}+d{\phi }^{4}.\end{array}\right.\end{eqnarray}$ The Hamiltonian function H (φ, z) of the system (11 ) is given by$ \begin{eqnarray}H(\phi ,z)=\displaystyle \frac{{z}^{2}}{2}-\displaystyle \frac{a{\phi }^{2}}{2}-\displaystyle \frac{b{\phi }^{3}}{3}-\displaystyle \frac{c{\phi }^{4}}{4}-\displaystyle \frac{d{\phi }^{5}}{5}.\end{eqnarray}$ If we consider $\vec{F}=(z,a\phi +b{\phi }^{2}+c{\phi }^{3}+d{\phi }^{4}),$ then$ \begin{eqnarray}\vec{{\rm{\nabla }}}\cdot \vec{F}=\displaystyle \frac{\partial }{\partial \phi }(z)+\displaystyle \frac{\partial }{\partial z}(a\phi +b{\phi }^{2}+c{\phi }^{3}+d{\phi }^{4})=0.\end{eqnarray}$ Equation (13 ) shows that the system (11 ) is a CDS.

To achieve different qualitative phase portrait profiles, we first obtain all equilibrium points of the dynamical system (11 ) by $\tfrac{{\rm{d}}\phi }{{\rm{d}}\xi }=0$ and $\tfrac{{\rm{d}}z}{{\rm{d}}\xi }=0$, which give$ \begin{eqnarray}z=0\,\mathrm{and}\,\phi (d{\phi }^{3}+c{\phi }^{2}+b\phi +a)=0,\end{eqnarray}$$ \begin{eqnarray}\Rightarrow \,z=0\,\mathrm{and}\,\phi ({\phi }^{3}+m{\phi }^{2}+n\phi +o)=0,\end{eqnarray}$ while $m=\tfrac{c}{d},\,n=\tfrac{b}{d}$ and $\,o=\tfrac{a}{d}$ . The CDS (11 ) has four equilibrium points at E0 (φ0 ,0), E1 (φ1 ,0), E2 (φ2 ,0), and E3 (φ3 ,0), while φ0 =0, and φ1, φ2, φ3 are real roots of the equation ${\phi }^{3}+m{\phi }^{2}+n\phi +o=0$ . The determinant of the Jacobian matrix at Ei (φi, 0), i =0, 1, 2, 3, of the CDS (11 ) is obtained as$ \begin{eqnarray}D({\phi }_{i},0)=-(a+2b{\phi }_{i}+3c{\phi }_{i}^{2}+4d{\phi }_{i}^{3}).\end{eqnarray}$ Using the concept of planar dynamical systems [78, 79], we attain saddle point and center about an equilibrium point of the CDS (11 ) for D <0 and D >0, respectively. In figure 1, we present all possible phase profiles of the CDS (11 ) depending on different possible sets of values of system parameters q α, γ, μc, β and V . Figures 1 (a) and (b) show phase profiles of the CDS (11 ) for $q=-0.17$ and q =−0.26, respectively, with $\alpha =0.9,$γ =0.01, μc =0.06, β =0.05 and V =0.5. The phase profile in figure 1 (a) shows saddle point at the origin and as we increase the value of q in the range −1<q <0 keeping other parameters fixed, we observe that the nature of stability at equilibrium points of the system (11 ) changes. Therefore, we observe a center at the origin in figure 1 (b). In this case, we find that there are equilibrium points featuring pairs of centers and saddle points, and four kind of orbits namely, nonlinear homoclinic and periodic orbits, supernonlinear periodic and homoclinic orbits of the CDS (11 ). Any particular orbit in phase space is associated to a particular wave. Therefore, IAWs in the EPI plasma system has nonlinear solitary and periodic IAW solutions corresponding to nonlinear homoclinic and periodic orbits. These types of nonlinear IAW solutions are already reported [80]. Also, we observe existence of supernonlinear solitary and periodic IAWs solutions corresponding to supernonlinear homoclinic and periodic orbit, respectively [71]. Figure 1 (c) show phase profiles of the CDS (11 ) for q =0.77 with α =0.62, γ =0.06, μc =0.03,β =0.07 and V =0.705. The phase portrait manifested in figure 1 (c) for the range 0<q <1 presents a saddle at the origin. Also, we perceive the existence of nonlinear solitary and periodic IAW solutions corresponding to nonlinear homoclinic and periodic orbits, and supernonlinear solitary and periodic IAWs solutions corresponding to supernonlinear homoclinic and periodic orbit, respectively. Figures 1 (d) and (e) display phase profiles of the CDS (11 ) for q =1.7 and q =1.8, respectively, with $\alpha =0.9,\,\gamma =0.01,$μc =0.06, β =0.01 and V =0.776 in the subextensive range of parameter q . It is evident from figures 1 (d) and (e) that the system (11 ) in the range of q >1 also shows the existence of nonlinear and supernonlinear periodic and solitary waves.

Figure 1.

New window|Download| PPT slide
Figure 1.Phase profiles of the CDS (11 ) for (a) q =−0.17 and (b) q =−0.26 with α =0.9, γ =0.01, μc =0.06, β =0.05 and V =0.5, (c) q =0.77 with α =0.62, γ =0.06, μc =0.03, β =0.07 and V =0.705, (d) q =1.7 and (e) q =1.8 with α =0.6, γ =0.01, μc =0.06, β =0.01 and V =0.776.


Now, we show the effects of system parameters on superperiodic wave solutions with respect to the distinct ranges of q . Hence, we consider parametric sets of figure 1 (a) for the range −1<q <0, figure 1 (c) for the range 0<q <1 and figure 1 (e) for the range q >1. In figures 2 (a)–(c), we display superperiodic wave solutions with effects of parameters q, β and γ, respectively, with parametric set of figure 1 (a) for the range −1<q <0. Similarly, we also show effects of q, β and γ on superperiodic wave solutions in figures 2 (d)–(f) associated to figure 1 (c) for the range 0<q <1, and in figures 2 (g)–(i) associated to figure 1 (e) for the range q >1. In this case, we observe that with increase in q, the amplitude of supernonlinear periodic IAWs decreases in the range −1<q <0(see figure 2 (a)) and increases in the ranges 0<q <1 and q >1 (see figures 2 (d) and (g)). With increase in β, the amplitude of supernonlinear periodic IAWs increases in the ranges −1<q <0 and 0<q <1 (see figures 2 (b) and (e)) and, decreases in the range q >1 (see figure 2 (h)). When γ increases, the amplitude of supernonlinear periodic IAWs increases in the range −1<q <0 (see figure 2 (c)) and decreases smoothly in the ranges 0<q <1 and q >1 (see figures 2 (f) and (i)).

Figure 2.

New window|Download| PPT slide
Figure 2.(a)–(c) Effects of q, β and γ, respectively, on superperiodic wave solutions corresponding to supernonlinear periodic trajectories shown in figures 1 (a) for the range −1<q <0. (d)–(f) Effects of q, β and γ, respectively, on superperiodic wave solutions associated with figures 1 (c) for the range 0<q <1. (g)–(i) Effects of q, β and γ, respectively, on superperiodic wave solutions associated with figures 1 (e) for the range q >1.


In figures 3 (a)–(c), we display nonlinear periodic wave solutions with effects of parameters q, β and γ, respectively, with parametric set of figure 1 (a) for the range −1<q <0. Similarly, we also show effects of q, β and γ on nonlinear periodic wave solutions in figures 3 (d)–(f) associated to figure 1 (c) for the range 0<q <1, and in figures 3 (g)–(i) associated to figure 1 (e) for the range q >1. Here, we see that with increase in q, the nonlinear periodic IAWs become smooth in the ranges −1<q <0 and q >1 (see figures 3 (a) and (g)) and spiky in the ranges 0<q <1 (see figure 3 (d)). With increase in β, the nonlinear periodic IAWs become smooth in the ranges −1<q <9 and 0<q <1 (see figures 3 (b) and (e)) and spiky in the range q >1 (see figure 3 (h)). However, when γ increases, the nonlinear periodic IAWs become spiky in the ranges −1<q <0 and 0<q <1 (see figures 3 (c) and (f)) and smooth in the range q >1 (see figure 3 (i)).

Figure 3.

New window|Download| PPT slide
Figure 3.(a)–(c) Effects of q, β and γ, respectively, on nonlinear periodic wave solutions corresponding to nonlinear periodic trajectories shown in figures 1 (a) for the range −1<q <0. (d)–(f) Effects of q, β and γ, respectively, on nonlinear periodic wave solutions associated with figures 1 (c) for the range 0<q <1. (g)–(i) Effects of q, β and γ, respectively, on nonlinear periodic wave solutions associated with figures 1 (e) for the range q >1.


4. Quasiperiodic and chaotic orbits for IAWs

Recently, the effect of the source perturbation on nonlinear waves under Gaussian shape is reported in the work [81]. But, the nonlinear perturbation as external forcing, can be of various types [8284]. In this study, we consider a nonlinear perturbation as f0 cos(ωξ). In presence of external periodic forcing f0 cos(ωξ), the system (11 ) can be expressed as:$ \begin{eqnarray}\left\{\begin{array}{l}\tfrac{{\rm{d}}\phi }{{\rm{d}}\xi }=z,\\ \tfrac{{\rm{d}}z}{{\rm{d}}\xi }=a\phi +b{\phi }^{2}+c{\phi }^{3}+d{\phi }^{4}+{f}_{0}\cos (\omega \xi )\end{array}\right.,\end{eqnarray}$ where f0 and ω denote strength and frequency of the external perturbation, respectively.

In figure 4, we display dynamical properties corresponding to the perturbed dynamical system (PDS) (17 ) by changing intensity of external periodic force (f0 ) and frequency (ω) considering different set of values from figures 1 (b)–(d). In figures 4 (a)–(c), we consider set of parameters same as figure 1 (b) with f0 =0.005 at initial condition (0.17, 0) by varying (a) ω =0.001, (b) ω =0.005 and (c) ω =0.05. We observe that orbits in phase spaces figures 4 (a)–(c) transit from regular pattern covering surface to irregular patterns. Hence, it is observed that dynamical motion of the PDS (17 ) transits slowly from quasiperiodic behavior (figure 4 (a)) to chaotic behavior (figure 4 (c)) as we increase frequency of external periodic force (ω). In figures 4 (d) and (e), we show phase space for (d) ω =0.5 and (e) ω =2 of the PDS (17 ) considering f0 =0.005 at initial condition (−0.004, 0) with set of parameters same as figure 1 (d). It is perceived that as we increase frequency (ω), orbits in phase spaces (figures 4 (d)–(e)) cover surface of space forming pattern of torus. Hence, we observe quasiperiodic behaviors in the subextensive range of parameter q . In figures 4 (f)–(g), we show phase space of the PDS (17 ) varying (f) f0 =0.005 and (g) f0 =0.01 for the range 0<q <1 and considering ω =0.08 with set of parameters same as figure 1 (c). It is observed that the PDS (17 ) exhibits two types of quasiperiodic behaviors. From figure 4 (f),it is transparent that there exists torus-like band forming quasiperiodic patterns. On the other hand, figure 4 (g) shows that there exist triple connected torus exhibiting complex quasiperiodic behaviors. However, in figures 4 (h)–(l), we show phase space of the perturbed dynamical system (17 ) varying (h) ω =0.01, (i) ω =0.1, (j) ω =0.57, (k) ω =0.6, and (l) ω =0.67 at initial condition (0.3, 0) for f0 =0.005 with the set of parameters same as figure 1 (c) in the range 0<q <1. It is observed that the PDS (17 ) exhibits continuous change in dynamical behaviors for different values of ω . In figure 4 (h), we display the phase space in which orbits filling around like a toroid exhibit quasiperiodic behavior. The quasiperiodic behaviors shown in figures 4 (i) and (j) are patterns of single and double twisted toruses, respectively. As we increase value of ω, the distortion of torus continues to give more complex structures, such as chaos in phase space (figure 4 (k)). In phase space figure 4 (l), orbits spread out steadily around like a toroid covering the surface densely.

Figure 4.

New window|Download| PPT slide
Figure 4.Phase spaces of system (17 ) with f0 =0.005 at initial condition (φ, z)=(0.17,0) for (a) ω =0.001, (b) ω =0.005 and (c) ω =0.05 while other parameters are kept the same as in figure 1 (b). Phase spaces (d) ω =0.5 and (e) ω =2 with f0 =0.005 at initial condition (φ, z)=(−0.004, 0) while other system parameters are the same as in figure 1 (d). Phase spaces (f) f0 =0.005 and (g) f0 =0.01 considering fixed value ω =0.08 while (h) ω =0.01, (i) ω =0.1, (j) ω =0.57, (k) ω =0.6, and (l) ω =0.67 at initial condition (φ, z)= (0.3, 0) while other system parameters are the same as in figure 1 (c).


To set up the existence of chaos, the main effectual approach is the alteration of the Lyapunov exponents. The Lyapunov exponent plays as an excellent measure for determining the rate of exponential convergence and divergence of nearby orbits in phase space. For any orbit ${\bf{X}}({\bf{t}})=({x}_{1}(t),{x}_{2}(t),\ldots ,{x}_{n}(t))$ of a continuous dynamical system with dimension n and consider in the phase space that two orbits beginning with initial conditions X0 and ${{\bf{X}}}_{0}^{{\prime} }={{\bf{X}}}_{0}+\delta {{\bf{X}}}_{0}$ develop with time as the vectors X (t) and X (t)+δX (t), respectively. Then, one can define the Lyapunov exponent [85] corresponding to that orbits as$ \begin{eqnarray}\lambda ({{\bf{X}}}_{0},\delta {\bf{X}})=\mathop{\mathrm{lim}}\limits_{t\to \infty }\displaystyle \frac{1}{t}\mathrm{log}\left(\displaystyle \frac{{\rm{d}}({{\bf{X}}}_{0},t)}{{\rm{d}}({{\bf{X}}}_{0},0)}\right).\end{eqnarray}$ Here, d(X0, t) is distance between orbits X (t) and ${{\bf{X}}}^{{\prime} }(t)$ . The main feature of Lyapunov exponent shows that for a chaotic dynamical system, at least one Lyapunov exponent is positive while, it is nonpositive for nonchaotic systems. Hence, the positive value of Lyapunov exponents indicates a chaotic motion. Using Wolf algorithm [86], we examine Lyapunov exponent for system (17 ). Here, we observe that the nonlinear dynamics of the dynamical system (17 ) varies effectively due to the impact of external force. Therefore, we investigate this change of system (17 ) by plotting the Lyapunov exponents with respect to frequency (ω) and the forcing parameter (f0 ) of an external periodic force in figure 5 corresponding to the phase space reported in figure 4 (c). The positive values of the Lyapunov exponents in figure 5 show the evidence of chaos for the different values of ω and f0 . However, the Lyapunov exponents at ω =0.05045, 0.05425 (in figure 5 (a)) and f0 =0.0048, 0.00492 (in figure 5 (b)) tend to be zero describing the presence of quasiperiodic orbits of the dynamical system (17 ). It is important to notice that the system (17 ) still preserves the conservative behavior in presence of external force.

Figure 5.

New window|Download| PPT slide
Figure 5.Plots of Lyapunov exponents with respect to (a) ω and (b) f0 corresponding to the phase space in figure 4 (c).


5. Coexisting orbits for IAWs

In this section, we investigate existence of coexisting orbits for IAWs in presence of the nonlinear controller f0 cos(ωξ). We observe that the system (17 ) can show coexisting periodic, quasiperiodic and chaotic orbits for IAWs.

In figure 6, existence of multistability for IAWs in the perturbed dynamical system (17 ) is shown with the help of numerical simulation. Figure 6 (a) shows coexistence of quasiperiodic (blue curve) and chaotic orbits (purple curve) for IAWs by changing initial condition from (0.17, 0) for the blue curve to (0.001, 0) for the purple curve, when other parameters are kept same as figure 4 (b). This shows that the system (17 ) is sensitive to initial condition. In figure 6 (b), by varying the initial condition from (−0.004, 0) to (0.5, 0), we observe that the system (17 ) transits from quasiperiodic to chaotic orbits for IAWs showing sensitive dependency to the initial condition. Therefore, figure 6 (b) depicts coexistence of quasiperiodic orbits (green curve) for IAWs at initial condition (−0.004, 0) and chaotic orbits (blue curve) for IAWs at initial condition (0.5, 0), while other parameters are same as figure 4 (d). Figure 6 (c) displays coexistence of periodic orbit for IAWs at initial condition (0.09, 0) (green curve) and on changing the initial condition to (0.3, 0), orbits in phase space fill torus spreading over the surface i.e. quasiperiodic orbits for IAWs is achieved (red curve) for parametric values same as figure 4 (h). Figure 6 (d) is enlarged view of the quasiperiodic orbits for IAWs shown in figure 6 (c) with the red curve.

Figure 6.

New window|Download| PPT slide
Figure 6.Phase space of the system (17 ) with ω =0.05 for (a) f0 =0.001 at initial conditions (φ, z)=(0.17, 0) for the blue curve and (φ, z)=(0.001, 0) for the purple curve while other parameters are the same as in figure 4 (b). (b) Phase space for f0 =0.5 at initial conditions (φ, z)=(−0.004, 0) for the green curve and (φ, z)=(0.5, 0) for the dark blue curve while other parameters are the same as in figure 4 (d). Phase space considering f0 =0.005 and ω =1 at (c) initial condition (φ, z)=(0.09, 0) for the green curve and (φ, z)=(0.3, 0) for the red curve while other parameters are the same as in figure 4 (h). (d) Enlarged view of the phase space shown with the red curve in figure 6 (c).


In figure 7, time series plots of IAWs are shown for the system (17 ) corresponding to coexisting orbits shown in phase spaces of figure 6 . Figure 7 (a) represents time series plot of coexisting quasiperiodic (dark blue color) and chaotic orbits (purple color) for IAWs shown in figure 6 (a). Figure 7 (b) shows time series plot of coexisting quasiperiodic (green color) and chaotic orbits (blue color) for IAWs presented in figure 6 (b). Time series plots in figures 7 (a) and (b) clearly show dynamical properties namely, quasiperiodic and chaotic orbits for IAWs existing in the perturbed system (17 ). However, figure 7 (c) displays time series plot of coexisting periodic (green curve) and quasiperiodic orbits (red color) for IAWs shown in figures 6 (c) and (d).

Figure 7.

New window|Download| PPT slide
Figure 7.(a) and (b) Time series plots corresponding to co-existing orbits for IAWs shown in figures 6 (a) and (b). (c) Time series plot corresponding to co-existing orbits for IAWs shown in figures 6 (c) (green curve) and (d) (red curve).


6. Conclusions

Superperiodicity, chaos and coexisting orbits of the nonlinear IAWs have been reported in a multi-component plasma containing fluid ions, Maxwellian hot positrons, cold and hot electrons with q -nonextensive distribution. The dynamics of IA superperiodic waves, which can be quantified by the phase portrait and the time series plots have been discussed. The significant impacts of system parameters like, nonextensive parameter (q), temperature ratio of cold and hot electrons (β) and the ratio of effective temperature and positron temperature (γ) have been observed on superperiodicity and nonlinear periodicity of IAWs. It has been encountered that superperiodic IAWs become smooth with increase in q and spiky with increase in β and γ for the range −1<q <0. The superperiodic IAWs have been spiky with increase in q and β, and less spiky with increase in γ for the range 0<q <1. However for the range q >1, the superperiodic IAWs have been spiky with increase with q and, smooth when β and γ increase. Furthermore, it has been observed that nonlinear periodic IAWs show similar features with increase in q, β and γ . The nonlinear periodic IAW has been smooth with increase in q and β while, has been sharp with increase in γ in the range −1<q <0. However, in the range 0<q <1, the nonlinear periodic IAW has been sharp and with increase in q and γ whereas, it has been smooth with increase in β . For the range q >1, it has been observed that the nonlinear periodic IAWs have been smooth for increase in q and γ, and sharp with increase in β . Effect of a nonlinear periodic perturbation has been presented on various types of orbits, like periodic, quasiperiodic and chaotic orbits for IAWs. The coexistence of such orbits for IAWs has been reported for suitable parametric region in their physical domain with different initial values.

Acknowledgments

The first author is thankful to SMIT and SMU for research support (6100/SMIT/R&D/Project/05/2018).


Reference By original order
By published year
By cited within times
By Impact factor

Gibbons G W Hawking S W Siklos S 1983 The Very Early Universe Cambridge Cambridge University Press
DOI:10.1017/S1323358000017744 [Cited within: 1]

Schlickeiser R Shukla P K 2003 Astrophys. J. 599 L57
DOI:10.1017/S1323358000017744 [Cited within: 1]

Michel F C 1991 Theory of Neutron Star Magnetospheres Chicago, IL University of Chicago Press
[Cited within: 1]

Miller H R 1998 Active Galactic Nuclei Georgia State University, Atlanta, GA 28–39 October 1987 Berlin Springer
[Cited within: 1]

Tandberg-Hansen E Emshie A G 1988 The Physics of Solar Flares Cambridge Cambridge University Press 124 p 124
[Cited within: 2]

Piran T 1999 Phys. Rep. 314 575 667 575–667
DOI:10.1016/S0370-1573(98)00127-6 [Cited within: 2]

Surko C M Lemvethal M Crane W S Passner A Wyocki F J Murphy T J Strachan J Rowan W L 1986 Rev. Sci. Instrum. 57 1862
DOI:10.1063/1.1139154 [Cited within: 1]

Tinkle M D Greaves R G Surko C M Spencer R L Mason G W 1994 Phys. Rev. Lett. 72 352
DOI:10.1103/PhysRevLett.72.352

Greaves R G Surko C M 1995 Phys. Rev. Lett. 75 3846
DOI:10.1103/PhysRevLett.75.3846 [Cited within: 1]

Berezhiani V I Tsintsadze L N Shukla P K 1992 J. Plasma Phys. 48 139
DOI:10.1017/S0022377800016421 [Cited within: 1]

Surko C M Murphy T J 1990 Phys. Fluids B 2 1372
DOI:10.1063/1.859558 [Cited within: 2]

Berezhiani V I El-Ashry M Y Mofiz U A 1994 Phys. Rev. E 50 448
DOI:10.1103/PhysRevE.50.448

Popel S I Vladimirov S V Shukla P K 1995 Phys. Plasmas 2 716
DOI:10.1063/1.871422

Yu M Y Shukla P K Stenflo L 1986 Astrophys. J. 309 L63
DOI:10.1086/184761 [Cited within: 1]

Kakati H Goswami K S 2000 Phys. Plasmas 7 808
DOI:10.1063/1.873875 [Cited within: 1]

Alinejad H Sobhanian S Mahmood J 2006 Phys. Plasmas 13 012304
DOI:10.1063/1.2158148

Tiwari R S Kaushik A Mishra M K 2007 Phys. Lett. A 365 335
DOI:10.1016/j.physleta.2007.01.020

Sabry R Moslem W M Haas F Ali S Shukla P K 2008 Phys. Plasmas 15 122308
DOI:10.1063/1.3037265

El-Shamy E F El-Bedwehy N A 2010 Phys. Lett. A 374 4425 4429 4425–9
DOI:10.1016/j.physleta.2010.09.002

Saha A Tamang J 2017 Phys. Plasmas 24 082101
DOI:10.1063/1.4994396 [Cited within: 1]

Dubinov A E Kolotkov D Y 2012 IEEE Trans Plasma Sci. 40 1429
DOI:10.1109/TPS.2012.2189026 [Cited within: 1]

Dubinov A E Kolotkov D Y 2012 High Energy Chem. 46 349
DOI:10.1134/S0018143912060033 [Cited within: 1]

Dubinov A E Kolotkov D Y 2012 Plasma Phys. Rep. 38 909
DOI:10.1134/S1063780X12100054 [Cited within: 1]

Dubinov A E Kolotkov D Y Sazonkin M A 2012 Plasma Phys. Rep. 38 833 844 833–44
DOI:10.1134/S1063780X12090036 [Cited within: 1]

Dubinov A E Sazonkin M A 2013 Phys. Wave Phenom. 21 118 128 118–28
DOI:10.3103/S1541308X13020039 [Cited within: 1]

Verheest F 2009 Phys Plasma 16 013704
DOI:10.1063/1.3059411 [Cited within: 1]

Baluku T K Hellberg M A Verheest F 2010 Europhys. Lett. 91 15001
DOI:10.1209/0295-5075/91/15001

Verheest F 2011 Phys Plasma 18 083701
DOI:10.1063/1.3610560

Das A Bandyopadhyay A Das K P 2012 J. Plasma Phys 78 149
DOI:10.1017/S002237781100050X [Cited within: 1]

Verheest F Hellberg M A Kourakis I 2013 Phys Plasmas 20 012302
DOI:10.1063/1.4775085 [Cited within: 1]

Verheest F Hellberg M A Kourakis I 2013 Phys Rev. 87 043107
DOI:10.1103/PhysRevE.87.043107

Verheest F Hellberg M A Kourakis I 2013 Phys Plasmas 20 082309
DOI:10.1063/1.4775085 [Cited within: 1]

Maharaja S K Bharuthram R Singh S V Lakhina G S 2013 Phys Plasmas 20 083705
DOI:10.1063/1.4818439 [Cited within: 1]

Rufai O R Bharuthram R Singh S V Lakhina G S 2014 Phys Plasmas 21 082304
DOI:10.1063/1.4891877 [Cited within: 1]

Lakhina G S Singh S V Kakad A P 2014 Phys. Plasmas 21 062311
DOI:10.1063/1.4884791 [Cited within: 1]

El-Wakil S A Abulwafa E M Elhanbaly A A 2017 Phys. Plasmas 24 073705
DOI:10.1063/1.4991406 [Cited within: 1]

Tsallis C 1988 J. Stat. Phys. 52 479
DOI:10.1007/BF01016429 [Cited within: 1]

Tsallis C 1994 New Trends in Magnetism, Magnetic Materials and Their Applications Moran-Lopez J L Sanchez J M New York Plenum 451 p 451
DOI:10.1007/978-1-4899-1334-0 [Cited within: 1]

Tsallis C 1995 Chaos Solitons Fractals 6 539
DOI:10.1016/0960-0779(95)80062-L [Cited within: 1]

Lima J A S Silva R Jr. Santos J 2000 Phys. Rev. E 61 3260
DOI:10.1103/PhysRevE.61.3260 [Cited within: 2]

Munoz V 2006 Nonlinear Process. Geophys. 13 237
DOI:10.5194/npg-13-237-2006 [Cited within: 1]

Latora V Rapisarda A Tsallis C 2001 Phys. Rev. E 64 056134
DOI:10.1103/PhysRevE.64.056134 [Cited within: 1]

Taruya A Sakagami M A 2003 Phys. Rev. Lett. 90 181101
DOI:10.1103/PhysRevLett.90.181101 [Cited within: 1]

Liu J M De Groot J S Matt J P Johnston T W Drake R P 1994 Phys. Rev. Lett. 72 2717
DOI:10.1103/PhysRevLett.72.2717 [Cited within: 1]

Du J 2004 Phys. Lett. A 329 262
DOI:10.1016/j.physleta.2004.07.010 [Cited within: 1]

Liyan L Du J 2008 Physica A 387 4821
DOI:10.1016/j.physa.2008.04.016

Saha A Chatterjee P 2015 Eur. Phys. J. D 69 203
DOI:10.1140/epjd/e2015-60115-7

Saha A Chatterjee P Wong C S 2015 Braz. J. Phys. 45 656
DOI:10.1007/s13538-015-0358-3

Tribeche M Djebarni L Amour L R 2010 Phys. Plasmas 17 042114
DOI:10.1063/1.3374429 [Cited within: 1]

Mukherjee S Palit S K Banerjee S Ariffin M Rondoni L Bhattacharya D 2015 Physica A 439 93 102 93–102
DOI:10.1016/j.physa.2015.07.030 [Cited within: 1]

Banerjee S Palit S K Mukherjee S Ariffin M Rondoni L 2016 Chaos 26 033105
DOI:10.1063/1.4941374 [Cited within: 1]

Banerjee S Saha P Chowdhury A R 2004 Chaos 14 347 357 347–57
DOI:10.1063/1.1755179 [Cited within: 1]

Theesar S J S Balasubramaniam P Banerjee S 2009 Chaos 19 013102
DOI:10.1063/1.3054711

Theesar S J S Balasubramaniam P Banerjee S 2010 Chaos 20 029902
DOI:10.1063/1.3054711

Theesar S J S Balasubramaniam P Banerjee S 2012 Chaos 22 013102
DOI:10.1063/1.3054711

Banerjee S Rondoni L Mukhopadhyay S 2011 Opt. Commun. 284 4623 4634 4623–34
DOI:10.1016/j.optcom.2011.06.009

Banerjee S Mukhopadhyay S Rondoni L 2012 Opt. Lasers Eng. 50 950 957 950–7
DOI:10.1016/j.optlaseng.2012.02.009

Rahim M F A Natiq H Fataf N A A Banerjee S 2019 Eur. Phys. J. Plus 134 499
[Cited within: 1]

Angeli D Ferrell J E Sontag E D 2004 Proc. Natl Acad. Sci. 101 1822 1827 1822–7
DOI:10.1073/pnas.0308265100 [Cited within: 1]

Arecchi F Meucc R Puccioni G Tredicce J 1982 Phys. Rev. Lett. 49 1217
DOI:10.1103/PhysRevLett.49.1217

Li C Sprott J C 2014 Int. J. Bifurcation Chaos 24 1450131
DOI:10.1142/S0218127414501314

Li C Sprott J C Hu W Xu Y 2017 Int. J. Bifurcation Chaos 27 1750160
DOI:10.1142/S0218127417501607

Wang M Deng Y Liao X Li Z Ma M Zeng Y 2019 Int. J. Non-Linear Mech. 111 149
DOI:10.1016/j.ijnonlinmec.2019.02.009

Natiq H Banerjee S Ariffin M Said M 2019 Chaos 29 011103
DOI:10.1063/1.5079886 [Cited within: 1]

Natiq H Banerjee S Misra A P Said M R M 2019 Chaos Solitons Fractals 122 58 68 58–68
DOI:10.1016/j.chaos.2019.03.009 [Cited within: 1]

Strogatz S H 2007 Nonlinear Dynamics and Chaos Boulder, CO Westview Press
DOI:10.1201/9780429492563 [Cited within: 2]

Samanta U K Saha A Chatterjee P 2013 Phys. Plasma 20 022111
DOI:10.1063/1.4791660 [Cited within: 1]

Saha A Pal N Chatterjee P 2014 Phys. Plasma 21 102101
DOI:10.1063/1.4896715 [Cited within: 1]

Selim M M El-Depsy A El-Shamy E F 2015 Astrophys. Space Sci. 360 66
DOI:10.1007/s10509-015-2574-2 [Cited within: 1]

Rufai O R Bharuthram R Singh S V Lakhina G S 2016 Adv. Spaces Res. 57 813 820 813–20
DOI:10.1016/j.asr.2015.11.021 [Cited within: 1]

Maharaj S K Bharuthram R Singh S V Lakhina G S 2013 Phys. Plasmas 20 083705
DOI:10.1063/1.4818439

Singh S V Lakhina G S 2015 Commun. Nonlinear Sci. Numer. Simulat. 23 274 281 274–81
DOI:10.1016/j.cnsns.2014.11.017

Mandi L Saha A Chatterjee P 2019 Adv. Space Res. 64 427 435 427–35
DOI:10.1016/j.asr.2019.04.028 [Cited within: 1]

Saha A Tamang J 2019 Adv. Space Res. 63 1596 1606 1596–606
DOI:10.1016/j.asr.2018.11.010 [Cited within: 2]

Prasad P K Sarkar S Saha A Mondal K K 2019 Braz. J. Phys. 49 698
DOI:10.1007/s13538-019-00697-y

Tamang J Saha A 2019 Z. Nat.forsch. A 74 499
DOI:10.1515/zna-2018-0476 [Cited within: 2]

Saha A Pradhan B Banerjee S 2020 Phys. Scr. 95 5
DOI:10.1088/1402-4896/ab7052 [Cited within: 1]

Sabry R 2009 Phys. Plasmas 16 072307
DOI:10.1063/1.3175693 [Cited within: 1]

Bains A S Tribeche M Gill T S 2011 Phys. Plasmas 18 022108
DOI:10.1063/1.3554658 [Cited within: 1]

Chow S N Hale J K 1982 Elements of Nonlinear Analysis, Methods of Bifurcation Theory New York Springer
DOI:10.1007/978-1-4613-8159-4 [Cited within: 1]

Saha A 2012 Commun. Nonlinear Sci. Numer. Simul. 17 3539
DOI:10.1016/j.cnsns.2012.01.005 [Cited within: 1]

Guckenheimer J Holmes P J 1983 Nonlinear Oscillations, Dynamical Systems and Bifurcations of Vector Fields New York Springer
DOI:10.1007/978-1-4612-1140-2 [Cited within: 2]

Saha A Chatterjee P 2014 Astrophys Space Sci. 350 631 636 631–6
DOI:10.1007/s10509-014-1796-z [Cited within: 1]

Sen A Tiwari S Mishra S Kaw P 2015 Adv. Space Res. 56 429
DOI:10.1016/j.asr.2015.03.021 [Cited within: 1]

Ali R Saha A Chatterjee P 2017 Phys. Plasmas 24 122106
DOI:10.1063/1.4994562 [Cited within: 1]

Chatterjee P Ali R Saha A 2018 Z. Nat.forsch. A 73 151 159 151–9
DOI:10.1515/zna-2017-0358

Zhen H Tian B Zhong H Sun W Li M 2013 Phys. Plasmas 20 082311
DOI:10.1063/1.4818508 [Cited within: 1]

Lakshmanan M Rajasekar S 2003 Nonlinear Dynamics Heidelberg Springer
DOI:10.1007/978-3-642-55688-3 [Cited within: 1]

Wolf A Swift J B Swinney H L Vastano J A 1985 Physica D 16 285 317 285–317
DOI:10.1016/0167-2789(85)90011-9 [Cited within: 1]

相关话题/Superperiodicity chaos coexisting